Multinational evaluation of genetic diversity indicators for the Kunming-Montreal Global Biodiversity Framework.

Under the recently adopted Kunming-Montreal Global Biodiversity Framework, 196 Parties committed to reporting the status of genetic diversity for all species. To facilitate reporting, three genetic diversity indicators were developed, two of which focus on processes contributing to genetic diversity conservation: maintaining genetically distinct populations and ensuring populations are large enough to maintain genetic diversity. The major advantage of these indicators is that they can be estimated with or without DNA-based data. However, demonstrating their feasibility requires addressing the methodological challenges of using data gathered from diverse sources, across diverse taxonomic groups, and for countries of varying socio-economic status and biodiversity levels. Here, we assess the genetic indicators for 919 taxa, representing 5271 populations across nine countries, including megadiverse countries and developing economies. Eighty-three percent of the taxa assessed had data available to calculate at least one indicator. Our results show that although the majority of species maintain most populations, 58% of species have populations too small to maintain genetic diversity. Moreover, genetic indicator values suggest that IUCN Red List status and other initiatives fail to assess genetic status, highlighting the critical importance of genetic indicators.

Evaluating pollution-related damage and restoration success in urban forests with participatory monitoring and digital tools.

Peri-urban forest monitoring requires indicators of vegetation damage. An example is the sacred fir (Abies religiosa) forests surrounding Mexico City, which have been heavily exposed to tropospheric ozone, a harmful pollutant, for over 4 decades. We developed a participatory monitoring system with which local community members and scientists generated data on ozone tree damage. Santa Rosa Xochiac rangers (13) used the digital tool KoboToolBox to record ozone damage to trees, tree height, tree ages, tree condition, tree position, and whether the tree had been planted. Thirty-five percent of the trees (n = 1765) had ozone damage. Younger trees had a lower percentage of foliage damaged by ozone than older trees (p < 0.0001), and asymptomatic trees tended to be younger (p < 0.0001). Symptomatic trees were taller than asymptomatic trees of the same age (R2 c  = 0.43, R2 m  = 0.27). Involving local communities facilitated forest monitoring and using digital technology improved data quality. This participatory system can be used to monitor forest condition change over time and thus aids restoration efforts driven by government or local communities' interests, facilitating local decision-making.

Evaluación del daño relacionado a la contaminación y del éxito de la restauración de los bosques urbanos con un monitoreo participativo y herramientas digitales Resumen El monitoreo de los bosques periurbanos requiere indicadores de daños en la vegetación. Un ejemplo son los bosques de abeto (Abies religiosa) que rodean la Ciudad de México, pues desde hace más de 4 décadas han estado expuestos al ozono troposférico, un contaminante nocivo. Desarrollamos un sistema de monitoreo participativo con el que miembros de la comunidad local y científicos generaron datos sobre los daños causados por el ozono en los árboles. Los guardabosques de Santa Rosa Xochiac (133) utilizaron la herramienta digital KoboToolBox para registrar los daños causados por el ozono en los árboles, su altura, edad, estado, posición y si eran árboles plantados. El 35% de los árboles (n = 1,765) presentó daños por ozono. Los árboles más jóvenes tenían un menor porcentaje de follaje dañado por el ozono que los árboles más viejos (p<0.0001), y los árboles asintomáticos tendían a ser más jóvenes (p<0.0001). Los árboles sintomáticos eran más altos que los asintomáticos de la misma edad (R2 c = 0.43, R2 m = 0.27). La participación de las comunidades locales facilitó el monitoreo forestal y el uso de tecnología digital mejoró la calidad de los datos. Este sistema participativo puede utilizarse para monitorear los cambios en el estado de los bosques a lo largo del tiempo y contribuir a los esfuerzos de restauración impulsados por el gobierno o las comunidades locales, facilitando la toma de decisiones a nivel local.

Contemporary evolution of maize landraces and their wild relatives influenced by gene flow with modern maize varieties.

