Parallel Evolution of Common Allelic Variants Confers Flowering Diversity in Capsella rubella

Molecular mechanisms of flowering time differentiation revealed by transcriptomic sequencing and de novo analysis in Chinese invasive populations of Ambrosia artemisiifolia

BACKGROUND

Ambrosia artemisiifolia is a highly invasive herb with deleterious effects on public health and agricultural systems. Flowering time in this species has been reported to vary along a latitudinal gradient, which may contribute to local adaptation and invasion success in China. However, the molecular basis for the flowering time differentiation remains unclear.

RESULTS

A common garden experiment confirmed a latitudinal gradient in flowering time among seven Chinese populations. Differentially expressed genes (DEGs) across sampling times and flowering time groups were identified through transcriptome sequencing and analyses of DGE and WGCNA, and were partially annotated to circadian rhythm, light response and hormone response through GO enrichment. By annotating to Flowering Interactive Database (FLOR-ID) and protein-protein interaction (PPI) databases, 53 candidate genes for flowering time differentiation were identified, with 23 of these genes linked to the photoperiod pathway. Additionally, 43 of 53 candidate genes exhibited expression correlated with latitude. Six genes, including FKF1, FT, FUL, MAF2, WNK4 and WNK5, were inferred to promote flowering, while 5 genes, FBH3, FLK, NCL(1), POL2A, and ZHD4, likely repress flowering, based on their expression patterns in relation to latitude and sampling times. Notably, NCL(1), FBH3, MAF2, and FLK may function differently in A. artemisiifolia compared to Arabidopsis thaliana.

CONCLUSIONS

This study identified key candidate genes related to the differentiation of flowering time in Chinese ragweed populations, providing valuable insights into molecular mechanisms of phenological adaptation and invasive success of ragweed.

The molecular basis of phenotypic evolution: beyond the usual suspects

It has been well documented that mutations in coding DNA or cis-regulatory elements underlie natural phenotypic variation in many organisms. However, the development of sophisticated functional tools in recent years in a wide range of traditionally non-model systems have revealed many 'unusual suspects' in the molecular bases of phenotypic evolution, including upstream open reading frames (uORFs), cryptic splice sites, and small RNAs. Furthermore, large-scale genome sequencing, especially long-read sequencing, has identified a cornucopia of structural variation underlying phenotypic divergence and elucidated the composition of supergenes that control complex multi-trait polymorphisms. In this review article we highlight recent studies that demonstrate this great diversity of molecular mechanisms producing adaptive genetic variation and the panoply of evolutionary paths leading to the 'grandeur of life'.

Strength in numbers: An isoform variety of homogalacturonan modifying enzymes may contribute to pollen tube fitness

Pectin is a major component of the cell wall in land plants. It plays crucial roles in cell wall assembly, cell growth, shaping, and signaling. The relative abundance of pectin in the cell wall is particularly high in rapidly growing organ regions and cell types. Homogalacturonan (HG), a polymer of 1,4-linked α-D-galacturonic acid, is a major pectin constituent in growing and dividing plant cells. In pollen tubes, an extremely rapidly growing cell type, HG is secreted at and inserted into the apical cell wall and is subject to further modification in muro by HG modifying enzymes (HGMEs). These enzymes, including pectin esterases and depolymerases, have multiple isoforms, some of which are specifically expressed in pollen. Given the importance of pectin chemistry for the fitness of pollen tubes, it is of interest to interrogate the potentially crucial roles these isoforms play in pollen germination and elongation. It is hypothesized that different HGME isoforms, through their action on apoplastic HG, may generate differential methylation and acetylation patterns endowing HG polysaccharides with specific, spatially and temporally varying properties that lead to a fine-tuned pattern of cell wall modification. In addition, these isoforms may be differentially activated and/or inhibited depending on the local conditions that may vary at subcellular resolution. In this Update we review the different HGME isoforms identified in recent years in Arabidopsis thaliana and postulate that the multiplicity of these isoforms may allow for specialized substrate recognition and conditional activation, leading to a sophisticated regulation scheme exemplified in the process that governs the dynamic properties of the cell wall in pollen tube growth.

Evolutionary Dynamics of FLC-like MADS-Box Genes in Brassicaceae

