In Silico Prediction of Deleteriousness for Nonsynonymous and Splice-Altering Single Nucleotide Variants in the Human Genome

Chromatin regulators in the TBX1 network confer risk for conotruncal heart defects in 22q11.2DS

Congenital heart disease (CHD) affecting the conotruncal region of the heart, occurs in 40-50% of patients with 22q11.2 deletion syndrome (22q11.2DS). This syndrome is a rare disorder with relative genetic homogeneity that can facilitate identification of genetic modifiers. Haploinsufficiency of TBX1, encoding a T-box transcription factor, is one of the main genes responsible for the etiology of the syndrome. We suggest that genetic modifiers of conotruncal defects in patients with 22q11.2DS may be in the TBX1 gene network. To identify genetic modifiers, we analyzed rare, predicted damaging variants in whole genome sequence of 456 cases with conotruncal defects and 537 controls, with 22q11.2DS. We then performed gene set approaches and identified chromatin regulatory genes as modifiers. Chromatin genes with recurrent damaging variants include EP400, KAT6A, KMT2C, KMT2D, NSD1, CHD7 and PHF21A. In total, we identified 37 chromatin regulatory genes, that may increase risk for conotruncal heart defects in 8.5% of 22q11.2DS cases. Many of these genes were identified as risk factors for sporadic CHD in the general population. These genes are co-expressed in cardiac progenitor cells with TBX1, suggesting that they may be in the same genetic network. The genes KAT6A, KMT2C, CHD7 and EZH2, have been previously shown to genetically interact with TBX1 in mouse models. Our findings indicate that disturbance of chromatin regulatory genes impact the TBX1 gene network serving as genetic modifiers of 22q11.2DS and sporadic CHD, suggesting that there are some shared mechanisms involving the TBX1 gene network in the etiology of CHD.

Possible Catch-Up Developmental Trajectories for Children with Mild Developmental Delay Caused by NAA15 Pathogenic Variants

Variants in NAA15 are closely related to neurodevelopmental disorders (NDDs). In this study, we investigated the spectrum and clinical features of NAA15 variants in a Chinese NDD cohort of 769 children. Four novel NAA15 pathogenic variants were detected by whole-exome sequencing, including three de novo variants and one maternal variant. The in vitro minigene splicing assay confirmed one noncanonical splicing variant (c.1410+5G>C), which resulted in abnormal mRNA splicing. All affected children presented mild developmental delay, and catch-up trajectories were noted in three patients based on their developmental scores at different ages. Meanwhile, the literature review also showed that half of the reported patients with NAA15 variants presented mild/moderate developmental delay or intellectual disability, and possible catch-up sign was indicated for three affected patients. Taken together, our study expanded the spectrum of NAA15 variants in NDD patients. The affected patients presented mild developmental delay, and possible catch-up developmental trajectories were suggested. Studying the natural neurodevelopmental trajectories of NDD patients with pathogenic variants and their benefits from physical rehabilitations are needed in the future for precise genetic counseling and clinical management.

The Somatic Mutation Paradigm in Congenital Malformations: Hirschsprung Disease as a Model

Patients with Hirschsprung disease (HSCR) do not always receive a genetic diagnosis after routine screening in clinical practice. One of the reasons for this could be that the causal mutation is not present in the cell types that are usually tested—whole blood, dermal fibroblasts or saliva—but is only in the affected tissue. Such mutations are called somatic, and can occur in a given cell at any stage of development after conception. They will then be present in all subsequent daughter cells. Here, we investigated the presence of somatic mutations in HSCR patients. For this, whole-exome sequencing and copy number analysis were performed in DNA isolated from purified enteric neural crest cells (ENCCs) and blood or fibroblasts of the same patient. Variants identified were subsequently validated by Sanger sequencing. Several somatic variants were identified in all patients, but causative mutations for HSCR were not specifically identified in the ENCCs of these patients. Larger copy number variants were also not found to be specific to ENCCs. Therefore, we believe that somatic mutations are unlikely to be identified, if causative for HSCR. Here, we postulate various modes of development following the occurrence of a somatic mutation, to describe the challenges in detecting such mutations, and hypothesize how somatic mutations may contribute to ‘missing heritability’ in developmental defects.

Identification of 22 susceptibility loci associated with testicular germ cell tumors

