Genomic analysis of 116 autism families strengthens known risk genes and highlights promising candidates

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition with a strong genetic component in which rare variants contribute significantly to risk. We performed whole genome and/or exome sequencing (WGS and WES) and SNP-array analysis to identify both rare sequence and copy number variants (SNVs and CNVs) in 435 individuals from 116 ASD families. We identified 37 rare potentially damaging de novo SNVs (pdSNVs) in the cases (n = 144). Interestingly, two of them (one stop-gain and one missense variant) occurred in the same gene, BRSK2. Moreover, the identification of 8 severe de novo pdSNVs in genes not previously implicated in ASD (AGPAT3, IRX5, MGAT5B, RAB8B, RAP1A, RASAL2, SLC9A1, YME1L1) highlighted promising candidates. Potentially damaging CNVs (pdCNVs) provided support to the involvement of inherited variants in PHF3, NEGR1, TIAM1 and HOMER1 in neurodevelopmental disorders (NDD), although mostly acting as susceptibility factors with incomplete penetrance. Interpretation of identified pdSNVs/pdCNVs according to the ACMG guidelines led to a molecular diagnosis in 19/144 cases, although this figure represents a lower limit and is expected to increase thanks to further clarification of the role of likely pathogenic variants in ASD/NDD candidate genes not yet established. In conclusion, our study highlights promising ASD candidate genes and contributes to characterize the allelic diversity, mode of inheritance and phenotypic impact of de novo and inherited risk variants in ASD/NDD genes.

Whole genome analysis of rare deleterious variants adds further evidence to BRSK2 and other risk genes in Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition with a strong genetic component in which rare variants contribute significantly to risk. We have performed whole genome and/or exome sequencing (WGS and WES) and SNP-array analysis to identify both rare sequence and copy number variants (SNVs and CNVs) in 435 individuals from 116 ASD families. We identified 37 rare potentially damaging de novo SNVs (pdSNVs) in cases (n = 144). Interestingly, two of them (one stop-gain and one missense variant) occurred in the same gene, BRSK2. Moreover, the identification of 9 severe de novo pdSNVs in genes not previously implicated in ASD (RASAL2, RAP1A, IRX5, SLC9A1, AGPAT3, MGAT3, RAB8B, MGAT5B, YME1L1), highlighted novel candidates. Potentially damaging CNVs (pdCNVs) provided support to the involvement of inherited variants in PHF3, NEGR1, TIAM1 and HOMER1 in neurodevelopmental disorders (NDD), although mostly acting as susceptibility factors with incomplete penetrance. Interpretation of identified pdSNVs/pdCNVs according to the ACMG guidelines led to a molecular diagnosis in 19/144 cases, but this figure represents a lower limit and is expected to increase thanks to further clarification of the role of likely pathogenic variants in new ASD/NDD candidates. In conclusion, our study strengthens the role of BRSK2 and other neurodevelopmental genes in ASD risk, highlights novel candidates and contributes to characterize the allelic diversity, mode of inheritance and phenotypic impact of de novo and inherited risk variants in ASD/NDD genes.

Clinical and genetic analysis further delineates the phenotypic spectrum of ALDH1A3-related anophthalmia and microphthalmia

Biallelic pathogenic variants in ALDH1A3 are responsible for approximately 11% of recessively inherited cases of severe developmental eye anomalies. Some individuals can display variable neurodevelopmental features, but the relationship to the ALDH1A3 variants remains unclear. Here, we describe seven unrelated families with biallelic pathogenic ALDH1A3 variants: four compound heterozygous and three homozygous. All affected individuals had bilateral anophthalmia/microphthalmia (A/M), three with additional intellectual or developmental delay, one with autism and seizures and three with facial dysmorphic features. This study confirms that individuals with biallelic pathogenic ALDH1A3 variants consistently manifest A/M, but additionally display neurodevelopmental features with significant intra- and interfamilial variability. Furthermore, we describe the first case with cataract and highlight the importance of screening ALDH1A3 variants in nonconsanguineous families with A/M.

Analysis of Fibroblast Growth Factor 14 (FGF14) structural variants reveals the genetic basis of the early onset nystagmus locus NYS4 and variable ataxia

