Transposable elements in Drosophila montana from harsh cold environments.

BACKGROUND: Substantial discoveries during the past century have revealed that transposable elements (TEs) can play a crucial role in genome evolution by affecting gene expression and inducing genetic rearrangements, among other molecular and structural effects. Yet, our knowledge on the role of TEs in adaptation to extreme climates is still at its infancy. The availability of long-read sequencing has opened up the possibility to identify and study potential functional effects of TEs with higher precision. In this work, we used Drosophila montana as a model for cold-adapted organisms to study the association between TEs and adaptation to harsh climates. RESULTS: Using the PacBio long-read sequencing technique, we de novo identified and manually curated TE sequences in five Drosophila montana genomes from eco-geographically distinct populations. We identified 489 new TE consensus sequences which represented 92% of the total TE consensus in D. montana. Overall, 11-13% of the D. montana genome is occupied by TEs, which as expected are non-randomly distributed across the genome. We identified five potentially active TE families, most of them from the retrotransposon class of TEs. Additionally, we found TEs present in the five analyzed genomes that were located nearby previously identified cold tolerant genes. Some of these TEs contain promoter elements and transcription binding sites. Finally, we detected TEs nearby fixed and polymorphic inversion breakpoints. CONCLUSIONS: Our research revealed a significant number of newly identified TE consensus sequences in the genome of D. montana, suggesting that non-model species should be studied to get a comprehensive view of the TE repertoire in Drosophila species and beyond. Genome annotations with the new D. montana library allowed us to identify TEs located nearby cold tolerant genes, and present at high population frequencies, that contain regulatory regions and are thus good candidates to play a role in D. montana cold stress response. Finally, our annotations also allow us to identify for the first time TEs present in the breakpoints of three D. montana inversions.

Comparative analysis of amphibian genomes: An emerging resource for basic and applied research.

Amphibians are the most threatened group of vertebrates and are in dire need of conservation intervention to ensure their continued survival. They exhibit unique features including a high diversity of reproductive strategies, permeable and specialized skin capable of producing toxins and antimicrobial compounds, multiple genetic mechanisms of sex determination and in some lineages, the ability to regenerate limbs and organs. Although genomic approaches would shed light on these unique traits and aid conservation, sequencing and assembly of amphibian genomes has lagged behind other taxa due to their comparatively large genome sizes. Fortunately, the development of long-read sequencing technologies and initiatives has led to a recent burst of new amphibian genome assemblies. Although growing, the field of amphibian genomics suffers from the lack of annotation resources, tools for working with challenging genomes and lack of high-quality assemblies in multiple clades of amphibians. Here, we analyse 51 publicly available amphibian genomes to evaluate their usefulness for functional genomics research. We report considerable variation in genome assembly quality and completeness and report some of the highest transposable element and repeat contents of any vertebrate. Additionally, we detected an association between transposable element content and climatic variables. Our analysis provides evidence of conserved genome synteny despite the long divergence times of this group, but we also highlight inconsistencies in chromosome naming and orientation across genome assemblies. We discuss sequencing gaps in the phylogeny and suggest key targets for future sequencing endeavours. Finally, we propose increased investment in amphibian genomics research to promote their conservation.

Mobilome impacts on physiology in the widely used non-toxic mutant Microcystis aeruginosa PCC 7806 ΔmcyB and toxic wildtype.

The Microcystis mobilome is a well-known but understudied component of this bloom-forming cyanobacterium. Through genomic and transcriptomic comparisons, we found five families of transposases that altered the expression of genes in the well-studied toxigenic type-strain, Microcystis aeruginosa PCC 7086, and a non-toxigenic genetic mutant, Microcystis aeruginosa PCC 7806 ΔmcyB. Since its creation in 1997, the ΔmcyB strain has been used in comparative physiology studies against the wildtype strain by research labs throughout the world. Some differences in gene expression between what were thought to be otherwise genetically identical strains have appeared due to insertion events in both intra- and intergenic regions. In our ΔmcyB isolate, a sulfate transporter gene cluster (sbp-cysTWA) showed differential expression from the wildtype, which may have been caused by the insertion of a miniature inverted repeat transposable element (MITE) in the sulfate-binding protein gene (sbp). Differences in growth in sulfate-limited media also were also observed between the two isolates. This paper highlights how Microcystis strains continue to "evolve" in lab conditions and illustrates the importance of insertion sequences / transposable elements in shaping genomic and physiological differences between Microcystis strains thought otherwise identical. This study forces the necessity of knowing the complete genetic background of isolates in comparative physiological experiments, to facilitate the correct conclusions (and caveats) from experiments.

