Chromosomal microarray analysis of 410 Han Chinese patients with autism spectrum disorder or unexplained intellectual disability and developmental delay

Copy number variants (CNVs) are recognized as a crucial genetic cause of neurodevelopmental disorders (NDDs). Chromosomal microarray analysis (CMA), the first-tier diagnostic test for individuals with NDDs, has been utilized to detect CNVs in clinical practice, but most reports are still from populations of European ancestry. To contribute more worldwide clinical genomics data, we investigated the genetic etiology of 410 Han Chinese patients with NDDs (151 with autism and 259 with unexplained intellectual disability (ID) and developmental delay (DD)) using CMA (Affymetrix) after G-banding karyotyping. Among all the NDD patients, 109 (26.6%) carried clinically relevant CNVs or uniparental disomies (UPDs), and 8 (2.0%) had aneuploidies (6 with trisomy 21 syndrome, 1 with 47,XXY, 1 with 47,XYY). In total, we found 129 clinically relevant CNVs and UPDs, including 32 CNVs in 30 ASD patients, and 92 CNVs and 5 UPDs in 79 ID/DD cases. When excluding the eight patients with aneuploidies, the diagnostic yield of pathogenic and likely pathogenic CNVs and UPDs was 20.9% for all NDDs (84/402), 3.3% in ASD (5/151), and 31.5% in ID/DD (79/251). When aneuploidies were included, the diagnostic yield increased to 22.4% for all NDDs (92/410), and 33.6% for ID/DD (87/259). We identified a de novo CNV in 14.9% (60/402) of subjects with NDDs. Interestingly, a higher diagnostic yield was observed in females (31.3%, 40/128) compared to males (16.1%, 44/274) for all NDDs (P = 4.8 × 10−4), suggesting that a female protective mechanism exists for deleterious CNVs and UPDs.

Most myopathic lamin variants aggregate: a functional genomics approach for assessing variants of uncertain significance

Hundreds of LMNA variants have been associated with several distinct disease phenotypes. However, genotype–phenotype relationships remain largely undefined and the impact for most variants remains unknown. We performed a functional analysis for 178 variants across five structural domains using two different overexpression models. We found that lamin A aggregation is a major determinant for skeletal and cardiac laminopathies. An in vitro solubility assay shows that aggregation-prone variants in the immunoglobulin-like domain correlate with domain destabilization. Finally, we demonstrate that myopathic-associated LMNA variants show aggregation patterns in induced pluripotent stem cell derived-cardiomyocytes (iPSC-CMs) in contrast to non-myopathic LMNA variants. Our data-driven approach (1) reveals that striated muscle laminopathies are predominantly protein misfolding diseases, (2) demonstrates an iPSC-CM experimental platform for characterizing laminopathic variants in human cardiomyocytes, and (3) supports a functional assay to aid in assessing pathogenicity for myopathic variants of uncertain significance.

The salivary metatranscriptome as an accurate diagnostic indicator of oral cancer

Despite advances in cancer treatment, the 5-year mortality rate for oral cancers (OC) is 40%, mainly due to the lack of early diagnostics. To advance early diagnostics for high-risk and average-risk populations, we developed and evaluated machine-learning (ML) classifiers using metatranscriptomic data from saliva samples (n = 433) collected from oral premalignant disorders (OPMD), OC patients (n = 71) and normal controls (n = 171). Our diagnostic classifiers yielded a receiver operating characteristics (ROC) area under the curve (AUC) up to 0.9, sensitivity up to 83% (92.3% for stage 1 cancer) and specificity up to 97.9%. Our metatranscriptomic signature incorporates both taxonomic and functional microbiome features, and reveals a number of taxa and functional pathways associated with OC. We demonstrate the potential clinical utility of an AI/ML model for diagnosing OC early, opening a new era of non-invasive diagnostics, enabling early intervention and improved patient outcomes.

De novo variants in H3-3A and H3-3B are associated with neurodevelopmental delay, dysmorphic features, and structural brain abnormalities

The histone H3 variant H3.3, encoded by two genes H3-3A and H3-3B, can replace canonical isoforms H3.1 and H3.2. H3.3 is important in chromatin compaction, early embryonic development, and lineage commitment. The role of H3.3 in somatic cancers has been studied extensively, but its association with a congenital disorder has emerged just recently. Here we report eleven de novo missense variants and one de novo stop-loss variant in H3-3A (n = 6) and H3-3B (n = 6) from Baylor Genetics exome cohort (n = 11) and Matchmaker Exchange (n = 1), of which detailed phenotyping was conducted for 10 individuals (H3-3A = 4 and H3-3B = 6) that showed major phenotypes including global developmental delay, short stature, failure to thrive, dysmorphic facial features, structural brain abnormalities, hypotonia, and visual impairment. Three variant constructs (p.R129H, p.M121I, and p.I52N) showed significant decrease in protein expression, while one variant (p.R41C) accumulated at greater levels than wild-type control. One H3.3 variant construct (p.R129H) was found to have stronger interaction with the chaperone death domain-associated protein 6.

Accurate detection of circulating tumor DNA using nanopore consensus sequencing

Levels of circulating tumor DNA (ctDNA) in liquid biopsies may serve as a sensitive biomarker for real-time, minimally-invasive tumor diagnostics and monitoring. However, detecting ctDNA is challenging, as much fewer than 5% of the cell-free DNA in the blood typically originates from the tumor. To detect lowly abundant ctDNA molecules based on somatic variants, extremely sensitive sequencing methods are required. Here, we describe a new technique, CyclomicsSeq, which is based on Oxford Nanopore sequencing of concatenated copies of a single DNA molecule. Consensus calling of the DNA copies increased the base-calling accuracy ~60×, enabling accurate detection of TP53 mutations at frequencies down to 0.02%. We demonstrate that a TP53-specific CyclomicsSeq assay can be successfully used to monitor tumor burden during treatment for head-and-neck cancer patients. CyclomicsSeq can be applied to any genomic locus and offers an accurate diagnostic liquid biopsy approach that can be implemented in clinical workflows.