Nystagmus (involuntary, rhythmical eye movements) can arise due to sensory eye defects, in association with neurological disorders or as an isolated condition. We identified a family with early onset nystagmus and additional neurological features carrying a partial duplication of FGF14, a gene associated with spinocerebellar ataxia type 27 (SCA27) and episodic ataxia. Detailed eye movement analysis revealed oculomotor anomalies strikingly similar to those reported in a previously described four-generation family with early onset nystagmus and linkage to a region on chromosome 13q31.3-q33.1 (NYS4). Since FGF14 lies within NYS4, we revisited the original pedigree using whole genome sequencing, identifying a 161 kb heterozygous deletion disrupting FGF14 and ITGBL1 in the affected individuals, suggesting an FGF14-related condition. Therefore, our study reveals the genetic variant underlying NYS4, expands the spectrum of pathogenic FGF14 variants, and highlights the importance of screening FGF14 in apparently isolated early onset nystagmus.

Genome-wide analyses of individual differences in quantitatively assessed reading- and language-related skills in up to 34,000 people

The use of spoken and written language is a fundamental human capacity. Individual differences in reading- and language-related skills are influenced by genetic variation, with twin-based heritability estimates of 30 to 80% depending on the trait. The genetic architecture is complex, heterogeneous, and multifactorial, but investigations of contributions of single-nucleotide polymorphisms (SNPs) were thus far underpowered. We present a multicohort genome-wide association study (GWAS) of five traits assessed individually using psychometric measures (word reading, nonword reading, spelling, phoneme awareness, and nonword repetition) in samples of 13,633 to 33,959 participants aged 5 to 26 y. We identified genome-wide significant association with word reading (rs11208009, P = 1.098 × 10-8) at a locus that has not been associated with intelligence or educational attainment. All five reading-/language-related traits showed robust SNP heritability, accounting for 13 to 26% of trait variability. Genomic structural equation modeling revealed a shared genetic factor explaining most of the variation in word/nonword reading, spelling, and phoneme awareness, which only partially overlapped with genetic variation contributing to nonword repetition, intelligence, and educational attainment. A multivariate GWAS of word/nonword reading, spelling, and phoneme awareness maximized power for follow-up investigation. Genetic correlation analysis with neuroimaging traits identified an association with the surface area of the banks of the left superior temporal sulcus, a brain region linked to the processing of spoken and written language. Heritability was enriched for genomic elements regulating gene expression in the fetal brain and in chromosomal regions that are depleted of Neanderthal variants. Together, these results provide avenues for deciphering the biological underpinnings of uniquely human traits.

Contribution of CACNA1H Variants in Autism Spectrum Disorder Susceptibility

Autism Spectrum Disorder (ASD) is a highly heterogeneous neuropsychiatric disorder with a strong genetic component. The genetic architecture is complex, consisting of a combination of common low-risk and more penetrant rare variants. Voltage-gated calcium channels (VGCCs or Ca_v) genes have been implicated as high-confidence susceptibility genes for ASD, in accordance with the relevant role of calcium signaling in neuronal function. In order to further investigate the involvement of VGCCs rare variants in ASD susceptibility, we performed whole genome sequencing analysis in a cohort of 105 families, composed of 124 ASD individuals, 210 parents and 58 unaffected siblings. We identified 53 rare inherited damaging variants in Ca_v genes, including genes coding for the principal subunit and genes coding for the auxiliary subunits, in 40 ASD families. Interestingly, biallelic rare damaging missense variants were detected in the CACNA1H gene, coding for the T-type Ca_v3.2 channel, in ASD probands from two different families. Thus, to clarify the role of these CACNA1H variants on calcium channel activity we performed electrophysiological analysis using whole-cell patch clamp technology. Three out of four tested variants were shown to mildly affect Ca_v3.2 channel current density and activation properties, possibly leading to a dysregulation of intracellular Ca²⁺ ions homeostasis, thus altering calcium-dependent neuronal processes and contributing to ASD etiology in these families. Our results provide further support for the role of CACNA1H in neurodevelopmental disorders and suggest that rare CACNA1H variants may be involved in ASD development, providing a high-risk genetic background.

Biallelic variants in the small optic lobe calpain CAPN15 are associated with congenital eye anomalies, deafness and other neurodevelopmental deficits

Microphthalmia, coloboma and cataract are part of a spectrum of developmental eye disorders in humans affecting ~12 per 100 000 live births. Currently, variants in over 100 genes are known to underlie these conditions. However, at least 40% of affected individuals remain without a clinical genetic diagnosis, suggesting variants in additional genes may be responsible. Calpain 15 (CAPN15) is an intracellular cysteine protease belonging to the non-classical small optic lobe (SOL) family of calpains, an important class of developmental proteins, as yet uncharacterized in vertebrates. We identified five individuals with microphthalmia and/or coloboma from four independent families carrying homozygous or compound heterozygous predicted damaging variants in CAPN15. Several individuals had additional phenotypes including growth deficits, developmental delay and hearing loss. We generated Capn15 knockout mice that exhibited similar severe developmental eye defects, including anophthalmia, microphthalmia and cataract, and diminished growth. We demonstrate widespread Capn15 expression throughout the brain and central nervous system, strongest during early development, and decreasing postnatally. Together, these findings demonstrate a critical role of CAPN15 in vertebrate developmental eye disorders, and may signify a new developmental pathway.

