Population genetic structure of two Medicago species shaped by distinct life form, mating system and seed dispersal

Back to the future: Using herbarium specimens to isolate nodule-associated bacteria

Herbarium specimens are increasingly being used as sources of information to understand the ecology and evolution of plants and their associated microbes. Most studies have used specimens as a source of genetic material using culture-independent approaches. We demonstrate that herbarium specimens can also be used to culture nodule-associated bacteria, opening the possibility of using specimens to understand plant-microbe interactions at new spatiotemporal scales. We used historic and contemporary nodules of a common legume, Medicago lupulina, to create a culture collection. We were able to recover historic bacteria in 15 genera from three specimens (collected in 1950, 2004, and 2015). This work is the first of its kind to isolate historic bacteria from herbarium specimens. Future work should include inoculating plants with historic strains to see if they produce nodules and if they affect plant phenotype and fitness. Although we were unable to recover any Ensifer, the main symbiont of Medicago lupulina, we recovered some other potential nodulating species, as well as many putative growth-promoting bacteria.

Evolutionary genomics of climatic adaptation and resilience to climate change in alfalfa

Given the escalating impact of climate change on agriculture and food security, gaining insights into the evolutionary dynamics of climatic adaptation and uncovering climate-adapted variation can empower the breeding of climate-resilient crops to face future climate change. Alfalfa (Medicago sativa subsp. sativa), the queen of forages, shows remarkable adaptability across diverse global environments, making it an excellent model for investigating species responses to climate change. In this study, we performed population genomic analyses using genome resequencing data from 702 accessions of 24 Medicago species to unravel alfalfa's climatic adaptation and genetic susceptibility to future climate change. We found that interspecific genetic exchange has contributed to the gene pool of alfalfa, particularly enriching defense and stress-response genes. Intersubspecific introgression between M. sativa subsp. falcata (subsp. falcata) and alfalfa not only aids alfalfa's climatic adaptation but also introduces genetic burden. A total of 1671 genes were associated with climatic adaptation, and 5.7% of them were introgressions from subsp. falcata. By integrating climate-associated variants and climate data, we identified populations that are vulnerable to future climate change, particularly in higher latitudes of the Northern Hemisphere. These findings serve as a clarion call for targeted conservation initiatives and breeding efforts. We also identified pre-adaptive populations that demonstrate heightened resilience to climate fluctuations, illuminating a pathway for future breeding strategies. Collectively, this study enhances our understanding about the local adaptation mechanisms of alfalfa and facilitates the breeding of climate-resilient alfalfa cultivars, contributing to effective agricultural strategies for facing future climate change.

Serpentine environment prevails over geographic distribution in shaping the genetic diversity of Medicago lupulina L

Edaphic conditions of serpentine soils, naturally rich in heavy metals, act as a strong selection pressure that shapes specific metal-tolerant ecotypes. Medicago lupulina L. (black medick) is not only a widespread plant species that prefers calcareous and dry soil types but also grows at the borders of serpentine formations. It can also be found in waste and disturbed habitats. This is a species with reported phytoremediation potential, however, there is no published data regarding the impact of the environment on the genetic distribution of this species. The aim of our research was to explore how selection pressure of serpentine soils affects genetic diversity of M. lupulina and to test heavy-metal accumulation capacity of this species. Specimens of 11 M. lupulina populations were collected from serpentine outcrops located in Central and Eastern Bosnia as well as from non-serpentine sites. Soil and plant samples were analyzed for the total contents of heavy metals using air-acetylene flame atomic absorption spectroscopy. Genetic diversity was analyzed using AFLP (Amplified Fragment Length Polymorphism) markers. Serpentine soils showed high nickel, cobalt, chromium and iron concentrations. Nickel and manganese concentrations in soil samples and plant material showed statistically significant correlation. Although plants in two populations show the ability to extract Ni, M. lupulina does not show hyperaccumulating properties. Despite severe selective pressure, genetic diversity in serpentine populations is not reduced. Analyses of intrapopulation and interpopulation genetic diversity showed significant genetic differentiation among populations which is not related to their geographic distance. Population from non-metalliferous soil showed clear separation from all other populations. Diversity data suggest that serpentine populations maintain genetic diversity by undetected mechanisms and that edaphic factors rather than geography influence genetic structure analyzed M. lupulina populations.

Genetic Diversity and Differentiation of Eleven Medicago Species from Campania Region Revealed by Nuclear and Chloroplast Microsatellites Markers

