Genome-wide characteristics of de novo mutations in autism

Proteome-Wide Assessment of Clustering of Missense Variants in Neurodevelopmental Disorders Versus Cancer

Missense de novo variants (DNVs) and missense somatic variants contribute to neurodevelopmental disorders (NDDs) and cancer, respectively. Proteins with statistical enrichment based on analyses of these variants exhibit convergence in the differing NDD and cancer phenotypes. Herein, the question of why some of the same proteins are identified in both phenotypes is examined through investigation of clustering of missense variation at the protein level. Our hypothesis is that missense variation is present in different protein locations in the two phenotypes leading to the distinct phenotypic outcomes. We tested this hypothesis in 1D protein space using our software CLUMP. Furthermore, we newly developed 3D-CLUMP that uses 3D protein structures to spatially test clustering of missense variation for proteome-wide significance. We examined missense DNVs in 39,883 parent-child sequenced trios with NDDs and missense somatic variants from 10,543 sequenced tumors covering five TCGA cancer types and two COSMIC pan-cancer aggregates of tissue types. There were 57 proteins with proteome-wide significant missense variation clustering in NDDs when compared to cancers and 79 proteins with proteome-wide significant missense clustering in cancers compared to NDDs. While our main objective was to identify differences in patterns of missense variation, we also identified a novel NDD protein BLTP2. Overall, our study is innovative, provides new insights into differential missense variation in NDDs and cancer at the protein-level, and contributes necessary information toward building a framework for thinking about prognostic and therapeutic aspects of these proteins.

Intellectual Disability and Behavioral Deficits Linked to CYFIP1 Missense Variants Disrupting Actin Polymerization

BACKGROUND

15q11.2 deletions and duplications have been linked to autism spectrum disorder, schizophrenia, and intellectual disability. Recent evidence suggests that dysfunctional CYFIP1 (cytoplasmic FMR1 interacting protein 1) contributes to the clinical phenotypes observed in individuals with 15q11.2 deletion/duplication syndrome. CYFIP1 plays crucial roles in neuronal development and brain connectivity, promoting actin polymerization and regulating local protein synthesis. However, information about the impact of single nucleotide variants in CYFIP1 on neurodevelopmental disorders is limited.

METHODS

Here, we report a family with 2 probands exhibiting intellectual disability, autism spectrum disorder, spastic tetraparesis, and brain morphology defects and who carry biallelic missense point mutations in the CYFIP1 gene. We used skin fibroblasts from one of the probands, the parents, and typically developing individuals to investigate the effect of the variants on the functionality of CYFIP1. In addition, we generated Drosophila knockin mutants to address the effect of the variants in vivo and gain insight into the molecular mechanism that underlies the clinical phenotype.

RESULTS

Our study revealed that the 2 missense variants are in protein domains responsible for maintaining the interaction within the wave regulatory complex. Molecular and cellular analyses in skin fibroblasts from one proband showed deficits in actin polymerization. The fly model for these mutations exhibited abnormal brain morphology and F-actin loss and recapitulated the core behavioral symptoms, such as deficits in social interaction and motor coordination.

CONCLUSIONS

Our findings suggest that the 2 CYFIP1 variants contribute to the clinical phenotype in the probands that reflects deficits in actin-mediated brain development processes.

Comprehensive systematic review and meta-analysis of the association between common genetic variants and autism spectrum disorder

BACKGROUND

Autism spectrum disorder (ASD) is neurodevelopmental disorder characterized by stereotyped behavior and deficits in communication and social interactions. To date, numerous studies have investigated the associations between genetic variants and ASD risk. However, the results of these published studies lack a clear consensus. In the present study, we performed a systematic review on the association between genetic variants and ASD risk. Meanwhile, we conducted a meta-analysis on available data to identify the association between the single nucleotide polymorphisms (SNPs) of candidate genes and ASD risk.

METHODS

We systematically searched public databases including English and Chinese from their inception to August 1, 2022. Two independent reviewers extracted data and assessed study quality. Odds ratio and 95 % confidence interval were used as effect indexes to evaluate the association between the SNPs of candidate genes and the risk of ASD. Heterogeneity was explored through subgroup, sensitivity, and meta-regression analyses. Publication bias was assessed by using Egger's and Begg's tests for funnel plot asymmetry. In addition, TSA analysis were performed to confirm the study findings.

RESULTS

We summarized 84 SNPs of 32 candidate genes from 81 articles included in the study. Subsequently, we analyzed 16 SNPs of eight genes by calculating pooled ORs, and identified eight significant SNPs of contactin associated protein 2 (CNTNAP2), methylentetrahydrofolate reductase (MTHFR), oxytocin receptor (OXTR), and vitamin D receptor (VDR). Results showed that seven SNPs, including the CNTNAP2 rs2710102 (homozygote, heterozygote, dominant and allelic models) and rs7794745 (heterozygote and dominant models), MTHFR C677T (homozygote, heterozygote, dominant, recessive and allelic models) and A1298C (dominant and allelic models), OXTR rs2254298 (homozygote and recessive models), VDR rs731236 (homozygote, dominant, recessive and allelic models) and rs2228570 (homozygote and recessive models), were showed to be correlated with an increased ASD risk. By contrast, the VDR rs7975232 was correlated with a decreased the risk of ASD under the homozygote and allelic models.

CONCLUSION

Our study summarized research evidence on the genetic variants of ASD and provides a broad and detailed overview of ASD risk genes. The C677T and A1298C polymorphisms of MTHFR, rs2710102 and rs7794745 polymorphisms of CNTNAP2, rs2254298 polymorphism of OXTR, and rs731236 and rs2228570 polymorphisms of VDR were genetic risk factors. The rs7975232 polymorphism of VDR was a genetic protective factor for ASD. Our study provides novel clues to clinicians and healthcare decision-makers to predict ASD susceptibility.

The Phenotypic variability of 16p11.2 distal BP2-BP3 deletion in a transgenerational family and in neurodevelopmentally ascertained samples

BACKGROUND

We present genomic and phenotypic findings of a transgenerational family consisting of three male offspring, each with a maternally inherited distal 220 kb deletion at locus 16p11.2 (BP2-BP3). Genomic analysis of all family members was prompted by a diagnosis of autism spectrum disorder (ASD) in the eldest child, who also presented with a low body mass index.

METHODS

All male offspring underwent extensive neuropsychiatric evaluation. Both parents were also assessed for social functioning and cognition. The family underwent whole-genome sequencing. Further data curation was undertaken from samples ascertained for neurodevelopmental disorders and congenital abnormalities.

RESULTS

On medical examination, both the second and third-born male offspring presented with obesity. The second-born male offspring met research diagnostic criteria for ASD at 8 years of age and presented with mild attention deficits. The third-born male offspring was only noted as having motor deficits and received a diagnosis of developmental coordination disorder. Other than the 16p11.2 distal deletion, no additional contributing variants of clinical significance were observed. The mother was clinically evaluated and noted as having a broader autism phenotype.

CONCLUSION

In this family, the phenotypes observed are most likely caused by the 16p11.2 distal deletion. The lack of other overt pathogenic mutations identified by genomic sequencing reinforces the variable expressivity that should be heeded in a clinical setting. Importantly, distal 16p11.2 deletions can present with a highly variable phenotype even within a single family. Our additional data curation provides further evidence on the variable clinical presentation among those with pathogenetic 16p11.2 (BP2-BP3) mutations.

De novo network analysis reveals autism causal genes and developmental links to co-occurring traits

Autism is a complex neurodevelopmental condition that manifests in various ways. Autism is often accompanied by other conditions, such as attention-deficit/hyperactivity disorder and schizophrenia, which can complicate diagnosis and management. Although research has investigated the role of specific genes in autism, their relationship with co-occurring traits is not fully understood. To address this, we conducted a two-sample Mendelian randomisation analysis and identified four genes located at the 17q21.31 locus that are putatively causal for autism in fetal cortical tissue (LINC02210, LRRC37A4P, RP11-259G18.1, and RP11-798G7.6). LINC02210 was also identified as putatively causal for autism in adult cortical tissue. By integrating data from expression quantitative trait loci, genes and protein interactions, we identified that the 17q21.31 locus contributes to the intersection between autism and other neurological traits in fetal cortical tissue. We also identified a distinct cluster of co-occurring traits, including cognition and worry, linked to the genetic loci at 3p21.1. Our findings provide insights into the relationship between autism and co-occurring traits, which could be used to develop predictive models for more accurate diagnosis and better clinical management.

