Distinct types of selection and genetic architecture shape molecular variation during the domestication of vegetable crops.

Humans select vegetable crops with desirable traits via a complex evolutionary process called domestication, generating a variety of cultivars worldwide. With advances in sequencing technologies, genomic scans for "signatures of selection" are widely used to identify target loci of selection. In the early phases of domestication, humans tended to favor similar sets of phenotypes in diverse crops, resulting in "domestication syndrome" and parallel evolution in multiple species. Subsequently, adaptation to distinct environments or different consumer preferences have diversified crop cultivars. Here, we review molecular and population genetic studies on genes affecting trait evolution during this complex process. We emphasize that, depending on interactions among different types of selection (directional selection within or divergent selection between groups), the genetic architecture of the target trait (Mendelian or polygenic), and the origin of the causal variant (new mutation or standing variation), the resulting molecular patterns of variation can be highly diverse. Situations in which the typical hard selective sweep model could be applied may be limited. Therefore, it is crucial to obtain a thorough understanding of the target species' historical, environmental, and ecological contexts.

The Past Contribution and Future Fate of Genetic Variants under Climate Change in an Island Population of Musa itinerans.

AbstractGenetic variation within species is crucial for sessile species to adapt to novel environments when facing dramatic climate changes. However, the debate continues whether standing ancestral variation adaptive to current environmental variability is sufficient to guarantee future suitability. Using wild banana Musa itinerans, we investigated the relative contribution of standing ancestral variation versus new mutations to environmental adaptation and inferred their future fate. On the continental island of Taiwan, local populations immigrated from the Southeast Asian continent during the ice age and have been isolated since then. This allows the classification of genetic variants into standing ancestral variation (polymorphic in Taiwan and the continent) and new mutations (polymorphic only in Taiwan). For temperature-related variables where Taiwan is mainly within the ancestral climatic range, standing ancestral variation had a slightly stronger association than new mutations. New mutations were more important for precipitation-related variables, where northeastern Taiwan had much more winter rainfall than most of continental Southeast Asia. Upon future climate change, new mutations showed higher genetic offset in regions of abrupt transition between allele frequency and local environments, suggesting their greater spatial heterogeneity of future vulnerability.

Transcriptome-wide association study coupled with eQTL analysis reveals the genetic connection between gene expression and flowering time in Arabidopsis.

Genome-wide association study (GWAS) has improved our understanding of complex traits, but challenges exist in distinguishing causation versus association caused by linkage disequilibrium. Instead, transcriptome-wide association studies (TWAS) detect direct associations between expression levels and phenotypic variations, providing an opportunity to better prioritize candidate genes. To assess the feasibility of TWAS, we investigated the association between transcriptomes, genomes, and various traits in Arabidopsis, including flowering time. The associated genes formerly known to regulate growth allometry or metabolite production were first identified by TWAS. Next, for flowering time, six TWAS-newly identified genes were functionally validated. Analysis of the expression quantitative trait locus (eQTL) further revealed a trans-regulatory hotspot affecting the expression of several TWAS-identified genes. The hotspot covers the FRIGIDA (FRI) gene body, which possesses multiple haplotypes differentially affecting the expression of downstream genes, such as FLOWERING LOCUS C (FLC) and SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1). We also revealed multiple independent paths towards the loss of function of FRI in natural accessions. Altogether, this study demonstrates the potential of combining TWAS with eQTL analysis to identify important regulatory modules of FRI-FLC-SOC1 for quantitative traits in natural populations.

Demographic history and distinct selection signatures of two domestication genes in mungbean.

Domestication is the long and complex process underlying the evolution of crops, in which artificial directional selection transformed wild progenitors into the desired form, affecting genomic variation and leaving traces of selection at targeted loci. However, whether genes controlling important domestication traits follow the same evolutionary pattern expected under the standard selective sweep model remains unclear. With whole-genome resequencing of mungbean (Vigna radiata), we investigated this issue by resolving its global demographic history and targeted dissection of the molecular footprints of genes underlying 2 key traits representing different stages of domestication. Mungbean originated in Asia, and the Southeast Asian wild population migrated to Australia about 50 thousand generations ago. Later in Asia, the cultivated form diverged from the wild progenitor. We identified the gene associated with the pod shattering resistance trait, VrMYB26a, with lower expression across cultivars and reduced polymorphism in the promoter region, reflecting a hard selective sweep. On the other hand, the stem determinacy trait was associated with VrDet1. We found that 2 ancient haplotypes of this gene have lower gene expression and exhibited intermediate frequencies in cultivars, consistent with selection favoring independent haplotypes in a soft selective sweep. In mungbean, contrasting signatures of selection were identified from the detailed dissection of 2 important domestication traits. The results suggest complex genetic architecture underlying the seemingly simple process of directional artificial selection and highlight the limitations of genome-scan methods relying on hard selective sweeps.

