Revisiting tandem repeats in psychiatric disorders from perspectives of genetics, physiology, and brain evolution

Multi-omic quantitative trait loci link tandem repeat size variation to gene regulation in human brain

Tandem repeat (TR) size variation is implicated in ~50 neurological disorders, yet its impact on gene regulation in the human brain remains largely unknown. In the present study, we quantified the impact of TR size variation on brain gene regulation across distinct molecular phenotypes, based on 4,412 multi-omics samples from 1,597 donors, including 1,586 newly sequenced ones. We identified ~2.2 million TR molecular quantitative trait loci (TR-xQTLs), linking ~139,000 unique TRs to nearby molecular phenotypes, including many known disease-risk TRs, such as the G2C4 expansion in C9orf72 associated with amyotrophic lateral sclerosis. Fine-mapping revealed ~18,700 TRs as potential causal variants. Our in vitro experiments further confirmed the causal and independent regulatory effects of three TRs. Additional colocalization analysis indicated the potential causal role of TR variation in brain-related phenotypes, highlighted by a 3'-UTR TR in NUDT14 linked to cortical surface area and a TG repeat in PLEKHA1, associated with Alzheimer's disease.

Genetic modifiers of repeat expansion disorders

Repeat expansion disorders (REDs) are monogenic diseases caused by a sequence of repetitive DNA expanding above a pathogenic threshold. A common feature of the REDs is a strong genotype-phenotype correlation in which a major determinant of age at onset (AAO) and disease progression is the length of the inherited repeat tract. Over a disease-gene carrier's life, the length of the repeat can expand in somatic cells, through the process of somatic expansion which is hypothesised to drive disease progression. Despite being monogenic, individual REDs are phenotypically variable, and exploring what genetic modifying factors drive this phenotypic variability has illuminated key pathogenic mechanisms that are common to this group of diseases. Disease phenotypes are affected by the cognate gene in which the expansion is found, the location of the repeat sequence in coding or non-coding regions and by the presence of repeat sequence interruptions. Human genetic data, mouse models and in vitro models have implicated the disease-modifying effect of DNA repair pathways via the mechanisms of somatic mutation of the repeat tract. As such, developing an understanding of these pathways in the context of expanded repeats could lead to future disease-modifying therapies for REDs.

The role of tandem repeat expansions in brain disorders

The human genome contains numerous genetic polymorphisms contributing to different health and disease outcomes. Tandem repeat (TR) loci are highly polymorphic yet under-investigated in large genomic studies, which has prompted research efforts to identify novel variations and gain a deeper understanding of their role in human biology and disease outcomes. We summarize the current understanding of TRs and their implications for human health and disease, including an overview of the challenges encountered when conducting TR analyses and potential solutions to overcome these challenges. By shedding light on these issues, this article aims to contribute to a better understanding of the impact of TRs on the development of new disease treatments.

Identification of a psychiatric risk gene NISCH at 3p21.1 GWAS locus mediating dendritic spine morphogenesis and cognitive function

BACKGROUND

Schizophrenia and bipolar disorder (BD) are believed to share clinical symptoms, genetic risk, etiological factors, and pathogenic mechanisms. We previously reported that single nucleotide polymorphisms spanning chromosome 3p21.1 showed significant associations with both schizophrenia and BD, and a risk SNP rs2251219 was in linkage disequilibrium with a human specific Alu polymorphism rs71052682, which showed enhancer effects on transcriptional activities using luciferase reporter assays in U251 and U87MG cells.

METHODS

CRISPR/Cas9-directed genome editing, real-time quantitative PCR, and public Hi-C data were utilized to investigate the correlation between the Alu polymorphism rs71052682 and NISCH. Primary neuronal culture, immunofluorescence staining, co-immunoprecipitation, lentiviral vector production, intracranial stereotaxic injection, behavioral assessment, and drug treatment were used to examine the physiological impacts of Nischarin (encoded by NISCH).

RESULTS

Deleting the Alu sequence in U251 and U87MG cells reduced mRNA expression of NISCH, the gene locates 180 kb from rs71052682, and Hi-C data in brain tissues confirmed the extensive chromatin contacts. These data suggested that the genetic risk of schizophrenia and BD predicted elevated NISCH expression, which was also consistent with the observed higher NISCH mRNA levels in the brain tissues from psychiatric patients compared with controls. We then found that overexpression of NISCH resulted in a significantly decreased density of mushroom dendritic spines with a simultaneously increased density of thin dendritic spines in primary cultured neurons. Intriguingly, elevated expression of this gene in mice also led to impaired spatial working memory in the Y-maze. Given that Nischarin is the target of anti-hypertensive agents clonidine and tizanidine, which have shown therapeutic effects in patients with schizophrenia and patients with BD in preliminary clinical trials, we demonstrated that treatment with those antihypertensive drugs could reduce NISCH mRNA expression and rescue the impaired working memory in mice.

CONCLUSIONS

