Gene Expression and Chromatin Modifications Associated with Maize Centromeres

A telomere-to-telomere cotton genome assembly reveals centromere evolution and a Mutator transposon-linked module regulating embryo development

Assembly of complete genomes can reveal functional genetic elements missing from draft sequences. Here we present the near-complete telomere-to-telomere and contiguous genome of the cotton species Gossypium raimondii. Our assembly identified gaps and misoriented or misassembled regions in previous assemblies and produced 13 centromeres, with 25 chromosomal ends having telomeres. In contrast to satellite-rich Arabidopsis and rice centromeres, cotton centromeres lack phased CENH3 nucleosome positioning patterns and probably evolved by invasion from long terminal repeat retrotransposons. In-depth expression profiling of transposable elements revealed a previously unannotated DNA transposon (MuTC01) that interacts with miR2947 to produce trans-acting small interfering RNAs (siRNAs), one of which targets the newly evolved LEC2 (LEC2b) to produce phased siRNAs. Systematic genome editing experiments revealed that this tripartite module, miR2947-MuTC01-LEC2b, controls the morphogenesis of complex folded embryos characteristic of Gossypium and its close relatives in the cotton tribe. Our study reveals a trans-acting siRNA-based tripartite regulatory pathway for embryo development in higher plants.

Characterization of the CsCENH3 protein and centromeric DNA profiles reveal the structures of centromeres in cucumber

Centromeres in eukaryotes mediate the accurate segregation of chromosomes during cell division. They serve as essential functional units of chromosomes and play a core role in the process of genome evolution. Centromeres are composed of satellite repeats and highly repetitive centromeric retrotransposons (CRs), which vary greatly even among closely related species. Cucumber (Cucumis sativus) is a globally cultivated and economically important vegetable and the only species in the Cucumis genus with seven pairs of chromosomes. Therefore, studying the centromeres of the Cucumis subgenus may yield valuable insights into its genome structure and evolution. Using chromatin immunoprecipitation (ChIP) techniques, we isolated centromeric DNA from cucumber reference line 9930. Our investigation into cucumber centromeres uncovered the centromeric satellite sequence, designated as CentCs, and the prevalence of Ty1/Copia long terminal repeat retrotransposons. In addition, active genes were identified in the CsCENH3 nucleosome regions with low transcription levels. To the best of our knowledge, this is the first time that characterization of centromeres has been achieved in cucumber. Meanwhile, our results on the distribution of CentCs and CsCRs in the subgenus Cucumis indicate that the content of centromeric repeats in the wild variants was significantly reduced compared with the cultivated cucumber. The results provide evidence for centromeric DNA amplification that occurred during the domestication process from wild to cultivated cucumber. Furthermore, these findings may offer new information for enhancing our understanding of phylogenetic relationships in the Cucumis genus.

The role of centromeric repeats and transcripts in kinetochore assembly and function

Centromeres are the chromosomal domains, where the kinetochore protein complex is formed, mediating proper segregation of chromosomes during cell division. Although the function of centromeres has remained conserved during evolution, centromeric DNA is highly variable, even in closely related species. In addition, the composition of the kinetochore complexes varies among organisms. Therefore, it is assumed that the centromeric position is determined epigenetically, and the centromeric histone H3 (CENH3) serves as an epigenetic marker. The loading of CENH3 onto centromeres depends on centromere-licensing factors, chaperones, and transcription of centromeric repeats. Several proteins that regulate CENH3 loading and kinetochore assembly interact with the centromeric transcripts and DNA in a sequence-independent manner. However, the functional aspects of these interactions are not fully understood. This review discusses the variability of centromeric sequences in different organisms and the regulation of their transcription through the RNA Pol II and RNAi machinery. The data suggest that the interaction of proteins involved in CENH3 loading and kinetochore assembly with centromeric DNA and transcripts plays a role in centromere, and possibly neocentromere, formation in a sequence-independent manner.

The transcribed centromeric gene OsMRPL15 is essential for pollen development in rice

