Transposable element influences on gene expression in plants

De novo Whole-Genome Assembly of the 10-Gigabase Fokienia Hodginsii Genome to Reveal Differential Epigenetic Events Between Callus and Xylem

Fokienia hodginsii (F. hodginsii), belonging to the genus Fokienia of the Cupressaceae. F. hodginsii has significant application value due to its wood properties and great research value in evolutionary studies as a gymnosperm. However, the genome of F. hodginsii remains unknown due to the large size of gymnosperms genome. Pacific Bioscience sequencing, Hi-C mapping, whole-genome Bisulfite Sequencing (BS-Seq), long-read isoform sequencing (Iso-Seq), direct RNA sequencing (DRS), quantitative proteomics, and metabonomics analysis are employed to facilitate genome assembly, gene annotation, and investigation into epigenetic mechanisms. In this study, the 10G F. hodginsii genome is assembled into 11 chromosomes. Furthermore, 50 521 protein-coding genes are annotated and determined that 65% of F. hodginsii genome comprises repetitive sequences. It is discovered that transposable element (TE)-including introns is associated with higher expression. The DNA methylome of F. hodginsii reveals that xylem has a higher DNA methylation level compared to callus. Moreover, DRS reveals the significant alterations in RNA full-length ratio, which potentially associated with poly(A) length (PAL) and alternative polyadenylation (APA). Finally, the morphology measurement and metabonomics analysis revealed the difference of 14 cultivars. In summary, the genomes and epigenetics datasets provide a molecular basis for callus formation in the gymnosperm family.

DANTE and DANTE_LTR: lineage-centric annotation pipelines for long terminal repeat retrotransposons in plant genomes

Long terminal repeat (LTR) retrotransposons constitute a predominant class of repetitive DNA elements in most plant genomes. With the increasing number of sequenced plant genomes, there is an ongoing demand for computational tools facilitating efficient annotation and classification of LTR retrotransposons in plant genome assemblies. Herein, we introduce DANTE, a computational pipeline for Domain-based ANnotation of Transposable Elements, designed for sensitive detection of these elements via their conserved protein domain sequences. The identified protein domains are subsequently inputted into the DANTE_LTR pipeline to annotate complete element sequences by detecting their structural features, such as LTRs, in adjacent genomic regions. Leveraging domain sequences allows for precise classification of elements into phylogenetic lineages, offering a more granular annotation compared with coarser conventional superfamily-based classification methods. The efficiency and accuracy of this approach were evidenced via annotation of LTR retrotransposons in 93 plant genomes. Results were benchmarked against several established pipelines, showing that DANTE_LTR is capable of identifying significantly more intact LTR retrotransposons. DANTE and DANTE_LTR are provided as user-friendly Galaxy tools accessible via a public server (https://repeatexplorer-elixir.cerit-sc.cz), installable on local Galaxy instances from the Galaxy tool shed or executable from the command line.

Unlocking plant genetics with telomere-to-telomere genome assemblies

Contiguous genome sequence assemblies will help us to realize the full potential of crop translational genomics. Recent advances in sequencing technologies, especially long-read sequencing strategies, have made it possible to construct gapless telomere-to-telomere (T2T) assemblies, thus offering novel insights into genome organization and function. Plant genomes pose unique challenges, such as a continuum of ancient to recent polyploidy and abundant highly similar and long repetitive elements. Owing to progress in sequencing approaches, for most crop plants, chromosome-scale reference genome assemblies are available, but T2T assembly construction remains challenging. Here we describe methods for haplotype-resolved, gapless T2T assembly construction in plants, including various crop species. We outline the impact of T2T assemblies in elucidating the roles of repetitive elements in gene regulation, as well as in pangenomics, functional genomics, genome-assisted breeding and targeted genome manipulation. In conjunction with sequence-enriched germplasm repositories, T2T assemblies thus hold great promise for basic and applied plant sciences.

Epigenetic and transcriptional consequences in the endosperm of chemically induced transposon mobilization in Arabidopsis

Genomic imprinting, an epigenetic phenomenon leading to parent-of-origin-specific gene expression, has independently evolved in the endosperm of flowering plants and the placenta of mammals-tissues crucial for nurturing embryos. While transposable elements (TEs) frequently colocalize with imprinted genes and are implicated in imprinting establishment, direct investigations of the impact of de novo TE transposition on genomic imprinting remain scarce. In this study, we explored the effects of chemically induced transposition of the Copia element ONSEN on genomic imprinting in Arabidopsis thaliana. Through the combination of chemical TE mobilization and doubled haploid induction, we generated a line with 40 new ONSEN copies. Our findings reveal a preferential targeting of maternally expressed genes (MEGs) for transposition, aligning with the colocalization of H2A.Z and H3K27me3 in MEGs-both previously identified as promoters of ONSEN insertions. Additionally, we demonstrate that chemically-induced DNA hypomethylation induces global transcriptional deregulation in the endosperm, leading to the breakdown of MEG imprinting. This study provides insights into the consequences of chemically induced TE remobilization in the endosperm, revealing that chemically-induced epigenome changes can have long-term consequences on imprinted gene expression.

Evolutionary History of the DD41D Family of Tc1/Mariner Transposons in Two Mayetiola Species

Tc1/mariner elements are ubiquitous in eukaryotic genomes including insects. They are diverse and divided into families and sub-families. The DD34D family including mauritiana and irritans subfamilies have already been identified in two closely related species of Cecidomyiids M. destructor and M. hordei. In the current study the de novo and similarity-based methods allowed the identification for the first time of seven consensuses in M. destructor and two consensuses in M. hordei belonging to DD41D family whereas the in vitro method allowed the amplification of two and three elements in these two species respectively. Most of identified elements accumulated different mutations and long deletions spanning the N-terminal region of the transposase. Phylogenetic analyses showed that the DD41D elements were clustered in two groups belonging to rosa and Long-TIR subfamilies. The age estimation of the last transposition events of the identified Tc1/mariner elements in M. destructor showed different evolutionary histories. Indeed, irritans elements have oscillated between periods of silencing and reappearance while rosa and mauritiana elements have shown regular activity with large recent bursts. The study of insertion sites showed that they are mostly intronic and that some recently transposed elements occurred in genes linked to putative DNA-binding domains and enzymes involved in metabolic chains. Thus, this study gave evidence of the existence of DD41D family in two Mayetiola species and an insight on their evolutionary history.

Ectopic enhancer-enhancer interactions as causal forces driving RNA-directed DNA methylation in gene regulatory regions

Cis-regulatory elements (CREs) are integral to the spatiotemporal and quantitative expression dynamics of target genes, thus directly influencing phenotypic variation and evolution. However, many of these CREs become highly susceptible to transcriptional silencing when in a transgenic state, particularly when organised as tandem repeats. We investigated the mechanism of this phenomenon and found that three of the six selected flower-specific CREs were prone to transcriptional silencing when in a transgenic context. We determined that this silencing was caused by the ectopic expression of non-coding RNAs (ncRNAs), which were processed into 24-nt small interfering RNAs (siRNAs) that drove RNA-directed DNA methylation (RdDM). Detailed analyses revealed that aberrant ncRNA transcription within the AGAMOUS enhancer (AGe) in a transgenic context was significantly enhanced by an adjacent CaMV35S enhancer (35Se). This particular enhancer is known to mis-activate the regulatory activities of various CREs, including the AGe. Furthermore, an insertion of 35Se approximately 3.5 kb upstream of the AGe in its genomic locus also resulted in the ectopic induction of ncRNA/siRNA production and de novo methylation specifically in the AGe, but not other regions, as well as the production of mutant flowers. This confirmed that interactions between the 35Se and AGe can induce RdDM activity in both genomic and transgenic states. These findings highlight a novel epigenetic role for CRE-CRE interactions in plants, shedding light on the underlying forces driving hypermethylation in transgenes, duplicate genes/enhancers, and repetitive transposons, in which interactions between CREs are inevitable.

Identification of Highly Repetitive Enhancers with Long-range Regulation Potential in Barley via STARR-seq

Enhancers are DNA sequences that can strengthen transcription initiation. However, the global identification of plant enhancers is complicated due to uncertainty in the distance and orientation of enhancers, especially in species with large genomes. In this study, we performed self-transcribing active regulatory region sequencing (STARR-seq) for the first time to identify enhancers across the barley genome. A total of 7323 enhancers were successfully identified, and among 45 randomly selected enhancers, over 75% were effective as validated by a dual-luciferase reporter assay system in the lower epidermis of tobacco leaves. Interestingly, up to 53.5% of the barley enhancers were repetitive sequences, especially transposable elements (TEs), thus reinforcing the vital role of repetitive enhancers in gene expression. Both the common active mark H3K4me3 and repressive mark H3K27me3 were abundant among the barley STARR-seq enhancers. In addition, the functional range of barley STARR-seq enhancers seemed much broader than that of rice or maize and extended to ±100 kb of the gene body, and this finding was consistent with the high expression levels of genes in the genome. This study specifically depicts the unique features of barley enhancers and provides available barley enhancers for further utilization.

ZmWRKY30 modulates drought tolerance in maize by influencing myo-inositol and reactive oxygen species homeostasis

Maize (Zea mays L.) is an important food crop with a wide range of uses in both industry and agriculture. Drought stress during its growth cycle can greatly reduce maize crop yield and quality. However, the molecular mechanisms underlying maize responses to drought stress remain unclear. In this work, a WRKY transcription factor-encoding gene, ZmWRKY30, from drought-treated maize leaves was screened out and characterized. ZmWRKY30 gene expression was induced by dehydration treatments. The ZmWRKY30 protein localized to the nucleus and displayed transactivation activity in yeast. Compared with wild-type (WT) plants, Arabidopsis lines overexpressing ZmWRKY30 exhibited a significantly enhanced drought stress tolerance, as evidenced by the improved survival rate, increased antioxidant enzyme activity by superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT), elevated proline content, and reduced lipid peroxidation recorded after drought stress treatment. In contrast, the mutator (Mu)-interrupted ZmWRKY30 homozygous mutant (zmwrky30) was more sensitive to drought stress than its null segregant (NS), characterized by the decreased survival rate, reduced antioxidant enzyme activity (SOD, POD, and CAT) and proline content, as well as increased malondialdehyde accumulation. RNA-Seq analysis further revealed that, under drought conditions, the knockout of the ZmWRKY30 gene in maize affected the expression of genes involved in reactive oxygen species (ROS), proline, and myo-inositol metabolism. Meanwhile, the zmwrky30 mutant exhibited significant downregulation of myo-inositol content in leaves under drought stress. Combined, our results suggest that ZmWRKY30 positively regulates maize responses to water scarcity. This work provides potential target genes for the breeding of drought-tolerant maize.