Widespread transposable element dysregulation in human aging brains with Alzheimer's disease.

INTRODUCTION: Transposable element (TE) dysregulation is associated with neuroinflammation in Alzheimer's disease (AD) brains. Yet, TE quantitative trait loci (teQTL) have not been well characterized in human aged brains with AD. METHODS: We leveraged large-scale bulk and single-cell RNA sequencing, whole-genome sequencing (WGS), and xQTL from three human AD brain biobanks to characterize TE expression dysregulation and experimentally validate AD-associated TEs using CRISPR interference (CRISPRi) assays in human induced pluripotent stem cell (iPSC)-derived neurons. RESULTS: We identified 26,188 genome-wide significant TE expression QTLs (teQTLs) in human aged brains. Subsequent colocalization analysis of teQTLs with AD genetic loci identified AD-associated teQTLs and linked locus TEs. Using CRISPRi assays, we pinpointed a neuron-specific suppressive role of the activated short interspersed nuclear element (SINE; chr11:47608036-47608220) on expression of C1QTNF4 via reducing neuroinflammation in human iPSC-derived neurons. DISCUSSION: We identified widespread TE dysregulation in human AD brains and teQTLs offer a complementary analytic approach to identify likely AD risk genes. HIGHLIGHTS: Widespread transposable element (TE) dysregulations are observed in human aging brains with degrees of neuropathology, apolipoprotein E (APOE) genotypes, and neuroinflammation in Alzheimer's disease (AD). A catalog of TE quantitative trait loci (teQTLs) in human aging brains was created using matched RNA sequencing and whole-genome sequencing data. CRISPR interference assays reveal that an upregulated intergenic TE from the MIR family (chr11: 47608036-47608220) suppresses expression of its nearest anti-inflammatory gene C1QTNF4 in human induced pluripotent stem cell-derived neurons.

Syntenic lncRNA locus exhibits DNA regulatory functions with sequence evolution.

Syntenic long non-coding RNAs (lncRNAs) often show limited sequence conservation across species, prompting concern in the field. This study delves into functional signatures of syntenic lncRNAs between humans and zebrafish. Syntenic lncRNAs are highly expressed in zebrafish, with ∼90 % located near protein-coding genes, either in sense or antisense orientation. During early zebrafish development and in human embryonic stem cells (H1-hESC), syntenic lncRNA loci are enriched with cis-regulatory repressor signatures, influencing the expression of development-associated genes. In later zebrafish developmental stages and specific human cell lines, these syntenic lncRNA loci function as enhancers or transcription start sites (TSS) for protein-coding genes. Analysis of transposable elements (TEs) in syntenic lncRNA sequences revealed intriguing patterns: human lncRNAs are enriched in simple repeat elements, while their zebrafish counterparts show enrichment in LTR elements. This sequence evolution likely arises from post-rearrangement mutations that enhance DNA elements or cis-regulatory functions. It may also contribute to vertebrate innovation by creating novel transcription factor binding sites within the locus. This study highlights the conserved functionality of syntenic lncRNA loci through DNA elements, emphasizing their conserved roles across species despite sequence divergence.

A T cell receptor specific for an HLA-A*03:01-restricted epitope in the endogenous retrovirus ERV-K-Env exhibits limited recognition of its cognate epitope.