Mexico is recognized as the center of origin and domestication of maize. Introduction of modern maize varieties (MVs) into Mexico raised concerns regarding the possible effects of gene flow from MVs into maize landraces (LRs) and their wild relatives (WRs), teosintes. However, after more than 60 y from the release of the first MVs, the impact of the sympatry with LRs and their WRs has not been explored with genetic data. In this work, we assessed changes in the genomes of 7 maize LRs and 2 WR subspecies from collections spanning over 70 y. We compared the genotypes obtained by genotyping by sequencing (GBS) for LRs and WRs before and after the adoption of MVs, and observed introgression from sympatric MVs into LRs and into the WR Zea mays ssp. mexicana sampled after the year 2000. We also found a decrease in the paired divergence index (FST ) between MV-LR and MV-WR over the same time frame. Moreover, we determined that LR genetic diversity increased after 2000, probably as a result of gene flow from MVs introduced in the 1990s. Our findings allowed us to identify ongoing changes in the domesticated and wild maize genetic pools, and concur with previous works that have evaluated short-term gene flow from MVs into LRs in other crops. Our approach represents a useful tool for tracking evolutionary change in wild and domesticated genetic resources, as well as for developing strategies for their conservation.

Genetic structure of Plasmodium vivax using the merozoite surface protein 1 icb5-6 fragment reveals new hybrid haplotypes in southern Mexico.

BACKGROUND: Plasmodium vivax is a protozoan parasite with an extensive worldwide distribution, being highly prevalent in Asia as well as in Mesoamerica and South America. In southern Mexico, P. vivax transmission has been endemic and recent studies suggest that these parasites have unique biological and genetic features. The msp1 gene has shown high rate of nucleotide substitutions, deletions, insertions, and its mosaic structure reveals frequent events of recombination, maybe between highly divergent parasite isolates. METHODS: The nucleotide sequence variation in the polymorphic icb5-6 fragment of the msp1 gene of Mexican and worldwide isolates was analysed. To understand how genotype diversity arises, disperses and persists in Mexico, the genetic structure and genealogical relationships of local isolates were examined. To identify new sequence hybrids and their evolutionary relationships with other P. vivax isolates circulating worldwide two haplotype networks were constructed questioning that two portions of the icb5-6 have different evolutionary history. RESULTS: Twelve new msp1 icb5-6 haplotypes of P. vivax from Mexico were identified. These nucleotide sequences show mosaic structure comprising three partially conserved and two variable subfragments and resulted into five different sequence types. The variable subfragment sV1 has undergone recombination events and resulted in hybrid sequences and the haplotype network allocated the Mexican haplotypes to three lineages, corresponding to the Sal I and Belem types, and other more divergent group. In contrast, the network from icb5-6 fragment but not sV1 revealed that the Mexican haplotypes belong to two separate lineages, none of which are closely related to Sal I or Belem sequences. CONCLUSIONS: These results suggest that the new hybrid haplotypes from southern Mexico were the result of at least three different recombination events. These rearrangements likely resulted from the recombination between haplotypes of highly divergent lineages that are frequently distributed in South America and Asia and diversified rapidly.

Genetic, morphological, geographical and ecological approaches reveal phylogenetic relationships in complex groups, an example of recently diverged pinyon pine species (Subsection Cembroides).

Elucidating phylogenetic relationships and species boundaries within complex taxonomic groups is challenging for intrinsic and extrinsic (i.e., technical) reasons. Mexican pinyon pines are a complex group whose phylogenetic relationships and species boundaries have been widely studied but poorly resolved, partly due to intrinsic ecological and evolutionary features such as low morphological and genetic differentiation caused by recent divergence, hybridization and introgression. Extrinsic factors such as limited sampling and difficulty in selecting informative molecular markers have also impeded progress. Some of the Mexican pinyon pines are of conservation concern but others may remain unprotected because the species boundaries have not been established. In this study we combined approaches to resolve the phylogenetic relationships in this complex group and to establish species boundaries in four recently diverged taxa: P. discolor, P. johannis, P. culminicola and P. cembroides. We performed phylogenetic analyses using the chloroplast markers matK and psbA-trnH as well as complete and partial chloroplast genomes of species of Subsection Cembroides. Additionally, we performed a phylogeographic analysis combining genetic data (18 chloroplast markers), morphological data and geographical data to define species boundaries in four recently diverged taxa. Ecological divergence was supported by differences in climate among localities for distinct genetic lineages. Whereas the phylogenetic analysis inferred with matK and psbA-trnH was unable to resolve the relationships in this complex group, we obtained a resolved phylogeny with the use of the chloroplast genomes. The resolved phylogeny was concordant with a haplotype network obtained using chloroplast markers. In species with potential for recent divergence, hybridization or introgression, nonhierarchical network-based approaches are probably more appropriate to protect against misclassification due to incomplete lineage sorting. The boundaries among genetic lineages were delimited by the inclusion of morphological, geographical and ecological data in the haplotype network. These multiple lines of evidence help to assign species boundaries in this complex group. P. johannis, P. discolor, P. culminicola and P. cembroides are different species based on their genetic, morphological and ecological niche differences. We suggest a reevaluation of the conservation status of these species considering the information generated in this study.