MADS-box genes encode transcription factors that play important roles in the development and evolution of plants. There are more than a dozen clades of MADS-box genes in angiosperms, of which those with functions in the specification of floral organ identity are especially well-known. From what has been elucidated in the model plant Arabidopsis thaliana, the clade of FLC-like MADS-box genes, comprising FLC-like genes sensu strictu and MAF-like genes, are somewhat special among the MADS-box genes of plants since FLC-like genes, especially MAF-like genes, show unusual evolutionary dynamics, in that they generate clusters of tandemly duplicated genes. Here, we make use of the latest genomic data of Brassicaceae to study this remarkable feature of the FLC-like genes in a phylogenetic context. We have identified all FLC-like genes in the genomes of 29 species of Brassicaceae and reconstructed the phylogeny of these genes employing a Maximum Likelihood method. In addition, we conducted selection analyses using PAML. Our results reveal that there are three major clades of FLC-like genes in Brassicaceae that all evolve under purifying selection but with remarkably different strengths. We confirm that the tandem arrangement of MAF-like genes in the genomes of Brassicaceae resulted in a high rate of duplications and losses. Interestingly, MAF-like genes also seem to be prone to transposition. Considering the role of FLC-like genes sensu lato (s.l.) in the timing of floral transition, we hypothesize that this rapid evolution of the MAF-like genes was a main contributor to the successful adaptation of Brassicaceae to different environments.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

An enhancing effect attributed to a nonsynonymous mutation in SOYBEAN SEED SIZE 1, a SPINDLY-like gene, is exploited in soybean domestication and improvement

Soybean (Glycine max) was domesticated from its wild relative Glycine soja. One-hundred-seed weight is one of the most important domesticated traits determining soybean yield; however, its underlying genetic basis remains elusive. We characterized a soybean seed size 1 (sss1) mutant featuring large seeds compared to its wild-type background. Positional cloning revealed that the candidate gene GmSSS1 encoded a SPINDLY homolog and was co-located in a well-identified quantitative trait locus (QTL)-rich region on chromosome 19. Knocking out GmSSS1 resulted in small seeds, while overexpressing GmSSS1/Gmsss1 induced large seeds. Modulating GmSSS1/Gmsss1 in transgenic plants can positively influence cell expansion and cell division. Relative to GmSSS1, one mutation leading to an E to Q substitution at the 182^nd residue in Gmsss1 conferred an enhancing effect on seed weight. GmSSS1 underwent diversification in wild-type and cultivated soybean, and the alleles encoding the Gmsss1-type substitution of 182^nd -Q, which originated along the central and downstream parts of the Yellow River, were selected and expanded during soybean domestication and improvement. We cloned the causative gene for the sss1 mutant, which is linked with a seed weight QTL, identified an elite allele of this gene for increasing seed weight, and provided new insights into soybean domestication and breeding.

The Boechera model system for evolutionary ecology

Model systems in biology expand the research capacity of individuals and the community. Closely related to Arabidopsis, the genus Boechera has emerged as an important ecological model owing to the ability to integrate across molecular, functional, and eco-evolutionary approaches. Boechera species are broadly distributed in relatively undisturbed habitats predominantly in western North America and provide one of the few experimental systems for identification of ecologically important genes through genome-wide association studies and investigations of selection with plants in their native habitats. The ecologically, evolutionarily, and agriculturally important trait of apomixis (asexual reproduction via seeds) is common in the genus, and field experiments suggest that abiotic and biotic environments shape the evolution of sex. To date, population genetic studies have focused on the widespread species B. stricta, detailing population divergence and demographic history. Molecular and ecological studies show that balancing selection maintains genetic variation in ~10% of the genome, and ecological trade-offs contribute to complex trait variation for herbivore resistance, flowering phenology, and drought tolerance. Microbiome analyses have shown that host genotypes influence leaf and root microbiome composition, and the soil microbiome influences flowering phenology and natural selection. Furthermore, Boechera offers numerous opportunities for investigating biological responses to global change. In B. stricta, climate change has induced a shift of >2 weeks in the timing of first flowering since the 1970s, altered patterns of natural selection, generated maladaptation in previously locally-adapted populations, and disrupted life history trade-offs. Here we review resources and results for this eco-evolutionary model system and discuss future research directions.

What can cold-induced transcriptomes of Arctic Brassicaceae tell us about the evolution of cold tolerance?

