Characterization of a Human-Specific Tandem Repeat Associated with Bipolar Disorder and Schizophrenia

Myelin Formation by Oligodendrocytes Is Enhanced Through Laminin-411 and Its Derived Peptide

In the central nervous system, oligodendrocytes (OLs) form myelin sheaths that accomplish the efficient transmission of nerve conduction for optimal motor and cognitive functions. OL development and differentiation are regulated by a variety of molecules, including extracellular matrix (ECM) proteins. ECM proteins are also useful as substrates for OL culture. However, the functions of ECM proteins in OL development and myelination remain unclear, and only a limited number of ECM proteins have been characterized and used in in vitro experiments. Here, we investigated the expression and function of laminin (LM) isoforms in OL differentiation and myelination. We found that LM α1, α2, and α4 chains were expressed around blood vessels at the stage of myelination in mice. Functional analyses using recombinant proteins of LM isoforms containing α1, α2, and α4 chains revealed that LM411 and LM411E8, the integrin binding domain of LM411, possessed significant activities in myelin membrane formation of OLs. Furthermore, the peptide A4G47 derived from LM411E8 promoted the activity, which provides evidence of the first peptide in OL myelin formation from ECM proteins. Our findings facilitate a better understanding of ECM functions in OL biology and the development of a new material in OL myelination.

Fast sequence alignment for centromeres with RaMA

The release of the first draft of the human pangenome has revolutionized genomic research by enabling access to complex regions like centromeres, composed of extra-long tandem repeats (ETRs). However, a significant gap remains as current methodologies are inadequate for producing sequence alignments that effectively capture genetic events within ETRs, highlighting a pressing need for improved alignment tools. Inspired by UniAligner, we developed a rare match aligner (RaMA), using rare matches as anchors and two-piece affine gap cost to generate complete pairwise alignment that better captures genetic evolution. RaMA also employs parallel computing and the wavefront algorithm to accelerate anchor discovery and sequence alignment, achieving up to 13.66 times faster processing using only 11% of UniAligner's memory. Downstream analysis of simulated data and the CHM13 and CHM1 higher-order repeat (HOR) arrays demonstrates that RaMA achieves more accurate alignments, effectively capturing true HOR structures. RaMA also introduces two methods for defining reliable alignment regions, further refining and enhancing the accuracy of centromeric alignment statistics.

Multisample motif discovery and visualization for tandem repeats

Tandem repeats (TRs) occupy a significant portion of the human genome and are a source of polymorphisms due to variations in sizes and motif compositions. Some of these variations have been associated with various neuropathological disorders, highlighting the clinical importance of assessing the motif structure of TRs. Moreover, assessing the TR motif variation can offer valuable insights into evolutionary dynamics and population structure. Previously, characterizations of TRs were limited by short-read sequencing technology, which lacks the ability to accurately capture the full TR sequences. As long-read sequencing becomes more accessible and can capture the full complexity of TRs, there is now also a need for tools to characterize and analyze TRs using long-read data across multiple samples. In this study, we present MotifScope, a novel algorithm for the characterization and visualization of TRs based on a de novo k-mer approach for motif discovery. Comparative analysis against established tools reveals that MotifScope can identify a greater number of motifs and more accurately represent the underlying repeat sequences. Moreover, MotifScope has been specifically designed to enable motif composition comparisons across assemblies of different individuals, as well as across long-read sequencing reads within an individual, through combined motif discovery and sequence alignment. We showcase potential applications of MotifScope in diverse fields, including population genetics, clinical settings, and forensic analyses.

Analysis of targeted and whole genome sequencing of PacBio HiFi reads for a comprehensive genotyping of gene-proximal and phenotype-associated Variable Number Tandem Repeats

Variable Number Tandem repeats (VNTRs) refer to repeating motifs of size greater than five bp. VNTRs are an important source of genetic variation, and have been associated with multiple Mendelian and complex phenotypes. However, the highly repetitive structures require reads to span the region for accurate genotyping. Pacific Biosciences HiFi sequencing spans large regions and is highly accurate but relatively expensive. Therefore, targeted sequencing approaches coupled with long-read sequencing have been proposed to improve efficiency and throughput. In this paper, we systematically explored the trade-off between targeted and whole genome HiFi sequencing for genotyping VNTRs. We curated a set of 10 , 787 gene-proximal (G-)VNTRs, and 48 phenotype-associated (P-)VNTRs of interest. Illumina reads only spanned 46% of the G-VNTRs and 71% of P-VNTRs, motivating the use of HiFi sequencing. We performed targeted sequencing with hybridization by designing custom probes for 9,999 VNTRs and sequenced 8 samples using HiFi and Illumina sequencing, followed by adVNTR genotyping. We compared these results against HiFi whole genome sequencing (WGS) data from 28 samples in the Human Pangenome Reference Consortium (HPRC). With the targeted approach only 4,091 (41%) G-VNTRs and only 4 (8%) of P-VNTRs were spanned with at least 15 reads. A smaller subset of 3,579 (36%) G-VNTRs had higher median coverage of at least 63 spanning reads. The spanning behavior was consistent across all 8 samples. Among 5,638 VNTRs with low-coverage ( < 15), 67% were located within GC-rich regions ( > 60%). In contrast, the 40X WGS HiFi dataset spanned 98% of all VNTRs and 49 (98%) of P-VNTRs with at least 15 spanning reads, albeit with lower coverage. Spanning reads were sufficient for accurate genotyping in both cases. Our findings demonstrate that targeted sequencing provides consistently high coverage for a small subset of low-GC VNTRs, but WGS is more effective for broad and sufficient sampling of a large number of VNTRs.

Aberrant Short Tandem Repeats: Pathogenicity, Mechanisms, Detection, and Roles in Neuropsychiatric Disorders

Short tandem repeat (STR) sequences are highly variable DNA segments that significantly contribute to human neurodegenerative disorders, highlighting their crucial role in neuropsychiatric conditions. This article examines the pathogenicity of abnormal STRs and classifies tandem repeat expansion disorders(TREDs), emphasizing their genetic characteristics, mechanisms of action, detection methods, and associated animal models. STR expansions exhibit complex genetic patterns that affect the age of onset and symptom severity. These expansions disrupt gene function through mechanisms such as gene silencing, toxic gain-of-function mutations leading to RNA and protein toxicity, and the generation of toxic peptides via repeat-associated non-AUG (RAN) translation. Advances in sequencing technologies-from traditional PCR and Southern blotting to next-generation and long-read sequencing-have enhanced the accuracy of STR variation detection. Research utilizing these technologies has linked STR expansions to a range of neuropsychiatric disorders, including autism spectrum disorders and schizophrenia, highlighting their contribution to disease risk and phenotypic expression through effects on genes involved in neurodevelopment, synaptic function, and neuronal signaling. Therefore, further investigation is essential to elucidate the intricate interplay between STRs and neuropsychiatric diseases, paving the way for improved diagnostic and therapeutic strategies.