Residual Effects of Transgenic Cotton on the Intestinal Microbiota of Dysdercus concinnus.

The interaction among plants, insects, and microbes (PIM) is a determinant factor for the assembly and functioning of natural and anthropic ecosystems. In agroecosystems, the relationships among PIM are based on the interacting taxa, environmental conditions, and agricultural management, including genetically modified (GM) organisms. Although evidence for the unintended effects of GM plants on non-target insects is increasingly robust, our knowledge remains limited regarding their impact on gut microbes and their repercussions on the host's ecology, especially in the wild. In this study, we compared the gut microbial community of Dysdercus concinnus bugs collected on wild cotton (Gossypium hirsutum), with and without insecticidal transgenes (cry1ab/ac), in its center of origin and diversity. By sequencing the V4-V5 region of 16S rRNA, we show differences in the diversity, structure, and topology of D. concinnus gut microbial interactions between specimens foraging cotton plants with and without transgenes. Identifying unintended residual effects of genetic engineering in natural ecosystems will provide first-line knowledge for informed decision-making to manage genetic, ecological, and evolutionary resources. Thus, determining which organisms interact with GM plants and how is the first step toward conserving natural ecosystems with evidence of transgenic introgression.

Multiple domestication events explain the origin of Gossypium hirsutum landraces in Mexico.

Several Mesoamerican crops constitute wild-to-domesticated complexes generated by multiple initial domestication events, and continuous gene flow among crop populations and between these populations and their wild relatives. It has been suggested that the domestication of cotton (Gossypium hirsutum) started in the northwest of the Yucatán Peninsula, from where it spread to other regions inside and outside of Mexico. We tested this hypothesis by assembling chloroplast genomes of 23 wild, landraces, and breeding lines (transgene-introgressed and conventional). The phylogenetic analysis showed that the evolutionary history of cotton in Mexico involves multiple events of introgression and genetic divergence. From this, we conclude that Mexican landraces arose from multiple wild populations. Our results also revealed that their structural and functional chloroplast organizations had been preserved. However, genetic diversity decreases as a consequence of domestication, mainly in transgene-introgressed (TI) individuals (π = 0.00020, 0.00001, 0.00016, 0, and 0, of wild, TI-wild, landraces, TI-landraces, and breeding lines, respectively). We identified homologous regions that differentiate wild from domesticated plants and indicate a relationship among the samples. A decrease in genetic diversity associated with transgene introgression in cotton was identified for the first time, and our outcomes are therefore relevant to both biosecurity and agrobiodiversity conservation.

Genomic diversity and population structure of teosinte (Zea spp.) and its conservation implications.

The wild species of the genus Zea commonly named teosintes, comprise nine different taxa, distributed from northern Mexico to Costa Rica. Although this genus of plants has been extensively studied from a morphological, ecogeographical and genetic point of view, most contributions have been limited to the study of a few populations and taxa. To understand the great variability that exists between and within teosinte species, it is necessary to include the vast majority of known populations. In this context, the objective of this work was to evaluate the diversity and genomic structure of 276 teosinte populations. Molecular analyzes were performed with 3,604 plants and with data from 33,929 SNPs. The levels of genetic diversity by taxonomic group show a marked difference between species, races and sections, where the highest values of genomic diversity was found in ssp. parviglumis and ssp. mexicana. The lower values were obtained for the Luxuriantes section as well as ssp. huehuetenagensis of the section Zea. The results of structure show that there is a great genetic differentiation in all the taxonomic groups considered. For ssp. parviglumis and mexicana, which are the taxa with the largest number of populations, a marked genomic differentiation was found that is consistent with their geographic distribution patterns. These results showed a loss of diversity in several teosinte populations, making a strong case for further collection, and ex situ and in situ conservation. Also, this study highlights the importance of integrating genomic diversity and structure for the applications of conservation and management.

