Subgenomic Stability of Progenitor Genomes During Repeated Allotetraploid Origins of the Same Grass Brachypodium hybridum

Both homeologous exchanges and homeologous expression bias are generally found in most allopolyploid species. Whether homeologous exchanges and homeologous expression bias differ between repeated allopolyploid speciation events from the same progenitor species remains unknown. Here, we detected a third independent and recent allotetraploid origin for the model grass Brachypodium hybridum. Our homeologous exchange with replacement analyses indicated the absence of significant homeologous exchanges in any of the three types of wild allotetraploids, supporting the integrity of their progenitor subgenomes and the immediate creation of the amphidiploids. Further homeologous expression bias tests did not uncover significant subgenomic dominance in different tissues and conditions of the allotetraploids. This suggests a balanced expression of homeologs under similar or dissimilar ecological conditions in their natural habitats. We observed that the density of transposons around genes was not associated with the initial establishment of subgenome dominance; rather, this feature is inherited from the progenitor genome. We found that drought response genes were highly induced in the two subgenomes, likely contributing to the local adaptation of this species to arid habitats in the third allotetraploid event. These findings provide evidence for the consistency of subgenomic stability of parental genomes across multiple allopolyploidization events that led to the same species at different periods. Our study emphasizes the importance of selecting closely related progenitor species genomes to accurately assess homeologous exchange with replacement in allopolyploids, thereby avoiding the detection of false homeologous exchanges when using less related progenitor species genomes.

Cell type-specific cytonuclear coevolution in three allopolyploid plant species

Cytonuclear disruption may accompany allopolyploid evolution as a consequence of the merger of different nuclear genomes in a cellular environment having only one set of progenitor organellar genomes. One path to reconcile potential cytonuclear mismatch is biased expression for maternal gene duplicates (homoeologs) encoding proteins that target to plastids and/or mitochondria. Assessment of this transcriptional form of cytonuclear coevolution at the level of individual cells or cell types remains unexplored. Using single-cell (sc-) and single-nucleus (sn-) RNAseq data from eight tissues in three allopolyploid species, we characterized cell type-specific variations of cytonuclear coevolutionary homoeologous expression and demonstrated the temporal dynamics of expression patterns across development stages during cotton fiber development. Our results provide unique insights into transcriptional cytonuclear coevolution in plant allopolyploids at the single-cell level.

Protein nonadditive expression and solubility contribute to heterosis in Arabidopsis hybrids and allotetraploids

Hybrid vigor or heterosis has been widely applied in agriculture and extensively studied using genetic and gene expression approaches. However, the biochemical mechanism underlying heterosis remains elusive. One theory suggests that a decrease in protein aggregation may occur in hybrids due to the presence of protein variants between parental alleles, but it has not been experimentally tested. Here, we report comparative analysis of soluble and insoluble proteomes in Arabidopsis intraspecific and interspecific hybrids or allotetraploids formed between A. thaliana and A. arenosa. Both allotetraploids and intraspecific hybrids displayed nonadditive expression (unequal to the sum of the two parents) of the proteins, most of which were involved in biotic and abiotic stress responses. In the allotetraploids, homoeolog-expression bias was not observed among all proteins examined but accounted for 17-20% of the nonadditively expressed proteins, consistent with the transcriptome results. Among expression-biased homoeologs, there were more A. thaliana-biased than A. arenosa-biased homoeologs. Analysis of the insoluble and soluble proteomes revealed more soluble proteins in the hybrids than their parents but not in the allotetraploids. Most proteins in ribosomal biosynthesis and in the thylakoid lumen, membrane, and stroma were in the soluble fractions, indicating a role of protein stability in photosynthetic activities for promoting growth. Thus, nonadditive expression of stress-responsive proteins and increased solubility of photosynthetic proteins may contribute to heterosis in Arabidopsis hybrids and allotetraploids and possibly hybrid crops.

Sequence capture: Obsolete or irreplaceable? A thorough validation across phylogenetic distances and its applicability to hybrids and allopolyploids

As whole-genome sequencing has become pervasive, some have suggested that reduced genomic representation approaches, for example, sequence capture, are becoming obsolete. In the present study, we argue that these techniques still provide excellent tools in terms of price and quality of data as well as in their ability to provide markers with specific features, as required, for example, in phylogenomics. A potential drawback of the wide-scale application of reduced representation approaches could be their drop in efficiency with increasing phylogenetic distance from the reference species. While some studies have focused on the degree and performance of reduced representation techniques in such situations, to our knowledge, none of them evaluated their applicability to inter-specific hybrids and polyploids. This highlights a significant gap in current knowledge since there is increasing evidence for the frequent occurrence of natural hybrids and polyploids, as well as for the major importance of both phenomena in evolution. The main aim of the present study was to carry out a thorough validation of SEQcap applicability to (1) a set of non-model taxa with a wide range of phylogenetic relatedness and (2) inter-specific hybrids of various ploidies and genomic compositions. Considering the latter point, we especially focused on mechanisms causing allelic bias and consequent allelic dropout, as these could have confounding effects with respect to the evolutionary genomic dynamics of hybrids, especially in asexuals, which virtually reproduce as a frozen F1 generation.

