Phylogenomic and morphological data reveal hidden patterns of diversity in the national tree of Brazil, Paubrasilia echinata.

PREMISE: Paubrasilia echinata (common names, pau brasil, brazilwood) is the national tree of Brazil and an endangered species endemic to the Brazilian Atlantic Forest. Over its wide distribution of 2000 km, its leaflets morphology exhibits extensive plasticity. Three morphotypes are commonly identified based on leaf size, but it is unclear if they represent distinct taxa or a single polymorphic species. This study aims to clarify the taxonomic position of the three morphotypes to inform conservation decisions. METHODS: A morphometric study of leaf characters of herbarium specimens was coupled with genetic analyses using genotype-by-sequencing data. We used maximum-likelihood and coalescent methods to evaluate the phylogenetic and population structure of the species. We compared these with a morphological dendrogram built from hierarchical clustering. RESULTS: Two of the three morphotypes formed separately evolving lineages, the third morphotype formed two geographically separate lineages, and northern trees with intermediate leaf morphology formed a separate fifth lineage. Leaflet size varied by over 35-fold, and although morphological clustering generally matched the genetic patterns, there were some overlaps, highlighting the cryptic diversity within this group. CONCLUSIONS: Our genetic and morphological results provide some evidence that cultivated trees from different states in Brazil seem to have a limited genetic origin and do not reflect the broader genetic and geographical diversity of the species. As a result, more care is likely needed to preserve the overall genomic diversity of this endangered and iconic species.

The role of oceanic currents in the dispersal and connectivity of the mangrove Rhizophora mangle on the Southwest Atlantic region.

Dispersal is a crucial mechanism to living beings, allowing them to reach new resources such that populations and species can occupy new environments. However, directly observing the dispersal mechanisms of widespread species can be costly or even impractical, which is the case for mangrove trees. The influence of ocean currents on mangrove dispersal is increasingly evident; however, few studies mechanistically relate the patterns of population distribution with the dispersal by oceanic currents under an integrated framework. Here, we evaluate the role of oceanic currents on connectivity of Rhizophora mangle along the Southwest Atlantic. We inferred population genetic structure and migration rates, simulated the displacement of propagules and tested our hypotheses with Mantel tests and redundancy analysis. We observed populations structured in two major groups, north and south, which is corroborated by other studies with Rhizophora and other coastal plants. Inferred recent migration rates do not indicate ongoing gene flow between sites. Conversely, long-term migration rates were low across groups and contrasting dispersal patterns within each one, which is consistent with long-distance dispersal events. Our hypothesis tests suggest that both isolation by distance and isolation by oceanography (derived from the oceanic currents) can explain the neutral genetic variation of R. mangle in the region. Our findings expand current knowledge of mangrove connectivity and highlight how the association of molecular methods with oceanographic simulations improve the interpretation of the dispersal process. This integrative approach is a cost- and time-efficient strategy to include dispersal and connectivity data into marine protected areas planning and management.

Geographical Patterns of Genetic Variation in Locoto Chile (Capsicum pubescens) in the Americas Inferred by Genome-Wide Data Analysis.

The locoto chile (Capsicum pubescens) is a regionally important food crop grown and marketed throughout the mid-highlands of South andCentral America, but little is known about its evolution and the diversity it harbours. An initial scan of genetic diversity and structure across its cultivation range was conducted, the first one using a genomic approach. The RAD-sequencing methodology was applied to a sampling of C. pubescens germplasm consisting of 67 accessions from different American countries, covering its range of distribution/cultivation on the continent. The RAD-seq SNP data obtained clustered the accessions into three major groups, with a high degree of admixture/reticulation among them. Moderate but significant differentiation and geographic structuration were found, depicting a south-north pattern in the distribution of genetic variation. The highest levels of diversity were found among central-western Bolivian individuals, while the lowest was found across Central America-Mexican germplasm. This study provides new genome-wide supported insights into the diversity and differentiation of C. pubescens, as well as a starting point for more efficient use of its genetic variation and germplasm conservation efforts. The findings also contribute to understanding the evolutionary history of C. pubescens, but further investigation is needed to disentangle its origin and diversification under domestication.