Land-Use Change and Management Intensification Is Associated with Shifts in Composition of Soil Microbial Communities and Their Functional Diversity in Coffee Agroecosystems.

Despite the central role of microorganisms in soil fertility, little understanding exists regarding the impact of management practices and soil microbial diversity on soil processes. Strong correlations among soil microbial composition, management practices, and microbially mediated processes have been previously shown. However, limited integration of the different parameters has hindered our understanding of agroecosystem functioning. Multivariate analyses of these systems allow simultaneous evaluation of the parameters and can lead to hypotheses on the microbial groups involved in specific nutrient transformations. In the present study, using a multivariate approach, we investigated the effect of microbial composition (16SrDNA sequencing) and soil properties in carbon mineralization (CMIN) (BIOLOG™, Hayward, CA, USA) across different management categories on coffee agroecosystems in Mexico. Results showed that (i) changes in soil physicochemical variables were related to management, not to region, (ii) microbial composition was associated with changes in management intensity, (iii) specific bacterial groups were associated with different management categories, and (iv) there was a broader utilization range of carbon sources in non-managed plots. The identification of specific bacterial groups, management practices, and soil parameters, and their correlation with the utilization range of carbon sources, presents the possibility to experimentally test hypotheses on the interplay of all these components and further our understanding of agroecosystem functioning and sustainable management.

Incorporating evolutionary and threat processes into crop wild relatives conservation.

Crop wild relatives (CWR) intra- and interspecific diversity is essential for crop breeding and food security. However, intraspecific genetic diversity, which is central given the idiosyncratic threats to species in landscapes, is usually not considered in planning frameworks. Here, we introduce an approach to develop proxies of genetic differentiation to identify conservation areas, applying systematic conservation planning tools that produce hierarchical prioritizations of the landscape. It accounts for: (i) evolutionary processes, including historical and environmental drivers of genetic diversity, and (ii) threat processes, considering taxa-specific tolerance to human-modified habitats, and their extinction risk status. Our analyses can be used as inputs for developing national action plans for the conservation and use of CWR. Our results also inform public policy to mitigate threat processes to CWR (like crops living modified organisms or agriculture subsidies), and could advise future research (e.g. for potential germplasm collecting). Although we focus on Mesoamerican CWR within Mexico, our methodology offers opportunities to effectively guide conservation and monitoring strategies to safeguard the evolutionary resilience of any taxa, including in regions of complex evolutionary histories and mosaic landscapes.

Ongoing ecological and evolutionary consequences by the presence of transgenes in a wild cotton population.

After 25 years of genetically modified cotton cultivation in Mexico, gene flow between transgenic individuals and their wild relatives represents an opportunity for analysing the impacts of the presence of novel genes in ecological and evolutionary processes in natural conditions. We show comprehensive empirical evidence on the physiological, metabolic, and ecological effects of transgene introgression in wild cotton, Gossypium hirsutum. We report that the expression of both the cry and cp4-epsps genes in wild cotton under natural conditions altered extrafloral nectar inducibility and thus, its association with different ant species: the dominance of the defensive species Camponotus planatus in Bt plants, the presence of cp4-epsps without defence role of Monomorium ebeninum ants, and of the invasive species Paratrechina longicornis in wild plants without transgenes. Moreover, we found an increase in herbivore damage to cp4-epsps plants. Our results reveal the influence of transgene expression on native ecological interactions. These findings can be useful in the design of risk assessment methodologies for genetically modified organisms and the in situ conservation of G. hirsutum metapopulations.

Evaluating the monophyly of Mammillaria series Supertextae (Cactaceae).

