Alicia Mastretta-Yanes

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference

Restriction site‐associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single‐nucleotide polymorphisms. As an empirical example, we use a double‐digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high‐altitude mountains in Mexico.

Gene duplication, population genomics and species-level differentiation within a tropical mountain shrub

Gene duplication leads to paralogy, which complicates the de novo assembly of genotyping-by-sequencing (GBS) data. The issue of paralogous genes is exacerbated in plants, because they are particularly prone to gene duplication events. Paralogs are normally filtered from GBS data before undertaking population genomics or phylogenetic analyses. However, gene duplication plays an important role in the functional diversification of genes and it can also lead to the formation of postzygotic barriers. Using populations and closely related species of a tropical mountain shrub, we examine 1) the genomic differentiation produced by putative orthologs, and 2) the distribution of recent gene duplication among lineages and geography. We find high differentiation among populations from isolated mountain peaks and species-level differentiation within what is morphologically described as a single species. The inferred distribution of paralogs among populations is congruent with taxonomy and shows that GBS could be used to examine recent gene duplication as a source of genomic differentiation of nonmodel species.

Genomic variation in recently collected maize landraces from Mexico

The present dataset comprises 36,931 SNPs genotyped in 46 maize landraces native to Mexico as well as the teosinte subspecies Zea maiz ssp. parviglumis and ssp. mexicana. These landraces were collected directly from farmers mostly between 2006 and 2010. We accompany these data with a short description of the variation within each landrace, as well as maps, principal component analyses and neighbor joining trees showing the distribution of the genetic diversity relative to landrace, geographical features and maize biogeography. High levels of genetic variation were detected for the maize landraces (HE = 0.234 to 0.318 (mean 0.311), while slightly lower levels were detected in Zea m. mexicana and Zea m. parviglumis (HE = 0.262 and 0.234, respectively). The distribution of genetic variation was better explained by environmental variables given by the interaction of altitude and latitude than by landrace identity. This dataset is a follow up product of the Global Native Maize Project, an initiative to update the data on Mexican maize landraces and their wild relatives, and to generate information that is necessary for implementing the Mexican Biosafety Law.

Domestication Genomics of the Open-Pollinated Scarlet Runner Bean (Phaseolus coccineus L.)

The runner bean is a legume species from Mesoamerica closely related to common bean (Phaseolus vulgaris). It is a perennial species, but it is usually cultivated in small-scale agriculture as an annual crop for its dry seeds and edible immature pods. Unlike the common bean, P. coccineus has received little attention from a genetic standpoint. In this work we aim to (1) provide information about the domestication history and domestication events of P. coccineus; (2) examine the distribution and level of genetic diversity in wild and cultivated Mexican populations of this species; and, (3) identify candidate loci to natural and artificial selection. For this, we generated genotyping by sequencing data (42,548 SNPs) from 242 individuals of P. coccineus and the domesticated forms of the closely related species P. vulgaris (20) and P. dumosus (35). Eight genetic clusters were detected, of which half corresponds to wild populations and the rest to domesticated plants. The cultivated populations conform a monophyletic clade, suggesting that only one domestication event occurred in Mexico, and that it took place around populations of the Trans-Mexican Volcanic Belt. No difference between wild and domesticated levels of genetic diversity was detected and effective population sizes are relatively high, supporting a weak genetic bottleneck during domestication. Most populations presented an excess of heterozygotes, probably due to inbreeding depression. One population of P. coccineus subsp. striatus had the greatest excess and seems to be genetically isolated despite being geographically close to other wild populations. Contrasting with previous studies, we did not find evidence of recent gene flow between wild and cultivated populations. Based on outlier detection methods, we identified 24 domestication-related SNPs, 13 related to cultivar diversification and eight under natural selection. Few of these SNPs fell within annotated loci, but the annotated domestication-related SNPs are highly expressed in flowers and pods. Our results contribute to the understanding of the domestication history of P. coccineus, and highlight how the genetic signatures of domestication can be substantially different between closely related species.

Finding a Needle in a Haystack: Distinguishing Mexican Maize Landraces Using a Small Number of SNPs

In Mexicos territory, the center of origin and domestication of maize (Zea mays), there is a large phenotypic diversity of this crop. This diversity has been classified into “landraces.” Previous studies have reported that genomic variation in Mexican maize is better explained by environmental factors, particularly those related with altitude, than by landrace. Still, landraces are extensively used by agronomists, who recognize them as stable and discriminatory categories for the classification of samples. In order to investigate the genomic foundation of maize landraces, we analyzed genomic data (35,909 SNPs from Illumina MaizeSNP50 BeadChip) obtained from 50 samples representing five maize landraces (Comiteco, Conejo, Tehua, Zapalote Grande, and Zapalote Chico), and searched for markers suitable for landrace assignment. Landrace clusters could not be identified taking all the genomic information, but they become manifest taking only a subset of SNPs with high FST among landraces. Discriminant analysis of principal components was conducted to classify samples using SNP data. Two classification analyses were done, first classifying samples by landrace and then by altitude category. Through this classification method, we identified 20 landrace-informative SNPs and 14 altitude-informative SNPs, with only 6 SNPs in common for both analyses. These results show that Mexican maize phenotypic diversity can be classified in landraces using a small number of genomic markers, given the fact that landrace genomic diversity is influenced by environmental factors as well as artificial selection due to bio-cultural practices.

Morphological and niche divergence of pinyon pines

Abstract The environmental variables that define a species ecological niche should be associated with the evolutionary patterns present in the adaptations that resulted from living in these conditions. Thus, when comparing across species, we can expect to find an association between phylogenetically independent phenotypic characters and ecological niche evolution. Few studies have evaluated how organismal phenotypes might mirror patterns of niche evolution if these phenotypes reflect adaptations. Doing so could contribute on the understanding of the origin and maintenance of phenotypic diversity observed in nature. Here, we show the pattern of niche evolution of the pinyon pine lineage (Pinus subsection Cembroides); then, we suggest morphological adaptations possibly related to niche divergence, and finally, we test for correlation between ecological niche and morphology. We demonstrate that niche divergence is the general pattern within the clade and that it is positively correlated with adaptation.

Evolutionary and food supply implications of ongoing maize domestication by Mexican campesinos

Maize evolution under domestication is a process that continues today. Case studies suggest that Mexican smallholder family farmers, known as campesinos, contribute importantly to this, but their significance has not been explicitly quantified and analysed as a whole. Here, we examine the evolutionary and food security implications of the scale and scope under which campesinos produce maize. We gathered official municipal-level data on maize production under rainfed conditions and identified campesino agriculture as occurring in municipalities with average yields of less than or equal to 3 t ha−1. Environmental conditions vary widely in those municipalities and are associated with a great diversity of maize races, representing 85.3% of native maize samples collected in the country. We estimate that in those municipalities, around 1.38 × 1011 genetically different individual plants are subjected to evolution under domestication each season. This implies that 5.24 × 108 mother plants contribute to the next generation with their standing genetic diversity and rare alleles. Such a large breeding population size also increases the total number of adaptive mutations that may appear and be selected for. We also estimate that campesino agriculture could potentially feed around 54.7 million people in Mexico. These analyses provide insights about the contributions of smallholder agriculture around the world.

Long-term in situ persistence of biodiversity in tropical sky islands revealed by landscape genomics

Tropical mountains are areas of high species richness and endemism. Two historical phenomena may have contributed to this: (i) fragmentation and isolation of habitats may have promoted the genetic differentiation of populations and increased the possibility of allopatric divergence and speciation and (ii) the mountain areas may have allowed long‐term population persistence during global climate fluctuations. These two phenomena have been studied using either species occurrence data or estimating species divergence times. However, only few studies have used intraspecific genetic data to analyse the mechanisms by which endemism may emerge at the microevolutionary scale. Here, we use landscape analysis of genomic SNP data sampled from two high‐elevation plant species from an archipelago of tropical sky islands (the Trans‐Mexican Volcanic Belt) to test for population genetic differentiation, synchronous demographic changes and habitat persistence. We show that genetic differentiation can be explained by the degree of glacial habitat connectivity among mountains and that mountains have facilitated the persistence of populations throughout glacial/interglacial cycles. Our results support the ongoing role of tropical mountains as cradles for biodiversity by uncovering cryptic differentiation and limits to gene flow.

An Initiative for the Study and Use of Genetic Diversity of Domesticated Plants and Their Wild Relatives

Domestication has been influenced by formal plant breeding since the onset of intensive agriculture and the Green Revolution. Despite providing food security for some regions, intensive agriculture has had substantial detrimental consequences for the environment and does not fulfill smallholder’s needs under most developing countries conditions. Therefore, it is necessary to look for alternative plant production techniques, effective for each environmental, socio-cultural, and economic conditions. This is particularly relevant for countries that are megadiverse and major centers of plant domestication and diversification. In this white paper, a Mexico-centered initiative is proposed, with two main objectives: (1) to study, understand, conserve, and sustainably use the genetic diversity of domesticated plants and their wild relatives, as well as the ongoing evolutionary processes that generate and maintain it; and (2) to strengthen food and forestry production in a socially fair and environmentally friendly way. To fulfill these objectives, the initiative focuses on the source of variability available for domestication (genetic diversity and functional genomics), the context in which domestication acts (breeding and production) and one of its main challenges (environmental change). Research on these components can be framed to target and connect both the theoretical understanding of the evolutionary processes, the practical aspects of conservation, and food and forestry production. The target, main challenges, problems to be faced and key research questions are presented for each component, followed by a roadmap for the consolidation of this proposal as a national initiative.

Climatic niche evolution is faster in sympatric than allopatric lineages of the butterfly genus Pyrgus

Understanding how speciation relates to ecological divergence has long fascinated biologists. It is assumed that ecological divergence is essential to sympatric speciation, as a mechanism to avoid competition and eventually lead to reproductive isolation, while divergence in allopatry is not necessarily associated with niche differentiation. The impact of the spatial context of divergence on the evolutionary rates of abiotic dimensions of the ecological niche has rarely been explored for an entire clade. Here, we compare the magnitude of climatic niche shifts between sympatric versus allopatric divergence of lineages in butterflies. By combining next-generation sequencing, parametric biogeography and ecological niche analyses applied to a genus-wide phylogeny of Palaearctic Pyrgus butterflies, we compare evolutionary rates along eight climatic dimensions across sister lineages that diverged in large-scale sympatry versus allopatry. In order to examine the possible effects of the spatial scale at which sympatry is defined, we considered three sets of biogeographic assignments, ranging from narrow to broad definition. Our findings suggest higher rates of niche evolution along all climatic dimensions for sister lineages that diverge in sympatry, when using a narrow delineation of biogeographic areas. This result contrasts with significantly lower rates of climatic niche evolution found in cases of allopatric speciation, despite the biogeographic regions defined here being characterized by significantly different climates. Higher rates in allopatry are retrieved when biogeographic areas are too widely defined—in such a case allopatric events may be recorded as sympatric. Our results reveal the macro-evolutionary significance of abiotic niche differentiation involved in speciation processes within biogeographic regions, and illustrate the importance of the spatial scale chosen to define areas when applying parametric biogeographic analyses.