Diversification of chiles (Capsicum, Solanaceae) through time and space: New insights from genome-wide RAD-seq data.

Capsicum L. encompasses 43 American species, including the five domesticated worldwide consumed sweet and hot chiles. This study presents new, updated and age-calibrated phylogenetic hypothesis on the genus evolution incorporating nearly all currently accepted Capsicum species. A new model of spatial diversification of Capsicum is proposed based on analyses of several thousands of genome-wide RAD-seq derived SNPs. Maximum likelihood approaches were used to reconstruct phylogenies and to estimate dates of all major splits. Ancestral ranges were inferred and diversification events were modeled in a time frame using a Bayesian approach. Nine clades corresponding to genetically and (mostly) geographically well-defined lineages, which diversified starting around mid-upper Miocene, were recovered with strong support. The Northern and Central Andes were inferred to represent the most likely ancestral range of the genus Capsicum. A few early vicariant and dispersal events were estimated to have driven the geographic divergence of the main Capsicum clades. Each lineage was inferred to have diversified within a distinct region of South America and expanded geographically to different extent. Extant species diversification was inferred to have begun at the beginning of Pliocene and continued through the Pleistocene. The Central Andes, represented mainly by the territory of present-day Bolivia, were proposed to play a central role in the diversification of lineages comprising domesticated capsicums and their wild allies. The genome-wide approach allowed for high resolution and support of deep phylogenetic nodes providing novel insights into the affinities of major lineages and clades as well as on the geographic expansion of Capsicum. This study provides the first dated evolutionary history of the genus encompassing most of the chile species diversity.

Genomic Traces of the Fruit Fly Anastrepha obliqua Associated with Its Polyphagous Nature.

Anastrepha obliqua (Macquart) (Diptera: Tephritidae) is an important pest in the neotropical region. It is considered a polyphagous insect, meaning it infests plants of different taxonomic families and readily colonizes new host plants. The change to new hosts can lead to diversification and the formation of host races. Previous studies investigating the effect of host plants on population structure and selection in Anastrepha obliqua have focused on the use of data from the mitochondrial DNA sequence and microsatellite markers of nuclear DNA, and there are no analyses at the genomic level. To better understand this issue, we used a pooled restriction site-associated DNA sequencing (pooled RAD-seq) approach to assess genomic differentiation and population structure across sympatric populations of Anastrepha obliqua that infest three host plants-Spondias purpurea (red mombin), Mangifera indica (mango) of the family Anacardiaceae and Averrhoa carambola (carambola) of the family Oxalidaceae-in sympatric populations of the species Anastrepha obliqua of Inter-Andean Valley of the Cauca River in southwestern Colombia. Our results show genomic differentiation of populations from carambola compared to mango and red mombin populations, but the genetic structure was mainly established by geography rather than by the host plant. On the other hand, we identified 54 SNPs in 23 sequences significantly associated with the use of the host plant. Of these 23 sequences, we identified 17 candidate genes and nine protein families, of which four protein families are involved in the nutrition of these flies. Future studies should investigate the adaptive processes undergone by phytophagous insects in the Neotropics, using fruit flies as a model and state-of-the-art molecular tools.

Investigating the origins and evolution of a glyphosate-resistant weed invasion in South America.

The global invasion, and subsequent spread and evolution of weeds provides unique opportunities to address fundamental questions in evolutionary and invasion ecology. Amaranthus palmeri is a widespread glyphosate-resistant (GR) weed in the USA. Since 2015, GR populations of A. palmeri have been confirmed in South America, raising questions about introduction pathways and the importance of pre- vs. post-invasion evolution of GR traits. We used RAD-sequencing genotyping to characterize genetic structure of populations from Brazil, Argentina, Uruguay and the USA. We also quantified gene copy number of the glyphosate target, 5-enolpyruvyl-3-shikimate phosphate synthase (EPSPS), and the presence of an extrachromosomal circular DNA (eccDNA) replicon known to confer glyphosate resistance in USA populations. Populations in Brazil, Argentina and Uruguay were only weakly differentiated (pairwise FST  ≤0.043) in comparison to USA populations (mean pairwise FST  =0.161, range =0.068-0.258), suggesting a single major invasion event. However, elevated EPSPS copy number and the EPSPS replicon were identified in all populations from Brazil and Uruguay, but only in a single Argentinean population. These observations are consistent with independent in situ evolution of glyphosate resistance in Argentina, followed by some limited recent migration of the eccDNA-based mechanism from Brazil to Argentina. Taken together, our results are consistent with an initial introduction of A. palmeri into South America sometime before the 1980s, and local evolution of GR in Argentina, followed by a secondary invasion of GR A. palmeri with the unique eccDNA-based mechanism from the USA into Brazil and Uruguay during the 2010s.