We identify a psychiatric risk gene NISCH at 3p21.1 GWAS locus influencing dendritic spine morphogenesis and cognitive function, and Nischarin may have potentials for future therapeutic development.

A novel hypervariable variable number tandem repeat in the dopamine transporter gene (SLC6A3)

The dopamine transporter gene, SLC6A3, has received substantial attention in genetic association studies of various phenotypes. Although some variable number tandem repeats (VNTRs) present in SLC6A3 have been tested in genetic association studies, results have not been consistent. VNTRs in SLC6A3 that have not been examined genetically were characterized. The Tandem Repeat Annotation Library was used to characterize the VNTRs of 64 unrelated long-read haplotype-phased SLC6A3 sequences. Sequence similarity of each repeat unit of the five VNTRs is reported, along with the correlations of SNP-SNP, SNP-VNTR, and VNTR-VNTR alleles across the gene. One of these VNTRs is a novel hyper-VNTR (hyVNTR) in intron 8 of SLC6A3, which contains a range of 3.4-133.4 repeat copies and has a consensus sequence length of 38 bp, with 82% G+C content. The 38-base repeat was predicted to form G-quadruplexes in silico and was confirmed by circular dichroism spectroscopy. In addition, this hyVNTR contains multiple putative binding sites for PRDM9, which, in combination with low levels of linkage disequilibrium around the hyVNTR, suggests it might be a recombination hotspot.

Mutation and selection processes regulating short tandem repeats give rise to genetic and phenotypic diversity across species

Short tandem repeats (STRs) are units of 1-6 bp that repeat in a tandem fashion in DNA. Along with single nucleotide polymorphisms and large structural variations, they are among the major genomic variants underlying genetic, and likely phenotypic, divergence. STRs experience mutation rates that are orders of magnitude higher than other well-studied genotypic variants. Frequent copy number changes result in a wide range of alleles, and provide unique opportunities for modulating complex phenotypes through variation in repeat length. While classical studies have identified key roles of individual STR loci, the advent of improved sequencing technology, high-quality genome assemblies for diverse species, and bioinformatics methods for genome-wide STR analysis now enable more systematic study of STR variation across wide evolutionary ranges. In this review, we explore mutation and selection processes that affect STR copy number evolution, and how these processes give rise to varying STR patterns both within and across species. Finally, we review recent examples of functional and adaptive changes linked to STRs.

Beyond monoamines: I. Novel targets and radiotracers for Positron emission tomography imaging in psychiatric disorders

With the emergence of positron emission tomography (PET) in the late 1970s, psychiatry had access to a tool capable of non-invasive assessment of human brain function. Early applications in psychiatry focused on identifying characteristic brain blood flow and metabolic derangements using radiotracers such as [¹⁵ O]H₂ O and [¹⁸ F]FDG. Despite the success of these techniques, it became apparent that more specific probes were needed to understand the neurochemical bases of psychiatric disorders. The first neurochemical PET imaging probes targeted sites of action of neuroleptic (dopamine D₂ receptors) and psychoactive (serotonin receptors) drugs. Based on the centrality of monoamine dysfunction in psychiatric disorders and the measured success of monoamine-enhancing drugs in treating them, the next 30 years witnessed the development of an armamentarium of PET radiopharmaceuticals and imaging methodologies for studying monoamines. Continued development of monoamine-enhancing drugs over this time however was less successful, realizing only modest gains in efficacy and tolerability. As patent protection for many widely prescribed and profitable psychiatric drugs lapsed, drug development pipelines shifted away from monoamines in search of novel targets with the promises of improved efficacy, or abandoned altogether. Over this period, PET radiopharmaceutical development activities closely paralleled drug development priorities resulting in the development of new PET imaging agents for non-monoamine targets. Part one of this review will briefly survey novel PET imaging targets with relevance to the field of psychiatry, which include the metabotropic glutamate receptor type 5 (mGluR5), purinergic P₂ X₇ receptor, type 1 cannabinoid receptor (CB₁ ), phosphodiesterase 10A (PDE10A), and describe radiotracers developed for these and other targets that have matured to human subject investigations. Current limitations of the targets and techniques will also be discussed.

Structural variation and eQTL analysis in two experimental populations of chickens divergently selected for feather-pecking behavior

Feather pecking (FP) is a damaging nonaggressive behavior in laying hens with a heritable component. Its occurrence has been linked to the immune system, the circadian clock, and foraging behavior. Furthermore, dysregulation of miRNA biogenesis, disturbance of the gamma-aminobutyric acid (GABAergic) system, as well as neurodevelopmental deficiencies are currently under debate as factors influencing the propensity for FP behavior. Past studies, which focused on the dissection of the genetic factors involved in FP, relied on single nucleotide polymorphisms (SNPs) and short insertions and deletions < 50 bp (InDels). These variant classes only represent a certain fraction of the genetic variation of an organism. Hence, we reanalyzed whole-genome sequencing data from two experimental populations, which have been divergently selected for FP behavior for over more than 15 generations, performed variant calling for structural variants (SVs) as well as tandem repeats (TRs), and jointly analyzed the data with SNPs and InDels. Genotype imputation and subsequent genome-wide association studies, in combination with expression quantitative trait loci analysis, led to the discovery of multiple variants influencing the GABAergic system. These include a significantly associated TR downstream of the GABA receptor subunit beta-3 (GABRB3) gene, two microRNAs targeting several GABA receptor genes, and dystrophin (DMD), a direct regulator of GABA receptor clustering. Furthermore, we found the transcription factor ETV1 to be associated with the differential expression of 23 genes, which points toward a role of ETV1, together with SMAD4 and KLF14, in the disturbed neurodevelopment of high-feather pecking chickens.