Approximating edit distances between complex tandem repeats efficiently

MOTIVATION

Extended tandem repeats (TRs) have been associated with 60 or more diseases over the past 30 years. Although most TRs have single repeat units (or motifs), complex TRs with different units have recently been correlated with some brain disorders. Of note, a population-scale analysis shows that complex TRs at one locus can be divergent, and different units are often expanded between individuals. To understand the evolution of high TR diversity, it is informative to visualize a phylogenetic tree. To do this, we need to measure the edit distance between pairs of complex TRs by considering duplication and contraction of units created by replication slippage. However, traditional rigorous algorithms for this purpose are computationally expensive.

RESULTS

We here propose an efficient heuristic algorithm to estimate the edit distance with duplication and contraction of units (EDDC, for short). We select a set of frequent units that occur in given complex TRs, encode each unit as a single symbol, compress a TR into an optimal series of unit symbols that partially matches the original TR with the minimum Levenshtein distance, and estimate the EDDC between a pair of complex TRs from their compressed forms. Using substantial synthetic benchmark datasets, we demonstrate that the estimated EDDC is highly correlated with the accurate EDDC, with a Pearson correlation coefficient of >0.983, while the heuristic algorithm achieves orders of magnitude performance speedup.

AVAILABILITY AND IMPLEMENTATION

The software program hEDDC that implements the proposed algorithm is available at https://github.com/Ricky-pon/hEDDC (DOI: 10.5281/zenodo.14732958).

ScatTR: Estimating the Size of Long Tandem Repeat Expansions from Short-Reads

Tandem repeats (TRs) are sequences of DNA where two or more base pairs are repeated back-to-back at specific locations in the genome. The expansions of TRs are implicated in over 50 conditions, including Friedreich's ataxia, autism, and cancer. However, accurately measuring the copy number of TRs is challenging, especially when their expansions are larger than the fragment sizes used in standard short-read genome sequencing. Here we introduce ScatTR, a novel computational method that leverages a maximum likelihood framework to estimate the copy number of large TR expansions from short-read sequencing data. ScatTR calculates the likelihood of different alignments between sequencing reads and reference sequences that represent various TR lengths and employs a Monte Carlo technique to find the best match. In simulated data, ScatTR outperforms state-of-the-art methods, particularly for TRs with longer motifs and those with lengths that greatly exceed typical sequencing fragment sizes. When applied to data from the 1000 Genomes Project, ScatTR detected potential large TR expansions that other methods missed, highlighting its ability to better identify genome-wide characterization of TR variation. ScatTR can be accessed via: https://github.com/g2lab/scattr .

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.

Novelty versus innovation of gene regulatory elements in human evolution and disease

It is not currently understood how much of human evolution is due to modifying existing functional elements in the genome versus forging novel elements from nonfunctional DNA. Many early experiments that aimed to assign genetic changes on the human lineage to their resulting phenotypic change have focused on mutations that modify existing elements. However, a number of recent studies have highlighted the potential ease and importance of forging novel gene regulatory elements from nonfunctional sequences on the human lineage. In this review, we distinguish gene regulatory element novelty from innovation. We propose definitions for these terms and emphasize their importance in studying the genetic basis of human uniqueness. We discuss why the forging of novel regulatory elements may have been less emphasized during the previous decades, and why novel regulatory elements are likely to play a significant role in both human adaptation and disease.

Genome-wide investigation of VNTR motif polymorphisms in 8,222 genomes: Implications for biological regulation and human traits

Variable number tandem repeat (VNTR) is a pervasive and highly mutable genetic feature that varies in both length and repeat sequence. Despite the well-studied copy-number variants, the functional impacts of repeat motif polymorphisms remain unknown. Here, we present the largest genome-wide VNTR polymorphism map to date, with over 2.5 million VNTR length polymorphisms (VNTR-LPs) and over 11 million VNTR motif polymorphisms (VNTR-MPs) detected in 8,222 high-coverage genomes. Leveraging the large-scale NyuWa cohort, we identified 2,982,456 (31.8%) NyuWa-specific VNTR-MPs, of which 95.3% were rare. Moreover, we found 1,937 out of 38,685 VNTRs that were associated with gene expression through VNTR-MPs in lymphoblastoid cell lines. Specifically, we clarified that the expansion of a likely causal motif could upregulate gene expression by improving the binding concentration of PU.1. We also explored the potential impacts of VNTR polymorphisms on phenotypic differentiation and disease susceptibility. This study expands our knowledge of VNTR-MPs and their functional implications.

How our brains are built: emerging approaches to understand human-specific features

Understanding what makes us uniquely human is a long-standing question permeating fields from genomics, neuroscience, and developmental biology to medicine. The discovery of human-specific genomic sequences has enabled a new understanding of the molecular features of human brain evolution. Advances in sequencing, computational, and in vitro screening approaches collectively reveal new roles of uniquely human sequences in regulating gene expression. Here, we review the landscape of human-specific loci and describe how emerging technologies are being used to understand their molecular functions and impact on brain development. We describe current challenges in the field and the potential of integrating new hypotheses and approaches to propel our understanding of the human brain.

Accurate Detection of Tandem Repeats from Error-Prone Sequences with EquiRep

A tandem repeat is a sequence of nucleotides that occurs as multiple contiguous and near-identical copies positioned next to each other. These repeats play critical roles in genetic diversity, gene regulation, and are strongly linked to various neurological and developmental disorders. While several methods exist for detecting tandem repeats, they often exhibit low accuracy when the repeat unit length increases or the number of copies is low. Furthermore, methods capable of handling highly mutated sequences remain scarce, highlighting a significant opportunity for improvement. We introduce EquiRep, a tool for accurate detection of tandem repeats from erroneous sequences. EquiRep estimates the likelihood of positions originating from the same position in the unit by self-alignment followed by a novel approach that refines the estimation. The built equivalent classes and the consecutive position information will be then used to build a weighted graph, and the cycle in this graph with maximum bottleneck weight while covering most nucleotide positions will be identified to reconstruct the repeat unit. We test EquiRep on simulated and real HOR and RCA datasets where it consistently outperforms or is comparable to state-of-the-art methods. EquiRep is robust to sequencing errors, and is able to make better predictions for long units and low frequencies which underscores its broad usability for studying tandem repeats.

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.