Pervasive admixture and the spread of a large-lipped form in a cichlid fish radiation.

Adaptive radiations have proven important for understanding the mechanisms and processes underlying biological diversity. The convergence of form and function, as well as admixture and adaptive introgression, are common in adaptive radiations. However, distinguishing between these two scenarios remains a challenge for evolutionary research. The Midas cichlid species complex (Amphilophus spp.) is a prime example of adaptive radiation, with phenotypic diversification occurring at various stages of genetic differentiation. One species, A. labiatus, has large fleshy lips, is associated with rocky lake substrates, and occurs patchily within Lakes Nicaragua and Managua. By contrast, the similar, but thin-lipped, congener, A. citrinellus, is more common and widespread. We investigated the evolutionary history of the large-lipped form, specifically regarding whether the trait has evolved independently in both lakes from ancestral thin-lipped populations, or via dispersal and/or admixture events. We collected samples from distinct locations in both lakes, and assessed differences in morphology and ecology. Using RAD-seq, we genotyped thousands of SNPs to measure population structure and divergence, demographic history, and admixture. We found significant between-species differences in ecology and morphology, local intraspecific differences in body shape and trophic traits, but only limited intraspecific variation in lip shape. Despite clear ecological differences, our genomic approach uncovered pervasive admixture between the species and low genomic differentiation, with species within lakes being genetically more similar than species between lakes. Taken together, our results suggest a single origin of large-lips, followed by pervasive admixture and adaptive introgression, with morphology being driven by local ecological opportunities, despite ongoing gene-flow.

Las radiaciones adaptativas han demostrado ser clave para entender los mecanismos y procesos responsables de la diversidad biológica. La convergencia en forma y función, así como la mezcla genética y la introgresión adaptativa, son algo común en las radiaciones adaptativas. Sin embargo, distinguir entre estos dos escenarios sigue siendo un desafío para la biología evolutiva. El complejo de especies del cíclido de Midas (Amphilophus spp.) es un ejemplo paradigmático de radiación adaptativa, con diversidad fenotípica en varias etapas de diferenciación genética. Una de las especies, A. labiatus, que tiene labios grandes y carnosos, se asocia a zonas rocosas de los lagos, y tiene una distribución irregular en los lagos Nicaragua y Managua. En cambio, A. citrinellus, es una especie similar pero con labios finos, más común y con una distribución más amplia. Investigamos la historia evolutiva de la especie de labios grandes y, en concreto, si este rasgo ha evolucionado de forma independiente en los dos grandes lagos nicaragüenses a partir de poblaciones ancestrales de labios finos, o si por el contrario se ha dispersado mediante migración y/o mezcla genética. Colectamos muestras de distintas localidades en ambos lagos y evaluamos las diferencias en morfología y ecología. Genotipamos miles de SNPs utilizando RAD-seq para medir la estructura genética, la divergencia, la historia demográfica y la mezcla genética de las poblaciones. Encontramos diferencias significativas entre especies en ecología y morfología, diferencias intraespecíficas locales en la forma del cuerpo y rasgos tróficos, pero sólo una limitada variación intraespecífica en la forma de los labios. A pesar de las claras diferencias ecológicas, el análisis genómico desveló una intensa mezcla genética entre especies, y una limitada diferenciación genómica, encontrando mayor semejanza genética entre especies dentro de un mismo lago, que entre especies de distintos lagos. Nuestros resultados sugieren un origen único de la especie de labios gruesos seguido de mezcla genética e introgresión adaptativa, e indican que la morfología habría sido modelada por las oportunidades ecológicas locales, a pesar del flujo génico.