The species belonging to the genus Medicago are considered a very important genetic resource at global level both for planet’s food security and for sustainable rangelands management. The checklist of the Italian flora (2021) includes a total number of 40 Medicago species for Italy, and 27 for Campania region, with a number of doubtful records or related to species no more found in the wild. In this study, 10 Medicago species native to Campania region, and one archaeophyte (M. sativa), identified by means of morphological diagnostic characters, were analyzed in a blind test to assay the efficacy of nine microsatellite markers (five cp-SSRs and four n-SSRs). A total number of 33 individuals from 6 locations were sampled and genotyped. All markers were polymorphic, 40 alleles were obtained with n-SSRs ranging from 8–12 alleles per locus with an average of 10 alleles per marker, PIC values ranged from 0.672 to 0.847, and the most polymorphic SSR was MTIC 564. The cp-SSRs markers were highly polymorphic too; PIC values ranged from 0.644 to 0.891 with an average of 0.776, the most polymorphic cp-SSR was CCMP10. 56 alleles were obtained with cp-SSRs ranging from 7 to 17 alleles per locus with an average of 11. AMOVA analysis with n-SSR markers highlighted a great level of genetic differentiation among the 11 species, with a statistically significant fixation index (FST). UPGMA clustering and Bayesian-based population structure analysis assigned these 11 species to two main clusters, but the distribution of species within clusters was not the same for the two analyses. In conclusion, our results demonstrated that the combination of the used SSRs well distinguished the 11 Medicago species. Moreover, our results demonstrated that the use of a limited number of SSRs might be considered for further genetic studies on other Medicago species.

Self-Fertilization, Inbreeding, and Yield in Alfalfa Seed Production

Selfing (self-pollination) is the ultimate form of inbreeding, or mating among close relatives. Selfing can create yield loss when inbreeding depression, defined as a lower survival and reproduction of inbred relative to outbred progeny, is present. To determine the impact of selfing in alfalfa (Medicago sativa L.), we quantified the selfing rate of 32 alfalfa seed production fields located in three regions, namely, the Pacific Northwest (PNW), the Central Valley of California (CEV), and the Imperial Valley of California (IMP). Selfing rates (the proportion of selfed seeds) varied between 5.3 and 30% with an average of 12.2% over the 32 seed production fields. In both the parents and their progeny, we observed an excess of heterozygotes relative to Hardy-Weinberg expectations. We detected notable levels of inbreeding in parents (0.231 ± 0.007 parental inbreeding coefficient) and progeny (0.229 ± 0.005). There were a 15% decrease in the number of seeds per stem (seed set) and a 13% decline in the number of seeds per pod in selfed relative to outcrossed stems, but negligible inbreeding depression for pods per raceme and seed weight. The number of racemes on selfed stems increased significantly in fields with greater selfing rates, supporting the presence of geitonogamous or among flower selfing. Despite the significant level of inbreeding depression, seed set did not decrease in fields with higher selfing rates, where the greater number of racemes on the selfed stems increased the seed set. The effects of the field selfing rate on the seed yield metrics were mostly indirect with direct effects of the number of racemes per stem. Available data indicate that the majority of selfing in alfalfa is pollinator-mediated, and thus, eliminating selfing in alfalfa seed production would require the selection of self-incompatible varieties, which, by eliminating inbreeding depression, would provide a 15% potential increase in seed yield and an increase in future hay yield.

Genetic variation in an ephemeral mudflat species: The role of the soil seed bank and dispersal in river and secondary anthropogenic habitats

Many ephemeral mudflat species, which rely on a soil seed bank to build up the next generation, are endangered in their natural habitat due to the widespread regulation of rivers. The aim of the present study was to elucidate the role of the soil seed bank and dispersal for the maintenance of genetic diversity in populations of near-natural river habitats and anthropogenic habitats created by traditional fish farming practices using Cyperus fuscus as a model. Using microsatellite markers, we found no difference in genetic diversity levels between soil seed bank and above-ground population and only moderate differentiation between the two fractions. One possible interpretation is the difference in short-term selection during germination under specific conditions (glasshouse versus field) resulting in an ecological filtering of genotypes out of the reservoir in the soil. River populations harbored significantly more genetic diversity than populations from the anthropogenic pond types. We suggest that altered levels and patterns of dispersal together with stronger selection pressures and historical bottlenecks in anthropogenic habitats are responsible for the observed reduction in genetic diversity. Dispersal is also supposed to largely prohibit genetic structure across Europe, although there is a gradient in private allelic richness from southern Europe (high values) to northern, especially north-western, Europe (low values), which probably relates to postglacial expansion out of southern and/or eastern refugia.