Centromeres consist of highly repetitive sequences that are challenging to map, clone, and sequence. Active genes exist in centromeric regions, but their biological functions are difficult to explore owing to extreme suppression of recombination in these regions. In this study, we used the CRISPR/Cas9 system to knock out the transcribed gene Mitochondrial Ribosomal Protein L15 (OsMRPL15), located in the centromeric region of rice (Oryza sativa) chromosome 8, resulting in gametophyte sterility. Osmrpl15 pollen was completely sterile, with abnormalities appearing at the tricellular stage including the absence of starch granules and disrupted mitochondrial structure. Loss of OsMRPL15 caused abnormal accumulation of mitoribosomal proteins and large subunit rRNA in pollen mitochondria. Moreover, the biosynthesis of several proteins in mitochondria was defective, and expression of mitochondrial genes was upregulated at the mRNA level. Osmrpl15 pollen contained smaller amounts of intermediates related to starch metabolism than wild-type pollen, while biosynthesis of several amino acids was upregulated, possibly to compensate for defective mitochondrial protein biosynthesis and initiate consumption of carbohydrates necessary for starch biosynthesis. These results provide further insight into how defects in mitoribosome development cause gametophyte male sterility.

Characterization of centromeric DNA of Gossypium anomalum reveals sequence-independent enrichment dynamics of centromeric repeats

Centromeres in eukaryotes are composed of highly repetitive DNAs, which evolve rapidly and are thought to achieve a favorable structure in mature centromeres. However, how the centromeric repeat evolves into an adaptive structure is largely unknown. We characterized the centromeric sequences of Gossypium anomalum through chromatin immunoprecipitation against CENH3 antibodies. We revealed that the G. anomalum centromeres contained only retrotransposon-like repeats but were depleted in long arrays of satellites. These retrotransposon-like centromeric repeats were present in the African-Asian and Australian lineage species, suggesting that they might have arisen in the common ancestor of these diploid species. Intriguingly, we observed a substantial increase and decrease in copy numbers among African-Asian and Australian lineages, respectively, for the retrotransposon-derived centromeric repeats without apparent structure or sequence variation in cotton. This result indicates that the sequence content is not a decisive aspect of the adaptive evolution of centromeric repeats or at least retrotransposon-like centromeric repeats. In addition, two active genes with potential roles in gametogenesis or flowering were identified in CENH3 nucleosome-binding regions. Our results provide new insights into the constitution of centromeric repetitive DNA and the adaptive evolution of centromeric repeats in plants.

An epigenetic basis of inbreeding depression in maize

Inbreeding depression is widespread across plant and animal kingdoms and may arise from the exposure of deleterious alleles and/or loss of overdominant alleles resulting from increased homozygosity, but these genetic models cannot fully explain the phenomenon. Here, we report epigenetic links to inbreeding depression in maize. Teosinte branched1/cycloidea/proliferating cell factor (TCP) transcription factors control plant development. During successive inbreeding among inbred lines, thousands of genomic regions across TCP-binding sites (TBS) are hypermethylated through the H3K9me2-mediated pathway. These hypermethylated regions are accompanied by decreased chromatin accessibility, increased levels of the repressive histone marks H3K27me2 and H3K27me3, and reduced binding affinity of maize TCP-proteins to TBS. Consequently, hundreds of TCP-target genes involved in mitochondrion, chloroplast, and ribosome functions are down-regulated, leading to reduced growth vigor. Conversely, random mating can reverse corresponding hypermethylation sites and TCP-target gene expression, restoring growth vigor. These results support a unique role of reversible epigenetic modifications in inbreeding depression.

Genome-wide mapping reveals R-loops associated with centromeric repeats in maize

R-loops are stable chromatin structures comprising a DNA:RNA hybrid and a displaced single-stranded DNA. R-loops have been implicated in gene expression and chromatin structure, as well as in replication blocks and genome instability. Here, we conducted a genome-wide identification of R-loops and identified more than 700,000 R-loop peaks in the maize (Zea mays) genome. We found that sense R-loops were mainly enriched in promoters and transcription termination sites and relatively less enriched in gene bodies, which is different from the main gene-body localization of sense R-loops in Arabidopsis and Oryza sativa. At the chromosome scale, maize R-loops were enriched in pericentromeric heterochromatin regions, and a significant portion of R-loops were derived from transposable elements. In centromeres, R-loops preferentially formed within the binding regions of the centromere-specific histone CENH3, and centromeric retrotransposons were strongly associated with R-loop formation. Furthermore, centromeric retrotransposon R-loops were observed by applying the single-molecule imaging technique of atomic force microscopy. These findings elucidate the fundamental character of R-loops in the maize genome and reveal the potential role of R-loops in centromeres.

