Assessing the regulatory potential of transposable elements using chromatin accessibility profiles of maize transposons

The regulatory potential of transposable elements in maize

The genomes of flowering plants consist largely of transposable elements (TEs), some of which modulate gene regulation and function. However, the repetitive nature of TEs and difficulty of mapping individual TEs by short-read sequencing have hindered our understanding of their regulatory potential. Here we show that long-read chromatin fibre sequencing (Fiber-seq) comprehensively identifies accessible chromatin regions (ACRs) and CpG methylation across the maize genome. We uncover stereotypical ACR patterns at young TEs that degenerate with evolutionary age, resulting in TE enhancers preferentially marked by a novel plant-specific epigenetic feature: simultaneous hyper-CpG methylation and chromatin accessibility. We show that TE ACRs are co-opted as gene promoters and that ACR-containing TEs can facilitate gene amplification. Lastly, we uncover a pervasive epigenetic signature-hypo-5mCpG methylation and diffuse chromatin accessibility-directing TEs to specific loci, including the loci that sparked McClintock's discovery of TEs.

The genetic architecture of cell type-specific cis regulation in maize

Gene expression and complex phenotypes are determined by the activity of cis-regulatory elements. However, an understanding of how extant genetic variants affect cis regulation remains limited. Here, we investigated the consequences of cis-regulatory diversity using single-cell genomics of more than 0.7 million nuclei across 172 Zea mays (maize) inbreds. Our analyses pinpointed cis-regulatory elements distinct to domesticated maize and revealed how historical transposon activity has shaped the cis-regulatory landscape. Leveraging population genetics principles, we fine-mapped about 22,000 chromatin accessibility-associated genetic variants with widespread cell type-specific effects. Variants in TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR-binding sites were the most prevalent determinants of chromatin accessibility. Finally, integrating chromatin accessibility-associated variants, organismal trait variation, and population differentiation revealed how local adaptation has rewired regulatory networks in unique cellular contexts to alter maize flowering.

Understanding the Regulation Activities of Transposons in Driving the Variation and Evolution of Polyploid Plant Genome

Transposon is the main component of the eukaryotic genome, and more and more plant genome data show that transposons are diverse in regulating genome structure, variation, function and evolution, with different transposition mechanisms in the genome. Hybridization and polyploidy play an important role in promoting plant speciation and evolution, and recent studies have shown that polyploidy is usually accompanied by the expansion of transposons, which affect the genome size and structure of polyploid plants. Transposons can insert into genes and intergenic regions, resulting in great differences in the overall genome structure of closely related plant species, and it can also capture gene segments in the genome to increase the copy number of genes. In addition, transposons influence the epigenetic modification state of the genome and regulate the expression of the gene, while plant phenotype, biological and abiotic stress response are also regulated by transposons. Overall, transposons play an important role in the plant genome, especially polyploid plant genome, adaptation and evolution.

The genomic and epigenomic landscapes of hemizygous genes across crops with contrasting reproductive systems

Hemizygous genes, which are present on only one of the two homologous chromosomes of diploid organisms, have been mainly studied in the context of sex chromosomes and sex-linked genes. However, these genes can also occur on the autosomes of diploid plants due to structural variants (SVs), such as a deletion/insertion of one allele, and this phenomenon largely unexplored in plants. Here, we investigated the genomic and epigenomic landscapes of hemizygous genes across 22 genomes with varying propagation histories: eleven clonal lineages, seven outcrossed samples, and four inbred and putatively homozygous genomes. We identified SVs leading to genic hemizygosity. As expected, very few genes (0.01 to 1.2%) were hemizygous in the homozygous genomes, representing negative controls. Hemizygosity was appreciable among outcrossed lineages, averaging 8.7% of genes, but consistently elevated for the clonal samples at 13.8% genes, likely reflecting heterozygous SV accumulation during clonal propagation. Compared to diploid genes, hemizygous genes were more often situated in centromeric than telomeric regions and experienced weaker purifying selection. They also had reduced levels of expression, averaging ~20% of the expression levels of diploid genes, violating the evolutionary model of dosage compensation. We also detected higher DNA methylation levels in hemizygous genes and transposable elements, which may contribute to their reduced expression. Finally, expression profiles showed that hemizygous genes were more specifically expressed in contexts related to fruit development, organ differentiation, and stress responses. Overall, hemizygous genes accumulate in clonally propagated lineages and display distinct genetic and epigenetic features compared to diploid genes, shedding unique insights into genetic studies and breeding programs of clonal crops.