Mammillaria (Cactaceae) taxonomy has been historically problematic due to the morphological variability and sympatry of the species. This has led to several proposals for infrageneric classification, including subgeneric, section and series categories. Mammillaria ser. Supertextae is one of 15 series and is made up of a variable set of species that are mainly distributed in southern Mexico and Central America. However, the phylogenetic relationships within M. ser. Supertextae and its relationship to other Mammillaria taxa are far from fully understood. Here we attempt to elucidate these relationships using complete terminal sampling and newly obtained chloroplast marker sequences and comparing them to Mammillaria species sequences from GenBank. Our phylogenetic analyses showed that M. ser. Supertextae comprises a well-supported monophyletic group that diverged approximately 2.1 Mya and has M. ser. Polyacanthae as its sister group; however, relationships within M. ser. Supertextae remain unresolved. The topology obtained within M. ser. Supertextae must also be interpreted under the distribution shared by these taxa, but it is difficult to differentiate ancestral polymorphisms from possible introgression, given the short time elapsed and the markers used. Our results show that the infrageneric units of M. haageana and M. albilanata can be considered independent evolutionary units. We also suggest that the relationship between M. haageana and M. albilanata is convoluted because their distribution overlaps (mainly towards southern Mexico), with genetic differences that possibly indicate they represent more than two taxonomic entities. One possible explanation is that there could still be gene flow between these taxa, and we might be witnessing an ongoing speciation process.

Host genotype explains rhizospheric microbial community composition: the case of wild cotton metapopulations (Gossypium hirsutum L.) in Mexico.

The rhizosphere provides several benefits to the plant host being a strong determinant for its health, growth and productivity. Nonetheless, the factors behind the assembly of the microbial communities associated with the rhizosphere such as the role of plant genotypes are not completely understood. In this study, we tested the role that intraspecific genetic variation has in rhizospheric microbial community assemblages, using genetically distinct wild cotton populations as a model of study. We followed a common garden experiment including five wild cotton populations, controlling for plant genotypes, environmental conditions and soil microbial community inoculum, to test for microbial differences associated with genetic variation of the plant hosts. Microbial communities of the treatments were characterized by culture-independent 16S rRNA gene amplicon sequencing with Illumina MiSeq platform. We analyzed microbial community diversity (alpha and beta), and diversity structure of such communities, determined by co-occurrence networks. Results show that different plant genotypes select for different and specific microbial communities from a common inoculum. Although we found common amplicon sequence variants (ASVs) to all plant populations (235), we also found unique ASVs for different populations that could be related to potential functional role of such ASVs in the rhizosphere.

The chloroplast genome of the pincushion cactus Mammilllaria haageana subsp. san-a ngelensis, a Mexican endangered species.

The genus Mammillaria occupies diverse habitats and exhibits diverse growth patterns and a large range of morphologies. Most of the species of this genus are used as ornamental plants and are subject to mass habitat loss. Due to these factors, they are being submitted to selective pressure that might affect conservational efforts and management plans. We obtained the 133 gene chloroplast genome as part of the project of sequencing the complete genome of pincushion cactus, including 88 protein-coding genes, 8 rRNA genes, and 37 tRNA genes. The phylogenetic tree indicates the pincushion cactus is a sister species of M. supertexta and M. huitzilopochtli.

A great diversity of Amerindian mitochondrial DNA ancestry is present in the Mexican mestizo population.

There are limited data on mitochondrial DNA (mtDNA) variation in the Mexican mestizo population. To examine the genetic diversity and matrilineal ancestry, the full mtDNA hypervariable regions I and II were sequenced in 270 unrelated mestizos from different regions of Mexico. A total of 202 different haplotypes were identified and the haplotype diversity was 0.9945. Amerindian haplotypes predominated in the sample with a proportion of 93.3%, followed by European (6.0%) and African haplotypes (0.7%). The frequency of the Amerindian haplogroups A2, B2, C1 and D1 was 51.1, 17.8, 18.5 and 5.9%, respectively. The frequency of Amerindian haplogroups was higher in the central region than in Mexico City, whereas it was the contrary for European haplogroups. This difference was accounted principally by the high frequency of B2 haplotypes in the central region. The minimum spanning network, the mismatch distribution and Tajima's D neutrality test suggest a population expansion for each Amerindian haplogroup, which could be initiated more recently for haplogroups A2 and D1. The present knowledge combined with other nuclear genetic markers will be essential in future association studies to correct for genetic substructure in mestizo populations.

In vitro performance in cotton plants with different genetic backgrounds: the case of Gossypium hirsutum in Mexico, and its implications for germplasm conservation.