Comparing DNA replication programs reveals large timing shifts at centromeres of endocycling cells in maize roots

Plant cells undergo two types of cell cycles–the mitotic cycle in which DNA replication is coupled to mitosis, and the endocycle in which DNA replication occurs in the absence of cell division. To investigate DNA replication programs in these two types of cell cycles, we pulse labeled intact root tips of maize (Zea mays) with 5-ethynyl-2’-deoxyuridine (EdU) and used flow sorting of nuclei to examine DNA replication timing (RT) during the transition from a mitotic cycle to an endocycle. Comparison of the sequence-based RT profiles showed that most regions of the maize genome replicate at the same time during S phase in mitotic and endocycling cells, despite the need to replicate twice as much DNA in the endocycle and the fact that endocycling is typically associated with cell differentiation. However, regions collectively corresponding to 2% of the genome displayed significant changes in timing between the two types of cell cycles. The majority of these regions are small with a median size of 135 kb, shift to a later RT in the endocycle, and are enriched for genes expressed in the root tip. We found larger regions that shifted RT in centromeres of seven of the ten maize chromosomes. These regions covered the majority of the previously defined functional centromere, which ranged between 1 and 2 Mb in size in the reference genome. They replicate mainly during mid S phase in mitotic cells but primarily in late S phase of the endocycle. In contrast, the immediately adjacent pericentromere sequences are primarily late replicating in both cell cycles. Analysis of CENH3 enrichment levels in 8C vs 2C nuclei suggested that there is only a partial replacement of CENH3 nucleosomes after endocycle replication is complete. The shift to later replication of centromeres and possible reduction in CENH3 enrichment after endocycle replication is consistent with a hypothesis that centromeres are inactivated when their function is no longer needed.

In traditional cell division, or mitosis, a cell’s genetic material is duplicated and then split between two daughter cells. In contrast, in some specialized cell types, the DNA is duplicated a second time without an intervening division step, resulting in cells that carry twice as much DNA. This phenomenon, which is called the endocycle, is common during plant development. At each step, DNA replication follows an ordered program in which highly compacted DNA is unraveled and replicated in sections at different times during the synthesis (S) phase. In plants, it is unclear whether traditional and endocycle programs are the same, especially since endocycling cells are typically in the process of differentiation. Using root tips of maize, we found that in comparison to replication in the mitotic cell cycle, there is a small portion of the genome whose replication in the endocycle is shifted in time, usually to later in S phase. Some of these regions are scattered around the genome and mostly coincide with active genes. However, the most prominent shifts occur in centromeres. The shift to later replication in centromeres is noteworthy because they orchestrate the process of separating duplicated chromosomes into daughter cells, a function that is not needed in the endocycle.

Meiotic Nuclear Architecture in Distinct Mole Vole Hybrids with Robertsonian Translocations: Chromosome Chains, Stretched Centromeres, and Distorted Recombination

Genome functioning in hybrids faces inconsistency. This mismatch is manifested clearly in meiosis during chromosome synapsis and recombination. Species with chromosomal variability can be a model for exploring genomic battles with high visibility due to the use of advanced immunocytochemical methods. We studied synaptonemal complexes (SC) and prophase I processes in 44-chromosome intraspecific (Ellobius tancrei × E. tancrei) and interspecific (Ellobius talpinus × E. tancrei) hybrid mole voles heterozygous for 10 Robertsonian translocations. The same pachytene failures were found for both types of hybrids. In the intraspecific hybrid, the chains were visible in the pachytene stage, then 10 closed SC trivalents formed in the late pachytene and diplotene stage. In the interspecific hybrid, as a rule, SC trivalents composed the SC chains and rarely could form closed configurations. Metacentrics involved with SC trivalents had stretched centromeres in interspecific hybrids. Linkage between neighboring SC trivalents was maintained by stretched centromeric regions of acrocentrics. This centromeric plasticity in structure and dynamics of SC trivalents was found for the first time. We assume that stretched centromeres were a marker of altered nuclear architecture in heterozygotes due to differences in the ancestral chromosomal territories of the parental species. Restructuring of the intranuclear organization and meiotic disturbances can contribute to the sterility of interspecific hybrids, and lead to the reproductive isolation of studied species.

