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.

Fostering relational autonomy in end-of-life care: a procedural approach and three-dimensional decision-making model.

Respect for patient autonomy is paramount in resolving ethical tensions in end-of-life care. The concept of relational autonomy has contributed to this debate; however, scholars often use this concept in a fragmented manner. This leads to partial answers on ascertaining patients' true wishes, meaningfully engaging patients' significant others, balancing interests among patients and significant others, and determining clinicians' obligations to change patients' unconventional convictions to enhance patient autonomy. A satisfactory solution based on relational autonomy must incorporate patients' competence (apart from decisional capacity), authenticity (their true desires or beliefs) and the involvement level of their significant others. To that end, we argue that John Christman's procedural approach to relational autonomy provides critical insights, such as the diachronic or socio-historical personhood, sustained critical reflection and his recent explication of the nature of asymmetrical relationships and helpful interlocutors. This study reviews Christman's account, proposes minor modifications and advocates for an integrated three-dimensional model for medical decision-making. Clarifying the relationship among the three elements promotes an ethical framework with a coherent understanding of relational autonomy. This model not only provides a descriptive and normative framework for end-of-life care practice but also reconsiders the nature of the clinician-patient relationship and its normative implications. We further present a case study to illustrate the merits of our proposed model. Altogether, our proposal will help navigate complex medical decision-making, foster trust and negotiate shared values between patients and their significant others, particularly in end-of-life care.

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.

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.

Understanding family dynamics in adult-to-adult living donor liver transplantation decision-making in Taiwan: Motivation, communication, and ambivalence.

This study investigated how multiple family members co-construct the ethical significance of living donor liver transplantation (LDLT) and how the family structure and dynamics of donor-recipient-caregiver relationships shape the communication and decision-making process within the sociocultural context of Taiwan. We conducted in-depth interviews with 36 participants from 13 families at the Chang Gung Memorial Hospital from February to December 2019. Interviews were audio recorded, transcribed, and translated from Mandarin to English and analyzed utilizing grounded theory and thematic analysis. Our findings revealed that LDLT in Taiwan is not merely a personal choice of either donor or recipient but essentially a collaborative process of family-centered medical decision-making, intertwined with socioeconomic conditions, cultural and social norms, gender roles, and the division of labor in the household. The notion of reciprocity and indebtedness for family members is a central theme underlying the entire process of motivation, deliberation, and decision-making, thereby reinforcing the naturalness of LDLT. However, our study highlighted that donors from families with traditional gender hierarchy and power imbalance may experience psychological and social vulnerabilities. Conclusively, we suggest that healthcare professionals should be sensitive to the contextual and relational factors involved in family dynamics and provide appropriate support to the ambivalent donors.

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.

Variation in Parasitoid Virulence of Tetrastichus brontispae during the Targeting of Two Host Beetles.

In host-parasitoid interactions, antagonistic relationship drives parasitoids to vary in virulence in facing different hosts, which makes these systems excellent models for stress-induced evolutionary studies. Venom compositions varied between two strains of Tetrastichus brontispae, Tb-Bl and Tb-On. Tb-Bl targets Brontispa longissima pupae as hosts, and Tb-On is a sub-population of Tb-Bl, which has been experimentally adapted to a new host, Octodonta nipae. Aiming to examine variation in parasitoid virulence of the two strains toward two hosts, we used reciprocal injection experiments to compare effect of venom/ovarian fluids from the two strains on cytotoxicity, inhibition of immunity and fat body lysis of the two hosts. We found that Tb-Onvenom was more virulent towards plasmatocyte spreading, granulocyte function and phenoloxidase activity than Tb-Blvenom. Tb-Blovary was able to suppress encapsulation and phagocytosis in both hosts; however, Tb-Onovary inhibition targeted only B. longissima. Our data suggest that the venom undergoes rapid evolution when facing different hosts, and that the wasp has good evolutionary plasticity.

Calmodulin Binding to Connexin 35: Specializations to Function as an Electrical Synapse.

Calmodulin binding is a nearly universal property of gap junction proteins, imparting a calcium-dependent uncoupling behavior that can serve in an emergency to decouple a stressed cell from its neighbors. However, gap junctions that function as electrical synapses within networks of neurons routinely encounter large fluctuations in local cytoplasmic calcium concentration; frequent uncoupling would be impractical and counterproductive. We have studied the properties and functional consequences of calmodulin binding to the electrical synapse protein Connexin 35 (Cx35 or gjd2b), homologous to mammalian Connexin 36 (Cx36 or gjd2). We find that specializations in Cx35 calmodulin binding sites make it relatively impervious to moderately high levels of cytoplasmic calcium. Calmodulin binding to a site in the C-terminus causes uncoupling when calcium reaches low micromolar concentrations, a behavior prevented by mutations that eliminate calmodulin binding. However, milder stimuli promote calcium/calmodulin-dependent protein kinase II activity that potentiates coupling without interference from calmodulin binding. A second calmodulin binding site in the end of the Cx35 cytoplasmic loop, homologous to a calmodulin binding site present in many connexins, binds calmodulin with very low affinity and stoichiometry. Together, the calmodulin binding sites cause Cx35 to uncouple only at extreme levels of intracellular calcium.

