Comparative genomics and molecular dynamics of DNA repeats in eukaryotes

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.

Chromosome mapping of retrotransposon AviRTE in a neotropical bird species: Trogon surrucura (Trogoniformes; Trogonidae)

Avian genomes are characterized as being more compact than other amniotes, with less diversity and density of transposable elements (TEs). In addition, birds usually show bimodal karyotypes, exhibiting a great variation in diploid numbers. Some species present unusually large sex chromosomes, possibly due to the accumulation of repetitive sequences. Avian retrotransposon-like element (AviRTE) is a long interspersed nuclear element (LINE) recently discovered in the genomes of birds and nematodes, and it is still poorly characterized in terms of chromosomal mapping and phylogenetic relationships. In this study, we mapped AviRTE isolated from the Trogon surrucura genome into the T. surrucura (TSU) karyotype. Furthermore, we analyzed the phylogenetic relationships of this LINE in birds and other vertebrates. Our results showed that the distribution pattern of AviRTE is not restricted to heterochromatic regions, with accumulation on the W chromosome of TSU, yet another species with an atypical sex chromosome and TE hybridization. The phylogenetic analysis of AviRTE sequences in birds agreed with the proposed phylogeny of species in most clades, and allowed the detection of this sequence in other species, expanding the distribution of the element.

A pilot study on the detection of microsatellite instability using long mononucleotide repeats in solid tumors

Microsatellite instability (MSI) status is a prognostic biomarker for immunotherapy in certain types of cancers, such as colorectal cancers (CRCs) and endometrial cancers (ECs). Tumors that are categorized as having high MSI (MSI-H) express high levels of neoantigens for immune recognition. The typical MSI test measures the length of short mononucleotide repeats (SMR) poly(A) 21-27; however, a limitation of this test is the difficulty in determining the shift size, particularly in endometrial cancer. To investigate an MSI detection assay with improved performance, the present study analyzed the use of poly(A) 40-44 mononucleotide repeats to detect the MSI status of 100 patients with either CRC (n=50) or EC (n=50). Capillary electrophoresis was used to evaluate five long mononucleotide repeat (LMR) markers, including poly(A) 40-A, 40-B, 40-C, 40-D and 44. The concordance rate of the LMR-MSI assay compared with an immunohistochemistry MSI detection assay was 96.0 and 95.1% for CRCs and ECs respectively, with the detection limit of the LMR-MSI assay demonstrated to be 2.5% MSI-H in HCT116 colorectal carcinoma cell lines. The LMR-MSI assay yielded a 95.1% concordance rate in ECs compared with that in the SMR-MSI test (87.8%). The LMR-MSI test identified a significantly higher mean shift size (13 bp) in MSI-H tumors compared with the SMR-MSI test (10 bp), in both EC and CRC tissue samples. Together, the present study suggested that the LMR-MSI test could potentially be a sensitive and practical technology for molecular laboratory testing, particularly in the use of immunotherapy for patients with CRCs and ECs.

Structure and Evolution of Ribosomal Genes of Insect Chromosomes

Currently, clusters of 45S and 5S ribosomal DNA (rDNA) have been studied in about 1000 and 100 species of the class Insecta, respectively. Although the number of insect species with known 45S rDNA clusters (also referred to as nucleolus-organizing regions, or NORs) constitutes less than 0.1 percent of the described members of this enormous group, certain conclusions can already be drawn. Since haploid karyotypes with single 45S and 5S rDNA clusters predominate in both basal and derived insect groups, this character state is apparently ancestral for the class Insecta in general. Nevertheless, the number, chromosomal location, and other characteristics of both 45S and 5S rDNA sites substantially vary across different species, and sometimes even within the same species. There are several main factors and molecular mechanisms that either maintain these parameters or alter them on the short-term and/or long-term scale. Chromosome structure (i.e., monocentric vs. holokinetic chromosomes), excessive numbers of rRNA gene copies per cluster, interactions with transposable elements, pseudogenization, and meiotic recombination are perhaps the most important among them.

Repetitive element expansions contribute to genome size gigantism in Pamphagidae: A comparative study (Orthoptera, Acridoidea)

Pamphagidae is a family of Acridoidea that inhabits the desert steppes of Eurasia and Africa. This study employed flow cytometry to estimate the genome size of eight species in the Pamphagidae. The results indicate that the genome size of the eight species ranged from 13.88 pg to 14.66 pg, with an average of 14.26 pg. This is the largest average genome size recorded for the Orthoptera families, as well as for the entire Insecta. Furthermore, the study explored the role of repetitive sequences in the genome, including their evolutionary dynamics and activity, using low-coverage next-generation sequencing data. The genome is composed of 14 different types of repetitive sequences, which collectively make up between 59.9% and 68.17% of the total genome. The Pamphagidae family displays high levels of transposable element (TE) activity, with the number of TEs increasing and accumulating since the family's emergence. The study found that the types of repetitive sequences contributing to the TE outburst events are similar across species. Additionally, the study identified unique repetitive elements for each species. The differences in repetitive sequences among the eight Pamphagidae species correspond to their phylogenetic relationships. The study sheds new light on genome gigantism in the Pamphagidae and provides insight into the correlation between genome size and repetitive sequences within the family.

Pericentromeric satellite RNAs as flexible protein partners in the regulation of nuclear structure

Pericentromeric heterochromatin is mainly composed of satellite DNA sequences. Although being historically associated with transcriptional repression, some pericentromeric satellite DNA sequences are transcribed. The transcription events of pericentromeric satellite sequences occur in highly flexible biological contexts. Hence, the apparent randomness of pericentromeric satellite transcription incites the discussion about the attribution of biological functions. However, pericentromeric satellite RNAs have clear roles in the organization of nuclear structure. Silencing pericentromeric heterochromatin depends on pericentromeric satellite RNAs, that, in a feedback mechanism, contribute to the repression of pericentromeric heterochromatin. Moreover, pericentromeric satellite RNAs can also act as scaffolding molecules in condensate subnuclear structures (e.g., nuclear stress bodies). Since the formation/dissociation of nuclear condensates provides cell adaptability, pericentromeric satellite RNAs can be an epigenetic platform for regulating (sub)nuclear structure. We review current knowledge about pericentromeric satellite RNAs that, irrespective of the meaning of biological function, should be functionally addressed in regular and disease settings. This article is categorized under: RNA Methods > RNA Analyses in Cells RNA in Disease and Development > RNA in Disease.

Genetic variation in recalcitrant repetitive regions of the Drosophila melanogaster genome

Many essential functions of organisms are encoded in highly repetitive genomic regions, including histones involved in DNA packaging, centromeres that are core components of chromosome segregation, ribosomal RNA comprising the protein translation machinery, telomeres that ensure chromosome integrity, piRNA clusters encoding host defenses against selfish elements, and virtually the entire Y chromosome. These regions, formed by highly similar tandem arrays, pose significant challenges for experimental and informatic study, impeding sequence-level descriptions essential for understanding genetic variation. Here, we report the assembly and variation analysis of such repetitive regions in Drosophila melanogaster, offering significant improvements to the existing community reference assembly. Our work successfully recovers previously elusive segments, including complete reconstructions of the histone locus and the pericentric heterochromatin of the X chromosome, spanning the Stellate locus to the distal flank of the rDNA cluster. To infer structural changes in these regions where alignments are often not practicable, we introduce landmark anchors based on unique variants that are putatively orthologous. These regions display considerable structural variation between different D. melanogaster strains, exhibiting differences in copy number and organization of homologous repeat units between haplotypes. In the histone cluster, although we observe minimal genetic exchange indicative of crossing over, the variation patterns suggest mechanisms such as unequal sister chromatid exchange. We also examine the prevalence and scale of concerted evolution in the histone and Stellate clusters and discuss the mechanisms underlying these observed patterns.

Comparative analysis, diversification, and functional validation of plant nucleotide-binding site domain genes

