Extraction-free clinical detection of SARS-CoV-2 virus from saline gargle samples using Hamilton STARlet liquid handler.

As part of the COVID-19 pandemic, clinical laboratories have been faced with massive increases in testing, resulting in sample collection systems, reagent, and staff shortages. We utilized self-collected saline gargle samples to optimize high throughput SARS-CoV-2 multiplex polymerase chain reaction (PCR) testing in order to minimize cost and technologist time. This was achieved through elimination of nucleic acid extraction and automation of sample handling on a widely available robotic liquid handler, Hamilton STARlet. A customized barcode scanning script for reading the sample ID by the Hamilton STARlet's software system was developed to allow primary tube sampling. Use of pre-frozen SARS-CoV-2 assay reaction mixtures reduced assay setup time. In both validation and live testing, the assay produced no false positive or false negative results. Of the 1060 samples tested during validation, 3.6% (39/1060) of samples required retesting as they were either single gene positive, had internal control failure or liquid aspiration error. Although the overall turnaround time was only slightly faster in the automated workflow (185 min vs 200 min), there was a 76% reduction in hands-on time, potentially reducing staff fatigue and burnout. This described process from sample self-collection to automated direct PCR testing significantly reduces the total burden on healthcare systems in terms of human resources and reagent requirements.

Automated Library Construction and Analysis for High-Throughput Nanopore Sequencing of SARS-CoV-2.

BACKGROUND: To support the implementation of high-throughput pipelines suitable for SARS-CoV-2 sequencing and analysis in a clinical laboratory, we developed an automated sample preparation and analysis workflow. METHODS: We used the established ARTIC protocol with approximately 400 bp amplicons sequenced on Oxford Nanopore's MinION. Sequences were analyzed using Nextclade, assigning both a clade and quality score to each sample. RESULTS: A total of 2179 samples on twenty-five 96-well plates were sequenced. Plates of purified RNA were processed within 12 h, sequencing required up to 24 h, and analysis of each pooled plate required 1 h. The use of samples with known threshold cycle (Ct) values enabled normalization, acted as a quality control check, and revealed a strong correlation between sample Ct values and successful analysis, with 85% of samples with Ct < 30 achieving a "good" Nextclade score. Less abundant samples responded to enrichment with the fraction of Ct > 30 samples achieving a "good" classification rising by 60% after addition of a post-ARTIC PCR normalization. Serial dilutions of 3 variant of concern samples, diluted from approximately Ct = 16 to approximately Ct = 50, demonstrated successful sequencing to Ct = 37. The sample set contained a median of 24 mutations per sample and a total of 1281 unique mutations with reduced sequence read coverage noted in some regions of some samples. A total of 10 separate strains were observed in the sample set, including 3 variants of concern prevalent in British Columbia in the spring of 2021. CONCLUSIONS: We demonstrated a robust automated sequencing pipeline that takes advantage of input Ct values to improve reliability.

Optimization of magnetic bead-based nucleic acid extraction for SARS-CoV-2 testing using readily available reagents.

The COVID-19 pandemic has highlighted the need for generic reagents and flexible systems in diagnostic testing. Magnetic bead-based nucleic acid extraction protocols using 96-well plates on open liquid handlers are readily amenable to meet this need. Here, one such approach is rigorously optimized to minimize cross-well contamination while maintaining sensitivity.

Transcription of foreign DNA in Escherichia coli.

Propagation of heterologous DNA in E. coli host cells is central to molecular biology. DNA constructs are often engineered for expression of recombinant protein in E. coli, but the extent of incidental transcription arising from natural regulatory sequences in cloned DNA remains underexplored. Here, we have used programmable microarrays and RT-PCR to measure, comprehensively, the transcription of H. influenzae, P. aeruginosa, and human DNA propagating in E. coli as bacterial artificial chromosomes. We find evidence that at least half of all H. influenzae genes are transcribed in E. coli. Highly transcribed genes are principally involved in energy metabolism, and their proximal promoter regions are significantly enriched with E. coli sigma(70) (also known as RpoD) binding sites. H. influenzae genes acquired from an ancient bacteriophage Mu insertion are also highly transcribed. Compared with H. influenzae, a smaller proportion of P. aeruginosa genes are transcribed in E. coli, and in E. coli there is punctuated transcription of human DNA. The presence of foreign DNA in E. coli disturbs the host transcriptional profile, with expression of the E. coli phage shock protein operon and the flagellar gene cluster being particularly strongly up-regulated. While cross-species transcriptional activation is expected to be enabling for horizontal gene transfer in bacteria, incidental expression of toxic genes can be problematic for DNA cloning. Ongoing characterization of cross-expression will help inform the design of biosynthetic gene clusters and synthetic microbial genomes.