Utilizing a comparative approach to assess genome evolution during diploidization in Artemisia tridentata, a keystone species of western North America

PREMISE

Polyploidization is often followed by diploidization. Diploidization is generally studied using synthetic polyploid lines and/or crop plants, but rarely using extant diploids or nonmodel plants such as Artemisia tridentata. This threatened western North American keystone species has a large genome compared to congeneric Artemisia species; dominated by diploid and tetraploid cytotypes, with multiple origins of tetraploids with genome size reduction.

METHODS

The genome of an A. tridentata sample was resequenced to study genome evolution and compared to that of A. annua, a diploid congener. Three diploid genomes of A. tridentata were compared to test for multiple diploidization events.

RESULTS

The A. tridentata genome had many chromosomal rearrangements relative to that of A. annua, while large-scale synteny of A. tridentata chromosome 3 and A. annua chromosome 4 was conserved. The three A. tridentata genomes had similar sizes (4.19-4.2 Gbp), heterozygosity (2.24-2.25%), and sequence (98.73-99.15% similarity) across scaffolds, and in k-mer analyses, similar patterns of diploid heterozygous k-mers (AB = 41%, 47%, and 47%), triploid heterozygous k-mers (AAB = 18-21%), and tetraploid k-mers (AABB = 13-17%). Biallelic SNPs were evenly distributed across scaffolds for all individuals. Comparisons of transposable element (TE) content revealed differential enrichment of TE clades.

CONCLUSIONS

Our findings suggest population-level TE differentiation after a shared polyploidization-to-diploidization event(s) and exemplify the complex processes of genome evolution. This research approached provides new resources for exploration of abiotic stress response, especially the roles of TEs in response pathways.

Evolutionary insights into the organization of chromatin structure and landscape of transcriptional regulation in plants

Development of complex traits necessitates the functioning and coordination of intricate regulatory networks involving multiple genes. Understanding 3D chromatin structure can facilitate insight into the regulation of gene expression by regulatory elements. This potential, of visualizing the role of chromatin organization in the evolution and function of regulatory elements, remains largely unexplored. Here, we describe new perspectives that arise from the dual considerations of sequence variation of regulatory elements and chromatin structure, with a special focus on whole-genome doubling or polyploidy. We underscore the significance of hierarchical chromatin organization in gene regulation during evolution. In addition, we describe strategies for exploring chromatin organization in future investigations of regulatory evolution in plants, enabling insights into the evolutionary influence of regulatory elements on gene expression and, hence, phenotypes.

ARGONAUTE1-binding Tudor domain proteins function in small interfering RNA production for RNA-directed DNA methylation

Transposable elements (TEs) contribute to plant evolution, development, and adaptation to environmental changes, but the regulatory mechanisms are largely unknown. RNA-directed DNA methylation (RdDM) is 1 TE regulatory mechanism in plants. Here, we identified that novel ARGONAUTE 1 (AGO1)-binding Tudor domain proteins Precocious dissociation of sisters C/E (PDS5C/E) are involved in 24-nt siRNA production to establish RdDM on TEs in Arabidopsis thaliana. PDS5 family proteins are subunits of the eukaryote-conserved cohesin complex. However, the double mutant lacking angiosperm-specific subfamily PDS5C and PDS5E (pds5c/e) exhibited different developmental phenotypes and transcriptome compared with those of the double mutant lacking eukaryote-conserved subfamily PDS5A and PDS5B (pds5a/b), suggesting that the angiosperm-specific PDS5C/E subfamily has a unique function in angiosperm plants. Proteome and imaging analyses revealed that PDS5C/E interact with AGO1. The pds5c/e double mutant had defects in 24-nt siRNA accumulation and CHH DNA methylation on TEs. In addition, some lncRNAs that accumulated in the pds5c/e mutant were targeted by AGO1-loading 21-nt miRNAs and 21-nt siRNAs. These results indicate that PDS5C/E and AGO1 participate in 24-nt siRNA production for RdDM in the cytoplasm. These findings indicate that angiosperm plants evolved a new regulator, the PDS5C/E subfamily, to control the increase in TEs during angiosperm evolution.

Molecular Mimicry of Transposable Elements in Plants

Transposable elements (TEs) are mobile DNA elements that are particularly abundant in the plant genomes. They have long been considered as junk DNA; however, a growing body of evidence suggests that TE insertions promote genetic diversity that is essential for the adaptive evolution of a species. Thus far, studies have mainly investigated the cis-acting regulatory roles of TEs generated by their insertions nearby or within the host genes. However, the trans-acting effects of TE-derived RNA and DNA remained obscure to date. TEs contain various regulatory elements within their sequences that can accommodate the binding of specific RNAs and proteins. Recently, it was suggested that some of these cellular regulators are shared between TEs and the host genes, and the competition for the common host factors underlies the fine-tuned developmental reprogramming. In this review, we will highlight and discuss the latest discoveries on the biological functions of plant TEs, with a particular focus on their competitive binding with specific developmental regulators.

Transposable elements contribute to the establishment of the glycine shuttle in Brassicaceae species

