Analysis of >3400 worldwide eggplant accessions reveals two independent domestication events and multiple migration-diversification routes.

Eggplant (Solanum melongena) is an important Solanaceous crop, widely cultivated and consumed in Asia, the Mediterranean basin, and Southeast Europe. Its domestication centers and migration and diversification routes are still a matter of debate. We report the largest georeferenced and genotyped collection to this date for eggplant and its wild relatives, consisting of 3499 accessions from seven worldwide genebanks, originating from 105 countries in five continents. The combination of genotypic and passport data points to the existence of at least two main centers of domestication, in Southeast Asia and the Indian subcontinent, with limited genetic exchange between them. The wild and weedy eggplant ancestor S. insanum shows admixture with domesticated S. melongena, similar to what was described for other fruit-bearing Solanaceous crops such as tomato and pepper and their wild ancestors. After domestication, migration and admixture of eggplant populations from different regions have been less conspicuous with respect to tomato and pepper, thus better preserving 'local' phenotypic characteristics. The data allowed the identification of misclassified and putatively duplicated accessions, facilitating genebank management. All the genetic, phenotypic, and passport data have been deposited in the Open Access G2P-SOL database, and constitute an invaluable resource for understanding the domestication, migration and diversification of this cosmopolitan vegetable.

Multi-environment association study highlights candidate genes for robust agronomic quantitative trait loci in a novel worldwide Capsicum core collection.

Investigating crop diversity through genome-wide association studies (GWAS) on core collections helps in deciphering the genetic determinants of complex quantitative traits. Using the G2P-SOL project world collection of 10 038 wild and cultivated Capsicum accessions from 10 major genebanks, we assembled a core collection of 423 accessions representing the known genetic diversity. Since complex traits are often highly dependent upon environmental variables and genotype-by-environment (G × E) interactions, multi-environment GWAS with a 10 195-marker genotypic matrix were conducted on a highly diverse subset of 350 Capsicum annuum accessions, extensively phenotyped in up to six independent trials from five climatically differing countries. Environment-specific and multi-environment quantitative trait loci (QTLs) were detected for 23 diverse agronomic traits. We identified 97 candidate genes potentially implicated in 53 of the most robust and high-confidence QTLs for fruit flavor, color, size, and shape traits, and for plant productivity, vigor, and earliness traits. Investigating the genetic architecture of agronomic traits in this way will assist the development of genetic markers and pave the way for marker-assisted selection. The G2P-SOL pepper core collection will be available upon request as a unique and universal resource for further exploitation in future gene discovery and marker-assisted breeding efforts by the pepper community.

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.

Global range expansion history of pepper (Capsicum spp.) revealed by over 10,000 genebank accessions.

Genebanks collect and preserve vast collections of plants and detailed passport information, with the aim of preserving genetic diversity for conservation and breeding. Genetic characterization of such collections has the potential to elucidate the genetic histories of important crops, use marker-trait associations to identify loci controlling traits of interest, search for loci undergoing selection, and contribute to genebank management by identifying taxonomic misassignments and duplicates. We conducted a genomic scan with genotyping by sequencing (GBS) derived single nucleotide polymorphisms (SNPs) of 10,038 pepper (Capsicum spp.) accessions from worldwide genebanks and investigated the recent history of this iconic staple. Genomic data detected up to 1,618 duplicate accessions within and between genebanks and showed that taxonomic ambiguity and misclassification often involve interspecific hybrids that are difficult to classify morphologically. We deeply interrogated the genetic diversity of the commonly consumed Capsicum annuum to investigate its history, finding that the kinds of peppers collected in broad regions across the globe overlap considerably. The method ReMIXTURE-using genetic data to quantify the similarity between the complement of peppers from a focal region and those from other regions-was developed to supplement traditional population genetic analyses. The results reflect a vision of pepper as a highly desirable and tradable cultural commodity, spreading rapidly throughout the globe along major maritime and terrestrial trade routes. Marker associations and possible selective sweeps affecting traits such as pungency were observed, and these traits were shown to be distributed nonuniformly across the globe, suggesting that human preferences exerted a primary influence over domesticated pepper genetic structure.

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.

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.

PCR-based assays for validation of single nucleotide polymorphism markers in rice and mungbean.