One of the best ex situ conservation strategies for wild germplasm is in vitro conservation of genetic banks. The success of in vitro conservation relies heavily on the micropropagation or performance of the species of interest. In the context of global change, crop production challenges and climate change, we face a reality of intensified crop production strategies, including genetic engineering, which can negatively impact biodiversity conservation. However, the possible consequences of transgene presence for the in vitro performance of populations and its implications for biodiversity conservation are poorly documented. In this study we analyzed experimental evidence of the potential effects of transgene presence on the in vitro performance of Gossypium hirsutum L. populations, representing the Mexican genetic diversity of the species, and reflect on the implications of such presence for ex situ genetic conservation of the natural variation of the species. We followed an experimental in vitro performance approach, in which we included individuals from different wild cotton populations as well as individuals from domesticated populations, in order to differentiate the effects of domestication traits dragged into the wild germplasm pool via gene flow from the effects of transgene presence. We evaluated the in vitro performance of five traits related to plant establishment (N = 300): propagation rate, leaf production rate, height increase rate, microbial growth and root development. Then we conducted statistical tests (PERMANOVA, Wilcoxon post-hoc tests, and NMDS multivariate analyses) to evaluate the differences in the in vitro performance of the studied populations. Although direct causality of the transgenes to observed phenotypes requires strict control of genotypes, the overall results suggest detrimental consequences for the in vitro culture performance of wild cotton populations in the presence of transgenes. This provides experimental, statistically sound evidence to support the implementation of transgene screening of plants to reduce time and economic costs in in vitro establishment, thus contributing to the overarching goal of germplasm conservation for future adaptation.

The Mating System of the Wild-to-Domesticated Complex of Gossypium hirsutum L. Is Mixed.

The domestication syndrome of many plants includes changes in their mating systems. The evolution of the latter is shaped by ecological and genetic factors that are particular to an area. Thus, the reproductive biology of wild relatives must be studied in their natural distribution to understand the mating system of a crop species as a whole. Gossypium hirsutum (upland cotton) includes both domesticated varieties and wild populations of the same species. Most studies on mating systems describe cultivated cotton as self-pollinated, while studies on pollen dispersal report outcrossing; however, the mating system of upland cotton has not been described as mixed and little is known about its wild relatives. In this study we selected two wild metapopulations for comparison with domesticated plants and one metapopulation with evidence of recent gene flow between wild relatives and the crop to evaluate the mating system of cotton's wild-to-domesticated complex. Using classic reproductive biology methods, our data demonstrate that upland cotton presents a mixed mating system throughout the complex. Given cotton's capacity for outcrossing, differences caused by the domestication process in cultivated individuals can have consequences for its wild relatives. This characterization of the diversity of the wild relatives in their natural distribution, as well as their interactions with the crop, will be useful to design and implement adequate strategies for conservation and biosecurity.

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

In Mexico's 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.

Domesticated, Genetically Engineered, and Wild Plant Relatives Exhibit Unintended Phenotypic Differences: A Comparative Meta-Analysis Profiling Rice, Canola, Maize, Sunflower, and Pumpkin.

Agronomic management of plants is a powerful evolutionary force acting on their populations. The management of cultivated plants is carried out by the traditional process of human selection or plant breeding and, more recently, by the technologies used in genetic engineering (GE). Even though crop modification through GE is aimed at specific traits, it is possible that other non-target traits can be affected by genetic modification due to the complex regulatory processes of plant metabolism and development. In this study, we conducted a meta-analysis profiling the phenotypic consequences of plant breeding and GE, and compared modified cultivars with wild relatives in five crops of global economic and cultural importance: rice, maize, canola, sunflower, and pumpkin. For these five species, we analyzed the literature with documentation of phenotypic traits that are potentially related to fitness for the same species in comparable conditions. The information was analyzed to evaluate whether the different processes of modification had influenced the phenotype in such a way as to cause statistical differences in the state of specific phenotypic traits or grouping of the organisms depending on their genetic origin [wild, domesticated with genetic engineering (domGE), and domesticated without genetic engineering (domNGE)]. In addition, we tested the hypothesis that, given that transgenic plants are a construct designed to impact, in many cases, a single trait of the plant (e.g., lepidopteran resistance), the phenotypic differences between domGE and domNGE would be either less (or inexistent) than between the wild and domesticated relatives (either domGE or domNGE). We conclude that (1) genetic modification (either by selective breeding or GE) can be traced phenotypically when comparing wild relatives with their domesticated relatives (domGE and domNGE) and (2) the existence and the magnitude of the phenotypic differences between domGE and domNGE of the same crop suggest consequences of genetic modification beyond the target trait(s).