Little is known about the evolution of cold tolerance in polar plant species and how they differ from temperate relatives. To gain insight into their biology and the evolution of cold tolerance, we compared the molecular basis of cold response in three Arctic Brassicaceae species. We conducted a comparative time series experiment to examine transcriptional responses to low temperature. RNA was sampled at 22°C, and after 3, 6, and 24 at 2°C. We then identified sets of genes that were differentially expressed in response to cold and compared them between species, as well as to published data from the temperate Arabidopsis thaliana. Most differentially expressed genes were species‐specific, but a significant portion of the cold response was also shared among species. Among thousands of differentially expressed genes, ~200 were shared among the three Arctic species and A. thaliana, while ~100 were exclusively shared among the three Arctic species. Our results show that cold response differs markedly between Arctic Brassicaceae species, but probably builds on a conserved basis found across the family. They also confirm that highly polygenic traits such as cold tolerance may show little repeatability in their patterns of adaptation.

Fine mapping of qPH9, a major quantitative trait locus, responsible for plant height in foxtail millet [Setaria italica (L.) P. Beauv.]

Plant height is vital for crop yield by influencing plant architecture and resistance to lodging. Although lots of quantitative trait loci (QTLs) controlling plant height had been mapped in foxtail millet, their contributions to phenotypic variation were generally small and mapping regions were relatively large, indicating the difficult application in molecular breeding using marker-assisted selection (MAS). In the present paper, a total of 23 QTLs involving in 15 traits were identified via a high-density Bin map containing 3024 Bin markers with an average distance of 0.48 cM through an F2 population. Among them, qPH9, with a large phenotypic variation explained (51.6%) related to plant height, was one of the major QTLs. Furthermore, qPH9 was repeatedly detected in multi-environments under field conditions using two new developed F2 populations from the same F1 plant, and was narrowed down to a smaller interval of 281 kb using 1024 recessive F2 individuals from the same F1 plant. Finally, we found that there was an extremely significant correlation between marker MRI1016 and plant height, and further speculated that Seita.9G088900 and Seita.9G089700 could be key candidates of qPH9. This study laid an important foundation for the cloning of qPH9 and molecular breeding of dwarf varieties via MAS.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01261-w.

Molecular Parallelism Underlies Convergent Highland Adaptation of Maize Landraces

Convergent phenotypic evolution provides some of the strongest evidence for adaptation. However, the extent to which recurrent phenotypic adaptation has arisen via parallelism at the molecular level remains unresolved, as does the evolutionary origin of alleles underlying such adaptation. Here, we investigate genetic mechanisms of convergent highland adaptation in maize landrace populations and evaluate the genetic sources of recurrently selected alleles. Population branch excess statistics reveal substantial evidence of parallel adaptation at the level of individual single-nucleotide polymorphism (SNPs), genes, and pathways in four independent highland maize populations. The majority of convergently selected SNPs originated via migration from a single population, most likely in the Mesoamerican highlands, while standing variation introduced by ancient gene flow was also a contributor. Polygenic adaptation analyses of quantitative traits reveal that alleles affecting flowering time are significantly associated with elevation, indicating the flowering time pathway was targeted by highland adaptation. In addition, repeatedly selected genes were significantly enriched in the flowering time pathway, indicating their significance in adapting to highland conditions. Overall, our study system represents a promising model to study convergent evolution in plants with potential applications to crop adaptation across environmental gradients.

Reinventing the wheel? Reassessing the roles of gene flow, sorting and convergence in repeated evolution

Biologists have long been intrigued by apparently predictable and repetitive evolutionary trajectories inferred across a variety of lineages and systems. In recent years, high-throughput sequencing analyses have started to transform our understanding of such repetitive shifts. While researchers have traditionally categorized such shifts as either "convergent" or "parallel," based on relatedness of the lineages involved, emerging genomic insights provide an opportunity to better describe the actual evolutionary mechanisms at play. A synthesis of recent genomic analyses confirms that convergence is the predominant driver of repetitive evolution among species, whereas repeated sorting of standing variation is the major driver of repeated shifts within species. However, emerging data reveal numerous notable exceptions to these expectations, with recent examples of de novo mutations underpinning convergent shifts among even very closely related lineages, while repetitive sorting processes have occurred among even deeply divergent taxa, sometimes via introgression. A number of very recent analyses have found evidence for both processes occurring on different scales within taxa. We suggest that the relative importance of convergent versus sorting processes depends on the interplay between gene flow among populations, and phylogenetic relatedness of the lineages involved.

Parallel evolution of apetalous lineages within the buttercup family (Ranunculaceae): outward expansion of AGAMOUS1, rather than disruption of APETALA3-3