Transposable elements (TEs) are often expressed at higher levels in tumor cells than normal cells, implicating these genomic regions as an untapped pool of tumor-associated antigens. In ovarian cancer (OC), protein from the TE ERV-K is frequently expressed by tumor cells. Here we determined whether the targeting of previously identified epitope in the envelope gene (env) of ERV-K resulted in target antigen specificity against cancer cells. We found that transducing healthy donor T cells with an ERV-K-Env-specific T cell receptor construct resulted in antigen specificity only when co-cultured with HLA-A*03:01 B lymphoblastoid cells. Furthermore, in vitro priming of several healthy donors with this epitope of ERV-K-Env did not result in target antigen specificity. These data suggest that the T cell receptor is a poor candidate for targeting this specific ERV-K-Env epitope and has limited potential as a T cell therapy for OC.

Characterization of hAT DNA transposon superfamily in the genome of Neotropical fish Apareiodon sp.

DNA transposons are diverse in fish genomes and have been described to generate genomic evolutionary novelties. hAT transposable element data are scarce in Teleostei genomes, making it challenging to conduct comparative genomic studies to understand their neutrality or function. This study aimed to perform a genomic and molecular characterization of hAT copies to assess the diversity of these elements and associate changes in these sequences to genomic and karyotypic novelties in Apareiodon sp. The data revealed that hAT TEs are highly abundant in the Apareiodon sp. genome, with few possibly autonomous copies. Highly conserved sequences with likely functional transposases were observed in nine hAT elements. A great diversity of hAT subgroups was observed, especially from Ac, Charlie, Blackjack, Tip100, hAT6, and hAT5, and a similar wave of hAT genomic invasion was identified in the genome for these six groups of hAT sequences. The data also revealed a distinct number of microsatellites within degenerated hAT copies. hAT sites were demonstrated to be dispersed in the Apareiodon sp. chromosomes and not involved in W chromosome-specific region differentiation. In conclusion, the genomic analysis revealed a great diversity of hAT elements, possible autonomous copies, and differentiation of degenerated transposable elements into tandem sequences.

Comparative Genomic Analyses of Colletotrichum lindemuthianum Pathotypes with Different Virulence Levels and Lifestyles.

Colletotrichum lindemuthianum is the most frequent pathogenic fungus of the common bean Phaseolus vulgaris. This filamentous fungus employs a hemibiotrophic nutrition/infection strategy, which is characteristic of many Colletotrichum species. Due to host-pathogen coevolution, C. lindemuthianum includes pathotypes with a diversity of virulence against differential common bean varieties. In this study, we performed comparative genomic analyses on three pathotypes with different virulence levels and a non-pathogenic pathotype, isolated from different geographical areas in Mexico. Our results revealed large genomes with high transposable element contents that have undergone expansions, generating intraspecific diversity. All the pathotypes exhibited a similar number of clusters of orthologous genes (COGs) and Gene Ontology (GO) terms. TFomes contain families that are typical in fungal genomes; however, they show different contents between pathotypes, mainly in transcription factors with the fungal-specific TF and Zn2Cys6 domains. Peptidase families mainly contain abundant serine peptidases, metallopeptidases, and cysteine peptidases. In the secretomes, the number of genes differed between the pathotypes, with a high percentage of candidate effectors. Both the virulence gene and CAZyme gene content for each pathotype was abundant and diverse, and the latter was enriched in hemicellulolytic enzymes. We provide new insights into the nature of intraspecific diversity among C. lindemuthianum pathotypes and the origin of their ability to rapidly adapt to genetic changes in its host and environmental conditions.

A reference genome for the Harpy Eagle reveals steady demographic decline and chromosomal rearrangements in the origin of Accipitriformes.