A temporal gradient of cytonuclear coordination of chaperonins and chaperones during RuBisCo biogenesis in allopolyploid plants

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCo), consisting of subunits encoded by nuclear and cytoplasmic genes, is a model for cytonuclear evolution in plant allopolyploids. To date, coordinated cytonuclear evolutionary responses of auxiliary cofactors involved in RuBisCo biogenesis remain unexplored. This study characterized and compared genomic and transcriptional cytonuclear coevolutionary responses of chaperonin/chaperones in RuBisCo folding and assembly processes across different allopolyploids. We discovered significant cytonuclear evolutionary responses in folding cofactors, with diminishing or attenuated responses later during assembly. Our results have general significance for understanding the unrecognized cytonuclear evolution of chaperonin/chaperone genes, structural and functional features of intermediate complexes, and the functioning stage of the Raf2 cofactor. Generally, the results reveal a hitherto unexplored dimension of allopolyploidy in plants.

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCo) has long been studied from many perspectives. As a multisubunit (large subunits [LSUs] and small subunits[SSUs]) protein encoded by genes residing in the chloroplast (rbcL) and nuclear (rbcS) genomes, RuBisCo also is a model for cytonuclear coevolution following allopolyploid speciation in plants. Here, we studied the genomic and transcriptional cytonuclear coordination of auxiliary chaperonin and chaperones that facilitate RuBisCo biogenesis across multiple natural and artificially synthesized plant allopolyploids. We found similar genomic and transcriptional cytonuclear responses, including respective paternal-to-maternal conversions and maternal homeologous biased expression, in chaperonin/chaperon-assisted folding and assembly of RuBisCo in different allopolyploids. One observation is about the temporally attenuated genomic and transcriptional cytonuclear evolutionary responses during early folding and later assembly process of RuBisCo biogenesis, which were established by long-term evolution and immediate onset of allopolyploidy, respectively. Our study not only points to the potential widespread and hitherto unrecognized features of cytonuclear evolution but also bears implications for the structural interaction interface between LSU and Cpn60 chaperonin and the functioning stage of the Raf2 chaperone.

SubPhaser: a robust allopolyploid subgenome phasing method based on subgenome-specific k-mers

With advanced sequencing technology, dozens of complex polyploid plant genomes have been characterized. However, for many polyploid species, their diploid ancestors are unknown or extinct, making it impossible to unravel the subgenomes and genome evolution directly. We developed a novel subgenome-phasing algorithm, SubPhaser, specifically designed for a neoallopolyploid or a homoploid hybrid. SubPhaser first searches for the subgenome-specific sequence (k-mer), then assigns homoeologous chromosomes into subgenomes, and further provides tools to annotate and investigate specific sequences. SubPhaser works well on neoallopolyploids and homoploid hybrids containing subgenome-specific sequences like wheat, but fails on autopolyploids lacking subgenome-specific sequences like alfalfa, indicating that SubPhaser can phase neoallopolyploid/homoploid hybrids with high accuracy, sensitivity and performance. This highly accurate, highly sensitive, ancestral data free chromosome phasing algorithm, SubPhaser, offers significant application value for subgenome phasing in neoallopolyploids and homoploid hybrids, and for the subsequent exploration of genome evolution and related genetic/epigenetic mechanisms.

Chromosome-specific painting unveils chromosomal fusions and distinct allopolyploid species in the Saccharum complex

Karyotypes provide key cytogenetic information on the phylogenetic relationships and evolutionary origins in related eukaryotic species. Despite our knowledge of the chromosome numbers of sugarcane and its wild relatives, the chromosome composition and evolution among the species in the Saccharum complex have been elusive owing to the complex polyploidy and the large numbers of chromosomes of these species. Oligonucleotide-based chromosome painting has become a powerful tool of cytogenetic studies especially for plant species with large numbers of chromosomes. We developed oligo-based chromosome painting probes for all 10 chromosomes in Saccharum officinarum (2n = 8x = 80). The 10 painting probes generated robust fluorescence in situ hybridization signals in all plant species within the Saccharum complex, including species in the genera Saccharum, Miscanthus, Narenga and Erianthus. We conducted comparative chromosome analysis using the same set of probes among species from four different genera within the Saccharum complex. Excitingly, we discovered several novel cytotypes and chromosome rearrangements in these species. We discovered that fusion from two different chromosomes is a common type of chromosome rearrangement associated with the species in the Saccharum complex. Such fusion events changed the basic chromosome number and resulted in distinct allopolyploids in the Saccharum complex.