Nucleotide-binding site (NBS) domain genes are one of the superfamily of resistance genes involved in plant responses to pathogens. The current study identified 12,820 NBS-domain-containing genes across 34 species covering from mosses to monocots and dicots. These identified genes are classified into 168 classes with several novel domain architecture patterns encompassing significant diversity among plant species. Several classical (NBS, NBS-LRR, TIR-NBS, TIR-NBS-LRR, etc.) and species-specific structural patterns (TIR-NBS-TIR-Cupin_1-Cupin_1, TIR-NBS-Prenyltransf, Sugar_tr-NBS etc.) were discovered. We observed 603 orthogroups (OGs) with some core (most common orthogroups; OG0, OG1, OG2, etc.) and unique (highly specific to species; OG80, OG82, etc.) OGs with tandem duplications. The expression profiling presented the putative upregulation of OG2, OG6, and OG15 in different tissues under various biotic and abiotic stresses in susceptible and tolerant plants to cotton leaf curl disease (CLCuD). The genetic variation between susceptible (Coker 312) and tolerant (Mac7) Gossypium hirsutum accessions identified several unique variants in NBS genes of Mac7 (6583 variants) and Coker312 (5173 variants). The protein-ligand and proteins-protein interaction showed a strong interaction of some putative NBS proteins with ADP/ATP and different core proteins of the cotton leaf curl disease virus. The silencing of GaNBS (OG2) in resistant cotton through virus-induced gene silencing (VIGS) demonstrated its putative role in virus tittering. The presented study will be further helpful in understanding the plant adaptation mechanism.

Role of the repeat expansion size in predicting age of onset and severity in RFC1 disease

RFC1 disease, caused by biallelic repeat expansion in RFC1, is clinically heterogeneous in terms of age of onset, disease progression and phenotype. We investigated the role of the repeat size in influencing clinical variables in RFC1 disease. We also assessed the presence and role of meiotic and somatic instability of the repeat. In this study, we identified 553 patients carrying biallelic RFC1 expansions and measured the repeat expansion size in 392 cases. Pearson's coefficient was calculated to assess the correlation between the repeat size and age at disease onset. A Cox model with robust cluster standard errors was adopted to describe the effect of repeat size on age at disease onset, on age at onset of each individual symptoms, and on disease progression. A quasi-Poisson regression model was used to analyse the relationship between phenotype and repeat size. We performed multivariate linear regression to assess the association of the repeat size with the degree of cerebellar atrophy. Meiotic stability was assessed by Southern blotting on first-degree relatives of 27 probands. Finally, somatic instability was investigated by optical genome mapping on cerebellar and frontal cortex and unaffected peripheral tissue from four post-mortem cases. A larger repeat size of both smaller and larger allele was associated with an earlier age at neurological onset [smaller allele hazard ratio (HR) = 2.06, P < 0.001; larger allele HR = 1.53, P < 0.001] and with a higher hazard of developing disabling symptoms, such as dysarthria or dysphagia (smaller allele HR = 3.40, P < 0.001; larger allele HR = 1.71, P = 0.002) or loss of independent walking (smaller allele HR = 2.78, P < 0.001; larger allele HR = 1.60; P < 0.001) earlier in disease course. Patients with more complex phenotypes carried larger expansions [smaller allele: complex neuropathy rate ratio (RR) = 1.30, P = 0.003; cerebellar ataxia, neuropathy and vestibular areflexia syndrome (CANVAS) RR = 1.34, P < 0.001; larger allele: complex neuropathy RR = 1.33, P = 0.008; CANVAS RR = 1.31, P = 0.009]. Furthermore, larger repeat expansions in the smaller allele were associated with more pronounced cerebellar vermis atrophy (lobules I-V β = -1.06, P < 0.001; lobules VI-VII β = -0.34, P = 0.005). The repeat did not show significant instability during vertical transmission and across different tissues and brain regions. RFC1 repeat size, particularly of the smaller allele, is one of the determinants of variability in RFC1 disease and represents a key prognostic factor to predict disease onset, phenotype and severity. Assessing the repeat size is warranted as part of the diagnostic test for RFC1 expansion.

Anthocyanin gene enrichment in the distal region of cotton chromosome A07: mechanisms of reproductive organ coloration

Introduction

The biosynthesis of secondary metabolites like anthocyanins is often governed by metabolic gene clusters (MGCs) in the plant ancestral genome. However, the existence of gene clusters specifically regulating anthocyanin accumulation in certain organs is not well understood.

Methods and results

In this study, we identify MGCs linked to the coloration of cotton reproductive organs, such as petals, spots, and fibers. Through genetic analysis and map-based cloning, we pinpointed key genes on chromosome A07, such as PCC/GhTT19, which is involved in anthocyanin transport, and GbBM and GhTT2-3A, which are associated with the regulation of anthocyanin and proanthocyanidin biosynthesis. Our results demonstrate the coordinated control of anthocyanin and proanthocyanidin pathways, highlighting the evolutionary significance of MGCs in plant adaptation. The conservation of these clusters in cotton chromosome A07 across species underscores their importance in reproductive development and color variation. Our study sheds light on the complex biosynthesis and transport mechanisms for plant pigments, emphasizing the role of transcription factors and transport proteins in pigment accumulation.

Discussion

This research offers insights into the genetic basis of color variation in cotton reproductive organs and the potential of MGCs to enhance our comprehension of plant secondary metabolism.

Structural polymorphism of the nucleic acids in pentanucleotide repeats associated with the neurological disorder CANVAS

Short tandem repeats are inherently unstable during DNA replication depending on repeat length, and the expansion of the repeat length in the human genome is responsible for repeat expansion disorders. Pentanucleotide AAGGG and ACAGG repeat expansions in intron 2 of the gene encoding replication factor C subunit 1 (RFC1) cause cerebellar ataxia, neuropathy, vestibular areflexia syndrome (CANVAS) and other phenotypes of late-onset cerebellar ataxia. Herein, we reveal the structural polymorphism of the RFC1 repeats associated with CANVAS in vitro. Single-stranded AAGGG repeat DNA formed a hybrid-type G-quadruplex, whereas its RNA formed a parallel-type G-quadruplex with three layers. The RNA of the ACAGG repeat formed hairpin structure comprising C-G and G-C base pairs with A:A and GA:AG mismatched repeats. Furthermore, both pathogenic repeat RNAs formed more rigid structures than those of the nonpathogenic repeat RNAs. These findings provide novel insights into the structural polymorphism of the RFC1 repeats, which may be closely related to the disease mechanism of CANVAS.

Effect of tokenization on transformers for biological sequences

MOTIVATION

Deep-learning models are transforming biological research, including many bioinformatics and comparative genomics algorithms, such as sequence alignments, phylogenetic tree inference, and automatic classification of protein functions. Among these deep-learning algorithms, models for processing natural languages, developed in the natural language processing (NLP) community, were recently applied to biological sequences. However, biological sequences are different from natural languages, such as English, and French, in which segmentation of the text to separate words is relatively straightforward. Moreover, biological sequences are characterized by extremely long sentences, which hamper their processing by current machine-learning models, notably the transformer architecture. In NLP, one of the first processing steps is to transform the raw text to a list of tokens. Deep-learning applications to biological sequence data mostly segment proteins and DNA to single characters. In this work, we study the effect of alternative tokenization algorithms on eight different tasks in biology, from predicting the function of proteins and their stability, through nucleotide sequence alignment, to classifying proteins to specific families.

RESULTS

We demonstrate that applying alternative tokenization algorithms can increase accuracy and at the same time, substantially reduce the input length compared to the trivial tokenizer in which each character is a token. Furthermore, applying these tokenization algorithms allows interpreting trained models, taking into account dependencies among positions. Finally, we trained these tokenizers on a large dataset of protein sequences containing more than 400 billion amino acids, which resulted in over a 3-fold decrease in the number of tokens. We then tested these tokenizers trained on large-scale data on the above specific tasks and showed that for some tasks it is highly beneficial to train database-specific tokenizers. Our study suggests that tokenizers are likely to be a critical component in future deep-network analysis of biological sequence data.