Rebuilding microbial genomes.

Engineered microbes are of great potential utility in biotechnology and basic research. In principle, a cell can be built from scratch by assembling small molecule sets with auto-catalytic properties. Alternatively, DNA can be isolated or directly synthesized and molded into a synthetic genome using existing genomic blueprints and molecular biology tools. Activating such a synthetic genome will yield a synthetic cell. Here we examine obstacles associated with this latter approach using a model system whereby a donor genome from H. influenzae is fragmented, and the pieces are then modified and reassembled stepwise in an E. coli host cell. There are obstacles associated with this strategy related to DNA transfer, DNA replication, cross-talk in gene regulation and compatibility of gene products between donor and host. Encouragingly, analysis of gene expression indicates widespread transcription of H. influenzae genes in E. coli, and analysis of gap locations in H. influenzae and other microbial genome assemblies reveals few genes routinely incompatible with E. coli. In conclusion, rebuilding and booting a genome remains a feasible and pragmatic approach to creating a synthetic microbial cell.

A physical map of the highly heterozygous Populus genome: integration with the genome sequence and genetic map and analysis of haplotype variation.

As part of a larger project to sequence the Populus genome and generate genomic resources for this emerging model tree, we constructed a physical map of the Populus genome, representing one of the few such maps of an undomesticated, highly heterozygous plant species. The physical map, consisting of 2802 contigs, was constructed from fingerprinted bacterial artificial chromosome (BAC) clones. The map represents approximately 9.4-fold coverage of the Populus genome, which has been estimated from the genome sequence assembly to be 485 +/- 10 Mb in size. BAC ends were sequenced to assist long-range assembly of whole-genome shotgun sequence scaffolds and to anchor the physical map to the genome sequence. Simple sequence repeat-based markers were derived from the end sequences and used to initiate integration of the BAC and genetic maps. A total of 2411 physical map contigs, representing 97% of all clones assigned to contigs, were aligned to the sequence assembly (JGI Populus trichocarpa, version 1.0). These alignments represent a total coverage of 384 Mb (79%) of the entire poplar sequence assembly and 295 Mb (96%) of linkage group sequence assemblies. A striking result of the physical map contig alignments to the sequence assembly was the co-localization of multiple contigs across numerous regions of the 19 linkage groups. Targeted sequencing of BAC clones and genetic analysis in a small number of representative regions showed that these co-aligning contigs represent distinct haplotypes in the heterozygous individual sequenced, and revealed the nature of these haplotype sequence differences.

Constructing large DNA segments by iterative clone recombination.

Methods for constructing large contiguous segments of DNA will be enabling for Synthetic Biology, where the assembly of genes encoding circuits, biosynthetic pathways or even whole microbial organisms is of interest. Currently, in vitro approaches to DNA synthesis are adequate for generating DNAs that are up to 10s of kbp in length, and in vivo recombination strategies are more suitable for building DNA constructs that are 100 kbp or larger. We have developed a vector system for efficient assembly of large DNA molecules by iterative in vivo recombination of fosmid clones. Two custom fosmid vectors have been built, pFOSAMP and pFOSKAN, that support antibiotic switching. Using this technique we rebuilt two non-contiguous regions of the Haemophilus influenzae genome as episomes in recombinogenic Escherichia coli host cells. These regions together comprise190 kbp, or 10.4% of the H. influenze genome.

Sequencing and analysis of 10,967 full-length cDNA clones from Xenopus laevis and Xenopus tropicalis reveals post-tetraploidization transcriptome remodeling.

