Chloroplast genome assembly of Serjania erecta Raldk: comparative analysis reveals gene number variation and selection in protein-coding plastid genes of Sapindaceae.

Serjania erecta Raldk is an essential genetic resource due to its anti-inflammatory, gastric protection, and anti-Alzheimer properties. However, the genetic and evolutionary aspects of the species remain poorly known. Here, we sequenced and assembled the complete chloroplast genome of S. erecta and used it in a comparative analysis within the Sapindaceae family. S. erecta has a chloroplast genome (cpDNA) of 159,297 bp, divided into a Large Single Copy region (LSC) of 84,556 bp and a Small Single Copy region (SSC) of 18,057 bp that are surrounded by two Inverted Repeat regions (IRa and IRb) of 28,342 bp. Among the 12 species used in the comparative analysis, S. erecta has the fewest long and microsatellite repeats. The genome structure of Sapindaceae species is relatively conserved; the number of genes varies from 128 to 132 genes, and this variation is associated with three main factors: (1) Expansion and retraction events in the size of the IRs, resulting in variations in the number of rpl22, rps19, and rps3 genes; (2) Pseudogenization of the rps2 gene; and (3) Loss or duplication of genes encoding tRNAs, associated with the duplication of trnH-GUG in X. sorbifolium and the absence of trnT-CGU in the Dodonaeoideae subfamily. We identified 10 and 11 mutational hotspots for Sapindaceae and Sapindoideae, respectively, and identified six highly diverse regions (tRNA-Lys - rps16, ndhC - tRNA-Val, petA - psbJ, ndhF, rpl32 - ccsA, and ycf1) are found in both groups, which show potential for the development of DNA barcode markers for molecular taxonomic identification of Serjania. We identified that the psaI gene evolves under neutrality in Sapindaceae, while all other chloroplast genes are under strong negative selection. However, local positive selection exists in the ndhF, rpoC2, ycf1, and ycf2 genes. The genes ndhF and ycf1 also present high nucleotide diversity and local positive selection, demonstrating significant potential as markers. Our findings include providing the first chloroplast genome of a member of the Paullinieae tribe. Furthermore, we identified patterns in variations in the number of genes and selection in genes possibly associated with the family's evolutionary history.

The complete chloroplast genome sequence of Eugenia klotzschiana O. Berg unveils the evolutionary dynamics in plastomes of Myrteae DC. Tribe (Myrtaceae).

Myrteae is the most diversified tribe in the Myrtaceae family and has great ecological and economic importance. Here, we performed the assembly and annotation of the chloroplast genome of Eugenia klotzschiana O. Berg and used this in a comparative analysis with other 13 species from the Myrteae tribe. The E. klotzschiana plastome exhibited a length of 158,977 bp and a very conserved structure and gene composition when compared with other Myrteae genomes. We identified 34 large repetitive sequences and 94 SSR repeats in E. klotzschiana plastome. The trnT-trnL, rpl32-trnL, ndhF-rpl32, psbE-petL, and ycf1 regions were identified as mutational hotspots. A negative selection signal was detected in 74 protein-coding genes while neutral evolution was detected in two genes (rps12 and psaI). Furthermore, 222 RNA editing sites were identified in the E. klotzschiana plastome. We also obtained a plastome-based Myrtales phylogenetic tree, including E. klotzschiana for the first time in a molecular phylogeny, recovering its sister relationship for all other Eugenia species. Our results illuminate how evolution shaped the chloroplast genome structure and composition in the Myrteae tribe, especially in the E. klotzschiana plastome.

Data of SSRs primers for high-throughput genotyping-by-sequencing (SSR-Seq) based on the partial genome assembly of Eugenia klotzschiana (Myrtaceae).

The neotropical fruit plant Eugenia klotzschiana Berg. is endemic from South America and occurs in the Brazilian savannah areas, a biome threatened by intensive agriculture. This species is a plant listed on the Brazilian list of Plants for the Future. The E. klotzschiana fruits have great nutritional value and antioxidant activity and are consumed in natura or processed into juice or jelly. However, their harvest is predominantly in native areas and needs further studies for large-scale commercialization. Nuclear genomic data and population genetic tools are still quite scarce for the species. Here, we provide data on the first partially assembled genome of E. klotzschiana (211 Mbp, ∼75.16% genome coverage, N50 = 3,407, and 46.8% BUSCO completeness), the raw Illumina sequencing reads, and two sets of primers for microsatellite (SSRs) high-throughput genotyping-by-sequencing (SSR-Seq) identified in the nuclear genome. These genomic resources are fundamental for this species conservation strategies and the development of a future breeding program.

Transferability and characterization of microsatellite markers from Byrsonima cydoniifolia A. Juss. (MALPIGHIACEAE) in seven related taxa from Cerrado biome reveal genetic relationships.

