Distribution Of Transposable Elements Research Articles

Decoding the Genomic Landscape of Pomegranate: A Genome-Wide Analysis of Transposable Elements and Their Structural Proximity to Functional Genes

Transposable elements (TEs) significantly drive dynamic changes that characterize genome evolution. However, understanding the variability associated with TE insertions among different cultivars remains challenging. The pomegranate (Punica granatum L.) has yet to be extensively studied regarding the roles of TEs in the diversification of cultivars. Herein, we explored the genome distribution of TEs and its potential functional implications among four pomegranate cultivars, ‘Bhagwa’, ‘Dabenzi’, ‘Taishanhong’ and ‘Tunisia’, whose genome sequences are available. A total of 8404 full-length TEs were isolated. The content of TEs varied among the cultivars, ranging from 41.67% of ‘Taishanhong’ to 52.45% of ‘Bhagwa’. In all cultivars, the Gypsy superfamily of retrotransposons accounted for a larger genome proportion than the Copia superfamily. Seventy-three full-length TEs were found at the same genomic loci in all four cultivars. By contrast, 947, 297, 311, and 874 TEs were found exclusively in ‘Bhagwa’, ‘Dabenzi’, ‘Taishanhong’, and ‘Tunisia’ cultivars, respectively. Phylogenetic clustering based on the presence of TE insertions in specific loci reflected the geographic origins of the cultivars. The insertion time profiles of LTR-REs were studied in the four cultivars. Shared elements across the four cultivars exhibited, on average, a more ancient insertion date than those exclusive to three, two, or one cultivars. The majority of TEs were located within 1000 bp from the nearest gene. This localization was observed for 57% of DNA TEs and 55% of long-terminal repeat retrotransposons (LTR-RE). More than 10% of TEs resulted inserted within genes. Concerning DNA TEs, 3.91% of insertions occurred in introns, while 2.42% occurred in exons. As to LTR-REs, 4% of insertions occurred in exons and 1.98% in introns. Functional analysis of the genes lying close to TEs was performed to infer if differences in TE insertion can affect the fruit quality. Two TE insertions were found close to two genes encoding 4-coumarate--CoA ligase, an enzyme involved in the phenylpropanoid pathway. Moreover, a TIR/Mariner element was found within the exon of a gene encoding anthocyanidin reductase in the ‘Tunisia’ genotype, crucial in the biosynthesis of flavan-3-ols and proanthocyanidins, strictly correlated with the nutraceutical properties of pomegranate. Although functional and metabolomic studies are essential to elucidate the consequences of TE insertions, these results contribute to advancing our comprehension of the role of TEs in pomegranate genomics, providing insights for crop breeding.

Relationship of MicroRNAs with Transposable Elements in the Type 1 Diabetes Development

The review article describes the involvement of epigenetic factors in type 1 diabetes mellitus (T1DM) etiopathogenesis. The disease is characterized by changes in expression of microRNAs that affect the transcription of genes involved in autoimmune reactions, destruction of beta cells and insulin production. However, the cause of the observed epigenetic changes is still unclear. In evolution, the sources of microRNA genes are transposable elements, which occupy up to 45 % of the entire human DNA sequence and are drivers of epigenetic regulation in ontogenesis. They are sources of transcription factor sequences and binding sites for them. Features of the genome distribution of transposable elements can cause changes in the number of 5’VNTR (variable number of tandem repeats) — repeats of insulin promoter region and HERV insertions into HLA genes, which affects their expression. Therefore, I assume that the cause of the development of type 1 diabetes mellitus may be an imbalance in transcription activation of transposons, which contributes to changes in the expression of specific microRNAs and protein-coding genes, and also contributes to autoimmune response development. Triggers for this may be individual features of genome distribution of transposons, viral infections and stress. An analysis of the scientific literature confirms my proposed mechanisms for T1DM development, since the global role of retroelements in hormonal regulation, the sensitivity of transposable elements to exogenous viral infections and stress, and HERV-W expression of the majority of patients with T1DM with activation of the autoimmune response have been proven. Analysis of the MDTE DB (miRNAs derived from transposable elements database) database showed the transposon origin of 12 T1DM-associated microRNAs (miR-192, miR-224, miR-31, miR-320c, miR-326, miR-340, miR-342, miR-44661, miR-548c, miR-652, miR-95), the use of which can become the basis for targeted therapy for T1DM.

