Expert Panel Curation of 113 Primary Mitochondrial Disease Genes for the Leigh Syndrome Spectrum.

OBJECTIVE: Primary mitochondrial diseases (PMDs) are heterogeneous disorders caused by inherited mitochondrial dysfunction. Classically defined neuropathologically as subacute necrotizing encephalomyelopathy, Leigh syndrome spectrum (LSS) is the most frequent manifestation of PMD in children, but may also present in adults. A major challenge for accurate diagnosis of LSS in the genomic medicine era is establishing gene-disease relationships (GDRs) for this syndrome with >100 monogenic causes across both nuclear and mitochondrial genomes. METHODS: The Clinical Genome Resource (ClinGen) Mitochondrial Disease Gene Curation Expert Panel (GCEP), comprising 40 international PMD experts, met monthly for 4 years to review GDRs for LSS. The GCEP standardized gene curation for LSS by refining the phenotypic definition, modifying the ClinGen Gene-Disease Clinical Validity Curation Framework to improve interpretation for LSS, and establishing a scoring rubric for LSS. RESULTS: The GDR with LSS across the nuclear and mitochondrial genomes was classified as definitive for 31 of 114 GDRs curated (27%), moderate for 38 (33%), limited for 43 (38%), and disputed for 2 (2%). Ninety genes were associated with autosomal recessive inheritance, 16 were maternally inherited, 5 were autosomal dominant, and 3 were X-linked. INTERPRETATION: GDRs for LSS were established for genes across both nuclear and mitochondrial genomes. Establishing these GDRs will allow accurate variant interpretation, expedite genetic diagnosis of LSS, and facilitate precision medicine, multisystem organ surveillance, recurrence risk counseling, reproductive choice, natural history studies, and determination of eligibility for interventional clinical trials. ANN NEUROL 2023;94:696-712.

Specifications of the ACMG/AMP standards and guidelines for mitochondrial DNA variant interpretation.

Mitochondrial DNA (mtDNA) variant pathogenicity interpretation has special considerations given unique features of the mtDNA genome, including maternal inheritance, variant heteroplasmy, threshold effect, absence of splicing, and contextual effects of haplogroups. Currently, there are insufficient standardized criteria for mtDNA variant assessment, which leads to inconsistencies in clinical variant pathogenicity reporting. An international working group of mtDNA experts was assembled within the Mitochondrial Disease Sequence Data Resource Consortium and obtained Expert Panel status from ClinGen. This group reviewed the 2015 American College of Medical Genetics and Association of Molecular Pathology standards and guidelines that are widely used for clinical interpretation of DNA sequence variants and provided further specifications for additional and specific guidance related to mtDNA variant classification. These Expert Panel consensus specifications allow for consistent consideration of the unique aspects of the mtDNA genome that directly influence variant assessment, including addressing mtDNA genome composition and structure, haplogroups and phylogeny, maternal inheritance, heteroplasmy, and functional analyses unique to mtDNA, as well as specifications for utilization of mtDNA genomic databases and computational algorithms.

HmtDB 2016: data update, a better performing query system and human mitochondrial DNA haplogroup predictor.

The HmtDB resource hosts a database of human mitochondrial genome sequences from individuals with healthy and disease phenotypes. The database is intended to support both population geneticists as well as clinicians undertaking the task to assess the pathogenicity of specific mtDNA mutations. The wide application of next-generation sequencing (NGS) has provided an enormous volume of high-resolution data at a low price, increasing the availability of human mitochondrial sequencing data, which called for a cogent and significant expansion of HmtDB data content that has more than tripled in the current release. We here describe additional novel features, including: (i) a complete, user-friendly restyling of the web interface, (ii) links to the command-line stand-alone and web versions of the MToolBox package, an up-to-date tool to reconstruct and analyze human mitochondrial DNA from NGS data and (iii) the implementation of the Reconstructed Sapiens Reference Sequence (RSRS) as mitochondrial reference sequence. The overall update renders HmtDB an even more handy and useful resource as it enables a more rapid data access, processing and analysis. HmtDB is accessible at http://www.hmtdb.uniba.it/.

