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.

Switch them off or not: selective rRNA gene repression in grasses.

Nucleolar dominance (ND) is selective epigenetic silencing of 35-48S rDNA loci. In allopolyploids, it is frequently manifested at the cytogenetic level by the inactivation of nucleolar organiser region(s) (NORs) inherited from one or several evolutionary ancestors. Grasses are ecologically and economically one of the most important land plant groups, which have frequently evolved through hybridisation and polyploidisation events. Here we review common and unique features of ND phenomena in this monocot family from cytogenetic, molecular, and genomic perspectives. We highlight recent advances achieved by using an allotetraploid model grass, Brachypodium hybridum, where ND commonly occurs at a population level, and we cover modern genomic approaches that decipher structural features of core arrays of NORs.

Multiple origins, one evolutionary trajectory: gradual evolution characterizes distinct lineages of allotetraploid Brachypodium.

The "genomic shock" hypothesis posits that unusual challenges to genome integrity such as whole genome duplication may induce chaotic genome restructuring. Decades of research on polyploid genomes have revealed that this is often, but not always the case. While some polyploids show major chromosomal rearrangements and derepression of transposable elements in the immediate aftermath of whole genome duplication, others do not. Nonetheless, all polyploids show gradual diploidization over evolutionary time. To evaluate these hypotheses, we produced a chromosome-scale reference genome for the natural allotetraploid grass Brachypodium hybridum, accession "Bhyb26." We compared 2 independently derived accessions of B. hybridum and their deeply diverged diploid progenitor species Brachypodium stacei and Brachypodium distachyon. The 2 B. hybridum lineages provide a natural timecourse in genome evolution because one formed 1.4 million years ago, and the other formed 140 thousand years ago. The genome of the older lineage reveals signs of gradual post-whole genome duplication genome evolution including minor gene loss and genome rearrangement that are missing from the younger lineage. In neither B. hybridum lineage do we find signs of homeologous recombination or pronounced transposable element activation, though we find evidence supporting steady post-whole genome duplication transposable element activity in the older lineage. Gene loss in the older lineage was slightly biased toward 1 subgenome, but genome dominance was not observed at the transcriptomic level. We propose that relaxed selection, rather than an abrupt genomic shock, drives evolutionary novelty in B. hybridum, and that the progenitor species' similarity in transposable element load may account for the subtlety of the observed genome dominance.

Cytomolecular Organisation of the Nuclear Genome.

Modern molecular cytogenetics allows many aspects of the nuclear genome structure, function, and evolution to be analysed within the topographic context of mitotic and meiotic chromosomes and interphase nuclei [...].

The Chromosome Number and rDNA Loci Evolution in Onobrychis (Fabaceae).

The evolution of chromosome number and ribosomal DNA (rDNA) loci number and localisation were studied in Onobrychis Mill. Diploid and tetraploid species, as well as two basic chromosome numbers, x = 7 and x = 8, were observed among analysed taxa. The chromosomal distribution of rDNA loci was presented here for the first time using fluorescence in situ hybridisation (FISH) with 5S and 35S rDNA probes. Onobrychis species showed a high polymorphism in the number and localisation of rDNA loci among diploids, whereas the rDNA loci pattern was very similar in polyploids. Phylogenetic relationships among the species, inferred from nrITS sequences, were used as a framework to reconstruct the patterns of basic chromosome number and rDNA loci evolution. Analysis of the evolution of the basic chromosome numbers allowed the inference of x = 8 as the ancestral number and the descending dysploidy and polyploidisation as the major mechanisms of the chromosome number evolution. Analyses of chromosomal patterns of rRNA gene loci in a phylogenetic context resulted in the reconstruction of one locus of 5S rDNA and one locus of 35S rDNA in the interstitial chromosomal position as the ancestral state in this genus.

Brachypodium: 20 years as a grass biology model system; the way forward?

It has been 20 years since Brachypodium distachyon was suggested as a model grass species, but ongoing research now encompasses the entire genus. Extensive Brachypodium genome sequencing programmes have provided resources to explore the determinants and drivers of population diversity. This has been accompanied by cytomolecular studies to make Brachypodium a platform to investigate speciation, polyploidisation, perenniality, and various aspects of chromosome and interphase nucleus organisation. The value of Brachypodium as a functional genomic platform has been underscored by the identification of key genes for development, biotic and abiotic stress, and cell wall structure and function. While Brachypodium is relevant to the biofuel industry, its impact goes far beyond that as an intriguing model to study climate change and combinatorial stress.

Corrigendum: To Be or Not to Be Expressed: The First Evidence of a Nucleolar Dominance Tissue-Specificity in Brachypodium hybridum.

[This corrects the article DOI: 10.3389/fpls.2021.768347.].

Plant Cell and Organism Development 2.0.