Byrsonima Rich. is one of the largest genera of the Malpighiaceae family with 97 species occurrence in Brazil and multiple potentialities, including pharmaceutical and food industries. In this study, 17 microsatellite markers characterized in Byrsonima cydoniifolia were tested for seven related taxa, all species are native to Brazil and four are endemic. Genomic DNA was extracted from leaves tissues and 17 microsatellite markers were used to cross-amplification of microsatellite regions. Polymorphism and genetic diversity were evaluated for B. intermedia, B. verbascifolia, B. laxiflora, B. subterranea, B. umbellata, B. linearifolia. from 16 individuals and for B. viminifolia from 14 individuals. Transferred microsatellite markers panels ranged from 11 (64.8%) in B. viminifolia to 6 (35.2%) in B. umbellata. The total number of alleles per locus ranged from 5 (B. linearifolia) to 8 (B. subterranea) alleles. B. umbellata showed lower values of observed and expected heterozygosity (HO = 0.312; HE = 0.436) and B. subterranea presented the highest values (HO = 0.687; HE = 0.778). A greater number of microsatellite markers should be developed for B. umbellata. The microsatellite marker panels transferred to the species B. intermedia, B. verbascifolia, B. laxiflora, B. subterranea, B. viminifolia and B. linearifolia are very informative, with a high combined probability of exclusion of paternity (Q ≥ 0.976) and the low combined probability of identity (I ≤ 9.91 × 10-6), potentially suitable for future genetic-population studies, supporting strategies for maintaining the genetic diversity and for exploration of Byrsonima species as genetic resources.

Genetic structure and chemical diversity in natural populations of Uncaria guianensis (Aubl.) J.F.Gmel. (Rubiaceae).

Uncaria guianensis is native to the Amazon and is used traditionally as an anti-inflammatory. Natural populations of the species have declined markedly in recent times because of strong anthropic pressure brought about by deforestation and indiscriminate collection. The aim of the present study was to assess the genetic and chemical diversity among eight natural populations of U. guianensis located in the Brazilian states of Acre, Amapá and Amazonas. A set of four primer combinations was employed in sequence-related amplified polymorphism (SRAP) amplifications of leaf DNA, and the fragments were analyzed in an LI-COR model 4300 DNA Analyzer. Genetic variability within the populations (81%) was substantially greater than that detected between them (19%). The highest percentage of polymorphic loci (90.21%) and the largest genetic variability were observed in the population located in Mazagão, Amapá. Genetic differentiation between populations was high (Fst = 0.188) and the studied populations formed three distinct genetic groups (K = 3). The population located in Assis Brasil, Acre, presented the highest average content of the mitraphylline (0.60 mg/g dry weight,). However, mitraphylline and isomitraphylline not detected in most individuals in the studied populations, and it is questionable whether they should be considered as chemical markers of the species. The genetic data confirm the urgent need for conservation programs for U. guianensis, and for further studies aimed at ascertaining the genetic basis and heritability of alkaloid accumulation.

Unravelling the genetic differentiation among varieties of the Neotropical savanna tree Hancornia speciosa Gomes.

Background and Aims: Spatial distribution of species genetic diversity is often driven by geographical distance (isolation by distance) or environmental conditions (isolation by environment), especially under climate change scenarios such as Quaternary glaciations. Here, we used coalescent analyses coupled with ecological niche modelling (ENM), spatially explicit quantile regression analyses and the multiple matrix regression with randomization (MMRR) approach to unravel the patterns of genetic differentiation in the widely distributed Neotropical savanna tree, Hancornia speciosa (Apocynaceae). Due to its high morphological differentiation, the species was originally classified into six botanical varieties by Monachino, and has recently been recognized as only two varieties by Flora do Brasil 2020. Thus, H. speciosa is a good biological model for learning about evolution of phenotypic plasticity under genetic and ecological effects, and predicting their responses to changing environmental conditions. Methods: We sampled 28 populations (777 individuals) of Monachino's four varieties of H. speciosa and used seven microsatellite loci to genotype them. Key Results: Bayesian clustering showed five distinct genetic groups (K = 5) with high admixture among Monachino's varieties, mainly among populations in the central area of the species geographical range. Genetic differentiation among Monachino's varieties was lower than the genetic differentiation among populations within varieties, with higher within-population inbreeding. A high historical connectivity among populations of the central Cerrado shown by coalescent analyses may explain the high admixture among varieties. In addition, areas of higher climatic suitability also presented higher genetic diversity in such a way that the wide historical refugium across central Brazil might have promoted the long-term connectivity among populations. Yet, FST was significantly related to geographic distances, but not to environmental distances, and coalescent analyses and ENM predicted a demographical scenario of quasi-stability through time. Conclusions: Our findings show that demographical history and isolation by distance, but not isolation by environment, drove genetic differentiation of populations. Finally, the genetic clusters do not support the two recently recognized botanical varieties of H. speciosa, but partially support Monachino's classification at least for the four sampled varieties, similar to morphological variation.

Demographic stability and high historical connectivity explain the diversity of a savanna tree species in the Quaternary.