C3 -C4 intermediate photosynthesis has evolved at least five times convergently in the Brassicaceae, despite this family lacking bona fide C4 species. The establishment of this carbon concentrating mechanism is known to require a complex suite of ultrastructural modifications, as well as changes in spatial expression patterns, which are both thought to be underpinned by a reconfiguration of existing gene-regulatory networks. However, to date, the mechanisms which underpin the reconfiguration of these gene networks are largely unknown. In this study, we used a pan-genomic association approach to identify genomic features that could confer differential gene expression towards the C3 -C4 intermediate state by analysing eight C3 species and seven C3 -C4 species from five independent origins in the Brassicaceae. We found a strong correlation between transposable element (TE) insertions in cis-regulatory regions and C3 -C4 intermediacy. Specifically, our study revealed 113 gene models in which the presence of a TE within a gene correlates with C3 -C4 intermediate photosynthesis. In this set, genes involved in the photorespiratory glycine shuttle are enriched, including the glycine decarboxylase P-protein whose expression domain undergoes a spatial shift during the transition to C3 -C4 photosynthesis. When further interrogating this gene, we discovered independent TE insertions in its upstream region which we conclude to be responsible for causing the spatial shift in GLDP1 gene expression. Our findings hint at a pivotal role of TEs in the evolution of C3 -C4 intermediacy, especially in mediating differential spatial gene expression.

The human PTGR1 gene expression is controlled by TE-derived Z-DNA forming sequence cooperating with miR-6867-5p

Z-DNA, a well-known non-canonical form of DNA involved in gene regulation, is often found in gene promoters. Transposable elements (TEs), which make up 45% of the human genome, can move from one location to another within the genome. TEs play various biological roles in host organisms, and like Z-DNA, can influence transcriptional regulation near promoter regions. MicroRNAs (miRNAs) are a class of small non-coding RNA molecules that play a critical role in the regulation of gene expression. Although TEs can generate Z-DNA and miRNAs can bind to Z-DNA, how these factors affect gene transcription has yet to be elucidated. Here, we identified potential Z-DNA forming sequence (ZFS), including TE-derived ZFS, in the promoter of prostaglandin reductase 1 (PTGR1) by data analysis. The transcriptional activity of these ZFS in PTGR1 was confirmed using dual-luciferase reporter assays. In addition, we discovered a novel ZFS-binding miRNA (miR-6867-5p) that suppressed PTGR1 expression by targeting to ZFS. In conclusion, these findings suggest that ZFS, including TE-derived ZFS, can regulate PTGR1 gene expression and that miR-6867-5p can suppress PTGR1 by interacting with ZFS.

Plant histone variants at the nexus of chromatin readouts, stress and development

Histones are crucial proteins that are involved in packaging the DNA as condensed chromatin inside the eukaryotic cell nucleus. Rather than being static packaging units, these molecules undergo drastic variations spatially and temporally to facilitate accessibility of DNA to replication, transcription as well as wide range of gene regulatory machineries. In addition, incorporation of paralogous variants of canonical histones in the chromatin is ascribed to specific functions. Given the peculiar requirement of plants to rapidly modulate gene expression levels on account of their sessile nature, histones and their variants serve as additional layers of gene regulation. This review summarizes the mechanisms and implications of distribution, modifications and differential incorporation of histones and their variants across plant genomes, and outlines emerging themes.

Quantifying the role of genome size and repeat content in adaptive variation and the architecture of flowering time in Amaranthus tuberculatus

Genome size variation, largely driven by repeat content, is poorly understood within and among populations, limiting our understanding of its significance for adaptation. Here we characterize intraspecific variation in genome size and repeat content across 186 individuals of Amaranthus tuberculatus, a ubiquitous native weed that shows flowering time adaptation to climate across its range and in response to agriculture. Sequence-based genome size estimates vary by up to 20% across individuals, consistent with the considerable variability in the abundance of transposable elements, unknown repeats, and rDNAs across individuals. The additive effect of this variation has important phenotypic consequences-individuals with more repeats, and thus larger genomes, show slower flowering times and growth rates. However, compared to newly-characterized gene copy number and polygenic nucleotide changes underlying variation in flowering time, we show that genome size is a marginal contributor. Differences in flowering time are reflected by genome size variation across sexes and marginally, habitats, while polygenic variation and a gene copy number variant within the ATP synthesis pathway show consistently stronger environmental clines than genome size. Repeat content nonetheless shows non-neutral distributions across the genome, and across latitudinal and environmental gradients, demonstrating the numerous governing processes that in turn influence quantitative genetic variation for phenotypes key to plant adaptation.

In Silico Verification of Predicted Potential Promoter Sequences in the Rice (Oryza sativa) Genome

The exact identification of promoter sequences remains a serious problem in computational biology, as the promoter prediction algorithms under development continue to produce false-positive results. Therefore, to fully assess the validity of predicted sequences, it is necessary to perform a comprehensive test of their properties, such as the presence of downstream transcribed DNA regions behind them, or chromatin accessibility for transcription factor binding. In this paper, we examined the promoter sequences of chromosome 1 of the rice Oryza sativa genome from the Database of Potential Promoter Sequences predicted using a mathematical algorithm based on the derivation and calculation of statistically significant promoter classes. In this paper TATA motifs and cis-regulatory elements were identified in the predicted promoter sequences. We also verified the presence of potential transcription start sites near the predicted promoters by analyzing CAGE-seq data. We searched for unannotated transcripts behind the predicted sequences by de novo assembling transcripts from RNA-seq data. We also examined chromatin accessibility in the region of the predicted promoters by analyzing ATAC-seq data. As a result of this work, we identified the predicted sequences that are most likely to be promoters for further experimental validation in an in vivo or in vitro system.

Advances in cis-element- and natural variation-mediated transcriptional regulation and applications in gene editing of major crops

Transcriptional regulation is crucial to control of gene expression. Both spatio-temporal expression patterns and expression levels of genes are determined by the interaction between cis-acting elements and trans-acting factors. Numerous studies have focused on the trans-acting factors that mediate transcriptional regulatory networks. However, cis-acting elements, such as enhancers, silencers, transposons, and natural variations in the genome, are also vital for gene expression regulation and could be utilized by clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9)-mediated gene editing to improve crop quality and yield. In this review, we discuss current understanding of cis-element-mediated transcriptional regulation in major crops, including rice (Oryza sativa), wheat (Triticum aestivum), and maize (Zea mays), as well as the latest advancements in gene editing techniques and their applications in crops to highlight prospective strategies for crop breeding.

Biotechnological Approaches for Enhancing Polyhydroxyalkanoates (PHAs) Production: Current and Future Perspectives

The benefits of biotechnology are not limited to genetic engineering, but it also displays its great impact on industrial uses of crops (e.g., biodegradable plastics). Polyhydroxyalkanoates (PHAs) make a diverse class of bio-based and biodegradable polymers naturally synthesized by numerous microorganisms. However, several C3 and C4 plants have also been genetically engineered to produce PHAs. The highest production yield of PHAs was obtained with a well-known C3 plant, Arabidopsis thaliana, upto 40% of the dry weight of the leaf. This review summarizes all biotechnological mechanisms that have been adopted with the goal of increasing PHAs production in bacteria and plant species alike. Moreover, the possibility of using some methods that have not been applied in bioplastic science are discussed with special attention to plants. These include producing PHAs in transgenic hairy roots and cell suspension cultures, making transformed bacteria and plants via transposons, constructing an engineered metabolon, and overexpressing of phaP and the ABC operon concurrently. Taken together, that biotechnology will be highly beneficial for reducing plastic pollution through the implementation of biotechnological strategies is taken for granted.

Young SINEs in pig genomes impact gene regulation, genetic diversity, and complex traits

