Structural insights into i-motif DNA structures in sequences from the insulin-linked polymorphic region.

The insulin-linked polymorphic region is a variable number of tandem repeats region of DNA in the promoter of the insulin gene that regulates transcription of insulin. This region is known to form the alternative DNA structures, i-motifs and G-quadruplexes. Individuals have different sequence variants of tandem repeats and although previous work investigated the effects of some variants on G-quadruplex formation, there is not a clear picture of the relationship between the sequence diversity, the DNA structures formed, and the functional effects on insulin gene expression. Here we show that different sequence variants of the insulin linked polymorphic region form different DNA structures in vitro. Additionally, reporter genes in cellulo indicate that insulin expression may change depending on which DNA structures form. We report the crystal structure and dynamics of an intramolecular i-motif, which reveal sequences within the loop regions forming additional stabilising interactions that are critical to formation of stable i-motif structures. The outcomes of this work reveal the detail in formation of stable i-motif DNA structures, with potential for rational based drug design for compounds to target i-motif DNA.

The pan-tandem repeat map highlights multiallelic variants underlying gene expression and agronomic traits in rice.

Tandem repeats (TRs) are genomic regions that tandemly change in repeat number, which are often multiallelic. Their characteristics and contributions to gene expression and quantitative traits in rice are largely unknown. Here, we survey rice TR variations based on 231 genome assemblies and the rice pan-genome graph. We identify 227,391 multiallelic TR loci, including 54,416 TR variations that are absent from the Nipponbare reference genome. Only 1/3 TR variations show strong linkage with nearby bi-allelic variants (SNPs, Indels and PAVs). Using 193 panicle and 202 leaf transcriptomic data, we reveal 485 and 511 TRs act as QTLs independently of other bi-allelic variations to nearby gene expression, respectively. Using plant height and grain width as examples, we identify and validate TRs contributions to rice agronomic trait variations. These findings would enhance our understanding of the functions of multiallelic variants and facilitate rice molecular breeding.

ModDotPlot-rapid and interactive visualization of tandem repeats.

MOTIVATION: A common method for analyzing genomic repeats is to produce a sequence similarity matrix visualized via a dot plot. Innovative approaches such as StainedGlass have improved upon this classic visualization by rendering dot plots as a heatmap of sequence identity, enabling researchers to better visualize multi-megabase tandem repeat arrays within centromeres and other heterochromatic regions of the genome. However, computing the similarity estimates for heatmaps requires high computational overhead and can suffer from decreasing accuracy. RESULTS: In this work, we introduce ModDotPlot, an interactive and alignment-free dot plot viewer. By approximating average nucleotide identity via a k-mer-based containment index, ModDotPlot produces accurate plots orders of magnitude faster than StainedGlass. We accomplish this through the use of a hierarchical modimizer scheme that can visualize the full 128 Mb genome of Arabidopsis thaliana in under 5 min on a laptop. ModDotPlot is bundled with a graphical user interface supporting real-time interactive navigation of entire chromosomes. AVAILABILITY AND IMPLEMENTATION: ModDotPlot is available at https://github.com/marbl/ModDotPlot.

DTDHM: detection of tandem duplications based on hybrid methods using next-generation sequencing data.