In the editorial summarising the first edition of the Special Issue on "Plant Cell and Organism Development", we listed the key features that make plants a unique and fascinating group of living organisms [...].

Migration without interbreeding: Evolutionary history of a highly selfing Mediterranean grass inferred from whole genomes.

Wild plant populations show extensive genetic subdivision and are far from the ideal of panmixia which permeates population genetic theory. Understanding the spatial and temporal scale of population structure is therefore fundamental for empirical population genetics - and of interest in itself, as it yields insights into the history and biology of a species. In this study we extend the genomic resources for the wild Mediterranean grass Brachypodium distachyon to investigate the scale of population structure and its underlying history at whole-genome resolution. A total of 86 accessions were sampled at local and regional scales in Italy and France, which closes a conspicuous gap in the collection for this model organism. The analysis of 196 accessions, spanning the Mediterranean from Spain to Iraq, suggests that the interplay of high selfing and seed dispersal rates has shaped genetic structure in B. distachyon. At the continental scale, the evolution in B. distachyon is characterized by the independent expansion of three lineages during the Upper Pleistocene. Today, these lineages may occur on the same meadow yet do not interbreed. At the regional scale, dispersal and selfing interact and maintain high genotypic diversity, thus challenging the textbook notion that selfing in finite populations implies reduced diversity. Our study extends the population genomic resources for B. distachyon and suggests that an important use of this wild plant model is to investigate how selfing and dispersal, two processes typically studied separately, interact in colonizing plant species.

To Be or Not to Be Expressed: The First Evidence of a Nucleolar Dominance Tissue-Specificity in Brachypodium hybridum.

Nucleolar dominance (ND) is an epigenetic, developmentally regulated phenomenon that describes the selective inactivation of 35S rDNA loci derived from one progenitor of a hybrid or allopolyploid. The presence of ND was documented in an allotetraploid grass, Brachypodium hybridum (genome composition DDSS), which is a polyphyletic species that arose from crosses between two putative ancestors that resembled the modern B. distachyon (DD) and B. stacei (SS). In this work, we investigated the developmental stability of ND in B. hybridum genotype 3-7-2 and compared it with the reference genotype ABR113. We addressed the question of whether the ND is established in generative tissues such as pollen mother cells (PMC). We examined condensation of rDNA chromatin by fluorescence in situ hybridization employing state-of-art confocal microscopy. The transcription of rDNA homeologs was determined by reverse-transcription cleaved amplified polymorphic sequence analysis. In ABR113, the ND was stable in all tissues analyzed (primary and adventitious root, leaf, and spikes). In contrast, the 3-7-2 individuals showed a strong upregulation of the S-genome units in adventitious roots but not in other tissues. Microscopic analysis of the 3-7-2 PMCs revealed extensive decondensation of the D-genome loci and their association with the nucleolus in meiosis. As opposed, the S-genome loci were always highly condensed and localized outside the nucleolus. These results indicate that genotype-specific loss of ND in B. hybridum occurs probably after fertilization during developmental processes. This finding supports our view that B. hybridum is an attractive model to study ND in grasses.

Changes in the Cell Wall Proteome of Leaves in Response to High Temperature Stress in Brachypodium distachyon.

High temperature stress leads to complex changes to plant functionality, which affects, i.a., the cell wall structure and the cell wall protein composition. In this study, the qualitative and quantitative changes in the cell wall proteome of Brachypodium distachyon leaves in response to high (40 °C) temperature stress were characterised. Using a proteomic analysis, 1533 non-redundant proteins were identified from which 338 cell wall proteins were distinguished. At a high temperature, we identified 46 differentially abundant proteins, and of these, 4 were over-accumulated and 42 were under-accumulated. The most significant changes were observed in the proteins acting on the cell wall polysaccharides, specifically, 2 over- and 12 under-accumulated proteins. Based on the qualitative analysis, one cell wall protein was identified that was uniquely present at 40 °C but was absent in the control and 24 proteins that were present in the control but were absent at 40 °C. Overall, the changes in the cell wall proteome at 40 °C suggest a lower protease activity, lignification and an expansion of the cell wall. These results offer a new insight into the changes in the cell wall proteome in response to high temperature.

3,4-Dehydro-L-proline Induces Programmed Cell Death in the Roots of Brachypodium distachyon.

