Chromosome number reduction in the sister clade of Carica papaya with concomitant genome size doubling

Repeat turnover meets stable chromosomes: repetitive DNA sequences mark speciation and gene pool boundaries in sugar beet and wild beets

Sugar beet and its wild relatives share a base chromosome number of nine and similar chromosome morphologies. Yet, interspecific breeding is impeded by chromosome and sequence divergence that is still not fully understood. Since repetitive DNAs are among the fastest evolving parts of the genome, we investigated, if repeatome innovations and losses are linked to chromosomal differentiation and speciation. We traced genome and chromosome-wide evolution across 13 beet species comprising all sections of the genera Beta and Patellifolia. For this, we combined short and long read sequencing, flow cytometry, and cytogenetics to build a comprehensive framework that spans the complete scale from DNA to chromosome to genome. Genome sizes and repeat profiles reflect the separation into three gene pools with contrasting evolutionary patterns. Among all repeats, satellite DNAs harbor most genomic variability, leading to fundamentally different centromere architectures, ranging from chromosomal uniformity in Beta and Patellifolia to the formation of patchwork chromosomes in Corollinae/Nanae. We show that repetitive DNAs are causal for the genome expansions and contractions across the beet genera, providing insights into the genomic underpinnings of beet speciation. Satellite DNAs in particular vary considerably between beet genomes, leading to the evolution of distinct chromosomal setups in the three gene pools, likely contributing to the barriers in beet breeding. Thus, with their isokaryotypic chromosome sets, beet genomes present an ideal system for studying the link between repeats, genomic variability, and chromosomal differentiation and provide a theoretical fundament for understanding barriers in any crop breeding effort.

Completing the TRB family: newly characterized members show ancient evolutionary origins and distinct localization, yet similar interactions

Telomere repeat binding proteins (TRBs) belong to a family of proteins possessing a Myb-like domain which binds to telomeric repeats. Three members of this family (TRB1, TRB2, TRB3) from Arabidopsis thaliana have already been described as associated with terminal telomeric repeats (telomeres) or short interstitial telomeric repeats in gene promoters (telo-boxes). They are also known to interact with several protein complexes: telomerase, Polycomb repressive complex 2 (PRC2) E(z) subunits and the PEAT complex (PWOs-EPCRs-ARIDs-TRBs). Here we characterize two novel members of the TRB family (TRB4 and TRB5). Our wide phylogenetic analyses have shown that TRB proteins evolved in the plant kingdom after the transition to a terrestrial habitat in Streptophyta, and consequently TRBs diversified in seed plants. TRB4-5 share common TRB motifs while differing in several others and seem to have an earlier phylogenetic origin than TRB1-3. Their common Myb-like domains bind long arrays of telomeric repeats in vitro, and we have determined the minimal recognition motif of all TRBs as one telo-box. Our data indicate that despite the distinct localization patterns of TRB1-3 and TRB4-5 in situ, all members of TRB family mutually interact and also bind to telomerase/PRC2/PEAT complexes. Additionally, we have detected novel interactions between TRB4-5 and EMF2 and VRN2, which are Su(z)12 subunits of PRC2.

Comparative genomics analysis of endangered wild Egyptian Moringa peregrina (Forssk.) Fiori plastome, with implications for the evolution of Brassicales order

Moringa is a mono-genus belonging to the Moringaceae family, which includes 13 species. Among them, Moringa peregrina is plant species native to the Arabian Peninsula, Southern Sinai in Egypt, and the Horn of Africa, and comprehensive studies on its nutritional, industrial, and medicinal values have been performed. Herein, we sequenced and analyzed the initial complete chloroplast genome of Moringa peregrina. Concurrently, we analyzed the new chloroplast genome along with 25 chloroplast genomes related to species representing eight families in the Brassicales order. The results indicate that the plastome sequence of M. peregrina consists of 131 genes, with an average GC content of 39.23%. There is a disparity in the IR regions of the 26 species ranging from 25,804 to 31,477 bp. Plastome structural variations generated 20 hotspot regions that could be considered prospective DNA barcode locations in the Brassicales order. Tandem repeats and SSR structures are reported as significant evidence of structural variations among the 26 tested specimens. Furthermore, selective pressure analysis was performed to estimate the substitution rate within the Moringaceae family, which revealing that the ndhA and accD genes are under positive selective pressure. The phylogenetic analysis of the Brassicales order produced an accurate monophyletic annotation cluster of the Moringaceae and Capparaceae species, offering unambiguous identification without overlapping groups between M. oleifera and M. peregrina, which are genetically strongly associated. Divergence time estimation suggests that the two Moringa species recently diversified, 0.467 Ma. Our findings highlight the first complete plastome of the Egyptian wild-type of M. peregrina, which can be used for determining plastome phylogenetic relationships and systematic evolution history within studies on the Moringaceae family.