Two-color fluorescent analysis of connexin 36 turnover: relationship to functional plasticity.

Gap junctions formed of connexin 36 (Cx36, also known as Gjd2) show tremendous functional plasticity on several time scales. Changes in connexin phosphorylation modify coupling in minutes through an order of magnitude, but recent studies also imply involvement of connexin turnover in regulating cell-cell communication. We utilized Cx36 with an internal HaloTag to study Cx36 turnover and trafficking in cultured cells. Irreversible, covalent pulse-chase labeling with fluorescent HaloTag ligands allowed clear discrimination of newly formed and pre-existing Cx36. Cx36 in junctional plaques turned over with a half-life of 3.1 h, and the turnover rate was unchanged by manipulations of protein kinase A (PKA) activity. In contrast, changes in PKA activity altered coupling within 20 min. New Cx36 in cargo vesicles was added directly to existing gap junctions and newly made Cx36 was not confined to points of addition, but diffused throughout existing gap junctions. Existing connexins also diffused into photobleached areas with a half-time of less than 2 s. In conclusion, studies of Cx36-HaloTag revealed novel features of connexin trafficking and demonstrated that phosphorylation-based changes in coupling occur on a different time scale than turnover.

Inactivation of the particulate methane monooxygenase (pMMO) in Methylococcus capsulatus (Bath) by acetylene.

Acetylene (HCCH) has a long history as a mechanism-based enzyme inhibitor and is considered an active-site probe of the particulate methane monooxygenase (pMMO). Here, we report how HCCH inactivates pMMO in Methylococcus capsulatus (Bath) by using high-resolution mass spectrometry and computational simulation. High-resolution MALDI-TOF MS of intact pMMO complexes has allowed us to confirm that the enzyme oxidizes HCCH to the ketene (C2H2O) intermediate, which then forms an acetylation adduct with the transmembrane PmoC subunit. LC-MS/MS analysis of the peptides derived from in-gel proteolytic digestion of the protein subunit identifies K196 of PmoC as the site of acetylation. No evidence is obtained for chemical modification of the PmoA or PmoB subunit. The inactivation of pMMO by a single adduct in the transmembrane PmoC domain is intriguing given the complexity of the structural fold of this large membrane-protein complex as well as the complicated roles played by the various metal cofactors in the enzyme catalysis. Computational studies suggest that the entry of hydrophobic substrates to, and migration of products from, the catalytic site of pMMO are controlled tightly within the transmembrane domain. Support of these conclusions is provided by parallel experiments with two related alkynes: propyne (CH3CCH) and trifluoropropyne (CF3CCH). Finally, we discuss the implication of these findings to the location of the catalytic site in pMMO.

[Study on release kinetics for gel of Aconitum brachypodum based on nonlinear mixed effect model].

OBJECTIVE: To establish the analytical method for the release kinetic (RK) of Aconitum Brachypodum gel based on the nonlinear mixed effect model (NLMEM), in order to rationally evaluate the drug release process and explain the release mechanism. METHOD: The zero-order kinetic model containing for non-corroded drug system with the random effect was taken as the base model. The fixed effect and random effect factors impacting the drug release were analyzed by PROC NLMIXED of SAS to establish the final typical model. Subsequently, 10 training subsets were randomly extracted from the primary data to respectively their RK models, calculate the corresponding predicted root-mean-square error and average relative error, and evaluate the model stability and prediction accuracy. RESULT: The burst effect F0 had a very significant effect on the RK model. Among the component factors, carbopol 940 showed an obvious effect on the inherence release speed constant k0 and the concentration gradient change constant a, with different variations on the basis of dosage range. The random effect factors of k0 and a had a significant impact. The final RK model was proved to be stable, effective and reliable in the cross validation. CONCLUSION: The drug release kinetic analysis method could be used to rationally evaluate the drug release process and explain the release mechanisms.

Evaluating shared decision making for dialysis initiation: A qualitative study on patient refusal of long-term dialysis in Taiwan.

