Comparative population genomics in Tabebuia alliance shows evidence of adaptation in Neotropical tree species.

The role of natural selection in shaping spatial patterns of genetic diversity in the Neotropics is still poorly understood. Here, we perform a genome scan with 24,751 probes targeting 11,026 loci in two Neotropical Bignoniaceae tree species: Handroanthus serratifolius from the seasonally dry tropical forest (SDTF) and Tabebuia aurea from savannas, and compared with the population genomics of H. impetiginosus from SDTF. OutFLANK detected 29 loci in 20 genes with selection signal in H. serratifolius and no loci in T. aurea. Using BayPass, we found evidence of selection in 335 loci in 312 genes in H. serratifolius, 101 loci in 92 genes in T. aurea, and 448 loci in 416 genes in H. impetiginosus. All approaches evidenced several genes affecting plant response to environmental stress and primary metabolic processes. The three species shared no SNPs with selection signal, but we found SNPs affecting the same gene in pair of species. Handroanthus serratifolius showed differences in allele frequencies at SNPs with selection signal among ecosystems, mainly between Caatinga/Cerrado and Atlantic Forest, while H. impetiginosus had one allele fixed across all populations, and T. aurea had similar allele frequency distribution among ecosystems and polymorphism across populations. Taken together, our results indicate that natural selection related to environmental stress shaped the spatial pattern of genetic diversity in the three species. However, the three species have different geographical distribution and niches, which may affect tolerances and adaption, and natural selection may lead to different signatures due to the differences in adaptive landscapes in different niches.

RNA-Seq reveals different responses to drought in Neotropical trees from savannas and seasonally dry forests.

BACKGROUND: Water is one of the main limiting factors for plant growth and crop productivity. Plants constantly monitor water availability and can rapidly adjust their metabolism by altering gene expression. This leads to phenotypic plasticity, which aids rapid adaptation to climate changes. Here, we address phenotypic plasticity under drought stress by analyzing differentially expressed genes (DEG) in four phylogenetically related neotropical Bignoniaceae tree species: two from savanna, Handroanthus ochraceus and Tabebuia aurea, and two from seasonally dry tropical forests (SDTF), Handroanthus impetiginosus and Handroanthus serratifolius. To the best of our knowledge, this is the first report of an RNA-Seq study comparing tree species from seasonally dry tropical forest and savanna ecosystems. RESULTS: Using a completely randomized block design with 4 species × 2 treatments (drought and wet) × 3 blocks (24 plants) and an RNA-seq approach, we detected a higher number of DEGs between treatments for the SDTF species H. serratifolius (3153 up-regulated and 2821 down-regulated under drought) and H. impetiginosus (332 and 207), than for the savanna species. H. ochraceus showed the lowest number of DEGs, with only five up and nine down-regulated genes, while T. aurea exhibited 242 up- and 96 down-regulated genes. The number of shared DEGs among species was not related to habitat of origin or phylogenetic relationship, since both T. aurea and H impetiginosus shared a similar number of DEGs with H. serratifolius. All four species shared a low number of enriched gene ontology (GO) terms and, in general, exhibited different mechanisms of response to water deficit. We also found 175 down-regulated and 255 up-regulated transcription factors from several families, indicating the importance of these master regulators in drought response. CONCLUSION: Our findings show that phylogenetically related species may respond differently at gene expression level to drought stress. Savanna species seem to be less responsive to drought at the transcriptional level, likely due to morphological and anatomical adaptations to seasonal drought. The species with the largest geographic range and widest edaphic-climatic niche, H. serratifolius, was the most responsive, exhibiting the highest number of DEG and up- and down-regulated transcription factors (TF).

Genome assemblies for two Neotropical trees: Jacaranda copaia and Handroanthus guayacan.