Transposable elements (TEs) are a major source of genetic polymorphisms and play a role in chromatin architecture, gene regulatory networks, and genomic evolution. However, their functional role in pigs and contributions to complex traits are largely unknown. We created a catalog of TEs (n = 3,087,929) in pigs and found that young SINEs were predominantly silenced by histone modifications, DNA methylation, and decreased accessibility. However, some transcripts from active young SINEs showed high tissue-specificity, as confirmed by analyzing 3570 RNA-seq samples. We also detected 211,067 dimorphic SINEs in 374 individuals, including 340 population-specific ones associated with local adaptation. Mapping these dimorphic SINEs to genome-wide associations of 97 complex traits in pigs, we found 54 candidate genes (e.g., ANK2 and VRTN) that might be mediated by TEs. Our findings highlight the important roles of young SINEs and provide a supplement for genotype-to-phenotype associations and modern breeding in pigs.

Vitis labrusca genome assembly reveals diversification between wild and cultivated grapevine genomes

Wild grapevines are important genetic resources in breeding programs to confer adaptive fitness traits and unique fruit characteristics, but the genetics underlying these traits, and their evolutionary origins, are largely unknown. To determine the factors that contributed to grapevine genome diversification, we performed comprehensive intragenomic and intergenomic analyses with three cultivated European (including the PN40024 reference genome) and two wild North American grapevine genomes, including our newly released Vitis labrusca genome. We found the heterozygosity of the cultivated grapevine genomes was twice as high as the wild grapevine genomes studied. Approximately 30% of V. labrusca and 48% of V. vinifera Chardonnay genes were heterozygous or hemizygous and a considerable number of collinear genes between Chardonnay and V. labrusca had different gene zygosity. Our study revealed evidence that supports gene gain-loss events in parental genomes resulted in the inheritance of hemizygous genes in the Chardonnay genome. Thousands of segmental duplications supplied source material for genome-specific genes, further driving diversification of the genomes studied. We found an enrichment of recently duplicated, adaptive genes in similar functional pathways, but differential retention of environment-specific adaptive genes within each genome. For example, large expansions of NLR genes were discovered in the two wild grapevine genomes studied. Our findings support variation in transposable elements contributed to unique traits in grapevines. Our work revealed gene zygosity, segmental duplications, gene gain-and-loss variations, and transposable element polymorphisms can be key driving forces for grapevine genome diversification.

Whole genome doubling-induced the enrichment of H3K27me3 in genes carrying specific TEs in Aegilops tauschii

Polyploidization plays important roles in the evolution and breeding of the common wheat. Aegilops tauschii, the D-genome progenitor of the common wheat, provides a valuable pool of resistance genes to multiple diseases. Extensive studies focus on the exploration of these genes for wheat improvement. However, few studies have unveiled alternations on genome-wide expression pattern and histone modifications induced by whole-genome doubling (WGD) process. In this study, we conducted transcriptome analysis for the diploid and tetraploid Ae. taushcii lines using the leaf and root tissues. Both lines tend to display similar tissue-specific pattern. Interestingly, we found that TEs located in genic regions were depleted of the repressive histone mark H3K27me3, whereas their adjacent chromatin was enriched with H3K27me3. The tetraploid line exhibited higher levels of H3K27me3 in those regions than the diploid line, particularly for genic regions associated with TEs of the long interspersed nuclear elements (LINEs), CACTA, PIF/Harbinger, Tc1/Mariner and unclassed DNA transposon. Surprisingly, the expression levels of these TEs cognate genes were negatively associated with the levels of H3K27me3 between the tetraploid and diploid lines, suggesting the five types of TEs located within genic regions might be involved in the regulation of the ploidy-related gene expression, possibly through differential enrichment of H3K27me3 in the genic regions. These findings will help to understand the potential role of specific types of TEs on transcription in response to WGD.

Methylome, transcriptome, and phenotype changes induced by temperature conditions experienced during sexual reproduction in Fragaria vesca

Temperature conditions experienced during embryogenesis and seed development may induce epigenetic changes that increase phenotypic variation in plants. Here we investigate if embryogenesis and seed development at two different temperatures (28 vs. 18°C) result in lasting phenotypic effects and DNA methylation changes in woodland strawberry (Fragaria vesca). Using five European ecotypes from Spain (ES12), Iceland (ICE2), Italy (IT4), and Norway (NOR2 and NOR29), we found statistically significant differences between plants from seeds produced at 18 or 28°C in three of four phenotypic features investigated under common garden conditions. This indicates the establishment of a temperature-induced epigenetic memory-like response during embryogenesis and seed development. The memory effect was significant in two ecotypes: in NOR2 flowering time, number of growth points and petiole length were affected, and in ES12 number of growth points was affected. This indicates that genetic differences between ecotypes in their epigenetic machinery, or other allelic differences, impact this type of plasticity. We observed statistically significant differences between ecotypes in DNA methylation marks in repetitive elements, pseudogenes, and genic elements. Leaf transcriptomes were also affected by embryonic temperature in an ecotype-specific manner. Although we observed significant and lasting phenotypic change in at least some ecotypes, there was considerable variation in DNA methylation between individual plants within each temperature treatment. This within-treatment variability in DNA methylation marks in F. vesca progeny may partly be a result of allelic redistribution from recombination during meiosis and subsequent epigenetic reprogramming during embryogenesis.

Long-read direct RNA sequencing reveals epigenetic regulation of chimeric gene-transposon transcripts in Arabidopsis thaliana

Transposable elements (TEs) are accumulated in both intergenic and intragenic regions in plant genomes. Intragenic TEs often act as regulatory elements of associated genes and are also co-transcribed with genes, generating chimeric TE-gene transcripts. Despite the potential impact on mRNA regulation and gene function, the prevalence and transcriptional regulation of TE-gene transcripts are poorly understood. By long-read direct RNA sequencing and a dedicated bioinformatics pipeline, ParasiTE, we investigated the transcription and RNA processing of TE-gene transcripts in Arabidopsis thaliana. We identified a global production of TE-gene transcripts in thousands of A. thaliana gene loci, with TE sequences often being associated with alternative transcription start sites or transcription termination sites. The epigenetic state of intragenic TEs affects RNAPII elongation and usage of alternative poly(A) signals within TE sequences, regulating alternative TE-gene isoform production. Co-transcription and inclusion of TE-derived sequences into gene transcripts impact regulation of RNA stability and environmental responses of some loci. Our study provides insights into TE-gene interactions that contributes to mRNA regulation, transcriptome diversity, and environmental responses in plants.

Whole genome sequencing and analysis of Armillaria gallica Jzi34 symbiotic with Gastrodia elata

BACKGROUND

Armillaria species are plant pathogens, but a few Armillaria species can establish a symbiotic relationship with Gastrodia elata, a rootless and leafless orchid, that is used as a Chinese herbal medicine. Armillaria is a nutrient source for the growth of G. elata. However, there are few reports on the molecular mechanism of symbiosis between Armillaria species and G. elata. The genome sequencing and analysis of Armillaria symbiotic with G. elata would provide genomic information for further studying the molecular mechanism of symbiosis.

RESULTS

The de novo genome assembly was performed with the PacBio Sequel platform and Illumina NovaSeq PE150 for the A. gallica Jzi34 strain, which was symbiotic with G. elata. Its genome assembly contained ~ 79.9 Mbp and consisted of 60 contigs with an N50 of 2,535,910 bp. There were only 4.1% repetitive sequences in the genome assembly. Functional annotation analysis revealed a total of 16,280 protein coding genes. Compared with the other five genomes of Armillaria, the carbohydrate enzyme gene family of the genome was significantly contracted, while it had the largest set of glycosyl transferase (GT) genes. It also had an expansion of auxiliary activity enzymes AA3-2 gene subfamily and cytochrome P450 genes. The synteny analysis result of P450 genes reveals that the evolutionary relationship of P450 proteins between A. gallica Jzi34 and other four Armillaria was complex.

CONCLUSIONS

These characteristics may be beneficial for establishing a symbiotic relationship with G. elata. These results explore the characteristics of A. gallica Jzi34 from a genomic perspective and provide an important genomic resource for further detailed study of Armillaria. This will help to further study the symbiotic mechanism between A. gallica and G. elata.

