Polyploidy and genome evolution in plants

Phylogenomics of the Andean Tetraploid Clade of the American Amaryllidaceae (Subfamily Amaryllidoideae): Unlocking a Polyploid Generic Radiation Abetted by Continental Geodynamics

One of the two major clades of the endemic American Amaryllidaceae subfam. Amaryllidoideae constitutes the tetraploid-derived (n = 23) Andean-centered tribes, most of which have 46 chromosomes. Despite progress in resolving phylogenetic relationships of the group with plastid and nrDNA, certain subclades were poorly resolved or weakly supported in those previous studies. Sequence capture using anchored hybrid enrichment was employed across 95 species of the clade along with five outgroups and generated sequences of 524 nuclear genes and a partial plastome. Maximum likelihood phylogenetic analyses were conducted on concatenated supermatrices, and coalescent-based species tree analyses were run on the gene trees, followed by hybridization network, age diversification and biogeographic analyses. The four tribes Clinantheae, Eucharideae, Eustephieae, and Hymenocallideae (sister to Clinantheae) are resolved in all analyses with > 90 and mostly 100% support, as are almost all genera within them. Nuclear gene supermatrix and species tree results were largely in concordance; however, some instances of cytonuclear discordance were evident. Hybridization network analysis identified significant reticulation in Clinanthus, Hymenocallis, Stenomesson and the subclade of Eucharideae comprising Eucharis, Caliphruria, and Urceolina. Our data support a previous treatment of the latter as a single genus, Urceolina, with the addition of Eucrosia dodsonii. Biogeographic analysis and penalized likelihood age estimation suggests an origin in the Cauca, Desert and Puna Neotropical bioprovinces for the complex in the mid-Oligocene, with more dispersals than vicariances in its history, but no extinctions. Hymenocallis represents the only instance of long-distance vicariance from the tropical Andean origin of its tribe Hymenocallideae. The absence of extinctions correlates with the lack of diversification rate shifts within the clade. The Eucharideae experienced a sudden lineage radiation ca. 10 Mya. We tie much of the divergences in the Andean-centered lineages to the rise of the Andes, and suggest that the Amotape-Huancabamba Zone functioned as both a corridor (dispersal) and a barrier to migration (vicariance). Several taxonomic changes are made. This is the largest DNA sequence data set to be applied within Amaryllidaceae to date.

Synthetic auxin herbicides: finding the lock and key to weed resistance

Synthetic auxin herbicides are designed to mimic indole-3-acetic acid (IAA), an integral plant hormone affecting cell growth, development, and tropism. In this review, we explore target site genes in the auxin signaling pathway including SCF^TIR1/AFB, Aux/IAA, and ARFs that are confirmed or proposed mechanisms for weed resistance to synthetic auxin herbicides. Resistance to auxin herbicides by metabolism, either by enhanced cytochrome P450 detoxification or by loss of pro-herbicide activation, is a major non-target-site resistance pathway. We speculate about potential fitness costs of resistance due to effects of resistance-conferring mutations, provide insight into the role of polyploidy in synthetic auxin resistance evolution, and address the genetic resources available for weeds. This knowledge will be the key to unlock the long-standing questions as to which components of the auxin signaling pathway are most likely to have a role in resistance evolution. We propose that an ambitious research effort into synthetic auxin herbicide/target site interactions is needed to 1) explain why some synthetic auxin chemical families have activity on certain dicot plant families but not others and 2) fully elucidate target-site cross-resistance patterns among synthetic auxin chemical families to guide best practices for resistance management.

Evolutionary dynamics of genome size in a radiation of woody plants

PREMISE

Plant genome size ranges widely, providing many opportunities to examine how genome size variation affects plant form and function. We analyzed trends in chromosome number, genome size, and leaf traits for the woody angiosperm clade Viburnum to examine the evolutionary associations, functional implications, and possible drivers of genome size.

METHODS

Chromosome counts and genome size estimates were mapped onto a Viburnum phylogeny to infer the location and frequency of polyploidization events and trends in genome size evolution. Genome size was analyzed with leaf anatomical and physiological data to evaluate the influence of genome size on plant function.

RESULTS

We discovered nine independent polyploidization events, two reductions in base chromosome number, and substantial variation in genome size with a slight trend toward genome size reduction in polyploids. We did not find strong relationships between genome size and the functional and morphological traits that have been highlighted at broader phylogenetic scales.

CONCLUSIONS

Polyploidization events were sometimes associated with rapid radiations, demonstrating that polyploid lineages can be highly successful. Relationships between genome size and plant physiological function observed at broad phylogenetic scales may be largely irrelevant to the evolutionary dynamics of genome size at smaller scales. The view that plants readily tolerate changes in ploidy and genome size, and often do so, appears to apply to Viburnum.

Evolutionary analysis of the Moringa oleifera genome reveals a recent burst of plastid to nucleus gene duplications

It is necessary to identify suitable alternative crops to ensure the nutritional demands of a growing global population. The genome of Moringa oleifera, a fast-growing drought-tolerant orphan crop with highly valuable agronomical, nutritional and pharmaceutical properties, has recently been reported. We model here gene family evolution in Moringa as compared with ten other flowering plant species. Despite the reduced number of genes in the compact Moringa genome, 101 gene families, grouping 957 genes, were found as significantly expanded. Expanded families were highly enriched for chloroplastidic and photosynthetic functions. Indeed, almost half of the genes belonging to Moringa expanded families grouped with their Arabidopsis thaliana plastid encoded orthologs. Microsynteny analysis together with modeling the distribution of synonymous substitutions rates, supported most plastid duplicated genes originated recently through a burst of simultaneous insertions of large regions of plastid DNA into the nuclear genome. These, together with abundant short insertions of plastid DNA, contributed to the occurrence of massive amounts of plastid DNA in the Moringa nuclear genome, representing 4.71%, the largest reported so far. Our study provides key genetic resources for future breeding programs and highlights the potential of plastid DNA to impact the structure and function of nuclear genes and genomes.

Fungal evolution: cellular, genomic and metabolic complexity

The question of how phenotypic and genomic complexity are inter-related and how they are shaped through evolution is a central question in biology that historically has been approached from the perspective of animals and plants. In recent years, however, fungi have emerged as a promising alternative system to address such questions. Key to their ecological success, fungi present a broad and diverse range of phenotypic traits. Fungal cells can adopt many different shapes, often within a single species, providing them with great adaptive potential. Fungal cellular organizations span from unicellular forms to complex, macroscopic multicellularity, with multiple transitions to higher or lower levels of cellular complexity occurring throughout the evolutionary history of fungi. Similarly, fungal genomes are very diverse in their architecture. Deep changes in genome organization can occur very quickly, and these phenomena are known to mediate rapid adaptations to environmental changes. Finally, the biochemical complexity of fungi is huge, particularly with regard to their secondary metabolites, chemical products that mediate many aspects of fungal biology, including ecological interactions. Herein, we explore how the interplay of these cellular, genomic and metabolic traits mediates the emergence of complex phenotypes, and how this complexity is shaped throughout the evolutionary history of Fungi.

Characterization and Comparative Analysis of RWP-RK Proteins from Arachis duranensis, Arachis ipaensis, and Arachis hypogaea

RWP-RK proteins are important factors involved in nitrate response and gametophyte development in plants, and the functions of RWP-RK proteins have been analyzed in many species. However, the characterization of peanut RWP-RK proteins is limited. In this study, we identified 16, 19, and 32 RWP-RK members from Arachis duranensis, Arachis ipaensis, and Arachis hypogaea, respectively, and investigated their evolution relationships. The RWP-RK proteins were classified into two groups, RWP-RK domain proteins and NODULE-INCEPTION-like proteins. Chromosomal distributions, gene structures, and conserved motifs of RWP-RK genes were compared among wild and cultivated peanuts. In addition, we identified 12 orthologous gene pairs from the two wild peanut species, 13 from A. duranensis and A. hypogaea, and 13 from A. ipaensis and A. hypogaea. One, one, and seventeen duplicated gene pairs were identified within the A. duranensis, A. ipaensis, and A. hypogaea genomes, respectively. Moreover, different numbers of cis-acting elements in the RWP-RK promoters were found in wild and cultivated species (87 in A. duranensis, 89 in A. ipaensis, and 92 in A. hypogaea), and as a result, many RWP-RK genes showed distinct expression patterns in different tissues. Our study will provide useful information for further functional and evolutionary analysis of the RWP-RK genes.