AVAILABILITY AND IMPLEMENTATION

Code, data, and trained tokenizers are available on https://github.com/technion-cs-nlp/BiologicalTokenizers.

Genome-wide analysis of cellulose synthase gene superfamily in Tectona grandis L.f

This study aimed to explore Cellulose synthase gene superfamily of teak, and its evolutionary relationship with homologous genes of other woody species. The incidence of evolutionary events like gene duplication and gene loss, influence of the selection pressure, and consequent adaptive functional divergence of the duplicated TgCes gene were assessed alongside it's role in wood coloration. This study identified 39 full-length non-redundant proteins belonging to CesA and Csl gene families. TgCesA and TgCsl proteins with Cellulose synthase domain repeats indicated tandem gene duplication and probable genetic variability, enabling local adaptation. Further, multi-domain protein (MYB-like DNA-binding domain and CesA domain) with maximum introns was also identified indicating gene fusion and formation of complex protein with novel functions. Phylogenetic analysis grouped the genes into seven subfamilies (CesA, CslA, CslC, CslD, CslE, CslG, and CslM) with each undergoing gene duplication and loss along their evolutionary history. Post-species gene duplications and probable neofunctionalization were identified in TgCesA and TgCsl gene families. Each subfamily was found to be under strong purifying selection with a few or no sites under positive selection. Functional divergence analysis further revealed site-specific selective constraints in CesA and Csl genes of the teak Cellulose synthase gene family. Furthermore, protein-protein interaction network analysis identified co-expression of Cellulose synthase gene with flavonoid 3',5'-hydroxylase (F3'5'H, CYP75A), involved in the biosynthesis of xylem anthocyanin compounds, probably responsible for wood coloration. This study thus offers a foundation for future research in wood formation and wood property traits specific to teak and its provenances.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03927-6.

Fossil-calibrated phylogenies of Southern cave wētā show dispersal and extinction confound biogeographic signal

The biota of continents and islands are commonly considered to have a source-sink relationship, but small islands can harbour distinctive taxa. The distribution of four monotypic genera of Orthoptera on young subantarctic islands indicates a role for long-distance dispersal and extinction. Phylogenetic relationships were inferred from whole mtDNA genomes and nuclear sequences (45S cassette; four histones). We used a fossil and one palaeogeographic event to calibrate molecular clock analysis. We confirm that neither the Australian nor Aotearoa-New Zealand Rhaphidophoridae faunas are monophyletic. The radiation of Macropathinae may have begun in the late Jurassic, but trans-oceanic dispersal is required to explain the current distribution of some lineages within this subfamily. Dating the most recent common ancestor of seven island endemic species with their nearest mainland relative suggests that each existed long before their island home was available. Time estimates from our fossil-calibrated molecular clock analysis suggest several lineages have not been detected on mainland New Zealand, Australia, or elsewhere most probably due to their extinction, providing evidence that patterns of extinction, which are not consistently linked to range size or lineage age, confound biogeographic signal.

TBP facilitates RNA Polymerase I transcription following mitosis

The TATA-box binding protein (TBP) is the sole transcription factor common in the initiation complexes of the three major eukaryotic RNA Polymerases (Pol I, II and III). Although TBP is central to transcription by the three RNA Pols in various species, the emergence of TBP paralogs throughout evolution has expanded the complexity in transcription initiation. Furthermore, recent studies have emerged that questioned the centrality of TBP in mammalian cells, particularly in Pol II transcription, but the role of TBP and its paralogs in Pol I transcription remains to be re-evaluated. In this report, we show that in murine embryonic stem cells TBP localizes onto Pol I promoters, whereas the TBP paralog TRF2 only weakly associates to the Spacer Promoter of rDNA, suggesting that it may not be able to replace TBP for Pol I transcription. Importantly, acute TBP depletion does not fully disrupt Pol I occupancy or activity on ribosomal RNA genes, but TBP binding in mitosis leads to efficient Pol I reactivation following cell division. These findings provide a more nuanced role for TBP in Pol I transcription in murine embryonic stem cells.

Genome-wide characterization of microsatellites in cobia Rachycentron canadum (Linnaeus, 1766): Survey and analysis of their abundance and diversity

The cobia Rachycentron canadum, mainly distributed in the warm waters of tropical and subtropical regions around the world, remains a fish of considerable economic importance. Detailed diversity and the number of microsatellite sequences in the cobia genome are still unintelligible. The primary aim of this work was to identify and quantify the miscellaneous SSR sequences in the cobia genome. More than 280,000 sequences were sequenced and screened using next-generation sequencing technology and microsatellite identification. Perfect mononucleotide repeats, dinucleotide microsatellites, and trinucleotide microsatellites contain (A)10 /(T)10 , (AC)6 /(TG)6 , and (AAT)5-32 as the largest number of motifs in each type of microsatellite, respectively. The tetranucleotide and pentanucleotide microsatellites (TTM and PTM) consist of the largest number of motifs of both (ATCT)5-32 and (TCAT)5-31 in TTMs, and (CTCTC)5-9 in PTMs, whereas the hexanucleotide microsatellites are rarely observed in the cobia genome. All c. 38000 sequences of composite microsatellites are extremely diverse, including compound (11.71%), interrupted compound (71.77%), complex (0.45%), and interrupted complex (16.07%). In this study, we developed a convenient and useful recording system for writing down and categorizing diverse composite microsatellite types. This system will provide great support for exploring repeat origins, evolutionary mechanisms, and the application of polymorphic microsatellites.

Rex1BD and the 14-3-3 protein control heterochromatin organization at tandem repeats by linking RNAi and HDAC

Tandem DNA repeats are often organized into heterochromatin that is crucial for genome organization and stability. Recent studies revealed that individual repeats within tandem DNA repeats can behave very differently. How DNA repeats are assembled into distinct heterochromatin structures remains poorly understood. Here, we developed a genome-wide genetic screen using a reporter gene at different units in a repeat array. This screen led to identification of a conserved protein Rex1BD required for heterochromatin silencing. Our structural analysis revealed that Rex1BD forms a four-helix bundle structure with a distinct charged electrostatic surface. Mechanistically, Rex1BD facilitates the recruitment of Clr6 histone deacetylase (HDAC) by interacting with histones. Interestingly, Rex1BD also interacts with the 14-3-3 protein Rad25, which is responsible for recruiting the RITS (RNA-induced transcriptional silencing) complex to DNA repeats. Our results suggest that coordinated action of Rex1BD and Rad25 mediates formation of distinct heterochromatin structure at DNA repeats via linking RNAi and HDAC pathways.

Mwhitiwhiti Aotearoa: Phylogeny and synonymy of the silent alpine grasshopper radiation of New Zealand (Orthoptera: Acrididae)

Aotearoa New Zealand has a fauna of endemic alpine grasshoppers, consisting of thirteen species distributed among four genera. The many re-classifications of species within this group and the presence of species complexes highlight the uncertainty that surrounds relationships within and between these genera. High-throughput Next Generation Sequencing was used to assemble the complete mitochondrial genomes, 45S ribosomal cassettes and histone sequences of New Zealands four endemic alpine genera: Alpinacris, Brachaspis, Paprides and Sigaus. Phylogenetic analysis of these molecular datasets, as individual genes, partitions and combinations returned a consistent topology that is incompatible with the current classification. The genera Sigaus, Alpinacris, and Paprides all exhibit paraphyly. A consideration of the pronotum, epiphallus and terminalia of adult specimens reveals species-specific differences, but fails to provide compelling evidence for species groups justifying distinct genera. In combination with phylogenetic, morphological and spatial evidence we propose a simplified taxonomy consisting of a single genus for the mwhitiwhiti Aotearoa species radiation.