MSeqDR mvTool: A mitochondrial DNA Web and API resource for comprehensive variant annotation, universal nomenclature collation, and reference genome conversion.

Accurate mitochondrial DNA (mtDNA) variant annotation is essential for the clinical diagnosis of diverse human diseases. Substantial challenges to this process include the inconsistency in mtDNA nomenclatures, the existence of multiple reference genomes, and a lack of reference population frequency data. Clinicians need a simple bioinformatics tool that is user-friendly, and bioinformaticians need a powerful informatics resource for programmatic usage. Here, we report the development and functionality of the MSeqDR mtDNA Variant Tool set (mvTool), a one-stop mtDNA variant annotation and analysis Web service. mvTool is built upon the MSeqDR infrastructure (https://mseqdr.org), with contributions of expert curated data from MITOMAP (https://www.mitomap.org) and HmtDB (https://www.hmtdb.uniba.it/hmdb). mvTool supports all mtDNA nomenclatures, converts variants to standard rCRS- and HGVS-based nomenclatures, and annotates novel mtDNA variants. Besides generic annotations from dbNSFP and Variant Effect Predictor (VEP), mvTool provides allele frequencies in more than 47,000 germline mitogenomes, and disease and pathogenicity classifications from MSeqDR, Mitomap, HmtDB and ClinVar (Landrum et al., 2013). mvTools also provides mtDNA somatic variants annotations. "mvTool API" is implemented for programmatic access using inputs in VCF, HGVS, or classical mtDNA variant nomenclatures. The results are reported as hyperlinked html tables, JSON, Excel, and VCF formats. MSeqDR mvTool is freely accessible at https://mseqdr.org/mvtool.php.

Inherited germline ATRX mutation in two brothers with ATR-X syndrome and osteosarcoma.

We report a family in which two brothers had an undiagnosed genetic disorder comprised of dysmorphic features, microcephaly, severe intellectual disability (non-verbal), mild anemia, and cryptorchidism. Both developed osteosarcoma. Trio exome sequencing (using blood samples from the younger brother and both parents) was performed and a nonsense NM_000489.4:c.7156C>T (p.Arg2386*) mutation in the ATRX gene was identified in the proband (hemizygous) and in the mother's peripheral blood DNA (heterozygous). The mother is healthy, does not exhibit any clinical manifestations of ATR-X syndrome and there was no family history of cancer. The same hemizygous pathogenic variant was confirmed in the affected older brother's skin tissue by subsequent Sanger sequencing. Chromosomal microarray studies of both brothers' osteosarcomas revealed complex copy number alterations consistent with the clinical diagnosis of osteosarcoma. Recently, somatic mutations in the ATRX gene have been observed as recurrent alterations in both osteosarcoma and brain tumors. However, it is unclear if there is any association between osteosarcoma and germline ATRX mutations, specifically in patients with constitutional ATR-X syndrome. This is the first report of osteosarcoma diagnosed in two males with ATR-X syndrome, suggesting a potential increased risk for cancer in patients with this disorder.

MSeqDR: A Centralized Knowledge Repository and Bioinformatics Web Resource to Facilitate Genomic Investigations in Mitochondrial Disease.

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.

Phy-Mer: a novel alignment-free and reference-independent mitochondrial haplogroup classifier.

MOTIVATION: All current mitochondrial haplogroup classification tools require variants to be detected from an alignment with the reference sequence and to be properly named according to the canonical nomenclature standards for describing mitochondrial variants, before they can be compared with the haplogroup determining polymorphisms. With the emergence of high-throughput sequencing technologies and hence greater availability of mitochondrial genome sequences, there is a strong need for an automated haplogroup classification tool that is alignment-free and agnostic to reference sequence. RESULTS: We have developed a novel mitochondrial genome haplogroup-defining algorithm using a k-mer approach namely Phy-Mer. Phy-Mer performs equally well as the leading haplogroup classifier, HaploGrep, while avoiding the errors that may occur when preparing variants to required formats and notations. We have further expanded Phy-Mer functionality such that next-generation sequencing data can be used directly as input. AVAILABILITY AND IMPLEMENTATION: Phy-Mer is publicly available under the GNU Affero General Public License v3.0 on GitHub (https://github.com/danielnavarrogomez/phy-mer). CONTACT: Xiaowu_Gai@meei.harvard.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Bone marrow transcriptome and epigenome profiles of equine common variable immunodeficiency patients unveil block of B lymphocyte differentiation.

Common variable immunodeficiency (CVID) is a late-onset humoral deficiency characterized by B lymphocyte dysfunction or loss, decreased immunoglobulin production, and recurrent bacterial infections. CVID is the most frequent human primary immunodeficiency but still presents challenges in the understanding of its etiology and treatment. CVID in equine patients manifests with a natural impairment of B lymphocyte differentiation, and is a unique model to identify genetic and epigenetic mechanisms of disease. Bone marrow transcriptome analyses revealed decreased expression of genes indicative of the pro-B cell differentiation stage, importantly PAX5 (p≤0.023). We hypothesized that aberrant epigenetic regulation caused PAX5 gene silencing, resulting in the late-onset and non-familial manifestation of CVID. A significant increase in PAX5 enhancer region methylation was identified in equine CVID patients by genome-wide reduced-representation bisulfite sequencing and bisulfite PCR sequencing (p=0.000). Thus, we demonstrate that integrating transcriptomics and epigenetics in CVID enlightens potential mechanisms of dysfunctional B lymphopoiesis or function.

Mitochondrial Disease Sequence Data Resource (MSeqDR): a global grass-roots consortium to facilitate deposition, curation, annotation, and integrated analysis of genomic data for the mitochondrial disease clinical and research communities.

Success rates for genomic analyses of highly heterogeneous disorders can be greatly improved if a large cohort of patient data is assembled to enhance collective capabilities for accurate sequence variant annotation, analysis, and interpretation. Indeed, molecular diagnostics requires the establishment of robust data resources to enable data sharing that informs accurate understanding of genes, variants, and phenotypes. The "Mitochondrial Disease Sequence Data Resource (MSeqDR) Consortium" is a grass-roots effort facilitated by the United Mitochondrial Disease Foundation to identify and prioritize specific genomic data analysis needs of the global mitochondrial disease clinical and research community. A central Web portal (https://mseqdr.org) facilitates the coherent compilation, organization, annotation, and analysis of sequence data from both nuclear and mitochondrial genomes of individuals and families with suspected mitochondrial disease. This Web portal provides users with a flexible and expandable suite of resources to enable variant-, gene-, and exome-level sequence analysis in a secure, Web-based, and user-friendly fashion. Users can also elect to share data with other MSeqDR Consortium members, or even the general public, either by custom annotation tracks or through the use of a convenient distributed annotation system (DAS) mechanism. A range of data visualization and analysis tools are provided to facilitate user interrogation and understanding of genomic, and ultimately phenotypic, data of relevance to mitochondrial biology and disease. Currently available tools for nuclear and mitochondrial gene analyses include an MSeqDR GBrowse instance that hosts optimized mitochondrial disease and mitochondrial DNA (mtDNA) specific annotation tracks, as well as an MSeqDR locus-specific database (LSDB) that curates variant data on more than 1300 genes that have been implicated in mitochondrial disease and/or encode mitochondria-localized proteins. MSeqDR is integrated with a diverse array of mtDNA data analysis tools that are both freestanding and incorporated into an online exome-level dataset curation and analysis resource (GEM.app) that is being optimized to support needs of the MSeqDR community. In addition, MSeqDR supports mitochondrial disease phenotyping and ontology tools, and provides variant pathogenicity assessment features that enable community review, feedback, and integration with the public ClinVar variant annotation resource. A centralized Web-based informed consent process is being developed, with implementation of a Global Unique Identifier (GUID) system to integrate data deposited on a given individual from different sources. Community-based data deposition into MSeqDR has already begun. Future efforts will enhance capabilities to incorporate phenotypic data that enhance genomic data analyses. MSeqDR will fill the existing void in bioinformatics tools and centralized knowledge that are necessary to enable efficient nuclear and mtDNA genomic data interpretation by a range of shareholders across both clinical diagnostic and research settings. Ultimately, MSeqDR is focused on empowering the global mitochondrial disease community to better define and explore mitochondrial diseases.

A-MYB (MYBL1) transcription factor is a master regulator of male meiosis.

The transcriptional regulation of mammalian meiosis is poorly characterized, owing to few genetic and ex vivo models. From a genetic screen, we identify the transcription factor MYBL1 as a male-specific master regulator of several crucial meiotic processes. Spermatocytes bearing a novel separation-of-function allele (Mybl1(repro9)) had subtle defects in autosome synapsis in pachynema, a high incidence of unsynapsed sex chromosomes, incomplete double-strand break repair on synapsed pachytene chromosomes and a lack of crossing over. MYBL1 protein appears in pachynema, and its mutation caused specific alterations in expression of diverse genes, including some translated postmeiotically. These data, coupled with chromatin immunoprecipitation (ChIP-chip) experiments and bioinformatic analysis of promoters, identified direct targets of MYBL1 regulation. The results reveal that MYBL1 is a master regulator of meiotic genes that are involved in multiple processes in spermatocytes, particularly those required for cell cycle progression through pachynema.

MSeqDR Quick-Mitome (QM): Combining Phenotype-Guided Variant Interpretation and Machine Learning Classifiers to Aid Primary Mitochondrial Disease Genetic Diagnosis.

The international Mitochondrial Disease Sequence Data Resource Consortium (MSeqDR) Quick-Mitome (QM) is a web-based platform enabling automated variant interpretation of whole-exome sequencing (WES) datasets for the genetic diagnosis of primary mitochondrial diseases (PMD). Designed specifically to address the unique dual genome nature of PMD etiologies, QM includes features for both nuclear and mitochondrial DNA (mtDNA) genome analysis. QM requires VCF variant lists, HPO ID clinical phenotypes, and pedigree files for multiple-sample VCF inputs. QM maps phenotypes to HPO terms before analysis. QM analysis requires 2 to 20 min for 100,000 variants on an 8-vCPU AWS server using Exomiser's "PASS_ONLY" mode for nuclear variants. QM ranks variants based on allele frequency, phenotype-gene association, functional impact, and inheritance mode. Variants are further annotated with multiple data sources such as OMIM, ClinVar, dbNSFP, gnoMAD, MITOMAP, and MSeqDR. In addition to standard Exomiser results, QM generates an Analysis Report and QM Integrated Report with add-on mtDNA-specific analyses, including haplogroup prediction with Phy-Mer, heteroplasmy calculation, and mvTool annotations. We developed the Mitochondrial Disease Variant (MDV) classifier using XGBoost to predict variant pathogenicity for PMD. The MDV classifier was trained on >120 features and performance benchmarking showed that it correctly classified >98% of nuclear gene variants as being pathogenic or benign, and predicted PMD-causing variants with 94% precision. The MSeqDR QM server is an open-access resource for phenotype-driven dual-genome analyses for PMD diagnosis by the global mitochondrial disease community. It is publicly available for non-commercial, non-clinical research use at https://mseqdr.org/quickmitome.php. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Standardizing clinical phenotypes into human phenotype ontology (HPO) terms as the phenotype input for Quick-Mitome (QM) Basic Protocol 2: Prepare the pedigree input for multiple-sample VCF Basic Protocol 3: Quick-Mitome (QM) analysis Basic Protocol 4: Reviewing and understanding the QM Integrated Report and Analysis Report.

Significance Associated with Phenotype Score Aids in Variant Prioritization for Exome Sequencing Analysis.

Several in silico annotation-based methods have been developed to prioritize variants in exome sequencing analysis. This study introduced a novel metric Significance Associated with Phenotypes (SAP) score, which generates a statistical score by comparing an individual's observed phenotypes against existing gene-phenotype associations. To evaluate the SAP score, a retrospective analysis was performed on 219 exomes. Among them, 82 family-based and 35 singleton exomes had at least one disease-causing variant that explained the patient's clinical features. SAP scores were calculated, and the rank of the disease-causing variant was compared with a known method, Exomiser. Using the SAP score, the known causative variant was ranked in the top 10 retained variants for 94% (77 of 82) of the family-based exomes and in first place for 73% of these cases. For singleton exomes, the SAP score analysis ranked the known pathogenic variants within the top 10 for 80% (28 of 35) of cases. The SAP score, which is independent of detected variants, demonstrates comparable performance with Exomiser, which considers both phenotype and variant-level evidence simultaneously. Among 102 cases with negative results or variants of uncertain significance, SAP score analysis revealed two cases with a potential new diagnosis based on rank. The SAP score, a phenotypic quantitative metric, can be used in conjunction with standard variant filtration and annotation to enhance variant prioritization in exome analysis.

A genome-wide association study identifies novel risk loci for type 2 diabetes.

Type 2 diabetes mellitus results from the interaction of environmental factors with a combination of genetic variants, most of which were hitherto unknown. A systematic search for these variants was recently made possible by the development of high-density arrays that permit the genotyping of hundreds of thousands of polymorphisms. We tested 392,935 single-nucleotide polymorphisms in a French case-control cohort. Markers with the most significant difference in genotype frequencies between cases of type 2 diabetes and controls were fast-tracked for testing in a second cohort. This identified four loci containing variants that confer type 2 diabetes risk, in addition to confirming the known association with the TCF7L2 gene. These loci include a non-synonymous polymorphism in the zinc transporter SLC30A8, which is expressed exclusively in insulin-producing beta-cells, and two linkage disequilibrium blocks that contain genes potentially involved in beta-cell development or function (IDE-KIF11-HHEX and EXT2-ALX4). These associations explain a substantial portion of disease risk and constitute proof of principle for the genome-wide approach to the elucidation of complex genetic traits.

Resilience to autosomal dominant Alzheimer's disease in a Reelin-COLBOS heterozygous man.

We characterized the world's second case with ascertained extreme resilience to autosomal dominant Alzheimer's disease (ADAD). Side-by-side comparisons of this male case and the previously reported female case with ADAD homozygote for the APOE3 Christchurch (APOECh) variant allowed us to discern common features. The male remained cognitively intact until 67 years of age despite carrying a PSEN1-E280A mutation. Like the APOECh carrier, he had extremely elevated amyloid plaque burden and limited entorhinal Tau tangle burden. He did not carry the APOECh variant but was heterozygous for a rare variant in RELN (H3447R, termed COLBOS after the Colombia-Boston biomarker research study), a ligand that like apolipoprotein E binds to the VLDLr and APOEr2 receptors. RELN-COLBOS is a gain-of-function variant showing stronger ability to activate its canonical protein target Dab1 and reduce human Tau phosphorylation in a knockin mouse. A genetic variant in a case protected from ADAD suggests a role for RELN signaling in resilience to dementia.

Multidrug-resistant Pseudomonas aeruginosa in healthcare facilities in Port-au-Prince, Haiti.

OBJECTIVES: Pseudomonas aeruginosa is a leading cause of opportunistic infections worldwide, particularly in healthcare settings, and frequently demonstrates resistance to commonly prescribed antimicrobials. Carbapenem resistance is prevalent worldwide, however there are currently limited data available from Haiti. The aim of this study was to characterise and document this phenotype in Port-au-Prince, Haiti, to further inform the need for appropriate infection control, empirical treatment guidelines and laboratory screening measures, both in Haiti and globally. METHODS: A total of 50 P. aeruginosa isolates were characterised by multilocus sequence typing (MLST) and antimicrobial susceptibility testing, of which 8 isolates were also subjected to whole-genome sequencing (WGS) to identify potential genetic correlations of phenotypic resistance. RESULTS: By MLST, 23 sequence types (STs) were identified, including 13 new STs. Nineteen isolates belonged to a single, previously characterised ST (ST654), all of which demonstrated a multidrug-resistant phenotype, including resistance to meropenem, imipenem and ceftazidime; two isolates were also resistant to colistin. WGS revealed the presence of genes encoding several previously characterised resistance determinants in ST654; notably ACC(6')-Ib3-cr and GES-7. Metallo-ß-lactamase genes (blaVIM-5) were also detected in three isolates. CONCLUSION: These findings confirm that drug-resistant clones of P. aeruginosa are present in Haiti, supporting the need for appropriate screening and control measures and confirming that drug-resistant micro-organisms pose a global threat. Further investigations are required to guide appropriate antimicrobial prescribing in this region.

Rapidly emerging SARS-CoV-2 B.1.1.7 sub-lineage in the United States of America with spike protein D178H and membrane protein V70L mutations.

The SARS-CoV-2 B.1.1.7 lineage is highly infectious and as of April 2021 accounted for 92% of COVID-19 cases in Europe and 59% of COVID-19 cases in the U.S. It is defined by the N501Y mutation in the receptor-binding domain (RBD) of the Spike (S) protein, and a few other mutations. These include two mutations in the N terminal domain (NTD) of the S protein, HV69-70del and Y144del (also known as Y145del due to the presence of tyrosine at both positions). We recently identified several emerging SARS-CoV-2 variants of concerns, characterized by Membrane (M) protein mutations, including I82T and V70L. We now identify a sub-lineage of B.1.1.7 that emerged through sequential acquisitions of M:V70L in November 2020 followed by a novel S:D178H mutation first observed in early February 2021. The percentage of B.1.1.7 isolates in the US that belong to this sub-lineage increased from 0.15% in February 2021 to 1.8% in April 2021. To date, this sub-lineage appears to be U.S.-specific with reported cases in 31 states, including Hawaii. As of April 2021, it constituted 36.8% of all B.1.1.7 isolates in Washington. Phylogenetic analysis and transmission inference with Nextstrain suggest this sub-lineage likely originated in either California or Washington. Structural analysis revealed that the S:D178H mutation is in the NTD of the S protein and close to two other signature mutations of B.1.1.7, HV69-70del and Y144del. It is surface exposed and may alter NTD tertiary configuration or accessibility, and thus has the potential to affect neutralization by NTD directed antibodies.

Emerging variants of concern in SARS-CoV-2 membrane protein: a highly conserved target with potential pathological and therapeutic implications.

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.

Aqueous Humor as a Liquid Biopsy for Retinoblastoma: Clear Corneal Paracentesis and Genomic Analysis.

There is significant potential clinical utility for the application of a liquid biopsy platform for retinoblastoma, given that direct tumor biopsy is prohibited in these patients. The aqueous humor (AH) forms in a separate compartment from the tumor but is enclosed within the same ocular space. Thus, it is an enriched source of eye-specific tumoral genomic information that can be used as a liquid biopsy or surrogate to tumor biopsy for this disease. This manuscript details a methodology for safely extracting the AH from retinoblastoma eyes via clear corneal paracentesis. Additionally, the steps for genomic analysis, including cell-free DNA isolation and purification, next-generation sequencing, somatic copy number alteration (SCNA) analysis, RB1 single nucleotide variant (SNV) mutation identification, and tumor fraction estimation are presented. The pre-analytical, analytical, and early clinical validity of the AH liquid biopsy platform have been evaluated; however, it is not without limitations. These are largely a consequence of the quantity of cell-free DNA that is required for certain steps of the assay. Compared to other blood-based liquid biopsy platforms currently under investigation for retinoblastoma, an AH-based platform is limited by the volume of biofluid (and thus the quantity of DNA) that can be extracted from the eye; the benefit is that AH is eye-specific. The platform discussed here is unique in that it detects circulating tumor DNA in the AH via two mechanisms (SCNAs and RB1 SNVs), yielding a higher sensitivity for identifying tumoral genomic information. The AH liquid biopsy has the potential for direct clinical application to precision oncology for retinoblastoma patients, with particular importance for patients with bilateral disease as the AH is specific to the tumors in each eye. There is ongoing research with applications of this platform to patients with other ocular tumors as well.

High Prevalence of SARS-CoV-2 Genetic Variation and D614G Mutation in Pediatric Patients With COVID-19.

BACKGROUND: The full spectrum of the disease phenotype and viral genotype of coronavirus disease 2019 (COVID-19) have yet to be thoroughly explored in children. Here, we analyze the relationships between viral genetic variants and clinical characteristics in children. METHODS: Whole-genome sequencing was performed on respiratory specimens collected for all SARS-CoV-2-positive children (n = 141) between March 13 and June 16, 2020. Viral genetic variations across the SARS-CoV-2 genome were identified and investigated to evaluate genomic correlates of disease severity. RESULTS: Higher viral load was detected in symptomatic patients (P = .0007) and in children <5 years old (P = .0004). Genomic analysis revealed a mean pairwise difference of 10.8 single nucleotide variants (SNVs), and the majority (55.4%) of SNVs led to an amino acid change in the viral proteins. The D614G mutation in the spike protein was present in 99.3% of the isolates. The calculated viral mutational rate of 22.2 substitutions/year contrasts the 13.5 substitutions/year observed in California isolates without the D614G mutation. Phylogenetic clade 20C was associated with severe cases of COVID-19 (odds ratio, 6.95; P = .0467). Epidemiological investigation revealed major representation of 3 of 5 major Nextstrain clades (20A, 20B, and 20C) consistent with multiple introductions of SARS-CoV-2 in Southern California. CONCLUSIONS: Genomic evaluation demonstrated greater than expected genetic diversity, presence of the D614G mutation, increased mutation rate, and evidence of multiple introductions of SARS-CoV-2 into Southern California. Our findings suggest a possible association of phylogenetic clade 20C with severe disease, but small sample size precludes a definitive conclusion. Our study warrants larger and multi-institutional genomic evaluation and has implications for infection control practices.

Inter-eye genomic heterogeneity in bilateral retinoblastoma via aqueous humor liquid biopsy.

Germline alterations in the RB1 tumor suppressor gene predispose patients to develop retinoblastoma (RB) in both eyes. While similar treatment is given for each eye, there is often a variable therapeutic response between the eyes. Herein, we use the aqueous humor (AH) liquid biopsy to evaluate the cell-free tumor DNA (ctDNA) from each eye in a patient with bilateral RB. Despite the same predisposing germline RB1 mutation, AH analysis identified a different somatic RB1 mutation as well as separate and distinct chromosomal alterations in each eye. The longitudinal alterations in tumor fraction (TFx) corresponded to therapeutic responses in each eye. This case demonstrates that bilateral RB tumors develop separate genomic alterations, which may play a role in tumorigenesis and prognosis for eye salvage. Identifying these inter-eye differences without the need for enucleated tumor tissue may help direct active management of RB, with particular usefulness in bilateral cases.