Chromatin Accessibility and Gene Expression Vary Between a New and Evolved Autopolyploid of Arabidopsis arenosa

Polyploids arise from whole-genome duplication (WGD) events, which have played important roles in genome evolution across eukaryotes. WGD can increase genome complexity, yield phenotypic novelty, and influence adaptation. Neo-polyploids have been reported to often show seemingly stochastic epigenetic and transcriptional changes, but this leaves open the question whether these changes persist in evolved polyploids. A powerful approach to address this is to compare diploids, neo-polyploids, and evolved polyploids of the same species. Arabidopsis arenosa is a species that allows us to do this-natural diploid and autotetraploid populations exist, while neo-tetraploids can be artificially generated. Here, we use ATAC-seq to assay local chromatin accessibility, and RNA-seq to study gene expression on matched leaf and petal samples from diploid, neo-tetraploid and evolved tetraploid A. arenosa. We found over 8,000 differentially accessible chromatin regions across all samples. These are largely tissue specific and show distinct trends across cytotypes, with roughly 70% arising upon WGD. Interestingly, only a small proportion is associated with expression changes in nearby genes. However, accessibility variation across cytotypes associates strongly with the number of nearby transposable elements. Relatively few genes were differentially expressed upon genome duplication, and ∼60% of these reverted to near-diploid levels in the evolved tetraploid, suggesting that most initial perturbations do not last. Our results provide new insights into how epigenomic and transcriptional mechanisms jointly respond to genome duplication and subsequent evolution of autopolyploids, and importantly, show that one cannot be directly predicted from the other.

The chromatin accessibility and transcriptomic landscape of the aging mice cochlea and the identification of potential functional super-enhancers in age-related hearing loss

BACKGROUND

Presbycusis, also referred to as age-related hearing loss (ARHL), is a condition that results from the cumulative effects of aging on an individual's auditory capabilities. Given the limited understanding of epigenetic mechanisms in ARHL, our research focuses on alterations in chromatin-accessible regions.

METHODS

We employed assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) in conjunction with unique identifier (UID) mRNA-seq between young and aging cochleae, and conducted integrated analysis as well as motif/TF-gene prediction. Additionally, the essential role of super-enhancers (SEs) in the development of ARHL was identified by comparative analysis to previous research. Meanwhile, an ARHL mouse model and an aging mimic hair cell (HC) model were established with a comprehensive identification of senescence phenotypes to access the role of SEs in ARHL progression.

RESULTS

The control cochlear tissue exhibited greater chromatin accessibility than cochlear tissue affected by ARHL. Furthermore, the levels of histone 3 lysine 27 acetylation were significantly depressed in both aging cochlea and aging mimic HEI-OC1 cells, highlighting the essential role of SEs in the development of ARHL. The potential senescence-associated super-enhancers (SASEs) of ARHL were identified, most of which exhibited decreased chromatin accessibility. The majority of genes related to the SASEs showed obvious decreases in mRNA expression level in aging HCs and was noticeably altered following treatment with JQ1 (a commonly used SE inhibitor).

CONCLUSION

The chromatin accessibility in control cochlear tissue was higher than that in cochlear tissue affected by ARHL. Potential SEs involved in ARHL were identified, which might provide a basis for future therapeutics targeting SASEs related to ARHL.

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.

MINI-AC: inference of plant gene regulatory networks using bulk or single-cell accessible chromatin profiles