The Genomes of the Allohexaploid Echinochloa crus-galli and Its Progenitors Provide Insights into Polyploidization-Driven Adaptation

The hexaploid species Echinochloa crus-galli is one of the most detrimental weeds in crop fields, especially in rice paddies. Its evolutionary history is similar to that of bread wheat, arising through polyploidization after hybridization between a tetraploid and a diploid species. In this study, we generated and analyzed high-quality genome sequences of diploid (E. haploclada), tetraploid (E. oryzicola), and hexaploid (E. crus-galli) Echinochloa species. Gene family analysis showed a significant loss of disease-resistance genes such as those encoding NB-ARC domain-containing proteins during Echinochloa polyploidization, contrary to their significant expansionduring wheat polyploidization, suggesting that natural selection might favor reduced investment in resistance in this weed to maximize its growth and reproduction. In contrast to the asymmetric patterns of genome evolution observed in wheat and other crops, no significant differences in selection pressure were detected between the subgenomes in E. oryzicola and E. crus-galli. In addition, distinctive differences in subgenome transcriptome dynamics during hexaploidization were observed between E. crus-galli and bread wheat. Collectively, our study documents genomic mechanisms underlying the adaptation of a major agricultural weed during polyploidization. The genomic and transcriptomic resources of three Echinochloa species and new insights into the polyploidization-driven adaptive evolution would be useful for future breeding cereal crops.

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.

Patterns and impacts of nonvertical evolution in eukaryotes: a paradigm shift

Evolution of eukaryotic species and their genomes has been traditionally understood as a vertical process in which genetic material is transmitted from parents to offspring along a lineage, and in which genetic exchange is restricted within species boundaries. However, mounting evidence from comparative genomics indicates that this paradigm is often violated. Horizontal gene transfer and mating between diverged lineages blur species boundaries and challenge the reconstruction of evolutionary histories of species and their genomes. Nonvertical evolution might be more restricted in eukaryotes than in prokaryotes, yet it is not negligible and can be common in certain groups. Recognition of such processes brings about the need to incorporate this complexity into our models, as well as to conceptually reframe eukaryotic diversity and evolution. Here, I review the recent work from genomics studies that supports the effects of nonvertical modes of evolution including introgression, hybridization, and horizontal gene transfer in different eukaryotic groups. I then discuss emerging patterns and effects, illustrated by specific examples, that support the conclusion that nonvertical processes are often at the root of important evolutionary transitions and adaptations. I will argue that a paradigm shift is needed to naturally accommodate nonvertical processes in eukaryotic evolution.

Cancer regeneration: Polyploid cells are the key drivers of tumor progression

In spite of significant advancements of therapies for initial eradication of cancers, tumor relapse remains a major challenge. It is for a long time known that polyploid malignant cells are a main source of resistance against chemotherapy and irradiation. However, therapeutic approaches targeting these cells have not been appropriately pursued which could partly be due to the shortage of knowledge on the molecular biology of cell polyploidy. On the other hand, there is a rising trend to appreciate polyploid/ multinucleated cells as key players in tissue regeneration. In this review, we suggest an analogy between the functions of polyploid cells in normal and malignant tissues and discuss the idea that cell polyploidy is an evolutionary conserved source of tissue regeneration also exploited by cancers as a survival factor. In addition, polyploid cells are highlighted as a promising therapeutic target to overcome drug resistance and relapse.

Evolutionary directions of single nucleotide substitutions and structural mutations in the chloroplast genomes of the family Calycanthaceae

BACKGROUND

Chloroplast genome sequence data is very useful in studying/addressing the phylogeny of plants at various taxonomic ranks. However, there are no empirical observations on the patterns, directions, and mutation rates, which are the key topics in chloroplast genome evolution. In this study, we used Calycanthaceae as a model to investigate the evolutionary patterns, directions and rates of both nucleotide substitutions and structural mutations at different taxonomic ranks.

RESULTS

There were 2861 polymorphic nucleotide sites on the five chloroplast genomes, and 98% of polymorphic sites were biallelic. There was a single-nucleotide substitution bias in chloroplast genomes. A → T or T → A (2.84%) and G → C or C → G (3.65%) were found to occur significantly less frequently than the other four transversion mutation types. Synonymous mutations kept balanced pace with nonsynonymous mutations, whereas biased directions appeared between transition and transversion mutations and among transversion mutations. Of the structural mutations, indels and repeats had obvious directions, but microsatellites and inversions were non-directional. Structural mutations increased the single nucleotide mutations rates. The mutation rates per site per year were estimated to be 0.14-0.34 × 10- 9 for nucleotide substitution at different taxonomic ranks, 0.64 × 10- 11 for indels and 1.0 × 10- 11 for repeats.

CONCLUSIONS

Our direct counts of chloroplast genome evolution events provide raw data for correctly modeling the evolution of sequence data for phylogenetic inferences.

Successful mating and hybridisation in two closely related flatworm species despite significant differences in reproductive morphology and behaviour

Reproductive traits are some of the fastest diverging characters and can serve as reproductive barriers. The free-living flatworm Macrostomum lignano, and its congener M. janickei are closely related, but differ substantially in their male intromittent organ (stylet) morphology. Here, we examine whether these morphological differences are accompanied by differences in behavioural traits, and whether these could represent barriers to successful mating and hybridization between the two species. Our data shows that the two species differ in many aspects of their mating behaviour. Despite these differences, the species mate readily with each other in heterospecific pairings. Although both species have similar fecundity in conspecific pairings, the heterospecific pairings revealed clear postmating barriers, as few heterospecific pairings produced F1 hybrids. These hybrids had a stylet morphology that was intermediate between that of the parental species, and they were fertile. Finally, using a mate choice experiment, we show that the nearly two-fold higher mating rate of M. lignano caused it to mate more with conspecifics, leading to assortative mating, while M. janickei ended up mating more with heterospecifics. Thus, while the two species can hybridize, the mating rate differences could possibly lead to higher fitness costs for M. janickei compared to M. lignano.

Gradual polyploid genome evolution revealed by pan-genomic analysis of Brachypodium hybridum and its diploid progenitors

Our understanding of polyploid genome evolution is constrained because we cannot know the exact founders of a particular polyploid. To differentiate between founder effects and post polyploidization evolution, we use a pan-genomic approach to study the allotetraploid Brachypodium hybridum and its diploid progenitors. Comparative analysis suggests that most B. hybridum whole gene presence/absence variation is part of the standing variation in its diploid progenitors. Analysis of nuclear single nucleotide variants, plastomes and k-mers associated with retrotransposons reveals two independent origins for B. hybridum, ~1.4 and ~0.14 million years ago. Examination of gene expression in the younger B. hybridum lineage reveals no bias in overall subgenome expression. Our results are consistent with a gradual accumulation of genomic changes after polyploidization and a lack of subgenome expression dominance. Significantly, if we did not use a pan-genomic approach, we would grossly overestimate the number of genomic changes attributable to post polyploidization evolution.

Genome-wide identification and characterization of nonspecific lipid transfer protein (nsLTP) genes in Arachis duranensis

Nonspecific lipid transfer proteins (nsLTPs) play vital roles in lipid metabolism, cell apoptosis and biotic and abiotic stresses in plants. However, the distribution of nsLTPs in Arachis duranensis has not been fully characterized. In this study, we identified 64 nsLTP genes in A. duranensis (designated AdLTPs), which were classified into six subfamilies and randomly distributed along nine chromosomes. Tandem and segmental duplication events were detected in the evolution of AdLTPs. The K_s and ω values differed significantly between Types 1 and D subfamilies, and eight AdLTPs were under positive selection. The expression levels of AdLTPs were changed after salinity, PEG, low-temperature and ABA treatments. Three AdLTPs were associated with resistance to nematode infection, and DOF and WRI1 transcription factors may regulate the AdLTP response to nematode infection. Our results may provide valuable genomic information for the breeding of peanut cultivars that are resistant to biotic and abiotic stresses.