Background: Tandem duplication (TD) is a common and important type of structural variation in the human genome. TDs have been shown to play an essential role in many diseases, including cancer. However, it is difficult to accurately detect TDs due to the uneven distribution of reads and the inherent complexity of next-generation sequencing (NGS) data. Methods: This article proposes a method called DTDHM (detection of tandem duplications based on hybrid methods), which utilizes NGS data to detect TDs in a single sample. DTDHM builds a pipeline that integrates read depth (RD), split read (SR), and paired-end mapping (PEM) signals. To solve the problem of uneven distribution of normal and abnormal samples, DTDHM uses the K-nearest neighbor (KNN) algorithm for multi-feature classification prediction. Then, the qualified split reads and discordant reads are extracted and analyzed to achieve accurate localization of variation sites. This article compares DTDHM with three other methods on 450 simulated datasets and five real datasets. Results: In 450 simulated data samples, DTDHM consistently maintained the highest F1-score. The average F1-score of DTDHM, SVIM, TARDIS, and TIDDIT were 80.0%, 56.2%, 43.4%, and 67.1%, respectively. The F1-score of DTDHM had a small variation range and its detection effect was the most stable and 1.2 times that of the suboptimal method. Most of the boundary biases of DTDHM fluctuated around 20 bp, and its boundary deviation detection ability was better than TARDIS and TIDDIT. In real data experiments, five real sequencing samples (NA19238, NA19239, NA19240, HG00266, and NA12891) were used to test DTDHM. The results showed that DTDHM had the highest overlap density score (ODS) and F1-score of the four methods. Conclusions: Compared with the other three methods, DTDHM achieved excellent results in terms of sensitivity, precision, F1-score, and boundary bias. These results indicate that DTDHM can be used as a reliable tool for detecting TDs from NGS data, especially in the case of low coverage depth and tumor purity samples.

Tandem repeats in the genome of Toxoplasma gondii display compositional bias that impacts in protein structure.

Repetitive elements in DNA sequences are a hallmark of Apicomplexan protozoa. A genome-wide screening for Tandem Repeats was conducted in Toxoplasma gondii and related Coccidian parasites with a novel strategy to assess compositional bias. A conserved pattern of GC skew and purine-pyrimidine bias was observed. Compositional bias was also present at the protein level. Glutamic acid was the most abundant amino acid in the purine (GA) rich cluster, while Serine prevailed in pyrimidine (CT) rich cluster. Purine rich repeats, and consequently glutamic acid abundance, correlated with high scores for intrinsically disordered protein regions/domains. Finally, variability was established for repetitive regions within a well-known rhoptry antigen (ROP1) and an uncharacterized hypothetical protein with similar features. The approach we present could be useful to identify potential antigens bearing repetitive elements.

Pindel-TD: A Tandem Duplication Detector Based on A Pattern Growth Approach.

Tandem duplication (TD) is a major type of structural variations (SVs) that plays an important role in novel gene formation and human diseases. However, TDs are often missed or incorrectly classified as insertions by most modern SV detection methods due to the lack of specialized operation on TD-related mutational signals. Herein, we developed a TD detection module for the Pindel tool, referred to as Pindel-TD, based on a TD-specific pattern growth approach. Pindel-TD is capable of detecting TDs with a wide size range at single nucleotide resolution. Using simulated and real read data from HG002, we demonstrated that Pindel-TD outperforms other leading methods in terms of precision, recall, F1-score, and robustness. Furthermore, by applying Pindel-TD to data generated from the K562 cancer cell line, we identified a TD located at the seventh exon of SAGE1, providing an explanation for its high expression. Pindel-TD is available for non-commercial use at https://github.com/xjtu-omics/pindel.

Immune function of a C-type lectin with long tandem repeats and abundant threonine in the ridgetail white prawn Exopalaemon carinicauda.

C-type lectins (CTLs) are an important class of pattern recognition receptors (PRRs) that exhibit structural and functional diversity in invertebrates. Repetitive DNA sequences are ubiquitous in eukaryotic genomes, representing distinct modes of genome evolution and promoting new gene generation. Our study revealed a new CTL that is composed of two long tandem repeats, abundant threonine, and one carbohydrate recognition domain (CRD) in Exopalaemon carinicauda and has been designated EcTR-CTL. The full-length cDNA of EcTR-CTL was 1242 bp long and had an open reading frame (ORF) of 999 bp that encoded a protein of 332 amino acids. The genome structure of EcTR-CTL contains 4 exons and 3 introns. The length of each repeat unit in EcTR-CTL was 198 bp, which is different from the short tandem repeats reported previously in prawns and crayfish. EcTR-CTL was abundantly expressed in the intestine and hemocytes. After Vibrio parahaemolyticus and white spot syndrome virus (WSSV) challenge, the expression level of EcTR-CTL in the intestine was upregulated. Knockdown of EcTR-CTL downregulated the expression of anti-lipopolysaccharide factor, crustin, and lysozyme during Vibrio infection. The recombinant CRD of EcTR-CTL (rCRD) could bind to bacteria, lipopolysaccharides, and peptidoglycans. Additionally, rCRD can directly bind to WSSV. These findings indicate that 1) CTLs with tandem repeats may be ubiquitous in crustaceans, 2) EcTR-CTL may act as a PRR to participate in the innate immune defense against bacteria via nonself-recognition and antimicrobial peptide regulation, and 3) EcTR-CTL may play a positive or negative role in the process of WSSV infection by capturing virions.