cis-Regulatory Elements in Plant Development, Adaptation, and Evolution

cis-Regulatory elements encode the genomic blueprints that ensure the proper spatiotemporal patterning of gene expression necessary for appropriate development and responses to the environment. Accumulating evidence implicates changes to gene expression as a major source of phenotypic novelty in eukaryotes, including acute phenotypes such as disease and cancer in mammals. Moreover, genetic and epigenetic variation affecting cis-regulatory sequences over longer evolutionary timescales has become a recurring theme in studies of morphological divergence and local adaptation. Here, we discuss the functions of and methods used to identify various classes of cis-regulatory elements, as well as their role in plant development and response to the environment. We highlight opportunities to exploit cis-regulatory variants underlying plant development and environmental responses for crop improvement efforts. Although a comprehensive understanding of cis-regulatory mechanisms in plants has lagged behind that in animals, we showcase several breakthrough findings that have profoundly influenced plant biology and shaped the overall understanding of transcriptional regulation in eukaryotes.

Plant polymerase IV sensitizes chromatin through histone modifications to preclude spread of silencing into protein-coding domains

Across eukaryotes, gene regulation is manifested via chromatin states roughly distinguished as heterochromatin and euchromatin. The establishment, maintenance, and modulation of the chromatin states is mediated using several factors including chromatin modifiers. However, factors that avoid the intrusion of silencing signals into protein-coding genes are poorly understood. Here we show that a plant specific paralog of RNA polymerase (Pol) II, named Pol IV, is involved in avoidance of facultative heterochromatic marks in protein-coding genes, in addition to its well-established functions in silencing repeats and transposons. In its absence, H3K27 trimethylation (me3) mark intruded the protein-coding genes, more profoundly in genes embedded with repeats. In a subset of genes, spurious transcriptional activity resulted in small(s) RNA production, leading to post-transcriptional gene silencing. We show that such effects are significantly pronounced in rice, a plant with a larger genome with distributed heterochromatin compared with Arabidopsis Our results indicate the division of labor among plant-specific polymerases, not just in establishing effective silencing via sRNAs and DNA methylation but also in influencing chromatin boundaries.

The power of retrotransposons in high-throughput genotyping and sequencing

The use of molecular markers has become an essential part of molecular genetics through their application in numerous fields, which includes identification of genes associated with targeted traits, operation of backcrossing programs, modern plant breeding, genetic characterization, and marker-assisted selection. Transposable elements are a core component of all eukaryotic genomes, making them suitable as molecular markers. Most of the large plant genomes consist primarily of transposable elements; variations in their abundance contribute to most of the variation in genome size. Retrotransposons are widely present throughout plant genomes, and replicative transposition enables them to insert into the genome without removing the original elements. Various applications of molecular markers have been developed that exploit the fact that these genetic elements are present everywhere and their ability to stably integrate into dispersed chromosomal localities that are polymorphic within a species. The ongoing development of molecular marker technologies is directly related to the deployment of high-throughput genotype sequencing platforms, and this research is of considerable significance. In this review, the practical application to molecular markers, which is a use of technology of interspersed repeats in the plant genome were examined using genomic sources from the past to the present. Prospects and possibilities are also presented.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Variation in mutation, recombination, and transposition rates in Drosophila melanogaster and Drosophila simulans

The rates of mutation, recombination, and transposition are core parameters in models of evolution. They impact genetic diversity, responses to ongoing selection, and levels of genetic load. However, even for key evolutionary model species such as Drosophila melanogaster and Drosophila simulans, few estimates of these parameters are available, and we have little idea of how rates vary between individuals, sexes, or populations. Knowledge of this variation is fundamental for parameterizing models of genome evolution. Here, we provide direct estimates of mutation, recombination, and transposition rates and their variation in a West African and a European population of D. melanogaster and a European population of D. simulans Across 89 flies, we observe 58 single-nucleotide mutations, 286 crossovers, and 89 transposable element (TE) insertions. Compared to the European D. melanogaster, we find the West African population has a lower mutation rate (1.67 × 10-9 site-1 gen-1 vs. 4.86 × 10-9 site-1 gen-1) and a lower transposition rate (8.99 × 10-5 copy-1 gen-1 vs. 23.36 × 10-5 copy-1 gen-1), but a higher recombination rate (3.44 cM/Mb vs. 2.06 cM/Mb). The European D. simulans population has a similar mutation rate to European D. melanogaster, but a significantly higher recombination rate and a lower, but not significantly different, transposition rate. Overall, we find paternal-derived mutations are more frequent than maternal ones in both species. Our study quantifies the variation in rates of mutation, recombination, and transposition among different populations and sexes, and our direct estimates of these parameters in D. melanogaster and D. simulans will benefit future studies in population and evolutionary genetics.

Maize PIMT2 repairs damaged 3-METHYLCROTONYL COA CARBOXYLASE in mitochondria, affecting seed vigor

PROTEIN l-ISOASPARTYL O-METHYLTRANSFERASE (PIMT) affects seed vigor by repairing damaged proteins. While PIMT is capable of isoaspartyl (isoAsp) repair in all proteins, those proteins most susceptible to isoAsp formation have not been well characterized, and the mechanisms by which PIMT affects seed vigor remain largely unknown. Using co-immunoprecipitation and LC-MS/MS, we found that maize (Zea mays) PIMT2 (ZmPIMT2) interacted predominantly with both subunits of maize 3-METHYLCROTONYL COA CARBOXYLASE (ZmMCC). ZmPIMT2 is specifically expressed in the maize embryo. Both mRNA and protein levels of ZmPIMT2 increased during seed maturation and declined during imbibition. Maize seed vigor was decreased in the zmpimt2 mutant line, while overexpression of ZmPIMT2 in maize and Arabidopsis thaliana increased seed vigor upon artificial aging. ZmPIMT2 was localized in the mitochondria, as determined by subcellular localization assays using maize protoplasts. ZmPIMT2 binding to ZmMCCα was confirmed by luciferase complementation tests in both tobacco (Nicotiana benthamiana) leaves and maize protoplasts. Knockdown of ZmMCCα decreased maize seed aging tolerance. Furthermore, overexpression of ZmPIMT2 decreased the accumulation of isoAsp of ZmMCCα protein in seed embryos that underwent accelerated aging treatment. Taken together, our results demonstrate that ZmPIMT2 binds ZmMCCα in mitochondria, repairs isoAsp damage, and positively affects maize seed vigor.

Miniature Inverted-Repeat Transposable Elements: Small DNA Transposons That Have Contributed to Plant MICRORNA Gene Evolution

Angiosperms form the largest phylum within the Plantae kingdom and show remarkable genetic variation due to the considerable difference in the nuclear genome size of each species. Transposable elements (TEs), mobile DNA sequences that can amplify and change their chromosome position, account for much of the difference in nuclear genome size between individual angiosperm species. Considering the dramatic consequences of TE movement, including the complete loss of gene function, it is unsurprising that the angiosperms have developed elegant molecular strategies to control TE amplification and movement. Specifically, the RNA-directed DNA methylation (RdDM) pathway, directed by the repeat-associated small-interfering RNA (rasiRNA) class of small regulatory RNA, forms the primary line of defense to control TE activity in the angiosperms. However, the miniature inverted-repeat transposable element (MITE) species of TE has at times avoided the repressive effects imposed by the rasiRNA-directed RdDM pathway. MITE proliferation in angiosperm nuclear genomes is due to their preference to transpose within gene-rich regions, a pattern of transposition that has enabled MITEs to gain further transcriptional activity. The sequence-based properties of a MITE results in the synthesis of a noncoding RNA (ncRNA), which, after transcription, folds to form a structure that closely resembles those of the precursor transcripts of the microRNA (miRNA) class of small regulatory RNA. This shared folding structure results in a MITE-derived miRNA being processed from the MITE-transcribed ncRNA, and post-maturation, the MITE-derived miRNA can be used by the core protein machinery of the miRNA pathway to regulate the expression of protein-coding genes that harbor homologous MITE insertions. Here, we outline the considerable contribution that the MITE species of TE have made to expanding the miRNA repertoire of the angiosperms.