Migration through a Major Andean Ecogeographic Disruption as a Driver of Genetic and Phenotypic Diversity in a Wild Tomato Species.

Evolutionary dynamics at the population level play a central role in creating the diversity of life on our planet. In this study, we sought to understand the origins of such population-level variation in mating systems and defensive acylsugar chemistry in Solanum habrochaites-a wild tomato species found in diverse Andean habitats in Ecuador and Peru. Using Restriction-site-Associated-DNA-Sequencing (RAD-seq) of 50 S. habrochaites accessions, we identified eight population clusters generated via isolation and hybridization dynamics of 4-6 ancestral populations. Detailed characterization of mating systems of these clusters revealed emergence of multiple self-compatible (SC) groups from progenitor self-incompatible populations in the northern part of the species range. Emergence of these SC groups was also associated with fixation of deleterious alleles inactivating acylsugar acetylation. The Amotape-Huancabamba Zone-a geographical landmark in the Andes with high endemism and isolated microhabitats-was identified as a major driver of differentiation in the northern species range, whereas large geographical distances contributed to population structure and evolution of a novel SC group in the central and southern parts of the range, where the species was also inferred to have originated. Findings presented here highlight the role of the diverse ecogeography of Peru and Ecuador in generating population differentiation, and enhance our understanding of the microevolutionary processes that create biological diversity.

Population Genomics in Rhamdia quelen (Heptapteridae, Siluriformes) Reveals Deep Divergence and Adaptation in the Neotropical Region.

Rhamdia quelen, a Neotropical fish with hybridization between highly divergent mitochondrial DNA (mtDNA) lineages, represents an interesting evolutionary model. Previous studies suggested that there might be demographic differences between coastal lagoons and riverine environments, as well as divergent populations that could be reproductively isolated. Here, we investigated the genetic diversity pattern of this taxon in the Southern Neotropical Basin system that includes the La Plata Basin, Patos-Merin lagoon basin and the coastal lagoons draining to the SW Atlantic Ocean, through a population genomics approach using 2b-RAD-sequencing-derived single nucleotide polymorphisms (SNPs). The genomic scan identified selection footprints associated with divergence and suggested local adaptation environmental drivers. Two major genomic clusters latitudinally distributed in the Northern and Southern basins were identified, along with consistent signatures of divergent selection between them. Population structure based on the whole set of loci and on the presumptive neutral vs. adaptive loci showed deep genomic divergence between the two major clusters. Annotation of the most consistent SNPs under divergent selection revealed some interesting candidate genes for further functional studies. Moreover, signals of adaptation to a coastal lagoon environment mediated by purifying selection were found. These new insights provide a better understanding of the complex evolutionary history of R. quelen in the southernmost basin of the Neotropical region.

Sympatric parallel diversification of major oak clades in the Americas and the origins of Mexican species diversity.

Oaks (Quercus, Fagaceae) are the dominant tree genus of North America in species number and biomass, and Mexico is a global center of oak diversity. Understanding the origins of oak diversity is key to understanding biodiversity of northern temperate forests. A phylogenetic study of biogeography, niche evolution and diversification patterns in Quercus was performed using 300 samples, 146 species. Next-generation sequencing data were generated using the restriction-site associated DNA (RAD-seq) method. A time-calibrated maximum likelihood phylogeny was inferred and analyzed with bioclimatic, soils, and leaf habit data to reconstruct the biogeographic and evolutionary history of the American oaks. Our highly resolved phylogeny demonstrates sympatric parallel diversification in climatic niche, leaf habit, and diversification rates. The two major American oak clades arose in what is now the boreal zone and radiated, in parallel, from eastern North America into Mexico and Central America. Oaks adapted rapidly to niche transitions. The Mexican oaks are particularly numerous, not because Mexico is a center of origin, but because of high rates of lineage diversification associated with high rates of evolution along moisture gradients and between the evergreen and deciduous leaf habits. Sympatric parallel diversification in the oaks has shaped the diversity of North American forests.