Three generation families: Analysis of de novo variants in autism

De novo variants (DNVs) analysis has proven to be a powerful approach to gene discovery in Autism Spectrum Disorder (ASD), which has not yet been shown in a Brazilian ASD cohort. The relevance of inherited rare variants has also been suggested, particularly in oligogenic models. We hypothesized that three-generation analyses of DNVs could provide new insights into the relevance of de novo and inherited variants across generations. To accomplish this goal, we performed whole-exome sequencing of 33 septet families composed of probands, parents, and grandparents (n = 231 individuals) and compared DNV rates (DNVr) between generations and those from two control cohorts. The DNVr in the probands (DNVr = 1.16) was marginally higher than in parents (DNVr = 0.60; p = 0.054), and in controls (DNVr = 0.68; p = 0.035, congenital heart disorder and DNVr = 0.70; p = 0.047, unaffected ASD siblings from Simons Simplex Collection). Moreover, most of the DNVs were found to have paternal origin in both generations (84.6%). Finally, we observed that 40% (6/15) of the DNVs in parents transmitted for probands are in ASD or ASD candidate genes, representing recently emerged risk variants to ASD in their families and suggest ZNF536, MSL2 and HDAC9 as ASD candidate genes. We did not observe an enrichment of risk variants nor sex bias of transmitted variants in the three generations, that can be due to sample size. These results further reinforce the relevance of de novo variants in ASD.

The landscape of somatic mutations in lymphoblastoid cell lines

Somatic mutations have important biological ramifications while exerting substantial rate, type, and genomic location heterogeneity. Yet, their sporadic occurrence makes them difficult to study at scale and across individuals. Lymphoblastoid cell lines (LCLs), a model system for human population and functional genomics, harbor large numbers of somatic mutations and have been extensively genotyped. By comparing 1,662 LCLs, we report that the mutational landscape of the genome varies across individuals in terms of the number of mutations, their genomic locations, and their spectra; this variation may itself be modulated by somatic trans-acting mutations. Mutations attributed to the translesion DNA polymerase η follow two different modes of formation, with one mode accounting for the hypermutability of the inactive X chromosome. Nonetheless, the distribution of mutations along the inactive X chromosome appears to follow an epigenetic memory of the active form.

Revisiting mutagenesis at non-B DNA motifs in the human genome

Non-B DNA structures formed by repetitive sequence motifs are known instigators of mutagenesis in experimental systems. Analyzing this phenomenon computationally in the human genome requires careful disentangling of intrinsic confounding factors, including overlapping and interrupted motifs and recurrent sequencing errors. Here, we show that accounting for these factors eliminates all signals of repeat-induced mutagenesis that extend beyond the motif boundary, and eliminates or dramatically shrinks the magnitude of mutagenesis within some motifs, contradicting previous reports. Mutagenesis not attributable to artifacts revealed several biological mechanisms. Polymerase slippage generates frequent indels within every variety of short tandem repeat motif, implicating slipped-strand structures. Interruption-correcting single nucleotide variants within short tandem repeats may originate from error-prone polymerases. Secondary-structure formation promotes single nucleotide variants within palindromic repeats and duplications within direct repeats. G-quadruplex motifs cause recurrent sequencing errors, whereas mutagenesis at Z-DNAs is conspicuously absent.

Evolution of the germline mutation rate across vertebrates

The germline mutation rate determines the pace of genome evolution and is an evolving parameter itself¹. However, little is known about what determines its evolution, as most studies of mutation rates have focused on single species with different methodologies². Here we quantify germline mutation rates across vertebrates by sequencing and comparing the high-coverage genomes of 151 parent-offspring trios from 68 species of mammals, fishes, birds and reptiles. We show that the per-generation mutation rate varies among species by a factor of 40, with mutation rates being higher for males than for females in mammals and birds, but not in reptiles and fishes. The generation time, age at maturity and species-level fecundity are the key life-history traits affecting this variation among species. Furthermore, species with higher long-term effective population sizes tend to have lower mutation rates per generation, providing support for the drift barrier hypothesis³. The exceptionally high yearly mutation rates of domesticated animals, which have been continually selected on fecundity traits including shorter generation times, further support the importance of generation time in the evolution of mutation rates. Overall, our comparative analysis of pedigree-based mutation rates provides ecological insights on the mutation rate evolution in vertebrates.

SCIP: software for efficient clinical interpretation of copy number variants detected by whole-genome sequencing