Complete loss of petals, or becoming apetalous, has occurred independently in many flowering plant lineages. However, the mechanisms underlying the parallel evolution of naturally occurring apetalous lineages remain largely unclear. Here, by sampling representatives of all nine apetalous genera/tribes of the family Ranunculaceae and conducting detailed morphological, expression, molecular evolutionary and functional studies, we investigate the mechanisms underlying parallel petal losses. We found that while non-expression/downregulation of the petal identity gene APETALA3-3 (AP3-3) is tightly associated with complete petal losses, disruptions of the AP3-3 orthologs were unlikely to be the real causes for the parallel evolution of apetalous lineages. We also found that, compared with their close petalous relatives, naturally occurring apetalous taxa usually bear slightly larger numbers of stamens, whereas the number of sepals remains largely unchanged, suggestive of petal-to-stamen rather than petal-to-sepal transformations. In addition, in the recently originated apetalous genus Enemion, the petal-to-stamen transformations have likely been caused by the mutations that led to the elevation and outward expansion of the expression of the C-function gene, AGAMOUS1 (AG1). Our results not only provide a general picture of parallel petal losses within the Ranunculaceae but also help understand the mechanisms underlying the independent originations of other apetalous lineages.

Less Is More, Natural Loss-of-Function Mutation Is a Strategy for Adaptation

Gene gain and loss are crucial factors that shape the evolutionary success of diverse organisms. In the past two decades, more attention has been paid to the significance of gene gain through gene duplication or de novo genes. However, gene loss through natural loss-of-function (LoF) mutations, which is prevalent in the genomes of diverse organisms, has been largely ignored. With the development of sequencing techniques, many genomes have been sequenced across diverse species and can be used to study the evolutionary patterns of gene loss. In this review, we summarize recent advances in research on various aspects of LoF mutations, including their identification, evolutionary dynamics in natural populations, and functional effects. In particular, we discuss how LoF mutations can provide insights into the minimum gene set (or the essential gene set) of an organism. Furthermore, we emphasize their potential impact on adaptation. At the genome level, although most LoF mutations are neutral or deleterious, at least some of them are under positive selection and may contribute to biodiversity and adaptation. Overall, we highlight the importance of natural LoF mutations as a robust framework for understanding biological questions in general.

Population genomics of parallel adaptation

Parallel evolution is one of the striking patterns in nature. The presence of repeated evolution of the same phenotypes, suites of traits, and adaptations suggests a strong role for natural selection in shaping biological diversity. The reasoning is straightforward: each instance of repeated evolution makes it less likely that these features evolved neutrally or due to stochastic forces in each population or species. With the growing sequencing capability, we are now poised to examine the genetic basis of parallel evolution in model and nonmodel systems. On pages 4102-4117 of this issue of Molecular Ecology, van Boheemen and Hodgins (2020) provide an exemplar study of this kind, using common ragweed (Ambrosia artemisiifolia; Figure 1a). Their study is noteworthy and ambitious in many respects, and we think will serve as a model for studying parallel adaptation, even in nonmodel species.

Evolutionary processes from the perspective of flowering time diversity

Although it is well appreciated that genetic studies of flowering time regulation have led to fundamental advances in the fields of molecular and developmental biology, the ways in which genetic studies of flowering time diversity have enriched the field of evolutionary biology have received less attention despite often being equally profound. Because flowering time is a complex, environmentally responsive trait that has critical impacts on plant fitness, crop yield, and reproductive isolation, research into the genetic architecture and molecular basis of its evolution continues to yield novel insights into our understanding of domestication, adaptation, and speciation. For instance, recent studies of flowering time variation have reconstructed how, when, and where polygenic evolution of phenotypic plasticity proceeded from standing variation and de novo mutations; shown how antagonistic pleiotropy and temporally varying selection maintain polymorphisms in natural populations; and provided important case studies of how assortative mating can evolve and facilitate speciation with gene flow. In addition, functional studies have built detailed regulatory networks for this trait in diverse taxa, leading to new knowledge about how and why developmental pathways are rewired and elaborated through evolutionary time.

Evolution of Yin and Yang isoforms of a chromatin remodeling subunit precedes the creation of two genes