Insights into genome plasticity and gene regulation in Orientia tsutsugamushi through genome-wide mining of microsatellite markers

Microsatellite markers are being used for molecular identification and characterization as well as estimation of evolution patterns due to their highly polymorphic nature. The repeats hold 40% of the entire genome of Orientia tsutsugamushi (OT), but not yet been characterized. Thus, we investigated the genome-wide presence of microsatellites within nine complete genomes of OT and analyzed their distribution pattern, composition, and complexity. The in-silico study revealed that the genome of OT enriched with microsatellites having a total of 126,187 SSRs and 10,374 cSSRs throughout the genome, of which 70% and 30% are represented within the coding and non-coding regions, respectively. The relative density (RD) and relative abundance (RA) of SSRs were 42-44.43/kb and 6.25-6.59/kb, while for cSSRs this value ranged from 7.06 to 8.1/kb and 0.50 to 0.55/kb, respectively. However, RA and RD were weakly correlated with genome size and incidence of microsatellites. The mononucleotide repeats (54.55%) were prevalent over di- (33.22%), tri- (11.88%), tetra- (0.27%), penta- (0.02%), hexanucleotide (0.04%) repeats, with poly (A/T) richness over poly (G/C). The motif composition of cSSRs revealed that maximum cSSRs were made up of two microsatellites having unique duplication patterns such as AT-x-AT and CG-x-CG. To our knowledge, this is the first study of microsatellites in the OT genome, where characterization of such variations in repeat sequences would be important in deciphering the origin, rate of mutation, and role of repeat sequences in the genome. More numbers of microsatellites represented within the coding region provide an insight into the genome plasticity that may interfere with gene regulation to mitigate host-pathogen interaction and evolution of the species.

Regulation and function of R-loops at repetitive elements

R-loops are atypical, three-stranded nucleic acid structures that contain a stretch of RNA:DNA hybrids and an unpaired, single stranded DNA loop. R-loops are physiological relevant and can act as regulators of gene expression, chromatin structure, DNA damage repair and DNA replication. However, unscheduled and persistent R-loops are mutagenic and can mediate replication-transcription conflicts, leading to DNA damage and genome instability if left unchecked. Detailed transcriptome analysis unveiled that 85% of the human genome, including repetitive regions, hold transcriptional activity. This anticipates that R-loops management plays a central role for the regulation and integrity of genomes. This function is expected to have a particular relevance for repetitive sequences that make up to 75% of the human genome. Here, we review the impact of R-loops on the function and stability of repetitive regions such as centromeres, telomeres, rDNA arrays, transposable elements and triplet repeat expansions and discuss their relevance for associated pathological conditions.

Characterization of the marsupial endogenous retrovirus walb with a focus on satellite DNA formation

The walbRep megasatellite DNA found in the red-necked wallaby was formed from the walb endogenous retrovirus. Our previous PCR experiments suggested the presence of walb and absence of walbRep in the genome of the tammar wallaby, which diverged from the red-necked wallaby 2-3 Mya. The results failed to exclude the possibility that certain walbRep sequences might have remained undetected owing to variation in the primer-annealing regions; therefore, the aforementioned suggestion was not confirmed. To obtain conclusive evidence, we analyzed the structure of walb sequences drawn from the tammar wallaby genome database recently updated to a chromosome-level assembly. All walb copies existed as separate DNA segments, not constituting tandem repeats. We concluded that walbRep was formed in the red-necked wallaby lineage after its divergence from the tammar wallaby. We also confirm the presence of a walb copy with an anomalistic, complex structure and propose a simple model for its generation mechanism.

Repetitive DNA sequence detection and its role in the human genome

Repetitive DNA sequences playing critical roles in driving evolution, inducing variation, and regulating gene expression. In this review, we summarized the definition, arrangement, and structural characteristics of repeats. Besides, we introduced diverse biological functions of repeats and reviewed existing methods for automatic repeat detection, classification, and masking. Finally, we analyzed the type, structure, and regulation of repeats in the human genome and their role in the induction of complex diseases. We believe that this review will facilitate a comprehensive understanding of repeats and provide guidance for repeat annotation and in-depth exploration of its association with human diseases.

Satellite RNAs: emerging players in subnuclear architecture and gene regulation

Satellite DNA is characterized by long, tandemly repeated sequences mainly found in centromeres and pericentromeric chromosomal regions. The recent advent of telomere-to-telomere sequencing data revealed the complete sequences of satellite regions, including centromeric α-satellites and pericentromeric HSat1-3, which together comprise ~ 5.7% of the human genome. Despite possessing constitutive heterochromatin features, these regions are transcribed to produce long noncoding RNAs with highly repetitive sequences that associate with specific sets of proteins to play various regulatory roles. In certain stress or pathological conditions, satellite RNAs are induced to assemble mesoscopic membraneless organelles. Specifically, under heat stress, nuclear stress bodies (nSBs) are scaffolded by HSat3 lncRNAs, which sequester hundreds of RNA-binding proteins. Upon removal of the stressor, nSBs recruit additional regulatory proteins, including protein kinases and RNA methylases, which modify the previously sequestered nSB components. The sequential recruitment of substrates and enzymes enables nSBs to efficiently regulate the splicing of hundreds of pre-mRNAs under limited temperature conditions. This review discusses the structural features and regulatory roles of satellite RNAs in intracellular architecture and gene regulation.

Microprobes for Label-Free Detection of Short Tandem Repeats: An Insight into Alleviating Secondary Structure Effects

Overgrowth of short tandem repeat sequences in our genes can cause various neurodegenerative disorders. Such repeat sequences are ideal targets for the label-free electrochemical detection of such potential expansions. However, their length- and sequence-dependent secondary structures may interfere with the interfacial charge transfer of a detection platform, making them complex targets. In addition, the gene contains sporadic repeats that may result in false-positive signals. Therefore, it is necessary to design a platform capable of mitigating these effects and ultimately enhancing the specificity of tandem repeats. Here, we analyzed three different backbones of nucleic acid microprobes [DNA, peptide nucleic acid, and lock-nucleic acid (LNA)] to detect in vitro transcribed RNA carrying CAG repeats, which are associated with Huntington's disease, based on the charge-transfer resistance of the interface. We found that the LNA microprobe can distinguish lengths down to the attomolar concentration level and alleviate the effect of secondary structures and sporadic repeats in the sequence, thus distinguishing the "tandem repeats" specifically. Additionally, the control experiments conducted with and without Mg2+ demonstrated the LNA microprobe to perform better in the presence of the divalent cation. The results suggest that the LNA-based platform may eventually lead to the development of a reliable and straightforward biosensor for genetic neurodegenerative disorders.

Are scattered microsatellites weak chromosomal markers? Guided mapping reveals new insights into Trachelyopterus (Siluriformes: Auchenipteridae) diversity

The scattered distribution pattern of microsatellites is a challenging problem in fish cytogenetics. This type of array hinders the identification of useful patterns and the comparison between species, often resulting in over-limited interpretations that only label it as "scattered" or "widely distributed". However, several studies have shown that the distribution pattern of microsatellites is non-random. Thus, here we tested whether a scattered microsatellite could have distinct distribution patterns on homeologous chromosomes of closely related species. The clustered sites of 18S and 5S rDNA, U2 snRNA and H3/H4 histone genes were used as a guide to compare the (GATA)n microsatellite distribution pattern on the homeologous chromosomes of six Trachelyopterus species: T. coriaceus and Trachelyopterus aff. galeatus from the Araguaia River basin; T. striatulus, T. galeatus and T. porosus from the Amazonas River basin; and Trachelyopterus aff. coriaceus from the Paraguay River basin. Most species had similar patterns of the (GATA)n microsatellite in the histone genes and 5S rDNA carriers. However, we have found a chromosomal polymorphism of the (GATA)n sequence in the 18S rDNA carriers of Trachelyopterus galeatus, which is in Hard-Weinberg equilibrium and possibly originated through amplification events; and a chromosome polymorphism in Trachelyopterus aff. galeatus, which combined with an inversion polymorphism of the U2 snRNA in the same chromosome pair resulted in six possible cytotypes, which are in Hardy-Weinberg disequilibrium. Therefore, comparing the distribution pattern on homeologous chromosomes across the species, using gene clusters as a guide to identify it, seems to be an effective way to further the analysis of scattered microsatellites in fish cytogenetics.