The role of gene presence-absence variations on genetic incompatibility in Asian rice.

Genetic incompatibilities are widespread between species. However, it remains unclear whether they all originated after population divergence as suggested by the Bateson-Dobzhansky-Muller model, and if not, what is their prevalence and distribution within populations. The gene presence-absence variations (PAVs) provide an opportunity for investigating gene-gene incompatibility. Here, we searched for the repulsion of coexistence between gene PAVs to identify the negative interaction of gene functions separately in two Oryza sativa subspecies. Many PAVs are involved in subspecies-specific negative epistasis and segregate at low-to-intermediate frequencies in focal subspecies but at low or high frequencies in the other subspecies. Incompatible PAVs are enriched in two functional groups, defense response and protein phosphorylation, which are associated with plant immunity and consistent with autoimmunity being a known mechanism of hybrid incompatibility in plants. Genes in the two enriched functional groups are older and seldom directly interact with each other. Instead, they interact with other younger gene PAVs with diverse functions. Our results illustrate the landscape of genetic incompatibility at gene PAVs in rice, where many incompatible pairs have already segregated as polymorphisms within subspecies, and many are novel negative interactions between older defense-related genes and younger genes with diverse functions.

Environment as a limiting factor of the historical global spread of mungbean.

While the domestication process has been investigated in many crops, the detailed route of cultivation range expansion and factors governing this process received relatively little attention. Here, using mungbean (Vigna radiata var. radiata) as a test case, we investigated the genomes of more than 1000 accessions to illustrate climatic adaptation's role in dictating the unique routes of cultivation range expansion. Despite the geographical proximity between South and Central Asia, genetic evidence suggests mungbean cultivation first spread from South Asia to Southeast, East and finally reached Central Asia. Combining evidence from demographic inference, climatic niche modeling, plant morphology, and records from ancient Chinese sources, we showed that the specific route was shaped by the unique combinations of climatic constraints and farmer practices across Asia, which imposed divergent selection favoring higher yield in the south but short-season and more drought-tolerant accessions in the north. Our results suggest that mungbean did not radiate from the domestication center as expected purely under human activity, but instead, the spread of mungbean cultivation is highly constrained by climatic adaptation, echoing the idea that human commensals are more difficult to spread through the south-north axis of continents.

Mungbean, also known as green gram, is an important crop plant in China, India, the Philippines and many other countries across Asia. Archaeological evidence suggests that humans first cultivated mungbeans from wild relatives in India over 4,000 years ago. However, it remains unclear how cultivation has spread to other countries and whether human activity alone dictated the route of the cultivated mungbean's expansion across Asia, or whether environmental factors, such as climate, also had an impact. To understand how a species of plant has evolved, researchers may collect specimens from the wild or from cultivated areas. Each group of plants of the same species they collect in a given location at a single point in time is known collectively as an accession. Ong et al. used a combination of genome sequencing, computational modelling and plant biology approaches to study more than 1,000 accessions of cultivated mungbean and trace the route of the crop's expansion across Asia. The data support the archaeological evidence that mungbean cultivation first spread from South Asia to Southeast Asia, then spread northwards to East Asia and afterwards to Central Asia. Computational modelling of local climates and the physical characteristics of different mungbean accessions suggest that the availability of water in the local area likely influenced the route. Specifically, accessions from arid Central Asia were better adapted to drought conditions than accessions from wetter South Asia. However, these drought adaptations decreased the yield of the plants, which may explain why the more drought tolerant accessions have not been widely grown in wetter parts of Asia. This study shows that human activity has not solely dictated where mungbean has been cultivated. Instead, both human activity and the various adaptations accessions evolved in response to their local environments shaped the route the crop took across Asia. In the future these findings may help plant breeders to identify varieties of mungbean and other crops with drought tolerance and other potentially useful traits for agriculture.

Modeling of Flowering Time in Vigna radiata with Artificial Image Objects, Convolutional Neural Network and Random Forest.

