Rare variants and the oligogenic architecture of autism

Rare variant association studies: Significance, methods, and applications in chronic pain studies

Rare genetic variants, characterized by their low frequency in a population, have emerged as essential components in the study of complex disease genetics. The biology of rare variants underscores their significance, as they can exert profound effects on phenotypic variation and disease susceptibility. Recent advancements in sequencing technologies have yielded the availability of large-scale sequencing data such as the UK Biobank whole-exome sequencing (WES) cohort empowered researchers to conduct rare variant association studies (RVASs). This review paper discusses the significance of rare variants, available methodologies, and applications. We provide an overview of RVASs, emphasizing their relevance in unraveling the genetic architecture of complex diseases with special focus on chronic pain and Arthritis. Additionally, we discuss the strengths and limitations of various rare variant association testing methods, outlining a typical pipeline for conducting rare variant association. This pipeline encompasses crucial steps such as quality control of WES data, rare variant annotation, and association testing. It serves as a comprehensive guide for researchers in the field of chronic pain diseases interested in rare variant association studies in large-scale sequencing datasets like the UK Biobank WES cohort. Lastly, we discuss how the identified variants can be further investigated through detailed experimental studies in animal models to elucidate their functional impact and underlying mechanisms.

Polygenicity in a box: Copy number variants, neural circuit development, and neurodevelopmental disorders

Clinically defined neurodevelopmental disorders (cd-NDDs), including Autistic Spectrum Disorder (ASD) and Schizophrenia (Scz), are primarily polygenic: Multiple risk genes distributed across the genome, in potentially infinite combinations, account for variable pathology. Polygenicity raises a fundamental question: Can "core" cd-NDD pathogenic mechanisms be identified given this genomic complexity? With the right models and analytic targets, a distinct class of polygenic mutations-Copy Number Variants (CNVs): contiguous gene deletions or duplications associated with cd-NDD risk-provide a singular opportunity to define cd-NDD pathology. CNVs orthologous to those that confer cd-NDD risk have been engineered in animals as well as human stem cells. Using these tools, one can determine how altered function of multiple genes cause serial stumbles over cell biological steps typically taken to build optimal "polygenic" neural circuits. Thus, cd-NDD pathology may be a consequence of polygenic deviations-stumbles-that exceed limits of adaptive variation for key developmental steps.

Phenotypic spectrum of dual diagnoses in developmental disorders

As more patients receive genome-wide sequencing, the number of individuals diagnosed with multiple monogenic conditions is increasing. We sought to investigate the relative phenotypic contribution of dual diagnoses using both manual curation and computational approaches. First, we computed 1,003,236 semantic similarity scores for all possible pairs of 1,417 genes in the Developmental Disorder Gene2Phenotype (DDG2P) database using Human Phenotype Ontology terms. Next, for 62 probands with two molecular diagnoses in the Deciphering Developmental Disorders study, we computed semantic similarity scores between the probands' phenotypes and DDG2P phenotypes associated with the two disorders and compared the results with manual attribution of proband phenotypes to none, one, or both of the genes. We found a spectrum of phenotypic similarity for dual diagnoses, both across all DDG2P genes and within dual diagnosed probands, from phenotypically distinct through blended to indistinguishable conditions. Pairwise semantic similarity scores between two DDG2P genes were a good predictor of the extent of phenotypic blending observed in probands. Dual diagnoses involving genes linked with synergistic phenotypes can result in more extreme presentations while those involving antagonistic phenotypes have spuriously high pairwise semantic similarity scores despite a potentially milder atypical presentation. We suggest that the phenotypic contribution of two molecular diagnoses may contain discrete, synergistic, or antagonistic elements. Conceptual recognition of this phenotypic spectrum is important for making a final clinico-molecular diagnosis and providing accurate genetic counseling.

Addressing Technical Pitfalls in Pursuit of Molecular Factors That Mediate Immunoglobulin Gene Regulation