Meiotic crossovers characterized by haplotype-specific chromosome painting in maize

Meiotic crossovers (COs) play a critical role in generating genetic variation and maintaining faithful segregation of homologous chromosomes during meiosis. We develop a haplotype-specific fluorescence in situ hybridization (FISH) technique that allows visualization of COs directly on metaphase chromosomes. Oligonucleotides (oligos) specific to chromosome 10 of maize inbreds B73 and Mo17, respectively, are synthesized and labeled as FISH probes. The parental and recombinant chromosome 10 in B73 x Mo17 F₁ hybrids and F₂ progenies can be unambiguously identified by haplotype-specific FISH. Analysis of 58 F₂ plants reveals lack of COs in the entire proximal half of chromosome 10. However, we detect COs located in regions very close to the centromere in recombinant inbred lines from an intermated B73 x Mo17 population, suggesting effective accumulation of COs in recombination-suppressed chromosomal regions through intermating and the potential to generate favorable allelic combinations of genes residing in these regions.

Meiotic crossovers (COs) are essential for proper chromosome segregation and generating novel combinations of alleles. Here, the authors develop haplotype-specific oligos on maize chromosome 10 for fluorescence in situ hybridization and analyze CO patterns in an intermated recombinant population derived from B73 and Mo17.

Transcriptional and epigenetic adaptation of maize chromosomes in Oat-Maize addition lines

By putting heterologous genomic regulatory systems into contact, chromosome addition lines derived from interspecific or intergeneric crosses allow the investigation of transcriptional regulation in new genomic environments. Here, we report the transcriptional and epigenetic adaptation of stably inherited alien maize chromosomes in two oat–maize addition (OMA) lines. We found that the majority of maize genes displayed maize-specific transcription in the oat genomic environment. Nevertheless, a quarter of the expressed genes encoded by the two maize chromosomes were differentially expressed genes (DEGs). Notably, highly conserved orthologs were more severely differentially expressed in OMAs than less conserved orthologs. Additionally, syntenic genes and highly abundant genes were over-represented among DEGs. Gene suppression was more common than activation among the DEGs; however, the genes in the former maize pericentromere, which expanded to become the new centromere in OMAs, were activated. Histone modifications (H3K4me3, H3K9ac and H3K27me3) were consistent with these transcriptome results. We expect that cis regulation is responsible for unchanged expression in OMA versus maize; and trans regulation is the predominant mechanism behind DEGs. The genome interaction identified here reveals the important consequences of interspecific/intergeneric crosses and potential mechanisms of plant evolution when genomic environments interact.

Comparison of Oryza sativa and Oryza brachyantha Genomes Reveals Selection-Driven Gene Escape from the Centromeric Regions

Centromeres are dynamic chromosomal regions, and the genetic and epigenetic environment of the centromere is often regarded as oppressive to protein-coding genes. Here, we used comparative genomic and phylogenomic approaches to study the evolution of centromeres and centromere-linked genes in the genus Oryza We report a 12.4-Mb high-quality BAC-based pericentromeric assembly for Oryza brachyantha, which diverged from cultivated rice (Oryza sativa) ∼15 million years ago. The synteny analyses reveal seven medium (>50 kb) pericentric inversions in O. sativa and 10 in O. brachyantha Of these inversions, three resulted in centromere movement (Chr1, Chr7, and Chr9). Additionally, we identified a potential centromere-repositioning event, in which the ancestral centromere on chromosome 12 in O. brachyantha jumped ∼400 kb away, possibly mediated by a duplicated transposition event (>28 kb). More strikingly, we observed an excess of syntenic gene loss at and near the centromeric regions (P < 2.2 × 10^-16). Most (33/47) of the missing genes moved to other genomic regions; therefore such excess could be explained by the selective loss of the copy in or near centromeric regions after gene duplication. The pattern of gene loss immediately adjacent to centromeric regions suggests centromere chromatin dynamics (e.g., spreading or microrepositioning) may drive such gene loss.

Proliferation of Regulatory DNA Elements Derived from Transposable Elements in the Maize Genome