Flowering time is an important target for breeders in developing new varieties adapted to changing conditions. In this work, a new approach is proposed in which the SNP markers influencing time to flowering in mung bean are selected as important features in a random forest model. The genotypic and weather data are encoded in artificial image objects, and a model for flowering time prediction is constructed as a convolutional neural network. The model uses weather data for only a limited time period of 5 days before and 20 days after planting and is capable of predicting the time to flowering with high accuracy. The most important factors for model solution were identified using saliency maps and a Score-CAM method. Our approach can help breeding programs harness genotypic and phenotypic diversity to more effectively produce varieties with a desired flowering time.

Ecological factors influence balancing selection on leaf chemical profiles of a wildflower.

Balancing selection is frequently invoked as a mechanism that maintains variation within and across populations. However, there are few examples of balancing selection operating on loci underpinning complex traits, which frequently display high levels of variation. We investigated mechanisms that may maintain variation in a focal polymorphism-leaf chemical profiles of a perennial wildflower (Boechera stricta, Brassicaceae)-explicitly interrogating multiple ecological and genetic processes including spatial variation in selection, antagonistic pleiotropy and frequency-dependent selection. A suite of common garden and greenhouse experiments showed that the alleles underlying variation in chemical profile have contrasting fitness effects across environments, implicating two ecological drivers of selection on chemical profile: herbivory and drought. Phenotype-environment associations and molecular genetic analyses revealed additional evidence of past selection by these drivers. Together, these data are consistent with balancing selection on chemical profile, probably caused by pleiotropic effects of secondary chemical biosynthesis genes on herbivore defence and drought response.

The ecological, genetic and genomic architecture of local adaptation and population differentiation in Boechera stricta.