An increased burden of rare exonic variants in NRXN1 microdeletion carriers is likely to enhance the penetrance for autism spectrum disorder

Autism spectrum disorder (ASD) is characterized by a complex polygenic background, but with the unique feature of a subset of cases (~15%‐30%) presenting a rare large‐effect variant. However, clinical interpretation in these cases is often complicated by incomplete penetrance, variable expressivity and different neurodevelopmental trajectories. NRXN1 intragenic deletions represent the prototype of such ASD‐associated susceptibility variants. From chromosomal microarrays analysis of 104 ASD individuals, we identified an inherited NRXN1 deletion in a trio family. We carried out whole‐exome sequencing and deep sequencing of mitochondrial DNA (mtDNA) in this family, to evaluate the burden of rare variants which may contribute to the phenotypic outcome in NRXN1 deletion carriers. We identified an increased burden of exonic rare variants in the ASD child compared to the unaffected NRXN1 deletion‐transmitting mother, which remains significant if we restrict the analysis to potentially deleterious rare variants only (P = 6.07 × 10⁻⁵). We also detected significant interaction enrichment among genes with damaging variants in the proband, suggesting that additional rare variants in interacting genes collectively contribute to cross the liability threshold for ASD. Finally, the proband's mtDNA presented five low‐level heteroplasmic mtDNA variants that were absent in the mother, and two maternally inherited variants with increased heteroplasmic load. This study underlines the importance of a comprehensive assessment of the genomic background in carriers of large‐effect variants, as penetrance modulation by additional interacting rare variants to might represent a widespread mechanism in neurodevelopmental disorders.

De Novo Missense Variants in FBXW11 Cause Diverse Developmental Phenotypes Including Brain, Eye, and Digit Anomalies

The identification of genetic variants implicated in human developmental disorders has been revolutionized by second-generation sequencing combined with international pooling of cases. Here, we describe seven individuals who have diverse yet overlapping developmental anomalies, and who all have de novo missense FBXW11 variants identified by whole exome or whole genome sequencing and not reported in the gnomAD database. Their phenotypes include striking neurodevelopmental, digital, jaw, and eye anomalies, and in one individual, features resembling Noonan syndrome, a condition caused by dysregulated RAS signaling. FBXW11 encodes an F-box protein, part of the Skp1-cullin-F-box (SCF) ubiquitin ligase complex, involved in ubiquitination and proteasomal degradation and thus fundamental to many protein regulatory processes. FBXW11 targets include β-catenin and GLI transcription factors, key mediators of Wnt and Hh signaling, respectively, critical to digital, neurological, and eye development. Structural analyses indicate affected residues cluster at the surface of the loops of the substrate-binding domain of FBXW11, and the variants are predicted to destabilize the protein and/or its interactions. In situ hybridization studies on human and zebrafish embryonic tissues demonstrate FBXW11 is expressed in the developing eye, brain, mandibular processes, and limb buds or pectoral fins. Knockdown of the zebrafish FBXW11 orthologs fbxw11a and fbxw11b resulted in embryos with smaller, misshapen, and underdeveloped eyes and abnormal jaw and pectoral fin development. Our findings support the role of FBXW11 in multiple developmental processes, including those involving the brain, eye, digits, and jaw.

New GJA8 variants and phenotypes highlight its critical role in a broad spectrum of eye anomalies

GJA8 encodes connexin 50 (Cx50), a transmembrane protein involved in the formation of lens gap junctions. GJA8 mutations have been linked to early onset cataracts in humans and animal models. In mice, missense mutations and homozygous Gja8 deletions lead to smaller lenses and microphthalmia in addition to cataract, suggesting that Gja8 may play a role in both lens development and ocular growth. Following screening of GJA8 in a cohort of 426 individuals with severe congenital eye anomalies, primarily anophthalmia, microphthalmia and coloboma, we identified four known [p.(Thr39Arg), p.(Trp45Leu), p.(Asp51Asn), and p.(Gly94Arg)] and two novel [p.(Phe70Leu) and p.(Val97Gly)] likely pathogenic variants in seven families. Five of these co-segregated with cataracts and microphthalmia, whereas the variant p.(Gly94Arg) was identified in an individual with congenital aphakia, sclerocornea, microphthalmia and coloboma. Four missense variants of unknown or unlikely clinical significance were also identified. Furthermore, the screening of GJA8 structural variants in a subgroup of 188 individuals identified heterozygous 1q21 microdeletions in five families with coloboma and other ocular and/or extraocular findings. However, the exact genotype-phenotype correlation of these structural variants remains to be established. Our data expand the spectrum of GJA8 variants and associated phenotypes, confirming the importance of this gene in early eye development.