Species-specific variation in mitochondrial genome tandem repeat polymorphisms in hares (Lepus spp., Lagomorpha, Leporidae) provides insight into their evolution.

The non-coding regions of the mitochondrial DNAs (mtDNAs) of hares, rabbits, and pikas (Lagomorpha) contain short (∼20 bp) and long (130-160 bp) tandem repeats, absent in related mammalian orders. In the presented study, we provide in-depth analysis for mountain hare (Lepus timidus) and brown hare (L. europaeus) mtDNA non-coding regions, together with a species- and population-level analysis of tandem repeat variation. Mountain hare short tandem repeats (SRs) as well as other analyzed hare species consist of two conserved 10 bp motifs, with only brown hares exhibiting a single, more variable motif. Long tandem repeats (LRs) also differ in sequence and copy number between species. Mountain hares have four to seven LRs, median value five, while brown hares exhibit five to nine LRs, median value six. Interestingly, introgressed mountain hare mtDNA in brown hares obtained an intermediate LR length distribution, with median copy number being the same as with conspecific brown hare mtDNA. In contrast, transfer of brown hare mtDNA into cultured mtDNA-less mountain hare cells maintained the original LR number, whereas the reciprocal transfer caused copy number instability, suggesting that cellular environment rather than the nuclear genomic background plays a role in the LR maintenance. Due to their dynamic nature and separation from other known conserved sequence elements on the non-coding region of hare mitochondrial genomes, the tandem repeat elements likely to represent signatures of ancient genetic rearrangements. clarifying the nature and dynamics of these rearrangements may shed light on the possible role of NCR repeated elements in mitochondria and in species evolution.

Characteristics of tandem repeat inheritance and sympathetic nerve involvement in GAA-FGF14 ataxia.

BACKGROUND: Intronic GAA repeat expansion ([GAA] ≥250) in FGF14 is associated with the late-onset neurodegenerative disorder, spinocerebellar ataxia 27B (SCA27B, GAA-FGF14 ataxia). We aim to determine the prevalence of the GAA repeat expansion in FGF14 in Chinese populations presenting late-onset cerebellar ataxia (LOCA) and evaluate the characteristics of tandem repeat inheritance, radiological features and sympathetic nerve involvement. METHODS: GAA-FGF14 repeat expansion was screened in an undiagnosed LOCA cohort (n = 664) and variations in repeat-length were analyzed in families of confirmed GAA-FGF14 ataxia patients. Brain magnetic resonance imaging (MRI) was used to evaluate the radiological feature in GAA-FGF14 ataxia patients. Clinical examinations and sympathetic skin response (SSR) recordings in GAA-FGF14 patients (n = 16) were used to quantify sympathetic nerve involvement. RESULTS: Two unrelated probands (2/664) were identified. Genetic screening for GAA-FGF14 repeat expansion was performed in 39 family members, 16 of whom were genetically diagnosed with GAA-FGF14 ataxia. Familial screening revealed expansion of GAA repeats in maternal transmissions, but contraction upon paternal transmission. Brain MRI showed slight to moderate cerebellar atrophy. SSR amplitude was lower in GAA-FGF14 patients in pre-symptomatic stage compared to healthy controls, and further decreased in the symptomatic stage. CONCLUSIONS: GAA-FGF14 ataxia was rare among Chinese LOCA cases. Parental gender appears to affect variability in GAA repeat number between generations. Reduced SSR amplitude is a prominent feature in GAA-FGF14 patients, even in the pre-symptomatic stage.