Sequencing of full-insert clones from full-length cDNA libraries from both Xenopus laevis and Xenopus tropicalis has been ongoing as part of the Xenopus Gene Collection Initiative. Here we present 10,967 full ORF verified cDNA clones (8049 from X. laevis and 2918 from X. tropicalis) as a community resource. Because the genome of X. laevis, but not X. tropicalis, has undergone allotetraploidization, comparison of coding sequences from these two clawed (pipid) frogs provides a unique angle for exploring the molecular evolution of duplicate genes. Within our clone set, we have identified 445 gene trios, each comprised of an allotetraploidization-derived X. laevis gene pair and their shared X. tropicalis ortholog. Pairwise dN/dS, comparisons within trios show strong evidence for purifying selection acting on all three members. However, dN/dS ratios between X. laevis gene pairs are elevated relative to their X. tropicalis ortholog. This difference is highly significant and indicates an overall relaxation of selective pressures on duplicated gene pairs. We have found that the paralogs that have been lost since the tetraploidization event are enriched for several molecular functions, but have found no such enrichment in the extant paralogs. Approximately 14% of the paralogous pairs analyzed here also show differential expression indicative of subfunctionalization.

A high-throughput screen identifying sequence and promiscuity characteristics of the loxP spacer region in Cre-mediated recombination.

BACKGROUND: Cre-loxP recombination refers to the process of site-specific recombination mediated by two loxP sequences and the Cre recombinase protein. Transgenic experiments exploit integrative recombination, where a donor plasmid carrying a loxP site and DNA of interest integrate into a recipient loxP site in a target genome. Unfortunately, integrative recombination is highly inefficient because the insert is flanked by two loxP sites, which themselves become targets for Cre and lead to subsequent excision of the insert. A small number of mutations have been discovered in parts of the loxP sequence, specifically the spacer and inverted repeat segments, that increase the efficiency of integrative recombination. In this study we introduce a high-throughput in vitro assay to rapidly detect novel loxP spacer mutants and describe the sequence characteristics of successful recombinants. RESULTS: We created synthetic loxP oligonucleotides that contained a combination of inverted repeat mutations (the lox66 and lox71 mutations) and mutant spacer sequences, degenerate at 6 of the 8 positions. After in vitro Cre recombination, 3,124 recombinant clones were identified by sequencing. Included in this set were 31 unique, novel, self-recombining sequences. Using network visualization tools, we recognized 12 spacer sets with restricted promiscuity. We observed that increased guanine content at all spacer positions save for position 8 resulted in increased recombination. Interestingly, recombination between identical spacers was not preferred over non-identical spacers. We also identified a set of 16 pairs of loxP spacers that reacted at least twice with another spacer, but not themselves. Further, neither the wild-type P1 phage loxP sequence nor any of the known loxP spacer mutants appeared to be kinetically favoured by Cre recombinase. CONCLUSION: This study approached loxP spacer mutant screening in an unbiased manner, assuming nothing about candidate loxP sites save for the conserved 4 and 5 spacer positions. Candidate sites were free to recombine with any other sequence in the pool of all possible sites. The subset of loxP sites identified here are candidates for in vivo serial recombination as they have already demonstrated limited promiscuity with other loxP spacer and stability in the presence of Cre.

Simple, robust methods for high-throughput nanoliter-scale DNA sequencing.

We have developed high-throughput DNA sequencing methods that generate high quality data from reactions as small as 400 nL, providing an approximate order of magnitude reduction in reagent use relative to standard protocols. Sequencing of clones from plasmid, fosmid, and BAC libraries yielded read lengths (PHRED20 bases) of 765 +/- 172 (n = 10,272), 621 +/- 201 (n = 1824), and 647 +/- 189 (n = 568), respectively. Implementation of these procedures at high-throughput genome centers could have a substantial impact on the amount of data that can be generated per unit cost.

Systematic recovery and analysis of full-ORF human cDNA clones.

