Forecasting effects of climate changes on the population genetic structure of Anacardium occidentale in the Cerrado biome, Brazil.

There has been a continuous interest in understanding the patterns of genetic diversity in natural populations because of the role of intraspecific genetic diversity as the basis of all evolutionary change and thus, its potential effects on population persistence when facing environmental changes. Here, we provided the first description of genetic diversity distribution and population structure of Anacardium occidentale L. (cashew) from the Brazilian Cerrado, one of the most economically important tropical crops in the world. We applied Bayesian clustering approaches (STRUCTURE and POPS) that allow predicting the effects of future climatic changes on the population genetic structure of A. occidentale. We identified distinct genetic groups corresponding to the southwestern, central, and northern regions of the species' range. The characterized genetic clusters will disappear under future climate change scenarios, leading to a homogenization of genetic variability across the landscape. Our findings suggest a high likelihood for the loss of genetic diversity, which in turn will reduce the evolutionary potential of the species to cope with predicted future climatic changes. Results from this study may help develop management strategies to conserve the genetic diversity and structure of natural cashew populations.

Complete Chloroplast Genomes of Pterodon emarginatus Vogel and Pterodon pubescens Benth: Comparative and Phylogenetic Analyses.

Background: The species Pterodon emarginatus and P. pubescens, popularly known as white sucupira or faveira, are native to the Cerrado biome and have the potential for medicinal use and reforestation. They are sister species with evolutionary proximity. Objective: Considering that the chloroplast genome exhibits a conserved structure and genes, the analysis of its sequences can contribute to the understanding of evolutionary, phylogenetic, and diversity issues. Methods: The chloroplast genomes of P. emarginatus and P. pubescens were sequenced on the Illumina MiSeq platform. The genomes were assembled based on the de novo strategy. We performed the annotation of the genes and the repetitive regions of the genomes. The nucleotide diversity and phylogenetic relationships were analyzed using the gene sequences of these species and others of the Leguminosae family, whose genomes are available in databases. Results: The complete chloroplast genome of P. emarginatus is 159,877 bp, and that of P. pubescens is 159,873 bp. The genomes of both species have circular and quadripartite structures. A total of 127 genes were predicted in both species, including 110 single-copy genes and 17 duplicated genes in the inverted regions. 141 microsatellite regions were identified in P. emarginatus and 140 in P. pubescens. The nucleotide diversity estimates of the gene regions in twenty-one species of the Leguminosae family were 0.062 in LSC, 0.086 in SSC, and 0.036 in IR. The phylogenetic analysis demonstrated the proximity between the genera Pterodon and Dipteryx, both from the clade Dipterygeae. Ten pairs of primers with potential for the development of molecular markers were designed. Conclusion: The genetic information obtained on the chloroplast genomes of P. emarginatus and P. pubescens presented here reinforces the similarity and evolutionary proximity between these species, with a similarity percentage of 99.8%.

A Pipeline for the Development of Microsatellite Markers using Next Generation Sequencing Data.

Background: Also known as Simple Sequence Repetitions (SSRs), microsatellites are profoundly informative molecular markers and powerful tools in genetics and ecology studies on plants. Objective: This research presents a workflow for developing microsatellite markers using genome skimming. Methods: The pipeline was proposed in several stages that must be performed sequentially: obtaining DNA sequences, identifying microsatellite regions, designing primers, and selecting candidate microsatellite regions to develop the markers. Our pipeline efficiency was analyzed using Illumina sequencing data from the non-model tree species Pterodon emarginatus Vog. Results: The pipeline revealed 4,382 microsatellite regions and drew 7,411 pairs of primers for P. emarginatus. However, a much larger number of microsatellite regions with the potential to develop markers were discovered from our pipeline. We selected 50 microsatellite regions with high potential for developing markers and organized 29 microsatellite regions in sets for multiplex PCR. Conclusion: The proposed pipeline is a powerful tool for fast and efficient development of microsatellite markers on a large scale in several species, especially nonmodel plant species.

Genetic variability and effective population size in Hymenaea stigonocarpa (Fabaceae) germplasm collection: tools for breeding programs and genetic conservation.