As cell wall proteins, the hydroxyproline-rich glycoproteins (HRGPs) take part in plant growth and various developmental processes. To fulfil their functions, HRGPs, extensins (EXTs) in particular, undergo the hydroxylation of proline by the prolyl-4-hydroxylases. The activity of these enzymes can be inhibited with 3,4-dehydro-L-proline (3,4-DHP), which enables its application to reveal the functions of the HRGPs. Thus, to study the involvement of HRGPs in the development of root hairs and roots, we treated seedlings of Brachypodium distachyon with 250 µM, 500 µM, and 750 µM of 3,4-DHP. The histological observations showed that the root epidermis cells and the cortex cells beneath them ruptured. The immunostaining experiments using the JIM20 antibody, which recognizes the EXT epitopes, demonstrated the higher abundance of this epitope in the control compared to the treated samples. The transmission electron microscopy analyses revealed morphological and ultrastructural features that are typical for the vacuolar-type of cell death. Using the TUNEL test (terminal deoxynucleotidyl transferase dUTP nick end labelling), we showed an increase in the number of nuclei with damaged DNA in the roots that had been treated with 3,4-DHP compared to the control. Finally, an analysis of two metacaspases' gene activity revealed an increase in their expression in the treated roots. Altogether, our results show that inhibiting the prolyl-4-hydroxylases with 3,4-DHP results in a vacuolar-type of cell death in roots, thereby highlighting the important role of HRGPs in root hair development and root growth.

Metabolomic Variation Aligns with Two Geographically Distinct Subpopulations of Brachypodium Distachyon before and after Drought Stress.

Brachypodium distachyon (Brachypodium) is a non-domesticated model grass that has been used to assess population level genomic variation. We have previously established a collection of 55 Brachypodium accessions that were sampled to reflect five different climatic regions of Turkey; designated 1a, 1c, 2, 3 and 4. Genomic and methylomic variation differentiated the collection into two subpopulations designated as coastal and central (respectively from regions 1a, 1c and the other from 2, 3 and 4) which were linked to environmental variables such as relative precipitation. Here, we assessed how far genomic variation would be reflected in the metabolomes and if this could be linked to an adaptive trait. Metabolites were extracted from eight-week-old seedlings from each accession and assessed using flow infusion high-resolution mass spectrometry (FIE-HRMS). Principal Component Analysis (PCA) of the derived metabolomes differentiated between samples from coastal and central subpopulations. The major sources of variation between seedling from the coastal and central subpopulations were identified. The central subpopulation was typified by significant increases in alanine, aspartate and glutamate metabolism and the tricarboxylic acid (TCA) cycle. Coastal subpopulation exhibited elevated levels of the auxin, indolacetic acid and rhamnose. The metabolomes of the seedling were also determined following the imposition of drought stress for seven days. The central subpopulation exhibited a metabolomic shift in response to drought, but no significant changes were seen in the coastal one. The drought responses in the central subpopulation were typified by changes in amino acids, increasing the glutamine that could be functioning as a stress signal. There were also changes in sugars that were likely to be an osmotic counter to drought, and changes in bioenergetic metabolism. These data indicate that genomic variation in our Turkish Brachypodium collection is largely reflected as distinctive metabolomes ("metabolotypes") through which drought tolerance might be mediated.

Allotetraploidization in Brachypodium May Have Led to the Dominance of One Parent's Metabolome in Germinating Seeds.

Seed germination is a complex process during which a mature seed resumes metabolic activity to prepare for seedling growth. In this study, we performed a comparative metabolomic analysis of the embryo and endosperm using the community standard lines of three annual Brachypodium species, i.e., B. distachyon (Bd) and B. stacei (Bs) and their natural allotetraploid B. hybridum (BdBs) that has wider ecological range than the other two species. We explored how far the metabolomic impact of allotetraploidization would be observable as over-lapping changes at 4, 12, and 24 h after imbibition (HAI) with water when germination was initiated. Metabolic changes during germination were more prominent in Brachypodium embryos than in the endosperm. The embryo and endosperm metabolomes of Bs and BdBs were similar, and those of Bd were distinctive. The Bs and BdBs embryos showed increased levels of sugars and the tricarboxylic acid cycle compared to Bd, which could have been indicative of better nutrient mobilization from the endosperm. Bs and BdBs also showed higher oxalate levels that could aid nutrient transfer through altered cellular events. In Brachypodium endosperm, the thick cell wall, in addition to starch, has been suggested to be a source of nutrients to the embryo. Metabolites indicative of sugar metabolism in the endosperm of all three species were not prominent, suggesting that mobilization mostly occurred prior to 4 HAI. Hydroxycinnamic and monolignol changes in Bs and BdBs were consistent with cell wall remodeling that arose following the release of nutrients to the respective embryos. Amino acid changes in both the embryo and endosperm were broadly consistent across the species. Taking our data together, the formation of BdBs may have maintained much of the Bs metabolome in both the embryo and endosperm during the early stages of germination. In the embryo, this conserved Bs metabolome appeared to include an elevated sugar metabolism that played a vital role in germination. If these observations are confirmed in the future with more Brachypodium accessions, it would substantiate the dominance of the Bs metabolome in BdBs allotetraploidization and the use of metabolomics to suggest important adaptive changes.