Modeling transposable elements dynamics during polyploidization in plants

In this work we study the proliferation of transposable elements (TEs) and the epigenetic response of plants during the process of polyploidization. Through a deterministic model, expanding on our previous work on TE proliferation under epigenetic regulation, we study the long-term TE distribution and TE stability in the subgenomes of both autopolyploids and allopolyploids. We also explore different small-interfering RNA (siRNA) action modes on the subgenomes, including a model where siRNAs are not directed to specific genomes and one where siRNAs are directed – i.e. more active – in subgenomes with higher TE loads. In the autopolyploid case, we find long-term stable equilbria that tend to equilibrate the number of active TEs between subgenomes. In the allopolyploid case, directed siRNA action is fundamental to avoid a ”winner takes all” outcome of the competition between the TEs in the different subgenomes. We also show that decaying oscillations in the number of TEs occur naturally in all cases, perhaps explaining some of the observed features of ’genomic shock’ after hybridization events, and that the balance in the dynamics of the different types of siRNA is determinant for the synchronization of these oscillations.

The genomic distribution of transposable elements is driven by spatially variable purifying selection.

It is widely accepted that the genomic distribution of transposable elements (TEs) mainly reflects the outcome of purifying selection and insertion bias (1). Nevertheless, the relative importance of these two evolutionary forces could not be tested thoroughly. Here, we introduce an experimental system, which allows separating purifying selection from TE insertion bias. We used experimental evolution to study the TE insertion patterns in Drosophila simulans founder populations harboring 1040 insertions of an active P-element. After 10 generations at a large population size, we detected strong selection against P-element insertions. The exception were P-element insertions in genomic regions for which a strong insertion bias has been proposed (2-4). Because recurrent P-element insertions cannot explain this pattern, we conclude that purifying selection, with variable strength along the chromosomes, is the major determinant of the genomic distribution of P-elements. Genomic regions with relaxed purifying selection against P-element insertions exhibit normal levels of purifying selection against base substitutions. This suggests that different types of purifying selection operate on base substitutions and P-element insertions. Our results highlight the power of experimental evolution to understand basic evolutionary processes, which are difficult to infer from patterns of natural variation alone.

Comparative analysis of repeat content in plant genomes, large and small.

The DNA Features pipeline is the analysis pipeline at EMBL-EBI that annotates repeat elements, including transposable elements. With Ensembl's goal to stay at the cutting edge of genome annotation, we proved that this pipeline needed an update. We then created a new analysis that allowed the Ensembl database to store the repeat classification from the PGSB repeat classification (Recat). This new dataset was then fetched using Perl scripts and used to prove that the pipeline modification induced a gain in sensitivity. Finally, we performed a comparative analysis of transposable element distribution in all plant species available, raising new questions about transposable elements in certain branches of the taxonomic tree.

Characterization of piggery farm waste-borne bacterial transposable elements associated with antimicrobial resistance phenotypes

Even though there is a link between antibiotic resistance and the presence of transposable elements few research has looked at the prevalence and distribution of transposable elements/ integrons in piggery farm samples. Present study identified the presence of six transposable elements namely Tn6763 (Accession number: OQ565300), Tn6764, (Accession number: OQ565299), Tn6765 (Accession number: OQ409902), Tn2003 (Accession number: OQ503494), Tn6072 (Accession number: OQ565298) and Tn6020 (Accession number: OQ503493) in piggery farm waste from India which are belongs to Enterobacteriaceae family. In a conjugative experiment, Klebsiella isolates carrying Tn6020 having the resistant phenotypes for nalidixic acid was used as donor cells while Escherichia coli DH5α Cells carrying chloramphenicol resistant plasmid was employed as recipient cells. Transconjugant bacterial colonies were shown to carry the Tn6020 transposable elements with both nalidixic acid (donor cell origin) and chloramphenicol (recipient cell origin) resistant antibiotic phenotypes. Given the presence of transposable elements in 21.4% of resistant Enterobacteriaceae strains, preventative measures are vital for avoiding the spread of mobile genetic resistance determinants in the piggery sector and to monitor their emergence.