The Harpy Eagle (Harpia harpyja) is an iconic species that inhabits forested landscapes in Neotropical regions, with decreasing population trends mainly due to habitat loss, and currently classified as vulnerable. Here, we report on a chromosome-scale genome assembly for a female individual combining long reads, optical mapping, and chromatin conformation capture reads. The final assembly spans 1.35 Gb, with N50scaffold equal to 58.1 Mb and BUSCO completeness of 99.7%. We built the first extensive transposable element (TE) library for the Accipitridae to date and identified 7,228 intact TEs. We found a burst of an unknown TE ~ 13-22 million years ago (MYA), coincident with the split of the Harpy Eagle from other Harpiinae eagles. We also report a burst of solo-LTRs and CR1 retrotransposons ~ 31-33 MYA, overlapping with the split of the ancestor to all Harpiinae from other Accipitridae subfamilies. Comparative genomics with other Accipitridae, the closely related Cathartidae and Galloanserae revealed major chromosome-level rearrangements at the basal Accipitriformes genome, in contrast to a conserved ancient genome architecture for the latter two groups. A historical demography reconstruction showed a rapid decline in effective population size over the last 20,000 years. This reference genome serves as a crucial resource for future conservation efforts towards the Harpy Eagle.

Characterization of active transposable elements and their new insertions in tuber propagated greater yam (Dioscorea alata).

BACKGROUND: Greater yam is a key staple crop grown in tropical and subtropical regions, while its asexual propagation mode had led to non-flowering mutations. How transposable elements contribute to its genetic variations is rarely analyzed. We used transcriptome and whole genome sequencing data to identify active transposable elements (TEs) and genetic variation caused by these active TEs. Our aim was to shed light on which TEs would lead to its intraspecies variation. RESULTS: Annotation of de novo assembly transcripts indicated that 0.8 - 0.9% of transcripts were TE related, with LTR retrotransposons (LTR-RTs) accounted for 65% TE transcripts. A large portion of these transcripts were non-autonomous TEs, which had incomplete functional domains. The majority of mapped transcripts were distributed in genic deficient regions, with 9% of TEs overlapping with genic regions. Moreover, over 90% TE transcripts exhibited low expression levels and insufficient reads coverage to support full-length structure assembly. Subfamily analysis of Copia and Gypsy, the two LTR-RTs revealed that a small number of subfamilies contained a significantly larger number of members, which play a key role in generating TE transcript. Based on resequencing data, 15,002 L-RT insertion loci were detected for active LTR-RT members. The insertion loci of LTR-RTs were highly divergent among greater yam accessions. CONCLUSIONS: This study showed the ongoing transcription and transpositions of TEs in greater yam, despite low transcription levels and incomplete proteins insufficient for autonomous transposition. While our research did not directly link these TEs to specific yam traits such as tuber yield and propagation mode, it lays a crucial foundation for further research on how these TE insertion polymorphisms (TIPs) might be related to variations in greater yam traits and its tuber propagation mode. Future research may explore the potential roles of TEs in trait variations, such as tuber yield and stress resistance, in greater yam.

Genomes of Aegilops umbellulata provide new insights into unique structural variations and genetic diversity in the U-genome for wheat improvement.

Aegilops umbellulata serve as an important reservoir for novel biotic and abiotic stress tolerance for wheat improvement. However, chromosomal rearrangements and evolutionary trajectory of this species remain to be elucidated. Here, we present a comprehensive investigation into Ae. umbellulata genome by generating a high-quality near telomere-to-telomere genome assembly of PI 554389 and resequencing 20 additional Ae. umbellulata genomes representing diverse geographical and phenotypic variations. Our analysis unveils complex chromosomal rearrangements, most prominently in 4U and 6U chromosomes, delineating a distinct evolutionary trajectory of Ae. umbellulata from wheat and its relatives. Furthermore, our data rectified the erroneous naming of chromosomes 4U and 6U in the past and highlighted multiple major evolutionary events that led to the present-day U-genome. Resequencing of diverse Ae. umbellulata accessions revealed high genetic diversity within the species, partitioning into three distinct evolutionary sub-populations and supported by extensive phenotypic variability in resistance against several races/pathotypes of five major wheat diseases. Disease evaluations indicated the presence of several novel resistance genes in the resequenced lines for future studies. Resequencing also resulted in the identification of six new haplotypes for Lr9, the first resistance gene cloned from Ae. umbellulata. The extensive genomic and phenotypic resources presented in this study will expedite the future genetic exploration of Ae. umbellulata, facilitating efforts aimed at enhancing resiliency and productivity in wheat.

