>>20034213
Medicinal Genomics MIP DNA Purification Kit used
https://store.medicinalgenomics.com/SenSATIVAx-DNA-Extraction-Kit-200-reactions_2
he CTAB/Chloroform/SPRI based DNA/RNA isolation methods are described above.
Cycling conditions
These conditions work for both qPCR and RT-qPCR. Note: The 50°C RT step can be skipped with qPCR. The MGC qPCR MasterMix kits used have a hot start enzyme which are unaffected by this 50°C step. For the sake of controlling RNA to DNA comparisons, we have put qPCR and RT-qPCR assays on the same plate and run the below program with the RT step included for all samples.
Cycling Conditions used for qPCR and RT-qPCR
Sequences of amplicons for gBlock Positive Controls. Ori = 106bp, Spike = 114bp.
Ori target
Spike target
Sequencing Data
Raw Illumina Reads RNA-seq
Read files are run through sha256 (Hash and stash) and etched onto the DASH blockchain. The sha256 hash of the read file is spent into the OP_RETURN of an immutable ledger. If the hash of the file doesn’t match the hash in these transactions, the file has been tampered with.
Megahit Assemblies
Pfizer Vial 1
Pfizer Vial 2
Moderna Vial 1
Moderna Vial 2
Illumina Reads mapped back to Megahit Assemblies
Q30 Filtered Illumina Reads (use these for transcriptional error rate estimates)
FastQ-Filter download: usagefastq-filter -e 0.001 -o output.fastq input.fastq
Q30 BAM files. Q30 Reads mapped against Megahit assemblies
IGVtools error by base on q30 reads
Fields = Position in contig, Positive stand (+)A, +C, +G, +T, +N, +Deletion, +Insertion, Negative strand -A, -C, -G, -T, -N, -Deletion, -Insertion
· Moderna Vial 1
· Moderna Vial 2
· Pfizer Vial 1
· Pfizer Vial 2
Analysis pipeline
Reads were demultiplexed and processed with
Trimgalore - Removes Illumina Sequencing adaptors.
Megahit- assembles reads into contigs.
Megahit for SARs-CoV-2
Samtools- generates BAM files for viewing in IGV.
Samtools stats used to calculate outie reads.
BWA-mem- Short read mapper used to align reads back to the assembled references.
SnapGene software- (www.snapgene.com)- Used to visualize and annotate expression vectors
IGV- Integrated Genome Viewer used to visualize Illumina sequencing reads.
RNase Treated Libraries-BAM files
contig specific BAM files were created using samtools
samtools view -h input.bam contig_name -O BAM contig.bam; samtools index contig.bam;
Samtools stats run on a each contig in each assembly.
for out_prefix in ls *.sort.bam | perl -pe "s/.sort.bam//"
; do mkdir -p ${out_prefix}-samtools-stats; for contig in samtools view -H ${out_prefix}.sort.bam | grep "^@SQ" | cut -f 2 | perl -pe "s/SN://"
; do echo "Now calculating stats for ${contig}/$out_prefix…"; samtools stats ${out_prefix}.sort.bam $contig ${out_prefix}-samtools-stats/${contig}-samtools-stats.txt; done; done
Author contributions
KJM- constructed the sequencing libraries, designed the qPCR assays, ran Qubit™ 3s and Agilent Tape Station™ and performed the analysis, drafted the manuscript.
YH-Optimized DNA isolations, Tape Station™ and qPCR results.
SM, LTK- assisted in demultiplexing and trimming the reads and assembly troubleshooting
Conflicts of interest- Authors of this paper are employees of Medicinal Genomics which manufacturers some of the qPCR and DNA isolation kits used in this study.