BACKGROUND: Single nucleotide polymorphism (SNP) markers are the method of choice for genetic analyses including diversity and quantitative trait loci (QTL) studies. Marker validation is essential for QTL studies, but the cost and workload are considerable when large numbers of markers need to be verified. Marker systems with low development costs would be most suitable for this task. RESULTS: We have tested allele specific polymerase chain reaction (PCR), tetra markers and a genotyping tool based on the single strand specific nuclease CEL-I to verify randomly selected SNP markers identified previously either with a SNP array or by genotyping by sequencing in rice and mungbean, respectively. The genotyping capacity of allele-specific PCR and tetra markers was affected by the sequence context surrounding the SNP; SNPs located in repeated sequences and in GC-rich stretches could not be correctly identified. In contrast, CEL-I digestion of mixed fragments produced from test and reference DNA reliably pinpointed the correct genotypes, yet scoring of the genotypes became complicated when multiple SNPs were present in the PCR fragments. A cost analysis showed that as long the sample number remains small, CEL-I genotyping is more cost-effective than tetra markers. CONCLUSIONS: CEL-I genotyping performed better in terms of genotyping accuracy and costs than tetra markers. The method is highly useful for validating SNPs in small to medium size germplasm panels.

Genomic and transcriptomic comparison of nucleotide variations for insights into bruchid resistance of mungbean (Vigna radiata [L.] R. Wilczek).

BACKGROUND: Mungbean (Vigna radiata [L.] R. Wilczek) is an important legume crop with high nutritional value in South and Southeast Asia. The crop plant is susceptible to a storage pest caused by bruchids (Callosobruchus spp.). Some wild and cultivated mungbean accessions show resistance to bruchids. Genomic and transcriptomic comparison of bruchid-resistant and -susceptible mungbean could reveal bruchid-resistant genes (Br) for this pest and give insights into the bruchid resistance of mungbean. RESULTS: Flow cytometry showed that the genome size varied by 61 Mb (mega base pairs) among the tested mungbean accessions. Next generation sequencing followed by de novo assembly of the genome of the bruchid-resistant recombinant inbred line 59 (RIL59) revealed more than 42,000 genes. Transcriptomic comparison of bruchid-resistant and -susceptible parental lines and their offspring identified 91 differentially expressed genes (DEGs) classified into 17 major and 74 minor bruchid-resistance-associated genes. We found 408 nucleotide variations (NVs) between bruchid-resistant and -susceptible lines in regions spanning 2 kb (kilo base pairs) of the promoters of 68 DEGs. Furthermore, 282 NVs were identified on exons of 148 sequence-changed-protein genes (SCPs). DEGs and SCPs comprised genes involved in resistant-related, transposable elements (TEs) and conserved metabolic pathways. A large number of these genes were mapped to a region on chromosome 5. Molecular markers designed for variants of putative bruchid-resistance-associated genes were highly diagnostic for the bruchid-resistant genotype. CONCLUSIONS: In addition to identifying bruchid-resistance-associated genes, we found that conserved metabolism and TEs may be modifier factors for bruchid resistance of mungbean. The genome sequence of a bruchid-resistant inbred line, candidate genes and sequence variations in promoter regions and exons putatively conditioning resistance as well as markers detecting these variants could be used for development of bruchid-resistant mungbean varieties.

Identification of single nucleotide polymorphism markers associated with resistance to bruchids (Callosobruchus spp.) in wild mungbean (Vigna radiata var. sublobata) and cultivated V. radiata through genotyping by sequencing and quantitative trait locus analysis.

BACKGROUND: Bruchid beetles are an important storage pest of grain legumes. Callosobruchus sp. infect mungbean (Vigna radiata) at low levels in the field, multiply during grain storage and can destroy seed stocks in a few months. Resistance against bruchid beetles has been found in wild mungbean V. radiata var. sublobata TC1966 and in cultivated mungbean line V2802. RESULTS: Bruchid resistance data were obtained from recombinant inbred line populations TC1966 (V. radiata var. sublobata) × NM92 (F12) and V2802 (V. radiata) × NM94 (F7). More than 6,000 single nucleotide polymorphic markers were generated through genotyping by sequencing (GBS) for each of these populations and were used to map bruchid resistance genes. One highly significant quantitative trait locus (QTL) associated with bruchid resistance was mapped to chromosome 5 on genetic maps of both populations, suggesting that TC1966 and V2802 contain the same resistance locus. Co-segregation of all markers associated with resistance indicated the presence of only one major resistance QTL on chromosome 5, while QTL analysis based on physical map positions of the markers suggested the presence of multiple QTLs on different chromosomes. The diagnostic capacity of the identified molecular markers located in the QTL to correctly predict resistance was up to 100 %. CONCLUSIONS: Molecular markers tightly linked to bruchid resistance loci of two different mungbean resistance sources were developed and validated. These markers are highly useful for developing resistant lines.

The AVRDC - The World Vegetable Center mungbean (Vigna radiata) core and mini core collections.