Environmental variation impacts trait expression and selection in the legume-rhizobium symbiosis

PREMISE

The ecological outcomes of mutualism are well known to shift across abiotic or biotic environments, but few studies have addressed how different environments impact evolutionary responses, including the intensity of selection on and the expression of genetic variance in key mutualism-related traits.

METHODS

We planted 30 maternal lines of the legume Medicago lupulina in four field common gardens and compared our measures of selection on and genetic variance in nodulation, a key trait reflecting legume investment in the symbiosis, with those from a previous greenhouse experiment using the same 30 M. lupulina lines.

RESULTS

We found that both the mean and genetic variance for nodulation were much greater in the greenhouse than in the field and that the form of selection on nodulation significantly differed across environments. We also found significant genotype-by-environment (G × E) effects for fitness-related traits that were generated by differences in the rank order of plant lines among environments.

CONCLUSIONS

Overall, our results suggest that the expression of genotypic variation and selection on nodulation differ across environments. In the field, significant rank-order changes for plant fitness potentially help maintain genetic variation in natural populations, even in the face of directional or stabilizing selection.

Rapid evolution in native plants cultivated for ecological restoration: not a general pattern

The growing number of restoration projects worldwide increases the demand for seed material of native species. To meet this demand, seeds are often produced through large-scale cultivation on specialised farms, using wild-collected seeds as the original sources. However, during cultivation, plants experience novel environmental conditions compared to those in natural populations, and there is a danger that the plants in cultivation are subject to unintended selection and lose their adaptation to natural habitats. Although the propagation methods are usually designed to maintain as much natural genetic diversity as possible, the effectiveness of these measures have never been tested. We obtained seed of five common grassland species from one of the largest native seed producers in Germany. For each species, the seeds were from multiple generations of seed production. We used AFLP markers and a common garden experiment to test for genetic and phenotypic changes during cultivation of these plants. The molecular markers detected significant evolutionary changes in three out of the five species and we found significant phenotypic changes in two species. The only species that showed substantial genetic and phenotypic changes was the short-lived and predominantly selfing Medicago lupulina, while in the other, mostly perennial and outcrossing species, the observed changes were mostly minor. Agricultural propagation of native seed material for restoration can cause evolutionary changes, at least in some species. We recommend caution, particularly in selfing and short-lived species, where evolution may be more rapid and effects may thus be more severe.

No evidence for adaptation to local rhizobial mutualists in the legume Medicago lupulina

Local adaptation is a common but not ubiquitous feature of species interactions, and understanding the circumstances under which it evolves illuminates the factors that influence adaptive population divergence. Antagonistic species interactions dominate the local adaptation literature relative to mutualistic ones, preventing an overall assessment of adaptation within interspecific interactions. Here, we tested whether the legume Medicago lupulina is adapted to the locally abundant species of mutualistic nitrogen-fixing rhizobial bacteria that vary in frequency across its eastern North American range. We reciprocally inoculated northern and southern M. lupulina genotypes with the northern (Ensifer medicae) or southern bacterium (E. meliloti) in a greenhouse experiment. Despite producing different numbers of root nodules (the structures in which the plants house the bacteria), neither northern nor southern plants produced more seeds, flowered earlier, or were more likely to flower when inoculated with their local rhizobia. We then used a pre-existing dataset to perform a genome scan for loci that showed elevated differentiation between field-collected plants that hosted different bacteria. None of the loci we identified belonged to the well-characterized suite of legume-rhizobia symbiosis genes, suggesting that the rhizobia do not drive genetic divergence between M. lupulina populations. Our results demonstrate that symbiont local adaptation has not evolved in this mutualism despite large-scale geographic variation in the identity of the interacting species.

Population Genetic Structure of Glycyrrhiza inflata B. (Fabaceae) Is Shaped by Habitat Fragmentation, Water Resources and Biological Characteristics

BACKGROUND

Habitat fragmentation, water resources and biological characteristics are important factors that shape the genetic structure and geographical distribution of desert plants. Analysis of the relationships between these factors and population genetic variation should help to determine the evolutionary potential and conservation strategies for genetic resources for desert plant populations. As a traditional Chinese herb, Glycyrrhiza inflata B. (Fabaceae) is restricted to the fragmented desert habitat in China and has undergone a dramatic decline due to long-term over-excavation. Determining the genetic structure of the G. inflata population and identifying a core collection could help with the development of strategies to conserve this species.

RESULTS