The expressed Ab repertoire is a critical determinant of immune-related phenotypes. Ab-encoding transcripts are distinct from other expressed genes because they are transcribed from somatically rearranged gene segments. Human Abs are composed of two identical H and L chain polypeptides derived from genes in IGH locus and one of two L chain loci. The combinatorial diversity that results from Ab gene rearrangement and the pairing of different H and L chains contributes to the immense diversity of the baseline Ab repertoire. During rearrangement, Ab gene selection is mediated by factors that influence chromatin architecture, promoter/enhancer activity, and V(D)J recombination. Interindividual variation in the composition of the Ab repertoire associates with germline variation in IGH, implicating polymorphism in Ab gene regulation. Determining how IGH variants directly mediate gene regulation will require integration of these variants with other functional genomic datasets. In this study, we argue that standard approaches using short reads have limited utility for characterizing regulatory regions in IGH at haplotype resolution. Using simulated and chromatin immunoprecipitation sequencing reads, we define features of IGH that limit use of short reads and a single reference genome, namely 1) the highly duplicated nature of the DNA sequence in IGH and 2) structural polymorphisms that are frequent in the population. We demonstrate that personalized diploid references enhance performance of short-read data for characterizing mappable portions of the locus, while also showing that long-read profiling tools will ultimately be needed to fully resolve functional impacts of IGH germline variation on expressed Ab repertoires.

Fetal Gene Regulatory Gene Deletions are Associated with Poor Cognition and Altered Cortical Morphology in Schizophrenia and Community-Based Samples

Schizophrenia spectrum disorders (SSDs) are characterized by substantial clinical and genetic heterogeneity. Multiple recurrent copy number variants (CNVs) increase risk for SSDs; however, how known risk CNVs and broader genome-wide CNVs influence clinical variability is unclear. The current study examined associations between borderline intellectual functioning or childhood-onset psychosis, known risk CNVs, and burden of deletions affecting genes in 18 previously validated neurodevelopmental gene-sets in 618 SSD individuals. CNV associations were assessed for replication in 235 SSD relatives and 583 controls, and 9,930 youth from the Adolescent Brain Cognitive Development (ABCD) Study. Known SSD- and neurodevelopmental disorder (NDD)-risk CNVs were associated with borderline intellectual functioning in SSD cases (odds ratios (OR) = 7.09 and 4.57, respectively); NDD-risk deletions were nominally associated with childhood-onset psychosis (OR = 4.34). Furthermore, deletion of genes involved in regulating gene expression during fetal brain development was associated with borderline intellectual functioning across SSD cases and non-cases (OR = 2.58), with partial replication in the ABCD cohort. Exploratory analyses of cortical morphology showed associations between fetal gene regulatory gene deletions and altered gray matter volume and cortical thickness across cohorts. Results highlight contributions of known risk CNVs to phenotypic variability in SSD and the utility of a neurodevelopmental framework for identifying mechanisms that influence phenotypic variability in SSDs, as well as the broader population, with implications for personalized medicine approaches to care.

Evidence for the additivity of rare and common variant burden throughout the spectrum of intellectual disability

Intellectual disability (ID) is a common disorder, yet there is a wide spectrum of impairment from mild to profoundly affected individuals. Mild ID is seen as the low extreme of the general distribution of intelligence, while severe ID is often seen as a monogenic disorder caused by rare, pathogenic, highly penetrant variants. To investigate the genetic factors influencing mild and severe ID, we evaluated rare and common variation in the Northern Finland Intellectual Disability cohort (n = 1096 ID patients), a cohort with a high percentage of mild ID (n = 550) and from a population bottleneck enriched in rare, damaging variation. Despite this enrichment, we found only a small percentage of ID was due to recessive Finnish-enriched variants (0.5%). A larger proportion was linked to dominant variation, with a significant burden of rare, damaging variation in both mild and severe ID. This rare variant burden was enriched in more severe ID (p = 2.4e-4), patients without a relative with ID (p = 4.76e-4), and in those with features associated with monogenic disorders. We also found a significant burden of common variants associated with decreased cognitive function, with no difference between mild and more severe ID. When we included common and rare variants in a joint model, the rare and common variants had additive effects in both mild and severe ID. A multimodel inference approach also found that common and rare variants together best explained ID status (ΔAIC = 16.8, ΔBIC = 10.2). Overall, we report evidence for the additivity of rare and common variant burden throughout the spectrum of intellectual disability.