Genomics of Evolutionary Novelty in Hybrids and Polyploids

It has long been recognized that hybridization and polyploidy are prominent processes in plant evolution. Although classically recognized as significant in speciation and adaptation, recognition of the importance of interspecific gene flow has dramatically increased during the genomics era, concomitant with an unending flood of empirical examples, with or without genome doubling. Interspecific gene flow is thus increasingly thought to lead to evolutionary innovation and diversification, via adaptive introgression, homoploid hybrid speciation and allopolyploid speciation. Less well understood, however, are the suite of genetic and genomic mechanisms set in motion by the merger of differentiated genomes, and the temporal scale over which recombinational complexity mediated by gene flow might be expressed and exposed to natural selection. We focus on these issues here, considering the types of molecular genetic and genomic processes that might be set in motion by the saltational event of genome merger between two diverged species, either with or without genome doubling, and how these various processes can contribute to novel phenotypes. Genetic mechanisms include the infusion of new alleles and the genesis of novel structural variation including translocations and inversions, homoeologous exchanges, transposable element mobilization and novel insertional effects, presence-absence variation and copy number variation. Polyploidy generates massive transcriptomic and regulatory alteration, presumably set in motion by disrupted stoichiometries of regulatory factors, small RNAs and other genome interactions that cascade from single-gene expression change up through entire networks of transformed regulatory modules. We highlight both these novel combinatorial possibilities and the range of temporal scales over which such complexity might be generated, and thus exposed to natural selection and drift.

A Reappraisal of the Evolutionary and Developmental Pathway of Apomixis and Its Genetic Control in Angiosperms

Apomixis sensu stricto (agamospermy) is asexual reproduction by seed. In angiosperms it represents an easy byway of life cycle renewal through gamete-like cells that give rise to maternal embryos without ploidy reduction (meiosis) and ploidy restitution (syngamy). The origin of apomixis still represents an unsolved problem, as it may be either evolved from sex or the other way around. This review deals with a reappraisal of the origin of apomixis in order to deepen knowledge on such asexual mode of reproduction which seems mainly lacking in the most basal angiosperm orders (i.e., Amborellales, Nymphaeales and Austrobaileyales, also known as ANA-grade), while it clearly occurs in different forms and variants in many unrelated families of monocots and eudicots. Overall findings strengthen the hypothesis that apomixis as a whole may have evolved multiple times in angiosperm evolution following different developmental pathways deviating to different extents from sexuality. Recent developments on the genetic control of apomixis in model species are also presented and adequately discussed in order to shed additional light on the antagonist theories of gain- and loss-of-function over sexuality.

Hybridization and the origin of new yeast lineages

Hybrids originate from the mating of two diverged organisms, resulting in novel lineages that have chimeric genomes. Hybrids may exhibit unique phenotypic traits that are not necessarily intermediate between those present in the progenitors. These unique traits may enable them to thrive in new environments. Many hybrid lineages have been discovered among yeasts in the Saccharomycotina, of which many have industrial or clinical relevance, but this might reflect a bias toward investigating species with relevance to humans. Hybridization has also been proposed to be at the root of the whole-genome duplication in the lineage leading to Saccharomyces cerevisiae. Thus, hybridization seems to have played a prominent role in the evolution of Saccharomycotina yeasts, although it is still unclear how common this evolutionary process has been during the evolution of this and other fungal clades. Similarly, the evolutionary aftermath of hybridization, including implications at the genomic, transcriptional, physiological or ecological levels, remains poorly understood. In this review, I survey recent findings from genomic analysis of yeast hybrids of industrial or clinical relevance, and discuss the evolutionary implications of genomic hybridization for the origin of new lineages, including when such hybridization results in a whole-genome duplication.

Polyploidization events shaped the transcription factor repertoires in legumes (Fabaceae)

Transcription factors (TFs) are essential for plant growth and development. Several legumes (e.g. soybean) are rich sources of protein and oil and have great economic importance. Here we report a phylogenomic analysis of TF families in legumes and their potential association with important traits (e.g. nitrogen fixation). We used TF DNA-binding domains to systematically screen the genomes of 15 leguminous and five non-leguminous species. Transcription factor orthologous groups (OGs) were used to estimate OG sizes in ancestral nodes using a gene birth-death model, which allowed the identification of lineage-specific expansions. The OG analysis and rate of synonymous substitutions show that major TF expansions are strongly associated with whole-genome duplication (WGD) events in the legume (approximately 58 million years ago) and Glycine (approximately 13 million years ago) lineages, which account for a large fraction of the Phaseolus vulgaris and Glycine max TF repertoires. Of the 3407 G. max TFs, 1808 and 676 have homeologs within single syntenic regions in Phaseolus vulgaris and Vitis vinifera, respectively. We found a trend for TFs expanded in legumes to be preferentially transcribed in roots and nodules, supporting their recruitment early in the evolution of nodulation in the legume clade. Some families also showed count differences between G. max and the wild soybean Glycine soja, including genes located within important quantitative trait loci. Our findings strongly support the roles of two WGDs in shaping the TF repertoires in the legume and Glycine lineages, and these are probably related to important aspects of legume and soybean biology.

Preferential Disomic Segregation and C. micrantha/C. medica Interspecific Recombination in Tetraploid 'Giant Key' Lime; Outlook for Triploid Lime Breeding

The triploid 'Tahiti' lime (C. x latifolia (Yu. Tanaka) Tanaka) naturally originated from a merger between a haploid ovule of lemon (C. x limon (L.) Burm) and a diploid pollen from a 'Mexican' lime (C. x aurantiifolia (Christm.) Swing). The very limited natural inter-varietal diversity and gametic sterility of C. latifolia requires a phylogenomic based reconstruction breeding strategy to insure its diversification. We developed a strategy based on interploid hybridization between diploid lemon and the doubled diploid 'Giant Key' lime. This lime is a doubled diploid of 'Mexican' lime, itself a natural interspecific F1 hybrid between C. medica L. and C. micrantha Wester. For an optimized breeding program, we analyzed the meiotic behavior of the allotetraploid lime, the genetic structure of its diploid gametes, the interspecific recombination between C. medica and C. micrantha, and constructed its genetic map. A population of 272 triploid hybrids was generated using 'Giant Key' lime as pollinator. One hundred fifty-eight SNPs diagnostic of C. micrantha, regularly distributed throughout the citrus genome were successfully developed and applied. The genetic structure of the diploid gametes was examined based on C. micrantha doses along the genome. The diploid gametes transmitted in average 91.17% of the parental interspecific C. medica/C. micrantha heterozygosity. Three chromosomes (2, 8, and 9) showed disomic segregation with high preferential pairing values, while the remaining chromosomes showed an intermediate inheritance with a preferential disomic trend. A total of 131 SNPs were assigned to nine linkage groups to construct the genetic map. It spanned 272.8 cM with a low average recombination rate (0.99 cM Mb-1) and high synteny and colinearity with the reference clementine genome. Our results confirmed that an efficient reconstruction breeding strategy for 'Tahiti' lime is possible, based on interploid hybridization using a doubled diploid of C. aurantiifolia. The tetraploid parent should be selected for favorable agronomic traits and its genetic value should be efficiently inherited by the progeny thanks to transmission of the high level of parental heterozygosity. However, it would require developing numerous progeny to overcome the linkage drag caused by the limited interspecific recombination associated with the predominant disomic inheritance.

3-D Nucleus Architecture in Oat × Maize Addition Lines

The nucleus architecture of hybrid crop plants is not a well-researched topic, yet it can have important implications for their genetic stability and usefulness in the successful expression of agronomically desired traits. In this work we studied the spatial distribution of introgressed maize chromatin in oat × maize addition lines with the number of added maize chromosomes varying from one to four. The number of chromosome additions was confirmed by genomic in situ hybridization (GISH). Maize chromosome-specific simple sequence repeat (SSR) markers were used to identify the added chromosomes. GISH on 3-D root and leaf nuclei was performed to assess the number, volume, and position of the maize-chromatin occupied regions. We revealed that the maize chromosome territory (CT) associations of varying degree prevailed in the double disomic lines, while CT separation was the most common distribution pattern in the double monosomic line. In all analyzed lines, the regions occupied by maize CTs were located preferentially at the nuclear periphery. A comparison between the tissues showed that the maize CTs in the leaf nuclei are positioned closer to the center of the nucleus than in the root nuclei. These findings shed more light on the processes that shape the nucleus architecture in hybrids.