Gene regulatory networks (GRNs) represent the interactions between transcription factors (TF) and their target genes. Plant GRNs control transcriptional programs involved in growth, development, and stress responses, ultimately affecting diverse agricultural traits. While recent developments in accessible chromatin (AC) profiling technologies make it possible to identify context-specific regulatory DNA, learning the underlying GRNs remains a major challenge. We developed MINI-AC (Motif-Informed Network Inference based on Accessible Chromatin), a method that combines AC data from bulk or single-cell experiments with TF binding site (TFBS) information to learn GRNs in plants. We benchmarked MINI-AC using bulk AC datasets from different Arabidopsis thaliana tissues and showed that it outperforms other methods to identify correct TFBS. In maize, a crop with a complex genome and abundant distal AC regions, MINI-AC successfully inferred leaf GRNs with experimentally confirmed, both proximal and distal, TF-target gene interactions. Furthermore, we showed that both AC regions and footprints are valid alternatives to infer AC-based GRNs with MINI-AC. Finally, we combined MINI-AC predictions from bulk and single-cell AC datasets to identify general and cell-type specific maize leaf regulators. Focusing on C4 metabolism, we identified diverse regulatory interactions in specialized cell types for this photosynthetic pathway. MINI-AC represents a powerful tool for inferring accurate AC-derived GRNs in plants and identifying known and novel candidate regulators, improving our understanding of gene regulation in plants.

Genomic variations combined with epigenetic modifications rewire open chromatin in rice

Cis-regulatory elements (CREs) fine-tune gene transcription in eukaryotes. CREs with sequence variations play vital roles in driving plant or crop domestication. However, how global sequence and structural variations (SVs) are responsible for multilevel changes between indica and japonica rice (Oryza sativa) is still not fully elucidated. To address this, we conducted multiomic studies using MNase hypersensitivity sequencing (MH-seq) in combination with RNA sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq), and bisulfite sequencing (BS-seq) between the japonica rice variety Nipponbare (NIP) and indica rice variety 93-11. We found that differential MNase hypersensitive sites (MHSs) exhibited some distinct intrinsic genomic sequence features between NIP and 93-11. Notably, through MHS-genome-wide association studies (GWAS) integration, we found that key sequence variations may be associated with differences of agronomic traits between NIP and 93-11, which is partly achieved by MHSs harboring CREs. In addition, SV-derived differential MHSs caused by transposable element (TE) insertion, especially by noncommon TEs among rice varieties, were associated with genes with distinct functions, indicating that TE-driven gene neo- or subfunctionalization is mediated by changes of chromatin openness. This study thus provides insights into how sequence and genomic SVs control agronomic traits of NIP and 93-11; it also provides genome-editing targets for molecular breeding aiming at improving favorable agronomic properties.

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.

Chromosomal-level genome assembly of the high-quality Xian/Indica rice (Oryza sativa L.) Xiangyaxiangzhan

The indica rice variety XYXZ carries elite traits including appearance and eating quality. Here, we report the de novo assembly of XYXZ using Illumine paired-end whole-genome shotgun sequencing and Nanopore sequencing. We annotated 39,722 protein-coding genes in the 395.04 Mb assembly. In comparison to other cultivars, XYXZ showed a larger gene size including the transcripts and introns, and more exons per gene. And hundreds of ultra-long genes were also detected. A total of 4362 complete LTRs were annotated, and among them, many were located next to or in protein-coding genes including several genes related to rice quality. We observed the different distributions of LTRs in these genes among XYXZ, Nipponbare, and R498, implying these LTRs might potentially affect expressions of the proximal genes and rice quality. Overall, This chromosome-length genome assembly of XYXZ provides a valuable resource for gene discovery, genetic variation and evolution, and the breeding of high-quality rice.

Annotation of the Turnera subulata (Passifloraceae) Draft Genome Reveals the S-Locus Evolved after the Divergence of Turneroideae from Passifloroideae in a Stepwise Manner

A majority of Turnera species (Passifloraceae) exhibit distyly, a reproductive system involving both self-incompatibility and reciprocal herkogamy. This system differs from self-incompatibility in Passiflora species. The genetic basis of distyly in Turnera is a supergene, restricted to the S-morph, and containing three S-genes. How supergenes and distyly evolved in Turnera, and the other Angiosperm families exhibiting distyly remain largely unknown. Unraveling the evolutionary origins in Turnera requires the generation of genomic resources and extensive phylogenetic analyses. Here, we present the annotated draft genome of the S-morph of distylous Turnera subulata. Our annotation allowed for phylogenetic analyses of the three S-genes' families across 56 plant species ranging from non-seed plants to eudicots. In addition to the phylogenetic analysis, we identified the three S-genes' closest paralogs in two species of Passiflora. Our analyses suggest that the S-locus evolved after the divergence of Passiflora and Turnera. Finally, to provide insights into the neofunctionalization of the S-genes, we compared expression patterns of the S-genes with close paralogs in Arabidopsis and Populus trichocarpa. The annotation of the T. subulata genome will provide a useful resource for future comparative work. Additionally, this work has provided insights into the convergent nature of distyly and the origin of supergenes.