Key Roles of CACNA1C/Cav1.2 and CALB1/Calbindin in Prefrontal Neurons Altered in Cognitive Disorders

Importance

The risk of mental disorders is consistently associated with variants in CACNA1C (L-type calcium channel Cav1.2) but it is not known why these channels are critical to cognition, and whether they affect the layer III pyramidal cells in the dorsolateral prefrontal cortex that are especially vulnerable in cognitive disorders.

Objective

To examine the molecular mechanisms expressed in layer III pyramidal cells in primate dorsolateral prefrontal cortices.

Design, Setting, and Participants

The design included transcriptomic analyses from human and macaque dorsolateral prefrontal cortex, and connectivity, protein expression, physiology, and cognitive behavior in macaques. The research was performed in academic laboratories at Yale, Harvard, Princeton, and the University of Pittsburgh. As dorsolateral prefrontal cortex only exists in primates, the work evaluated humans and macaques.

Main Outcomes and Measures

Outcome measures included transcriptomic signatures of human and macaque pyramidal cells, protein expression and interactions in layer III macaque pyramidal cells using light and electron microscopy, changes in neuronal firing during spatial working memory, and working memory performance following pharmacological treatments.

Results

Layer III pyramidal cells in dorsolateral prefrontal cortex coexpress a constellation of calcium-related proteins, delineated by CALB1 (calbindin), and high levels of CACNA1C (Cav1.2), GRIN2B (NMDA receptor GluN2B), and KCNN3 (SK3 potassium channel), concentrated in dendritic spines near the calcium-storing smooth endoplasmic reticulum. L-type calcium channels influenced neuronal firing needed for working memory, where either blockade or increased drive by β1-adrenoceptors, reduced neuronal firing by a mean (SD) 37.3% (5.5%) or 40% (6.3%), respectively, the latter via SK potassium channel opening. An L-type calcium channel blocker or β1-adrenoceptor antagonist protected working memory from stress.

Conclusions and Relevance

The layer III pyramidal cells in the dorsolateral prefrontal cortex especially vulnerable in cognitive disorders differentially express calbindin and a constellation of calcium-related proteins including L-type calcium channels Cav1.2 (CACNA1C), GluN2B-NMDA receptors (GRIN2B), and SK3 potassium channels (KCNN3), which influence memory-related neuronal firing. The finding that either inadequate or excessive L-type calcium channel activation reduced neuronal firing explains why either loss- or gain-of-function variants in CACNA1C were associated with increased risk of cognitive disorders. The selective expression of calbindin in these pyramidal cells highlights the importance of regulatory mechanisms in neurons with high calcium signaling, consistent with Alzheimer tau pathology emerging when calbindin is lost with age and/or inflammation.

Rare variation in non-coding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of non-coding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of non-coding regions: human accelerated regions (HARs), which show signatures of positive selection in humans; experimentally validated neural VISTA enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole-genome analysis of >16,600 samples and >4,900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly contribute, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in VEs near OTX1 and SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved non-coding regions in ASD risk and suggest potential mechanisms of how regulatory changes can modulate social behavior.

Comprehensive assessment of long-read sequencing platforms and calling algorithms for detection of copy number variation

Copy number variations (CNVs) play pivotal roles in disease susceptibility and have been intensively investigated in human disease studies. Long-read sequencing technologies offer opportunities for comprehensive structural variation (SV) detection, and numerous methodologies have been developed recently. Consequently, there is a pressing need to assess these methods and aid researchers in selecting appropriate techniques for CNV detection using long-read sequencing. Hence, we conducted an evaluation of eight CNV calling methods across 22 datasets from nine publicly available samples and 15 simulated datasets, covering multiple sequencing platforms. The overall performance of CNV callers varied substantially and was influenced by the input dataset type, sequencing depth, and CNV type, among others. Specifically, the PacBio CCS sequencing platform outperformed PacBio CLR and Nanopore platforms regarding CNV detection recall rates. A sequencing depth of 10x demonstrated the capability to identify 85% of the CNVs detected in a 50x dataset. Moreover, deletions were more generally detectable than duplications. Among the eight benchmarked methods, cuteSV, Delly, pbsv, and Sniffles2 demonstrated superior accuracy, while SVIM exhibited high recall rates.

Structural characterization of a polymorphic repeat at the CACNA1C schizophrenia locus

Genetic variation within intron 3 of the CACNA1C calcium channel gene is associated with schizophrenia and bipolar disorder, but analysis of the causal variants and their effect is complicated by a nearby variable-number tandem repeat (VNTR). Here, we used 155 long-read genome assemblies from 78 diverse individuals to delineate the structure and population variability of the CACNA1C intron 3 VNTR. We categorized VNTR sequences into 7 Types of structural alleles using sequence differences among repeat units. Only 12 repeat units at the 5' end of the VNTR were shared across most Types, but several Types were related through a series of large and small duplications. The most diverged Types were rare and present only in individuals with African ancestry, but the multiallelic structural polymorphism Variable Region 2 was present across populations at different frequencies, consistent with expansion of the VNTR preceding the emergence of early hominins. VR2 was in complete linkage disequilibrium with fine-mapped schizophrenia variants (SNPs) from genome-wide association studies (GWAS). This risk haplotype was associated with decreased CACNA1C gene expression in brain tissues profiled by the GTEx project. Our work suggests that sequence variation within a human-specific VNTR affects gene expression, and provides a detailed characterization of new alleles at a flagship neuropsychiatric locus.

Rediscovering tandem repeat variation in schizophrenia: challenges and opportunities

Tandem repeats (TRs) are prevalent throughout the genome, constituting at least 3% of the genome, and often highly polymorphic. The high mutation rate of TRs, which can be orders of magnitude higher than single-nucleotide polymorphisms and indels, indicates that they are likely to make significant contributions to phenotypic variation, yet their contribution to schizophrenia has been largely ignored by recent genome-wide association studies (GWAS). Tandem repeat expansions are already known causative factors for over 50 disorders, while common tandem repeat variation is increasingly being identified as significantly associated with complex disease and gene regulation. The current review summarizes key background concepts of tandem repeat variation as pertains to disease risk, elucidating their potential for schizophrenia association. An overview of next-generation sequencing-based methods that may be applied for TR genome-wide identification is provided, and some key methodological challenges in TR analyses are delineated.

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.

