Comprehensive de novo mutation discovery with HiFi long-read sequencing

BACKGROUND

Long-read sequencing (LRS) techniques have been very successful in identifying structural variants (SVs). However, the high error rate of LRS made the detection of small variants (substitutions and short indels < 20 bp) more challenging. The introduction of PacBio HiFi sequencing makes LRS also suited for detecting small variation. Here we evaluate the ability of HiFi reads to detect de novo mutations (DNMs) of all types, which are technically challenging variant types and a major cause of sporadic, severe, early-onset disease.

METHODS

We sequenced the genomes of eight parent-child trios using high coverage PacBio HiFi LRS (~ 30-fold coverage) and Illumina short-read sequencing (SRS) (~ 50-fold coverage). De novo substitutions, small indels, short tandem repeats (STRs) and SVs were called in both datasets and compared to each other to assess the accuracy of HiFi LRS. In addition, we determined the parent-of-origin of the small DNMs using phasing.

RESULTS

We identified a total of 672 and 859 de novo substitutions/indels, 28 and 126 de novo STRs, and 24 and 1 de novo SVs in LRS and SRS respectively. For the small variants, there was a 92 and 85% concordance between the platforms. For the STRs and SVs, the concordance was 3.6 and 0.8%, and 4 and 100% respectively. We successfully validated 27/54 LRS-unique small variants, of which 11 (41%) were confirmed as true de novo events. For the SRS-unique small variants, we validated 42/133 DNMs and 8 (19%) were confirmed as true de novo event. Validation of 18 LRS-unique de novo STR calls confirmed none of the repeat expansions as true DNM. Confirmation of the 23 LRS-unique SVs was possible for 19 candidate SVs of which 10 (52.6%) were true de novo events. Furthermore, we were able to assign 96% of DNMs to their parental allele with LRS data, as opposed to just 20% with SRS data.

CONCLUSIONS

HiFi LRS can now produce the most comprehensive variant dataset obtainable by a single technology in a single laboratory, allowing accurate calling of substitutions, indels, STRs and SVs. The accuracy even allows sensitive calling of DNMs on all variant levels, and also allows for phasing, which helps to distinguish true positive from false positive DNMs.

Large-scale analysis of de novo mutations identifies risk genes for female infertility characterized by oocyte and early embryo defects

BACKGROUND

Oocyte maturation arrest and early embryonic arrest are important reproductive phenotypes resulting in female infertility and cause the recurrent failure of assisted reproductive technology (ART). However, the genetic etiologies of these female infertility-related phenotypes are poorly understood. Previous studies have mainly focused on inherited mutations based on large pedigrees or consanguineous patients. However, the role of de novo mutations (DNMs) in these phenotypes remains to be elucidated.

RESULTS

To decipher the role of DNMs in ART failure and female infertility with oocyte and embryo defects, we explore the landscape of DNMs in 473 infertile parent-child trios and identify a set of 481 confident DNMs distributed in 474 genes. Gene ontology analysis reveals that the identified genes with DNMs are enriched in signaling pathways associated with female reproductive processes such as meiosis, embryonic development, and reproductive structure development. We perform functional assays on the effects of DNMs in a representative gene Tubulin Alpha 4a (TUBA4A), which shows the most significant enrichment of DNMs in the infertile parent-child trios. DNMs in TUBA4A disrupt the normal assembly of the microtubule network in HeLa cells, and microinjection of DNM TUBA4A cRNAs causes abnormalities in mouse oocyte maturation or embryo development, suggesting the pathogenic role of these DNMs in TUBA4A.

CONCLUSIONS

Our findings suggest novel genetic insights that DNMs contribute to female infertility with oocyte and embryo defects. This study also provides potential genetic markers and facilitates the genetic diagnosis of recurrent ART failure and female infertility.

Rapid survey of de novo mutations in naturally growing tree species following the March 2011 disaster in Fukushima: The effect of low-dose-rate radiation

The impact of low-dose-rate radiation on genetics is largely unknown, particularly in natural environments. The Fukushima Dai-ich Nuclear Power Plant disaster resulted in the creation of contaminated natural lands. In this study, de novo mutations (DNMs) in germ line cells were surveyed from double-digest RADseq fragments in Japanese cedar and flowering cherry trees exposed to ambient dose rates ranging from 0.08 to 6.86 μGy h^-1. These two species are among the most widely cultivated Japanese gymnosperm and angiosperm trees for forestry and horticultural purpose, respectively. For Japanese flowering cherry, open crossings were performed to produce seedlings, and only two candidate DNMs were detected from uncontaminated area. For Japanese cedar, the haploid megagametophytes were used as next generation samples. The use of megagametophytes from open crossing for next generation mutation screening had many advantages such as reducing exposure to radiation in contaminated areas because artificial crossings are not needed and the ease of data analysis owing to the haploid nature of megagametophytes. A direct comparison of the nucleotide sequences of parents and megagametophytes revealed an average of 1.4 candidate DNMs per megagametophyte sample (range: 0-40) after filtering procedures were optimized based on the validation of DNMs via Sanger sequencing. There was no relationship between the observed mutations and the ambient dose rate in the growing area or the concentration of ¹³⁷Cs in cedar branches. The present results also suggest that mutation rates differ among lineages and that the growing environment has a relatively large influence on these mutation rates. These results suggested there was no significant increase in the mutation rate of the germplasm of Japanese cedar and flowering cherry trees growing in the contaminated areas.

"Mutation blacklist" and "mutation whitelist" of SARS-CoV-2

