Does time matter? Intraspecific diversity of ribosomal RNA genes in lineages of the allopolyploid model grass Brachypodium hybridum with different evolutionary ages.

BACKGROUND: Polyploidisation often results in genome rearrangements that may involve changes in both the single-copy sequences and the repetitive genome fraction. In this study, we performed a comprehensive comparative analysis of repetitive DNA, with a particular focus on ribosomal DNA (rDNA), in Brachypodium hybridum (2n = 4x = 30, subgenome composition DDSS), an allotetraploid resulting from a natural cross between two diploid species that resemble the modern B. distachyon (2n = 10; DD) and B. stacei (2n = 20; SS). Taking advantage of the recurrent origin of B. hybridum, we investigated two genotypes, Bhyb26 and ABR113, differing markedly in their evolutionary age (1.4 and 0.14 Mya, respectively) and which resulted from opposite cross directions. To identify the origin of rDNA loci we employed cytogenetic and molecular methods (FISH, gCAPS and Southern hybridisation), phylogenetic and genomic approaches. RESULTS: Unlike the general maintenance of doubled gene dosage in B. hybridum, the rRNA genes showed a remarkable tendency towards diploidisation at both locus and unit levels. While the partial elimination of 35S rDNA units occurred in the younger ABR113 lineage, unidirectional elimination of the entire locus was observed in the older Bhyb26 lineage. Additionally, a novel 5S rDNA family was amplified in Bhyb26 replacing the parental units. The 35S and 5S rDNA units were preferentially eliminated from the S- and D-subgenome, respectively. Thus, in the more ancient B. hybridum lineage, Bhyb26, 5S and 35S rRNA genes are likely expressed from different subgenomes, highlighting the complexity of polyploid regulatory networks. CONCLUSION: Comparative analyses between two B. hybridum lineages of distinct evolutionary ages revealed that although the recent lineage ABR113 exhibited an additive pattern of rDNA loci distribution, the ancient lineage Bhyb26 demonstrated a pronounced tendency toward diploidisation manifested by the reduction in the number of both 35S and 5S loci. In conclusion, the age of the allopolyploid appears to be a decisive factor in rDNA turnover in B. hybridum.

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.

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.

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.

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 [...].

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.

The fate of 35S rRNA genes in the allotetraploid grass Brachypodium hybridum.

Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse-transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S-genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D-genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.

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.

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.

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.

Progressive refinement of the karyotyping of Brachypodium genomes.

Brachypodium distachyon is a weedy grass species that is firmly established as a model for the comparative and functional genomics of temperate cereals and grasses. Its simple, nuclear genome of five chromosomes contrasts it with other relatives of the genus with different, and usually higher, basic chromosome numbers and ploidy levels. This variation in karyotypic structure affords the possibility of reconstructing evolutionary pathways that have shaped the genome structure of extant species. This Tansley insight documents how key refinements in molecular cytogenetic approaches, from simple fluorescence in situ hybridization to comparative chromosome barcoding, have enabled genome structure studies and yielded valuable information about the drivers of karyotypic reorganization and evolution in the model grass genus Brachypodium.

Plant Cell and Organism Development.

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

From Single Cell to Plants: Mesophyll Protoplasts as a Versatile System for Investigating Plant Cell Reprogramming.

Plants are sessile organisms that have a remarkable developmental plasticity, which ensures their optimal adaptation to environmental stresses. Plant cell totipotency is an extreme example of such plasticity, whereby somatic cells have the potential to form plants via direct shoot organogenesis or somatic embryogenesis in response to various exogenous and/or endogenous signals. Protoplasts provide one of the most suitable systems for investigating molecular mechanisms of totipotency, because they are effectively single cell populations. In this review, we consider the current state of knowledge of the mechanisms that induce cell proliferation from individual, differentiated somatic plant cells. We highlight initial explant metabolic status, ploidy level and isolation procedure as determinants of successful cell reprogramming. We also discuss the importance of auxin signalling and its interaction with stress-regulated pathways in governing cell cycle induction and further stages of plant cell totipotency.

Comparative Genomics of Stenotrophomonas maltophilia and Stenotrophomonas rhizophila Revealed Characteristic Features of Both Species.

Although Stenotrophomonas maltophilia strains are efficient biocontrol agents, their field applications have raised concerns due to their possible threat to human health. The non-pathogenic Stenotrophomonas rhizophila species, which is closely related to S. maltophilia, has been proposed as an alternative. However, knowledge regarding the genetics of S. rhizophila is limited. Thus, the aim of the study was to define any genetic differences between the species and to characterise their ability to promote the growth of plant hosts as well as to enhance phytoremediation efficiency. We compared 37 strains that belong to both species using the tools of comparative genomics and identified 96 genetic features that are unique to S. maltophilia (e.g., chitin-binding protein, mechanosensitive channels of small conductance and KGG repeat-containing stress-induced protein) and 59 that are unique to S. rhizophila (e.g., glucosylglycerol-phosphate synthase, cold shock protein with the DUF1294 domain, and pteridine-dependent dioxygenase-like protein). The strains from both species have a high potential for biocontrol, which is mainly related to the production of keratinases (KerSMD and KerSMF), proteinases and chitinases. Plant growth promotion traits are attributed to the biosynthesis of siderophores, spermidine, osmoprotectants such as trehalose and glucosylglycerol, which is unique to S. rhizophila. In eight out of 37 analysed strains, the genes that are required to degrade protocatechuate were present. While our results show genetic differences between the two species, they had a similar growth promotion potential. Considering the information above, S. rhizophila constitutes a promising alternative for S. maltophilia for use in agricultural biotechnology.