Copy number variants (CNVs) represent major etiologic factors in rare genetic diseases. Current clinical CNV interpretation workflows require extensive back-and-forth with multiple tools and databases. This increases complexity and time burden, potentially resulting in missed genetic diagnoses. We present the Suite for CNV Interpretation and Prioritization (SCIP), a software package for the clinical interpretation of CNVs detected by whole-genome sequencing (WGS). The SCIP Visualization Module near-instantaneously displays all information necessary for CNV interpretation (variant quality, population frequency, inheritance pattern, and clinical relevance) on a single page-supported by modules providing variant filtration and prioritization. SCIP was comprehensively evaluated using WGS data from 1027 families with congenital cardiac disease and/or autism spectrum disorder, containing 187 pathogenic or likely pathogenic (P/LP) CNVs identified in previous curations. SCIP was efficient in filtration and prioritization: a median of just two CNVs per case were selected for review, yet it captured all P/LP findings (92.5% of which ranked 1st). SCIP was also able to identify one pathogenic CNV previously missed. SCIP was benchmarked against AnnotSV and a spreadsheet-based manual workflow and performed superiorly than both. In conclusion, SCIP is a novel software package for efficient clinical CNV interpretation, substantially faster and more accurate than previous tools (available at https://github.com/qd29/SCIP , a video tutorial series is available at https://bit.ly/SCIPVideos ).

Dendritic Spine in Autism Genetics: Whole-Exome Sequencing Identifying De Novo Variant of CTTNBP2 in a Quad Family Affected by Autism Spectrum Disorder

Autism spectrum disorder (ASD) affects around 1% of children with no effective blood test or cure. Recent studies have suggested that these are neurological disorders with a strong genetic basis and that they are associated with the abnormal formation of dendritic spines. Chromosome microarray (CMA) together with high-throughput sequencing technology has been used as a powerful tool to identify new candidate genes for ASD. In the present study, CMA was first used to scan for genome-wide copy number variants in a proband, and no clinically significant copy number variants were found. Whole-exome sequencing (WES) was used further for genetic testing of the whole quad family affected by ASD, including the proband, his non-autistic sister, and his parents. Sanger sequencing and MassARRAY-based validation were used to identify and confirm variants associated with ASD. WES yielded a 151-fold coverage depth for each sample. A total of 98.65% of the targeted whole-exome region was covered at >20-fold depth. A de novo variant in CTTNBP2, p.M115T, was identified. The CTTNBP2 gene belongs to a family of ankyrin repeat domain-containing proteins associated with dendritic spine formation. Although CTTNBP2 has been associated with ASD, limited studies have been developed to identify clinically relevant de novo mutations of CTTNBP2 in children with ASD; family-based WES successfully identified a clinically relevant mutation in the CTTNBP2 gene in a quad family affected by ASD. Considering the neuron-specific expression of CTTNBP2 and its role in dendritic spine formation, our results suggest a correlation between the CTTNBP2 mutation and ASD, providing genetic evidence for ASD spine pathology. Although the present study is currently insufficient to support the assertion that the de novo mutation M115T in CTTNBP2 directly causes the autism phenotype, our study provides support for the assertion that this mutation is a candidate clinically relevant variant in autism.

Fertility in the aging male: a systematic review

The negative effect of advanced female age on fertility and offspring health is well understood. In comparison, much less is known about the implications of male age on fertility, with many studies showing conflicting results. Nevertheless, increasing evidence suggests that advanced paternal age has negative effects on sperm parameters, reproductive success, and offspring health. Herein, we summarize the current body of knowledge on this controversial topic, with the belief that this review will serve as a resource for the clinicians providing fertility counseling to couples with older male partners.

Thyroid hormone elicits intergenerational epigenetic effects on adult social behavior and fetal brain expression of autism susceptibility genes

Genetic mutations identified in genome-wide association studies can only explain a small percentage of the cases of complex, highly heritable human conditions, including neurological and neurodevelopmental disorders. This suggests that intergenerational epigenetic effects, possibly triggered by environmental circumstances, may contribute to their etiology. We previously described altered DNA methylation signatures in the sperm of mice that experienced developmental overexposure to thyroid hormones as a result of a genetic defect in hormone clearance (DIO3 deficiency). Here we studied fetal brain gene expression and adult social behavior in genetically normal F2 generation descendants of overexposed mice. The brain of F2 generation E13.5 fetuses exhibited abnormal expression of genes associated with autism in humans, including Auts2, Disc1, Ldlr, Per2, Shank3, Oxtr, Igf1, Foxg1, Cd38, Grid2, Nrxn3, and Reln. These abnormal gene expression profiles differed depending on the sex of the exposed ancestor. In the three-chamber social box test, adult F2 generation males manifested significantly decreased interest in social interaction and social novelty, as revealed by decrease total time, distance traveled and time immobile in the area of interaction with novel strangers. F1 generation mice, compared to appropriate controls also exhibited altered profiles in fetal brain gene expression, although these profiles were substantially different to those in the F2 generation. Likewise adult F1 generation mice showed some abnormalities in social behavior that were sexually dimorphic and milder than those in F2 generation mice. Our results indicate that developmental overexposure to thyroid hormone causes intergenerational epigenetic effects impacting social behavior and the expression of autism-related genes during early brain development. Our results open the possibility that altered thyroid hormone states, by eliciting changes in the epigenetic information of the germ line, contribute to the susceptibility and the missing-but heriTables-etiology of complex neurodevelopmental conditions characterized by social deficits, including autism and schizophrenia.

Deciphering the impact of genetic variation on human polyadenylation using APARENT2

BACKGROUND

3'-end processing by cleavage and polyadenylation is an important and finely tuned regulatory process during mRNA maturation. Numerous genetic variants are known to cause or contribute to human disorders by disrupting the cis-regulatory code of polyadenylation signals. Yet, due to the complexity of this code, variant interpretation remains challenging.

RESULTS

We introduce a residual neural network model, APARENT2, that can infer 3'-cleavage and polyadenylation from DNA sequence more accurately than any previous model. This model generalizes to the case of alternative polyadenylation (APA) for a variable number of polyadenylation signals. We demonstrate APARENT2's performance on several variant datasets, including functional reporter data and human 3' aQTLs from GTEx. We apply neural network interpretation methods to gain insights into disrupted or protective higher-order features of polyadenylation. We fine-tune APARENT2 on human tissue-resolved transcriptomic data to elucidate tissue-specific variant effects. By combining APARENT2 with models of mRNA stability, we extend aQTL effect size predictions to the entire 3' untranslated region. Finally, we perform in silico saturation mutagenesis of all human polyadenylation signals and compare the predicted effects of [Formula: see text] million variants against gnomAD. While loss-of-function variants were generally selected against, we also find specific clinical conditions linked to gain-of-function mutations. For example, we detect an association between gain-of-function mutations in the 3'-end and autism spectrum disorder. To experimentally validate APARENT2's predictions, we assayed clinically relevant variants in multiple cell lines, including microglia-derived cells.

CONCLUSIONS

A sequence-to-function model based on deep residual learning enables accurate functional interpretation of genetic variants in polyadenylation signals and, when coupled with large human variation databases, elucidates the link between functional 3'-end mutations and human health.

Phasing of de novo mutations using a scaled-up multiple amplicon long-read sequencing approach

De novo mutations (DNMs) play an important role in severe genetic disorders that reduce fitness. To better understand their role in disease, it is important to determine the parent-of-origin and timing of mutational events that give rise to these mutations, especially in sex-specific developmental disorders such as male infertility. However, currently available short-read sequencing approaches are not ideally suited for phasing, as this requires long continuous DNA strands that span both the DNM and one or more informative single-nucleotide polymorphisms. To overcome these challenges, we optimized and implemented a multiplexed long-read sequencing approach using Oxford Nanopore technologies MinION platform. We focused on improving target amplification, integrating long-read sequenced data with high-quality short-read sequence data, and developing an anchored phasing computational method. This approach handled the inherent phasing challenges of long-range target amplification and the normal accumulation of sequencing error associated with long-read sequencing. In total, 77 of 109 DNMs (71%) were successfully phased and parent-of-origin identified. The majority of phased DNMs were prezygotic (90%), the accuracy of which is highlighted by an average mutant allele frequency of 49.6% and standard error of 0.84%. This study demonstrates the benefits of employing an integrated short-read and long-read sequencing approach for large-scale DNM phasing.

Non-coding de novo mutations in chromatin interactions are implicated in autism spectrum disorder

Three-dimensional chromatin interactions regulate gene expressions. The significance of de novo mutations (DNMs) in chromatin interactions remains poorly understood for autism spectrum disorder (ASD). We generated 813 whole-genome sequences from 242 Korean simplex families to detect DNMs, and identified target genes which were putatively affected by non-coding DNMs in chromatin interactions. Non-coding DNMs in chromatin interactions were significantly involved in transcriptional dysregulations related to ASD risk. Correspondingly, target genes showed spatiotemporal expressions relevant to ASD in developing brains and enrichment in biological pathways implicated in ASD, such as histone modification. Regarding clinical features of ASD, non-coding DNMs in chromatin interactions particularly contributed to low intelligence quotient levels in ASD probands. We further validated our findings using two replication cohorts, Simons Simplex Collection (SSC) and MSSNG, and showed the consistent enrichment of non-coding DNM-disrupted chromatin interactions in ASD probands. Generating human induced pluripotent stem cells in two ASD families, we were able to demonstrate that non-coding DNMs in chromatin interactions alter the expression of target genes at the stage of early neural development. Taken together, our findings indicate that non-coding DNMs in ASD probands lead to early neurodevelopmental disruption implicated in ASD risk via chromatin interactions.

Genome-wide rare variant score associates with morphological subtypes of autism spectrum disorder

Defining different genetic subtypes of autism spectrum disorder (ASD) can enable the prediction of developmental outcomes. Based on minor physical and major congenital anomalies, we categorize 325 Canadian children with ASD into dysmorphic and nondysmorphic subgroups. We develop a method for calculating a patient-level, genome-wide rare variant score (GRVS) from whole-genome sequencing (WGS) data. GRVS is a sum of the number of variants in morphology-associated coding and non-coding regions, weighted by their effect sizes. Probands with dysmorphic ASD have a significantly higher GRVS compared to those with nondysmorphic ASD (P = 0.03). Using the polygenic transmission disequilibrium test, we observe an over-transmission of ASD-associated common variants in nondysmorphic ASD probands (P = 2.9 × 10^-3). These findings replicate using WGS data from 442 ASD probands with accompanying morphology data from the Simons Simplex Collection. Our results provide support for an alternative genomic classification of ASD subgroups using morphology data, which may inform intervention protocols.

Predictive functional, statistical and structural analysis of CSNK2A1 and CSNK2B variants linked to neurodevelopmental diseases

Okur-Chung Neurodevelopmental Syndrome (OCNDS) and Poirier-Bienvenu Neurodevelopmental Syndrome (POBINDS) were recently identified as rare neurodevelopmental disorders. OCNDS and POBINDS are associated with heterozygous mutations in the CSNK2A1 and CSNK2B genes which encode CK2α, a serine/threonine protein kinase, and CK2β, a regulatory protein, respectively, which together can form a tetrameric enzyme called protein kinase CK2. A challenge in OCNDS and POBINDS is to understand the genetic basis of these diseases and the effect of the various CK2⍺ and CK2β mutations. In this study we have collected all variants available to date in CSNK2A1 and CSNK2B, and identified hotspots. We have investigated CK2⍺ and CK2β missense mutations through prediction programs which consider the evolutionary conservation, functionality and structure or these two proteins, compared these results with published experimental data on CK2α and CK2β mutants, and suggested prediction programs that could help predict changes in functionality of CK2α mutants. We also investigated the potential effect of CK2α and CK2β mutations on the 3D structure of the proteins and in their binding to each other. These results indicate that there are functional and structural consequences of mutation of CK2α and CK2β, and provide a rationale for further study of OCNDS and POBINDS-associated mutations. These data contribute to understanding the genetic and functional basis of these diseases, which is needed to identify their underlying mechanisms.

Krüppel-like factor 7 deficiency causes autistic-like behavior in mice via regulating Clock gene

BACKGROUND

Krüppel-like factor 7 (klf7), a transcription factor in the nervous system to regulate cell proliferation and differentiation, has been recently identified as a causal gene for autism spectrum disorder (ASD), but the mechanism behind remains unknown.

RESULT

To uncover this mechanism, in this study we characterized the involvement of klf7 in circadian rhythm by knocking down klf7 in N2A cells and examining the rhythmic expression of circadian genes, especially Clock gene. We constructed klf7^-/- mice and then investigated into klf7 regulation on the expression of rhythm genes in vivo as well as the use of melatonin to rescue the autism behavior. Our results illustrated that circadian rhythm was disrupted in klf7 knockdown cells and that klf7^-/- mice showed autism-like behavior. Also, we found that Clock gene was downregulated in the brain of these klf7^-/- mice and that the downstream rhythm genes of Clock were disturbed. Melatonin, as a circadian regulation drug, could regulate the expression level and amplitude of rhythm genes in klf7 knockout cells and further rescue the autistic behavior of klf7^-/- mice.

CONCLUSION

Klf7 deficiency causes ASD by disrupting circadian rhythm related genes to trigger rhythm oscillations. To treat ASD, maintaining circadian homeostasis is promising with the use of melatonin.

A rigorous in silico genomic interrogation at 1p13.3 reveals 16 autosomal dominant candidate genes in syndromic neurodevelopmental disorders

Genome-wide chromosomal microarray is extensively used to detect copy number variations (CNVs), which can diagnose microdeletion and microduplication syndromes. These small unbalanced chromosomal structural rearrangements ranging from 1 kb to 10 Mb comprise up to 15% of human mutations leading to monogenic or contiguous genomic disorders. Albeit rare, CNVs at 1p13.3 cause a variety of neurodevelopmental disorders (NDDs) including development delay (DD), intellectual disability (ID), autism, epilepsy, and craniofacial anomalies (CFA). Most of the 1p13.3 CNV cases reported in the pre-microarray era encompassed a large number of genes and lacked the demarcating genomic coordinates, hampering the discovery of positional candidate genes within the boundaries. In this study, we present four subjects with 1p13.3 microdeletions displaying DD, ID, autism, epilepsy, and CFA. In silico comparative genomic mapping with three previously reported subjects with CNVs and 22 unreported DECIPHER CNV cases has resulted in the identification of four different sub-genomic loci harboring five positional candidate genes for DD, ID, and CFA at 1p13.3. Most of these genes have pathogenic variants reported, and their interacting genes are involved in NDDs. RT-qPCR in various human tissues revealed a high expression pattern in the brain and fetal brain, supporting their functional roles in NDDs. Interrogation of variant databases and interacting protein partners led to the identification of another set of 11 potential candidate genes, which might have been dysregulated by the position effect of these CNVs at 1p13.3. Our studies define 1p13.3 as a genomic region harboring 16 NDD candidate genes and underscore the critical roles of small CNVs in in silico comparative genomic mapping for disease gene discovery. Our candidate genes will help accelerate the isolation of pathogenic heterozygous variants from exome/genome sequencing (ES/GS) databases.

De novo mutations, genetic mosaicism and human disease

Mosaicism—the existence of genetically distinct populations of cells in a particular organism—is an important cause of genetic disease. Mosaicism can appear as de novo DNA mutations, epigenetic alterations of DNA, and chromosomal abnormalities. Neurodevelopmental or neuropsychiatric diseases, including autism—often arise by de novo mutations that usually not present in either of the parents. De novo mutations might occur as early as in the parental germline, during embryonic, fetal development, and/or post-natally, through ageing and life. Mutation timing could lead to mutation burden of less than heterozygosity to approaching homozygosity. Developmental timing of somatic mutation attainment will affect the mutation load and distribution throughout the body. In this review, we discuss the timing of de novo mutations, spanning from mutations in the germ lineage (all ages), to post-zygotic, embryonic, fetal, and post-natal events, through aging to death. These factors can determine the tissue specific distribution and load of de novo mutations, which can affect disease. The disease threshold burden of somatic de novo mutations of a particular gene in any tissue will be important to define.

Contribution of copy number variations to the risk of severe eating disorders

AIM

Eating disorders (EDs) are complex, multifactorial psychiatric conditions. Previous studies identified pathogenic copy number variations associated with NDDs (NDD-CNVs) in ED patients. However, no statistical evidence for an association between NDD-CNVs and EDs has been demonstrated. Therefore, we examined whether NDD-CNVs confer risk for EDs.

METHODS

Using array comparative genomic hybridization (aCGH), we conducted a high-resolution CNV analysis of 71 severe female ED patients and 1045 female controls. According to the American College of Medical Genetics guidelines, we identified NDD-CNVs or pathogenic/likely pathogenic CNVs in NDD-linked loci. Gene set analysis was performed to examine the involvement of synaptic dysfunction in EDs. Clinical data were retrospectively examined for ED patients with NDD-CNVs.

RESULTS

Of the samples analyzed with aCGH, 70 severe ED patients (98.6%) and 1036 controls (99.1%) passed our quality control filtering. We obtained 189 and 2539 rare CNVs from patients and controls, respectively. NDD-CNVs were identified in 10.0% (7/70) of patients and 2.3% (24/1036) of controls. Statistical analysis revealed a significant association between NDD-CNVs and EDs (odds ratio = 4.69, P = 0.0023). NDD-CNVs in ED patients included 45,X and deletions at KATNAL2, DIP2A, PTPRT, RBFOX1, CNTN4, MACROD2, and FAM92B. Four of these genes were related to synaptic function. In gene set analysis, we observed a nominally significant enrichment of rare exonic CNVs in synaptic signaling in ED patients (odds ratio = 2.55, P = 0.0254).

CONCLUSION

Our study provides the first preliminary evidence that NDD-CNVs may confer risk for severe EDs. The pathophysiology may involve synaptic dysfunction.

Reactive Oxygen Species and Their Consequences on the Structure and Function of Mammalian Spermatozoa

Significance: Among the 200 or so cell types that comprise mammals, spermatozoa have an ambiguous relationship with the reactive oxygen species (ROS) inherent in the consumption of oxygen that supports aerobic metabolism. Recent Advances: In this review, we shall see that spermatozoa need the action of ROS to reach their structural and functional maturity, but that due to intrinsic unique characteristics, they are, perhaps more than any other cell type, susceptible to oxidative damage. Recent studies have improved our knowledge of how oxidative damage affects sperm structures and functions. The focus of this review will be on how genetic and epigenetic oxidative alterations to spermatozoa can have dramatic unintended consequences in terms of both the support and the suppression of sperm function. Critical Issues: Oxidative stress can have dramatic consequences not only for the spermatozoon itself, but also, and above all, on its primary objective, which is to carry out fertilization and to ensure, in part, that the embryonic development program should lead to a healthy progeny. Future Directions: Sperm oxidative DNA damage largely affects the integrity of the paternal genetic material to such an extent that the oocyte may have difficulties in correcting it. Diagnostic and therapeutic actions should be considered more systematically, especially in men with difficulties to conceive. Research is underway to determine whether the epigenetic information carried by spermatozoa is also subject to changes mediated by pro-oxidative situations. Antioxid. Redox Signal. 37, 481-500.

Rare variants and the oligogenic architecture of autism

Most large-scale genetic studies of autism have focused on the discovery of genes by proving an enrichment of de novo mutations (DNMs) in autism probands or characterizing polygenic risk based on the association of common variants. We present evidence in support of an oligogenic model where two or more ultrarare mutations of more modest effect are preferentially transmitted to children with autism. Such private gene-disruptive mutations are enriched in families where there are multiple affected individuals, emerged two or three generations ago, and map to genes not previously associated with autism. Although no single gene has reached statistical significance, this class of variation should be considered along with genetic and nongenetic factors to better explain the etiology of this complex trait.

Temporal analysis of enhancers during mouse cerebellar development reveals dynamic and novel regulatory functions

We have identified active enhancers in the mouse cerebellum at embryonic and postnatal stages which provides a view of novel enhancers active during cerebellar development. The majority of cerebellar enhancers have dynamic activity between embryonic and postnatal development. Cerebellar enhancers were enriched for neural transcription factor binding sites with temporally specific expression. Putative gene targets displayed spatially restricted expression patterns, indicating cell-type specific expression regulation. Functional analysis of target genes indicated that enhancers regulate processes spanning several developmental epochs such as specification, differentiation and maturation. We use these analyses to discover one novel regulator and one novel marker of cerebellar development: Bhlhe22 and Pax3, respectively. We identified an enrichment of de novo mutations and variants associated with autism spectrum disorder in cerebellar enhancers. Furthermore, by comparing our data with relevant brain development ENCODE histone profiles and cerebellar single-cell datasets we have been able to generalize and expand on the presented analyses, respectively. We have made the results of our analyses available online in the Developing Mouse Cerebellum Enhancer Atlas (https://goldowitzlab.shinyapps.io/developing_mouse_cerebellum_enhancer_atlas/), where our dataset can be efficiently queried, curated and exported by the scientific community to facilitate future research efforts. Our study provides a valuable resource for studying the dynamics of gene expression regulation by enhancers in the developing cerebellum and delivers a rich dataset of novel gene-enhancer associations providing a basis for future in-depth studies in the cerebellum.

Identifying Rare Genetic Variants of Immune Mediators as Risk Factors for Autism Spectrum Disorder

Autism spectrum disorder (ASD) affects more than 1% of children, and there is no viable pharmacotherapeutic agent to treat the core symptoms of ASD. Studies have shown that children with ASD show changes in their levels of immune response molecules. Our previous studies have shown that ASD is more common in children with folate receptor autoantibodies. We also found that children with ASD have abnormal gut immune function, which was characterized by a significant increase in the content of immunoglobulin A and an increase in gut-microbiota-associated epitope diversity. These studies suggest that the immune mechanism plays an important role in the occurrence of ASD. The present study aims to systematically assess gene mutations in immune mediators in patients with ASD. We collected genetic samples from 72 children with ASD (2−12 years old) and 107 healthy controls without ASD (20−78 years old). We used our previously-designed immune gene panel, which can capture cytokine and receptor genes, the coding regions of MHC genes, and genes of innate immunity. Target region sequencing (500×) and bioinformatics analytical methods were used to identify variants in immune response genes associated with patients with ASD. A total of 4 rare variants were found to be associated with ASD, including HLA-B: p.A93G, HLA-DQB1: p.S229N, LILRB2: p.R322H, and LILRB2: c.956-4C>T. These variants were present in 44.44% (32/72) of the ASD patients and were detected in 3.74% (4/107) of the healthy controls. We expect these genetic variants will serve as new targets for the clinical genetic assessment of ASD, and our findings suggest that immune abnormalities in children with ASD may have a genetic basis.

De novo mutations in children born after medical assisted reproduction

STUDY QUESTION

Are there more de novo mutations (DNMs) present in the genomes of children born through medical assisted reproduction (MAR) compared to spontaneously conceived children?

SUMMARY ANSWER

In this pilot study, no statistically significant difference was observed in the number of DNMs observed in the genomes of MAR children versus spontaneously conceived children.

WHAT IS KNOWN ALREADY

DNMs are known to play a major role in sporadic disorders with reduced fitness such as severe developmental disorders, including intellectual disability and epilepsy. Advanced paternal age is known to place offspring at increased disease risk, amongst others by increasing the number of DNMs in their genome. There are very few studies reporting on the effect of MAR on the number of DNMs in the offspring, especially when male infertility is known to be affecting the potential fathers. With delayed parenthood an ongoing epidemiological trend in the 21st century, there are more children born from fathers of advanced age and more children born through MAR every day.

STUDY DESIGN, SIZE, DURATION

This observational pilot study was conducted from January 2015 to March 2019 in the tertiary care centre at Radboud University Medical Center. We included a total of 53 children and their respective parents, forming 49 trios (mother, father and child) and two quartets (mother, father and two siblings). One group of children was born after spontaneous conception (n = 18); a second group of children born after IVF (n = 17) and a third group of children born after ICSI combined with testicular sperm extraction (ICSI-TESE) (n = 18). In this pilot study, we also subdivided each group by paternal age, resulting in a subgroup of children born to younger fathers (<35 years of age at conception) and older fathers (>45 years of age at conception).

PARTICIPANTS/MATERIALS, SETTING, METHODS

Whole-genome sequencing (WGS) was performed on all parent-offspring trios to identify DNMs. For 34 of 53 trios/quartets, WGS was performed twice to independently detect and validate the presence of DNMs. Quality of WGS-based DNM calling was independently assessed by targeted Sanger sequencing.

MAIN RESULTS AND THE ROLE OF CHANCE

No significant differences were observed in the number of DNMs per child for the different methods of conception, independent of parental age at conception (multi-factorial ANOVA, f(2) = 0.17, P-value = 0.85). As expected, a clear paternal age effect was observed after adjusting for method of conception and maternal age at conception (multiple regression model, t = 5.636, P-value = 8.97 × 10-7), with on average 71 DNMs in the genomes of children born to young fathers (<35 years of age) and an average of 94 DNMs in the genomes of children born to older fathers (>45 years of age).

LIMITATIONS, REASONS FOR CAUTION

This is a pilot study and other small-scale studies have recently reported contrasting results. Larger unbiased studies are required to confirm or falsify these results.

WIDER IMPLICATIONS OF THE FINDINGS

This pilot study did not show an effect for the method of conception on the number of DNMs per genome in offspring. Given the role that DNMs play in disease risk, this negative result is good news for IVF and ICSI-TESE born children, if replicated in a larger cohort.

STUDY FUNDING/COMPETING INTEREST(S)

This research was funded by the Netherlands Organisation for Scientific Research (918-15-667) and by an Investigator Award in Science from the Wellcome Trust (209451). The authors have no conflicts of interest to declare.

TRIAL REGISTRATION NUMBER

N/A.

Regionally defined proteomic profiles of human cerebral tissue and organoids reveal conserved molecular modules of neurodevelopment

Cerebral organoids have emerged as robust models for neurodevelopmental and pathological processes, as well as a powerful discovery platform for less-characterized neurobiological programs. Toward this prospect, we leverage mass-spectrometry-based proteomics to molecularly profile precursor and neuronal compartments of both human-derived organoids and mid-gestation fetal brain tissue to define overlapping programs. Our analysis includes recovery of precursor-enriched transcriptional regulatory proteins not found to be differentially expressed in previous transcriptomic datasets. To highlight the discovery potential of this resource, we show that RUVBL2 is preferentially expressed in the SOX2-positive compartment of organoids and that chemical inactivation leads to precursor cell displacement and apoptosis. To explore clinicopathological correlates of this cytoarchitectural disruption, we interrogate clinical datasets and identify rare de novo genetic variants involving RUVBL2 in patients with neurodevelopmental impairments. Together, our findings demonstrate how cell-type-specific profiling of organoids can help nominate previously unappreciated genes in neurodevelopment and disease.

Trio genome sequencing for developmental delay and pediatric heart conditions: A comparative microcost analysis

PURPOSE

Genome sequencing (GS) can aid clinical management of multiple pediatric conditions. Insurers require accurate cost information to inform funding and implementation decisions. The objective was to compare the laboratory workflows and microcosts of trio GS testing in children with developmental delay (DD) and in children with cardiac conditions.

METHODS

Cost items related to each step in trio GS (child and 2 parents) for both populations were identified and measured. Program costs over 5 years were estimated. Probabilistic and deterministic analyses were conducted.

RESULTS

The mean cost per trio GS was CAD$6634.11 (95% CI = 6352.29-6913.40) for DD and CAD$8053.10 (95% CI = 7699.30-8558.10) for cardiac conditions. The 5-year program cost was CAD$28.11 million (95% CI = 26.91-29.29) for DD and CAD$5.63 million (95% CI = 5.38-5.98) for cardiac conditions. Supplies constituted the largest cost component for both populations. The higher cost per sample for the population with cardiac conditions was due to the inclusion of pharmacogenomics, higher bioinformatics labor costs, and a more labor intensive case review.

CONCLUSION

This analysis indicated important variation in trio GS workflow and costs between pediatric populations in a single institution. Enhanced understanding of the clinical utility and costs of GS can inform harmonization and implementation decision-making.

Implications of Dosage Deficiencies in CTCF and Cohesin on Genome Organization, Gene Expression, and Human Neurodevelopment

Properly organizing DNA within the nucleus is critical to ensure normal downstream nuclear functions. CTCF and cohesin act as major architectural proteins, working in concert to generate thousands of high-intensity chromatin loops. Due to their central role in loop formation, a massive research effort has been dedicated to investigating the mechanism by which CTCF and cohesin create these loops. Recent results lead to questioning the direct impact of CTCF loops on gene expression. Additionally, results of controlled depletion experiments in cell lines has indicated that genome architecture may be somewhat resistant to incomplete deficiencies in CTCF or cohesin. However, heterozygous human genetic deficiencies in CTCF and cohesin have illustrated the importance of their dosage in genome architecture, cellular processes, animal behavior, and disease phenotypes. Thus, the importance of considering CTCF or cohesin levels is especially made clear by these heterozygous germline variants that characterize genetic syndromes, which are increasingly recognized in clinical practice. Defined primarily by developmental delay and intellectual disability, the phenotypes of CTCF and cohesin deficiency illustrate the importance of architectural proteins particularly in neurodevelopment. We discuss the distinct roles of CTCF and cohesin in forming chromatin loops, highlight the major role that dosage of each protein plays in the amplitude of observed effects on gene expression, and contrast these results to heterozygous mutation phenotypes in murine models and clinical patients. Insights highlighted by this comparison have implications for future research into these newly emerging genetic syndromes.

Molecular and Clinical Repercussions of GABA Transporter 1 Variants Gone Amiss: Links to Epilepsy and Developmental Spectrum Disorders

The human γ-aminobutyric acid (GABA) transporter 1 (hGAT-1) is the first member of the solute carrier 6 (SLC6) protein superfamily. GAT-1 (SLC6A1) is one of the main GABA transporters in the central nervous system. Its principal physiological role is retrieving GABA from the synapse into neurons and astrocytes, thus swiftly terminating neurotransmission. GABA is a key inhibitory neurotransmitter and shifts in GABAergic signaling can lead to pathological conditions, from anxiety and epileptic seizures to schizophrenia. Point mutations in the SLC6A1 gene frequently give rise to epilepsy, intellectual disability or autism spectrum disorders in the afflicted individuals. The mechanistic routes underlying these are still fairly unclear. Some loss-of-function variants impair the folding and intracellular trafficking of the protein (thus retaining the transporter in the endoplasmic reticulum compartment), whereas others, despite managing to reach their bona fide site of action at the cell surface, nonetheless abolish GABA transport activity (plausibly owing to structural/conformational defects). Whatever the molecular culprit(s), the physiological aftermath transpires into the absence of functional transporters, which in turn perturbs GABAergic actions. Dozens of mutations in the kin SLC6 family members are known to exhort protein misfolding. Such events typically elicit severe ailments in people, e.g., infantile parkinsonism-dystonia or X-linked intellectual disability, in the case of dopamine and creatine transporters, respectively. Flaws in protein folding can be rectified by small molecules known as pharmacological and/or chemical chaperones. The search for such apt remedies calls for a systematic investigation and categorization of the numerous disease-linked variants, by biochemical and pharmacological means in vitro (in cell lines and primary neuronal cultures) and in vivo (in animal models). We here give special emphasis to the utilization of the fruit fly Drosophila melanogaster as a versatile model in GAT-1-related studies. Jointly, these approaches can portray indispensable insights into the molecular factors underlying epilepsy, and ultimately pave the way for contriving efficacious therapeutic options for patients harboring pathogenic mutations in hGAT-1.

Environmental exposures associated with elevated risk for autism spectrum disorder may augment the burden of deleterious de novo mutations among probands

Although the full aetiology of autism spectrum disorder (ASD) is unknown, familial and twin studies demonstrate high heritability of 60-90%, indicating a predominant role of genetics in the development of the disorder. The genetic architecture of ASD consists of a complex array of rare and common variants of all classes of genetic variation usually acting additively to augment individual risk. The relative contribution of heredity in ASD persists despite selective pressures against the classic autistic phenotype; a phenomenon thought to be explained, in part, by the incidence of spontaneous (or de novo) mutations. Notably, environmental exposures attributed as salient risk factors for ASD may play a causal role in the emergence of deleterious de novo variations, with several ASD-associated agents having significant mutagenic potential. To explore this hypothesis, this review article assesses published epidemiological data with evidence derived from assays of mutagenicity, both in vivo and in vitro, to determine the likely role such agents may play in augmenting the genetic liability in ASD. Broadly, these exposures were observed to elicit genomic alterations through one or a combination of: (1) direct interaction with genetic material; (2) impaired DNA repair; or (3) oxidative DNA damage. However, the direct contribution of these factors to the ASD phenotype cannot be determined without further analysis. The development of comprehensive prospective birth cohorts in combination with genome sequencing is essential to forming a causal, mechanistic account of de novo mutations in ASD that links exposure, genotypic alterations, and phenotypic consequences.

Rare Does Not Mean Worthless: How Rare Diseases Have Shaped Neurodevelopment Research in the NGS Era

The advent of next-generation sequencing (NGS) is heavily changing both the diagnosis of human conditions and basic biological research. It is now possible to dig deep inside the genome of hundreds of thousands or even millions of people and find both common and rare genomic variants and to perform detailed phenotypic characterizations of both physiological organs and experimental models. Recent years have seen the introduction of multiple techniques using NGS to profile transcription, DNA and chromatin modifications, protein binding, etc., that are now allowing us to profile cells in bulk or even at a single-cell level. Although rare and ultra-rare diseases only affect a few people, each of these diseases represent scholarly cases from which a great deal can be learned about the pathological and physiological function of genes, pathways, and mechanisms. Therefore, for rare diseases, state-of-the-art investigations using NGS have double valence: their genomic cause (new variants) and the characterize the underlining the mechanisms associated with them (discovery of gene function) can be found. In a non-exhaustive manner, this review will outline the main usage of NGS-based techniques for the diagnosis and characterization of neurodevelopmental disorders (NDDs), under whose umbrella many rare and ultra-rare diseases fall.

A recurrent SHANK3 frameshift variant in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is genetically complex with ~100 copy number variants and genes involved. To try to establish more definitive genotype and phenotype correlations in ASD, we searched genome sequence data, and the literature, for recurrent predicted damaging sequence-level variants affecting single genes. We identified 18 individuals from 16 unrelated families carrying a heterozygous guanine duplication (c.3679dup; p.Ala1227Glyfs*69) occurring within a string of 8 guanines (genomic location [hg38]g.50,721,512dup) affecting SHANK3, a prototypical ASD gene (0.08% of ASD-affected individuals carried the predicted p.Ala1227Glyfs*69 frameshift variant). Most probands carried de novo mutations, but five individuals in three families inherited it through somatic mosaicism. We scrutinized the phenotype of p.Ala1227Glyfs*69 carriers, and while everyone (17/17) formally tested for ASD carried a diagnosis, there was the variable expression of core ASD features both within and between families. Defining such recurrent mutational mechanisms underlying an ASD outcome is important for genetic counseling and early intervention.

Mutational bias in spermatogonia impacts the anatomy of regulatory sites in the human genome

Mutation in the germline is the ultimate source of genetic variation, but little is known about the influence of germline chromatin structure on mutational processes. Using ATAC-seq, we profile the open chromatin landscape of human spermatogonia, the most proliferative cell type of the germline, identifying transcription factor binding sites (TFBSs) and PRDM9 binding sites, a subset of which will initiate meiotic recombination. We observe an increase in rare structural variant (SV) breakpoints at PRDM9-bound sites, implicating meiotic recombination in the generation of structural variation. Many germline TFBSs, such as NRF1, are also associated with increased rates of SV breakpoints, apparently independent of recombination. Singleton short insertions (≥5 bp) are highly enriched at TFBSs, particularly at sites bound by testis active TFs, and their rates correlate with those of structural variant breakpoints. Short insertions often duplicate the TFBS motif, leading to clustering of motif sites near regulatory regions in this male-driven evolutionary process. Increased mutation loads at germline TFBSs disproportionately affect neural enhancers with activity in spermatogonia, potentially altering neurodevelopmental regulatory architecture. Local chromatin structure in spermatogonia is thus pervasive in shaping both evolution and disease.

Prediction of genetic alteration of phospholipase C isozymes in brain disorders: Studies with deep learning

Genetic mutations leading to the development of various diseases, such as cancer, diabetes, and neurodegenerative disorders, can be attributed to multiple mechanisms and exposure to diverse environments. These disorders further increase gene mutation rates and affect the activity of translated proteins, both phenomena associated with cellular responses. Therefore, maintaining the integrity of genetic and epigenetic information is critical for disease suppression and prevention. With the advent of genome sequencing technologies, large-scale genomic data-based machine learning tools, including deep learning, have been used to predict and identify somatic inactivation or negative dominant expression of target genes in various diseases. Although deep learning studies have recently been highlighted for their ability to distinguish between the genetic information of diseases, conventional wisdom is also necessary to explain the correlation between genotype and phenotype. Herein, we summarize the current understanding of phosphoinositide-specific phospholipase C isozymes (PLCs) and an overview of their associations with genetic variation, as well as their emerging roles in several diseases. We also predicted and discussed new findings of cryptic PLC splice variants by deep learning and the clinical implications of the PLC genetic variations predicted using these tools.

Physiological Roles of Mammalian Transmembrane Adenylyl Cyclase Isoforms

Adenylyl cyclases (ACs) catalyze the conversion of ATP to the ubiquitous second messenger cAMP. Mammals possess nine isoforms of transmembrane ACs, dubbed AC1-9, that serve as major effector enzymes of G protein-coupled receptors. The transmembrane ACs display varying expression patterns across tissues, giving potential for them having a wide array of physiologic roles. Cells express multiple AC isoforms, implying that ACs have redundant functions. Furthermore, all transmembrane ACs are activated by Gα_s so it was long assumed that all ACs are activated by Gα_s-coupled GPCRs. AC isoforms partition to different microdomains of the plasma membrane and form prearranged signaling complexes with specific GPCRs that contribute to cAMP signaling compartments. This compartmentation allows for a diversity of cellular and physiological responses by enabling unique signaling events to be triggered by different pools of cAMP. Isoform specific pharmacological activators or inhibitors are lacking for most ACs, making knockdown and overexpression the primary tools for examining the physiological roles of a given isoform. Much progress has been made in understanding the physiological effects mediated through individual transmembrane ACs. GPCR-AC-cAMP signaling pathways play significant roles in regulating functions of every cell and tissue, so understanding each AC isoform's role holds potential for uncovering new approaches for treating a vast array of pathophysiological conditions.

Somatic Mosaicism and Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a genetically heterogenous neurodevelopmental disorder. In the early years of next-generation sequencing, de novo germline variants were shown to contribute to ASD risk. These germline mutations are present in all of the cells of an affected individual and can be detected in any tissue, including clinically accessible DNA sources such as blood or saliva. In recent years, studies have also implicated de novo somatic variants in ASD risk. These somatic mutations arise postzygotically and are present in only a subset of the cells of an affected individual. Depending on the developmental time and progenitor cell in which a somatic mutation occurs, it may be detectable in some tissues and not in others. Somatic mutations detectable at relatively low sequencing coverage in clinically accessible tissues are suggested to contribute to 3-5% of simplex ASD diagnoses, and "brain limited" somatic mutations have been identified in postmortem ASD brain tissue. Somatic mutations likely represent the genetic diagnosis in a proportion of otherwise unexplained individuals with ASD, and brain limited somatic mutations can be used as markers to discover risk genes, cell types, brain regions, and cellular pathways important for ASD pathogenesis and to potentially target for therapeutics.

Genome sequencing identifies rare tandem repeat expansions and copy number variants in Lennox-Gastaut syndrome

Epilepsies are a group of common neurological disorders with a substantial genetic basis. Despite this, the molecular diagnosis of epilepsies remains challenging due to its heterogeneity. Studies utilizing whole-genome sequencing may provide additional insights into genetic causes of epilepsies of unknown aetiology. Whole-genome sequencing was used to evaluate a cohort of adults with unexplained developmental and epileptic encephalopathies (n = 30), for whom prior genetic tests, including whole-exome sequencing in some cases, were negative or inconclusive. Rare single nucleotide variants, insertions/deletions, copy number variants and tandem repeat expansions were analysed. Seven pathogenic or likely pathogenic single nucleotide variants, and two pathogenic deleterious copy number variants were identified in nine patients (32.1% of the cohort). One of the copy number variants, identified in a patient with Lennox–Gastaut syndrome, was too small to be detected by chromosomal microarray techniques. We also identified two tandem repeat expansions with clinical implications in two other patients with Lennox–Gastaut syndrome: a CGG repeat expansion in the 5′untranslated region of DIP2B, and a CTG expansion in ATXN8OS (previously implicated in spinocerebellar ataxia type 8). Three patients had KCNA2 pathogenic variants. One of them died of sudden unexpected death in epilepsy. The other two patients had, in addition to a KCNA2 variant, a second de novo variant impacting potential epilepsy-relevant genes (KCNIP4 and UBR5). Overall, whole-genome sequencing provided a genetic explanation in 32.1% of the total cohort. This is also the first report of coding and non-coding tandem repeat expansions identified in patients with Lennox–Gastaut syndrome. This study demonstrates that using whole-genome sequencing, the examination of multiple types of rare genetic variation, including those found in the non-coding region of the genome, can help resolve unexplained epilepsies.

Prioritizing de novo autism risk variants with calibrated gene- and variant-scoring models

Whole-exome and whole-genome sequencing studies in autism spectrum disorder (ASD) have identified hundreds of thousands of exonic variants. Only a handful of them, primarily loss-of-function variants, have been shown to increase the risk for ASD, while the contributory roles of other variants, including most missense variants, remain unknown. New approaches that combine tissue-specific molecular profiles with patients' genetic data can thus play an important role in elucidating the functional impact of exonic variation and improve understanding of ASD pathogenesis. Here, we integrate spatio-temporal gene co-expression networks from the developing human brain and protein-protein interaction networks to first reach accurate prioritization of ASD risk genes based on their connectivity patterns with previously known high-confidence ASD risk genes. We subsequently integrate these gene scores with variant pathogenicity predictions to further prioritize individual exonic variants based on the positive-unlabeled learning framework with gene- and variant-score calibration. We demonstrate that this approach discriminates among variants between cases and controls at the high end of the prediction range. Finally, we experimentally validate our top-scoring de novo mutation NP_001243143.1:p.Phe309Ser in the sodium/potassium-transporting ATPase ATP1A3 to disrupt protein binding with different partners.

Rates of contributory de novo mutation in high and low-risk autism families

Autism arises in high and low-risk families. De novo mutation contributes to autism incidence in low-risk families as there is a higher incidence in the affected of the simplex families than in their unaffected siblings. But the extent of contribution in low-risk families cannot be determined solely from simplex families as they are a mixture of low and high-risk. The rate of de novo mutation in nearly pure populations of high-risk families, the multiplex families, has not previously been rigorously determined. Moreover, rates of de novo mutation have been underestimated from studies based on low resolution microarrays and whole exome sequencing. Here we report on findings from whole genome sequence (WGS) of both simplex families from the Simons Simplex Collection (SSC) and multiplex families from the Autism Genetic Resource Exchange (AGRE). After removing the multiplex samples with excessive cell-line genetic drift, we find that the contribution of de novo mutation in multiplex is significantly smaller than the contribution in simplex. We use WGS to provide high resolution CNV profiles and to analyze more than coding regions, and revise upward the rate in simplex autism due to an excess of de novo events targeting introns. Based on this study, we now estimate that de novo events contribute to 52-67% of cases of autism arising from low risk families, and 30-39% of cases of all autism.

Discovery of genomic variation across a generation

Over the past 30 years (the timespan of a generation), advances in genomics technologies have revealed tremendous and unexpected variation in the human genome and have provided increasingly accurate answers to long-standing questions of how much genetic variation exists in human populations and to what degree the DNA complement changes between parents and offspring. Tracking the characteristics of these inherited and spontaneous (or de novo) variations has been the basis of the study of human genetic disease. From genome-wide microarray and next-generation sequencing scans, we now know that each human genome contains over 3 million single nucleotide variants when compared with the ~ 3 billion base pairs in the human reference genome, along with roughly an order of magnitude more DNA—approximately 30 megabase pairs (Mb)—being ‘structurally variable’, mostly in the form of indels and copy number changes. Additional large-scale variations include balanced inversions (average of 18 Mb) and complex, difficult-to-resolve alterations. Collectively, ~1% of an individual’s genome will differ from the human reference sequence. When comparing across a generation, fewer than 100 new genetic variants are typically detected in the euchromatic portion of a child’s genome. Driven by increasingly higher-resolution and higher-throughput sequencing technologies, newer and more accurate databases of genetic variation (for instance, more comprehensive structural variation data and phasing of combinations of variants along chromosomes) of worldwide populations will emerge to underpin the next era of discovery in human molecular genetics.

Genetic architecture of autism spectrum disorder: Lessons from large-scale genomic studies

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with a strong genetic component. Recently developed genomic technologies, including microarray and next-generation sequencing (NGS), have enabled researchers to genetic analyses aimed at identifying genetic variations associated with ASD and to elucidate the genetic architecture of the disorder. Large-scale microarray, exome sequencing analyses, and robust statistical methods have resulted in successful gene discovery and identification of high-confidence ASD genes from among de novo and inherited variants. Efforts have been made to understand the genetic architecture of ASD using whole-genome sequencing and genome-wide association studies aimed at identifying noncoding mutations and common variants associated with ASD. In addition, the development of systems biology approaches has resulted in the integration of genetic findings with functional genomic datasets, thereby providing a unique insight into the functional convergence of ASD risk genes and their neurobiology. In this review, we summarize the latest findings of ASD genetic studies involving large cohorts and discuss their implications in ASD neurobiology and in clinical practice.

Epigenomic Convergence of Neural-Immune Risk Factors in Neurodevelopmental Disorder Cortex

Neurodevelopmental disorders (NDDs) affect 7-14% of all children in developed countries and are one of the leading causes of lifelong disability. Epigenetic modifications are poised at the interface between genes and environment and are predicted to reveal insight into NDD etiology. Whole-genome bisulfite sequencing was used to examine DNA cytosine methylation in 49 human cortex samples from 3 different NDDs (autism spectrum disorder, Rett syndrome, and Dup15q syndrome) and matched controls. Integration of methylation changes across NDDs with relevant genomic and genetic datasets revealed differentially methylated regions (DMRs) unique to each type of NDD but with shared regulatory functions in neurons and microglia. NDD DMRs were enriched within promoter regions and for transcription factor binding sites with identified methylation sensitivity. DMRs from all 3 disorders were enriched for ontologies related to nervous system development and genes with disrupted expression in brain from neurodevelopmental or neuropsychiatric disorders. Genes associated with NDD DMRs showed expression patterns indicating an important role for altered microglial function during brain development. These findings demonstrate an NDD epigenomic signature in human cortex that will aid in defining therapeutic targets and early biomarkers at the interface of genetic and environmental NDD risk factors.

HECW2-related disorder in four Japanese patients

The HECT, C2, and WW domain containing E3 ubiquitin protein ligase 2 gene (HECW2) is involved in protein ubiquitination. Several genes associated with protein ubiquitination have been linked to neurodevelopmental disorders. HECW2-related disorder has been established through the identification of de novo variants in HECW2 in patients with neurodevelopmental disorders with hypotonia, seizures, and absent language. Recently, we identified novel HECW2 variants in four Japanese patients with neurodevelopmental disorders. Regarding motor development, two of the patients cannot walk, whereas the other two can walk with an unsteady gait, owing to hypotonia. All HECW2 variants, including those that were previously reported, are missense, and no loss-of-function variants have been identified. Most of the identified variants are located around the HECT domain. These findings suggest that the dominant negative effects of missense variants around the HECT domain may be the mechanism underlying HECW2-related disorder.

Prospective and detailed behavioral phenotyping in DDX3X syndrome

BACKGROUND

DDX3X syndrome is a recently identified genetic disorder that accounts for 1-3% of cases of unexplained developmental delay and/or intellectual disability (ID) in females, and is associated with motor and language delays, and autism spectrum disorder (ASD). To date, the published phenotypic characterization of this syndrome has primarily relied on medical record review; in addition, the behavioral dimensions of the syndrome have not been fully explored.

METHODS

We carried out multi-day, prospective, detailed phenotyping of DDX3X syndrome in 14 females and 1 male, focusing on behavioral, psychological, and neurological measures. Three participants in this cohort were previously reported with limited phenotype information and were re-evaluated for this study. We compared results against population norms and contrasted phenotypes between individuals harboring either (1) protein-truncating variants or (2) missense variants or in-frame deletions.

RESULTS

Eighty percent (80%) of individuals met criteria for ID, 60% for ASD and 53% for attention-deficit/hyperactivity disorder (ADHD). Motor and language delays were common as were sensory processing abnormalities. The cohort included 5 missense, 3 intronic/splice-site, 2 nonsense, 2 frameshift, 2 in-frame deletions, and one initiation codon variant. Genotype-phenotype correlations indicated that, on average, missense variants/in-frame deletions were associated with more severe language, motor, and adaptive deficits in comparison to protein-truncating variants.

LIMITATIONS

Sample size is modest, however, DDX3X syndrome is a rare and underdiagnosed disorder.

CONCLUSION

This study, representing a first, prospective, detailed characterization of DDX3X syndrome, extends our understanding of the neurobehavioral phenotype. Gold-standard diagnostic approaches demonstrated high rates of ID, ASD, and ADHD. In addition, sensory deficits were observed to be a key part of the syndrome. Even with a modest sample, we observe evidence for genotype-phenotype correlations with missense variants/in-frame deletions generally associated with more severe phenotypes.

The germline mutational process in rhesus macaque and its implications for phylogenetic dating

BACKGROUND

Understanding the rate and pattern of germline mutations is of fundamental importance for understanding evolutionary processes.

RESULTS

Here we analyzed 19 parent-offspring trios of rhesus macaques (Macaca mulatta) at high sequencing coverage of ∼76× per individual and estimated a mean rate of 0.77 × 10-8de novo mutations per site per generation (95% CI: 0.69 × 10-8 to 0.85 × 10-8). By phasing 50% of the mutations to parental origins, we found that the mutation rate is positively correlated with the paternal age. The paternal lineage contributed a mean of 81% of the de novo mutations, with a trend of an increasing male contribution for older fathers. Approximately 3.5% of de novo mutations were shared between siblings, with no parental bias, suggesting that they arose from early development (postzygotic) stages. Finally, the divergence times between closely related primates calculated on the basis of the yearly mutation rate of rhesus macaque generally reconcile with divergence estimated with molecular clock methods, except for the Cercopithecoidea/Hominoidea molecular divergence dated at 58 Mya using our new estimate of the yearly mutation rate.

CONCLUSIONS

When compared to the traditional molecular clock methods, new estimated rates from pedigree samples can provide insights into the evolution of well-studied groups such as primates.

MicroRNA-29 is an essential regulator of brain maturation through regulation of CH methylation

Although embryonic brain development and neurodegeneration have received considerable attention, the events that govern postnatal brain maturation are less understood. Here, we identify the miR-29 family to be strikingly induced during the late stages of brain maturation. Brain maturation is associated with a transient, postnatal period of de novo non-CG (CH) DNA methylation mediated by DNMT3A. We examine whether an important function of miR-29 during brain maturation is to restrict the period of CH methylation via its targeting of Dnmt3a. Deletion of miR-29 in the brain, or knockin mutations preventing miR-29 to specifically target Dnmt3a, result in increased DNMT3A expression, higher CH methylation, and repression of genes associated with neuronal activity and neuropsychiatric disorders. These mouse models also develop neurological deficits and premature lethality. Our results identify an essential role for miR-29 in restricting CH methylation in the brain and illustrate the importance of CH methylation regulation for normal brain maturation.

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.

SPEN haploinsufficiency causes a neurodevelopmental disorder overlapping proximal 1p36 deletion syndrome with an episignature of X chromosomes in females

Deletion 1p36 (del1p36) syndrome is the most common human disorder resulting from a terminal autosomal deletion. This condition is molecularly and clinically heterogeneous. Deletions involving two non-overlapping regions, known as the distal (telomeric) and proximal (centromeric) critical regions, are sufficient to cause the majority of the recurrent clinical features, although with different facial features and dysmorphisms. SPEN encodes a transcriptional repressor commonly deleted in proximal del1p36 syndrome and is located centromeric to the proximal 1p36 critical region. Here, we used clinical data from 34 individuals with truncating variants in SPEN to define a neurodevelopmental disorder presenting with features that overlap considerably with those of proximal del1p36 syndrome. The clinical profile of this disease includes developmental delay/intellectual disability, autism spectrum disorder, anxiety, aggressive behavior, attention deficit disorder, hypotonia, brain and spine anomalies, congenital heart defects, high/narrow palate, facial dysmorphisms, and obesity/increased BMI, especially in females. SPEN also emerges as a relevant gene for del1p36 syndrome by co-expression analyses. Finally, we show that haploinsufficiency of SPEN is associated with a distinctive DNA methylation episignature of the X chromosome in affected females, providing further evidence of a specific contribution of the protein to the epigenetic control of this chromosome, and a paradigm of an X chromosome-specific episignature that classifies syndromic traits. We conclude that SPEN is required for multiple developmental processes and SPEN haploinsufficiency is a major contributor to a disorder associated with deletions centromeric to the previously established 1p36 critical regions.