We investigated the genetic variation of 25 G. inflata populations based on microsatellite markers. A high level of population genetic divergence (FST = 0.257), population bottlenecks, reduced gene flow and moderate genetic variation (HE = 0.383) were detected. The genetic distances between the populations significantly correlated with the geographical distances, and this suggests that habitat fragmentation has driven a special genetic structure of G. inflata in China through isolation by distance. STRUCTURE analysis showed that G. inflata populations were structured into three clusters and that the populations belonged to multiple water systems, which suggests that water resources were related to the genetic structure of G. inflata. In addition, the biological characteristics of the perennial species G. inflata, such as its long-lived seeds, asexual reproduction, and oasis ecology, may be related to its resistance to habitat fragmentation. A core collection of G. inflata, that included 57 accessions was further identified, which captured the main allelic diversity of G. inflata.

CONCLUSIONS

Recent habitat fragmentation has accelerated genetic divergence. The population genetic structure of G. inflata has been shaped by habitat fragmentation, water resources and biological characteristics. This genetic information and core collection will facilitate the conservation of wild germplasm and breeding of this Chinese medicinal plant.

Evaluating Agronomic Performance and Investigating Molecular Structure of Drought and Heat Tolerant Wild Alfalfa (Medicago sativa L.) Collection from the Southeastern Turkey

Drought is a major stress factor for agricultural production including alfalfa production. One way to counterbalance the yield losses is the introgression of drought tolerant germplasm into breeding programs. As an effort to exploit such germplasm, 16 individual plants were selected from the Southeastern Turkey from their natural habitat and clonally propagated in field trials with an ultimate goal to use the germplasm as parents for releasing a synthetic cultivar. Forage yield and forage quality traits were evaluated and molecular genetic diversity among genotypes were determined using inter simple sequence repeat markers. Genotypes showed a variation from growth habit to yield and quality traits indicating sufficient phenotypic variation for diverse breeding efforts (for grazing or harvesting) and long term selection schemes. A large amount of genetic variation was observed even with a limited number of marker and genotypes. However, no pattern of spatial genetic structure was observed for the scale of the study when genetic variation is linked to the geographic origin. We conclude that ex situ natural variation provides a wealth of germplasm that could be incorporated into breeding programs aiming to improve drought tolerance. We also suggest an extensive collection of seeds/plant tissue from unique plants with desirable traits rather than putting more efforts to create a spatial germplasm sampling efforts in narrow regions.

Standing genetic variation in host preference for mutualist microbial symbionts

Many models of mutualisms show that mutualisms are unstable if hosts lack mechanisms enabling preferential associations with mutualistic symbiotic partners over exploitative partners. Despite the theoretical importance of mutualism-stabilizing mechanisms, we have little empirical evidence to infer their evolutionary dynamics in response to exploitation by non-beneficial partners. Using a model mutualism-the interaction between legumes and nitrogen-fixing soil symbionts-we tested for quantitative genetic variation in plant responses to mutualistic and exploitative symbiotic rhizobia in controlled greenhouse conditions. We found significant broad-sense heritability in a legume host's preferential association with mutualistic over exploitative symbionts and selection to reduce frequency of associations with exploitative partners. We failed to detect evidence that selection will favour the loss of mutualism-stabilizing mechanisms in the absence of exploitation, as we found no evidence for a fitness cost to the host trait or indirect selection on genetically correlated traits. Our results show that genetic variation in the ability to preferentially reduce associations with an exploitative partner exists within mutualisms and is under selection, indicating that micro-evolutionary responses in mutualism-stabilizing traits in the face of rapidly evolving mutualistic and exploitative symbiotic bacteria can occur in natural host populations.

Reduced plant competition among kin can be explained by Jensen's inequality

Plants often compete with closely related individuals due to limited dispersal, leading to two commonly invoked predictions on competitive outcomes. Kin selection, from evolutionary theory, predicts that competition between relatives will likely be weaker. The niche partitioning hypothesis, from ecological theory, predicts that competition between close relatives will likely be stronger. We tested for evidence consistent with either of these predictions by growing an annual legume in kin and nonkin groups in the greenhouse. We grew plant groups in treatments of symbiotic nitrogen fixing bacteria differing in strain identity and composition to determine if differences in the microbial environment can facilitate or obscure plant competition patterns consistent with kin selection or niche partitioning. Nonkin groups had lower fitness than expected, based on fitness estimates of the same genotypes grown among kin. Higher fitness among kin groups was observed in mixtures of N-fixing bacteria strains compared to single inoculations of bacteria strains present in the soil, which increased fitness differences between kin and nonkin groups. Lower fitness in nonkin groups was likely caused by increased competitive asymmetry in nonkin groups due to genetic differences in plant size combined with saturating relationships with plant size and fitness- i.e. Jensen's inequality. Our study suggests that microbial soil symbionts alter competitive dynamics among kin and nonkin. Our study also suggests that kin groups can have higher fitness, as predicted by kin selection theory, through a commonly heritable trait (plant size), without requiring kin recognition mechanisms.