Testicular germ cell tumors (TGCT) are the most common tumor in young white men and have a high heritability. In this study, the international Testicular Cancer Consortium assemble 10,156 and 179,683 men with and without TGCT, respectively, for a genome-wide association study. This meta-analysis identifies 22 TGCT susceptibility loci, bringing the total to 78, which account for 44% of disease heritability. Men with a polygenic risk score (PRS) in the 95th percentile have a 6.8-fold increased risk of TGCT compared to men with median scores. Among men with independent TGCT risk factors such as cryptorchidism, the PRS may guide screening decisions with the goal of reducing treatment-related complications causing long-term morbidity in survivors. These findings emphasize the interconnected nature of two known pathways that promote TGCT susceptibility: male germ cell development within its somatic niche and regulation of chromosomal division and structure, and implicate an additional biological pathway, mRNA translation.

Developmental delay with hypotrophy associated with homozygous functionally relevant REV3L variant

REV3L encodes a catalytic subunit of DNA polymerase zeta (Pol zeta) which is essential for the tolerance of DNA damage by inducing translesion synthesis (TLS). So far, the only Mendelian disease associated with REV3L was Moebius syndrome (3 patients with dominant REV3L mutations causing monoallelic loss-of-function were reported). We describe a homozygous ultra-rare REV3L variant (T2753R) identified with whole exome sequencing in a child without Moebius syndrome but with developmental delay, hypotrophy, and dysmorphic features who was born to healthy parents (heterozygous carriers of the variant). The variant affects the amino acid adjacent to functionally important KKRY motif. By introducing an equivalent mutation (S1192R) into the REV3 gene in yeasts, we showed that, whereas it retained residual function, it caused clear dysfunction of TLS in the nucleus and instability of mitochondrial genetic information. In particular, the mutation increased UV sensitivity measured by cell survival, decreased both the spontaneous (P < 0.005) and UV-induced (P < 0.0001) mutagenesis rates of nuclear DNA and increased the UV-induced mutagenesis rates of mitochondrial DNA (P < 0.0005). We propose that our proband is the first reported case of a REV3L associated disease different from Moebius syndrome both in terms of clinical manifestations and inheritance (autosomal recessive rather than dominant). KEY MESSAGES: First description of a human recessive disorder associated with a REV3L variant. A study in yeast showed that the variant affected the enzymatic function of the protein. In particular, it caused increased UV sensitivity and abnormal mutagenesis rates.

The role of exome sequencing in newborn screening for inborn errors of metabolism

Public health newborn screening (NBS) programs provide population-scale ascertainment of rare, treatable conditions that require urgent intervention. Tandem mass spectrometry (MS/MS) is currently used to screen newborns for a panel of rare inborn errors of metabolism (IEMs)1-4. The NBSeq project evaluated whole-exome sequencing (WES) as an innovative methodology for NBS. We obtained archived residual dried blood spots and data for nearly all IEM cases from the 4.5 million infants born in California between mid-2005 and 2013 and from some infants who screened positive by MS/MS, but were unaffected upon follow-up testing. WES had an overall sensitivity of 88% and specificity of 98.4%, compared to 99.0% and 99.8%, respectively for MS/MS, although effectiveness varied among individual IEMs. Thus, WES alone was insufficiently sensitive or specific to be a primary screen for most NBS IEMs. However, as a secondary test for infants with abnormal MS/MS screens, WES could reduce false-positive results, facilitate timely case resolution and in some instances even suggest more appropriate or specific diagnosis than that initially obtained. This study represents the largest, to date, sequencing effort of an entire population of IEM-affected cases, allowing unbiased assessment of current capabilities of WES as a tool for population screening.

A novel splicing mutation of ARHGAP29 is associated with nonsyndromic cleft lip with or without cleft palate

BACKGROUND

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is one of the most common birth defects, and occurs in approximately 1/700 live births worldwide. The correlation between the ABCA4-ARHGAP29 region and NSCL/P was first identified by genome-wide association studies (GWAS), but few reports have examined NSCL/P caused by ARHGAP29 mutations in the Chinese population.

METHODS

We performed chromosome microarray analysis (CMA) for two consecutive abnormal fetuses and whole exome sequencing (WES) for the family, including 3 patients and 2 normal family members, Sanger sequencing and RT-PCR were used to confirm the mutation.

RESULTS

We identified a novel splice donor mutation (ARHGAP29 c.1920 + 1G > A) in two consecutive NSCL/P fetuses, and the variant was inherited from the mother and grandfather. The mutation caused abnormal skipping of exon 17, and the mRNA level of ARHGAP29 was significantly decreased compared to the wild type.