The lack of genomic resources for tropical canopy trees is impeding several research avenues in tropical forest biology. We present genome assemblies for two Neotropical hardwood species, Jacaranda copaia and Handroanthus (formerly Tabebuia) guayacan, that are model systems for research on tropical tree demography and flowering phenology. For each species, we combined Illumina short-read data with in vitro proximity-ligation (Chicago) libraries to generate an assembly. For Jacaranda copaia, we obtained 104X physical coverage and produced an assembly with N50/N90 scaffold lengths of 1.020/0.277 Mbp. For H. guayacan, we obtained 129X coverage and produced an assembly with N50/N90 scaffold lengths of 0.795/0.165 Mbp. J. copaia and H. guayacan assemblies contained 95.8% and 87.9% of benchmarking orthologs, although they constituted only 77.1% and 66.7% of the estimated genome sizes of 799 and 512 Mbp, respectively. These differences were potentially due to high repetitive sequence content (>59.31% and 45.59%) and high heterozygosity (0.5% and 0.8%) in each species. Finally, we compared each new assembly to a previously sequenced genome for Handroanthus impetiginosus using whole-genome alignment. This analysis indicated extensive gene duplication in H. impetiginosus since its divergence from H. guayacan.

Chloroplast genome assembly of Handroanthus impetiginosus: comparative analysis and molecular evolution in Bignoniaceae.

MAIN CONCLUSION: Bignoniaceae species have conserved chloroplast structure, with hotspots of nucleotide diversity. Several genes are under positive selection, and can be targets for evolutionary studies. Bignoniaceae is one of the most species-rich family of woody plants in Neotropical seasonally dry forests. Here we report the assembly of Handroanthus impetiginosus chloroplast genome and evolutionary comparative analyses of ten Bignoniaceae species representing the genera for which whole-genome chloroplast sequences were available. The chloroplast genome of H. impetiginosus is 159,462 bp in size and has a similar structure compared to the other nine species. The total number of genes was slightly variable amongst the Bignoniaceae, ranging from 124 in H. impetiginosus to 144 in Anemopaegma acutifolium. The inverted repeat (IR) size was variable, ranging from 24,657 bp (Tecomaria capensis) to 40,481 bp (A. acutifolium), due to the contraction and retraction at its boundaries. However, gene boundaries were very similar among the ten species. We found 98 forward and palindromic dispersed repeats, and 85 simple sequence repeats (SSRs). In general, chloroplast sequences were highly conserved, with few nucleotide diversity hotspots in the genes accD, clpP, rpoA, ycf1, ycf2. The phylogenetic analysis based on 77 coding genes was highly consistent with Angiosperm Phylogeny Group (APG) IV. Our results also indicate that most genes are under negative selection or neutral evolution. We found no evidence of branch-site selection, implying that H. impetiginosus is not evolving faster than the other species analyzed, notwithstanding we found site positive selection signal in several genes. These genes can provide targets for evolutionary studies in Bignoniaceae and Lamiales species.

Geology and climate influence rhizobiome composition of the phenotypically diverse tropical tree Tabebuia heterophylla.

Plant-associated microbial communities have diverse phenotypic effects on their hosts that are only beginning to be revealed. We hypothesized that morpho-physiological variations in the tropical tree Tabebuia heterophylla, observed on different geological substrates, arise in part due to microbial processes in the rhizosphere. We characterized the microbiota of the rhizosphere and soil communities associated with T. heterophylla trees in high and low altitude sites (with varying temperature and precipitation) of volcanic, karst and serpentine geologies across Puerto Rico. We sampled 6 areas across the island in three geological materials including volcanic, serpentine and karst soils. Collection was done in 2 elevations (>450m and 0-300m high), that included 3 trees for each site and 4 replicate soil samples per tree of both bulk and rhizosphere. Genomic DNA was extracted from 144 samples, and 16S rRNA V4 sequencing was performed on the Illumina MiSeq platform. Proteobacteria, Actinobacteria, and Verrucomicrobia were the most dominant phyla, and microbiomes clustered by geological substrate and elevation. Volcanic samples were enriched in Verrucomicrobia; karst was dominated by nitrogen-fixing Proteobacteria, and serpentine sites harbored the most diverse communities, with dominant Cyanobacteria. Sites with similar climates but differing geologies showed significant differences on rhizobiota diversity and composition demonstrating the importance of geology in shaping the rhizosphere microbiota, with implications for the plant's phenotype. Our study sheds light on the combined role of geology and climate in the rhizosphere microbial consortia, likely contributing to the phenotypic plasticity of the trees.