Regain of sex determination system and sexual reproduction ability in a synthetic octoploid male fish

Polyploids in vertebrates are generally associated with unisexual reproduction, but the direct consequences of polyploidy on sex determination system and reproduction mode remain unknown. Here, we synthesized a group of artificial octoploids between unisexual gynogenetic hexaploid Carassius gibelio and sexual tetraploid Carassius auratus. The synthetic octoploids were revealed to have more than 200 chromosomes, in which 50 chromosomes including the X/Y sex determination system were identified to transfer from sexual tetraploid C. auratus into the unisexual gynogenetic hexaploid C. gibelio. Significantly, a few synthetic octoploid males were found to be fertile, and one octoploid male was confirmed to regain sexual reproduction ability, which exhibits characteristics that are the same to sexual reproduction tetraploid males, such as 1:1 sex ratio occurrence, meiosis completion and euploid sperm formation in spermatogenesis, as well as normal embryo development and gene expression pattern during embryogenesis. Therefore, the current finding provides a unique case to explore the effect of sex determination system incorporation on reproduction mode transition from unisexual gynogenesis to sexual reproduction along with genome synthesis of recurrent polyploidy in vertebrates.

Polyploidization-driven differentiation of freezing tolerance in Solidago canadensis

Solidago canadensis, originating from the temperate region of North America, has expanded southward to subtropical regions through polyploidization. Here we investigated whether freezing tolerance of S. canadensis was weakened during expansion. Measurement of the temperature causing 50% ruptured cells (LT₅₀ ) in 35 S. canadensis populations revealed ploidy-related differentiation in freezing tolerance. Freezing tolerance was found to decrease with increasing ploidy. The polyploid populations of S. canadensis had lower ScICE1 gene expression levels but more ScICE1 gene copies than the diploids. Furthermore, more DNA methylation sites in the ScICE1 gene promoter were detected in the polyploids than in the diploids. The results suggest that promoter methylation represses the expression of multi-copy ScICE1 genes, leading to weaker freezing tolerance in polyploid S. canadensis compared to the diploids. The study provides empirical evidence that DNA methylation regulates expression of the gene copies and supports polyploidization-driven adaptation to new environments.

Comprehensive analysis of Translationally Controlled Tumor Protein (TCTP) provides insights for lineage-specific evolution and functional divergence

BACKGROUND

Translationally controlled tumor protein (TCTP) is a conserved, multifunctional protein involved in numerous cellular processes in eukaryotes. Although the functions of TCTP have been investigated sporadically in animals, invertebrates, and plants, few lineage-specific activities of this molecule, have been reported. An exception is in Arabidopsis thaliana, in which TCTP (AtTCTP1) functions in stomatal closuer by regulating microtubule stability. Further, although the development of next-generation sequencing technologies has facilitated the analysis of many eukaryotic genomes in public databases, inter-kingdom comparative analyses using available genome information are comparatively scarce.

METHODOLOGY

To carry out inter-kingdom comparative analysis of TCTP, TCTP genes were identified from 377 species. Then phylogenetic analysis, prediction of protein structure, molecular docking simulation and molecular dynamics analysis were performed to investigate the evolution of TCTP genes and their binding proteins.

RESULTS

A total of 533 TCTP genes were identified from 377 eukaryotic species, including protozoa, fungi, invertebrates, vertebrates, and plants. Phylogenetic and secondary structure analyses reveal lineage-specific evolution of TCTP, and inter-kingdom comparisons highlight the lineage-specific emergence of, or changes in, secondary structure elements in TCTP proteins from different kingdoms. Furthermore, secondary structure comparisons between TCTP proteins within each kingdom, combined with measurements of the degree of sequence conservation, suggest that TCTP genes have evolved to conserve protein secondary structures in a lineage-specific manner. Additional tertiary structure analysis of TCTP-binding proteins and their interacting partners and docking simulations between these proteins further imply that TCTP gene variation may influence the tertiary structures of TCTP-binding proteins in a lineage-specific manner.

CONCLUSIONS

Our analysis suggests that TCTP has undergone lineage-specific evolution and that structural changes in TCTP proteins may correlate with the tertiary structure of TCTP-binding proteins and their binding partners in a lineage-specific manner.

Gossypium Genomics: Trends, Scope, and Utilization for Cotton Improvement

Cotton (Gossypium spp.) is the most important natural fiber crop worldwide. The diversity of Gossypium species also provides an ideal model for investigating evolution and domestication of polyploids. However, the huge and complex cotton genome hinders genomic research. Technical advances in high-throughput sequencing and bioinformatics analysis have now largely overcome these obstacles, bringing about a new era of cotton genomics. Here, we review recent progress in Gossypium genomics based on whole genome sequencing, resequencing, and comparative genomics, which have provided insights about the genomic basis of fiber biogenesis and the landscape of cotton functional genomics. We address current challenges and present multidisciplinary genomics-enabled breeding strategies covering the breadth of high fiber yield, quality, and environmental resilience for future cotton breeding programs.

Homoeolog expression bias and expression level dominance (ELD) in four tissues of natural allotetraploid Brassica napus

Background

Allopolyploidy is widespread in angiosperms, and they can coordinate two or more different genomes through genetic and epigenetic modifications to exhibit stronger vigor and adaptability. To explore the changes in homologous gene expression patterns in the natural allotetraploid Brassica napus (A_nA_nC_nC_n) relative to its two diploid progenitors, B. rapa (A_rA_r) and B. oleracea (C_oC_o), after approximately 7500 years of domestication, the global gene pair expression patterns in four major tissues (stems, leaves, flowers and siliques) of these three species were analyzed using an RNA sequencing approach.

Results

The results showed that the ‘transcriptomic shock’ phenomenon was alleviated in natural B. napus after approximately 7500 years of natural domestication, and most differentially expressed genes (DEGs) in B. napus were downregulated relative to those in its two diploid progenitors. The KEGG analysis indicated that three pathways related to photosynthesis were enriched in both comparison groups (A_nA_nC_nC_n vs A_rA_r and A_nA_nC_nC_n vs C_oC_o), and these pathways were all downregulated in four tissues of B. napus. In addition, homoeolog expression bias and expression level dominance (ELD) in B. napus were thoroughly studied through analysis of expression levels of 27,609 B. rapa-B. oleracea orthologous gene pairs. The overwhelming majority of gene pairs (an average of 86.7%) in B. napus maintained their expression pattern in two diploid progenitors, and approximately 78.1% of the gene pairs showed expression bias with a preference toward the A subgenome. Overall, an average of 48, 29.7 and 22.3% homologous gene pairs exhibited additive expression, ELD and transgressive expression in B. napus, respectively. The ELD bias varies from tissue to tissue; specifically, more gene pairs in stems and siliques showed ELD-A, whereas the opposite was observed in leaves and flowers. More transgressive upregulation, rather than downregulation, was observed in gene pairs of B. napus.

Conclusions

In general, these results may provide a comprehensive understanding of the changes in homologous gene expression patterns in natural B. napus after approximately 7500 years of evolution and domestication and may enhance our understanding of allopolyploidy.

Phylotranscriptomic Analyses Reveal Asymmetrical Gene Duplication Dynamics and Signatures of Ancient Polyploidy in Mints

Ancient duplication events and retained gene duplicates have contributed to the evolution of many novel plant traits and, consequently, to the diversity and complexity within and across plant lineages. Although mounting evidence highlights the importance of whole-genome duplication (WGD; polyploidy) and its key role as an evolutionary driver, gene duplication dynamics and mechanisms, both of which are fundamental to our understanding of evolutionary process and patterns of plant diversity, remain poorly characterized in many clades. We use newly available transcriptomic data and a robust phylogeny to investigate the prevalence, occurrence, and timing of gene duplications in Lamiaceae (mints), a species-rich and chemically diverse clade with many ecologically, economically, and culturally important species. We also infer putative WGDs—an extreme mechanism of gene duplication—using large-scale data sets from synonymous divergence (K_S), phylotranscriptomic, and divergence time analyses. We find evidence for widespread but asymmetrical levels of gene duplication and ancient polyploidy in Lamiaceae that correlate with species richness, including pronounced levels of gene duplication and putative ancient WGDs (7–18 events) within the large subclade Nepetoideae and up to 10 additional WGD events in other subclades. Our results help disentangle WGD-derived gene duplicates from those produced by other mechanisms and illustrate the nonuniformity of duplication dynamics in mints, setting the stage for future investigations that explore their impacts on trait diversity and species diversification. Our results also provide a practical context for evaluating the benefits and limitations of transcriptome-based approaches to inferring WGD, and we offer recommendations for researchers interested in investigating ancient WGDs in other plant groups.