CONCLUSIONS

In this study, we successfully diagnosed the genetic cause of NSCL/P in a family and first report that the c.1920 + 1G > A mutation in ARHGAP29 is associated with NSCL/P. Our study enriches the genetic landscape of NSCL/P, extends the mutation spectrum of ARHGAP29, and provides a new direction for the diagnosis of NSCL/P in patients and its prenatal diagnosis in fetuses.

Human CRY1 variants associate with attention deficit/hyperactivity disorder

Attention deficit/hyperactivity disorder (ADHD) is a common and heritable phenotype frequently accompanied by insomnia, anxiety, and depression. Here, using a reverse phenotyping approach, we report heterozygous coding variations in the core circadian clock gene cryptochrome 1 in 15 unrelated multigenerational families with combined ADHD and insomnia. The variants led to functional alterations in the circadian molecular rhythms, providing a mechanistic link to the behavioral symptoms. One variant, CRY1Δ11 c.1657+3A>C, is present in approximately 1% of Europeans, therefore standing out as a diagnostic and therapeutic marker. We showed by exome sequencing in an independent cohort of patients with combined ADHD and insomnia that 8 of 62 patients and 0 of 369 controls carried CRY1Δ11. Also, we identified a variant, CRY1Δ6 c.825+1G>A, that shows reduced affinity for BMAL1/CLOCK and causes an arrhythmic phenotype. Genotype-phenotype correlation analysis revealed that this variant segregated with ADHD and delayed sleep phase disorder (DSPD) in the affected family. Finally, we found in a phenome-wide association study involving 9438 unrelated adult Europeans that CRY1Δ11 was associated with major depressive disorder, insomnia, and anxiety. These results defined a distinctive group of circadian psychiatric phenotypes that we propose to designate as "circiatric" disorders.

Gene panel for Mendelian strokes

BACKGROUND

Mendelian stroke causes nearly 7% of ischaemic strokes and is also an important aetiology of cryptogenic stroke. Identifying the genetic abnormalities in Mendelian strokes is important as it would facilitate therapeutic management and genetic counselling. Next-generation sequencing makes large-scale sequencing and genetic testing possible.

METHODS

A systematic literature search was conducted to identify causal genes of Mendelian strokes, which were used to construct a hybridization-based gene capture panel. Genetic variants for target genes were detected using Illumina HiSeq X10 and the Novaseq platform. The sensitivity and specificity were evaluated by comparing the results with Sanger sequencing.

RESULTS

53 suspected patients of Mendelian strokes were analysed using the panel of 181 causal genes. According to the American College of Medical Genetics and Genomics standard, 16 likely pathogenic/variants of uncertain significance genetic variants were identified. Diagnostic testing was conducted by comparing the consistency between the results of panel and Sanger sequencing. Both the sensitivity and specificity were 100% for the panel.

CONCLUSION

This panel provides an economical, time-saving and labour-saving method to detect causal mutations of Mendelian strokes.

Findings from a Genotyping Study of Over 1000 People with Inherited Retinal Disorders in Ireland

The Irish national registry for inherited retinal degenerations (Target 5000) is a clinical and scientific program to identify individuals in Ireland with inherited retinal disorders and to attempt to ascertain the genetic cause underlying the disease pathology. Potential participants first undergo a clinical assessment, which includes clinical history and analysis with multimodal retinal imaging, electrophysiology, and visual field testing. If suitable for recruitment, a sample is taken and used for genetic analysis. Genetic analysis is conducted by use of a retinal gene panel target capture sequencing approach. With over 1000 participants from 710 pedigrees now screened, there is a positive candidate variant detection rate of approximately 70% (495/710). Where an autosomal recessive inheritance pattern is observed, an additional 9% (64/710) of probands have tested positive for a single candidate variant. Many novel variants have also been detected as part of this endeavor. The target capture approach is an economic and effective means of screening patients with inherited retinal disorders. Despite the advances in sequencing technology and the ever-decreasing associated processing costs, target capture remains an attractive option as the data produced is easily processed, analyzed, and stored compared to more comprehensive methods. However, with decreasing costs of whole genome and whole exome sequencing, the focus will likely move towards these methods for more comprehensive data generation.

Biallelic DMXL2 mutations impair autophagy and cause Ohtahara syndrome with progressive course