Investigation of chimeric transcripts derived from LINE-1 and Alu retrotransposons in cerebellar tissues of individuals with autism spectrum disorder (ASD).

LINE-1 and Alu retrotransposons are components of the human genome and have been implicated in many human diseases. These elements can influence human transcriptome plasticity in various mechanisms. Chimeric transcripts derived from LINE-1 and Alu can also impact the human transcriptome, such as exonization and post-transcriptional modification. However, its specific role in ASD neuropathology remains unclear, particularly in the cerebellum tissues. We performed RNA-sequencing of post-mortem cerebellum tissues from ASD and unaffected individuals for transposable elements profiling and chimeric transcript identification. The majority of free transcripts of transposable elements were not changed in the cerebellum tissues of ASD compared with unaffected individuals. Nevertheless, we observed that chimeric transcripts derived from LINE-1 and Alu were embedded in the transcripts of differentially expressed genes in the cerebellum of ASD, and these genes were related to developments and abnormalities of the cerebellum. In addition, the expression levels of these genes were correlated with the significantly decreased thickness of the molecular layer in the cerebellum of ASD. We also found that global methylation and expression of LINE-1 and Alu elements were not changed in ASD, but observed in the ASD sub-phenotypes. Our findings showed associations between transposable elements and cerebellar abnormalities in ASD, particularly in distinct phenotypic subgroups. Further investigations using appropriate models are warranted to elucidate the structural and functional implications of LINE-1 and Alu elements in ASD neuropathology.

Three transposable elements exhibiting differential expression in pre-eclampsia overlap with enhancer regions.

Transposable elements (TEs) play a crucial role in placental development and dysfunction. Our study examined TE expression in pre-eclampsia (PE) using RNA-seq datasets. We identified differentially expressed TEs and explored the genomic location of the most significant TEs, investigating their possible regulatory roles. Notably, three TEs overlapped with putative enhancer regions, suggesting a potential regulatory impact on gene expression. These findings highlight the regulatory potential of TEs and their importance in placental development, supporting that TE dysregulation may contribute to PE pathogenesis.

Comparative Analysis of Transposable Elements in the Genomes of Citrus and Citrus-Related Genera.

Transposable elements (TEs) significantly contribute to the evolution and diversity of plant genomes. In this study, we explored the roles of TEs in the genomes of Citrus and Citrus-related genera by constructing a pan-genome TE library from 20 published genomes of Citrus and Citrus-related accessions. Our results revealed an increase in TE content and the number of TE types compared to the original annotations, as well as a decrease in the content of unclassified TEs. The average length of TEs per assembly was approximately 194.23 Mb, representing 41.76% (Murraya paniculata) to 64.76% (Citrus gilletiana) of the genomes, with a mean value of 56.95%. A significant positive correlation was found between genome size and both the number of TE types and TE content. Consistent with the difference in mean whole-genome size (39.83 Mb) between Citrus and Citrus-related genera, Citrus genomes contained an average of 34.36 Mb more TE sequences than Citrus-related genomes. Analysis of the estimated insertion time and half-life of long terminal repeat retrotransposons (LTR-RTs) suggested that TE removal was not the primary factor contributing to the differences among genomes. These findings collectively indicate that TEs are the primary determinants of genome size and play a major role in shaping genome structures. Principal coordinate analysis (PCoA) of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) identifiers revealed that the fragmented TEs were predominantly derived from ancestral genomes, while intact TEs were crucial in the recent evolutionary diversification of Citrus. Moreover, the presence or absence of intact TEs near the AdhE superfamily was closely associated with the bitterness trait in the Citrus species. Overall, this study enhances TE annotation in Citrus and Citrus-related genomes and provides valuable data for future genetic breeding and agronomic trait research in Citrus.

Centromeric transposable elements and epigenetic status drive karyotypic variation in the eastern hoolock gibbon.