Inferring the impacts of evolutionary history and ecological constraints on spore size and shape in the ferns

PREMISE

In the ferns, cell size has been explored with spores, which are largely uniform within species, produced in abundance, and durable. However, spore size and shape have been variously defined, and the relationship of these traits to genome size has not been well established. Here, we explore the variation in fern spore size and shape by ploidy level and genome size.

METHODS

Measurements of spore dimensions for two study sets of ferns, Polystichum and Adiantum, both including diploid and tetraploid taxa, provided the basis for computing estimates of shape and size as defined here. These traits were compared between and within ploidy levels and regressed with genome size estimates from flow cytometry analysis.

RESULTS

All size traits were strongly correlated with genome size; the shape trait was weakly correlated. Tetraploids were larger than diploids as expected; the spores of some closely related diploid species were also significantly different in size.

DISCUSSION

Researchers with access to a student-grade compound microscope can develop a valid estimate of relative genome size for ferns. These estimates provide enough resolution to infer ploidy level and explore the relationship between genome size, habitat, and physiological constraints for species within ploidy level.

Genome size variation at constant chromosome number is not correlated with repetitive DNA dynamism in Anacyclus (Asteraceae)

BACKGROUND AND AIMS

Changes in the amount of repetitive DNA (dispersed and tandem repeats) are considered the main contributors to genome size variation across plant species in the absence of polyploidy. However, the study of repeatome dynamism in groups showing contrasting genomic features and complex evolutionary histories is needed to determine whether other processes underlying genome size variation may have been overlooked. The main aim here was to elucidate which mechanism best explains genome size evolution in Anacyclus (Asteraceae).

METHODS

Using data from Illumina sequencing, we analysed the repetitive DNA in all species of Anacyclus, a genus with a reticulate evolutionary history, which displays significant genome size and karyotype diversity albeit presenting a stable chromosome number.

KEY RESULTS

By reconstructing ancestral genome size values, we inferred independent episodes of genome size expansions and contractions during the evolution of the genus. However, analysis of the repeatome revealed a similar DNA repeat composition across species, both qualitative and quantitative. Using comparative methods to study repeatome dynamics in the genus, we found no evidence for repeat activity causing genome size variation among species.

CONCLUSIONS

Our results, combined with previous cytogenetic data, suggest that genome size differences in Anacyclus are probably related to chromosome rearrangements involving losses or gains of chromosome fragments, possibly associated with homoploid hybridization. These could represent balanced rearrangements that do not disrupt gene dosage in merged genomes, for example via chromosome segment exchanges.

Has the Polyploid Wave Ebbed?

There was a wave of whole genome duplications (WGD) during and subsequent to the K-Pg interface, which was followed by an increase in the proportion of species that were polyploid. I consider why this wave of polyploid speciation has continued to rise through the divergent evolution of polyploid lineages, and through rounds of homoploid and heteroploid chromosomal change. I also consider why the polyploid speciation wave is likely to rise in the next millennium. I propose that the speed of polyploid genesis through ploidal increase and through diversification among polyploids likely will be greater than the speed of diploid speciation. The increase in polyploid diversity is expected to lag well behind episodes of WGD, owing to the very long period required for species diversification either by lineage splitting or additional rounds of polyploidy, in addition to the long period of genomic adjustment to higher ploidal levels in neopolyploids.

Comparative transcriptomic analyses of chlorogenic acid and luteolosides biosynthesis pathways at different flowering stages of diploid and tetraploid Lonicera japonica

The Flos Lonicerae Japonicae (FLJ), Lonicera japonica Thunb, belonging to the Caprifoliaceae family, is an economically important plant that is highly utilized in traditional Chinese medicine as well as in Japanese medicine. The flowers of these plants are rich in chlorogenic acid (CGA) and luteoloside. Our previous study revealed that tetraploid L. japonica has higher fresh/dry weight, phenolic acids and flavonoids contents than those of diploid plants. However, why tetraploid L. japonica can yield higher CGA and luteolosides than that in diploid and what is the difference in the molecular regulatory mechanism of these pathways between diploid and tetraploids remained unclear. Therefore, in the present study, we performed comprehensive transcriptome analyses of different flowering stages of diploid and tetraploid L. japonica. The CGA content of tetraploid was found higher than that of diploid at all the growth stages. While the luteoloside content of diploid was higher than that of tetraploid at S4 and S6 growth stages. We obtained a high-quality transcriptome assembly (N50 = 2,055 bp; Average length = 1,331 bp) compared to earlier studies. Differential expression analysis revealed that several important genes involving in plant hormone signal transduction, carbon metabolism, starch/sucrose metabolism and plant-pathogen interaction were upregulated in tetraploid compared with the diploid L. japonica, reflecting the higher adaptability and resistance of tetraploid species. Furthermore, by associating the phenotypic data and gene expression profiles, we were able to characterize the potential molecular regulatory mechanism of important biosynthetic pathways at different flowering stages. Overall, our work provides a foundation for further research on these important secondary metabolite pathways and their implications in traditional Chinese/Japanese medicine.

Hexaploid sweetpotato (Ipomoea batatas (L.) Lam.) may not be a true type to either auto- or allopolyploid

To better define the sweetpotato polyploidy, we sought to reconstruct phylogenies of its subgenomes based on hybridization networks that could trace reticulate lineages of differentiated homoeolog triplets of multiple single-copy genes. In search of such homoeolog triplets, we distinguished cDNA variants of 811 single-copy Conserved Ortholog Set II (COSII) genes from two sweetpotato clones into variation partitions specified by corresponding homologs from two I. trifida lines, I. tenuissima and I. littoralis using a phylogenetic partition method, and amplicon variants of the COSII-marker regions from 729 of these genes from two sweetpotato clones into putative homoeoallele groups using haplotype tree and the partition methods referenced by corresponding homologs from I. tenuissima. These analyses revealed partly or completely differentiated expressed-homoeologs and homoeologs from a majority of these genes with three important features. 1. Two variation types: the predominant interspecific variations (homoeoalleles), which are non-randomly clustered, differentially interspecifically conserved or sweetpotato-specific, and the minor intraspecific ones (alleles), which are randomly distributed mostly at non-interspecifically variable sites, and usually sweetpotato-specific. 2. A clear differentiation of cDNA variants of many COSII genes into the variation partition specified by I. tenuissima or I. littoralis from that by I. trifida. 3. Three species-homolog-specified and one sweetpotato-specific variation partitions among 293 different COSII cDNAs, and two or three out of the four partitions among cDNA variants of 306 COSII genes. We then constructed hybridization networks from two concatenations of 16 and 4 alignments of 8 homologous COSII cDNA regions each, which included three taxa of expressed homoeologs in a triple-partition combination from the 16 or 4 sweetpotato COSII genes and 5 taxa each of respective cDNA homologs from the three sweetpotato relatives and I. nil, and inferred a species tree embodying both networks. The species tree predicted close-relative origins of three partly differentiated sweetpotato subgenomes.

Polyploidy in gymnosperms - Insights into the genomic and evolutionary consequences of polyploidy in Ephedra

While polyploidization is recognized as a major evolutionary driver for ferns and angiosperms, little is known about its impact in gymnosperms, where polyploidy is much less frequent. We explore Ephedra to evaluate (i) the extent of genome size diversity in the genus and the influence polyploidy has had on the evolution of nuclear DNA contents, and (ii) identify where shifts in genome size and polyploidy have occurred both temporally and spatially. A phylogenetic framework of all Ephedra species together with genome sizes and karyotypes for 87% and 67% of them respectively, were used to explore ploidy evolution and its global distribution patterns. Polyploidy was shown to be extremely common, with 41 species (83%) being polyploid (up to 8×) or having polyploid cytotypes - the highest frequency and level reported for any gymnosperm. Genome size was also diverse, with values ranging ~5-fold (8.09-38.34 pg/1C) - the largest range for any gymnosperm family - and increasing in proportion to ploidy level (i.e. no genome downsizing). Our findings provide novel data which support the view that gymnosperms have a more conserved mode of genomic evolution compared with angiosperms.