In-Depth Genomic Analysis: The New Challenge in Congenital Heart Disease

The use of next-generation sequencing has provided new insights into the causes and mechanisms of congenital heart disease (CHD). Examinations of the whole exome sequence have detected detrimental gene variations modifying single or contiguous nucleotides, which are characterised as pathogenic based on statistical assessments of families and correlations with congenital heart disease, elevated expression during heart development, and reductions in harmful protein-coding mutations in the general population. Patients with CHD and extracardiac abnormalities are enriched for gene classes meeting these criteria, supporting a common set of pathways in the organogenesis of CHDs. Single-cell transcriptomics data have revealed the expression of genes associated with CHD in specific cell types, and emerging evidence suggests that genetic mutations disrupt multicellular genes essential for cardiogenesis. Metrics and units are being tracked in whole-genome sequencing studies.

An evolutionary perspective on complex neuropsychiatric disease

The forces of evolution-mutation, selection, migration, and genetic drift-shape the genetic architecture of human traits, including the genetic architecture of complex neuropsychiatric illnesses. Studying these illnesses in populations that are diverse in genetic ancestry, historical demography, and cultural history can reveal how evolutionary forces have guided adaptation over time and place. A fundamental truth of shared human biology is that an allele responsible for a disease in anyone, anywhere, reveals a gene critical to the normal biology underlying that condition in everyone, everywhere. Understanding the genetic causes of neuropsychiatric disease in the widest possible range of human populations thus yields the greatest possible range of insight into genes critical to human brain development. In this perspective, we explore some of the relationships between genes, adaptation, and history that can be illuminated by an evolutionary perspective on studies of complex neuropsychiatric disease in diverse populations.

Comparison of three bioinformatics tools in the detection of ASD candidate variants from whole exome sequencing data

Autism spectrum disorder (ASD) is a heterogenous multifactorial neurodevelopmental condition with a significant genetic susceptibility component. Thus, identifying genetic variations associated with ASD is a complex task. Whole-exome sequencing (WES) is an effective approach for detecting extremely rare protein-coding single-nucleotide variants (SNVs) and short insertions/deletions (INDELs). However, interpreting these variants' functional and clinical consequences requires integrating multifaceted genomic information. We compared the concordance and effectiveness of three bioinformatics tools in detecting ASD candidate variants (SNVs and short INDELs) from WES data of 220 ASD family trios registered in the National Autism Database of Israel. We studied only rare (< 1% population frequency) proband-specific variants. According to the American College of Medical Genetics (ACMG) guidelines, the pathogenicity of variants was evaluated by the InterVar and TAPES tools. In addition, likely gene-disrupting (LGD) variants were detected based on an in-house bioinformatics tool, Psi-Variant, that integrates results from seven in-silico prediction tools. Overall, 372 variants in 311 genes distributed in 168 probands were detected by these tools. The overlap between the tools was 64.1, 22.9, and 23.1% for InterVar-TAPES, InterVar-Psi-Variant, and TAPES-Psi-Variant, respectively. The intersection between InterVar and Psi-Variant (I ∩ P) was the most effective approach in detecting variants in known ASD genes (PPV = 0.274; OR = 7.09, 95% CI = 3.92-12.22), while the union of InterVar and Psi Variant (I U P) achieved the highest diagnostic yield (20.5%).Our results suggest that integrating different variant interpretation approaches in detecting ASD candidate variants from WES data is superior to each approach alone. The inclusion of additional criteria could further improve the detection of ASD candidate variants.

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.

Whole Genome Sequencing Revealed Inherited Rare Oligogenic Variants Contributing to Schizophrenia and Major Depressive Disorder in Two Families