Transposable elements are associated with genome-specific gene expression in bread wheat

INTRODUCTION

Recent studies in wheat emphasized the importance of TEs, which occupy ~85% of the wheat genome, as a major source of intraspecific genetic variation due to their recent activity and involvement in genomic rearrangements. The contribution of TEs to structural and functional variations in bread wheat genes is not fully understood.

METHODS

Here, publicly available RNA-Seq databases of bread wheat were integrated to identify TE insertions within gene bodies (exons\ introns) and assess the impact of TE insertions on gene expression variations of homoeologs gene groups. Overall, 70,818 homoeologs genes were analyzed: 55,170 genes appeared in each one of the three subgenomes (termed ABD), named triads; 12,640 genes appeared in two of the three subgenomes (in A and B only, termed AB; or in A and D only, termed AD; or in B and D only, termed BD);, named dyads; and 3,008 genes underwent duplication in one of the three subgenomes (two copies in: subgenome A, termed AABD; subgenome B, termed ABBD; or subgenome D, termed ABDD), named tetrads.

RESULTS

To this end, we found that ~36% of the 70,818 genes contained at least one TE insertion within the gene body, mostly in triads. Analysis of 14,258 triads revealed that the presence of TE insertion in at least one of the triad genes (7,439 triads) was associated with balanced expression (similar expression levels) between the homoeolog genes. TE insertions within the exon or in the untranslated regions (UTRs) of one or more of the homoeologs in a triad were significantly associated with homoeolog expression bias. Furthermore, we found a statistically significant correlation between the presence\absence of TEs insertions belonging to six TE superfamilies and 17 TE subfamilies and the suppression of a single homoeolog gene. A significant association was observed between the presence of TE insertions from specific superfamilies and the expression of genes that are associated with biotic and abiotic stress responses.

CONCLUSION

Our data strongly indicate that TEs might play a prominent role in controlling gene expression in a genome-specific manner in bread wheat.

A novel active transposon creates allelic variation through altered translation rate to influence protein abundance

Protein translation is tightly and precisely controlled by multiple mechanisms including upstream open reading frames (uORFs), but the origins of uORFs and their role in maize are largely unexplored. In this study, an active transposition event was identified during the propagation of maize inbred line B73. The transposon, which was named BTA for 'B73 active transposable element hAT', creates a novel dosage-dependent hypomorphic allele of the hexose transporter gene ZmSWEET4c through insertion within the coding sequence in the first exon, and results in reduced kernel size. The BTA insertion does not affect transcript abundance but reduces protein abundance of ZmSWEET4c, probably through the introduction of a uORF. Furthermore, the introduction of BTA sequence in the exon of other genes can regulate translation efficiency without affecting their mRNA levels. A transposon capture assay revealed 79 novel insertions for BTA and BTA-like elements. These insertion sites have typical euchromatin features, including low levels of DNA methylation and high levels of H3K27ac. A putative autonomous element that mobilizes BTA and BTA-like elements was identified. Together, our results suggest a transposon-based origin of uORFs and document a new role for transposable elements to influence protein abundance and phenotypic diversity by affecting the translation rate.

Expression pattern of resynthesized allotetraploid Capsella is determined by hybridization, not whole-genome duplication

Polyploidization, the process leading to the increase in chromosome sets, is a major evolutionary transition in plants. Whole-genome duplication (WGD) within the same species gives rise to autopolyploids, whereas allopolyploids result from a compound process with two distinct components: WGD and interspecific hybridization. To dissect the instant effects of WGD and hybridization on gene expression and phenotype, we created a series of synthetic hybrid and polyploid Capsella plants, including diploid hybrids, autotetraploids of both parental species, and two kinds of resynthesized allotetraploids with different orders of WGD and hybridization. Hybridization played a major role in shaping the relative expression pattern of the neo-allopolyploids, whereas WGD had almost no immediate effect on relative gene expression pattern but, nonetheless, still affected phenotypes. No transposable element-mediated genomic shock scenario was observed in either neo-hybrids or neo-polyploids. Finally, WGD and hybridization interacted and the distorting effects of WGD were less strong in hybrids. Whole-genome duplication may even improve hybrid fertility. In summary, while the initial relative gene expression pattern in neo-allotetraploids was almost entirely determined by hybridization, WGD only had trivial effects on relative expression patterns, both processes interacted and had a strong impact on physical attributes and meiotic behaviors.

Genome assembly of an Australian native grass species reveals a recent whole-genome duplication and biased gene retention of genes involved in stress response

BACKGROUND

The adaptive significance of polyploidy has been extensively debated, and chromosome-level genome assemblies of polyploids can provide insight into this. The Australian grass Bothriochloa decipiens belongs to the BCD clade, a group with a complex history of hybridization and polyploid. This is the first genome assembly and annotation of a species that belongs to this fascinating yet complex group.

FINDINGS

Using Illumina short reads, 10X Genomics linked reads, and Hi-C sequencing data, we assembled a highly contiguous genome of B. decipiens, with a total length of 1,218.22 Mb and scaffold N50 of 42.637 Mb. Comparative analysis revealed that the species experienced a relatively recent whole-genome duplication. We clustered the 20 major scaffolds, representing the 20 chromosomes, into the 2 subgenomes of the parental species using unique repeat signatures. We found evidence of biased fractionation and differences in the activity of transposable elements between the subgenomes prior to hybridization. Duplicates were enriched for genes involved in transcription and response to external stimuli, supporting a biased retention of duplicated genes following whole-genome duplication.

CONCLUSIONS

Our results support the hypotheses of a biased retention of duplicated genes following polyploidy and point to differences in repeat activity associated with subgenome dominance. B. decipiens is a widespread species with the ability to establish across many soil types, making it a prime candidate for climate change- resilient ecological restoration of Australian grasslands. This reference genome is a valuable resource for future population genomic research on Australian grasses.

Jack of all trades: Genome assembly of Wild Jack and comparative genomics of Artocarpus

Artocarpus (Moraceae), known as breadfruits for their diverse nutritious fruits, is prized for its high-quality timber, medicinal value, and economic importance. Breadfruits are native to Southeast Asia but have been introduced to other continents. The most commonly cultivated species are Artocarpus heterophyllus (Jackfruit) and Artocarpus altilis (Breadfruit). With numerous smaller but nutritionally comparable fruits on a larger tree, Artocarpus hirsutus, also called "Wild Jack" or "Ayani", is an elusive forest species endemic to Indian Western Ghats. In this study, we sequenced and assembled the whole genome of Artocarpus hirsutus sampled from the sacred groves of Coorg, India. To decipher demographic and evolutionary history, we compared our Wild Jack genome with previously published Jackfruit and Breadfruit genomes. Demographic history reconstruction indicates a stronger effect of habitat rather than phylogeny on the population histories of these plants. Repetitive genomic regions, especially LTR Copia, strongly affected the demographic trajectory of A. heterophyllus. Upon further investigation, we found a recent lineage-specific accumulation of LTR Copia in A. heterophyllus, which had a major contribution to its larger genome size. Several genes from starch, sucrose metabolism, and plant hormone signal transduction pathways, in Artocarpus species had signatures of selection and gene family evolution. Our comparative genomic framework provides important insights by incorporating endemic species such as the Wild Jack.

Pepper variome reveals the history and key loci associated with fruit domestication and diversification