RATIONALE: Previous studies have explored shared decision making (SDM) implementation to determine the renal replacement therapy modality; however, the SDM approach for dialysis initiation, especially when patients refuse physician suggestions for long-term dialysis, remains unclear. AIMS AND OBJECTIVES: This study aimed to explore physicians' responses towards patients' refusal of long-term dialysis during the SDM process and the thinking processes of both physicians and patients regarding dialysis refusal. METHOD: We conducted in-depth semi-structured interviews with 10 patients diagnosed with end-stage renal disease, each of whom refused long-term dialysis after physicians employed the SDM framework, and nine nephrologists at the Chang Gung Memorial Hospital, Taiwan, from March to May 2020. Interviews were audio-recorded, transcribed, and translated from Mandarin to English. They were then thematically analysed. RESULTS: Three main themes on dialysis initiation SDM implementation and the differences between physician and patient perceptions on patient treatment refusal were yielded. While the SDM approach for dialysis initiation developed by nephrologists in Taiwan respects patient decisions, physicians often actively persuade patients to undergo dialysis in case of treatment refusal. The motivation behind this approach is to promote the patient's best medical interests, particularly post-dialysis life quality, and to ensure a 'rational' medical decision is made. However, patients' perceptions of treatment refusal differ significantly from those of physicians, and their decision-making process is often iterative and based on comprehensive evaluation of immediate concerns beyond biomedical factors. CONCLUSIONS: Findings suggest that the current physician-led SDM approach for dialysis initiation characterises active persuasion with physicians' perspectives predominating the clinical encounter. To improve SDM implementation, we propose that physicians should acknowledge and understand patients' reasoning for dialysis refusal and the distinction between objective health and subjective well-being during the decision-making process.

Exploring the genetic makeup of Xanthomonas species causing bacterial spot in Taiwan: evidence of population shift and local adaptation.

The introduction of plant pathogens can quickly reshape disease dynamics in island agro-ecologies, representing a continuous challenge for local crop management strategies. Xanthomonas pathogens causing tomato bacterial spot were probably introduced in Taiwan several decades ago, creating a unique opportunity to study the genetic makeup and adaptive response of this alien population. We examined the phenotypic and genotypic identity of 669 pathogen entries collected across different regions of Taiwan in the last three decades. The analysis detected a major population shift, where X. euvesicatoria and X. vesicatoria races T1 and T2 were replaced by new races of X. perforans. After its introduction, race T4 quickly became dominant in all tomato-growing areas of the island. The genomic analysis of 317 global genomes indicates that the Xanthomonas population in Taiwan has a narrow genetic background, most likely resulting from a small number of colonization events. However, despite the apparent genetic uniformity, X. perforans race T4 shows multiple phenotypic responses in tomato lines. Additionally, an in-depth analysis of effector composition suggests diversification in response to local adaptation. These include unique mutations on avrXv3 which might allow the pathogen to overcome Xv3/Rx4 resistance gene. The findings underscore the dynamic evolution of a pathogen when introduced in a semi-isolated environment and provide insights into the potential management strategies for this important disease of tomato.

Landscape genomics reveals genetic signals of environmental adaptation of African wild eggplants.

Crop wild relatives (CWR) provide a valuable resource for improving crops. They possess desirable traits that confer resilience to various environmental stresses. To fully utilize crop wild relatives in breeding and conservation programs, it is important to understand the genetic basis of their adaptation. Landscape genomics associates environments with genomic variation and allows for examining the genetic basis of adaptation. Our study examined the differences in allele frequency of 15,416 single nucleotide polymorphisms (SNPs) generated through genotyping by sequencing approach among 153 accessions of 15 wild eggplant relatives and two cultivated species from Africa, the principal hotspot of these wild relatives. We also explored the correlation between these variations and the bioclimatic and soil conditions at their collection sites, providing a comprehensive understanding of the genetic signals of environmental adaptation in African wild eggplant. Redundancy analysis (RDA) results showed that the environmental variation explained 6% while the geographical distances among the collection sites explained 15% of the genomic variation in the eggplant wild relative populations when controlling for population structure. Our findings indicate that even though environmental factors are not the main driver of selection in eggplant wild relatives, it is influential in shaping the genomic variation over time. The selected environmental variables and candidate SNPs effectively revealed grouping patterns according to the environmental characteristics of sampling sites. Using four genotype-environment association methods, we detected 396 candidate SNPs (2.5% of the initial SNPs) associated with eight environmental factors. Some of these SNPs signal genes involved in pathways that help adapt to environmental stresses such as drought, heat, cold, salinity, pests, and diseases. These candidate SNPs will be useful for marker-assisted improvement and characterizing the germplasm of this crop for developing climate-resilient eggplant varieties. The study provides a model for applying landscape genomics to other crops' wild relatives.

Distinct neurogenic potential in the retinal margin and the pars plana of mammalian eye.