BACKGROUND: Large ex situ germplasm collections generally harbor a wide range of crop diversity. AVRDC--The World Vegetable Center is holding in trust the world's second largest mungbean (Vigna radiata) germplasm collection with more than 6,700 accessions. Screening large collections for traits of interest is laborious and expensive. To enhance the access of breeders to the diversity of the crop, mungbean core and mini core collections have been established. RESULTS: The core collection of 1,481 entries has been built by random selection of 20% of the accessions after geographical stratification and subsequent cluster analysis of eight phenotypic descriptors in the whole collection. Summary statistics, especially the low differences of means, equal variance of the traits in both the whole and core collection and the visual inspection of quantile-quantile plots comparing the variation of phenotypic traits present in both collections indicated that the core collection well represented the pattern of diversity of the whole collection. The core collection was genotyped with 20 simple sequence repeat markers and a mini core set of 289 accessions was selected, which depicted the allele and genotype diversity of the core collection. CONCLUSIONS: The mungbean core and mini core collections plus their phenotypic and genotypic data are available for distribution to breeders. It is expected that these collections will enhance the access to biodiverse mungbean germplasm for breeding.

Molecular markers associated with resistance to squash leaf curl China virus and tomato leaf curl New Delhi virus in tropical pumpkin (Cucurbita moschata Duchesne ex Poir.) breeding line AVPU1426.

Virus diseases are a major production constraint for pumpkin. Recessive resistance to squash leaf curl China virus and tomato leaf curl New Delhi virus has been mapped in Cucurbita moschata (Duchesne ex Poir.) breeding line AVPU1426 to chromosomes 7 and 8, respectively. Molecular markers tightly associated with the resistance loci have been developed and were able to correctly predict resistance and susceptibility with an accuracy of 99% for squash leaf curl China virus resistance and 94.34% for tomato leaf curl New Delhi virus in F2 and back cross populations derived from the original resistance source AVPU1426. The markers associated with resistance are recommended for use in marker-assisted breeding.

Current status, breeding strategies and future prospects for managing chilli leaf curl virus disease and associated begomoviruses in Chilli (Capsicum spp.).

Chilli leaf curl virus disease caused by begomoviruses, has emerged as a major threat to global chilli production, causing severe yield losses and economic harm. Begomoviruses are a highly successful and emerging group of plant viruses that are primarily transmitted by whiteflies belonging to the Bemisia tabaci complex. The most effective method for mitigating chilli leaf curl virus disease losses is breeding for host resistance to Begomovirus. This review highlights the current situation of chilli leaf curl virus disease and associated begomoviruses in chilli production, stressing the significant issues that breeders and growers confront. In addition, the various breeding methods used to generate begomovirus resistant chilli cultivars, and also the complicated connections between the host plant, vector and the virus are discussed. This review highlights the importance of resistance breeding, emphasising the importance of multidisciplinary approaches that combine the best of traditional breeding with cutting-edge genomic technologies. subsequently, the article highlights the challenges that must be overcome in order to effectively deploy begomovirus resistant chilli varieties across diverse agroecological zones and farming systems, as well as understanding the pathogen thus providing the opportunities for improving the sustainability and profitability of chilli production.

Critical Review of the Increasing Complexity of Access and Benefit-Sharing Policies of Genetic Resources for Genebank Curators and Plant Breeders-A Public and Private Sector Perspective.

Plant breeders develop competitive, high-yielding, resistant crop varieties that can cope with the challenges of biotic stresses and tolerate abiotic stresses, resulting in nutritious food for consumers worldwide. To achieve this, plant breeders need continuous and easy access to plant genetic resources (PGR) for trait screening, to generate new diversity that can be built into newly improved varieties. International agreements such as the Convention on Biological Diversity (CBD), the International Treaty on Plant Genetic Resources for Food and Agriculture (ITPGRFA) and the Nagoya Protocol recognised the sovereign rights of countries over their genetic resources. Under the CBD/Nagoya Protocol, countries are free to establish specific national legislations regulating germplasm access and benefit-sharing to be negotiated bilaterally. Consequently, access to PGR became increasingly restricted and cumbersome, resulting in a decrease in germplasm exchange. The ITPGRFA attempted to ease this situation by establishing a globally harmonised multilateral system (MLS). Unfortunately, the MLS is (still) restricted to a limited number of food and forage crops, with very few vegetable crops. Easy and continuous access to genetic diversity combined with equitable and fair sharing of derived benefits is a prerequisite to breeding new varieties. Facilitated access contributes to sustainable crop production and food and nutrition security; therefore, access to and, consequently, use of PGRFA needs to be improved. Thus, the authors recommend, among others, expanding the scope of the ITPGRFA to include all PGRFA and making them and all related information accessible under a Standard Material Transfer Agreement (SMTA) combined, if necessary, with a subscription system or a seed sales tax. Such a transparent, functional and efficient system would erase legal uncertainties and minimise transaction costs for conservers, curators and users of genetic resources, thus aiding plant breeders to fulfil their mission.

Plant Sample Collection and Shipment for Multi-omic Analyses and Phytosanitary Evaluation.