Genes can encode multiple isoforms, broadening their functions and providing a molecular substrate to evolve phenotypic diversity. Evolution of isoform function is a potential route to adapt to new environments. Here we show that de novo, beneficial alleles in the nurf-1 gene became fixed in two laboratory lineages of C. elegans after isolation from the wild in 1951, before methods of cryopreservation were developed. nurf-1 encodes an ortholog of BPTF, a large (>300 kD) multidomain subunit of the NURF chromatin remodeling complex. Using CRISPR-Cas9 genome editing and transgenic rescue, we demonstrate that in C. elegans, nurf-1 has split into two, largely non-overlapping isoforms (NURF-1.D and NURF-1.B, which we call Yin and Yang, respectively) that share only two of 26 exons. Both isoforms are essential for normal gametogenesis but have opposite effects on male/female gamete differentiation. Reproduction in hermaphrodites, which involves production of both sperm and oocytes, requires a balance of these opposing Yin and Yang isoforms. Transgenic rescue and genetic position of the fixed mutations suggest that different isoforms are modified in each laboratory strain. In a related clade of Caenorhabditis nematodes, the shared exons have duplicated, resulting in the split of the Yin and Yang isoforms into separate genes, each containing approximately 200 amino acids of duplicated sequence that has undergone accelerated protein evolution following the duplication. Associated with this duplication event is the loss of two additional nurf-1 transcripts, including the long-form transcript and a newly identified, highly expressed transcript encoded by the duplicated exons. We propose these lost transcripts are non-functional side products necessary to transcribe the Yin and Yang transcripts in the same cells. Our work demonstrates how gene sharing, through the production of multiple isoforms, can precede the creation of new, independent genes.

Transposable elements drive rapid phenotypic variation in Capsella rubella

The mechanisms underlying rapid adaptation to changing environments in species with reduced genetic variation, referred to as the “genetic paradox of invasion,” are unknown. We report that transposable elements (TEs) are highly enriched in the gene promoter regions of Capsella rubella compared with its outcrossing sister species Capsella grandiflora. We also show that a number of polymorphic TEs in C. rubella are associated with changes in gene expression. Frequent TE insertions at FLOWERING LOCUS C of C. rubella affect flowering-time variation, an important life history trait correlated with fitness. These results indicate that TE insertions drive rapid phenotypic variation, which could potentially help adapting to novel environments in species with limited genetic variation.

Rapid phenotypic changes in traits of adaptive significance are crucial for organisms to thrive in changing environments. How such phenotypic variation is achieved rapidly, despite limited genetic variation in species that experience a genetic bottleneck is unknown. Capsella rubella, an annual and inbreeding forb (Brassicaceae), is a great system for studying this basic question. Its distribution is wider than those of its congeneric species, despite an extreme genetic bottleneck event that severely diminished its genetic variation. Here, we demonstrate that transposable elements (TEs) are an important source of genetic variation that could account for its high phenotypic diversity. TEs are (i) highly enriched in C. rubella compared with its outcrossing sister species Capsella grandiflora, and (ii) 4.2% of polymorphic TEs in C. rubella are associated with variation in the expression levels of their adjacent genes. Furthermore, we show that frequent TE insertions at FLOWERING LOCUS C (FLC) in natural populations of C. rubella could explain 12.5% of the natural variation in flowering time, a key life history trait correlated with fitness and adaptation. In particular, we show that a recent TE insertion at the 3′ UTR of FLC affects mRNA stability, which results in reducing its steady-state expression levels, to promote the onset of flowering. Our results highlight that TE insertions can drive rapid phenotypic variation, which could potentially help with adaptation to changing environments in a species with limited standing genetic variation.

The Functional Change and Deletion of FLC Homologs Contribute to the Evolution of Rapid Flowering in Boechera stricta

Differences in the timing of vegetative-to-reproductive phase transition have evolved independently and repeatedly in different plant species. Due to their specific biological functions and positions in pathways, some genes are important targets of repeated evolution - independent mutations on these genes caused the evolution of similar phenotypes in distantly related organisms. While many studies have investigated these genes, it remains unclear how gene duplications influence repeated phenotypic evolution. Here we characterized the genetic architecture underlying a novel rapid-flowering phenotype in Boechera stricta and investigated the candidate genes BsFLC1 and BsFLC2. The expression patterns of BsFLC1 suggested its function in flowering time suppression, and the deletion of BsFLC1 is associated with rapid flowering and loss of vernalization requirement. In contrast, BsFLC2 did not appear to be associated with flowering and had accumulated multiple amino acid substitutions in the relatively short evolutionary timeframe after gene duplication. These non-synonymous substitutions greatly changed the physicochemical properties of the original amino acids, concentrated non-randomly near a protein-interacting domain, and had greater substitution rate than synonymous changes. Here we suggested that, after recent gene duplication of the FLC gene, the evolution of rapid phenology was made possible by the change of BsFLC2 expression pattern or protein sequences and the deletion of BsFLC1.