The conservation of plant genetic resources is essential for breeding programs. Regarding the native species of the Brazilian Cerrado biome, many studies have demonstrated their high potential for use in both medicines and foods. Hymenaea stigonocarpa, a tree with wide occurrence in the Cerrado, has economic importance, and due its extractive use, the establishment of a breeding program is relevant for sustainable use and conservation. Thus, the first germplasm collection of the species was installed at the Federal University of Goiás (UFG). To know the magnitude of genetic variability and how it was distributed in the collection, 353 individuals, distributed in 119 families from 24 subpopulations collected in the Cerrado biome, were genotyped using capillary electrophoresis. Nine pairs of microsatellite markers were genotyped. The UFG germplasm collection showed a high level of genetic diversity (mean [Formula: see text] = 0.554) at the evaluated loci. By Analysis of Molecular Variance (AMOVA), a significant genetic structure was detected (θP = 0.152, p < 0.01), which was expected since the subpopulations that originated the germplasm collection were collected in geographically distant locations. In addition, the germplasm collection had a population effective size of 54.9 and presented an allelic representation of 79.89% compared to 32 natural subpopulations. These results demonstrate that the germplasm collection preserves a high genetic diversity of H. stigonocarpa with a population effective size considered sufficient for the conduction of a breeding program.

Geographically weighted regression as a generalized Wombling to detect barriers to gene flow.

Barriers to gene flow play an important role in structuring populations, especially in human-modified landscapes, and several methods have been proposed to detect such barriers. However, most applications of these methods require a relative large number of individuals or populations distributed in space, connected by vertices from Delaunay or Gabriel networks. Here we show, using both simulated and empirical data, a new application of geographically weighted regression (GWR) to detect such barriers, modeling the genetic variation as a "local" linear function of geographic coordinates (latitude and longitude). In the GWR, standard regression statistics, such as R(2) and slopes, are estimated for each sampling unit and thus are mapped. Peaks in these local statistics are then expected close to the barriers if genetic discontinuities exist, capturing a higher rate of population differentiation among neighboring populations. Isolation-by-Distance simulations on a longitudinally warped lattice revealed that higher local slopes from GWR coincide with the barrier detected with Monmonier algorithm. Even with a relatively small effect of the barrier, the power of local GWR in detecting the east-west barriers was higher than 95 %. We also analyzed empirical data of genetic differentiation among tree populations of Dipteryx alata and Eugenia dysenterica Brazilian Cerrado. GWR was applied to the principal coordinate of the pairwise FST matrix based on microsatellite loci. In both simulated and empirical data, the GWR results were consistent with discontinuities detected by Monmonier algorithm, as well as with previous explanations for the spatial patterns of genetic differentiation for the two species. Our analyses reveal how this new application of GWR can viewed as a generalized Wombling in a continuous space and be a useful approach to detect barriers and discontinuities to gene flow.

Mantel test in population genetics.

The comparison of genetic divergence or genetic distances, estimated by pairwise FST and related statistics, with geographical distances by Mantel test is one of the most popular approaches to evaluate spatial processes driving population structure. There have been, however, recent criticisms and discussions on the statistical performance of the Mantel test. Simultaneously, alternative frameworks for data analyses are being proposed. Here, we review the Mantel test and its variations, including Mantel correlograms and partial correlations and regressions. For illustrative purposes, we studied spatial genetic divergence among 25 populations of Dipteryx alata ("Baru"), a tree species endemic to the Cerrado, the Brazilian savannas, based on 8 microsatellite loci. We also applied alternative methods to analyze spatial patterns in this dataset, especially a multivariate generalization of Spatial Eigenfunction Analysis based on redundancy analysis. The different approaches resulted in similar estimates of the magnitude of spatial structure in the genetic data. Furthermore, the results were expected based on previous knowledge of the ecological and evolutionary processes underlying genetic variation in this species. Our review shows that a careful application and interpretation of Mantel tests, especially Mantel correlograms, can overcome some potential statistical problems and provide a simple and useful tool for multivariate analysis of spatial patterns of genetic divergence.

Chloroplast genome characterization of Uncaria guianensis and Uncaria tomentosa and evolutive dynamics of the Cinchonoideae subfamily.