Perspective for Studying the Relationship of miRNAs with Transposable Elements.

Transposable elements are important sources of miRNA, long non-coding RNAs genes, and their targets in the composition of protein-coding genes in plants and animals. Therefore, the detection of expression levels of specific non-coding RNAs in various tissues and cells in normal and pathological conditions may indicate a programmed pattern of transposable elements' activation. This reflects the species-specific composition and distribution of transposable elements in genomes, which underlie gene regulation in every cell division, including during aging. TEs' expression is also regulated by epigenetic factors (DNA methylation, histone modifications), SIRT6, cytidine deaminases APOBEC3, APOBEC1, and other catalytic proteins, such as ERCC, TREX1, RB1, HELLS, and MEGP2. In evolution, protein-coding genes and their regulatory elements are derived from transposons. As part of non-coding regions and introns of genes, they are sensors for transcriptional and post-transcriptional control of expression, using miRNAs and long non-coding RNAs, that arose from transposable elements in evolution. Methods (Orbld, ncRNAclassifier) and databases have been created for determining the occurrence of miRNAs from transposable elements in plants (PlanTE-MIR DB, PlaNC-TE), which can be used to design epigenetic gene networks in ontogenesis. Based on the data accumulated in the scientific literature, the presence of 467 transposon-derived miRNA genes in the human genome has been reliably established. It was proposed to create an updated and controlled online bioinformatics database of miRNAs derived from transposable elements in healthy individuals, as well as expression changes of these miRNAs during aging and various diseases, such as cancer and difficult-to-treat diseases. The use of the information obtained can open new horizons in the management of tissue and organ differentiation to aging slow down. In addition, the created database could become the basis for clarifying the mechanisms of pathogenesis of various diseases (imbalance in the activity of transposable elements, reflected in changes in the expression of miRNAs) and designing their targeted therapy using specific miRNAs as targets. This article provides examples of the detection of transposable elements-derived miRNAs involved in the development of specific malignant neoplasms, aging, and idiopathic pulmonary fibrosis.

Transposable elements and their role in aging.

Transposable elements (TEs) are an important part of eukaryotic genomes. The role of somatic transposition in aging, carcinogenesis, and other age-related diseases has been determined. This review discusses the fundamental properties of TEs and their complex interactions with cellular processes, which are crucial for understanding the diverse effects of their activity on the genetics and epigenetics of the organism. The interactions of TEs with recombination, replication, repair, and chromosomal regulation; the ability of TEs to maintain a balance between their own activity and repression, the involvement of TEs in the creation of new or alternative genes, the expression of coding/non-coding RNA, and the role in DNA damage and modification of regulatory networks are reviewed. The contribution of the derepressed TEs to age-dependent effects in individual cells/tissues in different organisms was assessed. Conflicting information about TE activity under stress as well as theories of aging mechanisms related to TEs is discussed. On the one hand, transposition activity in response to stressors can lead to organisms acquiring adaptive innovations of great importance for evolution at the population level. On the other hand, the TE expression can cause decreased longevity and stress tolerance at the individual level. The specific features of TE effects on aging processes in germline and soma and the ways of their regulation in cells are highlighted. Recent results considering somatic mutations in normal human and animal tissues are indicated, with the emphasis on their possible functional consequences. In the context of aging, the correlation between somatic TE activation and age-related changes in the number of proteins required for heterochromatin maintenance and longevity regulation was analyzed. One of the original features of this review is a discussion of not only effects based on the TEs insertions and the associated consequences for the germline cell dynamics and somatic genome, but also the differences between transposon- and retrotransposon-mediated structural genome changes and possible phenotypic characteristics associated with aging and various age-related pathologies. Based on the analysis of published data, a hypothesis about the influence of the species-specific features of number, composition, and distribution of TEs on aging dynamics of different animal genomes was formulated.

Insight into weevil biology from a reference quality genome of the boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae).