Pepper (Capsicum spp.) is an important vegetable crop that provides a unique pungent sensation when eaten. Through construction of a pepper variome map, we examined the main groups that emerged during domestication and breeding of C. annuum, their relationships and temporal succession, and the molecular events underlying the main transitions. The results showed that the initial differentiation in fruit shape and pungency, increase in fruit weight, and transition from erect to pendent fruits, as well as the recent appearance of large, blocky, sweet fruits (bell peppers), were accompanied by strong selection/fixation of key alleles and introgressions in two large genomic regions. Furthermore, we identified Up, which encodes a BIG GRAIN protein involved in auxin transport, as a key domestication gene that controls erect vs pendent fruit orientation. The up mutation gained increased expression especially in the fruit pedicel through a 579-bp sequence deletion in its 5' upstream region, resulting in the phenotype of pendent fruit. The function of Up was confirmed by virus-induced gene silencing. Taken together, these findings constitute a cornerstone for understanding the domestication and differentiation of a key horticultural crop.

Genomic landscapes of bacterial transposons and their applications in strain improvement

Transposons are mobile genetic elements that can give rise to gene mutation and genome rearrangement. Due to their mobility, transposons have been exploited as genetic tools for modification of plants, animals, and microbes. Although a plethora of reviews have summarized families of transposons, the transposons from fermentation bacteria have not been systematically documented, which thereby constrain the exploitation for metabolic engineering and synthetic biology purposes. In this review, we summarize the transposons from the most used fermentation bacteria including Escherichia coli, Bacillus subtilis, Lactococcus lactis, Corynebacterium glutamicum, Klebsiella pneumoniae, and Zymomonas mobilis by literature retrieval and data mining from GenBank and KEGG. We also outline the state-of-the-art advances in basic research and industrial applications especially when allied with other genetic tools. Overall, this review aims to provide valuable insights for transposon-mediated strain improvement. KEY POINTS: • The transposons from the most-used fermentation bacteria are systematically summarized. • The applications of transposons in strain improvement are comprehensively reviewed.

Impact of LTR-Retrotransposons on Genome Structure, Evolution, and Function in Curcurbitaceae Species

Long terminal repeat (LTR)-retrotransposons (LTR-RTs) comprise a major portion of many plant genomes and may exert a profound impact on genome structure, function, and evolution. Although many studies have focused on these elements in an individual species, their dynamics on a family level remains elusive. Here, we investigated the abundance, evolutionary dynamics, and impact on associated genes of LTR-RTs in 16 species in an economically important plant family, Cucurbitaceae. Results showed that full-length LTR-RT numbers and LTR-RT content varied greatly among different species, and they were highly correlated with genome size. Most of the full-length LTR-RTs were amplified after the speciation event, reflecting the ongoing rapid evolution of these genomes. LTR-RTs highly contributed to genome size variation via species-specific distinct proliferations. The Angela and Tekay lineages with a greater evolutionary age were amplified in Trichosanthes anguina, whereas a recent activity burst of Reina and another ancient round of Tekay activity burst were examined in Sechium edule. In addition, Tekay and Retand lineages belonging to the Gypsy superfamily underwent a recent burst in Gynostemma pentaphyllum. Detailed investigation of genes with intronic and promoter LTR-RT insertion showed diverse functions, but the term of metabolism was enriched in most species. Further gene expression analysis in G.pentaphyllum revealed that the LTR-RTs within introns suppress the corresponding gene expression, whereas the LTR-RTs within promoters exert a complex influence on the downstream gene expression, with the main function of promoting gene expression. This study provides novel insights into the organization, evolution, and function of LTR-RTs in Cucurbitaceae genomes.

Exploiting the miniature inverted-repeat transposable elements insertion polymorphisms as an efficient DNA marker system for genome analysis and evolutionary studies in wheat and related species

Transposable elements (TEs) constitute ~80% of the complex bread wheat genome and contribute significantly to wheat evolution and environmental adaptation. We studied 52 TE insertion polymorphism markers to ascertain their efficiency as a robust DNA marker system for genetic studies in wheat and related species. Significant variation was found in miniature inverted-repeat transposable element (MITE) insertions in relation to ploidy with the highest number of "full site" insertions occurring in the hexaploids (32.6 ± 3.8), while the tetraploid and diploid progenitors had 22.3 ± 0.6 and 15.0 ± 3.5 "full sites," respectively, which suggested a recent rapid activation of these transposons after the formation of wheat. Constructed phylogenetic trees were consistent with the evolutionary history of these species which clustered mainly according to ploidy and genome types (SS, AA, DD, AABB, and AABBDD). The synthetic hexaploids sub-clustered near the tetraploid species from which they were re-synthesized. Preliminary genotyping in 104 recombinant inbred lines (RILs) showed predominantly 1:1 segregation for simplex markers, with four of these markers already integrated into our current DArT-and SNP-based linkage map. The MITE insertions also showed stability with no single excision observed. The MITE insertion site polymorphisms uncovered in this study are very promising as high-potential evolutionary markers for genomic studies in wheat.

Loss of pod strings in common bean is associated with gene duplication, retrotransposon insertion and overexpression of PvIND

Fruit development has been central in the evolution and domestication of flowering plants. In common bean (Phaseolus vulgaris), the principal global grain legume staple, two main production categories are distinguished by fibre deposition in pods: dry beans, with fibrous, stringy pods; and stringless snap/green beans, with reduced fibre deposition, which frequently revert to the ancestral stringy state. Here, we identify genetic and developmental patterns associated with pod fibre deposition. Transcriptional, anatomical, epigenetic and genetic regulation of pod strings were explored through RNA-seq, RT-qPCR, fluorescence microscopy, bisulfite sequencing and whole-genome sequencing. Overexpression of the INDEHISCENT ('PvIND') orthologue was observed in stringless types compared with isogenic stringy lines, associated with overspecification of weak dehiscence-zone cells throughout the pod vascular sheath. No differences in DNA methylation were correlated with this phenotype. Nonstringy varieties showed a tandemly direct duplicated PvIND and a Ty1-copia retrotransposon inserted between the two repeats. These sequence features are lost during pod reversion and are predictive of pod phenotype in diverse materials, supporting their role in PvIND overexpression and reversible string phenotype. Our results give insight into reversible gain-of-function mutations and possible genetic solutions to the reversion problem, of considerable economic value for green bean production.

Modeling chromatin state from sequence across angiosperms using recurrent convolutional neural networks

Accessible chromatin regions are critical components of gene regulation but modeling them directly from sequence remains challenging, especially within plants, whose mechanisms of chromatin remodeling are less understood than in animals. We trained an existing deep-learning architecture, DanQ, on data from 12 angiosperm species to predict the chromatin accessibility in leaf of sequence windows within and across species. We also trained DanQ on DNA methylation data from 10 angiosperms because unmethylated regions have been shown to overlap significantly with ACRs in some plants. The across-species models have comparable or even superior performance to a model trained within species, suggesting strong conservation of chromatin mechanisms across angiosperms. Testing a maize (Zea mays L.) held-out model on a multi-tissue chromatin accessibility panel revealed our models are best at predicting constitutively accessible chromatin regions, with diminishing performance as cell-type specificity increases. Using a combination of interpretation methods, we ranked JASPAR motifs by their importance to each model and saw that the TCP and AP2/ERF transcription factor (TF) families consistently ranked highly. We embedded the top three JASPAR motifs for each model at all possible positions on both strands in our sequence window and observed position- and strand-specific patterns in their importance to the model. With our publicly available across-species 'a2z' model it is now feasible to predict the chromatin accessibility and methylation landscape of any angiosperm genome.

Characterization of Transposon-Derived Accessible Chromatin Regions in Rice (Oryza Sativa)

