CNV analysis in Chinese children of mental retardation highlights a sex differentiation in parental contribution to de novo and inherited mutational burdens

Detection and characterization of copy-number variants from exome sequencing in the DDD study

Purpose

Structural variants such as multiexon deletions and duplications are an important cause of disease but are often overlooked in standard exome/genome sequencing analysis. We aimed to evaluate the detection of copy-number variants (CNVs) from exome sequencing (ES) in comparison with genome-wide low-resolution and exon-resolution chromosomal microarrays (CMAs) and to characterize the properties of de novo CNVs in a large clinical cohort.

Methods

We performed CNV detection using ES of 9859 parent-offspring trios in the Deciphering Developmental Disorders (DDD) study and compared them with CNVs detected from exon-resolution array comparative genomic hybridization in 5197 probands from the DDD study.

Results

Integrating calls from multiple ES-based CNV algorithms using random forest machine learning generated a higher quality data set than using individual algorithms. Both ES- and array comparative genomic hybridization-based approaches had the same sensitivity of 89% and detected the same number of unique pathogenic CNVs not called by the other approach. Of DDD probands prescreened with low-resolution CMAs, 2.6% had a pathogenic CNV detected by higher-resolution assays. De novo CNVs were strongly enriched in known DD-associated genes and exhibited no bias in parental age or sex.

Conclusion

ES-based CNV calling has higher sensitivity than low-resolution CMAs currently in clinical use and comparable sensitivity to exon-resolution CMA. With sufficient investment in bioinformatic analysis, exome-based CNV detection could replace low-resolution CMA for detecting pathogenic CNVs.

Human phenotype caused by biallelic KDM4B frameshift variant