Design and evaluation of a sequence capture system for genome-wide SNP genotyping in highly heterozygous plant genomes: a case study with a keystone Neotropical hardwood tree genome.

Targeted sequence capture coupled to high-throughput sequencing has become a powerful method for the study of genome-wide sequence variation. Following our recent development of a genome assembly for the Pink Ipê tree (Handroanthus impetiginosus), a widely distributed Neotropical timber species, we now report the development of a set of 24,751 capture probes for single-nucleotide polymorphisms (SNPs) characterization and genotyping across 18,216 distinct loci, sampling more than 10 Mbp of the species genome. This system identifies nearly 200,000 SNPs located inside or in close proximity to almost 14,000 annotated protein-coding genes, generating quality genotypic data in populations spanning wide geographic distances across the species native range. To provide recommendations for future developments of similar systems for highly heterozygous plant genomes we investigated issues such as probe design, sequencing coverage and bioinformatics, including the evaluation of the capture efficiency and a reassessment of the technical reproducibility of the assay for SNPs recall and genotyping precision. Our results highlight the value of a detailed probe screening on a preliminary genome assembly to produce reliable data for downstream genetic studies. This work should inspire and assist the development of similar genomic resources for other orphan crops and forest trees with highly heterozygous genomes.

Genome assembly of the Pink Ipê (Handroanthus impetiginosus, Bignoniaceae), a highly valued, ecologically keystone Neotropical timber forest tree.

Background: Handroanthus impetiginosus (Mart. ex DC.) Mattos is a keystone Neotropical hardwood tree widely distributed in seasonally dry tropical forests of South and Mesoamerica. Regarded as the "new mahogany," it is the second most expensive timber, the most logged species in Brazil, and currently under significant illegal trading pressure. The plant produces large amounts of quinoids, specialized metabolites with documented antitumorous and antibiotic effects. The development of genomic resources is needed to better understand and conserve the diversity of the species, to empower forensic identification of the origin of timber, and to identify genes for important metabolic compounds. Findings: The genome assembly covers 503.7 Mb (N50 = 81 316 bp), 90.4% of the 557-Mbp genome, with 13 206 scaffolds. A repeat database with 1508 sequences was developed, allowing masking of ∼31% of the assembly. Depth of coverage indicated that consensus determination adequately removed haplotypes assembled separately due to the extensive heterozygosity of the species. Automatic gene prediction provided 31 688 structures and 35 479 messenger RNA transcripts, while external evidence supported a well-curated set of 28 603 high-confidence models (90% of total). Finally, we used the genomic sequence and the comprehensive gene content annotation to identify genes related to the production of specialized metabolites. Conclusions: This genome assembly is the first well-curated resource for a Neotropical forest tree and the first one for a member of the Bignoniaceae family, opening exceptional opportunities to empower molecular, phytochemical, and breeding studies. This work should inspire the development of similar genomic resources for the largely neglected forest trees of the mega-diverse tropical biomes.

Two tropical conifers show strong growth and water-use efficiency responses to altered CO2 concentration.