The Mammalian Gene Collection (MGC) consortium (http://mgc.nci.nih.gov) seeks to establish publicly available collections of full-ORF cDNAs for several organisms of significance to biomedical research, including human. To date over 15,200 human cDNA clones containing full-length open reading frames (ORFs) have been identified via systematic expressed sequence tag (EST) analysis of a diverse set of cDNA libraries; however, further systematic EST analysis is no longer an efficient method for identifying new cDNAs. As part of our involvement in the MGC program, we have developed a scalable method for targeted recovery of cDNA clones to facilitate recovery of genes absent from the MGC collection. First, cDNA is synthesized from various RNAs, followed by polymerase chain reaction (PCR) amplification of transcripts in 96-well plates using gene-specific primer pairs flanking the ORFs. Amplicons are cloned into a sequencing vector, and full-length sequences are obtained. Sequences are processed and assembled using Phred and Phrap, and analyzed using Consed and a number of bioinformatics methods we have developed. Sequences are compared with the Reference Sequence (RefSeq) database, and validation of sequence discrepancies is attempted using other sequence databases including dbEST and dbSNP. Clones with identical sequence to RefSeq or containing only validated changes will become part of the MGC human gene collection. Clones containing novel splice variants or polymorphisms have also been identified. Our approach to clone recovery, applied at large scale, has the potential to recover many and possibly most of the genes absent from the MGC collection.

Development and application of a salmonid EST database and cDNA microarray: data mining and interspecific hybridization characteristics.

We report 80,388 ESTs from 23 Atlantic salmon (Salmo salar) cDNA libraries (61,819 ESTs), 6 rainbow trout (Oncorhynchus mykiss) cDNA libraries (14,544 ESTs), 2 chinook salmon (Oncorhynchus tshawytscha) cDNA libraries (1317 ESTs), 2 sockeye salmon (Oncorhynchus nerka) cDNA libraries (1243 ESTs), and 2 lake whitefish (Coregonus clupeaformis) cDNA libraries (1465 ESTs). The majority of these are 3' sequences, allowing discrimination between paralogs arising from a recent genome duplication in the salmonid lineage. Sequence assembly reveals 28,710 different S. salar, 8981 O. mykiss, 1085 O. tshawytscha, 520 O. nerka, and 1176 C. clupeaformis putative transcripts. We annotate the submitted portion of our EST database by molecular function. Higher- and lower-molecular-weight fractions of libraries are shown to contain distinct gene sets, and higher rates of gene discovery are associated with higher-molecular weight libraries. Pyloric caecum library group annotations indicate this organ may function in redox control and as a barrier against systemic uptake of xenobiotics. A microarray is described, containing 7356 salmonid elements representing 3557 different cDNAs. Analyses of cross-species hybridizations to this cDNA microarray indicate that this resource may be used for studies involving all salmonids.

Novel avian influenza H7N3 strain outbreak, British Columbia.

Genome sequences of chicken (low pathogenic avian influenza [LPAI] and highly pathogenic avian influenza [HPAI]) and human isolates from a 2004 outbreak of H7N3 avian influenza in Canada showed a novel insertion in the HA0 cleavage site of the human and HPAI isolate. This insertion likely occurred by recombination between the hemagglutination and matrix genes in the LPAI virus.

The Genome sequence of the SARS-associated coronavirus.

We sequenced the 29,751-base genome of the severe acute respiratory syndrome (SARS)-associated coronavirus known as the Tor2 isolate. The genome sequence reveals that this coronavirus is only moderately related to other known coronaviruses, including two human coronaviruses, HCoV-OC43 and HCoV-229E. Phylogenetic analysis of the predicted viral proteins indicates that the virus does not closely resemble any of the three previously known groups of coronaviruses. The genome sequence will aid in the diagnosis of SARS virus infection in humans and potential animal hosts (using polymerase chain reaction and immunological tests), in the development of antivirals (including neutralizing antibodies), and in the identification of putative epitopes for vaccine development.

An efficient strategy for large-scale high-throughput transposon-mediated sequencing of cDNA clones.

We describe an efficient high-throughput method for accurate DNA sequencing of entire cDNA clones. Developed as part of our involvement in the Mammalian Gene Collection full-length cDNA sequencing initiative, the method has been used and refined in our laboratory since September 2000. Amenable to large scale projects, we have used the method to generate >7 Mb of accurate sequence from 3695 candidate full-length cDNAs. Sequencing is accomplished through the insertion of Mu transposon into cDNAs, followed by sequencing reactions primed with Mu-specific sequencing primers. Transposon insertion reactions are not performed with individual cDNAs but rather on pools of up to 96 clones. This pooling strategy reduces the number of transposon insertion sequencing libraries that would otherwise be required, reducing the costs and enhancing the efficiency of the transposon library construction procedure. Sequences generated using transposon-specific sequencing primers are assembled to yield the full-length cDNA sequence, with sequence editing and other sequence finishing activities performed as required to resolve sequence ambiguities. Although analysis of the many thousands (22 785) of sequenced Mu transposon insertion events revealed a weak sequence preference for Mu insertion, we observed insertion of the Mu transposon into 1015 of the possible 1024 5mer candidate insertion sites.