Background and Aims: Cyclic glaciations were frequent throughout the Quaternary and this affected species distribution and population differentiation worldwide. The present study reconstructed the demographic history and dispersal routes of Eugenia dysenterica lineages and investigated the effects of Quaternary climate change on its spatial pattern of genetic diversity. Methods: A total of 333 individuals were sampled from 23 populations and analysed by sequencing four regions of the chloroplast DNA and the internal transcribed spacer of the nuclear DNA. The analyses were performed using a multi-model inference approach based on ecological niche modelling and statistical phylogeography. Key Results: Coalescent simulation showed that population stability through time is the most likely scenario. The palaeodistribution dynamics predicted by the ecological niche models revealed that the species was potentially distributed across a large area, extending over Central-Western Brazil through the last glaciation. The lineages of E. dysenterica dispersed from Central Brazil towards populations at the northern, western and south-eastern regions. A historical refugium through time may have favoured lineage dispersal and the maintenance of genetic diversity. Conclusions: The results suggest that the central region of the Cerrado biome is probably the centre of distribution of E. dysenterica and that the spatial pattern of its genetic diversity may be the outcome of population stability throughout the Quaternary. The lower genetic diversity in populations in the south-eastern Cerrado biome is probably due to local climatic instability during the Quaternary.

Mating system and pollen dispersal in Eugenia dysenterica (Myrtaceae) germplasm collection: tools for conservation and domestication.

Eugenia dysenterica DC. (Myrtaceae) is a perennial tree producing edible fruits and ornamental flowers of potential value widely distributed in Brazilian "Cerrados" (savannas), but available genetic resources and potential for future breeding programs must be evaluated. Here we evaluated the reproductive system and pollen-mediated gene flow in one generation of Eugenia dysenterica germplasm collection of Agronomy School, Federal University of Goiás (in Goiânia city, Central Brazil). We collected leaves from all adults from the germplasm collection (682 plants) and seeds (542) from 23 mother-trees. Genotypes were obtained for seven microsatellite loci. Genetic diversity was high and did not significantly differ between adults (H e = 0.777) and progeny arrays (H e = 0.617). Our results showed that E. dysenterica has an allogamous mating system in the germplasm collection (t m = 0.957), but with high and significant biparental inbreeding (t m - t s = 0.109). Because sibs are very close to each other, mating between closely related individuals is likely. Paternity correlation was also relatively high, indicating a 11.9 % probability that a randomly chosen pair of outcrossed progeny from the same array are full sibs. The maximum pollen dispersal distance (224 m), estimated using assignment test, corresponded to the boundaries of the orchard. We were able to assign the paternity to only 64 % of the 349 seeds analyzed, indicating potential pollen immigration to the germplasm collection. The variance effective population size estimated for one maternal family in the germplasm collection (N ev = 3.42) is very close to the theoretical maximum value for half-sibs (Nev = 4.0). Because E. dysenterica has a long life cycle and generation time, the maintenance of an effective population size of at least 100 in the germplasm collection is suggested, which can be achieved by maintaining a seed-trees number around 30 individuals.

Recovering species demographic history from multi-model inference: the case of a Neotropical savanna tree species.

BACKGROUND: Glaciations were recurrent throughout the Quaternary and potentially shaped species genetic structure worldwide by affecting population dynamics. Here, we implemented a multi-model inference approach to recover the distribution dynamics and demographic history of a Neotropical savanna tree, Tabebuia aurea (Bignoniaceae). Exploring different algorithms and paleoclimatic simulations, we used ecological niche modelling to generate alternative hypotheses of potential demographic changes through the last glacial cycle and estimated genetic parameters using coalescent modelling. RESULTS: Comparing predictions from demographic hypotheses with genetic parameters of modern populations, our findings revealed a likely scenario of population decline, with spatial displacement towards Northeast Brazil from the last glacial maximum to the mid-Holocene. Subsequently, populations expanded in response to the return of the climatically suitable conditions in Central-West Brazil. Nevertheless, a wide historical refugium across Central Brazil likely maintained large populations connected throughout time. The expected genetic signatures from such predicted distribution dynamics are also corroborated by spatial genetic structure observed in modern populations. CONCLUSION: By exploring uncertainties inherent in multiple working hypotheses, we have shown that multi-model inference is a fruitful and efficient approach to recover the nature, timing and geographical context of the Tabebuia aurea population dynamic in response to the Quaternary climate changes.

Paternity testing and behavioral ecology: a case study of jaguars (Panthera onca) in Emas National Park, Central Brazil

We used microsatellite loci to test the paternity of two male jaguars involved in an infanticide event recorded during a long-term monitoring program of this species. Seven microsatellite primers originally developed for domestic cats and previously selected for Panthera onca were used. In order to deal with uncertainty in the mother's genotypes for some of the loci, 10000 values of W were derived by simulation procedures. The male that killed the two cubs was assigned as the true sire. Although the reasons for this behavior remain obscure, it shows, in principle, a low recognition of paternity and kinship in the species. Since the two cubs were not very young, one possibility is that the adult male did not recognize the cubs and killed them for simple territorial reasons. Thus, ecological stress in this local population becomes a very plausible explanation for this infanticide, without further sociobiological implications.