Genomic regions free of nucleosomes, which are hypersensitive to DNase I digestion, are known as DNase I hypersensitive sites (DHSs) and frequently contain cis-regulatory DNA elements. To investigate their prevalence and characteristics in maize (Zea mays), we developed high-resolution genome-wide DHS maps using a modified DNase-seq technique. Maize DHSs exhibit depletion of nucleosomes and low levels of DNA methylation and are enriched with conserved noncoding sequences (CNSs). We developed a protoplast-based transient transformation assay to assess the potential gene expression enhancer and/or promoter functions associated with DHSs, which showed that more than 80% of DHSs overlapping with CNSs showed an enhancer function. Strikingly, nearly 25% of maize DHSs were derived from transposable elements (TEs), including both class I and class II transposons. Interestingly, TE-derived DHSs (teDHSs) homologous to retrotransposons were enriched with sequences related to the intrinsic cis-regulatory elements within the long terminal repeats of retrotransposons. We demonstrate that more than 80% of teDHSs can drive transcription of a reporter gene in protoplast assays. These results reveal the widespread occurrence of TE-derived cis-regulatory sequences and suggest that teDHSs play a major role in transcriptional regulation in maize.

Genomic Analysis of the DNA Replication Timing Program during Mitotic S Phase in Maize (Zea mays) Root Tips

All plants and animals must replicate their DNA, using a regulated process to ensure that their genomes are completely and accurately replicated. DNA replication timing programs have been extensively studied in yeast and animal systems, but much less is known about the replication programs of plants. We report a novel adaptation of the "Repli-seq" assay for use in intact root tips of maize (Zea mays) that includes several different cell lineages and present whole-genome replication timing profiles from cells in early, mid, and late S phase of the mitotic cell cycle. Maize root tips have a complex replication timing program, including regions of distinct early, mid, and late S replication that each constitute between 20 and 24% of the genome, as well as other loci corresponding to ∼32% of the genome that exhibit replication activity in two different time windows. Analyses of genomic, transcriptional, and chromatin features of the euchromatic portion of the maize genome provide evidence for a gradient of early replicating, open chromatin that transitions gradually to less open and less transcriptionally active chromatin replicating in mid S phase. Our genomic level analysis also demonstrated that the centromere core replicates in mid S, before heavily compacted classical heterochromatin, including pericentromeres and knobs, which replicate during late S phase.

Recurrent establishment of de novo centromeres in the pericentromeric region of maize chromosome 3

Centromeres can arise de novo from non-centromeric regions, which are often called "neocentromeres." Neocentromere formation provides the best evidence for the concept that centromere function is not determined by the underlying DNA sequences, but controlled by poorly understood epigenetic mechanisms. Numerous neocentromeres have been reported in several plant and animal species. However, it has been elusive how and why a specific chromosomal region is chosen to be a new centromere during the neocentromere activation events. We report recurrent establishment of neocentromeres in a pericentromeric region of chromosome 3 in maize (Zea mays). This latent region is located in the short arm and is only 2 Mb away from the centromere (Cen3) of chromosome 3. At least three independent neocentromere activation events, which were likely induced by different mechanisms, occurred within this latent region. We mapped the binding domains of CENH3, the centromere-specific H3 histone variant, of the three neocentromeres and analyzed the genomic and epigenomic features associated with Cen3, the de novo centromeres and an inactivated centromere derived from an ancestral chromosome. Our results indicate that lack of genes and transcription and a relatively high level of DNA methylation in this pericentromeric region may provide a favorable chromatin environment for neocentromere activation.

Repliscan: a tool for classifying replication timing regions

BACKGROUND

Replication timing experiments that use label incorporation and high throughput sequencing produce peaked data similar to ChIP-Seq experiments. However, the differences in experimental design, coverage density, and possible results make traditional ChIP-Seq analysis methods inappropriate for use with replication timing.

RESULTS

To accurately detect and classify regions of replication across the genome, we present Repliscan. Repliscan robustly normalizes, automatically removes outlying and uninformative data points, and classifies Repli-seq signals into discrete combinations of replication signatures. The quality control steps and self-fitting methods make Repliscan generally applicable and more robust than previous methods that classify regions based on thresholds.

CONCLUSIONS

Repliscan is simple and effective to use on organisms with different genome sizes. Even with analysis window sizes as small as 1 kilobase, reliable profiles can be generated with as little as 2.4x coverage.