Whole exome sequencing identifies deleterious rare variants in CCDC141 in familial self-limited delayed puberty

Developmental abnormalities of the gonadotropin-releasing hormone (GnRH) neuronal network result in a range of conditions from idiopathic hypogonadotropic hypogonadism to self-limited delayed puberty. We aimed to discover important underlying regulators of self-limited delayed puberty through interrogation of GnRH pathways. Whole exome sequencing (WES) data consisting of 193 individuals, from 100 families with self-limited delayed puberty, was analysed using a virtual panel of genes related to GnRH development and function (n = 12). Five rare predicted deleterious variants in Coiled-Coil Domain Containing 141 (CCDC141) were identified in 21 individuals from 6 families (6% of the tested cohort). Homology modeling predicted all five variants to be deleterious. CCDC141 mutant proteins showed atypical subcellular localization associated with abnormal distribution of acetylated tubulin, and expression of mutants resulted in a significantly delayed cell migration, demonstrated in transfected HEK293 cells. These data identify mutations in CCDC141 as a frequent finding in patients with self-limited delayed puberty. The mis-localization of acetylated tubulin and reduced cell migration seen with mutant CCDC141 suggests a role of the CCDC141-microtubule axis in GnRH neuronal migration, with heterozygous defects potentially impacting the timing of puberty.

Recurrent integration of human papillomavirus genomes at transcriptional regulatory hubs

Oncogenic human papillomavirus (HPV) genomes are often integrated into host chromosomes in HPV-associated cancers. HPV genomes are integrated either as a single copy or as tandem repeats of viral DNA interspersed with, or without, host DNA. Integration occurs frequently in common fragile sites susceptible to tandem repeat formation and the flanking or interspersed host DNA often contains transcriptional enhancer elements. When co-amplified with the viral genome, these enhancers can form super-enhancer-like elements that drive high viral oncogene expression. Here we compiled highly curated datasets of HPV integration sites in cervical (CESC) and head and neck squamous cell carcinoma (HNSCC) cancers, and assessed the number of breakpoints, viral transcriptional activity, and host genome copy number at each insertion site. Tumors frequently contained multiple distinct HPV integration sites but often only one “driver” site that expressed viral RNA. As common fragile sites and active enhancer elements are cell-type-specific, we mapped these regions in cervical cell lines using FANCD2 and Brd4/H3K27ac ChIP-seq, respectively. Large enhancer clusters, or super-enhancers, were also defined using the Brd4/H3K27ac ChIP-seq dataset. HPV integration breakpoints were enriched at both FANCD2-associated fragile sites and enhancer-rich regions, and frequently showed adjacent focal DNA amplification in CESC samples. We identified recurrent integration “hotspots” that were enriched for super-enhancers, some of which function as regulatory hubs for cell-identity genes. We propose that during persistent infection, extrachromosomal HPV minichromosomes associate with these transcriptional epicenters and accidental integration could promote viral oncogene expression and carcinogenesis.

Geometric network analysis provides prognostic information in patients with high grade serous carcinoma of the ovary treated with immune checkpoint inhibitors

Network analysis methods can potentially quantify cancer aberrations in gene networks without introducing fitted parameters or variable selection. A new network curvature-based method is introduced to provide an integrated measure of variability within cancer gene networks. The method is applied to high-grade serous ovarian cancers (HGSOCs) to predict response to immune checkpoint inhibitors (ICIs) and to rank key genes associated with prognosis. Copy number alterations (CNAs) from targeted and whole-exome sequencing data were extracted for HGSOC patients (n = 45) treated with ICIs. CNAs at a gene level were represented on a protein–protein interaction network to define patient-specific networks with a fixed topology. A version of Ollivier–Ricci curvature was used to identify genes that play a potentially key role in response to immunotherapy and further to stratify patients at high risk of mortality. Overall survival (OS) was defined as the time from the start of ICI treatment to either death or last follow-up. Kaplan–Meier analysis with log-rank test was performed to assess OS between the high and low curvature classified groups. The network curvature analysis stratified patients at high risk of mortality with p = 0.00047 in Kaplan–Meier analysis in HGSOC patients receiving ICI. Genes with high curvature were in accordance with CNAs relevant to ovarian cancer. Network curvature using CNAs has the potential to be a novel predictor for OS in HGSOC patients treated with immunotherapy.

The diagnostic trajectory of infants and children with clinical features of genetic disease

We characterized US pediatric patients with clinical indicators of genetic diseases, focusing on the burden of disease, utilization of genetic testing, and cost of care. Curated lists of diagnosis, procedure, and billing codes were used to identify patients with clinical indicators of genetic disease in healthcare claims from Optum’s de-identified Clinformatics® Database (13,076,038 unique patients). Distinct cohorts were defined to represent permissive and conservative estimates of the number of patients. Clinical phenotypes suggestive of genetic diseases were observed in up to 9.4% of pediatric patients and up to 44.7% of critically-ill infants. Compared with controls, patients with indicators of genetic diseases had higher utilization of services (e.g., mean NICU length of stay of 31.6d in a cohort defined by multiple congenital anomalies or neurological presentations compared with 10.1d for patients in the control population (P < 0.001)) and higher overall costs. Very few patients received any genetic testing (4.2–8.4% depending on cohort criteria). These results highlight the substantial proportion of the population with clinical features associated with genetic disorders and underutilization of genetic testing in these populations.