New variant and expression studies provide further insight into the genotype-phenotype correlation in YAP1-related developmental eye disorders

YAP1, which encodes the Yes-associated protein 1, is part of the Hippo pathway involved in development, growth, repair and homeostasis. Nonsense YAP1 mutations have been shown to co-segregate with autosomal dominantly inherited coloboma. Therefore, we screened YAP1 for variants in a cohort of 258 undiagnosed UK patients with developmental eye disorders, including anophthalmia, microphthalmia and coloboma. We identified a novel 1 bp deletion in YAP1 in a boy with bilateral microphthalmia and bilateral chorioretinal coloboma. This variant is located in the coding region of all nine YAP1 spliceforms, and results in a frameshift and subsequent premature termination codon in each. The variant is predicted to result in the loss of part of the transactivation domain of YAP1, and sequencing of cDNA from the patient shows it does not result in nonsense mediated decay. To investigate the role of YAP1 in human eye development, we performed in situ hybridisation utilising human embryonic tissue, and observed expression in the developing eye, neural tube, brain and kidney. These findings help confirm the role of YAP1 and the Hippo developmental pathway in human eye development and its associated anomalies and demonstrate its expression during development in affected organ systems.

Multidisciplinary investigation links backward-speech trait and working memory through genetic mutation

Case studies of unusual traits can provide unique snapshots of the effects of modified systems. In this study, we report on an individual from a Serbian family with the ability to rapidly, accurately and voluntarily speak backwards. We consider psychological, neural and genetic correlates of this trait to identify specific relevant neural mechanisms and new molecular pathways for working memory and speech-related tasks. EEG data suggest that the effect of word reversal precedes semantic integration of visually presented backward-words, and that event-related potentials above the frontal lobe are affected by both word reversal and the maintenance of backward-words in working memory. fMRI revealed that the left fusiform gyrus may facilitate the production of backward-speech. Exome sequencing identified three novel coding variants of potential significance in the RIC3, RIPK1 and ZBED5 genes. Taken together, our data suggest that, in this individual, the ability to speak backwards is afforded by an extraordinary working memory capacity. We hypothesise that this is served by cholinergic projections from the basal forebrain to the frontal cortex and supported by visual semantic loops within the left fusiform gyrus and that these neural processes may be mediated by a genetic mutation in RIC3; a chaperone for nicotinic acetylcholine receptors.

Genome-wide analysis identifies a role for common copy number variants in specific language impairment

An exploratory genome-wide copy number variant (CNV) study was performed in 127 independent cases with specific language impairment (SLI), their first-degree relatives (385 individuals) and 269 population controls. Language-impaired cases showed an increased CNV burden in terms of the average number of events (11.28 vs 10.01, empirical P=0.003), the total length of CNVs (717 vs 513 Kb, empirical P=0.0001), the average CNV size (63.75 vs 51.6 Kb, empirical P=0.0005) and the number of genes spanned (14.29 vs 10.34, empirical P=0.0007) when compared with population controls, suggesting that CNVs may contribute to SLI risk. A similar trend was observed in first-degree relatives regardless of affection status. The increased burden found in our study was not driven by large or de novo events, which have been described as causative in other neurodevelopmental disorders. Nevertheless, de novo CNVs might be important on a case-by-case basis, as indicated by identification of events affecting relevant genes, such as ACTR2 and CSNK1A1, and small events within known micro-deletion/-duplication syndrome regions, such as chr8p23.1. Pathway analysis of the genes present within the CNVs of the independent cases identified significant overrepresentation of acetylcholine binding, cyclic-nucleotide phosphodiesterase activity and MHC proteins as compared with controls. Taken together, our data suggest that the majority of the risk conferred by CNVs in SLI is via common, inherited events within a ‘common disorder–common variant' model. Therefore the risk conferred by CNVs will depend upon the combination of events inherited (both CNVs and SNPs), the genetic background of the individual and the environmental factors.

A CTNNA3 compound heterozygous deletion implicates a role for αT-catenin in susceptibility to autism spectrum disorder

Background