​Plant centromeres​

Plant centromeres, which are determined epigenetically by centromeric histone 3 (CENH3) have revealed surprising structural diversity, ranging from the canonical monocentric seen in vertebrates, to polycentric, and holocentric. Normally stable, centromeres can change position over evolutionary times or upon genomic stress, such as when chromosomes are broken. At the DNA level, centromeres can be based on single copy DNA or more commonly on repeats. Rapid evolution of centromeric sequences and of CENH3 protein remains a mystery, as evidence of co-adaptation is lacking. Epigenetic differences between parents can trigger uniparental centromere failure and genome elimination, contributing to postzygotic hybridization barriers.​.

[Reciprocity between active transcription and histone methylation]

In the nucleus of eukaryotic cells, the chromatin states dictated by the different combinations of histone post-translational modifications, such as the methylation of lysine residues, are an integral part of the multitude of epigenomes involved in the fine tuning of all genome functions, and in particular transcription. Over the last decade, an increasing number of factors have been identified as regulators involved in the establishment, reading or erasure of histone methylations. Their characterization in model organisms such as Arabidopsis has thus unraveled their fundamental roles in the control and regulation of essential developmental processes such as the floral transition, cell differentiation, gametogenesis, and/or the response/adaptation of plants to environmental stresses. In this review, we will focus on the methylation of histones functioning as a mark of activate transcription and we will try to highlight, based on recent findings, the more or less direct links between this mark and gene expression. Thus, we will discuss the different mechanisms allowing the dynamics and the integration of the chromatin states resulting from the different histone methylations in connection with the transcriptional machinery of the RNA polymerase II.

Isolation and characterization of centromeric repetitive DNA sequences in Saccharum spontaneum

Sugarcane (Saccharum hybrids spp.) is the most important sugar crop that accounts for ~75% of the world’s sugar production. Recently, a whole-genome sequencing project was launched on the wild species S. spontaneum. To obtain information on the DNA composition of the repeat-enriched region of the centromere, we conducted a genome-wide analysis of the DNA sequences associated with CenH3 (a mutant of histone H3 located in eukaryote centromeres) using chromatin immunoprecipitation followed by sequencing (ChIP-seq) method. We demonstrate that the centromeres contain mainly SCEN-like single satellite repeat (Ss1) and several Ty3/gypsy retrotransposon-related repeats (Ss166, Ss51, and Ss68). Ss1 dominates in the centromeric regions and spans up to 500 kb. In contrast, the Ty3/gypsy retrotransposon-related repeats are either clustered spanning over a short range, or dispersed in the centromere regions. Interestingly, Ss1 exhibits a chromosome-specific enrichment in the wild species S. spontaneum and S. robustum, but not in the domesticated species S. officinarum and modern sugarcane cultivars. This finding suggests an autopolyploid genome identity of S. spontaneum with a high level of homology among its eight sub-genomes. We also conducted a genome-wide survey of the repetitive DNAs in S. spontaneum following a similarity-based sequence clustering strategy. These results provide insight into the composition of sugarcane genome as well as the genome assembly of S. spontaneum.

Dynamic chromatin changes associated with de novo centromere formation in maize euchromatin

The inheritance and function of centromeres are not strictly dependent on any specific DNA sequence, but involve an epigenetic component in most species. CENH3, a centromere histone H3 variant, is one of the best-described epigenetic factors in centromere identity, but the chromatin features required during centromere formation have not yet been revealed. We previously identified two de novo centromeres on Zea mays (maize) minichromosomes derived from euchromatic sites with high-density gene distributions but low-density transposon distributions. The distribution of gene location and gene expression in these sites indicates that transcriptionally active regions can initiate de novo centromere formation, and CENH3 seeding shows a preference for gene-free regions or regions with no gene expression. The locations of the expressed genes detected were at relatively hypomethylated loci, and the altered gene expression resulted from de novo centromere formation, but not from the additional copy of the minichromosome. The initial overall DNA methylation level of the two de novo regions was at a low level, but increased substantially to that of native centromeres after centromere formation. These results illustrate the dynamic chromatin changes during euchromatin-originated de novo centromere formation, which provides insight into the mechanism of de novo centromere formation and regulation of subsequent consequences.