Down, then up: non-parallel genome size changes and a descending chromosome series in a recent radiation of the Australian allotetraploid plant species, Nicotiana section Suaveolentes (Solanaceae)

BACKGROUND AND AIMS

The extent to which genome size and chromosome numbers evolve in concert is little understood, particularly after polyploidy (whole-genome duplication), when a genome returns to a diploid-like condition (diploidization). We study this phenomenon in 46 species of allotetraploid Nicotiana section Suaveolentes (Solanaceae), which formed <6 million years ago and radiated in the arid centre of Australia.

METHODS

We analysed newly assessed genome sizes and chromosome numbers within the context of a restriction site-associated nuclear DNA (RADseq) phylogenetic framework.

KEY RESULTS

RADseq generated a well-supported phylogenetic tree, in which multiple accessions from each species formed unique genetic clusters. Chromosome numbers and genome sizes vary from n = 2x = 15 to 24 and 2.7 to 5.8 pg/1C nucleus, respectively. Decreases in both genome size and chromosome number occur, although neither consistently nor in parallel. Species with the lowest chromosome numbers (n = 15-18) do not possess the smallest genome sizes and, although N. heterantha has retained the ancestral chromosome complement, n = 2x = 24, it nonetheless has the smallest genome size, even smaller than that of the modern representatives of ancestral diploids.

CONCLUSIONS

The results indicate that decreases in genome size and chromosome number occur in parallel down to a chromosome number threshold, n = 20, below which genome size increases, a phenomenon potentially explained by decreasing rates of recombination over fewer chromosomes. We hypothesize that, more generally in plants, major decreases in genome size post-polyploidization take place while chromosome numbers are still high because in these stages elimination of retrotransposons and other repetitive elements is more efficient. Once such major genome size change has been accomplished, then dysploid chromosome reductions take place to reorganize these smaller genomes, producing species with small genomes and low chromosome numbers such as those observed in many annual angiosperms, including Arabidopsis.

Genomes, repeatomes and interphase chromosome organization in the meadowfoam family (Limnanthaceae, Brassicales)

The meadowfoam family (Limnanthaceae) is one of the smallest and genomically underexplored families of the Brassicales. The Limnanthaceae harbor about seven species in the genus Limnanthes (meadowfoam) and Floerkea proserpinacoides (false mermaidweed), all native to North America. Because all Limnanthes and Floerkea species have only five chromosome pairs, i.e., a chromosome number rare in Brassicales and shared with Arabidopsis thaliana (Arabidopsis), we examined the Limnanthaceae genomes as a potential model system. Using low-coverage whole-genome sequencing data, we reexamined phylogenetic relationships and characterized the repeatomes of Limnanthaceae genomes. Phylogenies based on complete chloroplast and 35S rDNA sequences corroborated the sister relationship between Floerkea and Limnanthes and two major clades in the latter genus. The genome size of Limnanthaceae species ranges from 1.5 to 2.1 Gb, apparently due to the large increase in DNA repeats, which constitute 60-70% of their genomes. Repeatomes are dominated by long terminal repeat retrotransposons, while tandem repeats represent only less than 0.5% of the genomes. The average chromosome size in Limnanthaceae species (340-420 Mb) is more than 10 times larger than in Arabidopsis (32 Mb). A three-dimensional fluorescence in situ hybridization analysis demonstrated that the five chromosome pairs in interphase nuclei of Limnanthes species adopt the Rabl-like configuration.

Interstitial Arabidopsis-Type Telomeric Repeats in Asteraceae