The formation and evolution of centromeric satellite repeats in Saccharum species.

Centromeres in eukaryotes are composed of tandem DNAs and retrotransposons. However, centromeric repeats exhibit considerable diversity, even among closely related species, and their origin and evolution are largely unknown. We conducted a genome-wide characterization of the centromeric sequences in sugarcane (Saccharum officinarum). Four centromeric tandem repeat sequences, So1, So103, So137 and So119, were isolated. So1 has a monomeric length of 137 bp, typical of a centromeric satellite, and has evolved four variants. However, these So1 variants had distinct centromere distributions and some were unique to an individual centromere. The distributions of the So1 variants were unexpectedly consistent among the Saccharum species that had different basic chromosome numbers or ploidy levels, thus suggesting evolutionary stability for approximately 7 million years in sugarcane. So103, So137 and So119 had unusually longer monomeric lengths that ranged from 327 to 1371 bp and lacked translational phasing on the CENH3 nucleosomes. Moreover, So103, So137 and So119 seemed to be highly similar to retrotransposons, which suggests that they originated from these mobile elements. Notably, all three repeats were flanked by direct repeats, and formed extrachromosomal circular DNAs (eccDNAs). The presence of circular molecules for these retrotransposon-derived centromeric satellites suggests an eccDNA-mediated centromeric satellite formation pathway in sugarcane.

Defining the Cell Wall, Cell Cycle and Chromatin Landmarks in the Responses of Brachypodium distachyon to Salinity.

Excess salinity is a major stress that limits crop yields. Here, we used the model grass Brachypodium distachyon (Brachypodium) reference line Bd21 in order to define the key molecular events in the responses to salt during germination. Salt was applied either throughout the germination period ("salt stress") or only after root emergence ("salt shock"). Germination was affected at ≥100 mM and root elongation at ≥75 mM NaCl. The expression of arabinogalactan proteins (AGPs), FLA1, FLA10, FLA11, AGP20 and AGP26, which regulate cell wall expansion (especially FLA11), were mostly induced by the "salt stress" but to a lesser extent by "salt shock". Cytological assessment using two AGP epitopes, JIM8 and JIM13 indicated that "salt stress" increases the fluorescence signals in rhizodermal and exodermal cell wall. Cell division was suppressed at >75 mM NaCl. The cell cycle genes (CDKB1, CDKB2, CYCA3, CYCB1, WEE1) were induced by "salt stress" in a concentration-dependent manner but not CDKA, CYCA and CYCLIN-D4-1-RELATED. Under "salt shock", the cell cycle genes were optimally expressed at 100 mM NaCl. These changes were consistent with the cell cycle arrest, possibly at the G1 phase. The salt-induced genomic damage was linked with the oxidative events via an increased glutathione accumulation. Histone acetylation and methylation and DNA methylation were visualized by immunofluorescence. Histone H4 acetylation at lysine 5 increased strongly whereas DNA methylation decreased with the application of salt. Taken together, we suggest that salt-induced oxidative stress causes genomic damage but that it also has epigenetic effects, which might modulate the cell cycle and AGP expression gene. Based on these landmarks, we aim to encourage functional genomics studies on the responses of Brachypodium to salt.

Genetic and Methylome Variation in Turkish Brachypodium Distachyon Accessions Differentiate Two Geographically Distinct Subpopulations.

Brachypodium distachyon (Brachypodium) is a non-domesticated model grass species that can be used to test if variation in genetic sequence or methylation are linked to environmental differences. To assess this, we collected seeds from 12 sites within five climatically distinct regions of Turkey. Seeds from each region were grown under standardized growth conditions in the UK to preserve methylated sequence variation. At six weeks following germination, leaves were sampled and assessed for genomic and DNA methylation variation. In a follow-up experiment, phenomic approaches were used to describe plant growth and drought responses. Genome sequencing and population structure analysis suggested three ancestral clusters across the Mediterranean, two of which were geographically separated in Turkey into coastal and central subpopulations. Phenotypic analyses showed that the coastal subpopulation tended to exhibit relatively delayed flowering and the central, increased drought tolerance as indicated by reduced yellowing. Genome-wide methylation analyses in GpC, CHG and CHH contexts also showed variation which aligned with the separation into coastal and central subpopulations. The climate niche modelling of both subpopulations showed a significant influence from the "Precipitation in the Driest Quarter" on the central subpopulation and "Temperature of the Coldest Month" on the coastal subpopulation. Our work demonstrates genetic diversity and variation in DNA methylation in Turkish accessions of Brachypodium that may be associated with climate variables and the molecular basis of which will feature in ongoing analyses.

Plant Cell and Organism Development.

Plants represent a unique and fascinating group of living organisms [...].