Genetic Diversity and Structure Revealed by Genomic Microsatellite Markers of Mytilus unguiculatus in the Coast of China Sea

The hard-shelled mussel Mytilus unguiculatus plays an important role in mussel aquaculture in China due to its characteristic and nutritive value. In this study, 10 microsatellite loci are used to study the genetic diversity and genetic structure of seven location populations of M. unguiculatus in coastal areas of China. The results of amplification and genotyping show that the observed heterozygosity (H_o) and the expected heterozygosity (H_e) are 0.61~0.71 and 0.72~0.83, respectively. M. unguiculatus has high genetic diversity. The inbreeding index (F_IS) of M. unguiculatus is significantly positive (F_IS: 0.14~0.19), indicating that inbreeding might exist within populations. The genetic structure of M. unguiculatus is weak within populations from the East China Sea All results showed that genetic differences existed between the Qingdao population from the Yellow Sea and other populations from the East China Sea. It does not detect a population bottleneck event or expansion event in the populations. The results from this study can be used to provide important insights in genetic management units and sustainable utilization of M. unguiculatus resources and provide a better understand of genetic structure of marine bivalve with similar planktonic larval stage in the China Sea.

Overexpression of Pericentromeric HSAT2 DNA Increases Expression of EMT Markers in Human Epithelial Cancer Cell Lines

Pericentromeric tandemly repeated DNA of human satellites 1, 2, and 3 (HS1, HS2, and HS3) is actively transcribed in some cells. However, the functionality of the transcription remains obscure. Studies in this area have been hampered by the absence of a gapless genome assembly. The aim of our study was to map a transcript that we have previously described as HS2/HS3 on chromosomes using a newly published gapless genome assembly T2T-CHM13, and create a plasmid overexpressing the transcript to assess the influence of HS2/HS3 transcription on cancer cells. We report here that the sequence of the transcript is tandemly repeated on nine chromosomes (1, 2, 7, 9, 10, 16, 17, 22, and Y). A detailed analysis of its genomic localization and annotation in the T2T-CHM13 assembly revealed that the sequence belonged to HSAT2 (HS2) but not to the HS3 family of tandemly repeated DNA. The transcript was found on both strands of HSAT2 arrays. The overexpression of the HSAT2 transcript increased the transcription of the genes encoding the proteins involved in the epithelial-to-mesenchymal transition, EMT (SNAI1, ZEB1, and SNAI2), and the genes that mark cancer-associated fibroblasts (VIM, COL1A1, COL11A1, and ACTA2) in cancer cell lines A549 and HeLa. Co-transfection of the overexpression plasmid and antisense nucleotides eliminated the transcription of EMT genes observed after HSAT2 overexpression. Antisense oligonucleotides also decreased transcription of the EMT genes induced by tumor growth factor beta 1 (TGFβ1). Thus, our study suggests HSAT2 lncRNA transcribed from the pericentromeric tandemly repeated DNA is involved in EMT regulation in cancer cells.

Clonal hematopoiesis driven by chromosome 1q/MDM4 trisomy defines a canonical route toward leukemia in Fanconi anemia

Fanconi anemia (FA) patients experience chromosome instability, yielding hematopoietic stem/progenitor cell (HSPC) exhaustion and predisposition to poor-prognosis myeloid leukemia. Based on a longitudinal cohort of 335 patients, we performed clinical, genomic, and functional studies in 62 patients with clonal evolution. We found a unique pattern of somatic structural variants and mutations that shares features of BRCA-related cancers, the FA-hallmark being unbalanced, microhomology-mediated translocations driving copy-number alterations. Half the patients developed chromosome 1q gain, driving clonal hematopoiesis through MDM4 trisomy downmodulating p53 signaling later followed by secondary acute myeloid lukemia genomic alterations. Functionally, MDM4 triplication conferred greater fitness to murine and human primary FA HSPCs, rescued inflammation-mediated bone marrow failure, and drove clonal dominance in FA mouse models, while targeting MDM4 impaired leukemia cells in vitro and in vivo. Our results identify a linear route toward secondary leukemogenesis whereby early MDM4-driven downregulation of basal p53 activation plays a pivotal role, opening monitoring and therapeutic prospects.

Evolution of satDNAs on holocentric chromosomes: insights from hemipteran insects of the genus Mahanarva

Satellite DNAs (satDNAs) constitute one of the main components of eukaryote genomes and are involved in chromosomal organization and diversification. Although largely studied, little information was gathered about their evolution on holocentric species, i.e., diffuse centromeres, which, due to differences in repeat organization, could result in different evolutionary patterns. Here, we combined bioinformatics and cytogenetic approaches to evaluate the evolution of the satellitomes in Mahanarva holocentric insects. In two species, de novo identification revealed a high number of satDNAs, 110 and 113, with an extreme monomer length range of 18-4228 bp. The overall abundance of satDNAs was observed to be 6.67% in M. quadripunctata and 1.98% in M. spectabilis, with different abundances for the shared satDNAs. Chromosomal mapping of the most abundant repeats of M. quadripunctata and M. spectabilis on other Mahanarva reinforced the dynamic nature of satDNAs. Variable patterns of chromosomal distribution for the satDNAs were noticed, with the occurrence of clusters on distinct numbers of chromosomes and at different positions and the occurrence of scattered signals or nonclustered satDNAs. Altogether, our data demonstrated the high dynamism of satDNAs in Mahanarva with the involvement of this genomic fraction in chromosome diversification of the genus. The general characteristics and patterns of evolution of satDNAs are similar to those observed on monocentric chromosomes, suggesting that the differential organization of genome compartments observed on holocentric chromosomes compared with monocentric chromosomes does not have a large impact on the evolution of satDNAs. Analysis of the satellitomes of other holocentric species in a comparative manner will shed light on this issue.

Discovery of Macrocycle-Based HPK1 Inhibitors for T-Cell-Based Immunotherapy

Hematopoietic progenitor kinase 1 (HPK1) is a negative regulator of T-cell activation, and targeting HPK1 is considered a promising strategy for improving responses to antitumor immune therapies. The biggest challenge of HPK1 inhibitor design is to achieve a higher selectivity to GLK, an HPK1 homology protein as a positive regulator of T-cell activation. Herein, we report the design of a series of macrocycle-based HPK1 inhibitors via a conformational constraint strategy. The identified candidate compound 5i exhibited HPK1 inhibition with an IC₅₀ value of 0.8 nM and 101.3-fold selectivity against GLK. Compound 5i also displayed good oral bioavailability (F = 27-49%) in mice and beagles and favorable metabolic stability (T_1/2 > 186.4 min) in human liver microsomes. More importantly, compound 5i demonstrated a clear synergistic effect with anti-PD-1 in both MC38 (MSI) and CT26 (MSS) syngeneic tumor mouse models. These results showed that compound 5i has a great potential in immunotherapy.

Synthetical lethality of Werner helicase and mismatch repair deficiency is mediated by p53 and PUMA in colon cancer