Incipient sympatric speciation in Midas cichlid fish from the youngest and one of the smallest crater lakes in Nicaragua due to differential use of the benthic and limnetic habitats?

Understanding how speciation can occur without geographic isolation remains a central objective in evolutionary biology. Generally, some form of disruptive selection and assortative mating are necessary for sympatric speciation to occur. Disruptive selection can arise from intraspecific competition for resources. If this competition leads to the differential use of habitats and variation in relevant traits is genetically determined, then assortative mating can be an automatic consequence (i.e., habitat isolation). In this study, we caught Midas cichlid fish from the limnetic (middle of the lake) and benthic (shore) habitats of Crater Lake Asososca Managua to test whether some of the necessary conditions for sympatric speciation due to intraspecific competition and habitat isolation are given. Lake As. Managua is very small (<900 m in diameter), extremely young (maximally 1245 years of age), and completely isolated. It is inhabited by, probably, only a single endemic species of Midas cichlids, Amphilophus tolteca. We found that fish from the limnetic habitat were more elongated than fish collected from the benthic habitat, as would be predicted from ecomorphological considerations. Stable isotope analyses confirmed that the former also exhibit a more limnetic lifestyle than the latter. Furthermore, split-brood design experiments in the laboratory suggest that phenotypic plasticity is unlikely to explain much of the observed differences in body elongation that we observed in the field. Yet, neutral markers (microsatellites) did not reveal any genetic clustering in the population. Interestingly, demographic inferences based on RAD-seq data suggest that the apparent lack of genetic differentiation at neutral markers could simply be due to a lack of time, as intraspecific competition may only have begun a few hundred generations ago.

Adaptive genetic variation distinguishes Chilean blue mussels (Mytilus chilensis) from different marine environments.

Chilean mussel populations have been thought to be panmictic with limited genetic structure. Genotyping-by-sequencing approaches have enabled investigation of genomewide variation that may better distinguish populations that have evolved in different environments. We investigated neutral and adaptive genetic variation in Mytilus from six locations in southern Chile with 1240 SNPs obtained with RAD-seq. Differentiation among locations with 891 neutral SNPs was low (FST = 0.005). Higher differentiation was obtained with a panel of 58 putative outlier SNPs (FST = 0.114) indicating the potential for local adaptation. This panel identified clusters of genetically related individuals and demonstrated that much of the differentiation (~92%) could be attributed to the three major regions and environments: extreme conditions in Patagonia, inner bay influenced by aquaculture (Reloncaví), and outer bay (Chiloé Island). Patagonia samples were most distinct, but additional analysis carried out excluding this collection also revealed adaptive divergence between inner and outer bay samples. The four locations within Reloncaví area were most similar with all panels of markers, likely due to similar environments, high gene flow by aquaculture practices, and low geographical distance. Our results and the SNP markers developed will be a powerful tool supporting management and programs of this harvested species.

Genomic tests of the species-pump hypothesis: Recent island connectivity cycles drive population divergence but not speciation in Caribbean crickets across the Virgin Islands.

Harnessing the power of genomic scans, we test the debated "species pump" hypothesis that implicates repeated cycles of island connectivity and isolation as drivers of divergence. This question has gone understudied given the limited resolution of past molecular markers for studying such dynamic phenomena. With an average of 32,000 SNPs from the genome of 136 individuals from 10 populations of a Caribbean flightless ground cricket species (Amphiacusta sanctaecrucis) and a complementary set of statistical approaches, we infer a stepping-stone colonization model and high levels of genetic differentiation across the Virgin Islands, which have been periodically interconnected until 8 ka. Estimates of divergence times from models based on the site frequency spectrum coincide with a period of repeated connection and fragmentation of the islands at 75-130 ka. These results are consistent with a role of island connectivity cycles in promoting genomic divergence and indicate that the genetic distinctiveness of island populations has persisted despite subsequent and extended interisland connections identified from bathymetric data. We discuss these findings in the broader context of Caribbean biogeography, and more specifically why high levels of genomic divergence across the Virgin Islands associated with repeated connectivity cycles do not actually translate into species diversification.