Chromosome-Level Genome Assembly Provides New Insights into Genome Evolution and Tuberous Root Formation of Potentilla anserina

Potentilla anserina is a perennial stoloniferous plant with edible tuberous roots in Rosaceae, served as important food and medicine sources for Tibetans in the Qinghai-Tibetan Plateau (QTP), China, over thousands of years. However, a lack of genome information hindered the genetic study. Here, we presented a chromosome-level genome assembly using single-molecule long-read sequencing, and the Hi-C technique. The assembled genome was 454.28 Mb, containing 14 chromosomes, with contig N50 of 2.14 Mb. A total of 46,495 protein-coding genes, 169.74 Mb repeat regions, and 31.76 Kb non-coding RNA were predicted. P. anserina diverged from Potentilla micrantha ∼28.52 million years ago (Mya). Furthermore, P. anserina underwent a recent tetraploidization ∼6.4 Mya. The species-specific genes were enriched in Starch and sucrose metabolism and Galactose metabolism pathways. We identified the sub-genome structures of P. anserina, with A sub-genome was larger than B sub-genome and closer to P. micrantha phylogenetically. Despite lacking significant genome-wide expression dominance, the A sub-genome had higher homoeologous gene expression in shoot apical meristem, flower and tuberous root. The resistance genes was contracted in P. anserina genome. Key genes involved in starch biosynthesis were expanded and highly expressed in tuberous roots, which probably drives the tuber formation. The genomics and transcriptomics data generated in this study advance our understanding of the genomic landscape of P. anserina, and will accelerate genetic studies and breeding programs.

Identification of genome compositions in allopolyploid species of the genus Elymus (Poaceae: Triticeae) in the Asian part of Russia by CAPS analysis

The genus Elymus L., together with wheat, rye, and barley, belongs to the tribe Triticeae. Apart from its economic value, this tribe is characterized by abundance of polyploid taxa formed in the course of remote hybridization. Single-copy nuclear genes are convenient markers for identification of source genomes incorporated into polyploids. In the present work, a CAPS-marker is developed to distinguish basic St, H, and Y genomes comprising polyploid genomes of Asiatic species of the genus Elymus. The test is based on electrophoretic analysis of restriction patterns of a PCR-amplified fragment of the gene coding for beta-amylase. There are about 50 Elymus species in Russia, and most of them are supposed to possess one of three haplome combinations, StH, StY and StHY. Boreal StH-genomic species endemic for Russia are the least studied. On the basis of nucleotide sequences from public databases, TaqI restrictase was selected, as it produced patterns of restriction fragments specific for St, H, and Y haplomes easily recognizable in agarose gel. A sample of 68 accessions belonging to 32 species was analyzed. In 15 species, the earlier known genomic constitutions were confirmed, but in E. kamoji this assay failed to reveal the presence of H genome. This unusual H genome was suggested to originate from a different Hordeum species. In 16 species, genomic constitutions were identified for the first time. Fifteen accessions from Asian Russia possessed the genomic constitution StStHH, and E. amurensis, phylogenetically close to the StY-genomic species E. ciliaris, had the genomic constitution StStYY. It is inferred that the center of species diversity of the StH-genomic group is shifted to the north as compared to the center of origin of StY-genomic species, confined to China. Key words: Elymus; taxonomy; allopolyploids; genome constitution; CAPS markers.

Inherited allelic variants and novel karyotype changes influence fertility and genome stability in Brassica allohexaploids

Synthetic allohexaploid Brassica hybrids (2n = AABBCC) do not exist naturally, but can be synthesized by crosses between diploid and/or allotetraploid Brassica species. Using these hybrids, we aimed to identify how novel allohexaploids restore fertility and normal meiosis after formation. Chromosome inheritance, genome structure, fertility and meiotic behaviour were assessed in three segregating allohexaploid populations derived from the cross (B. napus × B. carinata) × B. juncea using a combination of molecular marker genotyping, phenotyping and cytogenetics. Plants with unbalanced A-C translocations in one direction (where a C-genome chromosome fragment replaces an A-genome fragment) but not the other (where an A-genome fragment replaces a C-genome fragment) showed significantly reduced fertility across all populations. Genomic regions associated with fertility contained several meiosis genes with putatively causal mutations inherited from the parents (copies of SCC2 in the A genome, PAIR1/PRD3, PRD1 and ATK1/KATA/KIN14a in the B genome, and MSH2 and SMC1/TITAN8 in the C genome). Reduced seed fertility associated with the loss of chromosome fragments from only one subgenome following homoeologous exchanges could comprise a mechanism for biased genome fractionation in allopolyploids. Pre-existing meiosis gene variants present in allotetraploid parents may help to stabilize meiosis in novel allohexaploids.