Uncaria species are used in traditional medicine and are considered of high therapeutic value and economic importance. This work describes the assembly and annotation of the chloroplast genomes of U. guianensis and U. tomentosa, as well as a comparative analysis. The genomes were sequenced on MiSeq Illumina, assembled with NovoPlasty, and annotated using CHLOROBOX GeSeq. Addictionaly, comparative analysis were performed with six species from NCBI databases and primers were designed in Primer3 for hypervariable regions based on the consensus sequence of 16 species of the Rubiaceae family and validated on an in-silico PCR in OpenPrimeR. The genome size of U. guianensis and U. tomentosa was 155,505 bp and 156,390 bp, respectively. Both Species have 131 genes and GC content of 37.50%. The regions rpl32-ccsA, ycf1, and ndhF-ccsA showed the three highest values of nucleotide diversity within the species of the Rubiaceae family and within the Uncaria genus, these regions were trnH-psbA, psbM-trnY, and rps16-psbK. Our results indicates that the primer of the region ndhA had an amplification success for all species tested and can be promising for usage in the Rubiaceae family. The phylogenetic analysis recovered a congruent topology to APG IV. The gene content and the chloroplast genome structure of the analyzed species are conserved and most of the genes are under negative selection. We provide the cpDNA of Neotropical Uncaria species, an important genomic resource for evolutionary studies of the group.

A geographical genetics framework for inferring homing reproductive behavior in fishes.

One of the most intriguing patterns of migration and gene flow that affects genetic structure is the reproductive homing behavior of fishes, wherein the adults return to the areas in which they were spawned. Here we reviewed the literature on homing behavior in fish and propose an analytical framework for testing hypotheses regarding this behavior and its effects on the genetic structure of fish in an explicit geographical context, using a geographical genetics toolbox. Although disentangling the many potential causes underlying genetic population structure and unambiguously demonstrating that the homing behavior causes these genetic patterns is difficult, our framework allows the successive testing of homing behavior with increasing levels of complexity based on the following: (1) establishment of population structures among waterheads; (2) patterns of genetic variability throughout the adult migratory pool; (3) analyses of the non-migratory adult pool; and (4) comparisons among successive generations. We expect that the framework presented here will help delineating the appropriate uses of different sampling designs to make inferences regarding homing behavior and illustrate the limits imposed by the interpretation of different types of genetic data. More importantly, we hope this framework enables researchers to understand how a particular dataset can be utilized in a broader context as an ongoing part of a larger research program and thus guide future research by developing better and more integrated sampling designs.

BarcodingGO: A problem-based approach to teach concepts related to environmental-DNA and bioinformatics.

In this paper, we propose and describe a new approach, named BarcodingGO, to teach environmental DNA and bioinformatics concepts to undergraduate or graduate students in molecular biology-related fields. The learning pipeline proposed here aims to solve a simulated environmental monitoring problem, in which a biodiversity survey of a particular region is needed to assess the impact of an environmental disaster. Biological surveys, in the context of environmental DNA studies, are performed by analyzing the DNA released by organisms living in a specific environment. We proposed a scenario in which quick response (QR) codes represented a given environmental DNA, and they were positioned in a scattered pattern across two regions of the classroom (representing pre and post scenarios for a particular environmental disaster). The QR codes redirect to a page that contained a fictional representation of an animal or a plant. Students then survey the region's biodiversity using QR code scanning applications on their cell phones by "capturing" these organisms as an analogy to the Pokémon GO game of the international Pokémon franchise. We believe this method (or even an adaptation of it) can be an essential tool to engage students in molecular biology classes. Moreover, this approach can help to teach how modern genomics and bioinformatics tools can be applied to solve real problems in conservation biology.

Genetic analysis of a local population of Oryza glumaepatula using SSR markers: implications for management and conservation programs.

Knowledge of natural diversity and population structures of wild species, which might be related to cultivated species, is fundamental for conservation and breeding purposes. In this study, a genetic characterization of a large population of Oryza glumaepatula, occurring in a 10 km(2) area located at Tamengo Basin (Paraguay River, Brazil), was performed using SSR markers. This population is annually dragged from the river to permit navigation; one goal of this study was to examine the impact of this removal on genetic variability. From 18 polymorphic SSR markers, a total of 190 alleles were detected in a sample of 126 individuals, with an average of 10.3 alleles/locus, and a H(e) of 0.67. The five QTL-related markers showed an average H(e) value of 0.56, while the remaining 13 markers detected an average estimate of 0.70. An apparent outcrossing rate of 30%, a high proportion of alleles at low frequencies (56%), and the presence of exclusive alleles (9.5%) were found, with strong evidence of the establishment of individuals from different populations upstream in the Paraguay River. For conservation purposes, the river drag has no effect on the population. However, periodical seed collection from the Corumbá population can preserve part of the genetic variability present in upstream populations reducing the need for upriver collecting expeditions.

A review of techniques for spatial modeling in geographical, conservation and landscape genetics.