Synthetic lethality is a powerful approach for targeting oncogenic drivers in cancer. Recent studies revealed that cancer cells with microsatellite instability (MSI) require Werner (WRN) helicase for survival; however, the underlying mechanism remains unclear. In this study, we found that WRN depletion strongly induced p53 and its downstream apoptotic target PUMA in MSI colorectal cancer (CRC) cells. p53 or PUMA deletion abolished apoptosis induced by WRN depletion in MSI CRC cells. Importantly, correction of MSI abrogated the activation of p53/PUMA and cell killing, while induction of MSI led to sensitivity in isogenic CRC cells. Rare p53-mutant MSI CRC cells are resistant to WRN depletion due to lack of PUMA induction, which could be restored by wildtype (WT) p53 knock in or reconstitution. WRN depletion or treatment with the RecQ helicase inhibitor ML216 suppressed in vitro and in vivo growth of MSI CRCs in a p53/PUMA-dependent manner. ML216 treatment was efficacious in MSI CRC patient-derived xenografts. Interestingly, p53 gene remains WT in the majority of MSI CRCs. These results indicate a critical role of p53/PUMA-mediated apoptosis in the vulnerability of MSI CRCs to WRN loss, and support WRN as a promising therapeutic target in p53-WT MSI CRCs.

Time-resolved microfluidics unravels individual cellular fates during double-strand break repair

BACKGROUND

Double-strand break repair (DSBR) is a highly regulated process involving dozens of proteins acting in a defined order to repair a DNA lesion that is fatal for any living cell. Model organisms such as Saccharomyces cerevisiae have been used to study the mechanisms underlying DSBR, including factors influencing its efficiency such as the presence of distinct combinations of microsatellites and endonucleases, mainly by bulk analysis of millions of cells undergoing repair of a broken chromosome. Here, we use a microfluidic device to demonstrate in yeast that DSBR may be studied at a single-cell level in a time-resolved manner, on a large number of independent lineages undergoing repair.

RESULTS

We used engineered S. cerevisiae cells in which GFP is expressed following the successful repair of a DSB induced by Cas9 or Cpf1 endonucleases, and different genetic backgrounds were screened to detect key events leading to the DSBR efficiency. Per condition, the progenies of 80-150 individual cells were analyzed over 24 h. The observed DSBR dynamics, which revealed heterogeneity of individual cell fates and their contributions to global repair efficacy, was confronted with a coupled differential equation model to obtain repair process rates. Good agreement was found between the mathematical model and experimental results at different scales, and quantitative comparisons of the different experimental conditions with image analysis of cell shape enabled the identification of three types of DSB repair events previously not recognized: high-efficacy error-free, low-efficacy error-free, and low-efficacy error-prone repair.

CONCLUSIONS

Our analysis paves the way to a significant advance in understanding the complex molecular mechanism of DSB repair, with potential implications beyond yeast cell biology. This multiscale and multidisciplinary approach more generally allows unique insights into the relation between in vivo microscopic processes within each cell and their impact on the population dynamics, which were inaccessible by previous approaches using molecular genetics tools alone.

A PCR-based method for the diagnosis of Enterobius vermicularis in stool samples, specifically designed for clinical application

BACKGROUND

Enterobius vermicularis (E. vermicularis) is a nematode that infects up to 200 million people worldwide, despite effective medications being available. Conventional diagnostic tests are hindered by low sensitivity and poor patient compliance. Furthermore, no biomolecular techniques are available for clinical application. The aim of this study was to develop a procedure specifically designed for clinical application to detect E. vermicularis by means of PCR.

MATERIALS AND METHODS

Two subject groups were taken into account: a group of 27 infected patients and a control group of 27 healthy subjects. A nested-PCR was performed on fecal samples to detect E. vermicularis. Due to the intrinsic difficulties of the fecal matrix, several countermeasures were adopted to ensure the efficient performance of the method: (a) a large amount of feces for the extraction process (20 g instead of 200 mg); (b) a combination of chemical and physical treatments to grind the fecal matrix; (c) an additional purification process for the negative samples after the first nested-PCR; and (d) the selection of a very specific target region for the PCR.

RESULTS

Due to the lack of overlap with other organisms, a sequence of the 5S ribosomal DNA (rDNA) spacer region including the tract SL1 was chosen to design appropriate external and internal primers. The first nested-PCR detected E.vermicularis in 19/27 samples from infected patients. After further purification, 5/8 of the negative samples resulted positive at the second PCR. Conversely, all the samples from healthy controls resulted negative to both PCRs. Sensitivity and specificity of the method were, respectively, 88.9% and 100%.

CONCLUSION

The results prove the high diagnostic accuracy of the proposed method, addressing and overcoming the challenges posed by both conventional tests and PCR-based approaches. Therefore, the method can be proposed for clinical application.

Evolution of the ancestral mammalian karyotype and syntenic regions

Decrypting the rearrangements that drive mammalian chromosome evolution is critical to understanding the molecular bases of speciation, adaptation, and disease susceptibility. Using 8 scaffolded and 26 chromosome-scale genome assemblies representing 23/26 mammal orders, we computationally reconstructed ancestral karyotypes and syntenic relationships at 16 nodes along the mammalian phylogeny. Three different reference genomes (human, sloth, and cattle) representing phylogenetically distinct mammalian superorders were used to assess reference bias in the reconstructed ancestral karyotypes and to expand the number of clades with reconstructed genomes. The mammalian ancestor likely had 19 pairs of autosomes, with nine of the smallest chromosomes shared with the common ancestor of all amniotes (three still conserved in extant mammals), demonstrating a striking conservation of synteny for ∼320 My of vertebrate evolution. The numbers and types of chromosome rearrangements were classified for transitions between the ancestral mammalian karyotype, descendent ancestors, and extant species. For example, 94 inversions, 16 fissions, and 14 fusions that occurred over 53 My differentiated the therian from the descendent eutherian ancestor. The highest breakpoint rate was observed between the mammalian and therian ancestors (3.9 breakpoints/My). Reconstructed mammalian ancestor chromosomes were found to have distinct evolutionary histories reflected in their rates and types of rearrangements. The distributions of genes, repetitive elements, topologically associating domains, and actively transcribed regions in multispecies homologous synteny blocks and evolutionary breakpoint regions indicate that purifying selection acted over millions of years of vertebrate evolution to maintain syntenic relationships of developmentally important genes and regulatory landscapes of gene-dense chromosomes.

Megasatellite formation and evolution in vertebrate genes

Since formation of the first proto-eukaryotes, gene repertoire and genome complexity have significantly increased. Among genetic elements responsible for this increase are tandem repeats. Here we describe a genome-wide analysis of large tandem repeats, called megasatellites, in 58 vertebrate genomes. Two bursts occurred, one after the radiation between Agnatha and Gnathostomata fishes and the second one in therian mammals. Megasatellites are enriched in subtelomeric regions and frequently encoded in genes involved in transcription regulation, intracellular trafficking, and cell membrane metabolism, reminiscent of what is observed in fungus genomes. The presence of many introns within young megasatellites suggests that an exon-intron DNA segment is first duplicated and amplified before accumulation of mutations in intronic parts partially erases the megasatellite in such a way that it becomes detectable only in exons. Our results suggest that megasatellite formation and evolution is a dynamic and still ongoing process in vertebrate genomes.

Mutated KLF4(K409Q) in meningioma binds STRs and activates FGF3 gene expression

Krüppel-like factor 4 (KLF4) is a transcription factor that has been proven necessary for both induction and maintenance of pluripotency and self-renewal. Whole-genome sequencing defined a unique mutation in KLF4 (KLF4K409Q) in human meningiomas. However, the molecular mechanism of this tumor-specific KLF4 mutation is unknown. Using genome-wide high-throughput and focused quantitative transcriptional approaches in human cell lines, primary meningeal cells, and meningioma tumor tissue, we found that a change in the evolutionarily conserved DNA-binding domain of KLF4 alters its DNA recognition preference, resulting in a shift in downstream transcriptional activity. In the KLF4K409Q-specific targets, the normally silent fibroblast growth factor 3 (FGF3) is activated. We demonstrated a neomorphic function of KLF4K409Q in stimulating FGF3 transcription through binding to its promoter and in using short tandem repeats (STRs) located within the locus as enhancers.