The boll weevil, Anthonomus grandis grandis Boheman, is one of the most historically impactful insects due to its near destruction of the US cotton industry in the early 20th century. Contemporary efforts to manage this insect primarily use pheromone baited traps for detection and organophosphate insecticides for control, but this strategy is not sustainable due to financial and environmental costs. We present a high-quality boll weevil genome assembly, consisting of 306 scaffolds with approximately 24,000 annotated genes, as a first step in the identification of gene targets for novel pest control. Gene content and transposable element distribution are similar to those found in other Curculionidae genomes; however, this is the most contiguous and only assembly reported to date for a member in the species-rich genus Anthonomus. Transcriptome profiles across larval, pupal, and adult life stages led to identification of several genes and gene families that could present targets for novel control strategies.

Transposable Elements in Bats Show Differential Accumulation Patterns Determined by Class and Functionality

Bat genomes are characterized by a diverse transposable element (TE) repertoire. In particular, the genomes of members of the family Vespertilionidae contain both active retrotransposons and active DNA transposons. Each TE type is characterized by a distinct pattern of accumulation over the past ~40 million years. Each also exhibits its own target site preferences (sometimes shared with other TEs) that impact where they are likely to insert when mobilizing. Therefore, bats provide a great resource for understanding the diversity of TE insertion patterns. To gain insight into how these diverse TEs impact genome structure, we performed comparative spatial analyses between different TE classes and genomic features, including genic regions and CpG islands. Our results showed a depletion of all TEs in the coding sequence and revealed patterns of species- and element-specific attraction in the transcript. Trends of attraction in the distance tests also suggested significant TE activity in regions adjacent to genes. In particular, the enrichment of small, non-autonomous TE insertions in introns and near coding regions supports the hypothesis that the genomic distribution of TEs is the product of a balance of the TE insertion preference in open chromatin regions and the purifying selection against TEs within genes.

Development and application of universal ND-FISH probes for detecting P-genome chromosomes based on Agropyron cristatum transposable elements.

Fluorescence in situ hybridization (FISH) is a basic tool that is widely used in cytogenetic research. The detection efficiency of conventional FISH is limited due to its time-consuming nature. Oligonucleotide (oligo) probes with fluorescent labels have been applied in non-denaturing FISH (ND-FISH) assays, which greatly streamline experimental processes and save costs and time. Agropyron cristatum, which contains one basic genome, "P," is a vital wild relative for wheat improvement. However, oligo probes for detecting P-genome chromosomes based on ND-FISH assays have not been reported. In this study, according to the distribution of transposable elements (TEs) in Triticeae genomes, 94 oligo probes were designed based on three types of A. cristatum sequences. ND-FISH validation showed that 12 single oligo probes generated a stable and obvious hybridization signal on whole P chromosomes in the wheat background. To improve signal intensity, mixed probes (Oligo-pAc) were prepared by using the 12 successful probes and validated in the diploid accession A. cristatum Z1842, a small segmental translocation line and six allopolyploid wild relatives containing the P genome. The signals of Oligo-pAc covered the entire chromosomes of A. cristatum and were more intense than those of single probes. The results indicate that Oligo-pAc can replace conventional genomic in situ hybridization (GISH) probes to identify P chromosomes or segments in non-P-genome backgrounds. Finally, we provide a rapid and efficient method specifically for detecting P chromosomes in wheat backgrounds by combining the Oligo-pAc probe with the Oligo-pSc119.2-1 and Oligo-pTa535-1 probes, which can replace conventional sequential GISH/FISH assays. Altogether, we developed a set of oligo probes based on the ND-FISH assays to identify P-genome chromosomes, which can promote utilization of A. cristatum in wheat improvement programs.

Transposon- and Genome Dynamics in the Fungal Genus Neurospora: Insights from Nearly Gapless Genome Assemblies