Unlike many other vertebrates, a healthy mammalian retina does not grow throughout life and lacks a ciliary margin zone capable of actively generating new neurons. The isolation of stem-like cells from the ciliary epithelium has led to speculation that the mammalian retina and/or surrounding tissues may retain neurogenic potential capable of responding to retinal damage. Using genetically altered mouse lines with varying degrees of retinal ganglion cell loss, we show that the retinal margin responds to ganglion cell loss by prolonging specific neurogenic activity, as characterized by increased numbers of Atoh7(LacZ)-expressing cells. The extent of neurogenic activity correlated with the degree of ganglion cell deficiency. In the pars plana, but not the retinal margin, cells remain proliferative into adulthood, marking the junction of pars plana and retinal margin as a niche capable of producing proliferative cells in the mammalian retina and a potential cellular source for retinal regeneration.

Scenario- and discussion-based approach for teaching preclinical medical students the socio-philosophical aspects of psychiatry.

BACKGROUND: This study used a scenario- and discussion-based approach to teach preclinical medical students the socio-philosophical aspects of psychiatry and qualitatively evaluated the learning outcomes in a medical humanities course in Taiwan. METHODS: The seminar session focused on three hypothetical psychiatry cases. Students discussed the cases in groups and were guided by facilitators from multiple disciplines and professions. At the end of the semester, students submitted a narrative report comprising their reflections on the cases and discussions. The authors utilized content analysis to categorize students' narratives into three facets, namely, the philosophical, social and individual. RESULTS: In total 163 preclinical medical students participated in the class; 150 of them mentioned the scenario-based lesson in their reports; 33.3% of these reports discussed the case at the philosophical dimension (n = 50), 45.3% at the social dimension (n = 68), and 26.6% at the individual dimension (n = 40). Four major themes emerged: (1) a psychiatric diagnosis has far-reaching consequences for an individual's life, (2) the social structure affects how patients experience psychiatric disorders, (3) students related personal experience or those of friends and family to understand psychiatric disorders, and (4) medical humanities are of particular importance in psychiatric education. CONCLUSIONS: This study demonstrated that the scenario-based discussions led by a multidisciplinary team of facilitators can benefit medical students with limited clinical experience to contemplate the socio-philosophical aspects of psychiatry. The authors suggest that this pedagogical model during preclinical education should be encouraged.

Sources of Calcium at Connexin 36 Gap Junctions in the Retina.

Synaptic plasticity is a fundamental feature of the CNS that controls the magnitude of signal transmission between communicating cells. Many electrical synapses exhibit substantial plasticity that modulates the degree of coupling within groups of neurons, alters the fidelity of signal transmission, or even reconfigures functional circuits. In several known examples, such plasticity depends on calcium and is associated with neuronal activity. Calcium-driven signaling is known to promote potentiation of electrical synapses in fish Mauthner cells, mammalian retinal AII amacrine cells, and inferior olive neurons, and to promote depression in thalamic reticular neurons. To measure local calcium dynamics in situ, we developed a transgenic mouse expressing a GCaMP calcium biosensor fused to Connexin 36 (Cx36) at electrical synapses. We examined the sources of calcium for activity-dependent plasticity in retina slices using confocal or Super-Resolution Radial Fluctuations imaging. More than half of Cx36-GCaMP gap junctions responded to puffs of glutamate with transient increases in fluorescence. The responses were strongly dependent on NMDA receptors, in keeping with known activity-dependent signaling in some amacrine cells. We also found that some responses depended on the activity of voltage-gated calcium channels, representing a previously unrecognized source of calcium to control retinal electrical synaptic plasticity. The high prevalence of calcium signals at electrical synapses in response to glutamate application indicates that a large fraction of electrical synapses has the potential to be regulated by neuronal activity. This provides a means to tune circuit connectivity dynamically based on local activity.

Neurobeachin controls the asymmetric subcellular distribution of electrical synapse proteins.

The subcellular positioning of synapses and their specialized molecular compositions form the fundamental basis of neural circuits. Like chemical synapses, electrical synapses are constructed from an assortment of adhesion, scaffolding, and regulatory molecules, yet little is known about how these molecules localize to specific neuronal compartments. Here, we investigate the relationship between the autism- and epilepsy-associated gene Neurobeachin, the neuronal gap junction channel-forming Connexins, and the electrical synapse scaffold ZO1. Using the zebrafish Mauthner circuit, we find Neurobeachin localizes to the electrical synapse independently of ZO1 and Connexins. By contrast, we show Neurobeachin is required postsynaptically for the robust localization of ZO1 and Connexins. We demonstrate that Neurobeachin binds ZO1 but not Connexins. Finally, we find Neurobeachin is required to restrict electrical postsynaptic proteins to dendrites, but not electrical presynaptic proteins to axons. Together, the results reveal an expanded understanding of electrical synapse molecular complexity and the hierarchical interactions required to build neuronal gap junctions. Further, these findings provide novel insight into the mechanisms by which neurons compartmentalize the localization of electrical synapse proteins and provide a cell biological mechanism for the subcellular specificity of electrical synapse formation and function.

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.