Paramutation at the maize pl1 locus is associated with RdDM activity at distal tandem repeats.

Exceptions to Mendelian inheritance often highlight novel chromosomal behaviors. The maize Pl1-Rhoades allele conferring plant pigmentation can display inheritance patterns deviating from Mendelian expectations in a behavior known as paramutation. However, the chromosome features mediating such exceptions remain unknown. Here we show that small RNA production reflecting RNA polymerase IV function within a distal downstream set of five tandem repeats is coincident with meiotically-heritable repression of the Pl1-Rhoades transcription unit. A related pl1 haplotype with three, but not one with two, repeat units also displays the trans-homolog silencing typifying paramutations. 4C interactions, CHD3a-dependent small RNA profiles, nuclease sensitivity, and polyadenylated RNA levels highlight a repeat subregion having regulatory potential. Our comparative and mutant analyses show that transcriptional repression of Pl1-Rhoades correlates with 24-nucleotide RNA production and cytosine methylation at this subregion indicating the action of a specific DNA-dependent RNA polymerase complex. These findings support a working model in which pl1 paramutation depends on trans-chromosomal RNA-directed DNA methylation operating at a discrete cis-linked and copy-number-dependent transcriptional regulatory element.

Genome comparison of long-circulating field CnmeGV isolates from the same region.

Cnaphalocrocis medinalis granulovirus (CnmeGV), belonging to Betabaculovirus cnamedinalis, can infect the rice pest, the rice leaf roller. In 1979, a CnmeGV isolate, CnmeGV-EP, was collected from Enping County, China. In 2014, we collected another CnmeGV isolate, CnmeGV-EPDH3, at the same location and obtained the complete virus genome sequence using Illumina and ONT sequencing technologies. By combining these two virus isolates, we updated the genome annotation of CnmeGV and conducted an in-depth analysis of its genome features. CnmeGV genome contains abundant tandem repeat sequences, and the repeating units in the homologous regions (hrs) exhibit overlapping and nested patterns. The genetic variations within EPDH3 population show the high stability of CnmeGV genome, and tandem repeats are the only region of high genetic variation in CnmeGV genome replication. Some defective viral genomes formed by recombination were found within the population. Comparison analysis of the two virus isolates collected from Enping showed that the proteins encoded by the CnmeGV-specific genes were less conserved relative to the baculovirus core genes. At the genomic level, there are a large number of SNPs and InDels between the two virus isolates, especially in and around the bro genes and hrs. Additionally, we discovered that CnmeGV acquired a segment of non-ORF sequence from its host, which does not provide any new proteins but rather serves as redundant genetic material integrated into the viral genome. Furthermore, we observed that the host's transposon piggyBac has inserted into some virus genes. Together, dsDNA viruses could acquire non-coding genetic material from their hosts to expand the size of their genomes. These findings provide new insights into the evolution of dsDNA viruses.

A new tandem repeat-enriched lncRNA XLOC_008672 promotes gastric carcinogenesis by regulating G3BP1 expression.

Aberrant expression of forkhead box transcription factor 1 (FOXM1) plays critical roles in a variety of human malignancies and predicts poor prognosis. However, little is known about the crosstalk between FOXM1 and long noncoding RNAs (lncRNAs) in tumorigenesis. The present study identifies a previously uncharacterized lncRNA XLOC_008672 in gastric cancer (GC), which is regulated by FOXM1 and possesses multiple copies of tandem repetitive sequences. LncRNA microarrays are used to screen differentially expressed lncRNAs in FOXM1 knockdown GC cells, and then the highest fold downregulation lncRNA XLOC_008672 is screened out. Sequence analysis reveals that the new lncRNA contains 62 copies of 37-bp tandem repeats. It is transcriptionally activated by FOXM1 and functions as a downstream effector of FOXM1 in GC cells through in vitro and in vivo functional assays. Elevated expression of XLOC_008672 is found in GC tissues and indicates worse prognosis. Mechanistically, XLOC_008672 can bind to small nuclear ribonucleoprotein polypeptide A (SNRPA), thereby enhancing mRNA stability of Ras-GTPase-activating protein SH3 domain-binding protein 1 (G3BP1) and, consequently, facilitating GC cell proliferation and migration. Our study discovers a new uncharacterized lncRNA XLOC_008672 involved in GC carcinogenesis and progression. Targeting FOXM1/XLOC_008672/SNRPA/G3BP1 signaling axis might be a promising therapeutic strategy for GC.