Expression patterns and functional divergence of homologous genes accompanied by polyploidization in cotton (Gossypium hirsutum L.)

Naturally allotetraploid cotton has been widely used as an ideal model to investigate gene expression remodeling as a consequence of polyploidization. However, the global gene pattern variation during early fiber development was unknown. In this study, through RNA-seq technology, we comprehensively investigated the expression patterns of homologous genes between allotetraploid cotton (G. hirsutum) and its diploid progenitors (G. arboreum and G. raimondii) at the fiber early development stage. In tetraploid cotton, genes showed expression level dominance (ELD) bias toward the A genome. This phenomenon was explained by the up-/downregulation of the homologs from the nondominant progenitor (D genome). Gene ontology (GO) enrichment results indicated that the ELD-A genes might be a prominent cause responsible for fiber property change through regulating the fatty acid biosynthesis/metabolism and microtubule procession, and the ELD-D genes might be involved in transcription regulation and stress inducement. In addition, the number and proportion of completely A- and D-subfunctionalized gene were similar at different fiber development stages. However, for neofunctionalization, the number and proportion of reactivated D-derived genes were greater than those of A at 3 and 5 DPA. Eventually, we found that some homologous genes belonging to several specific pathways might create novel asymmetric transcripts between two subgenomes during polyploidization and domestication process, further making the fiber property meet the human demands. Our study identified determinate pathways and their involved genes between allotetraploid cotton and their progenitors at early fiber development stages, providing new insights into the mechanism of cotton fiber evolution.

Homeologs of Brassica SOC1, a central regulator of flowering time, are differentially regulated due to partitioning of evolutionarily conserved transcription factor binding sites in promoters

Evolution of Brassica genome post-polyploidization reveals asymmetrical genome fractionation and copy number variation. Herein, we describe the impact of promoter divergence among SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) homeologs on expression and function in Brassica spp. SOC1, a regulated floral pathway integrator, is conserved as 3 redundant homeologs in diploid Brassicas. Even with high sequence identity within coding regions (92.8-100%), the spatio-temporal expression patterns of 9 SOC1 homologs in B. juncea and B. nigra indicates regulatory divergence. While LF and MF2 SOC1 homeologs are upregulated during floral transition, MF1 is barely expressed. Also, MF2 homeolog levels do not decline post-flowering, unlike LF. To investigate the underlying source of divergence, we analyzed the sequence and phylogeny of all reported (22) and isolated (21) upstream regions of Brassica SOC1. Full length upstream regions (4712-19189 bp) reveal 5 ubiquitously conserved ancestral Blocks, harboring binding sites of 18 TFs (TFBSs) characterized in Arabidopsis thaliana. The orthologs of these TFBSs are differentially conserved among Brassica SOC1 homeologs, imparting expression divergence. No crucial TFBSs are exclusively lost from LF_SOC1 promoter, while MF1_SOC1 has lost NF-Y binding site crucial for SOC1 activation by CONSTANS. MF2_SOC1 homeologs have lost important TFBSs (SEP3, AP1 and SMZ), responsible for SOC1 repression post-flowering. BjuAALF_SOC1 promoter (proximal 2 kb) shows ubiquitous reporter expression in B. juncea cv. Varuna transgenics, while BjuAAMF1_SOC1 promoter shows absence of reporter expression, validating the impact of TFBS divergence. Conservation of the original primary protein sequence is discovered in B. rapa homeologs (46) of 18 TFs. Co-regulation pattern of these TFs appeared similar for B. rapa LF and MF2 SOC1 homeologs; MF1 shows significant variation. Strong regulatory association is recorded for AP1, AP2, SEP3, FLC and CONSTANS/NF-Y, highlighting their importance in homeolog-specific SOC1 regulation. Correlation of B. juncea AP1, AP2 and FLC expression with SOC1 homeologs also complies with the TFBS differences. We thus conclude that redundant SOC1 loci contribute differentially to cumulative expression of SOC1 due to divergent selection of ancestral TFBSs.

Ecological setup, ploidy diversity, and reproductive biology of Paspalum modestum, a promising wetland forage grass from South America

With ever-rising demand for food, forage breeding for intensification of cattle production is also taking a leap. In South America, cattle production systems are displaced to marginal areas poorly exploited with cultivated pastures yet with high potential for growing stocking rates. This places the need for using native genetic resources to breed locally adapted plant genotypes that benefits from better forage quality, yield, and lesser threat to the local biodiversity. Paspalum modestum Mez is a grass species that produces quality forage and grows in marginal areas like estuaries and floodplains, suitable for introduction in breeding programs. In this study we characterize the species' reproductive biology and ecological preferences needed beforehand any improvement. P. modestum plants found in nature are commonly diploids, rarely triploids, and tetraploids. Chromosome associations during meiosis in polyploids indicate they are autopolyploids. While diploids are sexual self-sterile, analyses of embryology, gamete fertility and experimental crossings show tetraploids are self-compatible facultative apomicts, highly fertile and have a high proportion of sexuality compared to other apomictic species. Ecological niche analysis and species distribution modelling show mean annual temperature and precipitation as main ecological drivers and a wide geographical area of climatic suitability where P. modestum can grow and be exploited as a forage grass. Our study points to P. modestum as a native plant resource appropriate for breeding waterlogging tolerant ecotypes and genotypes of high biomass production adapted to low flow areas in the Subtropics of Brazil, Paraguay, Uruguay and Argentina.

Preferential gene retention increases the robustness of cold regulation in Brassicaceae and other plants after polyploidization

Cold stress profoundly affects plant growth and development and is a key factor affecting the geographic distribution and evolution of plants. Plants have evolved adaptive mechanisms to cope with cold stress. Here, through the genomic analysis of Arabidopsis, three Brassica species and 17 other representative species, we found that both cold-related genes (CRGs) and their collinearity were preferentially retained after polyploidization followed by genome instability, while genome-wide gene sets exhibited a variety of other expansion mechanisms. The cold-related regulatory network was increased in Brassicaceae genomes, which were recursively affected by polyploidization. By combining our findings regarding the selective retention of CRGs from this ecological genomics study with the available knowledge of cold-induced chromosome doubling, we hypothesize that cold stress may have contributed to the success of polyploid plants through both increasing polyploidization and selectively maintaining CRGs during evolution. This hypothesis requires further biological and ecological exploration to obtain solid supporting evidence, which will potentially contribute to understanding the generation of polyploids and to the field of ecological genomics.

The modular nature of protein evolution: domain rearrangement rates across eukaryotic life

BACKGROUND

Modularity is important for evolutionary innovation. The recombination of existing units to form larger complexes with new functionalities spares the need to create novel elements from scratch. In proteins, this principle can be observed at the level of protein domains, functional subunits which are regularly rearranged to acquire new functions.

RESULTS

In this study we analyse the mechanisms leading to new domain arrangements in five major eukaryotic clades (vertebrates, insects, fungi, monocots and eudicots) at unprecedented depth and breadth. This allows, for the first time, to directly compare rates of rearrangements between different clades and identify both lineage specific and general patterns of evolution in the context of domain rearrangements. We analyse arrangement changes along phylogenetic trees by reconstructing ancestral domain content in combination with feasible single step events, such as fusion or fission. Using this approach we explain up to 70% of all rearrangements by tracing them back to their precursors. We find that rates in general and the ratio between these rates for a given clade in particular, are highly consistent across all clades. In agreement with previous studies, fusions are the most frequent event leading to new domain arrangements. A lineage specific pattern in fungi reveals exceptionally high loss rates compared to other clades, supporting recent studies highlighting the importance of loss for evolutionary innovation. Furthermore, our methodology allows us to link domain emergences at specific nodes in the phylogenetic tree to important functional developments, such as the origin of hair in mammals.