Plant sample preparation for analyses is a fundamental step in high-throughput omics strategies. Especially for plant metabolomics, quenching of hydrolytic enzymes able to affect metabolite concentrations is crucial for the accuracy of results. Given that DNA is usually less labile than metabolites, most sampling and shipment procedures able to preserve the metabolome are also suitable for preventing the degradation of plant DNA or of DNA of pathogens in the plant tissue. In this article, we describe all the steps of sample collection, shipment (including the phytosanitary issues of moving plant samples), and processing for combined genomics and metabolomics from a single sample, as well as the protocols used in our laboratories for downstream approaches for crop plants, allowing collection of multi-omic datasets in large experimental setups. The protocols have been adjusted to apply to both freeze-dried and fresh-frozen material to allow the processing of crop plant samples that will require long-distance transport. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Preparation of freeze-dried leaf disks for multiplexed PCR or DArT-Seq genotyping Basic Protocol 2: Medium-throughput preparation of pathogen-free nucleic acids for most genotyping-resequencing applications or pathogen detection Alternate Protocol: Low-throughput extraction of high-quality DNA for resequencing using commercial kits Support Protocol: DNA quality control Basic Protocol 3: Preparation of freeze-dried plant material for metabolomics Basic Protocol 4: Preparation of fresh-frozen plant material for metabolomics Basic Protocol 5: Preparation and shipment of metabolite extracts for metabolomic analyses Basic Protocol 6: Sample shipping and long-term storage.

Corrigendum: Molecular markers associated to two non-allelic genic male sterility genes in peppers (Capsicum annuum L.).

[This corrects the article DOI: 10.3389/fpls.2018.01343.].

Genetics, genomics, and breeding of black gram [Vigna mungo (L.) Hepper].

Black gram [Vigna mungo (L.) Hepper] is a highly nutritious grain legume crop, mainly grown in South and Southeast Asia, with the largest area in India, where the crop is challenged by several biotic and abiotic stresses leading to significant yield losses. Improving genetic gains to increase on-farm yields is the primary goal of black gram breeding programs. This could be achieved by developing varieties resistant to major diseases like mungbean yellow mosaic disease, urdbean leaf crinkle virus, Cercospora leaf spot, anthracnose, powdery mildew, and insect pests such as whitefly, cowpea aphids, thrips, stem flies, and bruchids. Along with increasing on-farm yields, incorporating market-preferred traits ensures the adoption of improved varieties. Black gram breeding programs rely upon a limited number of parental lines, leading to a narrow genetic base of the developed varieties. For accelerating genetic gain, there is an urgent need to include more diverse genetic material for improving traits for better adaptability and stress resistance in breeding populations. The present review summarizes the importance of black gram, the major biotic and abiotic stresses, available genetic and genomic resources, major traits for potential crop improvement, their inheritance, and the breeding approaches being used in black gram for the development of new varieties.

Genome wide association analysis for grain micronutrients and anti-nutritional traits in mungbean [Vigna radiata (L.) R. Wilczek] using SNP markers.

Mungbean is an important food grain legume for human nutrition and nutritional food due to its nutrient-dense seed, liked palatability, and high digestibility. However, anti-nutritional factors pose a significant risk to improving nutritional quality for bio-fortification. In the present study, genetic architecture of grain micronutrients (grain iron and zinc concentration) and anti-nutritional factors (grain phytic acid and tannin content) in association mapping panel of 145 diverse mungbean were evaluated. Based on all four parameters genotypes PUSA 1333 and IPM 02-19 were observed as desired genotypes as they had high grain iron and zinc concentration but low grain phytic acid and tannin content. The next generation sequencing (NGS)-based genotyping by sequencing (GBS) identified 14,447 genome-wide SNPs in a diverse selected panel of 127 mungbean genotypes. Population admixture analysis revealed the presence of four different ancestries among the genotypes and LD decay of ∼57.6 kb kb physical distance was noted in mungbean chromosomes. Association mapping analysis revealed that a total of 20 significant SNPs were shared by both GLM and Blink models associated with grain micronutrient and anti-nutritional factor traits, with Blink model identifying 35 putative SNPs. Further, this study identified the 185 putative candidate genes. Including potential candidate genes Vradi07g30190, Vradi01g09630, and Vradi09g05450 were found to be associated with grain iron concentration, Vradi10g04830 with grain zinc concentration, Vradi08g09870 and Vradi01g11110 with grain phytic acid content and Vradi04g11580 and Vradi06g15090 with grain tannin content. Moreover, two genes Vradi07g15310 and Vradi09g05480 showed significant variation in protein structure between native and mutated versions. The identified SNPs and candidate genes are potential powerful tools to provide the essential information for genetic studies and marker-assisted breeding program for nutritional improvement in mungbean.

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.