Genomic Organization of Microsatellites and LINE-1-like Retrotransposons: Evolutionary Implications for Ctenomys minutus (Rodentia: Ctenomyidae) Cytotypes

Simple Summary

In animals, several species contain substantial chromosomal and genomic variation among their populations, but as to what could have driven such diversification is still a puzzle for most cases. Here, we used molecular cytogenetic analysis to expose the main genomic elements involved in the population variation observed in the Neotropical underground rodents of the genus Ctenomys (Rodentia: Ctenomyidae), which harbor the most significant chromosomal variation among mammals (2n = 10 to 2n = 70). These data provide evidence for a correlation between repetitive genomic content and localization of evolutionary breakpoint regions (EBRs) and highlight their direct impact in promoting chromosomal rearrangements.

Abstract

The Neotropical underground rodents of the genus Ctenomys (Rodentia: Ctenomyidae) comprise about 65 species, which harbor the most significant chromosomal variation among mammals (2n = 10 to 2n = 70). Among them, C. minutus stands out with 45 different cytotypes already identified, among which, seven parental ones, named A to G, are parapatrically distributed in the coastal plains of Southern Brazil. Looking for possible causes that led to such extensive karyotype diversification, we performed chromosomal mapping of different repetitive DNAs, including microsatellites and long interspersed element-1 (LINE-1) retrotransposons in the seven parental cytotypes. Although microsatellites were found mainly in the centromeric and telomeric regions of the chromosomes, different patterns occur for each cytotype, thus revealing specific features. Likewise, the LINE-1-like retrotransposons also showed a differential distribution for each cytotype, which may be linked to stochastic loss of LINE-1 in some populations. Here, microsatellite motifs (A)₃₀, (C)₃₀, (CA)₁₅, (CAC)₁₀, (CAG)₁₀, (CGG)₁₀, (GA)₁₅, and (GAG)₁₀ could be mapped to fusion of chromosomes 20/17, fission and inversion in the short arm of chromosome 2, fusion of chromosomes 23/19, and different combinations of centric and tandem fusions of chromosomes 22/24/16. These data provide evidence for a correlation between repetitive genomic content and localization of evolutionary breakpoints and highlight their direct impact in promoting chromosomal rearrangements.

Breaking satellite silence: human satellite II RNA expression in ovarian cancer

Multiple cancer types demonstrate abnormal expression of repetitive RNA sequences as a form of epigenetic instability. There is growing interest in understanding the role of repetitive RNAs in cancer pathogenesis and immunogenicity and in their potential role as diagnostic or therapeutic biomarkers. In this issue of the JCI, Porter and colleagues report on satellite RNA in a subset of ovarian cancers. The authors found that high expression of human satellite (HSAT) repeats — but not other families of repeats — was associated with an immunosuppressive phenotype in ovarian cancer cell lines and tumor samples. Further induction of HSAT RNA levels in vitro, surprisingly, leads to innate immune activation, suggesting a potential therapeutic strategy. This work highlights the expanding role of repetitive RNAs in tumor biology and the need to better define specific classes of repetitive elements expressed in cancer — as well as their role in tumorigenesis, tumor immunity, and the host response to cancer.

Role of an FNIP Repeat Domain-Containing Protein Encoded by Megavirus Baoshan during Viral Infection

FNIP repeat domain-containing protein (FNIP protein) is a little-studied atypical leucine-rich repeat domain-containing protein found in social amoebae and mimiviruses. Here, a recently reported mimivirus of lineage C, Megavirus baoshan, was analyzed for FNIP protein genes. A total of 82 FNIP protein genes were identified, each containing up to 26 copies of the FNIP repeat, and mostly having an F-box domain at the N terminus. Both nucleotide and amino acid sequences of FNIP repeat were highly conserved. Most of the FNIP protein genes clustered together tandemly in groups of two to 14 genes. Nearly all FNIP protein genes shared similar expression patterns and were expressed 4 to 9 h postinfection. A typical viral FNIP protein, Mb0983, was selected for functional analysis. Protein interactome analysis identified two small GTPases, Rap1B and Rab7A, that interacted with Mb0983 in cytoplasm. The overexpression of Mb0983 in Acanthamoeba castellanii accelerated the degradation of Rap1B and Rab7A during viral infection. Mb0983 also interacted with host SKP1 and cullin-1, which were conserved components of the SKP1-cullin-1-F-box protein (SCF)-type ubiquitin E3 ligase complex. Deletion of the F-box domain of Mb0983 not only abolished its interaction with SKP1 and cullin-1 but also returned the speed of Rap1B and Rab7A degradation to normal in infected A. castellanii. These results suggested that Mb0983 is a part of the SCF-type ubiquitin E3 ligase complex and plays a role in the degradation of Rap1B and Rab7A. They also implied that other viral F-box-containing FNIP proteins might have similar effects on various host proteins. IMPORTANCE Megavirus baoshan encodes 82 FNIP proteins, more than any other reported mimiviruses. Their genetic and transcriptional features suggest that they are important for virus infection and adaption. Since most mimiviral FNIP proteins have the F-box domain, they were predicted to be involved in protein ubiquitylation. FNIP protein Mb0983 interacted with host SKP1 and cullin-1 through the F-box domain, supporting the idea that it is a part of the SCF-type ubiquitin E3 ligase complex. The substrates of Mb0983 for degradation were identified as the host small GTPases Rap1B and Rab7A. Combining the facts of the presence of a large number of FNIP genes in megavirus genomes, the extremely high expression level of the viral ubiquitin gene, and the reported observation that 35% of megavirus-infected amoeba cells died without productive infection, it is likely that megavirus actively explores the host ubiquitin-proteasome pathway in infection and that viral FNIP proteins play roles in the process.

A Structural Potential of Rare Trinucleotide Repeat Tracts in RNA

Among types of trinucleotide repeats, there is some disproportion in the frequency of their occurrence in the human exome. This research presents new data describing the folding and thermodynamic stability of short, tandem RNA repeats of 23 types, focusing on the rare, yet poorly analyzed ones. UV-melting experiments included the presence of PEG or potassium and magnesium ions to determine their effect on the stability of RNA repeats structures. Rare repeats predominantly stayed single-stranded but had the potential for base pairing with other partially complementary repeat tracts. A coexistence of suitably complementary repeat types in a single RNA creates opportunities for interaction in the context of the secondary structure of RNA. We searched the human transcriptome for model RNAs in which different, particularly rare trinucleotide repeats coexist and selected the GABRA4 and CHIC1 RNAs to study intramolecular interactions between the repeat tracts that they contain. In vitro secondary structure probing results showed that the UAA and UUG repeat tracts, present in GABRA4 3' UTR, form a double helix, which separates one of its structural domains. For the RNA CHIC1 ORF fragment containing four short AGG repeat tracts and the CGU tract, we proved the formation of quadruplexes that blocked reverse transcription.

Complete chloroplast genome of Gentianopsis barbata and comparative analysis with related species from Gentianaceae

Gentianopsis barbata is an essential medicinal plant in China with high ornamental and medicinal values. Unfortunately, the study of the chloroplast genome of this plant still has a gap. This study sequenced and characterized the complete chloroplast genome of G. barbata. The complete chloroplast genome of G. barbata is a typical circular structure with 151,123 bp. It consists of a large single-copy region (82,690 bp) and a small single-copy region (17,887 bp) separated by a pair of inverted repeats (25,273 bp), which covers 78 protein-coding genes, 30 tRNAs, and 4 rRNAs. Repeat analysis showed the highest frequency of palindrome. Thirty-seven simple sequence repeats were identified, most of which were single nucleotides. The bayesian inference tree, maximum likelihood tree, and neighbor joining tree suggested that G. barbata is grouped with Gentianopsis grandis and Gentianopsis paludosa. The divergence time analysis showed that G. barbata diverged at 1.243 Mya. Comparative chloroplast analysis can reveal interspecific diversity, and regions with high variation can be used to develop molecular markers applicable to various research areas. Our results provide new insight into plastome evolution and valuable resource for further studies on G. barbata.

