The ability to capture and sequence large contiguous DNA fragments represents a significant advancement towards the comprehensive characterization of complex genomic regions. While emerging sequencing platforms are capable of producing several kilobases-long reads, the fragment sizes generated by current DNA target enrichment technologies remain a limiting factor, producing DNA fragments generally shorter than 1 kbp. The DNA enrichment methodology described herein, Region-Specific Extraction (RSE), produces DNA segments in excess of 20 kbp in length. Coupling this enrichment method to appropriate sequencing platforms will significantly enhance the ability to generate complete and accurate sequence characterization of any genomic region without the need for reference-based assembly. RSE is a long-range DNA target capture methodology that relies on the specific hybridization of short (20-25 base) oligonucleotide primers to selected sequence motifs within the DNA target region. These capture primers are then enzymatically extended on the 3’-end, incorporating biotinylated nucleotides into the DNA. Streptavidin-coated beads are subsequently used to pull-down the original, long DNA template molecules via the newly synthesized, biotinylated DNA that is bound to them. We demonstrate the accuracy, simplicity and utility of the RSE method by capturing and sequencing a 4 Mbp stretch of the major histocompatibility complex (MHC). Our results show an average depth of coverage of 164X for the entire MHC. This depth of coverage contributes significantly to a 99.94 % total coverage of the targeted region and to an accuracy that is over 99.99 %. RSE represents a cost-effective target enrichment method capable of producing sequencing templates in excess of 20 kbp in length. The utility of our method has been proven to generate superior coverage across the MHC as compared to other commercially available methodologies, with the added advantage of producing longer sequencing templates amenable to DNA sequencing on recently developed platforms. Although our demonstration of the method does not utilize these DNA sequencing platforms directly, our results indicate that the capture of long DNA fragments produce superior coverage of the targeted region.

Because direct tumor biopsy is prohibited for retinoblastoma (RB), eye-specific molecular biomarkers are not used in clinical practice for RB. Recently, we demonstrated that the aqueous humor (AH) is a rich liquid biopsy source of cell-free tumor DNA. Herein, we detail clinically-relevant molecular biomarkers from the first year of prospective validation data. Seven eyes from 6 RB patients who had AH sampled at diagnosis and throughout therapy with ≥12 months of follow-up were included. Cell-free DNA (cfDNA) from each sample was isolated and sequenced to assess genome-wide somatic copy number alterations (SCNAs), followed by targeted resequencing for pathogenic variants using a RB1 and MYCN custom hybridization panel. Tumoral genomic information was detected in 100% of diagnostic AH samples. Of the seven diagnostic AH samples, 5/7 were positive for RB SCNAs. Mutational analysis identified RB1 variants in 5/7 AH samples, including the 2 samples in which no SCNAs were detected. Two eyes failed therapy and required enucleation; both had poor prognostic biomarkers (chromosome 6p gain or MYCN amplification) present in the AH at the time of diagnosis. In the context of previously established pre-analytical, analytical, and clinical validity, this provides evidence for larger, prospective studies to further establish the clinical utility of the AH liquid biopsy and its applications to precision oncology for RB.

Advancing the clinical utility of whole-exome sequencing (WES) for patients with suspected genetic disorders is largely driven by bioinformatics approaches that streamline data processing and analysis. Herein, we describe our experience with implementing a semiautomated and phenotype-driven WES diagnostic workflow, incorporating both the DRAGEN pipeline and the Exomiser variant prioritization tool, at an academic children's hospital with an ethnically diverse pediatric patient population. We achieved a 41% molecular diagnostic rate for 66 duo-, quad-, or trio-WES cases, and 28% for 40 singleton-WES cases. Preliminary results were returned to ordering physicians within 1 wk for 12 of 38 (32%) probands with positive findings, which were instrumental in guiding the appropriate clinical management for a variety of patients, especially in critical care settings. The semiautomated and streamlined WES workflow also enabled us to identify novel variants in candidate disease genes in patients with developmental delay and autism and immune disorders and cancer, including ANK2, BPTF, BCL11A, FOXN1, PLAA, ATRX, DNAJC21, and RAD50. Together, we demonstrated the implementation of a streamlined WES workflow that was successfully applied for both clinical and research purposes.

Mutations in the SARS-CoV-2 Membrane (M) gene are relatively uncommon. The M gene encodes the most abundant viral structural protein, and is implicated in multiple viral functions, including initial attachment to the host cell via heparin sulphate proteoglycan, viral protein assembly in conjunction with the N and E genes, and enhanced glucose transport. We have identified a recent spike in the frequency of reported SARS-CoV-2 genomes carrying M gene mutations. This is associated with emergence of a new sub-B.1 clade, B.1.I82T, defined by the previously unreported M:I82T mutation within TM3, the third of three membrane spanning helices implicated in glucose transport. The frequency of this mutation increased in the USA from 0.014% in October 2020 to 1.62% in February 2021, a 116-fold change. While constituting 0.7% of the isolates overall, M:I82T sub-B.1 lineage accounted for 14.4% of B.1 lineage isolates in February 2021, similar to the rapid initial increase previously seen with the B.1.1.7 and B.1.429 lineages, which quickly became the dominant lineages in Europe and California over a period of several months. A similar increase in incidence was also noted in another related mutation, V70L, also within the TM2 transmembrane helix. These M mutations are associated with younger patient age (4.6 to 6.3 years). The rapid emergence of this B.1.I82T clade, recently named Pangolin B.1.575 lineage, suggests that this M gene mutation is more biologically fit, perhaps related to glucose uptake during viral replication, and should be included in ongoing genomic surveillance efforts and warrants further evaluation for potentially increased pathogenic and therapeutic implications.

MSeqDR is the Mitochondrial Disease Sequence Data Resource, a centralized and comprehensive genome and phenome bioinformatics resource built by the mitochondrial disease community to facilitate clinical diagnosis and research investigations of individual patient phenotypes, genomes, genes, and variants. A central Web portal (https://mseqdr.org) integrates community knowledge from expert‐curated databases with genomic and phenotype data shared by clinicians and researchers. MSeqDR also functions as a centralized application server for Web‐based tools to analyze data across both mitochondrial and nuclear DNA, including investigator‐driven whole exome or genome dataset analyses through MSeqDR‐Genesis. MSeqDR‐GBrowse genome browser supports interactive genomic data exploration and visualization with custom tracks relevant to mtDNA variation and mitochondrial disease. MSeqDR‐LSDB is a locus‐specific database that currently manages 178 mitochondrial diseases, 1,363 genes associated with mitochondrial biology or disease, and 3,711 pathogenic variants in those genes. MSeqDR Disease Portal allows hierarchical tree‐style disease exploration to evaluate their unique descriptions, phenotypes, and causative variants. Automated genomic data submission tools are provided that capture ClinVar compliant variant annotations. PhenoTips will be used for phenotypic data submission on deidentified patients using human phenotype ontology terminology. The development of a dynamic informed patient consent process to guide data access is underway to realize the full potential of these resources.