A large portion of nuclear DNA is composed of transposable element (TE) sequences, whose transposition is controlled by diverse host defense strategies in order to maintain genomic integrity. One such strategy is the fungal-specific Repeat-Induced Point mutation (RIP) that hyper-mutates repetitive DNA sequences. While RIP is found across Fungi, it has been shown to vary in efficiency. The filamentous ascomycete Neurospora crassa has been a pioneer in the study of RIP, but data on TEs and RIP from other species in the genus is limited. In this study, we investigated 18 nearly gapless genome assemblies of ten Neurospora species, which diverged from a common ancestor about 7 MYA, to determine and compare genome-wide TE distribution and their associated RIP patterns. Four of these assemblies, generated by PacBio technology, represent new genomic datasets. We showed that the TE contents between 8.7-18.9% covary with genome sizes that range between 37.8-43.9 Mb. Degraded copies of Long Terminal Repeat (LTR) retrotransposons were abundant among the identified TEs, and these are distributed across the genome at varying frequencies. In all investigated Neurospora genomes, TE sequences had signs of numerous C-to-T substitutions, suggesting that RIP occurred in all species, and accordingly, RIP signatures correlated with TE-dense regions in all genomes. In conclusion, essentially gapless genome assemblies allowed us to identify TEs in Neurospora genomes, and reveal that TEs contribute to genome size variation in this group. Our study suggests that TEs and RIP are highly correlated in each examined Neurospora species, and hence, the pattern of interaction is conserved over the investigated evolutionary timescale. Finally, with our results, we verify that RIP signatures can be used to facilitate the identification of TE-rich region in the genome. The comprehensive genomic dataset of Neurospora is a rich resource for further in-depth analyses of fungal genomes by the community.

Controlling for Variable Transposition Rate with an Age-Adjusted Site Frequency Spectrum.

Recognition of the important role of transposable elements (TEs) in eukaryotic genomes quickly led to a burgeoning literature modeling and estimating the effects of selection on TEs. Much of the empirical work on selection has focused on analyzing the site frequency spectrum (SFS) of TEs. But TE evolution differs from standard models in a number of ways that can impact the power and interpretation of the SFS. For example, rather than mutating under a clock-like model, transposition often occurs in bursts which can inflate particular frequency categories compared with expectations under a standard neutral model. If a TE burst has been recent, the excess of low-frequency polymorphisms can mimic the effect of purifying selection. Here, we investigate how transposition bursts affect the frequency distribution of TEs and the correlation between age and allele frequency. Using information on the TE age distribution, we propose an age-adjusted SFS to compare TEs and neutral polymorphisms to more effectively evaluate whether TEs are under selective constraints. We show that our approach can minimize instances of false inference of selective constraint, remains robust to simple demographic changes, and allows for a correct identification of even weak selection affecting TEs which experienced a transposition burst. The results presented here will help researchers working on TEs to more reliably identify the effects of selection on TEs without having to rely on the assumption of a constant transposition rate.

Zebrafish transposable elements show extensive diversification in age, genomic distribution, and developmental expression.

There is considerable interest in understanding the effect of transposable elements (TEs) on embryonic development. Studies in humans and mice are limited by the difficulty of working with mammalian embryos and by the relative scarcity of active TEs in these organisms. The zebrafish is an outstanding model for the study of vertebrate development, and over half of its genome consists of diverse TEs. However, zebrafish TEs remain poorly characterized. Here we describe the demography and genomic distribution of zebrafish TEs and their expression throughout embryogenesis using bulk and single-cell RNA sequencing data. These results reveal a highly dynamic genomic ecosystem comprising nearly 2000 distinct TE families, which vary in copy number by four orders of magnitude and span a wide range of ages. Longer retroelements tend to be retained in intergenic regions, whereas short interspersed nuclear elements (SINEs) and DNA transposons are more frequently found nearby or within genes. Locus-specific mapping of TE expression reveals extensive TE transcription during development. Although two-thirds of TE transcripts are likely driven by nearby gene promoters, we still observe stage- and tissue-specific expression patterns in self-regulated TEs. Long terminal repeat (LTR) retroelements are most transcriptionally active immediately following zygotic genome activation, whereas DNA transposons are enriched among transcripts expressed in later stages of development. Single-cell analysis reveals several endogenous retroviruses expressed in specific somatic cell lineages. Overall, our study provides a valuable resource for using zebrafish as a model to study the impact of TEs on vertebrate development.