Great apes have maintained a stable karyotype with few large-scale rearrangements; in contrast, gibbons have undergone a high rate of chromosomal rearrangements coincident with rapid centromere turnover. Here we characterize assembled centromeres in the Eastern hoolock gibbon, Hoolock leuconedys (HLE), finding a diverse group of transposable elements (TEs) that differ from the canonical alpha satellites found across centromeres of other apes. We find that HLE centromeres contain a CpG methylation centromere dip region, providing evidence this epigenetic feature is conserved in the absence of satellite arrays; nevertheless, we report a variety of atypical centromeric features, including protein-coding genes and mismatched replication timing. Further, large structural variations define HLE centromeres and distinguish them from other gibbons. Combined with differentially methylated TEs, topologically associated domain boundaries, and segmental duplications at chromosomal breakpoints, we propose that a "perfect storm" of multiple genomic attributes with propensities for chromosome instability shaped gibbon centromere evolution.

Identification of novel myelodysplastic syndromes prognostic subgroups by integration of inflammation, cell-type composition, and immune signatures in the bone marrow.

Mutational profiles of myelodysplastic syndromes (MDS) have established that a relatively small number of genetic aberrations, including SF3B1 and SRSF2 spliceosome mutations, lead to specific phenotypes and prognostic subgrouping. We performed a multi-omics factor analysis (MOFA) on two published MDS cohorts of bone marrow mononuclear cells (BMMNCs) and CD34 + cells with three data modalities (clinical, genotype, and transcriptomics). Seven different views, including immune profile, inflammation/aging, retrotransposon (RTE) expression, and cell-type composition, were derived from these modalities to identify the latent factors with significant impact on MDS prognosis. SF3B1 was the only mutation among 13 mutations in the BMMNC cohort, indicating a significant association with high inflammation. This trend was also observed to a lesser extent in the CD34 + cohort. Interestingly, the MOFA factor representing the inflammation shows a good prognosis for MDS patients with high inflammation. In contrast, SRSF2 mutant cases show a granulocyte-monocyte progenitor (GMP) pattern and high levels of senescence, immunosenescence, and malignant myeloid cells, consistent with their poor prognosis. Furthermore, MOFA identified RTE expression as a risk factor for MDS. This work elucidates the efficacy of our integrative approach to assess the MDS risk that goes beyond all the scoring systems described thus far for MDS.

Natural Transposable Element Insertions Contribute to Host Fitness in Model Yeasts.

Transposable elements (TEs) are ubiquitous in the eukaryote genomes, but their evolutionary and functional significance remains largely obscure and contentious. Here, we explore the evolution and functional impact of TEs in two model unicellular eukaryotes, the fission yeast Schizosaccharomyces pombe and the budding yeast Saccharomyces cerevisiae, which diverged around 330 to 420 million years ago. We analyze the distribution of LTR retrotransposons (LTR-RTs, the only TE order identified in both species) and their solo-LTR derivatives in 35 strains of S. pombe and 128 strains of S. cerevisiae. We find that natural LTR-RT and solo-LTR insertions exhibit high presence-absence polymorphism among individuals in both species. Population genetics analyses show that solo-LTR insertions experienced functional constraints similar to synonymous sites of host genes in both species, indicating a majority of solo-LTR insertions might have evolved in a neutral manner. When knocking out nine representative solo-LTR insertions separately in the S. pombe strain 972h- and 12 representative solo-LTR insertions separately in the S. cerevisiae strain S288C, we find that one solo-LTR insertion in S. pombe has a significant effect on the fitness and transcriptome of its host. Together, our findings indicate that a fraction of natural TE insertions likely shape their host transcriptomes and thereby contribute to their host fitness, with implications for understanding the functional significance of TEs in eukaryotes.

Impact of Chromosomal Fusion and Transposable Elements on the Genomic Evolution and Genetic Diversity of Ilex Species.