Graph-based pan-genomes: increased opportunities in plant genomics

Due to the development of sequencing technology and the great reduction in sequencing costs, an increasing number of plant genomes have been assembled, and numerous genomes have revealed large amounts of variations. However, a single reference genome does not allow the exploration of species diversity, and therefore the concept of pan-genome was developed. A pan-genome is a collection of all sequences available for a species, including a large number of consensus sequences, large structural variations, and small variations including single nucleotide polymorphisms and insertions/deletions. A simple linear pan-genome does not allow these structural variations to be intuitively characterized, so graph-based pan-genomes have been developed. These pan-genomes store sequence and structural variation information in the form of nodes and paths to store and display species variation information in a more intuitive manner. The key role of graph-based pan-genomes is to expand the coordinate system of the linear reference genome to accommodate more regions of genetic diversity. Here, we review the origin and development of graph-based pan-genomes, explore their application in plant research, and further highlight the application of graph-based pan-genomes for future plant breeding.

Canalization of genome-wide transcriptional activity in Arabidopsis thaliana accessions by MET1-dependent CG methylation

BACKGROUND

Despite its conserved role on gene expression and transposable element (TE) silencing, genome-wide CG methylation differs substantially between wild Arabidopsis thaliana accessions.

RESULTS

To test our hypothesis that global reduction of CG methylation would reduce epigenomic, transcriptomic, and phenotypic diversity in A. thaliana accessions, we knock out MET1, which is required for CG methylation, in 18 early-flowering accessions. Homozygous met1 mutants in all accessions suffer from common developmental defects such as dwarfism and delayed flowering, in addition to accession-specific abnormalities in rosette leaf architecture, silique morphology, and fertility. Integrated analysis of genome-wide methylation, chromatin accessibility, and transcriptomes confirms that MET1 inactivation greatly reduces CG methylation and alters chromatin accessibility at thousands of loci. While the effects on TE activation are similarly drastic in all accessions, the quantitative effects on non-TE genes vary greatly. The global expression profiles of accessions become considerably more divergent from each other after genome-wide removal of CG methylation, although a few genes with diverse expression profiles across wild-type accessions tend to become more similar in mutants. Most differentially expressed genes do not exhibit altered chromatin accessibility or CG methylation in cis, suggesting that absence of MET1 can have profound indirect effects on gene expression and that these effects vary substantially between accessions.

CONCLUSIONS

Systematic analysis of MET1 requirement in different A. thaliana accessions reveals a dual role for CG methylation: for many genes, CG methylation appears to canalize expression levels, with methylation masking regulatory divergence. However, for a smaller subset of genes, CG methylation increases expression diversity beyond genetically encoded differences.

The evolution and function of transposons in epigenetic regulation in response to the environment

Transposable elements (TEs) make up a major proportion of plant genomes. Despite their prevalence genome-wide, TEs are often tossed aside as "junk DNA" since they rarely cause phenotypes, and epigenetic mechanisms silence TEs to prevent them from causing deleterious mutations through movement. While this bleak picture of TEs in genomes is true on average, a growing number of examples across many plant species point to TEs as drivers of phenotypic diversity and novel stress responses. Examples of TE-influenced phenotypes illustrate the many ways that novel transposition events can alter local gene expression and how this relates to potential variation in plant responses to environmental stress. Since TE families and insertions at the locus level lack evolutionary conservation, advancements in the field will require TE experts across diverse species to identify and utilize TE variation in their own systems as a means of crop improvement.

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.

Asymmetric subgenomic chromatin architecture impacts on gene expression in resynthesized and natural allopolyploid Brassica napus