Schizophrenia and affective disorder are two major complex mental disorders with high heritability. Evidence shows that rare variants with significant clinical impacts contribute to the genetic liability of these two disorders. Also, rare variants associated with schizophrenia and affective disorders are highly personalized; each patient may carry different variants. We used whole genome sequencing analysis to study the genetic basis of two families with schizophrenia and major depressive disorder. We did not detect de novo, autosomal dominant, or recessive pathogenic or likely pathogenic variants associated with psychiatric disorders in these two families. Nevertheless, we identified multiple rare inherited variants with unknown significance in the probands. In family 1, with singleton schizophrenia, we detected four rare variants in genes implicated in schizophrenia, including p.Arg1627Trp of LAMA2, p.Pro1338Ser of CSMD1, p.Arg691Gly of TLR4, and Arg182X of AGTR2. The p.Arg691Gly of TLR4 was inherited from the father, while the other three were inherited from the mother. In family 2, with two affected sisters diagnosed with major depressive disorder, we detected three rare variants shared by the two sisters in three genes implicated in affective disorders, including p.Ala4551Gly of FAT1, p.Val231Leu of HOMER3, and p.Ile185Met of GPM6B. These three rare variants were assumed to be inherited from their parents. Prompted by these findings, we suggest that these rare inherited variants may interact with each other and lead to psychiatric conditions in these two families. Our observations support the conclusion that inherited rare variants may contribute to the heritability of psychiatric disorders.

The association between polygenic risk scores for mental disorders and social cognition: A scoping review

People with mental disorders, such as psychosis or autism spectrum disorder (ASD), often present impairments in social cognition (SC), which may cause significant difficulties in real-world functioning. SC deficits are seen also in unaffected relatives, indicating a genetic substratum. The present review evaluated the evidence on the association between SC and the polygenic risk score (PRS), a single metric of the molecular genetic risk to develop a specific disorder. In July 2022, we conducted systematic searches in Scopus and PubMed following the PRISMA-ScR guidelines. We selected original articles written in English reporting results on the association between PRSs for any mental disorder and domains of SC either in people with mental disorders or controls. The search yielded 244 papers, of which 13 were selected for inclusion. Studies tested mainly PRSs for schizophrenia, ASD, and attention-deficit hyperactivity disorder. Emotion recognition was the most investigated domain of SC. Overall, evidence revealed that currently available PRSs for mental disorders do not explain variation in SC performances. To enhance the understanding of mechanisms underlying SC in mental disorders, future research should focus on the development of transdiagnostic PRSs, study their interaction with environmental risk factors, and standardize outcome measurement.

A systematic review of non-coding RNA genes with differential expression profiles associated with autism spectrum disorders

AIMS

To identify differential expression of shorter non-coding RNA (ncRNA) genes associated with autism spectrum disorders (ASD).

BACKGROUND

ncRNA are functional molecules that derive from non-translated DNA sequence. The HUGO Gene Nomenclature Committee (HGNC) have approved ncRNA gene classes with alignment to the reference human genome. One subset is microRNA (miRNA), which are highly conserved, short RNA molecules that regulate gene expression by direct post-transcriptional repression of messenger RNA. Several miRNA genes are implicated in the development and regulation of the nervous system. Expression of miRNA genes in ASD cohorts have been examined by multiple research groups. Other shorter classes of ncRNA have been examined less. A comprehensive systematic review examining expression of shorter ncRNA gene classes in ASD is timely to inform the direction of research.

METHODS

We extracted data from studies examining ncRNA gene expression in ASD compared with non-ASD controls. We included studies on miRNA, piwi-interacting RNA (piRNA), small NF90 (ILF3) associated RNA (snaR), small nuclear RNA (snRNA), small nucleolar RNA (snoRNA), transfer RNA (tRNA), vault RNA (vtRNA) and Y RNA. The following electronic databases were searched: Cochrane Library, EMBASE, PubMed, Web of Science, PsycINFO, ERIC, AMED and CINAHL for papers published from January 2000 to May 2022. Studies were screened by two independent investigators with a third resolving discrepancies. Data was extracted from eligible papers.

RESULTS

Forty-eight eligible studies were included in our systematic review with the majority examining miRNA gene expression alone. Sixty-four miRNA genes had differential expression in ASD compared to controls as reported in two or more studies, but often in opposing directions. Four miRNA genes had differential expression in the same direction in the same tissue type in at least 3 separate studies. Increased expression was reported in miR-106b-5p, miR-155-5p and miR-146a-5p in blood, post-mortem brain, and across several tissue types, respectively. Decreased expression was reported in miR-328-3p in bloods samples. Seven studies examined differential expression from other classes of ncRNA, including piRNA, snRNA, snoRNA and Y RNA. No individual ncRNA genes were reported in more than one study. Six studies reported differentially expressed snoRNA genes in ASD. A meta-analysis was not possible because of inconsistent methodologies, disparate tissue types examined, and varying forms of data presented.