Advances in the discovery and analyses of human tandem repeats

Long-read sequencing platforms provide unparalleled access to the structure and composition of all classes of tandemly repeated DNA from STRs to satellite arrays. This review summarizes our current understanding of their organization within the human genome, their importance with respect to disease, as well as the advances and challenges in understanding their genetic diversity and functional effects. Novel computational methods are being developed to visualize and associate these complex patterns of human variation with disease, expression, and epigenetic differences. We predict accurate characterization of this repeat-rich form of human variation will become increasingly relevant to both basic and clinical human genetics.

Effect of rs1108580 of DBH and rs1006737 of CACNA1C on Cognition and Tardive Dyskinesia in a North Indian Schizophrenia Cohort

Genetic perturbations in dopamine neurotransmission and calcium signaling pathways are implicated in the etiology of schizophrenia. We aimed to test the association of a functional splice variant each in Dopamine β-Hydroxylase (DBH; rs1108580) and Calcium voltage-gated channel subunit alpha1 C (CACNA1C; rs1006737) genes in these pathways with schizophrenia (506 cases, 443 controls); Abnormal Involuntary Movement Scale (AIMS) scores in subjects assessed for tardive dyskinesia (76 TD-positive, 95 TD-negative) and Penn Computerized Neurocognitive Battery (PennCNB) scores (334 cases, 234 controls). The effect of smoking status and SNP genotypes on AIMS scores were assessed using ANOVA; health status and SNP genotypes on three performance functions of PennCNB cognitive domains were assessed by ANCOVA with age and sex as covariates. Association with Positive and Negative Syndrome Scale (PANSS) scores in the TD cohort and cognitive scores in healthy controls of the cognition cohort were tested by linear regression. None of the markers were associated with schizophrenia. Smoking status [F(2, 139) = 10.6; p = 5 × 10-5], rs1006737 [F(2, 139) = 7.1; p = 0.001], TD status*smoking [F(2, 139) = 8.0; p = 5.0 × 10-4] and smoking status*rs1006737 [F(4, 139) = 2.7; p = 0.03] had an effect on AIMS score. Furthermore, rs1006737 was associated with orofacial [F(2, 139) = 4.6; p = 0.01] and limb-truncal TD [(F(2, 139) = 3.8; p = 0.02]. Main effect of rs1108580 on working memoryprocessing speed [F(2, 544) = 3.8; p = 0.03] and rs1006737 on spatial abilityefficiency [F(1, 550) = 9.4; p = 0.02] was identified. Health status*rs1006737 interaction had an effect on spatial memoryprocessing speed [F(1, 550) = 6.9; p = 0.01]. Allelic/genotypic association (p = 0.01/0.03) of rs1006737 with disorganized/concrete factor and allelic association of rs1108580 (p = 0.04) with a depressive factor of PANSS was observed in the TD-negative subcohort. Allelic association of rs1006737 with sensorimotor dexterityaccuracy (p = 0.03), attentionefficiency (p = 0.05), and spatial abilityefficiency (p = 0.02); allelic association of rs1108580 with face memoryaccuracy (p = 0.05) and emotionefficiency (p = 0.05); and allelic/genotypic association with emotionaccuracy (p = 0.003/0.009) were observed in healthy controls of the cognition cohort. These association findings may have direct implications for personalized medicine and cognitive remediation.

Human-specific genetics: new tools to explore the molecular and cellular basis of human evolution

Our ancestors acquired morphological, cognitive and metabolic modifications that enabled humans to colonize diverse habitats, develop extraordinary technologies and reshape the biosphere. Understanding the genetic, developmental and molecular bases for these changes will provide insights into how we became human. Connecting human-specific genetic changes to species differences has been challenging owing to an abundance of low-effect size genetic changes, limited descriptions of phenotypic differences across development at the level of cell types and lack of experimental models. Emerging approaches for single-cell sequencing, genetic manipulation and stem cell culture now support descriptive and functional studies in defined cell types with a human or ape genetic background. In this Review, we describe how the sequencing of genomes from modern and archaic hominins, great apes and other primates is revealing human-specific genetic changes and how new molecular and cellular approaches - including cell atlases and organoids - are enabling exploration of the candidate causal factors that underlie human-specific traits.

Rare variation in noncoding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of noncoding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of noncoding regions: Human Accelerated Regions (HARs), which show signatures of positive selection in humans; experimentally validated neural Vista Enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole genome analysis of >16,600 samples and >4900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in a VE near SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved noncoding regions in ASD risk and suggest potential mechanisms of how changes in regulatory regions can modulate social behavior.

A landscape of complex tandem repeats within individual human genomes

Markedly expanded tandem repeats (TRs) have been correlated with ~60 diseases. TR diversity has been considered a clue toward understanding missing heritability. However, haplotype-resolved long TRs remain mostly hidden or blacked out because their complex structures (TRs composed of various units and minisatellites containing >10-bp units) make them difficult to determine accurately with existing methods. Here, using a high-precision algorithm to determine complex TR structures from long, accurate reads of PacBio HiFi, an investigation of 270 Japanese control samples yields several genome-wide findings. Approximately 322,000 TRs are difficult to impute from the surrounding single-nucleotide variants. Greater genetic divergence of TR loci is significantly correlated with more events of younger replication slippage. Complex TRs are more abundant than single-unit TRs, and a tendency for complex TRs to consist of <10-bp units and single-unit TRs to be minisatellites is statistically significant at loci with ≥500-bp TRs. Of note, 8909 loci with extended TRs (>100b longer than the mode) contain several known disease-associated TRs and are considered candidates for association with disorders. Overall, complex TRs and minisatellites are found to be abundant and diverse, even in genetically small Japanese populations, yielding insights into the landscape of long TRs.

vamos: variable-number tandem repeats annotation using efficient motif sets

Roughly 3% of the human genome is composed of variable-number tandem repeats (VNTRs): arrays of motifs at least six bases. These loci are highly polymorphic, yet current approaches that define and merge variants based on alignment breakpoints do not capture their full diversity. Here we present a method vamos: VNTR Annotation using efficient Motif Sets that instead annotates VNTR using repeat composition under different levels of motif diversity. Using vamos we estimate 7.4-16.7 alleles per locus when applied to 74 haplotype-resolved human assemblies, compared to breakpoint-based approaches that estimate 4.0-5.5 alleles per locus.

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.

Comparative physiological anthropogeny: exploring molecular underpinnings of distinctly human phenotypes

Anthropogeny is a classic term encompassing transdisciplinary investigations of the origins of the human species. Comparative anthropogeny is a systematic comparison of humans and other living nonhuman hominids (so-called "great apes"), aiming to identify distinctly human features in health and disease, with the overall goal of explaining human origins. We begin with a historical perspective, briefly describing how the field progressed from the earliest evolutionary insights to the current emphasis on in-depth molecular and genomic investigations of "human-specific" biology and an increased appreciation for cultural impacts on human biology. While many such genetic differences between humans and other hominids have been revealed over the last two decades, this information remains insufficient to explain the most distinctive phenotypic traits distinguishing humans from other living hominids. Here we undertake a complementary approach of "comparative physiological anthropogeny," along the lines of the preclinical medical curriculum, i.e., beginning with anatomy and considering each physiological system and in each case considering genetic and molecular components that are relevant. What is ultimately needed is a systematic comparative approach at all levels from molecular to physiological to sociocultural, building networks of related information, drawing inferences, and generating testable hypotheses. The concluding section will touch on distinctive considerations in the study of human evolution, including the importance of gene-culture interactions.

TRviz: a Python library for decomposing and visualizing tandem repeat sequences

Summary

TRviz is an open-source Python library for decomposing, encoding, aligning and visualizing tandem repeat (TR) sequences. TRviz takes a collection of alleles (TR containing sequences) and one or more motifs as input and generates a plot showing the motif composition of the TR sequences.

Availability and implementation

TRviz is an open-source Python library and freely available at https://github.com/Jong-hun-Park/trviz. Detailed documentation is available at https://trviz.readthedocs.io.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Interplay Between Polymorphic Short Tandem Repeats and Gene Expression Variation in Caenorhabditis elegans

Short tandem repeats (STRs) have orders of magnitude higher mutation rates than single nucleotide variants (SNVs) and have been proposed to accelerate evolution in many organisms. However, only few studies have addressed the impact of STR variation on phenotypic variation at both the organismal and molecular levels. Potential driving forces underlying the high mutation rates of STRs also remain largely unknown. Here, we leverage the recently generated expression and STR variation data among wild Caenorhabditis elegans strains to conduct a genome-wide analysis of how STRs affect gene expression variation. We identify thousands of expression STRs (eSTRs) showing regulatory effects and demonstrate that they explain missing heritability beyond SNV-based expression quantitative trait loci. We illustrate specific regulatory mechanisms such as how eSTRs affect splicing sites and alternative splicing efficiency. We also show that differential expression of antioxidant genes and oxidative stresses might affect STR mutations systematically using both wild strains and mutation accumulation lines. Overall, we reveal the interplay between STRs and gene expression variation by providing novel insights into regulatory mechanisms of STRs and highlighting that oxidative stress could lead to higher STR mutation rates.

The motif composition of variable number tandem repeats impacts gene expression

Understanding the impact of DNA variation on human traits is a fundamental question in human genetics. Variable number tandem repeats (VNTRs) make up ∼3% of the human genome but are often excluded from association analysis owing to poor read mappability or divergent repeat content. Although methods exist to estimate VNTR length from short-read data, it is known that VNTRs vary in both length and repeat (motif) composition. Here, we use a repeat-pangenome graph (RPGG) constructed on 35 haplotype-resolved assemblies to detect variation in both VNTR length and repeat composition. We align population-scale data from the Genotype-Tissue Expression (GTEx) Consortium to examine how variations in sequence composition may be linked to expression, including cases independent of overall VNTR length. We find that 9422 out of 39,125 VNTRs are associated with nearby gene expression through motif variations, of which only 23.4% are accessible from length. Fine-mapping identifies 174 genes to be likely driven by variation in certain VNTR motifs and not overall length. We highlight two genes, CACNA1C and RNF213, that have expression associated with motif variation, showing the utility of RPGG analysis as a new approach for trait association in multiallelic and highly variable loci.

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.

Leveraging primate-specific genomic information for genetic studies of complex diseases

Genomic changes specific to higher primates are regarded as primate-specific genomic information (PSI). Using PSI to inform genetic studies is highly desirable but hampered by three factors: heterogeneity among PSI studies, lack of integrated profiles of the identified PSI elements and dearth of relevant functional information. We report a database of 19,767 PSI elements collated from nine types of brain-related studies, which form 19,473 non-overlapping PSI regions that distribute unevenly but jointly cover only 0.81% of the genome. About 2.5% of the PSI regions colocalized with variants identified in genome-wide association studies, with disease loci more likely colocalized than quantitative trait loci (p = 1.6 × 10−5), particularly in regions without obvious regulatory roles. We further showed an LRP4 exemplar region with PSI elements orchestrated with common and rare disease variants and other functional elements. Our results render PSI elements as a valuable source to inform genetic studies of complex diseases.

The genetic architecture of schizophrenia: review of large-scale genetic studies

Schizophrenia is a complex and often chronic psychiatric disorder with high heritability. Diagnosis of schizophrenia is still made clinically based on psychiatric symptoms; no diagnostic tests or biomarkers are available. Pathophysiology-based diagnostic scheme and treatments are also not available. Elucidation of the pathogenesis is needed for development of pathology-based diagnostics and treatments. In the past few decades, genetic research has made substantial advances in our understanding of the genetic architecture of schizophrenia. Rare copy number variations (CNVs) and rare single-nucleotide variants (SNVs) detected by whole-genome CNV analysis and whole-genome/-exome sequencing analysis have provided the great advances. Common single-nucleotide polymorphisms (SNPs) detected by large-scale genome-wide association studies have also provided important information. Large-scale genetic studies have been revealed that both rare and common genetic variants play crucial roles in this disorder. In this review, we focused on CNVs, SNVs, and SNPs, and discuss the latest research findings on the pathogenesis of schizophrenia based on these genetic variants. Rare variants with large effect sizes can provide mechanistic hypotheses. CRISPR-based genetics approaches and induced pluripotent stem cell technology can facilitate the functional analysis of these variants detected in patients with schizophrenia. Recent advances in long-read sequence technology are expected to detect variants that cannot be detected by short-read sequence technology. Various studies that bring together data from common variant and transcriptomic datasets provide biological insight. These new approaches will provide additional insight into the pathophysiology of schizophrenia and facilitate the development of pathology-based therapeutics.

Gene-Environment Interactions in Repeat Expansion Diseases: Mechanisms of Environmentally Induced Repeat Instability

Short tandem repeats (STRs) are units of 1-6 base pairs that occur in tandem repetition to form a repeat tract. STRs exhibit repeat instability, which generates expansions or contractions of the repeat tract. Over 50 diseases, primarily affecting the central nervous system and muscles, are characterized by repeat instability. Longer repeat tracts are typically associated with earlier age of onset and increased disease severity. Environmental exposures are suspected to play a role in the pathogenesis of repeat expansion diseases. Here, we review the current knowledge of mechanisms of environmentally induced repeat instability in repeat expansion diseases. The current evidence demonstrates that environmental factors modulate repeat instability via DNA damage and induction of DNA repair pathways, with distinct mechanisms for repeat expansion and contraction. Of particular note, oxidative stress is a key mediator of environmentally induced repeat instability. The preliminary evidence suggests epigenetic modifications as potential mediators of environmentally induced repeat instability. Future research incorporating an array of environmental exposures, new human cohorts, and improved model systems, with a continued focus on cell-types, tissues, and critical windows, will aid in identifying mechanisms of environmentally induced repeat instability. Identifying environmental modulators of repeat instability and their mechanisms of action will inform preventions, therapies, and public health measures.

Variable number tandem repeats (VNTRs) as modifiers of breast cancer risk in carriers of BRCA1 185delAG

Despite substantial efforts in identifying both rare and common variants affecting disease risk, in the majority of diseases, a large proportion of unexplained genetic risk remains. We propose that variable number tandem repeats (VNTRs) may explain a proportion of the missing genetic risk. Herein, in a pilot study with a retrospective cohort design, we tested whether VNTRs are causal modifiers of breast cancer risk in 347 female carriers of the BRCA1 185delAG pathogenic variant, an important group given their high risk of developing breast cancer. We performed targeted-capture to sequence VNTRs, called genotypes with adVNTR, tested the association of VNTRs and breast cancer risk using Cox regression models, and estimated the effect size using a retrospective likelihood approach. Of 303 VNTRs that passed quality control checks, 4 VNTRs were significantly associated with risk to develop breast cancer at false discovery rate [FDR] < 0.05 and an additional 4 VNTRs had FDR < 0.25. After determining the specific risk alleles, there was a significantly earlier age at diagnosis of breast cancer in carriers of the risk alleles compared to those without the risk alleles for seven of eight VNTRs. One example is a VNTR in exon 2 of LINC01973 with a per-allele hazard ratio of 1.58 (1.07-2.33) and 5.28 (2.79-9.99) for the homozygous risk-allele genotype. Results from this first systematic study of VNTRs demonstrate that VNTRs may explain a proportion of the unexplained genetic risk for breast cancer.

Influencing cognitive performance via social interactions: a novel therapeutic approach for brain disorders based on neuroanatomical mapping?

Many psychiatric and neurological disorders present deficits in both the social and cognitive domain. In this perspectives article, we provide an overview and the potential of the existence of an extensive neurobiological substrate underlying the close relationship between these two domains. By mapping the rodent brain regions involved in the social and/or cognitive domain, we show that the vast majority of brain regions involved in the cognitive domain are also involved in the social domain. The identified neuroanatomical overlap has an evolutionary basis, as complex social behavior requires cognitive skills, and aligns with the reported functional interactions of processes underlying cognitive and social performance. Based on the neuroanatomical mapping, recent (pre-)clinical findings, and the evolutionary perspective, we emphasize that the social domain requires more focus as an important treatment target and/or biomarker, especially considering the presently limited treatment strategies for these disorders.

Rare tandem repeat expansions associate with genes involved in synaptic and neuronal signaling functions in schizophrenia

Tandem repeat expansions (TREs) are associated with over 60 monogenic disorders and have recently been implicated in complex disorders such as cancer and autism spectrum disorder. The role of TREs in schizophrenia is now emerging. In this study, we have performed a genome-wide investigation of TREs in schizophrenia. Using genome sequence data from 1154 Swedish schizophrenia cases and 934 ancestry-matched population controls, we have detected genome-wide rare (<0.1% population frequency) TREs that have motifs with a length of 2-20 base pairs. We find that the proportion of individuals carrying rare TREs is significantly higher in the schizophrenia group. There is a significantly higher burden of rare TREs in schizophrenia cases than in controls in genic regions, particularly in postsynaptic genes, in genes overlapping brain expression quantitative trait loci, and in brain-expressed genes that are differentially expressed between schizophrenia cases and controls. We demonstrate that TRE-associated genes are more constrained and primarily impact synaptic and neuronal signaling functions. These results have been replicated in an independent Canadian sample that consisted of 252 schizophrenia cases of European ancestry and 222 ancestry-matched controls. Our results support the involvement of rare TREs in schizophrenia etiology.

Approaches to long-read sequencing in a clinical setting to improve diagnostic rate

Over the past decade, advances in genetic testing, particularly the advent of next-generation sequencing, have led to a paradigm shift in the diagnosis of molecular diseases and disorders. Despite our present collective ability to interrogate more than 90% of the human genome, portions of the genome have eluded us, resulting in stagnation of diagnostic yield with existing methodologies. Here we show how application of a new technology, long-read sequencing, has the potential to improve molecular diagnostic rates. Whole genome sequencing by long reads was able to cover 98% of next-generation sequencing dead zones, which are areas of the genome that are not interpretable by conventional industry-standard short-read sequencing. Through the ability of long-read sequencing to unambiguously call variants in these regions, we discovered an immunodeficiency due to a variant in IKBKG in a subject who had previously received a negative genome sequencing result. Additionally, we demonstrate the ability of long-read sequencing to detect small variants on par with short-read sequencing, its superior performance in identifying structural variants, and thirdly, its capacity to determine genomic methylation defects in native DNA. Though the latter technical abilities have been demonstrated, we demonstrate the clinical application of this technology to successfully identify multiple types of variants using a single test.

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.

CACNA1C (CaV1.2) and other L-type calcium channels in the pathophysiology and treatment of psychiatric disorders: Advances from functional genomics and pharmacoepidemiology

A role for voltage-gated calcium channels (VGCCs) in psychiatric disorders has long been postulated as part of a broader involvement of intracellular calcium signalling. However, the data were inconclusive and hard to interpret. We review three areas of research that have markedly advanced the field. First, there is now robust genomic evidence that common variants in VGCC subunit genes, notably CACNA1C which encodes the L-type calcium channel (LTCC) Ca_V1.2 subunit, are trans-diagnostically associated with psychiatric disorders including schizophrenia and bipolar disorder. Rare variants in these genes also contribute to the risk. Second, pharmacoepidemiological evidence supports the possibility that calcium channel blockers, which target LTCCs, might have beneficial effects on the onset or course of these disorders. This is especially true for calcium channel blockers that are brain penetrant. Third, long-range sequencing is revealing the repertoire of full-length LTCC transcript isoforms. Many novel and abundant CACNA1C isoforms have been identified in human and mouse brain, including some which are enriched compared to heart or aorta, and predicted to encode channels with differing functional and pharmacological properties. These isoforms may contribute to the molecular mechanisms of genetic association to psychiatric disorders. They may also enable development of therapeutic agents that can preferentially target brain LTCC isoforms and be of potential value for psychiatric indications.

Identifying enhancer properties associated with genetic risk for complex traits using regulome-wide association studies

Genetic risk for complex traits is strongly enriched in non-coding genomic regions involved in gene regulation, especially enhancers. However, we lack adequate tools to connect the characteristics of these disruptions to genetic risk. Here, we propose RWAS (Regulome Wide Association Study), a new application of the MAGMA software package to identify the characteristics of enhancers that contribute to genetic risk for disease. RWAS involves three steps: (i) assign genotyped SNPs to cell type- or tissue-specific regulatory features (e.g., enhancers); (ii) test associations of each regulatory feature with a trait of interest for which genome-wide association study (GWAS) summary statistics are available; (iii) perform enhancer-set enrichment analyses to identify quantitative or categorical features of regulatory elements that are associated with the trait. These steps are implemented as a novel application of MAGMA, a tool originally developed for gene-based GWAS analyses. Applying RWAS to interrogate genetic risk for schizophrenia, we discovered a class of risk-associated AT-rich enhancers that are active in the developing brain and harbor binding sites for multiple transcription factors with neurodevelopmental functions. RWAS utilizes open-source software, and we provide a comprehensive collection of annotations for tissue-specific enhancer locations and features, including their evolutionary conservation, AT content, and co-localization with binding sites for hundreds of TFs. RWAS will enable researchers to characterize properties of regulatory elements associated with any trait of interest for which GWAS summary statistics are available.

Antidepressants and Antipsychotics in Human Pregnancy: Transfer Across the Placenta and Opportunities for Modeling Studies

There is limited information about the transfer of antidepressants and antipsychotics across the human placenta. The objective of the current review was to systematically screen the scientific literature using relevant keywords to collect quantitative data on placental transfer of these drugs in humans and to give an overview of current modeling approaches used in this context. The collected data encompassed clinically measured fetal:maternal (F:M) concentration ratios (ie, the ratio between drug concentrations measured in the umbilical cord and drug concentrations measured in the mother) and transfer data obtained from ex vivo cotyledon perfusion experiments. These data were found for 18 antidepressants and some of their pharmacologically active metabolites, and for 10 antipsychotics and the metabolites thereof. Based on the collected data, similar maternal and fetal exposure could be observed for only a few compounds (eg, norfluoxetine and desvenlafaxine), whereas for most drugs (eg, paroxetine, sertraline, and quetiapine), fetal exposure appeared to be on average lower than maternal exposure. Venlafaxine appeared to be an exception in that the data indicated equivalent or higher concentrations in the umbilical cord than in the mother. Physiologically based pharmacokinetic (PBPK) models were sporadically used to investigate maternal pharmacokinetics of antidepressants or antipsychotics (eg, for sertraline, aripiprazole, and olanzapine), although without explicitly addressing fetal drug exposure. It is recommended that PBPK modeling is applied more frequently to these drugs. Although no substitute for clinical studies, these tools can help to better understand pregnancy-induced pharmacokinetic changes and ultimately contribute to a more evidence-based pharmacotherapy of depression and psychosis in pregnant subjects.

A phenome-wide association study identifies effects of copy-number variation of VNTRs and multicopy genes on multiple human traits

The human genome contains tens of thousands of large tandem repeats and hundreds of genes that show common and highly variable copy-number changes. Due to their large size and repetitive nature, these variable number tandem repeats (VNTRs) and multicopy genes are generally recalcitrant to standard genotyping approaches and, as a result, this class of variation is poorly characterized. However, several recent studies have demonstrated that copy-number variation of VNTRs can modify local gene expression, epigenetics, and human traits, indicating that many have a functional role. Here, using read depth from whole-genome sequencing to profile copy number, we report results of a phenome-wide association study (PheWAS) of VNTRs and multicopy genes in a discovery cohort of ∼35,000 samples, identifying 32 traits associated with copy number of 38 VNTRs and multicopy genes at 1% FDR. We replicated many of these signals in an independent cohort and observed that VNTRs showing trait associations were significantly enriched for expression QTLs with nearby genes, providing strong support for our results. Fine-mapping studies indicated that in the majority (∼90%) of cases, the VNTRs and multicopy genes we identified represent the causal variants underlying the observed associations. Furthermore, several lie in regions where prior SNV-based GWASs have failed to identify any significant associations with these traits. Our study indicates that copy number of VNTRs and multicopy genes contributes to diverse human traits and suggests that complex structural variants potentially explain some of the so-called "missing heritability" of SNV-based GWASs.

DeepRepeat: direct quantification of short tandem repeats on signal data from nanopore sequencing

Despite recent improvements in basecalling accuracy, nanopore sequencing still has higher error rates on short-tandem repeats (STRs). Instead of using basecalled reads, we developed DeepRepeat which converts ionic current signals into red-green-blue channels, thus transforming the repeat detection problem into an image recognition problem. DeepRepeat identifies and accurately quantifies telomeric repeats in the CHM13 cell line and achieves higher accuracy in quantifying repeats in long STRs than competing methods. We also evaluate DeepRepeat on genome-wide or candidate region datasets from seven different sources. In summary, DeepRepeat enables accurate quantification of long STRs and complements existing methods relying on basecalled reads.

Advancing genomic technologies and clinical awareness accelerates discovery of disease-associated tandem repeat sequences

Expansions of gene-specific DNA tandem repeats (TRs), first described in 1991 as a disease-causing mutation in humans, are now known to cause >60 phenotypes, not just disease, and not only in humans. TRs are a common form of genetic variation with biological consequences, observed, so far, in humans, dogs, plants, oysters, and yeast. Repeat diseases show atypical clinical features, genetic anticipation, and multiple and partially penetrant phenotypes among family members. Discovery of disease-causing repeat expansion loci accelerated through technological advances in DNA sequencing and computational analyses. Between 2019 and 2021, 17 new disease-causing TR expansions were reported, totaling 63 TR loci (>69 diseases), with a likelihood of more discoveries, and in more organisms. Recent and historical lessons reveal that properly assessed clinical presentations, coupled with genetic and biological awareness, can guide discovery of disease-causing unstable TRs. We highlight critical but underrecognized aspects of TR mutations. Repeat motifs may not be present in current reference genomes but will be in forthcoming gapless long-read references. Repeat motif size can be a single nucleotide to kilobases/unit. At a given locus, repeat motif sequence purity can vary with consequence. Pathogenic repeats can be "insertions" within nonpathogenic TRs. Expansions, contractions, and somatic length variations of TRs can have clinical/biological consequences. TR instabilities occur in humans and other organisms. TRs can be epigenetically modified and/or chromosomal fragile sites. We discuss the expanding field of disease-associated TR instabilities, highlighting prospects, clinical and genetic clues, tools, and challenges for further discoveries of disease-causing TR instabilities and understanding their biological and pathological impacts-a vista that is about to expand.

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

Genome-wide association studies (GWASs) have revealed substantial genetic components comprised of single nucleotide polymorphisms (SNPs) in the heritable risk of psychiatric disorders. However, genetic risk factors not covered by GWAS also play pivotal roles in these illnesses. Tandem repeats, which are likely functional but frequently overlooked by GWAS, may account for an important proportion in the "missing heritability" of psychiatric disorders. Despite difficulties in characterizing and quantifying tandem repeats in the genome, studies have been carried out in an attempt to describe impact of tandem repeats on gene regulation and human phenotypes. In this review, we have introduced recent research progress regarding the genomic distribution and regulatory mechanisms of tandem repeats. We have also summarized the current knowledge of the genetic architecture and biological underpinnings of psychiatric disorders brought by studies of tandem repeats. These findings suggest that tandem repeats, in candidate psychiatric risk genes or in different levels of linkage disequilibrium (LD) with psychiatric GWAS SNPs and haplotypes, may modulate biological phenotypes related to psychiatric disorders (e.g., cognitive function and brain physiology) through regulating alternative splicing, promoter activity, enhancer activity and so on. In addition, many tandem repeats undergo tight natural selection in the human lineage, and likely exert crucial roles in human brain evolution. Taken together, the putative roles of tandem repeats in the pathogenesis of psychiatric disorders is strongly implicated, and using examples from previous literatures, we wish to call for further attention to tandem repeats in the post-GWAS era of psychiatric disorders.

Calcium imaging reveals depressive- and manic-phase-specific brain neural activity patterns in a murine model of bipolar disorder: a pilot study

Brain pathological features during manic/hypomanic and depressive episodes in the same patients with bipolar disorder (BPD) have not been described precisely. The study aimed to investigate depressive and manic-phase-specific brain neural activity patterns of BPD in the same murine model to provide information guiding investigation of the mechanism of phase switching and tailored prevention and treatment for patients with BPD. In vivo two-photon imaging was used to observe brain activity alterations in the depressive and manic phases in the same murine model of BPD. Two-photon imaging showed significantly reduced Ca²⁺ activity in temporal cortex pyramidal neurons in the depression phase in mice exposed to chronic unpredictable mild stress (CUMS), but not in the manic phase in mice exposed to CUMS and ketamine. Total integrated calcium values correlated significantly with immobility times. Brain Ca²⁺ hypoactivity was observed in the depression and manic phases in the same mice exposed to CUMS and ketamine relative to naïve controls. The novel object recognition preference ratio correlated negatively with the immobility time in the depression phase and the total distance traveled in the manic phase. With recognition of its limitations, this study revealed brain neural activity impairment indicating that intrinsic emotional network disturbance is a mechanism of BPD and that brain neural activity is associated with cognitive impairment in the depressive and manic phases of this disorder. These findings are consistent with those from macro-imaging studies of patients with BPD. The observed correlation of brain neural activity with the severity of depressive, but not manic, symptoms need to be investigated further.

Missing heritability may be hiding in repeats

[Figure: see text].

Characterizing nucleotide variation and expansion dynamics in human-specific variable number tandem repeats

There are more than 55,000 variable number tandem repeats (VNTRs) in the human genome, notable for both their striking polymorphism and mutability. Despite their role in human evolution and genomic variation, they have yet to be studied collectively and in detail, partially owing to their large size, variability, and predominant location in noncoding regions. Here, we examine 467 VNTRs that are human-specific expansions, unique to one location in the genome, and not associated with retrotransposons. We leverage publicly available long-read genomes, including from the Human Genome Structural Variant Consortium, to ascertain the exact nucleotide composition of these VNTRs and compare their composition of alleles. We then confirm repeat unit composition in more than 3000 short-read samples from the 1000 Genomes Project. Our analysis reveals that these VNTRs contain highly structured repeat motif organization, modified by frequent deletion and duplication events. Although overall VNTR compositions tend to remain similar between 1000 Genomes Project superpopulations, we describe a notable exception with substantial differences in repeat composition (in PCBP3), as well as several VNTRs that are significantly different in length between superpopulations (in ART1, PROP1, DYNC2I1, and LOC102723906). We also observe that most of these VNTRs are expanded in archaic human genomes, yet remain stable in length between single generations. Collectively, our findings indicate that repeat motif variability, repeat composition, and repeat length are all informative modalities to consider when characterizing VNTRs and their contribution to genomic variation.

Profiling variable-number tandem repeat variation across populations using repeat-pangenome graphs

Variable number tandem repeats (VNTRs) are composed of consecutive repetitive DNA with hypervariable repeat count and composition. They include protein coding sequences and associations with clinical disorders. It has been difficult to incorporate VNTR analysis in disease studies that use short-read sequencing because the traditional approach of mapping to the human reference is less effective for repetitive and divergent sequences. In this work, we solve VNTR mapping for short reads with a repeat-pangenome graph (RPGG), a data structure that encodes both the population diversity and repeat structure of VNTR loci from multiple haplotype-resolved assemblies. We develop software to build a RPGG, and use the RPGG to estimate VNTR composition with short reads. We use this to discover VNTRs with length stratified by continental population, and expression quantitative trait loci, indicating that RPGG analysis of VNTRs will be critical for future studies of diversity and disease.