CONCLUSIONS

Our results demonstrate that domain rearrangements are based on a canonical set of mutational events with rates which lie within a relatively narrow and consistent range. In addition, gained knowledge about these rates provides a basis for advanced domain-based methodologies for phylogenetics and homology analysis which complement current sequence-based methods.

A Phylogenomic Perspective on Evolution and Discordance in the Alpine-Arctic Plant Clade Micranthes (Saxifragaceae)

The increased availability of large phylogenomic datasets is often accompanied by difficulties in disentangling and harnessing the data. These difficulties may be enhanced for species resulting from reticulate evolution and/or rapid radiations producing large-scale discordance. As a result, there is a need for methods to investigate discordance, and in turn, use this conflict to inform and aid in downstream analyses. Therefore, we drew upon multiple analytical tools to investigate the evolution of Micranthes (Saxifragaceae), a clade of primarily arctic-alpine herbs impacted by reticulate and rapid radiations. To elucidate the evolution of Micranthes we sought near-complete taxon sampling with multiple accessions per species and assembled extensive nuclear (518 putatively single copy loci) and plastid (95 loci) datasets. In addition to a robust phylogeny for Micranthes, this research shows that genetic discordance presents a valuable opportunity to develop hypotheses about its underlying causes, such as hybridization, polyploidization, and range shifts. Specifically, we present a multi-step approach that incorporates multiple checks points for paralogy, including reciprocally blasting targeted genes against transcriptomes, running paralogy checks during the assembly step, and grouping genes into gene families to look for duplications. We demonstrate that a thorough assessment of discordance can be a source of evidence for evolutionary processes that were not adequately captured by a bifurcating tree model, and helped to clarify processes that have structured the evolution of Micranthes.

Sequential Paleotetraploidization shaped the carrot genome

BACKGROUND

Carrot (Daucus carota subsp. carota L.) is an important root crop with an available high-quality genome. The carrot genome is thought to have undergone recursive paleo-polyploidization, but the extent, occurrences, and nature of these events are not clearly defined.

RESULTS

Using a previously published comparative genomics pipeline, we reanalysed the carrot genome and characterized genomic fractionation, as well as gene loss and retention, after each of the two tetraploidization events and inferred a dominant and sensitive subgenome for each event. In particular, we found strong evidence of two sequential tetraploidization events, with one (Dc-α) approximately 46-52 million years ago (Mya) and the other (Dc-β) approximately 77-87 Mya, both likely allotetraploidization in nature. The Dc-β event was likely common to all Apiales plants, occurring around the divergence of Apiales-Bruniales and after the divergence of Apiales-Asterales, likely playing an important role in the derivation and divergence of Apiales species. Furthermore, we found that rounds of polyploidy events contributed to the expansion of gene families responsible for plastidial methylerythritol phosphate (MEP), the precursor of carotenoid accumulation, and shaped underlying regulatory pathways. The alignment of orthologous and paralogous genes related to different events of polyploidization and speciation constitutes a comparative genomics platform for studying Apiales, Asterales, and many other related species.

CONCLUSIONS

Hierarchical inference of homology revealed two tetraploidization events that shaped the carrot genome, which likely contributed to the successful establishment of Apiales plants and the expansion of MEP, upstream of the carotenoid accumulation pathway.

Genetic Contribution of Paleopolyploidy to Adaptive Evolution in Angiosperms

Ancient whole-genome duplications (WGDs or polyploidy) are prevalent in plants, and some WGDs occurred during the timing of severe global environmental changes. It has been suggested that WGDs may have contributed to plant adaptation. However, this still lacks empirical evidence at the genetic level to support the hypothesis. Here, we investigated the survivors of gene duplicates from multiple ancient WGD events on the major branches of angiosperm phylogeny, and aimed to explore genetic evidence supporting the significance of polyploidy. Duplicated genes co-retained from three waves of independent WGDs (∼120 million years ago [Ma], ∼66, and <20 Ma) were investigated in 25 selected species. Gene families functioning in low temperature and darkness were commonly retained gene duplicates after the eight independently occurring WGDs in many lineages around the Cretaceous-Paleocene boundary, when the global cooling and darkness were the two main stresses. Moreover, the commonly retained duplicates could be key factors which may have contributed to the robustness of the critical stress-related pathways. In addition, genome-wide transcription factors (TFs) functioning in stresses tend to retain duplicates after waves of WGDs, and the coselected gene duplicates in many lineages may play critical roles during severe environmental stresses. Collectively, these results shed new light on the significant contribution of paleopolyploidy to plant adaptation during global environmental changes in the evolutionary history of angiosperms.

Current Challenges in Studying Alternative Splicing in Plants: The Case of Physcomitrella patens SR Proteins

To colonize different terrestrial habitats, early land plants had to overcome the challenge of coping with harsh new environments. Alternative splicing - an RNA processing mechanism through which splice sites are differentially recognized, originating multiple transcripts and potentially different proteins from the same gene - can be key for plant stress tolerance. Serine/arginine-rich (SR) proteins constitute an evolutionarily conserved family of major alternative splicing regulators that in plants subdivides into six subfamilies. Despite being well studied in animals and a few plant species, such as the model angiosperm Arabidopsis thaliana and the crop Oryza sativa, little is known of these splicing factors in early land plants. Establishing the whole complement of SR proteins in different species is essential to understand the functional and evolutionary significance of alternative splicing. An in silico search for SR proteins in the extant moss Physcomitrella patens revealed inconsistencies both in the published data and available databases, likely arising from automatic annotation lacking adequate manual curation. These misannotations interfere with the description not only of the number and subfamily classification of Physcomitrella SR proteins but also of their domain architecture, potentially hindering the elucidation of their molecular functions. We therefore advise caution when looking into P. patens genomic resources. Our systematic survey nonetheless confidently identified 16 P. patens SR proteins that fall into the six described subfamilies and represent counterparts of well-established members in Arabidopsis and rice. Intensified research efforts should disclose whether SR proteins were already determining alternative splicing modulation and stress tolerance in early land plants.

Ploidy-Dependent Effects of Light Stress on the Mode of Reproduction in the Ranunculus auricomus Complex (Ranunculaceae)

Polyploidy in angiosperms is an influential factor to trigger apomixis, the reproduction of asexual seeds. Apomixis is usually facultative, which means that both sexual and apomictic seeds can be formed by the same plant. Environmental abiotic stress, e.g. light stress, can change the frequency of apomixis. Previous work suggested effects of stress treatments on meiosis and megasporogenesis. We hypothesized that polyploidy would alter the stress response and hence reproductive phenotypes of different cytotypes. The main aims of this research were to explore with prolonged photoperiods, whether polyploidy alters proportions of sexual ovule and sexual seed formation under light stress conditions. We used three facultative apomictic, pseudogamous cytotypes of the Ranunculus auricomus complex (diploid, tetraploid, and hexaploid). Stress treatments were applied by extended light periods (16.5 h) and control (10 h) in climate growth chambers. Proportions of apomeiotic vs. meiotic development in the ovule were evaluated with clearing methods, and mode of seed formation was examined by single seed flow cytometric seed screening (ssFCSS). We further studied pollen stainability to understand effects of pollen quality on seed formation. Results revealed that under extended photoperiod, all cytotypes produced significantly more sexual ovules than in the control, with strongest effects on diploids. The stress treatment affected neither the frequency of seed set nor the proportion of sexual seeds nor pollen quality. Successful seed formation appears to be dependent on balanced maternal: paternal genome contributions. Diploid cytotypes had mostly sexual seed formation, while polyploid cytotypes formed predominantly apomictic seeds. Pollen quality was in hexaploids better than in diploids and tetraploids. These findings confirm our hypothesis that megasporogenesis is triggered by light stress treatments. Comparisons of cytotypes support the hypothesis that ovule development in polyploid plants is less sensitive to prolonged photoperiods and responds to a lesser extent with sexual ovule formation. Polyploids may better buffer environmental stress, which releases the potential for aposporous ovule development from somatic cells, and may facilitate the establishment of apomictic seed formation.

Evolution of small RNA expression following hybridization and allopolyploidization: insights from Spartina species (Poaceae, Chloridoideae)

