Advances in genotyping microsatellite markers through sequencing and consequences of scoring methods for Ceratonia siliqua (Leguminosae)

Phenotypic and genetic characterization of an Avena sativa L. germplasm collection of diverse origin: implications for food-oat breeding in Chile

Oats are known for their nutritional value and also for their beneficial properties on human health, such as the reduction of cholesterol levels and risk of coronary heart disease; they are an important export product for Chile. During the last decade (2010-2022) over 90% of the oat cultivated area in Chile has been covered with Avena sativa L. cv. Supernova INIA. This lack of genetic diversity in a context of climate change could limit the long-term possibility of growing oats in Chile. The present study is a phenotypic and genetic analysis of 132 oat cultivars and pure lines of diverse origin that can be considered as potential breeding material. The germplasm was evaluated for 28 traits and analyzed with 14 SSR markers. The effects of genotypes on phenotype were significant over all traits (P ≤ 0.05). Most traits exhibited moderate to high broad-sense heritability with exceptions such as yield (H2 = 0.27) and hulls staining (H2 = 0.32). Significant undesirable correlations between traits were generally of small biological importance, which is auspicious for achieving breeding objectives. Some of the heritability data and correlations provided here have not been previously reported. The overall phenotypic diversity was high (H' = 0.68 ± 0.18). The germplasm was grouped into three phenotypic clusters, differing in their qualities for breeding. Twenty-six genotypes outperforming Supernova INIA were identified for breeding of conventional food-oats. The genetic diversity of the germplasm was moderate on average (He = 0.58 ± 0.03), varying between 0.32 (AM22) and 0.77 (AME178). Two genetic subpopulations supported by the Structure algorithm exhibited a genetic distance of 0.24, showing low divergence of the germplasm. The diversity and phenotypic values found in this collection of oat genotypes are promising with respect to obtaining genetic gain in the short term in breeding programs. However, the similar genetic diversity, higher phenotypic diversity, and better phenotypic performance of the germplasm created in Chile compared to foreign germplasm suggest that germplasm harboring new genetic diversity will be key to favor yield and quality in new oat cultivars in the long term.

Unmasking microsatellite deceptiveness and debunking hybridization with SNPs in four marine copepod species of Calanus

Interspecific hybridization events are on the rise in natural systems due to climate change disrupting species barriers. Across taxa, microsatellites have long been the molecular markers of choice to identify admixed individuals. However, with the advent of high-throughput sequencing easing the generation of genome-wide datasets, incorrect reports of hybridization resulting from microsatellite technical artefacts have been uncovered in a growing number of taxa. In the marine zooplankton genus Calanus (Copepoda), whose species are used as climate change indicators, microsatellite markers have suggested hybridization between C. finmarchicus and C. glacialis, while other nuclear markers (InDels) never detected any admixed individuals, leaving the scientific community divided. Here, for the first time, we investigated the potential for hybridization among C. finmarchicus, C. glacialis, C. helgolandicus and C. hyperboreus using two large and independent SNP datasets. These were derived firstly from a protocol of target-capture applied to 179 individuals collected from 17 sites across the North Atlantic and Arctic Oceans, including sympatric areas, and second from published RNA sequences. All SNP-based analyses were congruent in showing that Calanus species are distinct and do not appear to hybridize. We then thoroughly re-assessed the microsatellites showing hybrids, with the support of published transcriptomes, and identified technical issues plaguing eight out of 10 microsatellites, including size homoplasy, paralogy, potential for null alleles and even two primer pairs targeting the same locus. Our study illustrates how deceptive microsatellites can be when applied to the investigation of hybridization.

Joint analysis of microsatellites and flanking sequences enlightens complex demographic history of interspecific gene flow and vicariance in rear-edge oak populations

Inference of recent population divergence requires fast evolving markers and necessitates to differentiate shared genetic variation caused by ancestral polymorphism and gene flow. Theoretical research shows that the use of compound marker systems integrating linked polymorphisms with different mutational dynamics, such as a microsatellite and its flanking sequences, can improve estimation of population structure and inference of demographic history, especially in the case of complex population dynamics. However, empirical application in natural populations has so far been limited by lack of suitable methods for data collection. A solution comes from the development of sequence-based microsatellite genotyping which we used to study molecular variation at 36 sequenced nuclear microsatellites in seven Quercus canariensis and four Q. faginea rear-edge populations across Algeria. We aim to decipher their taxonomic relationship, past evolutionary history and recent demographic trajectory. First, we compare the estimation of population genetics parameters and simulation-based inference of demographic history from microsatellite sequence alone, flanking sequence alone or the combination of linked microsatellite and flanking sequence variation. Second, we apply random forest approximate Bayesian computation to identify which of these sequence types is most informative. Whereas analysing microsatellite variation alone indicates recent interspecific gene flow, additional information gained by integrating nucleotide variation in flanking sequences, by reducing homoplasy, suggests ancient interspecific gene flow followed by drift in isolation instead. The weight of each polymorphism in the inference also demonstrates the value of linked variations with contrasted mutation dynamic to improve estimation of both demographic and mutational parameters.

Genome-wide footprints in the carob tree (Ceratonia siliqua) unveil a new domestication pattern of a fruit tree in the Mediterranean

Intense research efforts over the last two decades have renewed our understanding of plant phylogeography and domestication in the Mediterranean basin. Here we aim to investigate the evolutionary history and the origin of domestication of the carob tree (Ceratonia siliqua), which has been cultivated for millennia for food and fodder. We used >1000 microsatellite genotypes to delimit seven carob evolutionary units (CEUs). We investigated genome‐wide diversity and evolutionary patterns of the CEUs with 3557 single nucleotide polymorphisms generated by restriction‐site associated DNA sequencing (RADseq). To address the complex wild vs. cultivated status of sampled trees, we classified 56 sampled populations across the Mediterranean basin as wild, seminatural or cultivated. Nuclear and cytoplasmic loci were identified from RADseq data and separated for analyses. Phylogenetic analyses of these genomic‐wide data allowed us to resolve west‐to‐east expansions from a single long‐term refugium probably located in the foothills of the High Atlas Mountains near the Atlantic coast. Our findings support multiple origins of domestication with a low impact on the genetic diversity at range‐wide level. The carob was mostly domesticated from locally selected wild genotypes and scattered long‐distance westward dispersals of domesticated varieties by humans, concomitant with major historical migrations by Romans, Greeks and Arabs. Ex situ efforts to preserve carob genetic resources should prioritize accessions from both western and eastern populations, with emphasis on the most differentiated CEUs situated in southwest Morocco, south Spain and eastern Mediterranean. Our study highlights the relevance of wild and seminatural habitats in the conservation of genetic resources for cultivated trees.

Microsatellites and petal morphology reveal new patterns of admixture in Orchis hybrid zones

PREMISE

The genetic structure of hybrid zones provides insight into the potential for gene flow to occur between plant taxa. Four closely related European orchid species (Orchis anthropophora, O. militaris, O. purpurea, and O. simia) hybridize when they co-occur. We aimed to characterize patterns of hybridization in O. militaris-O. purpurea, O. purpurea-O. simia, and O. anthropophora-O. simia hybrid zones using molecular and morphological data.

METHODS

We used 11 newly isolated nuclear microsatellites to genotype 695 individuals collected from seven hybrid zones and six allopatric parental populations in France. Geometric morphometric analysis was conducted using 15 labellum landmarks to capture the main aspects of petal shape.

RESULTS

Backcrossing was asymmetric toward O. militaris in multiple O. militaris-O. purpurea hybrid zones. Hybrids in O. purpurea-O. simia and O. anthropophora-O. simia hybrid zones were largely limited to F1 and F2 generations, but further admixture had occurred. These patterns were reflected in labellum geometric morphometric data, which correlated strongly with nuclear microsatellite data in all three species combinations.

CONCLUSIONS

The coexistence of parental and admixed individuals in these Orchis hybrid zones implies they are likely to be tension zones being maintained by a balance between gene flow into the hybrid zone and selection acting against admixed individuals. The pattern of admixture in the three species combinations suggests intrinsic selection acting on the hybrids is weaker in more closely related taxa.

Revealing environmentally driven population dynamics of an Arctic diatom using a novel microsatellite PoolSeq barcoding approach

Ecological stability under environmental change is determined by both interspecific and intraspecific processes. Particularly for planktonic microorganisms, it is challenging to follow intraspecific dynamics over space and time. We propose a new method, microsatellite PoolSeq barcoding (MPB), for tracing allele frequency changes in protist populations. We successfully applied this method to experimental community incubations and field samples of the diatom Thalassiosira hyalina from the Arctic, a rapidly changing ecosystem. Validation of the method found compelling accuracy in comparison with established genotyping approaches within different diversity contexts. In experimental and environmental samples, we show that MPB can detect meaningful patterns of population dynamics, resolving allelic stability and shifts within a key diatom species in response to experimental treatments as well as different bloom phases and years. Through our novel MPB approach, we produced a large dataset of populations at different time-points and locations with comparably little effort. Results like this can add insights into the roles of selection and plasticity in natural protist populations under stable experimental but also variable field conditions. Especially for organisms where genotype sampling remains challenging, MPB holds great potential to efficiently resolve eco-evolutionary dynamics and to assess the mechanisms and limits of resilience to environmental stressors.

Pod Morphology, Primary and Secondary Metabolite Profiles in Non-grafted and Grafted Carob Germplasm Are Configured by Agro-Environmental Zone, Genotype, and Growing Season

Carob is a predominantly rainfed tree crop of high nutritive value and a long history of adaptation to the edaphoclimatic stress conditions of the Mediterranean. However, declining attention to the carob tree in recent decades has aggravated genetic erosion. The extant in situ germplasm varies both in terms of pod morphology and composition, reflecting the genetic and physiological divide chiefly among grafted and non-grafted material, and possibly the impact of variable agro-environments. Accordingly, the present study aimed to establish a systematic categorization of the genetic and phenotypic diversity encountered across carob germplasm identified in situ throughout Cyprus, a historical center of production and genetic diversity for the species. Linking pod morphology, primary and secondary metabolite profiles with genotyped source material originating in different agro-environments and crop seasons would provide a framework for interpreting (a) the interaction of these factors in configuring carob pod physicochemical constitution, and (b) the relative stability of phenotypic traits against environmental and seasonal variation. Microsatellite analysis discriminated 36 genotypes out of the 124 trees located in nine traditional agro-environmental zones and revealed low genetic diversity within the grafted germplasm. Two landraces were identified: "Tillyria," which is widespread and predominant, and "Kountourka," which is mainly localized to the northeastern peninsula of Karpasia. Morphological traits, such as seeds-to-pod weight ratio, pod width and thickness were principally under genetic control. Contrarily, compositional traits, particularly total phenolic content-including condensed tannins, in vitro antioxidant capacity and to a lesser extent gallic acid, organic acids and minerals were under agro-environmental control. Agro-environmental zone also modulated principally fructose and glucose; sucrose was modulated equally by genotype and agro-environment, while total sugars were under genetic control. Statistically significant differences between seasons were detected for all traits except for the seeds-to-pod weight ratio, pod length and width. Hierarchical cluster analysis corroborates that Cyprus may be divided into two major agro-environmental zones modulating the compositional properties of the carob pulp. The present study provides a comprehensive insight into the extant carob genetic resources of Cyprus and advances our understanding of how genetic, agro-environmental and seasonal factors interact in shaping carob pod morphology and composition.

Fast sequence-based microsatellite genotyping development workflow

Application of high-throughput sequencing technologies to microsatellite genotyping (SSRseq) has been shown to remove many of the limitations of electrophoresis-based methods and to refine inference of population genetic diversity and structure. We present here a streamlined SSRseq development workflow that includes microsatellite development, multiplexed marker amplification and sequencing, and automated bioinformatics data analysis. We illustrate its application to five groups of species across phyla (fungi, plant, insect and fish) with different levels of genomic resource availability. We found that relying on previously developed microsatellite assay is not optimal and leads to a resulting low number of reliable locus being genotyped. In contrast, de novo ad hoc primer designs gives highly multiplexed microsatellite assays that can be sequenced to produce high quality genotypes for 20-40 loci. We highlight critical upfront development factors to consider for effective SSRseq setup in a wide range of situations. Sequence analysis accounting for all linked polymorphisms along the sequence quickly generates a powerful multi-allelic haplotype-based genotypic dataset, calling to new theoretical and analytical frameworks to extract more information from multi-nucleotide polymorphism marker systems.

Microsatellite marker development in the crop wild relative Linum bienne using genome skimming

PREMISE

Nuclear microsatellite markers were developed for Linum bienne, the sister species of the crop L. usitatissimum, to provide molecular genetic tools for the investigation of L. bienne genetic diversity and structure.

METHODS AND RESULTS

Fifty microsatellite loci were identified in L. bienne by means of genome skimming, and 44 loci successfully amplified. Of these, 16 loci evenly spread across the L. usitatissimum reference nuclear genome were used for genotyping six L. bienne populations. Excluding one monomorphic locus, the number of alleles per locus ranged from two to 12. Four out of six populations harbored private alleles. The levels of expected and observed heterozygosity were 0.076 to 0.667 and 0.000 to 1.000, respectively. All 16 loci successfully cross-amplified in L. usitatissimum.

CONCLUSIONS

The 16 microsatellite loci developed here can be used for population genetic studies in L. bienne, and 28 additional loci that successfully amplified are available for further testing.