Growing evidence indicates that transposons or transposable elements (TEs)-derived accessible chromatin regions (ACRs) play essential roles in multiple biological processes by interacting with trans-acting factors. However, the function of TE-derived ACRs in the regulation of gene expression in the rice genome has not been well characterized. In this study, we examined the chromatin dynamics in six types of rice tissues and found that ~8% of ACRs were derived from TEs and exhibited distinct levels of accessibility and conservation as compared to those without TEs. TEs exhibited a TE subtype-dependent impact on ACR formation, which can be mediated by changes in the underlying DNA methylation levels. Moreover, we found that tissue-specific TE-derived ACRs might function in the tissue development through the modulation of nearby gene expression. Interestingly, many genes in domestication sweeps were found to overlap with TE-derived ACRs, suggesting their potential functions in the rice domestication. In addition, we found that the expression divergence of 1070 duplicate gene pairs were associated with TE-derived ACRs and had distinct distributions of TEs and ACRs around the transcription start sites (TSSs), which may experience different selection pressures. Thus, our study provides some insights into the biological implications of TE-derived ACRs in the rice genome. Our results imply that these ACRs are likely involved in the regulation of tissue development, rice domestication and functional divergence of duplicated 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.

Hsa-miR-422a Originated from Short Interspersed Nuclear Element Increases ARID5B Expression by Collaborating with NF-E2

MicroRNAs (miRNAs) are a class of small non-coding RNAs that regulate the expression of target messenger RNA (mRNA) complementary to the 3' untranslated region (UTR) at the post-transcriptional level. Hsa-miR-422a, which is commonly known as miRNA derived from transposable element (MDTE), was derived from short interspersed nuclear element (SINE). Through expression analysis, hsa-miR-422a was found to be highly expressed in both the small intestine and liver of crab-eating monkey. AT-Rich Interaction Domain 5 B (ARID5B) was selected as the target gene of hsa-miR-422a, which has two binding sites in both the exon and 3'UTR of ARID5B. To identify the interaction between hsa-miR-422a and ARID5B, a dual luciferase assay was conducted in HepG2 cell line. The luciferase activity of cells treated with the hsa-miR-422a mimic was upregulated and inversely downregulated when both the hsa-miR-422a mimic and inhibitor were administered. Nuclear factor erythroid-2 (NF-E2) was selected as the core transcription factor (TF) via feed forward loop analysis. The luciferase expression was downregulated when both the hsa-miR-422a mimic and siRNA of NF-E2 were treated, compared to the treatment of the hsa-miR-422a mimic alone. The present study suggests that hsa-miR-422a derived from SINE could bind to the exon region as well as the 3'UTR of ARID5B. Additionally, hsa-miR-422a was found to share binding sites in ARID5Bwith several TFs, including NF-E2. The hsa-miR-422a might thus interact with TF to regulate the expression of ARID5B, as demonstrated experimentally. Altogether, hsa-miR-422a acts as a super enhancer miRNA of ARID5Bby collaborating with TF and NF-E2.

Neuronal role of taxi is imperative for flight in Drosophila melanogaster

Extensive studies in Drosophila have led to the elucidation of the roles of many molecular players involved in the sensorimotor coordination of flight. However, the identification and characterisation of new players can add novel perspectives to the process. In this paper, we show that the extant mutant, jumper, is a hypermorphic allele of the taxi/delilah gene, which encodes a transcription factor. The defective flight of jumper flies results from the insertion of an I-element in the 5'-UTR of taxi gene, leading to an over-expression of the taxi. We also show that the molecular lesion responsible for the taxi¹ allele results from a 25 bp deletion leading to a shift in the reading frame at the C-terminus of the taxi coding sequence. Thus, the last 20 residues are replaced by 32 disparate residues in taxi¹. Both taxi¹, a hypomorphic allele, and the CRISPR-Cas9 knock-out (taxi^KO) null allele, show a defective flight phenotype. Electrophysiological studies show taxi hypermorphs, hypomorphs, and knock out flies show abnormal neuronal firing. We further show that neuronal-specific knock-down or over-expression of taxi cause a defect in the brain's inputs to the flight muscles, leading to reduced flight ability. Through transcriptomic analysis of the taxi^KO fly head, we have identified several putative targets of Taxi that may play important roles in flight. In conclusion, from molecularly characterising jumper to establishing Taxi's role during Drosophila flight, our work shows that the forward genetics approach still can lead to the identification of novel molecular players required for neuronal transmission.

Variation in leaf transcriptome responses to elevated ozone corresponds with physiological sensitivity to ozone across maize inbred lines

We examine the impact of sustained elevated ozone concentration on the leaf transcriptome of 5 diverse maize inbred genotypes, which vary in physiological sensitivity to ozone (B73, Mo17, Hp301, C123, and NC338), using long reads to assemble transcripts and short reads to quantify expression of these transcripts. More than 99% of the long reads, 99% of the assembled transcripts, and 97% of the short reads map to both B73 and Mo17 reference genomes. Approximately 95% of the genes with assembled transcripts belong to known B73-Mo17 syntenic loci and 94% of genes with assembled transcripts are present in all temperate lines in the nested association mapping pan-genome. While there is limited evidence for alternative splicing in response to ozone stress, there is a difference in the magnitude of differential expression among the 5 genotypes. The transcriptional response to sustained ozone stress in the ozone resistant B73 genotype (151 genes) was modest, while more than 3,300 genes were significantly differentially expressed in the more sensitive NC338 genotype. There is the potential for tandem duplication in 30% of genes with assembled transcripts, but there is no obvious association between potential tandem duplication and differential expression. Genes with a common response across the 5 genotypes (83 genes) were associated with photosynthesis, in particular photosystem I. The functional annotation of genes not differentially expressed in B73 but responsive in the other 4 genotypes (789) identifies reactive oxygen species. This suggests that B73 has a different response to long-term ozone exposure than the other 4 genotypes. The relative magnitude of the genotypic response to ozone, and the enrichment analyses are consistent regardless of whether aligning short reads to: long read assembled transcripts; the B73 reference; the Mo17 reference. We find that prolonged ozone exposure directly impacts the photosynthetic machinery of the leaf.

Role of Transposable Elements in Genome Stability: Implications for Health and Disease

Most living organisms have in their genome a sizable proportion of DNA sequences capable of mobilization; these sequences are commonly referred to as transposons, transposable elements (TEs), or jumping genes. Although long thought to have no biological significance, advances in DNA sequencing and analytical technologies have enabled precise characterization of TEs and confirmed their ubiquitous presence across all forms of life. These findings have ignited intense debates over their biological significance. The available evidence now supports the notion that TEs exert major influence over many biological aspects of organismal life. Transposable elements contribute significantly to the evolution of the genome by giving rise to genetic variations in both active and passive modes. Due to their intrinsic nature of mobility within the genome, TEs primarily cause gene disruption and large-scale genomic alterations including inversions, deletions, and duplications. Besides genomic instability, growing evidence also points to many physiologically important functions of TEs, such as gene regulation through cis-acting control elements and modulation of the transcriptome through epigenetic control. In this review, we discuss the latest evidence demonstrating the impact of TEs on genome stability and the underling mechanisms, including those developed to mitigate the deleterious impact of TEs on genomic stability and human health. We have also highlighted the potential therapeutic application of TEs.

Chromosome-level genome assembly of the aquatic plant Nymphoides indica reveals transposable element bursts and NBS-LRR gene family expansion shedding light on its invasiveness

Nymphoides indica, an aquatic plant, is an invasive species that causes both ecological and economic damage in North America and elsewhere. However, the lack of genomic data of N. indica limits the in-depth analysis of this invasive species. Here, we report a chromosome-level genome assembly of nine pseudochromosomes of N. indica with a total size of ∼ 520 Mb. More than half of the N. indica genome consists of transposable elements (TEs), and a higher density of TEs around genes may play a significant role in response to an ever-changing environment by regulating the nearby gene. Additionally, our analysis revealed that N. indica only experienced a gamma (γ) whole-genome triplication event. Functional enrichment of the N. indica-specific and expanded gene families highlighted genes involved in the responses to hypoxia and plant-pathogen interactions, which may strengthen the ability to adapt to external challenges and improve ecological fitness. Furthermore, we identified 160 members of the nucleotide-binding site and leucine-rich repeat gene family, which may be linked to the defence response. Collectively, the high-quality N. indica genome reported here opens a novel avenue to understand the evolution and rapid invasion of Nymphoides spp.