A genome-wide spectrum of tandem repeat expansions in 338,963 humans.

The Genome Aggregation Database (gnomAD), widely recognized as the gold-standard reference map of human genetic variation, has largely overlooked tandem repeat (TR) expansions, despite the fact that TRs constitute ∼6% of our genome and are linked to over 50 human diseases. Here, we introduce the TR-gnomAD (https://wlcb.oit.uci.edu/TRgnomAD), a biobank-scale reference of 0.86 million TRs derived from 338,963 whole-genome sequencing (WGS) samples of diverse ancestries (39.5% non-European samples). TR-gnomAD offers critical insights into ancestry-specific disease prevalence using disparities in TR unit number frequencies among ancestries. Moreover, TR-gnomAD is able to differentiate between common, presumably benign TR expansions, which are prevalent in TR-gnomAD, from those potentially pathogenic TR expansions, which are found more frequently in disease groups than within TR-gnomAD. Together, TR-gnomAD is an invaluable resource for researchers and physicians to interpret TR expansions in individuals with genetic diseases.

UBTF tandem duplications in pediatric myelodysplastic syndrome and acute myeloid leukemia: implications for clinical screening and diagnosis.

Recent genomic studies in adult and pediatric acute myeloid leukemia (AML) demonstrated recurrent in-frame tandem duplications (TD) in exon 13 of upstream binding transcription factor (UBTF). These alterations, which account for approximately 4.3% of AML in childhood and about 3% in adult AML aged <60 years of age, are subtype-defining and associated with poor outcomes. Here, we provide a comprehensive investigation into the clinicopathological features of UBTF-TD myeloid neoplasms in childhood, including 89 unique pediatric AML and 6 myelodysplastic syndrome (MDS) cases harboring a tandem duplication in exon 13 of UBTF. We demonstrate that UBTF-TD myeloid tumors are associated with dysplastic features, low bone marrow blast infiltration, and low white blood cell count. Furthermore, using bulk and single-cell analyses, we confirm that UBTF-TD is an early and clonal event associated with a distinct transcriptional profile, whereas the acquisition of FLT3 or WT1 mutations is associated with more stem cell-like programs. Lastly, we report rare duplications within exon 9 of UBTF that phenocopy exon 13 duplications, expanding the spectrum of UBTF alterations in pediatric myeloid tumors. Collectively, we comprehensively characterize pediatric AML and MDS with UBTF-TD, and highlight key clinical and pathologic features that distinguish this new entity from other molecular subtypes of AML.

Genome-wide enhancer-associated tandem repeats are expanded in cardiomyopathy.

BACKGROUND: Cardiomyopathy is a clinically and genetically heterogeneous heart condition that can lead to heart failure and sudden cardiac death in childhood. While it has a strong genetic basis, the genetic aetiology for over 50% of cardiomyopathy cases remains unknown. METHODS: In this study, we analyse the characteristics of tandem repeats from genome sequence data of unrelated individuals diagnosed with cardiomyopathy from Canada and the United Kingdom (n = 1216) and compare them to those found in the general population. We perform burden analysis to identify genomic and epigenomic features that are impacted by rare tandem repeat expansions (TREs), and enrichment analysis to identify functional pathways that are involved in the TRE-associated genes in cardiomyopathy. We use Oxford Nanopore targeted long-read sequencing to validate repeat size and methylation status of one of the most recurrent TREs. We also compare the TRE-associated genes to those that are dysregulated in the heart tissues of individuals with cardiomyopathy. FINDINGS: We demonstrate that tandem repeats that are rarely expanded in the general population are predominantly expanded in cardiomyopathy. We find that rare TREs are disproportionately present in constrained genes near transcriptional start sites, have high GC content, and frequently overlap active enhancer H3K27ac marks, where expansion-related DNA methylation may reduce gene expression. We demonstrate the gene silencing effect of expanded CGG tandem repeats in DIP2B through promoter hypermethylation. We show that the enhancer-associated loci are found in genes that are highly expressed in human cardiomyocytes and are differentially expressed in the left ventricle of the heart in individuals with cardiomyopathy. INTERPRETATION: Our findings highlight the underrecognized contribution of rare tandem repeat expansions to the risk of cardiomyopathy and suggest that rare TREs contribute to ∼4% of cardiomyopathy risk. FUNDING: Government of Ontario (RKCY), The Canadian Institutes of Health Research PJT 175329 (RKCY), The Azrieli Foundation (RKCY), SickKids Catalyst Scholar in Genetics (RKCY), The University of Toronto McLaughlin Centre (RKCY, SM), Ted Rogers Centre for Heart Research (SM), Data Sciences Institute at the University of Toronto (SM), The Canadian Institutes of Health Research PJT 175034 (SM), The Canadian Institutes of Health Research ENP 161429 under the frame of ERA PerMed (SM, RL), Heart and Stroke Foundation of Ontario & Robert M Freedom Chair in Cardiovascular Science (SM), Bitove Family Professorship of Adult Congenital Heart Disease (EO), Canada Foundation for Innovation (SWS, JR), Canada Research Chair (PS), Genome Canada (PS, JR), The Canadian Institutes of Health Research (PS).

Characterization of the diversity of barn owl's mitochondrial genome reveals high copy number variations in the control region.

Mitochondria are known to play an essential role in the cell. These organelles contain their own DNA, which is divided in a coding and non-coding region (NCR). While much of the NCR's function is unknown, tandem repeats have been observed in several vertebrates, with extreme intra-individual, intraspecific and interspecific variation. Taking advantage of a new complete reference for the mitochondrial genome of the Afro-European Barn Owl (Tyto alba), as well as 172 whole genome-resequencing; we (i) describe the reference mitochondrial genome with a special focus on the repeats in the NCR, (ii) quantify the variation in number of copies between individuals, and (iii) explore the possible factors associated with the variation in the number of repetitions. The reference mitochondrial genome revealed a long (256bp) and a short (80bp) tandem repeat in the NCR region. The re-sequenced genomes showed a great variation in number of copies between individuals, with 4 to 38 copies of the Long and 6 to 135 copies of the short repeat. Among the factors associated with this variation between individuals, the tissue used for extraction was the most significant. The exact mechanisms of the formations of these repeats are still to be discovered and understanding them will help explain the maintenance of the polymorphism in the number of copies, as well as their interactions with the metabolism, the aging and health of the individuals.

Structural peculiarities of tandem repeats and their clinical significance.

While it is well established that a mere 2% of human DNA nucleotides are involved in protein coding, the remainder of the DNA plays a vital role in the preservation of normal cellular genetic function. A significant proportion of tandem repeats (TRs) are present in non-coding DNA. TRs - specific sequences of nucleotides that entail numerous repetitions of a given fragment. In this study, we employed our novel algorithm grounded in finite automata theory, which we refer to as Dafna, to investigate for the first time the likelihood of these nucleotide sequences forming non-canonical DNA structures (NS). Such structures include G-quadruplexes, i-motifs, hairpins, and triplexes. The tandem repeats under consideration in our research encompassed sequences containing 1 to 6 nucleotides per repeated fragment. For comparison, we employed a set of randomly generated sequences of the same length (60 nucleotides) as a benchmark. The outcomes of our research exposed a disparity between the potential for NS formation in random sequences and tandem repeats. Our findings affirm that the propensity of DNA and RNA to form NS is closely tied to various genetic disorders, including Huntington's disease, Fragile X syndrome, and Friedreich's ataxia. In the concluding discussion, we present a proposal for a new therapeutic mechanism to address these diseases. This novel approach revolves around the ability of specific nucleic acid fragments to form multiple types of NS.