Transposable Elements in the Genome of the Lichen-Forming Fungus Umbilicaria pustulata and Their Distribution in Different Climate Zones along Elevation.

Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Drift and natural selection are important forces shaping TE distribution and accumulation. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host's ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold temperate. We trace the occurrence of the newly identified TEs in populations along three elevation gradients using a Pool-Seq approach to identify TE insertions of potential adaptive significance. We found that TEs cover 21.26% of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. We identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function. This pioneering study of the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution.

Comparative analysis of transposable elements provides insights into genome evolution in the genus Camelus

BackgroundTransposable elements (TEs) are common features in eukaryotic genomes that are known to affect genome evolution critically and to play roles in gene regulation. Vertebrate genomes are dominated by TEs, which can reach copy numbers in the hundreds of thousands. To date, details regarding the presence and characteristics of TEs in camelid genomes have not been made available.ResultsWe conducted a genome-wide comparative analysis of camelid TEs, focusing on the identification of TEs and elucidation of transposition histories in four species: Camelus dromedarius, C. bactrianus, C. ferus, and Vicugna pacos. Our TE library was created using both de novo structure-based and homology-based searching strategies (https://github.com/kacst-bioinfo-lab/TE_ideintification_pipeline). Annotation results indicated a similar proportion of each genomes comprising TEs (35–36%). Class I LTR retrotransposons comprised 16–20% of genomes, and mostly consisted of the endogenous retroviruses (ERVs) groups ERVL, ERVL-MaLR, ERV_classI, and ERV_classII. Non-LTR elements comprised about 12% of genomes and consisted of SINEs (MIRs) and the LINE superfamilies LINE1, LINE2, L3/CR1, and RTE clades. Least represented were the Class II DNA transposons (2%), consisting of hAT-Charlie, TcMar-Tigger, and Helitron elements and comprising about 1–2% of each genome.ConclusionsThe findings of the present study revealed that the distribution of transposable elements across camelid genomes is approximately similar. This investigation presents a characterization of TE content in four camelid to contribute to developing a better understanding of camelid genome architecture and evolution.

Comparative genomic and transcriptomic analyses of transposable elements in polychaetous annelids highlight LTR retrotransposon diversity and evolution

BackgroundWith the expansion of high throughput sequencing, we now have access to a larger number of genome-wide studies analyzing the Transposable elements (TEs) composition in a wide variety of organisms. However, genomic analyses often remain too limited in number and diversity of species investigated to study in depth the dynamics and evolutionary success of the different types of TEs among metazoans. Therefore, we chose to investigate the use of transcriptomes to describe the diversity of TEs in phylogenetically related species by conducting the first comparative analysis of TEs in two groups of polychaetes and evaluate the diversity of TEs that might impact genomic evolution as a result of their mobility.ResultsWe present a detailed analysis of TEs distribution in transcriptomes extracted from 15 polychaetes depending on the number of reads used during assembly, and also compare these results with additional TE scans on associated low-coverage genomes. We then characterized the clades defined by 1021 LTR-retrotransposon families identified in 26 species. Clade richness was highly dependent on the considered superfamily. Copia elements appear rare and are equally distributed in only three clades, GalEa, Hydra and CoMol. Among the eight BEL/Pao clades identified in annelids, two small clades within the Sailor lineage are new for science. We characterized 17 Gypsy clades of which only 4 are new; the C-clade largely dominates with a quarter of the families. Finally, all species also expressed for the majority two distinct transcripts encoding PIWI proteins, known to be involved in control of TEs mobilities.ConclusionsThis study shows that the use of transcriptomes assembled from 40 million reads was sufficient to access to the diversity and proportion of the transposable elements compared to those obtained by low coverage sequencing. Among LTR-retrotransposons Gypsy elements were unequivocally dominant but results suggest that the number of Gypsy clades, although high, may be more limited than previously thought in metazoans. For BEL/Pao elements, the organization of clades within the Sailor lineage appears more difficult to establish clearly. The Copia elements remain rare and result from the evolutionary consistent success of the same three clades.

The pattern of DNA methylation alteration, and its association with the changes of gene expression and alternative splicing during phosphate starvation in tomato

SUMMARYDNA methylation is changed and associates with gene expression alterations in plant response to phosphate starvation (Pi–), a common stress that impacts plant growth and productivity. However, in the horticultural model species Solanum lycopersicum (tomato), the dynamics of DNA methylation and its relationship with changes in gene transcription and alternative splicing (AS) under Pi– are unknown. Here, we performed integrative methylome and transcriptome analyses of tomato seedlings under Pi‐deficient and ‐sufficient conditions. We found Pi– caused a slight increase in the overall methylation level, with millions of differentially methylated cytosines (DmCs) and a few hundred differentially methylated regions (DMRs). We also identified thousands of differentially expressed (DE) and differential AS (DAS) genes induced by Pi–, and found that DmCs were more abundant in non‐expressed genes than in DE or DAS genes. Moreover, DNA methylation alterations weakly correlated with transcription changes but not with DAS events, and hyper‐CHH‐DMRs overlapping with transposable elements (TEs) were enriched in a subset of Pi starvation response (PSR) genes. We propose that changes in DNA methylation may be associated with the differential expression of some PSR genes, but that most of these changes probably control the expression of nearby TEs, rather than directly affecting the transcription or AS of PSR genes. Besides, the pattern of methylation changes upon Pi– may largely be shaped by TE distributions. Together, our study provides comprehensive insights into the association of DNA methylation with gene transcription and AS under Pi– in tomato and may contribute to unveiling novel roles of epigenetic mechanisms in plant stress response.

Chromosome-scale assembly reveals asymmetric paleo-subgenome evolution and targets for the acceleration of fungal resistance breeding in the nut crop, pecan

Pecan (Carya illinoinensis) is a tree nut crop of worldwide economic importance that is rich in health-promoting factors. However, pecan production and nut quality are greatly challenged by environmental stresses such as the outbreak of severe fungal diseases. Here, we report a high-quality, chromosome-scale genome assembly of the controlled-cross pecan cultivar ‘Pawnee’ constructed by integrating Nanopore sequencing and Hi-C technologies. Phylogenetic and evolutionary analyses reveal two whole-genome duplication (WGD) events and two paleo-subgenomes in pecan and walnut. Time estimates suggest that the recent WGD event and considerable genome rearrangements in pecan and walnut account for expansions in genome size and chromosome number after the divergence from bayberry. The two paleo-subgenomes differ in size and protein-coding gene sets. They exhibit uneven ancient gene loss, asymmetrical distribution of transposable elements (especially LTR/Copia and LTR/Gypsy), and expansions in transcription factor families (such as the extreme pecan-specific expansion in the far-red impaired response 1 family), which are likely to reflect the long evolutionary history of species in the Juglandaceae. A whole-genome scan of resequencing data from 86 pecan scab-associated core accessions identified 47 chromosome regions containing 185 putative candidate genes. Significant changes were detected in the expression of candidate genes associated with the chitin response pathway under chitin treatment in the scab-resistant and scab-susceptible cultivars ‘Excell’ and ‘Pawnee’. These findings enable us to identify key genes that may be important susceptibility factors for fungal diseases in pecan. The high-quality sequences are valuable resources for pecan breeders and will provide a foundation for the production and quality improvement of tree nut crops.

Transposable elements and introgression introduce genetic variation in the invasive ant Cardiocondyla obscurior

Introduced populations of invasive organisms have to cope with novel environmental challenges, while having reduced genetic variation caused by founder effects. The mechanisms associated with this "genetic paradox of invasive species" has received considerable attention, yet few studies have examined the genomic architecture of invasive species. Populations of the heart node ant Cardiocondyla obscurior belong to two distinct lineages, a New World lineage so far only found in Latin America and a more globally distributed Old World lineage. In the present study, we use population genomic approaches to compare populations of the two lineages with apparent divergent invasive potential. We find that the strong genetic differentiation of the two lineages began at least 40,000 generations ago and that activity of transposable elements (TEs) has contributed significantly to the divergence of both lineages, possibly linked to the very unusual genomic distribution of TEs in this species. Furthermore, we show that introgression from the Old World lineage is a dominant source of genetic diversity in the New World lineage, despite the lineages' strong genetic differentiation. Our study uncovers mechanisms underlying novel genetic variation in introduced populations of C.obscurior that could contribute to the species' adaptive potential.