Although asymmetric subgenomic epigenetic modification and gene expression have been revealed in the successful establishment of allopolyploids, the changes in chromatin accessibility and their relationship with epigenetic modifications and gene expression are poorly understood. Here, we synthetically analyzed chromatin accessibility, four epigenetic modifications and gene expression in natural allopolyploid Brassica napus, resynthesized allopolyploid B. napus, and diploid progenitors (B. rapa and B. oleracea). “Chromatin accessibility shock” occurred in both allopolyploidization and natural evolutionary processes, and genic accessible chromatin regions (ACRs) increased after allopolyploidization. ACRs associated with H3K27me3 modifications were more accessible than those with H3K27ac or H3K4me3. Although overall chromatin accessibility may be defined by H3K27me3, the enrichment of H3K4me3 and H3K27ac and depletion of DNA methylation around transcriptional start sites up-regulated gene expression. Moreover, we found that subgenome C_n exhibited higher chromatin accessibility than A_n, which depended on the higher chromatin accessibility of C_n-unique genes but not homologous genes.

Changes in chromatin accessibility occuring during the process of allopolyploidization of Brassica napus are analysed using ATAC and ChIPseq, with differences in asymmetric chromatin accessibility between subgenomes of B. napus investigated.

Pangenomics in crop improvement-from coding structural variations to finding regulatory variants with pangenome graphs

Since the first reported crop pangenome in 2014, advances in high-throughput and cost-effective DNA sequencing technologies facilitated multiple such studies including the pangenomes of oilseed rape (Brassica napus L.), soybean [Glycine max (L.) Merr.], rice (Oryza sativa L.), wheat (Triticum aestivum L.), and barley (Hordeum vulgare L.). Compared with single-reference genomes, pangenomes provide a more accurate representation of the genetic variation present in a species. By combining the genomic data of multiple accessions, pangenomes allow for the detection and annotation of complex DNA polymorphisms such as structural variations (SVs), one of the major determinants of genetic diversity within a species. In this review we summarize the current literature on crop pangenomics, focusing on their application to find candidate SVs involved in traits of agronomic interest. We then highlight the potential of pangenomes in the discovery and functional characterization of noncoding regulatory sequences and their variations. We conclude with a summary and outlook on innovative data structures representing the complete content of plant pangenomes including annotations of coding and noncoding elements and outcomes of transcriptomic and epigenomic experiments.

The arches and spandrels of maize domestication, adaptation, and improvement

People living in the Balsas River basin in southwest México domesticated maize from the bushy grass teosinte. Nine thousand years later, in 2021, Ms. Deb Haaland - a member of the Pueblo of Laguna tribe of New Mexico - wore a dress adorned with a cornstalk when she was sworn in as the Secretary of Interior of the United States of America. This choice of garment highlights the importance of the coevolution of maize and the farmers who, through careful selection over thousands of years, domesticated maize and adapted the physiology and shoot architecture of maize to fit local environments and growth habits. Some traits such as tillering were directly selected on (arches), and others such as tassel size are the by-products (spandrels) of maize evolution. Here, we review current knowledge of the underlying cellular, developmental, physiological, and metabolic processes that were selected by farmers and breeders, which have positioned maize as a top global staple crop.

Intergenic spaces: a new frontier to improving plant health

To more sustainably mitigate the impact of crop diseases on plant health and productivity, there is a need for broader spectrum, long-lasting resistance traits. Defense response (DR) genes, located throughout the genome, participate in cellular and system-wide defense mechanisms to stave off infection by diverse pathogens. This multigenic resistance avoids rapid evolution of a pathogen to overcome host resistance. DR genes reside within resistance-associated quantitative trait loci (QTL), and alleles of DR genes in resistant varieties are more active during pathogen attack relative to susceptible haplotypes. Differential expression of DR genes results from polymorphisms in their regulatory regions, that includes cis-regulatory elements such as transcription factor binding sites as well as features that influence epigenetic structural changes to modulate chromatin accessibility during infection. Many of these elements are found in clusters, known as cis-regulatory modules (CRMs), which are distributed throughout the host genome. Regulatory regions involved in plant-pathogen interactions may also contain pathogen effector binding elements that regulate DR gene expression, and that, when mutated, result in a change in the plants' response. We posit that CRMs and the multiple regulatory elements that comprise them are potential targets for marker-assisted breeding for broad-spectrum, durable disease resistance.