Tandem repeats of telomeric-like motifs at intra-chromosomal regions, known as interstitial telomeric repeats (ITR), have drawn attention as potential markers of structural changes, which might convey information about evolutionary relationships if preserved through time. Building on our previous work that reported outstanding ITR polymorphisms in the genus Anacyclus, we undertook a survey across 132 Asteraceae species, focusing on the six most speciose subfamilies and considering all the ITR data published to date. The goal was to assess whether the presence, site number, and chromosomal location of ITRs convey any phylogenetic signal. We conducted fluorescent in situ hybridization (FISH) using an Arabidopsis-type telomeric sequence as a probe on karyotypes obtained from mitotic chromosomes. FISH signals of ITR sites were detected in species of subfamilies Asteroideae, Carduoideae, Cichorioideae, Gymnarhenoideae, and Mutisioideae, but not in Barnadesioideae. Although six small subfamilies have not yet been sampled, altogether, our results suggest that the dynamics of ITR formation in Asteraceae cannot accurately trace the complex karyological evolution that occurred since the early diversification of this family. Thus, ITRs do not convey a reliable signal at deep or shallow phylogenetic levels and cannot help to delimitate taxonomic categories, a conclusion that might also hold for other important families such as Fabaceae.

Interstitial Telomeric-like Repeats (ITR) in Seed Plants as Assessed by Molecular Cytogenetic Techniques: A Review

The discovery of telomeric repeats in interstitial regions of plant chromosomes (ITRs) through molecular cytogenetic techniques was achieved several decades ago. However, the information is scattered and has not been critically evaluated from an evolutionary perspective. Based on the analysis of currently available data, it is shown that ITRs are widespread in major evolutionary lineages sampled. However, their presence has been detected in only 45.6% of the analysed families, 26.7% of the sampled genera, and in 23.8% of the studied species. The number of ITR sites greatly varies among congeneric species and higher taxonomic units, and range from one to 72 signals. ITR signals mostly occurs as homozygous loci in most species, however, odd numbers of ITR sites reflecting a hemizygous state have been reported in both gymnosperm and angiosperm groups. Overall, the presence of ITRs appears to be poor predictors of phylogenetic and taxonomic relatedness at most hierarchical levels. The presence of ITRs and the number of sites are not significantly associated to the number of chromosomes. The longitudinal distribution of ITR sites along the chromosome arms indicates that more than half of the ITR presences are between proximal and terminal locations (49.5%), followed by proximal (29.0%) and centromeric (21.5%) arm regions. Intraspecific variation concerning ITR site number, chromosomal locations, and the differential presence on homologous chromosome pairs has been reported in unrelated groups, even at the population level. This hypervariability and dynamism may have likely been overlooked in many lineages due to the very low sample sizes often used in cytogenetic studies.

Evolutionary diversification of cytokinin-specific glucosyltransferases in angiosperms and enigma of missing cis-zeatin O-glucosyltransferase gene in Brassicaceae

In the complex process of homeostasis of phytohormones cytokinins (CKs), O-glucosylation catalyzed by specific O-glucosyltransferases represents one of important mechanisms of their reversible inactivation. The CK O-glucosyltransferases belong to a highly divergent and polyphyletic multigene superfamily of glycosyltransferases, of which subfamily 1 containing UDP-glycosyltransferases (UGTs) is the largest in the plant kingdom. It contains recently discovered O and P subfamilies present in higher plant species but not in Arabidopsis thaliana. The cis-zeatin O-glucosyltransferase (cisZOG) genes belong to the O subfamily encoding a stereo-specific O-glucosylation of cis-zeatin-type CKs. We studied different homologous genes, their domains and motifs, and performed a phylogenetic reconstruction to elucidate the plant evolution of the cisZOG gene. We found that the cisZOG homologs do not form a clear separate clade, indicating that diversification of the cisZOG gene took place after the diversification of the main angiosperm families, probably within genera or closely related groups. We confirmed that the gene(s) from group O is(are) not present in A. thaliana and is(are) also missing in the family Brassicaceae. However, cisZOG or its metabolites are found among Brassicaceae clade, indicating that remaining genes from other groups (UGT73-group D and UGT85-group G) are able, at least in part, to substitute the function of group O lost during evolution. This study is the first detailed evolutionary evaluation of relationships among different plant ZOGs within angiosperms.

Karyotype evolution in Phalaris (Poaceae): The role of reductional dysploidy, polyploidy and chromosome alteration in a wide-spread and diverse genus