KDM4B (MIM*609765, NM_015015.3, formerly JMJD2B) encodes a histone demethylase and regulates gene expression via demethylation, mainly of H3K9 tri-methylation. Heterozygous KDM4B loss-of-function variants cause autosomal dominant intellectual developmental disorder 65 (MIM#619320), which is characterized by global developmental delay, intellectual disability, language and gross motor delays, structural brain anomalies, characteristic facial features, and clinodactyly. Although the majority of reported patients have de novo pathogenic variants, some patients inherit pathogenic variants from affected parents. To our knowledge, only 23 patients with heterozygous KDM4B variants have been reported to date, and there are no reports of patients with biallelic KDM4B pathogenic variants. Herein, we report a female patient with a biallelic KDM4B frameshift variant (NM_015015.3: c.1384_1394delinsGGG, p.(Leu462Glyfs*43)) located at exon 12 of 23 protein-coding exons, which is thought to be subject to nonsense-mediated mRNA decay and no protein production. She presented developmental and language delays and a hypotonic and characteristic face. The patient's phenotype was more obvious than that of her mother, who is heterozygous for the same variant. Although declining birth rate (embryonic lethality in male mice) in homozygous knockout mice has been demonstrated, our report suggests that homozygous KDM4B frameshift variants can be viable in humans at least female.

Rare copy number variants in males and females with childhood attention-deficit/hyperactivity disorder

While childhood attention-deficit/hyperactivity disorder (ADHD) is more prevalent in males than females, genetic contributors to this effect have not been established. Here, we explore sex differences in the contribution of common and/or rare genetic variants to ADHD. Participants were from the Adolescent Brain and Cognitive Development study (N = 1253 youth meeting DSM-5 criteria for ADHD [mean age = 11.46 years [SD = 0.87]; 31% female] and 5577 unaffected individuals [mean age = 11.42 years [SD = 0.89]; 50% female], overall 66% White, non-Hispanic (WNH), 19% Black/African American, and 15% other races. Logistic regression tested for interactions between sex (defined genotypically) and both rare copy number variants (CNV) and polygenic (common variant) risk in association with ADHD. There was a significant interaction between sex and the presence of a CNV deletion larger than 200 kb, both in the entire cohort (β = -0.74, CI = [-1.27 to -0.20], FDR-corrected p = 0.048) and, at nominal significance levels in the WNH ancestry subcohort (β = -0.86, CI = [-1.51 to -0.20], p = 0.010). Additionally, the number of deleted genes interacted with sex in association with ADHD (whole cohort. β = -0.13, CI = [-0.23 to -0.029], FDR-corrected p = 0.048; WNH. β = -0.17, CI = [-0.29 to -0.050], FDR-corrected p = 0.044) as did the total length of CNV deletions (whole cohort. β = -0.12, CI = [-0.19 to -0.044], FDR-corrected p = 0.028; WNH. β = -0.17, CI = [-0.28 to -0.061], FDR-corrected p = 0.034). This sex effect was driven by increased odds of childhood ADHD for females but not males in the presence of CNV deletions. No similar sex effect was found for CNV duplications or polygenic risk scores. The association between CNV deletions and ADHD was partially mediated by measures of cognitive flexibility. In summary, CNV deletions were associated with increased odds for childhood ADHD in females, but not males.

Urine Organic Acids as Metabolic Indicators for Global Developmental Delay/Intellectual Disability in Chinese Children

Objective: The purpose of this study was to search for differential metabolites in urine organic acids, and to characterize metabolic features that can be used to identify metabolites for exploration of global developmental delay (GDD)/intellectual disability (ID) etiology and pathogenesis. Methods: We screened positive test results that could explain GDD/ID from 1,253 cases, and the major differential metabolites in 132 urine organic acids from the 1,230 cases with negative results (863 GDD cases, 367 ID cases), and 100 typically developing children (TD). Non-supervisory principal component analysis and orthogonal partial least squares discriminant analysis were used to develop models to distinguish GDD/ID from TD children, and to detect major differential metabolites. Results: We get 23 positive results that could identify the cause of GDD/ID from 1253 cases diagnosed with GDD/ID. Among 1,230 negative results, we get the differential metabolites of the GDD group and the ID group had the same trend compared with the TD group. Twenty four differential metabolites were obtained from the GDD group, and 25 from the ID group (VIP > 1.0, p < 0.01). These differential metabolites were mainly related to the following pathways: the synthesis and degradation of ketone bodies, citrate cycle, alanine, aspartate and glutamate metabolism, pyrimidine metabolism, butanoate metabolism, pyruvate metabolism, fatty acid biosynthesis, valine, leucine and isoleucine degradation. Conclusion: The use of metabolomics research methods to detect urine organic acids of children with GDD/ID can discover differential metabolites, which might be valuable for future research on the etiology, pathogenesis, prognosis and possible interventions of GDD/ID. The significantly altered differential metabolites indicators could therefore be potential diagnostic biomarkers for GDD/ID.

Chromosomal variants accumulate in genomes of the spontaneous aborted fetuses revealed by chromosomal microarray analysis

Spontaneous abortion is an impeding factor for the success rates of human assistant reproductive technology (ART). Causes of spontaneous abortion include not only the pregnant mothers’ health conditions and lifestyle habits, but also the fetal development potential. Evidences had shown that fetal chromosome aneuploidy is associated with fetal spontaneous abortion, however, it is still not definite that whether other genome variants, like copy number variations (CNVs) or loss of heterozygosity (LOHs) is associated with the spontaneous abortion. To assess the relationship between the fetal genome variants and abortion during ART, a chromosomal microarray data including chromosomal information of 184 spontaneous aborted fetuses, 147 adult female patients and 78 adult male patients during ART were collected. We firstly analyzed the relationship of fetal aneuploidy with maternal ages and then compared the numbers and lengths of CNVs (< 4Mbp) and LOHs among adults and aborted fetuses. In addition to the already known association between chromosomal aneuploidy and maternal ages, from the chromosomal microarray data we found that the numbers and the accumulated lengths of short CNVs and LOHs in the aborted fetuses were significantly larger or longer than those in adults. Our findings indicated that the increased numbers and accumulated lengths of CNVs or LOHs might be associated with the spontaneous abortion during ART.

Rare and de novo variants in 827 congenital diaphragmatic hernia probands implicate LONP1 as candidate risk gene

Congenital diaphragmatic hernia (CDH) is a severe congenital anomaly that is often accompanied by other anomalies. Although the role of genetics in the pathogenesis of CDH has been established, only a small number of disease-associated genes have been identified. To further investigate the genetics of CDH, we analyzed de novo coding variants in 827 proband-parent trios and confirmed an overall significant enrichment of damaging de novo variants, especially in constrained genes. We identified LONP1 (lon peptidase 1, mitochondrial) and ALYREF (Aly/REF export factor) as candidate CDH-associated genes on the basis of de novo variants at a false discovery rate below 0.05. We also performed ultra-rare variant association analyses in 748 affected individuals and 11,220 ancestry-matched population control individuals and identified LONP1 as a risk gene contributing to CDH through both de novo and ultra-rare inherited largely heterozygous variants clustered in the core of the domains and segregating with CDH in affected familial individuals. Approximately 3% of our CDH cohort who are heterozygous with ultra-rare predicted damaging variants in LONP1 have a range of clinical phenotypes, including other anomalies in some individuals and higher mortality and requirement for extracorporeal membrane oxygenation. Mice with lung epithelium-specific deletion of Lonp1 die immediately after birth, most likely because of the observed severe reduction of lung growth, a known contributor to the high mortality in humans. Our findings of both de novo and inherited rare variants in the same gene may have implications in the design and analysis for other genetic studies of congenital anomalies.

Non-Coding Variants in Cancer: Mechanistic Insights and Clinical Potential for Personalized Medicine

The cancer genome is characterized by extensive variability, in the form of Single Nucleotide Polymorphisms (SNPs) or structural variations such as Copy Number Alterations (CNAs) across wider genomic areas. At the molecular level, most SNPs and/or CNAs reside in non-coding sequences, ultimately affecting the regulation of oncogenes and/or tumor-suppressors in a cancer-specific manner. Notably, inherited non-coding variants can predispose for cancer decades prior to disease onset. Furthermore, accumulation of additional non-coding driver mutations during progression of the disease, gives rise to genomic instability, acting as the driving force of neoplastic development and malignant evolution. Therefore, detection and characterization of such mutations can improve risk assessment for healthy carriers and expand the diagnostic and therapeutic toolbox for the patient. This review focuses on functional variants that reside in transcribed or not transcribed non-coding regions of the cancer genome and presents a collection of appropriate state-of-the-art methodologies to study them.

The diagnostic yield of intellectual disability: combined whole genome low-coverage sequencing and medical exome sequencing

Background

Intellectual disability (ID) is a heterogeneous neurodevelopmental disorder with a complex genetic underpinning in its etiology. Chromosome microarray (CMA) is recommended as the first-tier diagnostic test for ID due to high detection rate of copy number variation (CNV).

Methods

To identify an appropriate clinical detection scheme for ID in Han Chinese patients, whole genome low-coverage sequencing was performed as the first-tier diagnostic test, and medical exome sequencing (MES) as the second-tier diagnostic test for patients with negative results of CNVs.

Results

A total of 19 pathogenic CNVs in 16/95(16.84%) ID patients and 10 pathogenic single-nucleotide variations (SNVs), including 6 novel mutations in 8/95(8.42%) ID patients were identified on whom no pathogenic CNVs were discovered. The detection rate of CNVs in ID with multiple congenital anomalies (MCA) subgroup was significantly higher than ID with autism spectrum disorders and other IDs subgroups. And the single-nucleotide variations showed a higher occurrence rate in the other IDs subgroup.

Conclusions

There were differences in the diagnostic yields of different variation types among the three ID subgroups. Our findings provided a new perspective on appropriate clinical detection scheme in different ID subgroups based on statistically significant differences among the three ID subgroups. The application of whole genome low-coverage sequencing as the first-tier diagnostic test for ID with MCA subgroup and MES as the first-tier diagnostic test for other ID subgroup was considered as an efficient clinical detection scheme.

Copy number variations in candidate genomic regions confirm genetic heterogeneity and parental bias in Hirschsprung disease

Background

Hirschsprung Disease (HSCR) is a congenital defect of the intestinal innervations characterized by complex inheritance. Many susceptibility genes including RET, the major HSCR gene, and several linked regions and associated loci have been shown to contribute to disease pathogenesis. Nonetheless, a proportion of patients still remains unexplained. Copy Number Variations (CNVs) have already been involved in HSCR, and for this reason we performed Comparative Genomic Hybridization (CGH), using a custom array with high density probes.

Results

A total of 20 HSCR candidate regions/genes was tested in 55 sporadic patients and four patients with already known chromosomal aberrations. Among 83 calls, 12 variants were experimentally validated, three of which involving the HSCR crucial genes SEMA3A/3D, NRG1, and PHOX2B. Conversely RET involvement in HSCR does not seem to rely on the presence of CNVs while, interestingly, several gains and losses did co-occur with another RET defect, thus confirming that more than one predisposing event is necessary for HSCR to develop. New loci were also shown to be involved, such as ALDH1A2, already found to play a major role in the enteric nervous system. Finally, all the inherited CNVs were of maternal origin.

Conclusions

Our results confirm a wide genetic heterogeneity in HSCR occurrence and support a role of candidate genes in expression regulation and cell signaling, thus contributing to depict further the molecular complexity of the genomic regions involved in the Enteric Nervous System development. The observed maternal transmission bias for HSCR associated CNVs supports the hypothesis that in females these variants might be more tolerated, requiring additional alterations to develop HSCR disease.

Experience of chromosomal microarray applied in prenatal and postnatal settings in Hong Kong

Chromosomal microarray (CMA) is recommended as a first tier investigation for patients with developmental delay (DD), intellectual disability (ID), autistic spectrum disorder (ASD), and multiple congenital anomalies (MCA). It is widely used in the prenatal and postnatal settings for detection of chromosomal aberrations. This is a retrospective review of all array comparative genomic hybridization (aCGH/ array CGH) findings ascertained in two major prenatal and postnatal genetic diagnostic centers in Hong Kong from June 2012 to December 2017. Medical records were reviewed for cases with pathogenic and variants of uncertain clinical significance (VUS). Classification of copy number variants (CNVs) was based on current knowledge and experience by August 2018. The aims of this review are to study the diagnostic yield of array CGH application in prenatal and postnatal settings in Hong Kong and to describe the spectrum of abnormalities found. Prenatal indications included abnormal ultrasound findings, positive Down syndrome screening, abnormal noninvasive prenatal test results, advanced maternal age and family history of chromosomal or genetic abnormalities. Postnatal indications included unexplained DD, ID, ASD, and MCA. A total of 1,261 prenatal subjects and 3,096 postnatal patients were reviewed. The prenatal diagnostic yield of pathogenic CNV and VUS (excluding those detectable by karyotype) was 3.5%. The postnatal diagnostic yield of pathogenic CNV was 15.2%. The detection rates for well-defined microdeletion and microduplication syndromes were 4.6% in prenatal and 6.1% (1 in 16 index patients) in postnatal cases, respectively. Chromosomes 15, 16, and 22 accounted for over 21 and 25% of pathogenic CNVs detected in prenatal and postnatal cohorts, respectively. This review provides the first large scale overview of genomic imbalance of mostly Chinese patients in prenatal and postnatal settings.

De novo variants in congenital diaphragmatic hernia identify MYRF as a new syndrome and reveal genetic overlaps with other developmental disorders

Congenital diaphragmatic hernia (CDH) is a severe birth defect that is often accompanied by other congenital anomalies. Previous exome sequencing studies for CDH have supported a role of de novo damaging variants but did not identify any recurrently mutated genes. To investigate further the genetics of CDH, we analyzed de novo coding variants in 362 proband-parent trios including 271 new trios reported in this study. We identified four unrelated individuals with damaging de novo variants in MYRF (P = 5.3x10(-8)), including one likely gene-disrupting (LGD) and three deleterious missense (D-mis) variants. Eight additional individuals with de novo LGD or missense variants were identified from our other genetic studies or from the literature. Common phenotypes of MYRF de novo variant carriers include CDH, congenital heart disease and genitourinary abnormalities, suggesting that it represents a novel syndrome. MYRF is a membrane associated transcriptional factor highly expressed in developing diaphragm and is depleted of LGD variants in the general population. All de novo missense variants aggregated in two functional protein domains. Analyzing the transcriptome of patient-derived diaphragm fibroblast cells suggest that disease associated variants abolish the transcription factor activity. Furthermore, we showed that the remaining genes with damaging variants in CDH significantly overlap with genes implicated in other developmental disorders. Gene expression patterns and patient phenotypes support pleiotropic effects of damaging variants in these genes on CDH and other developmental disorders. Finally, functional enrichment analysis implicates the disruption of regulation of gene expression, kinase activities, intra-cellular signaling, and cytoskeleton organization as pathogenic mechanisms in CDH.

Paternally inherited cis-regulatory structural variants are associated with autism

The genetic basis of autism spectrum disorder (ASD) is known to consist of contributions from de novo mutations in variant-intolerant genes. We hypothesize that rare inherited structural variants in cis-regulatory elements (CRE-SVs) of these genes also contribute to ASD. We investigated this by assessing the evidence for natural selection and transmission distortion of CRE-SVs in whole genomes of 9274 subjects from 2600 families affected by ASD. In a discovery cohort of 829 families, structural variants were depleted within promoters and untranslated regions, and paternally inherited CRE-SVs were preferentially transmitted to affected offspring and not to their unaffected siblings. The association of paternal CRE-SVs was replicated in an independent sample of 1771 families. Our results suggest that rare inherited noncoding variants predispose children to ASD, with differing contributions from each parent.

Clinical Application of Chromosome Microarray Analysis in Han Chinese Children with Neurodevelopmental Disorders

Chromosome microarray analysis (CMA) is a cost-effective molecular cytogenetic technique that has been used as a first-line diagnostic test in neurodevelopmental disorders in the USA since 2011. The impact of CMA results on clinical practice in China is not yet well studied, so we aimed to better evaluate this phenomenon. We analyzed the CMA results from 434 patients in our clinic, and characterized their molecular diagnoses, clinical features, and follow-up clinical actions based on these results. The overall diagnostic yield for our patients was 13.6% (59 out of 434). This gave a detection rate of 14.7% for developmental delay/intellectual disability (DD/ID, 38/259) and 12% for autism spectrum disorders (ASDs, 21/175). Thirty-three recurrent (n ≥ 2) variants were found, distributed at six chromosomal loci involving known chromosome syndromes (such as DiGeorge, Williams Beuren, and Angelman/Prader-Willi syndromes). The spectrum of positive copy number variants in our study was comparable to that reported in Caucasian populations, but with specific characteristics. Parental origin tests indicated an effect involving a significant maternal transmission bias to sons. The majority of patients with positive results (94.9%) had benefits, allowing earlier diagnosis (36/59), prioritized full clinical management (28/59), medication changes (7/59), a changed prognosis (30/59), and prenatal genetic counseling (15/59). Our results provide information on de novo mutations in Chinese children with DD/ID and/or ASDs. Our data showed that microarray testing provides immediate clinical utility for patients. It is expected that the personalized medical care of children with developmental disabilities will lead to improved outcomes in long-term developmental potential. We advocate using the diagnostic yield of clinically actionable results to evaluate CMA as it provides information of both clinical validity and clinical utility.

CAGI4 SickKids clinical genomes challenge: A pipeline for identifying pathogenic variants

Compared with earlier more restricted sequencing technologies, identification of rare disease variants using whole-genome sequence has the possibility of finding all causative variants, but issues of data quality and an overwhelming level of background variants complicate the analysis. The CAGI4 SickKids clinical genome challenge provided an opportunity to assess the landscape of variants found in a difficult set of 25 unsolved rare disease cases. To address the challenge, we developed a three-stage pipeline, first carefully analyzing data quality, then classifying high-quality gene-specific variants into seven categories, and finally examining each candidate variant for compatibility with the often complex phenotypes of these patients for final prioritization. Variants consistent with the phenotypes were found in 24 out of the 25 cases, and in a number of these, there are prioritized variants in multiple genes. Data quality analysis suggests that some of the selected variants are likely incorrect calls, complicating interpretation. The data providers followed up on three suggested variants with Sanger sequencing, and in one case, a prioritized variant was confirmed as likely causative by the referring physician, providing a diagnosis in a previously intractable case.