Differential expression of mi-RNAs targeting developmental processes and progressive downregulation of repeat-associated siRNAs following genome merger and genome duplication in the context of allopolyploid speciation in Spartina. The role of small RNAs on gene expression regulation and genome stability is arousing increased interest and is being explored in various plant systems. In spite of prominence of reticulate evolution and polyploidy that affects the evolutionary history of all plant lineages, very few studies analysed RNAi mechanisms with this respect. Here, we explored small RNAs diversity and expression in the context of recent allopolyploid speciation, using the Spartina system, which offers a unique opportunity to explore the immediate changes following hybridization and genome duplication. Small RNA-Seq analyses were conducted on hexaploid parental species (S. alterniflora and S. maritima), their F1 hybrid S. x townsendii, and the neoallododecaploid S. anglica. We identified 594 miRNAs, 2197 miRNA-target genes, and 3730 repeat-associated siRNAs (mostly targeting Class I/Copia-Ivana- Copia-SIRE and LINEs elements). For both mi- and ra-siRNAs, we detected differential expression patterns following genome merger and genome duplication. These misregulations include non-additive expression of miRNAs in the F1 hybrid and additional changes in the allopolyploid targeting developmental processes. Expression of repeat-associated siRNAs indicates a strengthen of transposable element repression during the allopolyploidization process. Altogether, these results confirm the central role small RNAs play in shaping regulatory changes in naturally formed recent allopolyploids.

Relaxed purifying selection in autopolyploids drives transposable element over-accumulation which provides variants for local adaptation

Polyploidization is frequently associated with increased transposable element (TE) content. However, what drives TE dynamics following whole genome duplication (WGD) and the evolutionary implications remain unclear. Here, we leverage whole-genome resequencing data available for ~300 individuals of Arabidopsis arenosa, a well characterized natural diploid-autotetraploid plant species, to address these questions. Based on 43,176 TE insertions we detect in these genomes, we demonstrate that relaxed purifying selection rather than transposition bursts is the main driver of TE over-accumulation after WGD. Furthermore, the increased pool of TE insertions in tetraploids is especially enriched within or near environmentally responsive genes. Notably, we show that the major flowering-time repressor gene FLC is disrupted by a TE insertion specifically in the rapid-cycling tetraploid lineage that colonized mainland railways. Together, our findings indicate that tetrasomy leads to an enhanced accumulation of genic TE insertions, some of which likely contribute to local adaptation.

Why transposable elements (TEs) accumulate in polyploids and the evolutionary implications remain unclear. Here, the authors show that following whole genome duplication, relaxed purifying selection is the main driver of TE over-accumulation, which provides variants for rapid local adaptation.

Comparatively Barcoded Chromosomes of Brachypodium Perennials Tell the Story of Their Karyotype Structure and Evolution

The Brachypodium genus is an informative model system for studying grass karyotype organization. Previous studies of a limited number of species and reference chromosomes have not provided a comprehensive picture of the enigmatic phylogenetic relationships in the genus. Comparative chromosome barcoding, which enables the reconstruction of the evolutionary history of individual chromosomes and their segments, allowed us to infer the relationships between putative ancestral karyotypes of extinct species and extant karyotypes of current species. We used over 80 chromosome-specific BAC (bacterial artificial chromosome) clones derived from five reference chromosomes of B. distachyon as probes against the karyotypes of twelve accessions representing five diploid and polyploid Brachypodium perennials. The results showed that descending dysploidy is common in Brachypodium and occurs primarily via nested chromosome fusions. Brachypodiumdistachyon was rejected as a putative ancestor for allotetraploid perennials and B. stacei for B. mexicanum. We propose two alternative models of perennial polyploid evolution involving either the incorporation of a putative x = 5 ancestral karyotype with different descending dysploidy patterns compared to B. distachyon chromosomes or hybridization of two x = 9 ancestors followed by genome doubling and descending dysploidy. Details of the karyotype structure and evolution in several Brachypodium perennials are revealed for the first time.

Interaction among ploidy, breeding system and lineage diversification

If particular traits consistently affect rates of speciation and extinction, broad macroevolutionary patterns can be interpreted as consequences of selection at high levels of the biological hierarchy. Identifying traits associated with diversification rates is difficult because of the wide variety of characters under consideration and the statistical challenges of testing for associations from comparative phylogenetic data. Ploidy (diploid vs polyploid states) and breeding system (self-incompatible vs self-compatible states) are both thought to be drivers of differential diversification in angiosperms. We fit 29 diversification models to extensive trait and phylogenetic data in Solanaceae and investigate how speciation and extinction rate differences are associated with ploidy, breeding system, and the interaction between these traits. We show that diversification patterns in Solanaceae are better explained by breeding system and an additional unobserved factor, rather than by ploidy. We also find that the most common evolutionary pathway to polyploidy in Solanaceae occurs via direct breakdown of self-incompatibility by whole genome duplication, rather than indirectly via breakdown followed by polyploidization. Comparing multiple stochastic diversification models that include complex trait interactions alongside hidden states enhances our understanding of the macroevolutionary patterns in plant phylogenies.

TaCER1-1A is involved in cuticular wax alkane biosynthesis in hexaploid wheat and responds to plant abiotic stresses

To protect above-ground plant organs from excessive water loss, their surfaces are coated by waxes. The genes involved in wax formation have been investigated in detail in Arabidopsis but scarcely in crop species. Here, we aimed to isolate and characterize a CER1 enzyme responsible for formation of the very long-chain alkanes present in high concentrations especially during late stages of wheat development. On the basis of comparative wax and transcriptome analyses of various wheat organs, we selected TaCER1-1A as a primary candidate and demonstrated that it was located to the endoplasmic reticulum, the subcellular compartment for wax biosynthesis. A wheat nullisomic-tetrasomic substitution line lacking TaCER1-1A had significantly reduced amounts of C₃₃ alkane, whereas rice plants overexpressing TaCER1-1A showed substantial increases of C₂₅ -C₃₃ alkanes relative to wild type control. Similarly, heterologous expression of TaCER1-1A in Arabidopsis wild type and the cer1 mutant resulted in increased levels of unbranched alkanes, iso-branched alkanes and alkenes. Finally, the expression of TaCER1-1A was found activated by abiotic stresses and abscisic acid treatment, resulting in increased production of alkanes in wheat. Taken together, our results demonstrate that TaCER1-1A plays an important role in wheat wax alkane biosynthesis and involved in responding to drought and other environmental stresses.

Associations between genomic ancestry, genome size and capitula morphology in the invasive meadow knapweed hybrid complex (Centaurea × moncktonii) in eastern North America

Plant invasions are prime opportunities for studying hybridization and the nature of species boundaries, but hybrids also complicate the taxonomic treatment and management of introduced taxa. In this study, we use population genomics to estimate the extent of genomic admixture and test for its association with morphology and genome size in a hybrid complex of knapweeds invasive to North America: meadow knapweed (Centaurea × moncktonii) and its parental species (C. jacea and C. nigra). We sampled 20 populations from New York and Vermont, USA, and used genotyping by sequencing to identify single nucleotide polymorphisms in order to estimate genome-wide ancestry and classify individuals into hybrid genotype classes. We then tested for association between degree of genomic introgression and variation in a subset of traits diagnostic for the parental taxa, namely capitula morphology and monoploid genome size. Genomic clustering revealed two clearly defined lineages, as well as many admixed individuals forming a continuous gradation of introgression. Individual assignments to hybrid genotype classes revealed many advanced generation intercrosses and backcrosses, suggesting introgression has been extensive and unimpeded by strong reproductive barriers between taxa. Variation in capitula traits between the two unadmixed, presumed parental, lineages exhibited continuous, and in some cases transgressive, segregation among introgressed hybrids. Genome size was also divergent between lineages, although advanced generation hybrids had smaller genomes relative to additive expectations. Our study demonstrates deep introgression between the porous genomes of a hybrid invasive species complex. In addition to strong associations among genomic ancestry, genome size and morphology, hybrids expressed more extreme phenotypic values for capitula traits and genome size, indicating transgressive segregation, as well as a bias towards smaller genomes, possibly due to genomic downsizing. Future studies will apply these results to experimentally test how introgression, transgressive segregation and genome size reduction interact to confer invasiveness.