A change in SHATTERPROOF protein lies at the origin of a fruit morphological novelty and a new strategy for seed dispersal in medicago genus

Angiosperms are the most diverse and numerous group of plants, and it is generally accepted that this evolutionary success owes in part to the diversity found in fruits, key for protecting the developing seeds and ensuring seed dispersal. Although studies on the molecular basis of morphological innovations are few, they all illustrate the central role played by transcription factors acting as developmental regulators. Here, we show that a small change in the protein sequence of a MADS-box transcription factor correlates with the origin of a highly modified fruit morphology and the change in seed dispersal strategies that occurred in Medicago, a genus belonging to the large legume family. This protein sequence modification alters the functional properties of the protein, affecting the affinities for other protein partners involved in high-order complexes. Our work illustrates that variation in coding regions can generate evolutionary novelties not based on gene duplication/subfunctionalization but by interactions in complex networks, contributing also to the current debate on the relative importance of changes in regulatory or coding regions of master regulators in generating morphological novelties.

Cotyledon damage affects seed number through final plant size in the annual grassland species Medicago lupulina

BACKGROUND AND AIMS

The effects of cotyledon damage on seedling growth and survival are relatively well established, but little is known about the effects on aspects of plant fitness such as seed number and size. Here the direct and indirect mechanisms linking cotyledon damage and plant fitness in the annual species Medicago lupulina are examined.

METHODS

Growth and reproductive traits, including mature plant size, time to first flowering, flower number, seed number and individual seed mass were monitored in M. lupulina plants when zero, one or two cotyledons were removed at 7 d old. Structural equation modelling (SEM) was used to examine the mechanisms linking cotyledon damage to seed number and seed mass.

KEY RESULTS

Cotyledon damage reduced seed number but not individual seed mass. The primary mechanism was a reduction in plant biomass with cotyledon damage that in turn reduced seed number primarily through a reduction in flower numbers. Although cotyledon damage delayed flower initiation, it had little effect on seed number. Individual seed mass was not affected by cotyledon removal, but there was a trade-off between seed number and seed mass.

CONCLUSIONS

It is shown how a network of indirect mechanisms link damage to cotyledons and fitness in M. lupulina. Cotyledon damage had strong direct effects on both plant size and flowering phenology, but an analysis of the causal relationships among plant traits and fitness components showed that a reduction in plant size associated with cotyledon damage was an important mechanism influencing fitness.

Genetic connectivity and inter-population seed dispersal of Banksia hookeriana at the landscape scale

BACKGROUND AND AIMS

Landscape genetics combines approaches from population genetics and landscape ecology, increasing the scope for conceptual advances in biology. Banksia hookeriana comprises clusters of individuals located on dune crests (geographical populations) physically separated by uninhabitable swales, with local extinctions common through frequent fire and/or severe drought.

METHODS

A landscape genetics approach was used to explore landscape-scale genetic connectivity and structure among geographical populations of B. hookeriana on 18 physically separated dunes located within a heterogeneous landscape of 3 x 5 km. These geographical populations were separated by approx. 0.1 to >1 km of unsuitable intervening swale habitat. Using 11 highly variable microsatellite loci, we utilized a Bayesian approach to identify genetic discontinuities within and between these geographical populations. Population allocation tests were then used to detect inter-dune seed dispersal inferred from assignment of individuals to a source population other than that from which they were collected.

KEY RESULTS

For the modal number of genetically distinct clusters (n = 17 genetic populations), two coincided with the geographical (dune) populations, eight spanned two to four geographical populations, and the remaining seven were spread among various parts of the sampled dunes, so that most geographical populations were spatially defined mosaics of individuals (subpopulations) belonging to two or more genetic populations. We inferred 25 inter-dune immigrants among the 582 individuals assessed, with an average distance between sink and source dunes of 1.1 km, and a maximum of 3.3 km.

CONCLUSIONS

The results show that genetic structure in an apparently strongly spatially structured landscape is not solely dependent on landscape structure, and that many physically defined geographical populations were genetic mosaics. More strikingly, there were physically separated individuals and groups of individuals that were part of the same genetically defined populations. We attribute this mismatch between spatially and genetically defined population structure to the varying closeness of the dunes and the ability of seeds to disperse long distances.