CONCLUSION

There is limited but promising evidence associating the expression of certain miRNA genes and ASD, although the studies are of variable methodological quality and the results are largely inconsistent. There is emerging evidence associating differential expression of snoRNA genes in ASD. It is not currently possible to say whether the reports of differential expression in ncRNA may relate to ASD aetiology, a response to shared environmental factors linked to ASD such as sleep and nutrition, other molecular functions, human diversity, or chance findings. To improve our understanding of any potential association, we recommend improved and standardised methodologies and reporting of raw data. Further high-quality research is required to shine a light on possible associations, which may yet yield important information.

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.

The impact of paternal age on new mutations and disease in the next generation

Advanced paternal age is associated with an increased risk of fathering children with genetic disorders and other adverse reproductive consequences. However, the mechanisms underlying this phenomenon remain largely unexplored. In this review, we focus on the impact of paternal age on de novo mutations that are an important contributor to genetic disease and can be studied both indirectly through large-scale sequencing studies and directly in the tissue in which they predominantly arise-the aging testis. We discuss the recent data that have helped establish the origins and frequency of de novo mutations, and highlight experimental evidence about the close link between new mutations, parental age, and genetic disease. We then focus on a small group of rare genetic conditions, the so-called "paternal age effect" disorders that show a strong association between paternal age and disease prevalence, and discuss the underlying mechanism ("selfish selection") and implications of this process in more detail. More broadly, understanding the causes and consequences of paternal age on genetic risk has important implications both for individual couples and for public health advice given that the average age of fatherhood is steadily increasing in many developed nations.

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.

Two Genetic Mechanisms in Two Siblings with Intellectual Disability, Autism Spectrum Disorder, and Psychosis

Intellectual disability (ID) and autism spectrum disorder (ASD) are complex neurodevelopmental disorders with high heritability. To search for the genetic deficits in two siblings affected with ID and ASD in a family, we first performed a genome-wide copy number variation (CNV) analysis using chromosomal microarray analysis (CMA). We found a 3.7 Mb microdeletion at 22q13.3 in the younger sister. This de novo microdeletion resulted in the haploinsufficiency of SHANK3 and several nearby genes involved in neurodevelopment disorders. Hence, she was diagnosed with Phelan–McDermid syndrome (PMS, OMIM#606232). We further performed whole-genome sequencing (WGS) analysis in this family. We did not detect pathogenic mutations with significant impacts on the phenotypes of the elder brother. Instead, we identified several rare, likely pathogenic variants in seven genes implicated in neurodevelopmental disorders: KLHL17, TDO2, TRRAP, EIF3F, ATP10A, DICER1, and CDH15. These variants were transmitted from his unaffected parents, indicating these variants have only moderate clinical effects. We propose that these variants worked together and led to the clinical phenotypes in the elder brother. We also suggest that the combination of multiple genes with moderate effects is part of the genetic mechanism of neurodevelopmental disorders.

Non-Mendelian inheritance patterns and extreme deviation rates of CGG repeats in autism

As expansions of CGG short tandem repeats (STRs) are established as the genetic etiology of many neurodevelopmental disorders, we aimed to elucidate the inheritance patterns and role of CGG STRs in autism-spectrum disorder (ASD). By genotyping 6063 CGG STR loci in a large cohort of trios and quads with an ASD-affected proband, we determined an unprecedented rate of CGG repeat length deviation across a single generation. Although the concept of repeat length being linked to deviation rate was solidified, we show how shorter STRs display greater degrees of size variation. We observed that CGG STRs did not segregate by Mendelian principles but with a bias against longer repeats, which appeared to magnify as repeat length increased. Through logistic regression, we identified 19 genes that displayed significantly higher rates and degrees of CGG STR expansion within the ASD-affected probands (P < 1 × 10^-5). This study not only highlights novel repeat expansions that may play a role in ASD but also reinforces the hypothesis that CGG STRs are specifically linked to human cognition.