Over the past two years, scientists throughout the world have completed more than 6 million SARS-CoV-2 genome sequences. Today, the number of SARS-CoV-2 genomes exceeds the total number of all other viral genomes. These genomes are a record of the evolution of SARS-CoV-2 in the human host, and provide information on the emergence of mutations. In this study, analysis of these sequenced genomes identified 296,728 de novo mutations (DNMs), and found that six types of base substitutions reached saturation in the sequenced genome population. Based on this analysis, a "mutation blacklist" of SARS-CoV-2 was compiled. The loci on the "mutation blacklist" are highly conserved, and these mutations likely have detrimental effects on virus survival, replication, and transmission. This information is valuable for SARS-CoV-2 research on gene function, vaccine design, and drug development. Through association analysis of DNMs and viral transmission rates, we identified 185 DNMs that positively correlated with the SARS-CoV-2 transmission rate, and these DNMs where classified as the "mutation whitelist" of SARS-CoV-2. The mutations on the "mutation whitelist" are beneficial for SARS-CoV-2 transmission and could therefore be used to evaluate the transmissibility of new variants. The occurrence of mutations and the evolution of viruses are dynamic processes. To more effectively monitor the mutations and variants of SARS-CoV-2, we built a SARS-CoV-2 mutation and variant monitoring and pre-warning system (MVMPS), which can monitor the occurrence and development of mutations and variants of SARS-CoV-2, as well as provide pre-warning for the prevention and control of SARS-CoV-2 (https://www.omicx.cn/). Additionally, this system could be used in real-time to update the "mutation whitelist" and "mutation blacklist" of SARS-CoV-2.

De novo mutations within metabolism networks of amino acid/protein/energy in Chinese autistic children with intellectual disability

BACKGROUND

Autism spectrum disorder (ASD) is often accompanied by intellectual disability (ID). Despite extensive studies, however, the genetic basis for this comorbidity is still not clear. In this study, we tried to develop an analyzing pipeline for de novo mutations and possible pathways related to ID phenotype in ASD. Whole-exome sequencing (WES) was performed to screen de novo mutations and candidate genes in 79 ASD children together with their parents (trios). The de novo altering genes and relative pathways which were associated with ID phenotype were analyzed. The connection nodes (genes) of above pathways were selected, and the diagnostic value of these selected genes for ID phenotype in the study population was also evaluated.

RESULTS

We identified 89 de novo mutant genes, of which 34 genes were previously reported to be associated with ASD, including double hits in the EGF repeats of NOTCH1 gene (p.V999M and p.S1027L). Interestingly, of these 34 genes, 22 may directly affect intelligence quotient (IQ). Further analyses revealed that these IQ-related genes were enriched in protein synthesis, energy metabolism, and amino acid metabolism, and at least 9 genes (CACNA1A, ALG9, PALM2, MGAT4A, PCK2, PLEKHA1, PSME3, ADI1, and TLE3) were involved in all these three pathways. Seven patients who harbored these gene mutations showed a high prevalence of a low IQ score (< 70), a non-verbal language, and an early diagnostic age (< 4 years). Furthermore, our panel of these 9 genes reached a 10.2% diagnostic rate (5/49) in early diagnostic patients with a low IQ score and also reached a 10% diagnostic yield in those with both a low IQ score and non-verbal language (4/40).

CONCLUSION

We found some new genetic disposition for ASD accompanied with intellectual disability in this study. Our results may be helpful for etiologic research and early diagnoses of intellectual disability in ASD. Larger population studies and further mechanism studies are warranted.

[Exploring the association between de novo mutations and non-syndromic cleft lip with or without palate based on whole exome sequencing of case-parent trios]

OBJECTIVE

To explore the association between de novo mutations (DNM) and non-syndromic cleft lip with or without palate (NSCL/P) using case-parent trio design.

METHODS

Whole-exome sequencing was conducted for twenty-two NSCL/P trios and Genome Analysis ToolKit (GATK) was used to identify DNM by comparing the alleles of the cases and their parents. Information of predictable functions was annotated to the locus with SnpEff. Enrichment analysis for DNM was conducted to test the difference between the actual number and the expected number of DNM, and to explore whether there were genes with more DNM than expected. NSCL/P-related genes indicated by previous studies with solid evidence were selected by literature reviewing. Protein-protein interactions analysis was conducted among the genes with protein-altering DNM and NSCL/P-related genes. R package "denovolyzeR" was used for the enrichment analysis (Bonferroni correction: P=0.05/n, n is the number of genes in the whole genome range). Protein-protein interactions among genes with DNM and genes with solid evidence on the risk factors of NSCL/P were predicted depending on the information provided by STRING database.

RESULTS

A total of 339 908 SNPs were qualified for the subsequent analysis after quality control. The number of high confident DNM identified by GATK was 345. Among those DNM, forty-four DNM were missense mutations, one DNM was nonsense mutation, two DNM were splicing site mutations, twenty DNM were synonymous mutations and others were located in intron or intergenic regions. The results of enrichment analysis showed that the number of protein-altering DNM on the exome regions was larger than expected (P < 0.05), and five genes (KRTCAP2, HMCN2, ANKRD36C, ADGRL2 and DIPK2A) had more DNM than expected (P < 0.05/(2×19 618)). Protein-protein interaction analysis was conducted among forty-six genes with protein-altering DNM and thirteen genes associated with NSCL/P selected by literature reviewing. Six pairs of interactions occurred between the genes with DNM and known NSCL/P-related genes. The score measuring the confidence level of the predicted interaction between RGPD4 and SUMO1 was 0.868, which was higher than the scores for other pairs of genes.

CONCLUSION

Our study provided novel insights into the development of NSCL/P and demonstrated that functional analyses of genes carrying DNM were warranted to understand the genetic architecture of complex diseases.

Protein phosphatase 2A - structure, function and role in neurodevelopmental disorders

Neurodevelopmental disorders (NDDs), including intellectual disability (ID), autism and schizophrenia, have high socioeconomic impact, yet poorly understood etiologies. A recent surge of large-scale genome or exome sequencing studies has identified a multitude of mostly de novo mutations in subunits of the protein phosphatase 2A (PP2A) holoenzyme that are strongly associated with NDDs. PP2A is responsible for at least 50% of total Ser/Thr dephosphorylation in most cell types and is predominantly found as trimeric holoenzymes composed of catalytic (C), scaffolding (A) and variable regulatory (B) subunits. PP2A can exist in nearly 100 different subunit combinations in mammalian cells, dictating distinct localizations, substrates and regulatory mechanisms. PP2A is well established as a regulator of cell division, growth, and differentiation, and the roles of PP2A in cancer and various neurodegenerative disorders, such as Alzheimer's disease, have been reviewed in detail. This Review summarizes and discusses recent reports on NDDs associated with mutations of PP2A subunits and PP2A-associated proteins. We also discuss the potential impact of these mutations on the structure and function of the PP2A holoenzymes and the etiology of NDDs.

Rare variants discovery by extensive whole-genome sequencing of the Han Chinese population in Taiwan: Applications to cardiovascular medicine

Introduction

A population-specific genomic reference is important for research and clinical practice, yet it remains unavailable for Han Chinese (HC) in Taiwan.

Objectives

We report the first whole genome sequencing (WGS) database of HC (1000 Taiwanese genome (1KTW-WGS)) and demonstrate several applications to cardiovascular medicine.

Methods

Whole genomes of 997 HC were sequenced to at least 30X depth. A total of 20,117 relatively healthy HC individuals were genotyped using a customized Axiom GWAS array. We performed a genome-wide genotype imputation technique using IMPUTE2.

Results

We identified 26.7 million single-nucleotide variants (SNVs) and 4.2 million insertions-deletions. Of the SNVs, 16.1% were novel relative to dbSNP (build 152), and 34.2% were novel relative to gnomAD. A total of 18,450 healthy HC individuals were genotyped using a customized Genome-Wide Association Study (GWAS) array. We identified hypertension-associated variants and developed a hypertension prediction model based on the correlation between the WGS data and GWAS data (combined clinical and genetic models, AUC 0.887), and also identified 3 novel hyperlipidemia-associated variants. Each individual carried an average of 16.42 (SD = 3.72) disease-causing variants. Additionally, we established an online SCN5A (an important cardiac gene) database that can be used to explore racial differences. Finally, pharmacogenetics studies identified HC population-specific SNVs in genes (CYP2C9 and VKORC1) involved in drug metabolism and blood clotting.

Conclusion

This research demonstrates the benefits of constructing a population-specific genomic reference database for precision medicine.

Targeted sequencing and integrative analysis of 3,195 Chinese patients with neurodevelopmental disorders prioritized 26 novel candidate genes

Neurodevelopmental disorders (NDDs) are a set of complex disorders characterized by diverse and co-occurring clinical symptoms. The genetic contribution in patients with NDDs remains largely unknown. Here, we sequence 519 NDD-related genes in 3,195 Chinese probands with neurodevelopmental phenotypes and identify 2,522 putative functional mutations consisting of 137 de novo mutations (DNMs) in 86 genes and 2,385 rare inherited mutations (RIMs) with 22 X-linked hemizygotes in 13 genes, 2 homozygous mutations in 2 genes and 23 compound heterozygous mutations in 10 genes. Furthermore, the DNMs of 16,807 probands with NDDs are retrieved from public datasets and combine in an integrated analysis with the mutation data of our Chinese NDD probands by taking 3,582 in-house controls of Chinese origin as background. We prioritize 26 novel candidate genes. Notably, six of these genes - ITSN1, UBR3, CADM1, RYR3, FLNA, and PLXNA3 - preferably contribute to autism spectrum disorders (ASDs), as demonstrated by high co-expression and/or interaction with ASD genes confirmed via rescue experiments in a mouse model. Importantly, these genes are differentially expressed in the ASD cortex in a significant manner and involved in ASD-associated networks. Together, our study expands the genetic spectrum of Chinese NDDs, further facilitating both basic and translational research.

Antenatal antioxidants to avert autism?

Paternally derived de novo mutations (DNMs) caused by oxidative stress (OS) have been implicated in the development of autism spectrum disorders (ASDs). Whether preconception antioxidant supplementation can reduce the incidence of ASDs by reducing OS is an area of uncertainty and potentially important future scientific investigation.

Genetics of Childhood-Onset Schizophrenia: 2019 Update

The genetic architecture of schizophrenia is complex and highly polygenic. This article discusses key findings from genetic studies of childhood-onset schizophrenia (COS) and the more common adult-onset schizophrenia (AOS), including studies of familial aggregation and common, rare, and copy number variants. Extant literature suggests that COS is a rare variant of AOS involving greater familial aggregation of schizophrenia spectrum disorders and a potentially higher occurrence of pathogenic copy number variants. The direct utility of genetics to clinical practice for COS is currently limited; however, identifying common pathways through which risk genes affect brain function offers promise for novel interventions.

Detection of de novo genetic variants in Mayer-Rokitansky-Küster-Hauser syndrome by whole genome sequencing

Objective

The aim of this study was to use whole genome sequencing (WGS) help detect de novo mutations or pathogenic genes of Mayer-Rokitansky-Küster-Hauser syndrome type 1(MRKH syndrome type 1).

Study design

This was a case-parent trios study. Nine unrelated probands, with MRKH syndrome type 1 and their parents were enrolled. The enrollment, sequencing process, establishment of the de novo mutations detecting procedure and experiment part were performed over a 2-year period.

Results

we detected 632 de novo single nucleotide variants (SNVs), 267 de novo small insertions/deletions (indels), 39 de novo structural variations (SVs) and 28 de novo copy number alterations (CNAs). Three novel damaging coding de novo SNVs with three damaging coding de novo genes (PIK3CD, SLC4A10 and TNK2) were revealed. Two SNVs were annotated of the promoter region of gene NBPF10 and 3'UTR of NOTCH2NL, potentially contributing to the pathogenesis of MRKH.

Conclusion

We identified five de novo mutations in BAZ2B, KLHL18, PIK3CD, SLC4A10 and TNK2 by performing WGS, the functional involvement of all deleterious mutations in MRKH candidate genes of the trios warrant further study. WGS may complement conventional array to capture the complete landscape of the genome in MRKH.

Failure to attune to language predicts autism in high risk infants

Young humans are typically sensitive to evolutionarily important aspects of information in the surrounding environment in a way that makes us thrive. Seeking to probe the putative disruptions of this process in infancy, I examined the statistical character of head movements in 52 9-10 mo-old infants, half at high familial risk (HR) for Autism Spectrum Disorders (ASD), who underwent an fMRI scan while listening to words spoken with alternating stress patterns on syllables. Relative to low risk (LR) infants, HR infants, in particular those showing the least rapid receptive language progress, had significantly lower noise-to-signal levels and increased symmetry. A comparison of patterns during a native language and a sleep scan revealed the most atypical ordering of signatures on the 3 tasks in a subset of HR infants, suggesting that the biological mechanism of language development is least acquisitive in those HR infants who go on to develop ASD in toddlerhood.

A structural look at GABAA receptor mutations linked to epilepsy syndromes

Understanding the genetic variation in GABA_A receptor subunit genes (GABRs), GABRA1-6, GABRB1-3, GABRG1-3 and GABRD, in individuals affected by epilepsy may improve the diagnosis and treatment of epilepsy syndromes through identification of disease-associated variants. However, the lack of functional analysis and validation of many novel and previously reported familial and de novo mutations have made it challenging to address meaningful gene associations with epilepsy syndromes. GABA_A receptors belong to the Cys-loop receptor family. Even though GABA_A receptor mutant residues are widespread among different GABRs, their frequent occurrence in important structural domains that share common functional features suggests associations between structure and function.

Inherited and multiple de novo mutations in autism/developmental delay risk genes suggest a multifactorial model

Background

We previously performed targeted sequencing of autism risk genes in probands from the Autism Clinical and Genetic Resources in China (ACGC) (phase I). Here, we expand this analysis to a larger cohort of patients (ACGC phase II) to better understand the prevalence, inheritance, and genotype-phenotype correlations of likely gene-disrupting (LGD) mutations for autism candidate genes originally identified in cohorts of European descent.

Methods

We sequenced 187 autism candidate genes in an additional 784 probands and 85 genes in 599 probands using single-molecule molecular inversion probes. We tested the inheritance of potentially pathogenic mutations, performed a meta-analysis of phase I and phase II data and combined our results with existing exome sequence data to investigate the phenotypes of carrier parents and patients with multiple hits in different autism risk genes.

Results

We validated recurrent, LGD, de novo mutations (DNMs) in 13 genes. We identified a potential novel risk gene (ZNF292), one novel gene with recurrent LGD DNMs (RALGAPB), as well as genes associated with macrocephaly (GIGYF2 and WDFY3). We identified the transmission of private LGD mutations in genes predominantly associated with DNMs and showed that parental carriers tended to share milder autism-related phenotypes. Patients that carried DNMs in two or more candidate genes show more severe phenotypes.

Conclusions

We identify new risk genes and transmission of deleterious mutations in genes primarily associated with DNMs. The fact that parental carriers show milder phenotypes and patients with multiple hits are more severe supports a multifactorial model of risk.

Genomic landscapes of Chinese sporadic autism spectrum disorders revealed by whole-genome sequencing

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with considerable clinical and genetic heterogeneity. In this study, we identified all classes of genomic variants from whole-genome sequencing (WGS) dataset of 32 Chinese trios with ASD, including de novo mutations, inherited variants, copy number variants (CNVs) and genomic structural variants. A higher mutation rate (Poisson test, P < 2.2 × 10^-16) in exonic (1.37 × 10^-8) and 3'-UTR regions (1.42 × 10^-8) was revealed in comparison with that of whole genome (1.05 × 10^-8). Using an integrated model, we identified 87 potentially risk genes (P < 0.01) from 4832 genes harboring various rare deleterious variants, including CHD8 and NRXN2, implying that the disorders may be in favor to multiple-hit. In particular, frequent rare inherited mutations of several microcephaly-associated genes (ASPM, WDR62, and ZNF335) were found in ASD. In chromosomal structure analyses, we found four de novo CNVs and one de novo chromosomal rearrangement event, including a de novo duplication of UBE3A-containing region at 15q11.2-q13.1, which causes Angelman syndrome and microcephaly, and a disrupted TNR due to de novo chromosomal translocation t(1; 5)(q25.1; q33.2). Taken together, our results suggest that abnormalities of centrosomal function and chromatin remodeling of the microcephaly-associated genes may be implicated in pathogenesis of ASD. Adoption of WGS as a new yet efficient technique to illustrate the full genetic spectrum in complex disorders, such as ASD, could provide novel insights into pathogenesis, diagnosis and treatment.

A 73,128 bp de novo deletion encompassing the OPN1LW/OPN1MW gene cluster in sporadic Blue Cone Monochromacy: a case report

Background

Blue Cone Monochromacy (BCM) is a rare congenital cone dysfunction disorder with X-linked recessive mode of inheritance. BCM is caused by mutations at the OPN1LW/MW cone opsin gene cluster including deletions of the locus control region (LCR) and/or parts of the gene cluster. We aimed at investigating the clinical presentation, genetic cause and inheritance underlying a sporadic case of BCM.

Case presentation

We report a 24-year-old male presenting with congenital photophobia, nystagmus and colour vision abnormalities. There was no history of retinal dystrophy in the family. Clinical diagnosis of BCM was supported by genetic investigations of the patient and his family members. Molecular genetic analysis of the OPN1LW/OPN1MW gene cluster revealed a novel deletion of about 73 kb in the patient encompassing the LCR. The deletion was absent in the X-chromosomes of both the mother and transmitting grandfather.

Conclusions

The present report provides the clinical findings and the genetic basis underlying a sporadic BCM case which is caused by a de novo deletion within the OPN1LW/MW gene cluster originating from the mother’s germline due to Alu-repeat mediated recombination. This is the first report of a de novo deletion resulting in BCM, highlighting the importance to consider BCM and perform genetic testing for this condition in male patients with cone dysfunction also in the absence of a positive family history.

Discovering human germ cell mutagens with whole genome sequencing: Insights from power calculations reveal the importance of controlling for between-family variability

Mutations in germ cells pose potential genetic risks to offspring. However, de novo mutations are rare events that are spread across the genome and are difficult to detect. Thus, studies in this area have generally been under-powered, and no human germ cell mutagen has been identified. Whole Genome Sequencing (WGS) of human pedigrees has been proposed as an approach to overcome these technical and statistical challenges. WGS enables analysis of a much wider breadth of the genome than traditional approaches. Here, we performed power analyses to determine the feasibility of using WGS in human families to identify germ cell mutagens. Different statistical models were compared in the power analyses (ANOVA and multiple regression for one-child families, and mixed effect model sampling between two to four siblings per family). Assumptions were made based on parameters from the existing literature, such as the mutation-by-paternal age effect. We explored two scenarios: a constant effect due to an exposure that occurred in the past, and an accumulating effect where the exposure is continuing. Our analysis revealed the importance of modeling inter-family variability of the mutation-by-paternal age effect. Statistical power was improved by models accounting for the family-to-family variability. Our power analyses suggest that sufficient statistical power can be attained with 4-28 four-sibling families per treatment group, when the increase in mutations ranges from 40 to 10% respectively. Modeling family variability using mixed effect models provided a reduction in sample size compared to a multiple regression approach. Much larger sample sizes were required to detect an interaction effect between environmental exposures and paternal age. These findings inform study design and statistical modeling approaches to improve power and reduce sequencing costs for future studies in this area.

Reproductive genetics and the aging male

PURPOSE

To examine current evidence of the known effects of advanced paternal age on sperm genetic and epigenetic changes and associated birth defects and diseases in offspring.

METHODS

Review of published PubMed literature.

RESULTS

Advanced paternal age (> 40 years) is associated with accumulated damage to sperm DNA and mitotic and meiotic quality control mechanisms (mismatch repair) during spermatogenesis. This in turn causes well-delineated abnormalities in sperm chromosomes, both numerical and structural, and increased sperm DNA fragmentation (3%/year of age) and single gene mutations (relative risk, RR 10). An increase in related abnormalities in offspring has also been described, including miscarriage (RR 2) and fetal loss (RR 2). There is also a significant increase in rare, single gene disorders (RR 1.3 to 12) and congenital anomalies (RR 1.2) in offspring. Current research also suggests that autism, schizophrenia, and other forms of "psychiatric morbidity" are more likely in offspring (RR 1.5 to 5.7) with advanced paternal age. Genetic defects related to faulty sperm quality control leading to single gene mutations and epigenetic alterations in several genetic pathways have been implicated as root causes.

CONCLUSIONS

Advanced paternal age is associated with increased genetic and epigenetic risk to offspring. However, the precise age at which risk develops and the magnitude of the risk are poorly understood or may have gradual effects. Currently, there are no clinical screenings or diagnostic panels that target disorders associated with advanced paternal age. Concerned couples and care providers should pursue or recommend genetic counseling and prenatal testing regarding specific disorders.

Comprehensive Analyses of Tissue-Specific Networks with Implications to Psychiatric Diseases

Recent advances in genome sequencing and "omics" technologies are opening new opportunities for improving diagnosis and treatment of human diseases. The precision medicine initiative in particular aims at developing individualized treatment options that take into account individual variability in genes and environment of each person. Systems biology approaches that group genes, transcripts and proteins into functionally meaningful networks will play crucial role in the future of personalized medicine. They will allow comparison of healthy and disease-affected tissues and organs from the same individual, as well as between healthy and disease-afflicted individuals. However, the field faces a multitude of challenges ranging from data integration to statistical and combinatorial issues in data analyses. This chapter describes computational approaches developed by us and the others to tackle challenges in tissue-specific network analyses, with the main focus on psychiatric diseases.

Two de novo variations identified by massively parallel sequencing in 13 Chinese families with children diagnosed with autism spectrum disorder

Autism spectrum disorder (ASD) is a genetically heterogeneous neurodevelopmental disorder characterized by impairments in social interaction and communication, and by restricted and repetitive behaviors. The genetic architecture of ASD has been elucidated, including chromosomal rearrangements, de novo or inherited rare variants, and copy number variants. However, the genetic mechanism of Chinese families with ASD children is explored rarely. To identify genetic pathogenesis, we performed massively parallel sequencing on 13 Chinese ASD trio families, and found two de novo variations. The novel de novo splice alteration c.664 + 2T > G in the DEAF1 gene and the novel de novo missense mutation c.95 C > T in the AADAT gene associated with ASD may be important clues for exploring the etiology of this disorder.

Integrated Bayesian analysis of rare exonic variants to identify risk genes for schizophrenia and neurodevelopmental disorders

BACKGROUND

Integrating rare variation from trio family and case-control studies has successfully implicated specific genes contributing to risk of neurodevelopmental disorders (NDDs) including autism spectrum disorders (ASD), intellectual disability (ID), developmental disorders (DDs), and epilepsy (EPI). For schizophrenia (SCZ), however, while sets of genes have been implicated through the study of rare variation, only two risk genes have been identified.

METHODS

We used hierarchical Bayesian modeling of rare-variant genetic architecture to estimate mean effect sizes and risk-gene proportions, analyzing the largest available collection of whole exome sequence data for SCZ (1,077 trios, 6,699 cases, and 13,028 controls), and data for four NDDs (ASD, ID, DD, and EPI; total 10,792 trios, and 4,058 cases and controls).

RESULTS

For SCZ, we estimate there are 1,551 risk genes. There are more risk genes and they have weaker effects than for NDDs. We provide power analyses to predict the number of risk-gene discoveries as more data become available. We confirm and augment prior risk gene and gene set enrichment results for SCZ and NDDs. In particular, we detected 98 new DD risk genes at FDR < 0.05. Correlations of risk-gene posterior probabilities are high across four NDDs (ρ>0.55), but low between SCZ and the NDDs (ρ<0.3). An in-depth analysis of 288 NDD genes shows there is highly significant protein-protein interaction (PPI) network connectivity, and functionally distinct PPI subnetworks based on pathway enrichment, single-cell RNA-seq cell types, and multi-region developmental brain RNA-seq.

CONCLUSIONS

We have extended a pipeline used in ASD studies and applied it to infer rare genetic parameters for SCZ and four NDDs ( https://github.com/hoangtn/extTADA ). We find many new DD risk genes, supported by gene set enrichment and PPI network connectivity analyses. We find greater similarity among NDDs than between NDDs and SCZ. NDD gene subnetworks are implicated in postnatally expressed presynaptic and postsynaptic genes, and for transcriptional and post-transcriptional gene regulation in prenatal neural progenitor and stem cells.

Recurrent de novo mutations in neurodevelopmental disorders: properties and clinical implications

Next-generation sequencing (NGS) is now more accessible to clinicians and researchers. As a result, our understanding of the genetics of neurodevelopmental disorders (NDDs) has rapidly advanced over the past few years. NGS has led to the discovery of new NDD genes with an excess of recurrent de novo mutations (DNMs) when compared to controls. Development of large-scale databases of normal and disease variation has given rise to metrics exploring the relative tolerance of individual genes to human mutation. Genetic etiology and diagnosis rates have improved, which have led to the discovery of new pathways and tissue types relevant to NDDs. In this review, we highlight several key findings based on the discovery of recurrent DNMs ranging from copy number variants to point mutations. We explore biases and patterns of DNM enrichment and the role of mosaicism and secondary mutations in variable expressivity. We discuss the benefit of whole-genome sequencing (WGS) over whole-exome sequencing (WES) to understand more complex, multifactorial cases of NDD and explain how this improved understanding aids diagnosis and management of these disorders. Comprehensive assessment of the DNM landscape across the genome using WGS and other technologies will lead to the development of novel functional and bioinformatics approaches to interpret DNMs and drive new insights into NDD biology.

Whole exome sequencing coupled with unbiased functional analysis reveals new Hirschsprung disease genes

Background

Hirschsprung disease (HSCR), which is congenital obstruction of the bowel, results from a failure of enteric nervous system (ENS) progenitors to migrate, proliferate, differentiate, or survive within the distal intestine. Previous studies that have searched for genes underlying HSCR have focused on ENS-related pathways and genes not fitting the current knowledge have thus often been ignored. We identify and validate novel HSCR genes using whole exome sequencing (WES), burden tests, in silico prediction, unbiased in vivo analyses of the mutated genes in zebrafish, and expression analyses in zebrafish, mouse, and human.

Results

We performed de novo mutation (DNM) screening on 24 HSCR trios. We identify 28 DNMs in 21 different genes. Eight of the DNMs we identified occur in RET, the main HSCR gene, and the remaining 20 DNMs reside in genes not reported in the ENS. Knockdown of all 12 genes with missense or loss-of-function DNMs showed that the orthologs of four genes (DENND3, NCLN, NUP98, and TBATA) are indispensable for ENS development in zebrafish, and these results were confirmed by CRISPR knockout. These genes are also expressed in human and mouse gut and/or ENS progenitors. Importantly, the encoded proteins are linked to neuronal processes shared by the central nervous system and the ENS.

Conclusions

Our data open new fields of investigation into HSCR pathology and provide novel insights into the development of the ENS. Moreover, the study demonstrates that functional analyses of genes carrying DNMs are warranted to delineate the full genetic architecture of rare complex diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1174-6) contains supplementary material, which is available to authorized users.

[Effects of father's age on fertility, results of ART and health of children]

Many studies exist on the impact of female age on fertility, success of assisted reproductive technologies and on obstetric, fetal and neonatal adverse outcomes. Late paternity seems commonplace especially in the media… But there are reliable scientific data which confirm decline of fertility related to male age but also an increased risk of genetic diseases for the offspring. The objective of this article is to make a synthesis of the literature on this subject.

Truncating de novo mutations in the Krüppel-type zinc-finger gene ZNF148 in patients with corpus callosum defects, developmental delay, short stature, and dysmorphisms

BACKGROUND

Krüppel-type zinc finger genes (ZNF) constitute a large yet relatively poorly characterized gene family. ZNF genes encode proteins that recognize specific DNA motifs in gene promotors. They act as transcriptional co-activators or -repressors via interaction with chromatin remodeling proteins and other transcription factors. Only few ZNF genes are currently linked to human disorders and identification of ZNF gene-associated human diseases may help understand their function. Here we provide genetic, statistical, and clinical evidence to support association of ZNF148 with a new intellectual disability (ID) syndrome disorder.

METHODS

Routine diagnostic exome sequencing data were obtained from 2172 patients with ID and/or multiple congenital anomalies.

RESULTS

In a cohort of 2172 patient-parent trios referred for routine diagnostic whole exome sequencing for ID and/or multiple congenital anomalies (MCA) in the period 2012-2016, four patients were identified who carried de novo heterozygous nonsense or frameshift mutations in the ZNF148 gene. This was the only ZNF gene with recurrent truncating de novo mutations in this cohort. All mutations resulted in premature termination codons in the last exon of ZNF148. The number of the de novo truncating mutations in the ZNF148 gene was significantly enriched (p = 5.42 × 10^-3). The newly described ZNF148-associated syndrome is characterized by underdevelopment of the corpus callosum, mild to moderate developmental delay and ID, variable microcephaly or mild macrocephaly, short stature, feeding problems, facial dysmorphisms, and cardiac and renal malformations.

CONCLUSIONS

We propose ZNF148 as a gene involved in a newly described ID syndrome with a recurrent phenotype and postulate that the ZNF148 is a hitherto unrecognized but crucial transcription factor in the development of the corpus callosum. Our study illustrates the advantage of whole exome sequencing in a large cohort using a parent-offspring trio approach for identifying novel genes involved in rare human diseases.

Effects of a Mutation in the HSPE1 Gene Encoding the Mitochondrial Co-chaperonin HSP10 and Its Potential Association with a Neurological and Developmental Disorder

We here report molecular investigations of a missense mutation in the HSPE1 gene encoding the HSP10 subunit of the HSP60/ HSP10 chaperonin complex that assists protein folding in the mitochondrial matrix. The mutation was identified in an infant who came to clinical attention due to infantile spasms at 3 months of age. Clinical exome sequencing revealed heterozygosity for a HSPE1 NM_002157.2:c.217C>T de novo mutation causing replacement of leucine with phenylalanine at position 73 of the HSP10 protein. This variation has never been observed in public exome sequencing databases or the literature. To evaluate whether the mutation may be disease-associated we investigated its effects by in vitro and ex vivo studies. Our in vitro studies indicated that the purified mutant protein was functional, yet its thermal stability, spontaneous refolding propensity, and resistance to proteolytic treatment were profoundly impaired. Mass spectrometric analysis of patient fibroblasts revealed barely detectable levels of HSP10-p.Leu73Phe protein resulting in an almost 2-fold decrease of the ratio of HSP10 to HSP60 subunits. Amounts of the mitochondrial superoxide dismutase SOD2, a protein whose folding is known to strongly depend on the HSP60/HSP10 complex, were decreased to approximately 20% in patient fibroblasts in spite of unchanged SOD2 transcript levels. As a likely consequence, mitochondrial superoxide levels were increased about 2-fold. Although, we cannot exclude other causative or contributing factors, our experimental data support the notion that the HSP10-p.Leu73Phe mutation could be the cause or a strong contributing factor for the disorder in the described patient.

VARPRISM: incorporating variant prioritization in tests of de novo mutation association

Background

Patients with certain genetic diseases, such as autism spectrum disorder, have increased rates of de novo mutations within some protein-coding genes.

Results

We introduce the VARiant PRIoritization SuM (VARPRISM), a software package which incorporates functional variant prioritization information to improve the power to detect de novo mutations influencing disease risk. VARPRISM evaluates the consequence of any given exonic mutation on the protein sequence to estimate the likelihood that the mutation is benign or damaging and conducts a likelihood ratio test on the gene level. We analyzed the Simons Simplex Collection of 2508 parent-offspring autism trios using VARPRISM, replicating 44 genes previously implicated in autism susceptibility and identifying 20 additional candidate genes, including MYO1E, KCND3, PDCD1, DLX3, and TSPAN4 (false discovery rate < 0.3).

Conclusion

By incorporating functional predictions, VARPRISM improved the statistical power to identify de novo mutations increasing disease risks. VARPRISM is available at http://www.hufflab.org/software/VARPRISM.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0341-9) contains supplementary material, which is available to authorized users.

High-specificity detection of rare alleles with Paired-End Low Error Sequencing (PELE-Seq)

BACKGROUND

Polymorphic loci exist throughout the genomes of a population and provide the raw genetic material needed for a species to adapt to changes in the environment. The minor allele frequencies of rare Single Nucleotide Polymorphisms (SNPs) within a population have been difficult to track with Next-Generation Sequencing (NGS), due to the high error rate of standard methods such as Illumina sequencing.

RESULTS

We have developed a wet-lab protocol and variant-calling method that identifies both sequencing and PCR errors, called Paired-End Low Error Sequencing (PELE-Seq). To test the specificity and sensitivity of the PELE-Seq method, we sequenced control E. coli DNA libraries containing known rare alleles present at frequencies ranging from 0.2-0.4 % of the total reads. PELE-Seq had higher specificity and sensitivity than standard libraries. We then used PELE-Seq to characterize rare alleles in a Caenorhabditis remanei nematode worm population before and after laboratory adaptation, and found that minor and rare alleles can undergo large changes in frequency during lab-adaptation.

CONCLUSION

We have developed a method of rare allele detection that mitigates both sequencing and PCR errors, called PELE-Seq. PELE-Seq was evaluated using control E. coli populations and was then used to compare a wild C. remanei population to a lab-adapted population. The PELE-Seq method is ideal for investigating the dynamics of rare alleles in a broad range of reduced-representation sequencing methods, including targeted amplicon sequencing, RAD-Seq, ddRAD, and GBS. PELE-Seq is also well-suited for whole genome sequencing of mitochondria and viruses, and for high-throughput rare mutation screens.

Evolutionary Genetics of the Cavefish Astyanax mexicanus

Blind and depigmented fish belonging to the species Astyanax mexicanus are outstanding models for evolutionary genetics. During their evolution in the darkness of caves, they have undergone a number of changes at the morphological, physiological, and behavioral levels, but they can still breed with their river-dwelling conspecifics. The fertile hybrids between these two morphotypes allow forward genetic approaches, from the search of quantitative trait loci to the identification of the mutations underlying the evolution of troglomorphism. We review here the past 30years of evolutionary genetics on Astyanax: from the first crosses and the discovery of convergent evolution of different Astyanax cavefish populations to the most recent evolutionary transcriptomics and genomics studies that have provided researchers with potential candidate genes to be tested using functional genetic approaches. Although significant progress has been made and some genes have been identified, cavefish have not yet fully revealed the secret of their adaptation to the absence of light. In particular, the genetic determinism of their loss of eyes seems complex and still puzzles researchers. We also discuss future research directions, including searches for the origin of cave alleles and searches for selection genome-wide, as well as the necessary but missing information on the timing of cave colonization by surface fish.

The importance of de novo mutations for pediatric neurological disease--It is not all in utero or birth trauma

The advent of next generation sequencing (NGS, which consists of massively parallel sequencing to perform TGS (total genome sequencing) or WES (whole exome sequencing)) has abundantly discovered many causative mutations in patients with pediatric neurological disease. A surprisingly high number of these are de novo mutations which have not been inherited from either parent. For epilepsy, autism spectrum disorders, and neuromotor disorders, including cerebral palsy, initial estimates put the frequency of causative de novo mutations at about 15% and about 10% of these are somatic. There are some shared mutated genes between these three classes of disease. Studies of copy number variation by comparative genomic hybridization (CGH) proceded the NGS approaches but they also detect de novo variation which is especially important for ASDs. There are interesting differences between the mutated genes detected by CGS and NGS. In summary, de novo mutations cause a very significant proportion of pediatric neurological disease.

Paternal and maternal ages have contrasting associations with self-reported schizophrenia liability

BACKGROUND

Older paternal age predicts schizophrenia diagnosis in offspring. If this relationship reflects a pathogenic process, paternal age should predict the expression of subclinical schizophrenia liability (schizotypy). We hypothesized that paternal and maternal ages predict positive, negative, and disorganized features of schizotypy, that family history of psychosis moderates the relationship of paternal age with schizotypy, and that stress sensitivity mediates the relationship of maternal age with schizotypy.

METHOD

Two studies are reported, each of undergraduates (n=500 and n=211) who completed the Schizotypal Personality Questionnaire. The second was designed to replicate and extend the first and included assessment of stress sensitivity.

RESULTS

In Study 1, older paternal age and younger maternal age predicted greater positive schizotypy (β=.13 and β=-.19, respectively). Parental ages did not predict negative or disorganized features and family history did not moderate the paternal age association. In Study 2, the same pattern of associations between parental ages and schizotypy components was observed. Additionally, stress sensitivity partially mediated the association of maternal age with positive schizotypy whereas it did not contribute to the paternal age association.

CONCLUSION

The association between older paternal age and schizophrenia extends to self-reported positive features of schizophrenia liability, consistent with the notion that this relationship arises from a pathogenic process, such as de novo mutations. Importantly, younger maternal age was an equally potent predictor of positive schizotypy, with its association partially mediated by stress sensitivity.

De novo missense variants in PPP2R5D are associated with intellectual disability, macrocephaly, hypotonia, and autism

Protein phosphatase 2A (PP2A) is a heterotrimeric protein serine/threonine phosphatase and is involved in a broad range of cellular processes. PPP2R5D is a regulatory B subunit of PP2A and plays an important role in regulating key neuronal and developmental regulation processes such as PI3K/AKT and glycogen synthase kinase 3 beta (GSK3β)-mediated cell growth, chromatin remodeling, and gene transcriptional regulation. Using whole-exome sequencing (WES), we identified four de novo variants in PPP2R5D in a total of seven unrelated individuals with intellectual disability (ID) and other shared clinical characteristics, including autism spectrum disorder, macrocephaly, hypotonia, seizures, and dysmorphic features. Among the four variants, two have been previously reported and two are novel. All four amino acids are highly conserved among the PP2A subunit family, and all change a negatively charged acidic glutamic acid (E) to a positively charged basic lysine (K) and are predicted to disrupt the PP2A subunit binding and impair the dephosphorylation capacity. Our data provides further support for PPP2R5D as a genetic cause of ID.

Paternal age of schizophrenia probands and endophenotypic differences from unaffected siblings

We evaluated the discrepancy of endophenotypic performance between probands with schizophrenia and unaffected siblings by paternal age at proband birth, a possible marker for de novo mutations. Pairs of schizophrenia probands and unaffected siblings (N=220 pairs) were evaluated on 11 neuropsychological or neurophysiological endophenotypes previously identified as heritable. For each endophenotype, the sibling-minus-proband differences were transformed to standardized scores. Then for each pair, the average discrepancy was calculated from its standardized scores. We tested the hypothesis that the discrepancy is associated with paternal age, controlling for the number of endophenotypes shared between proband and his or her sibling, and proband age, which were both associated with paternal age. The non-significant association between the discrepancy and paternal age was in the opposite direction from the hypothesis. Of the 11 endophenotypes only sensori-motor dexterity was significant, but in the opposite direction. Eight other endophenotypes were also in the opposite direction, but not significant. The results did not support the hypothesized association of increased differences between sibling/proband pairs with greater paternal age. A possible explanation is that the identification of heritable endophenotypes was based on samples for which schizophrenia was attributable to inherited rather than de novo/non-inherited causes.