Autism spectrum disorder (ASD) is a highly heritable, neurodevelopmental condition showing extreme genetic heterogeneity. While it is well established that rare genetic variation, both de novo and inherited, plays an important role in ASD risk, recent studies also support a rare recessive contribution.

Methods

We identified a compound heterozygous deletion intersecting the CTNNA3 gene, encoding αT-catenin, in a proband with ASD and moderate intellectual disability. The deletion breakpoints were mapped at base-pair resolution, and segregation analysis was performed. We compared the frequency of CTNNA3 exonic deletions in 2,147 ASD cases from the Autism Genome Project (AGP) study versus the frequency in 6,639 controls. Western blot analysis was performed to get a quantitative characterisation of Ctnna3 expression during early brain development in mouse.

Results

The CTNNA3 compound heterozygous deletion includes a coding exon, leading to a putative frameshift and premature stop codon. Segregation analysis in the family showed that the unaffected sister is heterozygote for the deletion, having only inherited the paternal deletion. While the frequency of CTNNA3 exonic deletions is not significantly different between ASD cases and controls, no homozygous or compound heterozygous exonic deletions were found in a sample of over 6,000 controls. Expression analysis of Ctnna3 in the mouse cortex and hippocampus (P0-P90) provided support for its role in the early stage of brain development.

Conclusion

The finding of a rare compound heterozygous CTNNA3 exonic deletion segregating with ASD, the absence of CTNNA3 homozygous exonic deletions in controls and the high expression of Ctnna3 in both brain areas analysed implicate CTNNA3 in ASD susceptibility.

A deletion involving CD38 and BST1 results in a fusion transcript in a patient with autism and asthma

CD38 encodes a ligand in the oxytocin signaling pathway. Some single nucleotide polymorphisms in this gene have been associated with low serum oxytocin levels in autism spectrum disorder (ASD) patients. Oxytocin disruption has been hypothesized to account for features of ASD, including impaired communication and social behavior, based on animal studies. Recent human studies have shown administration of oxytocin improving emotion recognition, promoting social behavior, and improving auditory processing of social stimuli in ASD patients. In addition to its role in oxytocin signaling, CD38 is involved in the regulation of calcium concentration in airway smooth muscle with impairment of CD38 being implicated in airway diseases like asthma. While a number of studies have implicated rare chromosomal deletions and duplications in helping determine genetic risk for autism, there are to our knowledge no reports describing rearrangements involving CD38 or deletions in patients with ASD. Here, we present two sisters diagnosed with autism and with features of regression-previously acquired speech lost in the second year of life. The younger sister, who also had asthma, inherited a maternal deletion of 4p15.32 that results in a BST1-CD38 fusion transcript. Their mother's deletion was mosaic and she was not affected. Although further work is required to assess functional consequences of the fusion transcript, we hypothesize that the proband's deletion may have served as a risk factor for autism that, when combined with other susceptibility variants, resulted in a more severe presentation than her sister.

Homozygous microdeletion of exon 5 in ZNF277 in a girl with specific language impairment

Specific language impairment (SLI), an unexpected failure to develop appropriate language skills despite adequate non-verbal intelligence, is a heterogeneous multifactorial disorder with a complex genetic basis. We identified a homozygous microdeletion of 21,379 bp in the ZNF277 gene (NM_021994.2), encompassing exon 5, in an individual with severe receptive and expressive language impairment. The microdeletion was not found in the proband's affected sister or her brother who had mild language impairment. However, it was inherited from both parents, each of whom carries a heterozygous microdeletion and has a history of language problems. The microdeletion falls within the AUTS1 locus, a region linked to autistic spectrum disorders (ASDs). Moreover, ZNF277 is adjacent to the DOCK4 and IMMP2L genes, which have been implicated in ASD. We screened for the presence of ZNF277 microdeletions in cohorts of children with SLI or ASD and panels of control subjects. ZNF277 microdeletions were at an increased allelic frequency in SLI probands (1.1%) compared with both ASD family members (0.3%) and independent controls (0.4%). We performed quantitative RT-PCR analyses of the expression of IMMP2L, DOCK4 and ZNF277 in individuals carrying either an IMMP2L_DOCK4 microdeletion or a ZNF277 microdeletion. Although ZNF277 microdeletions reduce the expression of ZNF277, they do not alter the levels of DOCK4 or IMMP2L transcripts. Conversely, IMMP2L_DOCK4 microdeletions do not affect the expression levels of ZNF277. We postulate that ZNF277 microdeletions may contribute to the risk of language impairments in a manner that is independent of the autism risk loci previously described in this region.