BACKGROUND AND AIMS: Conifers dominated wet lowland tropical forests 100 million years ago (MYA). With a few exceptions in the Podocarpaceae and Araucariaceae, conifers are now absent from this biome. This shift to angiosperm dominance also coincided with a large decline in atmospheric CO2 concentration (ca). We compared growth and physiological performance of two lowland tropical angiosperms and conifers at ca levels representing pre-industrial (280 ppm), ambient (400 ppm) and Eocene (800 ppm) conditions to explore how differences in ca affect the growth and water-use efficiency (WUE) of seedlings from these groups. METHODS: Two conifers (Araucaria heterophylla and Podocarpus guatemalensis) and two angiosperm trees (Tabebuia rosea and Chrysophyllum cainito) were grown in climate-controlled glasshouses in Panama. Growth, photosynthetic rates, nutrient uptake, and nutrient use and water-use efficiencies were measured. KEY RESULTS: Podocarpus seedlings showed a stronger (66 %) increase in relative growth rate with increasing ca relative to Araucaria (19 %) and the angiosperms (no growth enhancement). The response of Podocarpus is consistent with expectations for species with conservative growth traits and low mesophyll diffusion conductance. While previous work has shown limited stomatal response of conifers to ca, we found that the two conifers had significantly greater increases in leaf and whole-plant WUE than the angiosperms, reflecting increased photosynthetic rate and reduced stomatal conductance. Foliar nitrogen isotope ratios (δ15N) and soil nitrate concentrations indicated a preference in Podocarpus for ammonium over nitrate, which may impact nitrogen uptake relative to nitrate assimilators under high ca SIGNIFICANCE: Podocarps colonized tropical forests after angiosperms achieved dominance and are now restricted to infertile soils. Although limited to a single species, our data suggest that higher ca may have been favourable for podocarp colonization of tropical South America 60 MYA, while plasticity in photosynthetic capacity and WUE may help account for their continued persistence under large changes in ca since the Eocene.

Coalescent Simulation and Paleodistribution Modeling for Tabebuia rosealba Do Not Support South American Dry Forest Refugia Hypothesis.

Studies based on contemporary plant occurrences and pollen fossil records have proposed that the current disjunct distribution of seasonally dry tropical forests (SDTFs) across South America is the result of fragmentation of a formerly widespread and continuously distributed dry forest during the arid climatic conditions associated with the Last Glacial Maximum (LGM), which is known as the modern-day dry forest refugia hypothesis. We studied the demographic history of Tabebuia rosealba (Bignoniaceae) to understand the disjunct geographic distribution of South American SDTFs based on statistical phylogeography and ecological niche modeling (ENM). We specifically tested the dry forest refugia hypothesis; i.e., if the multiple and isolated patches of SDTFs are current climatic relicts of a widespread and continuously distributed dry forest during the LGM. We sampled 235 individuals across 18 populations in Central Brazil and analyzed the polymorphisms at chloroplast (trnS-trnG, psbA-trnH and ycf6-trnC intergenic spacers) and nuclear (ITS nrDNA) genomes. We performed coalescence simulations of alternative hypotheses under demographic expectations from two a priori biogeographic hypotheses (1. the Pleistocene Arc hypothesis and, 2. a range shift to Amazon Basin) and other two demographic expectances predicted by ENMs (3. expansion throughout the Neotropical South America, including Amazon Basin, and 4. retraction during the LGM). Phylogenetic analyses based on median-joining network showed haplotype sharing among populations with evidence of incomplete lineage sorting. Coalescent analyses showed smaller effective population sizes for T. roseoalba during the LGM compared to the present-day. Simulations and ENM also showed that its current spatial pattern of genetic diversity is most likely due to a scenario of range retraction during the LGM instead of the fragmentation from a once extensive and largely contiguous SDTF across South America, not supporting the South American dry forest refugia hypothesis.

Development and characterization of microsatellite loci for Tabebuia cassinoides (Bignoniaceae).

Tabebuia cassinoides (Lam.) DC., popularly known as caxeta, is a tree species that belongs to the plant family Bignoniaceae. This species is endemic to the Brazilian Atlantic Forest and is widely exploited commercially. To date, little is known about its genetic structure, preventing the establishment of adequate management plans for this taxon. The objective of this study was to construct a microsatellite-enriched genomic library for T. cassinoides to select polymorphic loci, and standardize polymerase chain reaction amplification conditions. Of the 15 loci examined, 5 were polymorphic. The number of alleles per locus ranged from 2 to 8, with a mean of 4.4. The microsatellite loci described here represent the basis for detailed population genetic studies of this species, which will greatly contribute for the development of better conservation strategies for this taxon.