In Vitro Tissue Culture in Brachypodium: Applications and Challenges.

Brachypodium distachyon has become an excellent model for plant breeding and bioenergy grasses that permits many fundamental questions in grass biology to be addressed. One of the constraints to performing research in many grasses has been the difficulty with which they can be genetically transformed and the generally low frequency of such transformations. In this review, we discuss the contribution that transformation techniques have made in Brachypodium biology as well as how Brachypodium could be used to determine the factors that might contribute to transformation efficiency. In particular, we highlight the latest research on the mechanisms that govern the gradual loss of embryogenic potential in a tissue culture and propose using B. distachyon as a model for other recalcitrant monocots.

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.

Centromeric DNA characterization in the model grass Brachypodium distachyon provides insights on the evolution of the genus.

Brachypodium distachyon is a well-established model monocot plant, and its small and compact genome has been used as an accurate reference for the much larger and often polyploid genomes of cereals such as Avena sativa (oats), Hordeum vulgare (barley) and Triticum aestivum (wheat). Centromeres are indispensable functional units of chromosomes and they play a core role in genome polyploidization events during evolution. As the Brachypodium genus contains about 20 species that differ significantly in terms of their basic chromosome numbers, genome size, ploidy levels and life strategies, studying their centromeres may provide important insight into the structure and evolution of the genome in this interesting and important genus. In this study, we isolated the centromeric DNA of the B. distachyon reference line Bd21 and characterized its composition via the chromatin immunoprecipitation of the nucleosomes that contain the centromere-specific histone CENH3. We revealed that the centromeres of Bd21 have the features of typical multicellular eukaryotic centromeres. Strikingly, these centromeres contain relatively few centromeric satellite DNAs; in particular, the centromere of chromosome 5 (Bd5) consists of only ~40 kb. Moreover, the centromeric retrotransposons in B. distachyon (CRBds) are evolutionarily young. These transposable elements are located both within and adjacent to the CENH3 binding domains, and have similar compositions. Moreover, based on the presence of CRBds in the centromeres, the species in this study can be grouped into two distinct lineages. This may provide new evidence regarding the phylogenetic relationships within the Brachypodium genus.

Ribosomal DNA loci derived from Brachypodium stacei are switched off for major parts of the life cycle of Brachypodium hybridum.

Nucleolar dominance is an epigenetic phenomenon that occurs in some plant and animal allopolyploids and hybrids, whereby only one ancestral set of 35S rRNA genes retains the ability to form the nucleolus while the rDNA loci derived from the other progenitor are transcriptionally silenced. There is substantial evidence that nucleolar dominance is regulated developmentally. This study focuses upon the establishment and/or maintenance of nucleolar dominance during different stages of development in the model grass allotetraploid Brachypodium hybridum. Fluorescence in situ hybridization with a 25S rDNA probe to cells in three-dimensional cytogenetic preparations showed that nucleolar dominance is present not only in root meristematic and differentiated cells of this species, but also in male meiocytes at prophase I, tetrads of microspores, and different embryonic tissues. The inactive state of Brachypodium stacei-originated rDNA loci was confirmed by silver staining. Only B. distachyon-derived 35S rDNA loci formed nucleoli in the aforementioned tissues, whereas B. stacei-like loci remained highly condensed and thus transcriptionally suppressed. The establishment of nucleolar dominance during earlier stages of B. hybridum embryo development cannot be ruled out. However, we propose that gradual pseudogenization of B. stacei-like loci in the evolution of the allotetraploid seems to be more likely.

Detecting Brachypodium distachyon Chromosomes Bd4 and Bd5 in MH- and X-Ray-Induced Micronuclei Using mcFISH.

Micronuclei are biomarkers of genotoxic effects and chromosomal instability. They are formed when chromosome fragments or whole chromosomes fail to disjoin into daughter nuclei. We present qualitative and quantitative analyses of the involvement of specific chromosome regions of chromosomes Bd4 and Bd5 in the formation of micronuclei of Brachypodium distachyon root tip cells following maleic hydrazide (MH) treatment and X-radiation. This is visualised by cytomolecular approaches using bacterial artificial chromosome (BAC)-based multicolour fluorescence in situ hybridisation (mcFISH) in combination with 5S and 25S rDNA probes. The results showed that the long arm of submetacentric chromosome Bd4 forms micronuclei at twice the frequency of its short arm, suggesting that the former is more prone to double-strand breaks (DSBs). In contrast, no difference was observed in the frequency of micronuclei derived from the long and short arms of submetacentric chromosome Bd5. Interestingly, the proximal region of the short arm of Bd5 is more prone to DSBs than its distal part. This demonstrates that 5S rDNA and 35S rDNA loci are not "hot spots" for DNA breaks after the application of these mutagens.