Karyotype characteristics can provide valuable information on genome evolution and speciation, in particular in taxa with varying basic chromosome numbers and ploidy levels. Due to its worldwide distribution, remarkable variability in morphological traits and the fact that ploidy change plays a key role in its evolution, the canary grass genus Phalaris (Poaceae) is an excellent study system to investigate the role of chromosomal changes in species diversification and expansion. Phalaris comprises diploid species with two basic chromosome numbers of x = 6 and 7 as well as polyploids based on x = 7. To identify distinct karyotype structures and to trace chromosome evolution within the genus, we apply fluorescence in situ hybridisation (FISH) of 5S and 45S rDNA probes in four diploid and four tetraploid Phalaris species of both basic numbers. The data agree with a dysploid reduction from x = 7 to x = 6 as the result of reciprocal translocations between three chromosomes of an ancestor with a diploid chromosome complement of 2n = 14. We recognize three different genomes in the genus: (1) the exclusively Mediterranean genome A based on x = 6, (2) the cosmopolitan genome B based on x = 7 and (3) a genome C based on x = 7 and with a distribution in the Mediterranean and the Middle East. Both auto- and allopolyploidy of genomes B and C are suggested for the formation of tetraploids. The chromosomal divergence observed in Phalaris can be explained by the occurrence of dysploidy, the emergence of three different genomes, and the chromosome rearrangements accompanied by karyotype change and polyploidization. Mapping the recognized karyotypes on the existing phylogenetic tree suggests that genomes A and C are restricted to sections Phalaris and Bulbophalaris, respectively, while genome B occurs across all taxa with x = 7.

Cytogenetic comparison of heteromorphic and homomorphic sex chromosomes in Coccinia (Cucurbitaceae) points to sex chromosome turnover

Our understanding of the evolution of plant sex chromosomes is increasing rapidly due to high-throughput sequencing data and phylogenetic and molecular-cytogenetic approaches that make it possible to infer the evolutionary direction and steps leading from homomorphic to heteromorphic sex chromosomes. Here, we focus on four species of Coccinia, a genus of 25 dioecious species, including Coccinia grandis, the species with the largest known plant Y chromosome. Based on a phylogeny for the genus, we selected three species close to C. grandis to test the distribution of eight repetitive elements including two satellites, and several plastid and mitochondrial probes, that we had previously found to have distinct accumulation patterns in the C. grandis genome. Additionally, we determined C-values and performed immunostaining experiments with (peri-)centromere-specific antibodies on two species (for comparison with C. grandis). In spite of no microscopic chromosomal heteromorphism, single pairs of chromosomes in male cells of all three species accumulate some of the very same repeats that are enriched on the C. grandis Y chromosome, pointing to either old (previous) sex chromosomes or incipient (newly arising) ones, that is, to sex chromosome turnover. A 144-bp centromeric satellite repeat (CgCent) that characterizes all C. grandis chromosomes except the Y is highly abundant in all centromeric regions of the other species, indicating that the centromeric sequence of the Y chromosome diverged very recently.

Development of plants resistant to Papaya leaf distortion mosaic virus by intergeneric hybridization between Carica papaya and Vasconcellea cundinamarcensis

In this study, we confirmed that Vasconcellea cundinamarcensis resists Papaya leaf distortion mosaic virus (PLDMV), and used it to produce intergeneric hybrids with Carica papaya. From the cross between C. papaya and V. cundinamarcensis, we obtained 147 seeds with embryos. Though C. papaya is a monoembryonic plant, multiple embryos were observed in all 147 seeds. We produced 218 plants from 28 seeds by means of embryo-rescue culture. All plants had pubescence on their petioles and stems characteristic of V. cundinamarcensis. Flow cytometry and PCR of 28 plants confirmed they were intergeneric hybrids. To evaluate virus resistance, mechanical inoculation of PLDMV was carried out. The test showed that 41 of 134 intergeneric hybrid plants showed no symptoms and were resistant. The remaining 93 hybrids showed necrotic lesions on the younger leaves than the inoculated leaves. In most of the 93 hybrids, the necrotic lesions enclosed the virus and prevented further spread. These results suggest that the intergeneric hybrids will be valuable material for PLDMV-resistant papaya breeding.