Differential local adaptation restricts gene flow between populations inhabiting distinct environments, resulting in isolation by adaptation. In addition to the statistical inferences of genotype-environment associations, an integrative approach is needed to investigate the effect of local adaptation on population divergence at the ecological, genetic and genomic scale. Here, we combine reciprocal transplant, genome-environment association and QTL mapping to investigate local adaptation in Boechera stricta (Drummond's rockcress). With reciprocal transplant experiment, we found local genetic groups exhibit phenotypic characteristics corresponding to the distinct selection forces from different water availability. At the genetic level, the local allele of a major fitness QTL confers higher and sturdier flowering stalks, maximizing the fecundity fitness component under sufficient water supply, and its genetic variation is associated with precipitation across the landscape. At the genomewide scale, we further showed that multiple loci associated with precipitation are highly differentiated between genetic groups, suggesting that local adaptation has a widespread effect on reducing gene flow. This study provides one of the few comprehensive examples demonstrating how local adaptation facilitates population divergence at the trait, gene and genome level.

Dynamical climatic model for time to flowering in Vigna radiata.

BACKGROUND: Phenology data collected recently for about 300 accessions of Vigna radiata (mungbean) is an invaluable resource for investigation of impacts of climatic factors on plant development. RESULTS: We developed a new mathematical model that describes the dynamic control of time to flowering by daily values of maximal and minimal temperature, precipitation, day length and solar radiation. We obtained model parameters by adaptation to the available experimental data. The models were validated by cross-validation and used to demonstrate that the phenology of adaptive traits, like flowering time, is strongly predicted not only by local environmental factors but also by plant geographic origin and genotype. CONCLUSIONS: Of local environmental factors maximal temperature appeared to be the most critical factor determining how faithfully the model describes the data. The models were applied to forecast time to flowering of accessions grown in Taiwan in future years 2020-2030.

Genome-wide association study in accessions of the mini-core collection of mungbean (Vigna radiata) from the World Vegetable Gene Bank (Taiwan).

BACKGROUND: Mungbean (Vigna radiata (L.) R. Wilczek, or green gram) is important tropical and sub-tropical legume and a rich source of dietary protein and micronutrients. In this study we employ GWAS to examine the genetic basis of variation in several important traits in mungbean, using the mini-core collection established by the World Vegetable Center, which includes 296 accessions that represent the major market classes. This collection has been grown in a common field plot in southern European part of Russia in 2018. RESULTS: We used 5041 SNPs in 293 accessions that passed strict filtering for genetic diversity, linkage disequilibrium, population structure and GWAS analysis. Polymorphisms were distributed among all chromosomes, but with variable density. Linkage disequilibrium decayed in approximately 105 kb. Four distinct subgroups were identified within 293 accessions with 70% of accessions attributed to one of the four populations. By performing GWAS on the mini-core collection we have found several loci significantly associated with two important agronomical traits. Four SNPs associated with possibility of maturation in Kuban territory of Southern Russia in 2018 were identified within a region of strong linkage which contains genes encoding zinc finger A20 and an AN1 domain stress-associated protein. CONCLUSIONS: The core collection of mungbean established by the World Vegetable Center is a valuable resource for mungbean breeding. The collection has been grown in southern European part of Russia in 2018 under incidental stresses caused by abnormally hot weather and different photoperiod. We have found several loci significantly associated with color of hypocotyl and possibility of maturation under these stressful conditions. SNPs associated with possibility of maturation localize to a region on chromosome 2 with strong linkage, in which genes encoding zinc finger A20 and AN1 domain stress associated protein (SAP) are located. Phenotyping of WorldVeg collection for maturation traits in temperate climatic locations is important as phenology remains a critical breeding target for mungbean. As demand rises for mungbean, production in temperate regions with shorter growing seasons becomes crucial to keep up with needs. Uncovering SNPs for phenology traits will speed breeding efforts.

Long-read bitter gourd (Momordica charantia) genome and the genomic architecture of nonclassic domestication.

The genetic architecture of quantitative traits is determined by both Mendelian and polygenic factors, yet classic examples of plant domestication focused on selective sweep of newly mutated Mendelian genes. Here we report the chromosome-level genome assembly and the genomic investigation of a nonclassic domestication example, bitter gourd (Momordica charantia), an important Asian vegetable and medicinal plant of the family Cucurbitaceae. Population resequencing revealed the divergence between wild and South Asian cultivars about 6,000 y ago, followed by the separation of the Southeast Asian cultivars about 800 y ago, with the latter exhibiting more extreme trait divergence from wild progenitors and stronger signs of selection on fruit traits. Unlike some crops where the largest phenotypic changes and traces of selection happened between wild and cultivar groups, in bitter gourd large differences exist between two regional cultivar groups, likely reflecting the distinct consumer preferences in different countries. Despite breeding efforts toward increasing female flower proportion, a gynoecy locus exhibits complex patterns of balanced polymorphism among haplogroups, with potential signs of selective sweep within haplogroups likely reflecting artificial selection and introgression from cultivars back to wild accessions. Our study highlights the importance to investigate such nonclassic example of domestication showing signs of balancing selection and polygenic trait architecture in addition to classic selective sweep in Mendelian factors.

The climatic association of population divergence and future extinction risk of Solanum pimpinellifolium.

Under intraspecific differentiation driven by differential climatic adaptation, it may be expected that intraspecific genetic groups occur at distinct environments. Populations occupying different niches may therefore differ in their ability to cope with climate change. Here, we addressed this hypothesis with a wild tomato, Solanum pimpinellifolium. This species is distributed from the west side of Andes to the coastal region in Peru and Ecuador and occupies a wide environmental diversity. This environmental diversity is related to the genetic structure of the species providing an ideal material to investigate the isolation by environment hypothesis. While previous hypothesis stated that S. pimpinellifolium originated from northern Peru and migrated northwards and southwards, our results support that S. pimpinellifolium originated from Ecuador and expanded to northern and southern Peru, and during this process, the niche space of S. pimpinellifolium became more associated with cold and drought. We further predicted its fate under anthropogenic climate change. According to our predictions, the northern group will maintain its current extent or even expand to the entire western region of Ecuador. In contrast, we predicted low habitat suitability for the southern group which could potentially lead to the shrinkage of its distribution. In conclusion, we revealed the distinct fates among the differentiated populations driven by environment under global warming conditions.

Quantitative Trait Locus Analysis of Leaf Morphology Indicates Conserved Shape Loci in Grapevine.

Leaf shape in plants plays important roles in water use, canopy structure, and physiological tolerances to abiotic stresses; all important traits for the future development and sustainability of grapevine cultivation. Historically, researchers have used ampelography, the study of leaf shape in grapevines, to differentiate Vitis species and cultivars based on finite leaf attributes. However, ampelographic measurements have limitations and new methods for quantifying shape are now available. We paired an analysis of finite trait attributes with a 17-point landmark survey and generalized Procrustes analysis (GPA) to reconstruct grapevine leaves digitally from five interspecific hybrid mapping families. Using the reconstructed leaves, we performed three types of quantitative trait loci (QTL) analyses to determine the genetic architecture that defines leaf shape. In the first analysis, we compared several important ampelographic measurements as finite trait QTL. In the second and third analyses, we identified significant shape variation via principal components analysis (PCA) and using a multivariate least squares interval mapping (MLSIM) approach. In total, we identified 271 significant QTL across the three measures of leaf shape and identified specific QTL hotspots in the grape genome which appear to drive major aspects of grapevine leaf shape.

Correction to: Ancient polymorphisms contribute to genome-wide variation by long-term balancing selection and divergent sorting in Boechera stricta.

Following publication of the original article [1], the authors reported that the Availability of data and materials section required updating. The updated text reads as follows.

Ancient polymorphisms contribute to genome-wide variation by long-term balancing selection and divergent sorting in Boechera stricta.

BACKGROUND: Genomic variation is widespread, and both neutral and selective processes can generate similar patterns in the genome. These processes are not mutually exclusive, so it is difficult to infer the evolutionary mechanisms that govern population and species divergence. Boechera stricta is a perennial relative of Arabidopsis thaliana native to largely undisturbed habitats with two geographic and ecologically divergent subspecies. Here, we delineate the evolutionary processes driving the genetic diversity and population differentiation in this species. RESULTS: Using whole-genome re-sequencing data from 517 B. stricta accessions, we identify four genetic groups that diverged around 30-180 thousand years ago, with long-term small effective population sizes and recent population expansion after the Last Glacial Maximum. We find three genomic regions with elevated nucleotide diversity, totaling about 10% of the genome. These three regions of elevated nucleotide diversity show excess of intermediate-frequency alleles, higher absolute divergence (dXY), and lower relative divergence (FST) than genomic background, and significant enrichment in immune-related genes, reflecting long-term balancing selection. Scattered across the genome, we also find regions with both high FST and dXY among the groups, termed FST-islands. Population genetic signatures indicate that FST-islands with elevated divergence, which have experienced directional selection, are derived from divergent sorting of ancient polymorphisms. CONCLUSIONS: Our results suggest that long-term balancing selection on disease resistance genes may have maintained ancestral haplotypes across different geographical lineages, and unequal sorting of balanced polymorphisms may have generated genomic regions with elevated divergence. This study highlights the importance of ancestral balanced polymorphisms as crucial components of genome-wide variation.

On the postglacial spread of human commensal Arabidopsis thaliana: journey to the East.

With more sequenced genomes, our understanding of the demographic history of Arabidopsis thaliana is rapidly expanding. However, no-one has yet compiled previous data to investigate patterns of genetic variation across Eurasia. While sub-Saharan accessions have been reported to be the most divergent group, in the nuclear genome we found accessions from Yunnan, China to be genetically closest to the sub-Saharan group. In chloroplast, several deeply diverged haplogroups exist only in Eurasia, and African populations have lower variation in many haplogroups that they share with the Eurasian populations. These patterns cannot be easily explained by a single out-of-Africa event suggested previously. For more recent demographic history, we dated the nonrelict expansion to 10 ka. In the Chinese Yangtze nonrelicts, we found clear traces of gene flow with local relicts, and genes under strong selection were enriched for traces of relict introgression, especially those related to biotic and immune responses. The results suggest the ability of nonrelicts to obtain locally adaptive alleles through admixture with relicts is important for the expansion across environmental gradients of Eurasia. Our re-analyses provide another model for the early history as well as elucidating factors contributing to the recent demographic turnover event of this species.

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.

Ecological causes and consequences of flower color polymorphism in a self-pollinating plant (Boechera stricta).

Intraspecific variation in flower color is often attributed to pollinator-mediated selection, yet this mechanism cannot explain flower color polymorphisms in self-pollinating species. Indirect selection mediated via biotic and abiotic stresses could maintain flower color variation in these systems. The selfing forb, Boechera stricta, typically displays white flowers, but some individuals produce purple flowers. We quantified environmental correlates of flower color in natural populations. To disentangle plasticity from genotypic variation, we performed a multiyear field experiment in five gardens. In controlled conditions, we evaluated herbivore preferences and the effects of drought stress and soil pH on flower color expression. In natural populations, purple-flowered individuals experienced lower foliar herbivory than did their white-flowered counterparts. This pattern also held in the common gardens. Additionally, low-elevation environments induced pigmented flowers (plasticity), and the likelihood of floral pigmentation decreased with source elevation of maternal families (genetic cline). Viability selection favored families with pigmented flowers. In the laboratory, herbivores exerted greater damage on tissue derived from white- vs purple-flowered individuals. Furthermore, drought induced pigmentation in white-flowered lineages, and white-flowered plants had a fecundity advantage in the well-watered control. Flower color variation in selfing species is probably maintained by herbivory, drought stress, and other abiotic factors that vary spatially.