Esposito et al. identify biallelic loss-of-function mutations in DMXL2, encoding a v-ATPase regulatory protein, in three sibling pairs exhibiting Ohtahara syndrome with a progressive course. Patient-derived fibroblasts and Dmxl2-silenced mouse hippocampal neurons show defective lysosomal function and autophagy, resulting in the latter in impaired neuronal development and synapse formation.

Future Preventive Gene Therapy of Polygenic Diseases from a Population Genetics Perspective

With the accumulation of scientific knowledge of the genetic causes of common diseases and continuous advancement of gene-editing technologies, gene therapies to prevent polygenic diseases may soon become possible. This study endeavored to assess population genetics consequences of such therapies. Computer simulations were used to evaluate the heterogeneity in causal alleles for polygenic diseases that could exist among geographically distinct populations. The results show that although heterogeneity would not be easily detectable by epidemiological studies following population admixture, even significant heterogeneity would not impede the outcomes of preventive gene therapies. Preventive gene therapies designed to correct causal alleles to a naturally-occurring neutral state of nucleotides would lower the prevalence of polygenic early- to middle-age-onset diseases in proportion to the decreased population relative risk attributable to the edited alleles. The outcome would manifest differently for late-onset diseases, for which the therapies would result in a delayed disease onset and decreased lifetime risk; however, the lifetime risk would increase again with prolonging population life expectancy, which is a likely consequence of such therapies. If the preventive heritable gene therapies were to be applied on a large scale, the decreasing frequency of risk alleles in populations would reduce the disease risk or delay the age of onset, even with a fraction of the population receiving such therapies. With ongoing population admixture, all groups would benefit over generations.

MoBiDiC Prioritization Algorithm, a Free, Accessible, and Efficient Pipeline for Single-Nucleotide Variant Annotation and Prioritization for Next-Generation Sequencing Routine Molecular Diagnosis

Interpretation of next-generation sequencing constitutes the main limitation of molecular diagnostics. In diagnosing myopathies and muscular dystrophies, another issue is efficiency in predicting the pathogenicity of variants identified in large genes, especially TTN; current in silico prediction tools show limitations in predicting and ranking the numerous variants of such genes. We propose a variant-prioritization tool, the MoBiDiCprioritization algorithm (MPA). MPA is based on curated interpretation of data on previously reported variants, biological assumptions, and splice and missense predictors, and is used to prioritize all types of single-nucleotide variants. MPA was validated by comparing its sensitivity and specificity to those of dbNSFP database prediction tools, using a data set composed of DYSF, DMD, LMNA, NEB, and TTN variants extracted from expert-reviewed and ExAC databases. MPA obtained the best annotation rates for missense and splice variants. As MPA aggregates the results from several predictors, individual predictor errors are counterweighted, improving the sensitivity and specificity of missense and splice variant predictions. We propose a sequential use of MPA, beginning with the selection of variants with higher scores and followed by, in the absence of candidate pathologic variants, consideration of variants with lower scores. We provide scripts and documentation for free academic use and a validated annotation pipeline scaled for panel and exome sequencing to prioritize single-nucleotide variants from a VCF file.

Combined genetic approaches yield a 48% diagnostic rate in a large cohort of French hearing-impaired patients

Hearing loss is the most common sensory disorder and because of its high genetic heterogeneity, implementation of Massively Parallel Sequencing (MPS) in diagnostic laboratories is greatly improving the possibilities of offering optimal care to patients. We present the results of a two-year period of molecular diagnosis that included 207 French families referred for non-syndromic hearing loss. Our multi-step strategy involved (i) DFNB1 locus analysis, (ii) MPS of 74 genes, and (iii) additional approaches including Copy Number Variations, in silico analyses, minigene studies coupled when appropriate with complete gene sequencing, and a specific assay for STRC. This comprehensive screening yielded an overall diagnostic rate of 48%, equally distributed between DFNB1 (24%) and the other genes (24%). Pathogenic genotypes were identified in 19 different genes, with a high prevalence of GJB2, STRC, MYO15A, OTOF, TMC1, MYO7A and USH2A. Involvement of an Usher gene was reported in 16% of the genotyped cohort. Four de novo variants were identified. This study highlights the need to develop several molecular approaches for efficient molecular diagnosis of hearing loss, as this is crucial for genetic counselling, audiological rehabilitation and the detection of syndromic forms.