The genus Ilex belongs to the sole family and is the single genus within the order Aquifoliales, exhibiting significant phenotypic diversity. However, the genetic differences underlying these phenotypic variations have rarely been studied. In this study, collinearity analyses of three Ilex genomes, Ilex latifolia Thunb., Ilex polyneura (Hand.-Mazz.) S. Y. Hu, and Ilex asprella Champ. ex Benth., indicated a recent fusion event contributing to the reduction of chromosomes in I. asprella. Comparative genome analyses showed slight differences in gene annotation among the three species, implying a minimal disruption of genes following chromosomal fusion in I. asprella. Comprehensive annotation of transposable elements (TEs) revealed that TEs constitute a significant portion of the Ilex genomes, with LTR transposons being predominant. TEs exhibited an inverse relationship with gene density, potentially influencing gene regulation and chromosomal architecture. TE insertions were shown to affect the conformation and binding sites of key genes such as 7-deoxyloganetin glucosyltransferase and transmembrane kinase (TMK) genes, highlighting potential functional impacts. The structural variations caused by TE insertions suggest significant roles in the evolutionary dynamics, leading to either loss or gain of gene function. This study underscores the importance of TEs in shaping the genomic landscape and evolutionary trajectories of Ilex species.

3t-seq: automatic gene expression analysis of single-copy genes, transposable elements, and tRNAs from RNA-seq data.

RNA sequencing is the gold-standard method to quantify transcriptomic changes between two conditions. The overwhelming majority of data analysis methods available are focused on polyadenylated RNA transcribed from single-copy genes and overlook transcripts from repeated sequences such as transposable elements (TEs). These self-autonomous genetic elements are increasingly studied, and specialized tools designed to handle multimapping sequencing reads are available. Transfer RNAs are transcribed by RNA polymerase III and are essential for protein translation. There is a need for integrated software that is able to analyze multiple types of RNA. Here, we present 3t-seq, a Snakemake pipeline for integrated differential expression analysis of transcripts from single-copy genes, TEs, and tRNA. 3t-seq produces an accessible report and easy-to-use results for downstream analysis starting from raw sequencing data and performing quality control, genome mapping, gene expression quantification, and statistical testing. It implements three methods to quantify TEs expression and one for tRNA genes. It provides an easy-to-configure method to manage software dependencies that lets the user focus on results. 3t-seq is released under MIT license and is available at https://github.com/boulardlab/3t-seq.

Repeatome evolution across space and time: Unravelling repeats dynamics in the plant genus Erythrostemon Klotzsch (Leguminosae Juss).

Fluctuations in genomic repetitive fractions (repeatome) are known to impact several facets of evolution, such as ecological adaptation and speciation processes. Therefore, investigating the divergence of repetitive elements can provide insights into an important evolutionary force. However, it is not clear how the different repetitive element clades are impacted by the different factors such as ecological changes and/or phylogeny. To discuss this, we used the Neotropical legume genus Erythrostemon (Caesalpinioideae) as a model, given its ancient origin (~33 Mya), lineage-specific niche conservatism, macroecological heterogeneity, and disjunct distribution in Meso- and South American (MA and SA respectively) lineages. We performed a comparative repeatomic analysis of 18 Erythrostemon species to test the impact of environmental variables over repeats diversification. Overall, repeatome composition was diverse, with high abundances of satDNAs and Ty3/gypsy-Tekay transposable elements, predominantly in the MA and SA lineages respectively. However, unexpected repeatome profiles unrelated to the phylogeny/biogeography were found in a few MA (E. coccineus, E. pannosus and E. placidus) and SA (E. calycinus) species, related to reticulate evolution and incongruence between nuclear and plastid topology, suggesting ancient hybridizations. The plesiomorphic Tekay and satDNA pattern was altered in the MA-sensu stricto subclade with a striking genomic differentiation (expansion of satDNA and retraction of Tekay) associated with the colonization of a new environment in Central America around 20 Mya. Our data reveal that the current species-specific Tekay pool was the result of two bursts of amplification probably in the Miocene, with distinct patterns for the MA and SA repeatomes. This suggests a strong role of the Tekay elements as modulators of the genome-environment interaction in Erythrostemon, providing macroevolutionary insights about mechanisms of repeatome differentiation and plant diversification across space and time.