Most evolutionary processes occur in a spatial context and several spatial analysis techniques have been employed in an exploratory context. However, the existence of autocorrelation can also perturb significance tests when data is analyzed using standard correlation and regression techniques on modeling genetic data as a function of explanatory variables. In this case, more complex models incorporating the effects of autocorrelation must be used. Here we review those models and compared their relative performances in a simple simulation, in which spatial patterns in allele frequencies were generated by a balance between random variation within populations and spatially-structured gene flow. Notwithstanding the somewhat idiosyncratic behavior of the techniques evaluated, it is clear that spatial autocorrelation affects Type I errors and that standard linear regression does not provide minimum variance estimators. Due to its flexibility, we stress that principal coordinate of neighbor matrices (PCNM) and related eigenvector mapping techniques seem to be the best approaches to spatial regression. In general, we hope that our review of commonly used spatial regression techniques in biology and ecology may aid population geneticists towards providing better explanations for population structures dealing with more complex regression problems throughout geographic space.

The complete mitochondrial genome of the aquatic coralsnake Micrurus surinamensis (Reptilia, Serpentes, Elapidae).

In this study, we report the first complete mitochondrial genome sequence of the Aquatic Coralsnake Micrurus surinamensis. The mitochondrial genome lengthis 17,375 bp, comprising 13 protein-coding genes, 2 rRNA (12S and 16S) and 22 tRNA, as well as two typical control regions. Phylogenetic analysis based upon 13 protein-coding genes showed clusters based on terrestrial and marine species.

Phylogeography of Tibouchina papyrus (Pohl) Toledo (Melastomataceae), an endangered tree species from rocky savannas, suggests bidirectional expansion due to climate cooling in the Pleistocene.

Many endemic species present disjunct geographical distribution; therefore, they are suitable models to test hypotheses about the ecological and evolutionary mechanisms involved in the origin of disjunct distributions in these habitats. We studied the genetic structure and phylogeography of Tibouchina papyrus (Melastomataceae), endemic to rocky savannas in Central Brazil, to test hypothesis of vicariance and dispersal in the origin of the disjunct geographical distribution. We sampled 474 individuals from the three localities where the species is reported: Serra dos Pirineus, Serra Dourada, and Serra de Natividade. Analyses were based on the polymorphisms at cpDNA and on nuclear microsatellite loci. To test for vicariance and dispersal we constructed a median-joining network and performed an analysis of molecular variance (AMOVA). We also tested population bottleneck and estimated demographic parameters and time to most recent common ancestor (TMRCA) using coalescent analyses. A remarkable differentiation among populations was found. No significant effect of population expansion was detected and coalescent analyses showed a negligible gene flow among populations and an ancient coalescence time for chloroplast genome. Our results support that the disjunct distribution of T. papyrus may represent a climatic relict. With an estimated TMRCA dated from ∼836.491 ± 107.515 kyr BP (before present), we hypothesized that the disjunct distribution may be the outcome of bidirectional expansion of the geographical distribution favored by the drier and colder conditions that prevailed in much of Brazil during the Pre-Illinoian glaciation, followed by the retraction as the climate became warmer and moister.

Landscape conservation genetics of Dipteryx alata ("baru" tree: Fabaceae) from Cerrado region of central Brazil.

In this paper random amplified polymorphic DNA (RAPD) was used to evaluate the degree of among-population differentiation and associated spatial patterns of genetic divergence for Dipteryx alata Vogel populations from Cerrado region of central Brazil, furnishing support for future programs of conservation of this species. We analyzed patterns of genetic and spatial population structure using 45 RAPD loci scored for 309 trees, sampled from five different regions with two populations each. Genetic structure analysis suggested that panmixia null hypothesis can be rejected, with significant among-population components of 15%. Hierarchical partition by Analysis of Molecular Variance (AMOVA) shows that 5% of genetic variation is within regions, whereas 10% of variation is among regions, and these results were confirmed by a Bayesian analyses on HICKORY. The Mantel correlogram revealed that this divergence is spatially structured, so that local populations situated at short geographic distances could not be considered independent units for conservation and management. However, genetic discontinuities among populations were found in the northwest and southeast parts of the study area, corresponding to regions of recent socio-economic expansion and high population density, respectively. Taking both geographic distances and genetic discontinuities into account it is possible to establish a group of population to be conserved, covering most of D. alata geographic distribution and congruent with previously established priority areas for conservation in the Cerrado region.