Structural variants and tandem repeats in the founder individuals of four F2 pig crosses and implications to F2 GWAS results

BACKGROUND

Structural variants and tandem repeats are relevant sources of genomic variation that are not routinely analyzed in genome wide association studies mainly due to challenging identification and genotyping. Here, we profiled these variants via state-of-the-art strategies in the founder animals of four F2 pig crosses using whole-genome sequence data (20x coverage). The variants were compared at a founder level with the commonly screened SNPs and small indels. At the F2 level, we carried out an association study using imputed structural variants and tandem repeats with four growth and carcass traits followed by a comparison with a previously conducted SNPs and small indels based association study.

RESULTS

A total of 13,201 high confidence structural variants and 103,730 polymorphic tandem repeats (with a repeat length of 2-20 bp) were profiled in the founders. We observed a moderate to high (r from 0.48 to 0.57) level of co-localization between SNPs or small indels and structural variants or tandem repeats. In the association step 56.56% of the significant variants were not in high LD with significantly associated SNPs and small indels identified for the same traits in the earlier study and thus presumably not tagged in case of a standard association study. For the four growth and carcass traits investigated, many of the already proposed candidate genes in our previous studies were confirmed and additional ones were identified. Interestingly, a common pattern on how structural variants or tandem repeats regulate the phenotypic traits emerged. Many of the significant variants were embedded or nearby long non-coding RNAs drawing attention to their functional importance. Through which specific mechanisms the identified long non-coding RNAs and their associated structural variants or tandem repeats contribute to quantitative trait variation will need further investigation.

CONCLUSIONS

The current study provides insights into the characteristics of structural variants and tandem repeats and their role in association studies. A systematic incorporation of these variants into genome wide association studies is advised. While not of immediate interest for genomic prediction purposes, this will be particularly beneficial for elucidating biological mechanisms driving the complex trait variation.

L1 Retrotransposons: A Potential Endogenous Regulator for Schizophrenia

The long interspersed nuclear elements 1 (LINE-1/L1s) are the only active autonomous retrotransposons found in humans which can integrate anywhere in the human genome. They can expand the genome and thus bring good or bad effects to the host cells which really depends on their integration site and associated polymorphism. LINE-1 retrotransposition has been found participating in various neurological disorders such as autism spectrum disorder, Alzheimer’s disease, major depression disorder, post-traumatic stress disorder and schizophrenia. Despite the recent progress, the roles and pathological mechanism of LINE-1 retrotransposition in schizophrenia and its heritable risks, particularly, contribution to “missing heritability” are yet to be determined. Therefore, this review focuses on the potentially etiological roles of L1s in the development of schizophrenia, possible therapeutic choices and unaddressed questions in order to shed lights on the future research.

Longitudinal Influences of DRD4 Polymorphism and Early Maternal Caregiving on Personality Development and Problem Behavior in Middle Childhood and Adolescence

Most studies examining gene-environment effects on self-regulation focus on outcomes early childhood or adulthood. However, only a few studies investigate longitudinal effects during middle childhood and adolescence and compare two domains of early caregiving. In a longitudinal follow-up with a sample of N = 87, we studied the effects of differences in the DRD4 tandem repeat polymorphisms and two domains of early maternal caregiving quality on children's personality development using Block's California Child Q-Set (CCQ) at age six and age 12 and on problem behavior at ages six and seven. Early maternal regulation quality predicted later ego-resiliency and aggressiveness. In addition, significant gene-environment interactions revealed that children with the 7+ DRD4 tandem repeat polymorphism and poor maternal regulation quality in infancy showed lower scores in ego-resiliency and higher scores in ego-undercontrol and CCQ aggressiveness. In contrast, children who had experienced effective maternal regulation in infancy showed a comparable level in personality traits and problem behavior as the DRD4 7- group independent of the levels of maternal regulatory behavior. Similarly, longitudinal caregiving × DRD4 interactions were found for behavior problems in middle childhood, especially for oppositional-aggression, inattentive-hyperactivity, and social competence. Early caregiving effects were only found for maternal regulation quality, but not for maternal responsiveness. Effective early maternal regulation in infancy can moderate the negative effect of DRD4 7+ on children's self-regulation in middle childhood and adolescence. However, maternal responsiveness has no comparable effects. It seems relevant to consider several dimensions of early caregiving and to also measure the environment in more detail in gene-environment studies.