Effect of Hybridization on Somatic Mutations and Genomic Rearrangements in Plants

Hybridization has been routinely practiced in agriculture to enhance the crop yield. Principally, it can cause hybrid vigor where hybrid plants display increased size, biomass, fertility, and resistance to diseases, when compared to their parents. During hybridization, hybrid offspring receive a genomic shock due to mixing of distant parental genomes, which triggers a myriad of genomic rearrangements, e.g., transpositions, genome size changes, chromosomal rearrangements, and other effects on the chromatin. Recently, it has been reported that, besides genomic rearrangements, hybridization can also alter the somatic mutation rates in plants. In this review, we provide in-depth insights about hybridization triggered genomic rearrangements and somatic mutations in plants.

Adaptation of the Pivotal-Differential Genome Pattern for the Induction of Intergenomic Chromosome Recombination in Hybrids of Synthetic Amphidiploids within Triticeae Tribe

A pivotal-differential evolution pattern is when two allopolyploids share a common genome, which is called pivotal, and differ with respect to the other genome or genomes, called differential. This feature induces the intergenomic recombination between chromosomes of differential genomes, which can lead to speciation. Our study is a cytomolecular insight into this mechanism which was adapted for the induction of intergenomic chromosome recombination in hybrids of synthetic amphidiploids Aegilops biuncialis × S. cereale (UUMMRR) and triticale (AABBRR) where R-genome was pivotal. We observed chromosome recombination events which were induced by both: (1) random chromosome fragmentation and non-homologous chromosome end joining at mitosis of root meristem cells and (2) intergenomic chromosome associations at meiosis of pollen mother cells (PMCs) of F₁ hybrids. Reciprocal chromosome translocations were identified in six F₁ plants and 15 plants of F₂ generation using fluorescence in situ hybridization (FISH) with DNA clones (pTa-86, pTa-k374, pTa-465, pTa-535, pTa-k566, and pTa-713). We observed signals of pTa-86, pTa-535, and pTa-k566 probes in several chromosome breakpoints. The comparison of the DNA clone sequences distinguished a number of common motifs, which can be considered as characteristics of chromosome breakpoint loci. Immunodetection of synaptonemal complex proteins and genomic in situ hybridization analysis at meiosis of PMCs of F₁ hybrids showed, that the homologous pairing of pivotal R-genome chromosomes is crucial for the fertility of F₁ hybrids, however, these chromosomes can be also involved in the intergeneric recombination.

The evolutionary history of sea barley (Hordeum marinum) revealed by comparative physical mapping of repetitive DNA

BACKGROUND AND AIMS

Hordeum marinum is a species complex that includes the diploid subspecies marinum and both diploid and tetraploid forms of gussoneanum. Their relationships, the rank of the taxa and the origin of the polyploid forms remain points of debate. The present work reports a comparative karyotype analysis of six H. marinum accessions representing all taxa and cytotypes.

METHODS

Karyotypes were determined by analysing the chromosomal distribution of several tandemly repeated sequences, including the Triticeae cloned probes pTa71, pTa794, pAs1 and pSc119·2 and the simple sequence repeats (SSRs) (AG)10, (AAC)5, (AAG)5, (ACT)5 and (ATC)5.

KEY RESULTS

The identification of each chromosome pair in all subspecies and cytotypes is reported for the first time. Homologous relationships are also established. Wide karyotypic differences were detected within marinum accessions. Specific chromosomal markers characterized and differentiated the genomes of marinum and diploid gussoneanum. Two subgenomes were detected in the tetraploids. One of these had the same chromosome complement as diploid gussoneanum; the second subgenome, although similar to the chromosome complement of diploid H. marinum sensu lato, appeared to have no counterpart in the marinum accessions analysed here.

CONCLUSIONS

The tetraploid forms of gussoneanum appear to have come about through a cross between a diploid gussoneanum progenitor and a second, related-but unidentified-diploid ancestor. The results reveal the genome structure of the different H. marinum taxa and demonstrate the allopolyploid origin of the tetraploid forms of gussoneanum.