Functional variability in adhesion and flocculation of yeast megasatellite genes

Megasatellites are large tandem repeats found in all fungal genomes but especially abundant in the opportunistic pathogen Candida glabrata. They are encoded in genes involved in cell-cell interactions, either between yeasts or between yeast and human cells. In the present work, we have been using an iterative genetic system to delete several Candida glabrata megasatellite-containing genes and found that 2 of them were positively involved in adhesion to epithelial cells, whereas 3 genes negatively controlled adhesion. Two of the latter, CAGL0B05061g or CAGL0A04851g, were also negative regulators of yeast-to-yeast adhesion, making them central players in controlling Candida glabrata adherence properties. Using a series of synthetic Saccharomyces cerevisiae strains in which the FLO1 megasatellite was replaced by other tandem repeats of similar length but different sequences, we showed that the capacity of a strain to flocculate in liquid culture was unrelated to its capacity to adhere to epithelial cells or to invade agar. Finally, to understand how megasatellites were initially created and subsequently expanded, an experimental evolution system was set up, in which modified yeast strains containing different megasatellite seeds were grown in bioreactors for more than 200 generations and selected for their ability to sediment at the bottom of the culture tube. Several flocculation-positive mutants were isolated. Functionally relevant mutations included general transcription factors as well as a 230-kbp segmental duplication.

Genetic and Epigenetic Interplay Define Disease Onset and Severity in Repeat Diseases

Repeat diseases, such as fragile X syndrome, myotonic dystrophy, Friedreich ataxia, Huntington disease, spinocerebellar ataxias, and some forms of amyotrophic lateral sclerosis, are caused by repetitive DNA sequences that are expanded in affected individuals. The age at which an individual begins to experience symptoms, and the severity of disease, are partially determined by the size of the repeat. However, the epigenetic state of the area in and around the repeat also plays an important role in determining the age of disease onset and the rate of disease progression. Many repeat diseases share a common epigenetic pattern of increased methylation at CpG islands near the repeat region. CpG islands are CG-rich sequences that are tightly regulated by methylation and are often found at gene enhancer or insulator elements in the genome. Methylation of CpG islands can inhibit binding of the transcriptional regulator CTCF, resulting in a closed chromatin state and gene down regulation. The downregulation of these genes leads to some disease-specific symptoms. Additionally, a genetic and epigenetic interplay is suggested by an effect of methylation on repeat instability, a hallmark of large repeat expansions that leads to increasing disease severity in successive generations. In this review, we will discuss the common epigenetic patterns shared across repeat diseases, how the genetics and epigenetics interact, and how this could be involved in disease manifestation. We also discuss the currently available stem cell and mouse models, which frequently do not recapitulate epigenetic patterns observed in human disease, and propose alternative strategies to study the role of epigenetics in repeat diseases.

Satellitome Analysis and Transposable Elements Comparison in Geographically Distant Populations of Spodoptera frugiperda

Spodoptera frugiperda (fall armyworm) is a member of the superfamily Noctuoidea that accounts for more than a third of all Lepidoptera and includes a considerable number of agricultural and forest pest species. Spodoptera frugiperda is a polyphagous species that is a significant agricultural pest worldwide, emphasizing its economic importance. Spodoptera frugiperda's genome size, assembly, phylogenetic classification, and transcriptome analysis have all been previously described. However, the different studies reported different compositions of repeated DNA sequences that occupied the whole assembled genome, and the Spodoptera frugiperda genome also lacks the comprehensive study of dynamic satellite DNA. We conducted a comparative analysis of repetitive DNA across geographically distant populations of Spodoptera frugiperda, particularly satellite DNA, using publicly accessible raw genome data from eight different geographical regions. Our results showed that most transposable elements (TEs) were commonly shared across all geographically distant samples, except for the Maverick and PIF/Harbinger elements, which have divergent repeat copies. The TEs age analysis revealed that most TEs families consist of young copies 1-15 million years old; however, PIF/Harbinger has some older/degenerated copies of 30-35 million years old. A total of seven satellite DNA families were discovered, accounting for approximately 0.65% of the entire genome of the Spodoptera frugiperda fall armyworm. The repeat profiling analysis of satellite DNA families revealed differential read depth coverage or copy numbers. The satellite DNA families range in size from the lowest 108 bp SfrSat06-108 families to the largest (1824 bp) SfrSat07-1824 family. We did not observe a statistically significant correlation between monomer length and K2P divergence, copy number, or abundance of each satellite family. Our findings suggest that the satellite DNA families identified in Spodoptera frugiperda account for a considerable proportion of the genome's repetitive fraction. The satellite DNA families' repeat profiling revealed a point mutation along the reference sequences. Limited TEs differentiation exists among geographically distant populations of Spodoptera frugiperda.

Satellitome analysis illuminates the evolution of ZW sex chromosomes of Triportheidae fishes (Teleostei: Characiformes)

Satellites are an abundant source of repetitive DNAs that play an essential role in the chromosomal organization and are tightly linked with the evolution of sex chromosomes. Among fishes, Triportheidae stands out as the only family where almost all species have a homeologous ZZ/ZW sex chromosomes system. While the Z chromosome is typically conserved, the W is always smaller, with variations in size and morphology between species. Here, we report an analysis of the satellitome of Triportheus auritus (TauSat) by integrating genomic and chromosomal data, with a special focus on the highly abundant and female-biased satDNAs. In addition, we investigated the evolutionary trajectories of the ZW sex chromosomes in the Triportheidae family by mapping satDNAs in selected representative species of this family. The satellitome of T. auritus comprised 53 satDNA families of which 24 were also hybridized by FISH. Most satDNAs differed significantly between sexes, with 19 out of 24 being enriched on the W chromosome of T. auritus. The number of satDNAs hybridized into the W chromosomes of T. signatus and T. albus decreased to six and four, respectively, in accordance with the size of their W chromosomes. No TauSat probes produced FISH signals on the chromosomes of Agoniates halecinus. Despite its apparent conservation, our results indicate that each species differs in the satDNA accumulation on the Z chromosome. Minimum spanning trees (MSTs), generated for three satDNA families with different patterns of FISH mapping data, revealed different homogenization rates between the Z and W chromosomes. These results were linked to different levels of recombination between them. The most abundant satDNA family (TauSat01) was exclusively hybridized in the centromeres of all 52 chromosomes of T. auritus, and its putative role in the centromere evolution was also highlighted. Our results identified a high differentiation of both ZW chromosomes regarding satellites composition, highlighting their dynamic role in the sex chromosomes evolution.

Neurodegenerative diseases associated with non-coding CGG tandem repeat expansions

Non-coding CGG repeat expansions cause multiple neurodegenerative disorders, including fragile X-associated tremor/ataxia syndrome, neuronal intranuclear inclusion disease, oculopharyngeal myopathy with leukodystrophy, and oculopharyngodistal myopathy. The underlying genetic causes of several of these diseases have been identified only in the past 2-3 years. These expansion disorders have substantial overlapping clinical, neuroimaging and histopathological features. The shared features suggest common mechanisms that could have implications for the development of therapies for this group of diseases - similar therapeutic strategies or drugs may be effective for various neurodegenerative disorders induced by non-coding CGG expansions. In this Review, we provide an overview of clinical and pathological features of these CGG repeat expansion diseases and consider the likely pathological mechanisms, including RNA toxicity, CGG repeat-associated non-AUG-initiated translation, protein aggregation and mitochondrial impairment. We then discuss future research needed to improve the identification and diagnosis of CGG repeat expansion diseases, to improve modelling of these diseases and to understand their pathogenesis. We also consider possible therapeutic strategies. Finally, we propose that CGG repeat expansion diseases may represent manifestations of a single underlying neuromyodegenerative syndrome in which different organs are affected to different extents depending on the gene location of the repeat expansion.