Development of genomic tools in a widespread tropical tree, Symphonia globulifera L.f.: a new low-coverage draft genome, SNP and SSR markers

Chromosome-Level Genome Assembly and Genetic Maker System of the Endangered Largemouth Bronze Gudgeon (Coreius guichenoti) with Focus on Conservation Applications

The largemouth bronze gudgeon (Coreius guichenoti), an endemic fish species, is distributed in the upper Yangtze River drainage. Due to anthropogenetic factors such as water pollution, overfishing, and dam construction, the wild populations of C. guichenoti have dramatically declined in recent decades. In this study, we generated a reference chromosomal-level genome assembly of C. guichenoti on the basis of PacBio HiFi sequencing and Hi-C scaffolding technologies. The final genome assembly was 1.10 Gb in length (contig N50: 28.64 Mb; scaffold N50: 42.39 Mb) with 25 chromosomes. The completeness score of the C. guichenoti genome was 96.4%, and high synteny was detected compared with Danio rerio and Ictalurus punctatus genomes. A total of 24 325 PCGs were annotated for the C. guichenoti genome. Comparative genomics analysis identified 986 expanded gene families in C. guichenoti, which were significantly enriched in GO items associated with the development and interaction of sperm and egg as well as immunity. Furthermore, positively selected genes (PSGs) detected in C. guichenoti were mainly associated with DNA repair, ATP binding, mitochondrion, and lipid homeostasis. Based on the reference genome and resequencing data, the polymorphic microsatellite (SSR) loci were comprehensively analyzed for C. guichenoti, and the top 15 tetra-nucleotide SSR loci were selected for the construction of the genetic maker system after validation through PCR and genotyping. All of these 15 tetra-nucleotide SSR loci without Hardy-Weinberg equilibrium (HWE) deviation showed high polymorphism, good amplification stability, and selective neutrality. The PID (sibs) curves revealed that the subset of four tetra-nucleotide SSR loci (cgui1, cgui5, cgui3, cgui13) was sufficient for accurate identification of C. guichenoti individuals (PIDsib < 0.01). These 15 tetra-nucleotide SSR loci could also serve as genetic markers in subsequent parentage identification and genetic diversity analysis. The chromosome-level genome assembly and findings laid solid foundations for molecular breeding, genomic research, and biological conservation of C. guichenoti.

Selection in space and time: Individual tree growth is adapted to tropical forest gap dynamics

Tropical forest dynamics are driven by growth and survival strategies of tree species in relation to treefall gaps; however, the ecological and evolutionary roles of intraspecific variation in the response to forest gaps remain unexplored. Here, we associated genomic data of three related tree species of the genus Symphonia in a French Guiana forest with (1) each individual tree's growth potential, and (2) with its light and competition environment estimated based on 33 years of forest monitoring in plots covering 120 ha. We show that individual trees within species have genetically determined growth strategies that are adapted to the local light and competition environments, which are shaped by the time since the last treefall. Within species, fast-growing genotypes are more frequent in light-enriched environments and slow-growing genotypes in more shaded environments. Forest gap dynamics is thus a strong selection driver that shapes adaptive strategies and maintains genetic variation within tropical tree species.

Topography drives microgeographic adaptations of closely related species in two tropical tree species complexes

Closely related tree species that grow in sympatry are abundant in rainforests. However, little is known of the eco-evolutionary processes that govern their niches and local coexistence. We assessed genetic species delimitation in closely related sympatric species belonging to two Neotropical tree species complexes and investigated their genomic adaptation to a fine-scale topographic gradient with associated edaphic and hydrologic features. Combining LiDAR-derived topography, tree inventories, and single nucleotide polymorphisms (SNPs) from gene capture experiments, we explored genome-wide population genetic structure, covariation of environmental variables, and genotype-environment association to assess microgeographic adaptations to topography within the species complexes Symphonia (Clusiaceae), and Eschweilera (Lecythidaceae) with three species per complex and 385 and 257 individuals genotyped, respectively. Within species complexes, closely related tree species had different realized optima for topographic niches defined through the topographic wetness index or the relative elevation, and species displayed genetic signatures of adaptations to these niches. Symphonia species were genetically differentiated along water and nutrient distribution particularly in genes responding to water deprivation, whereas Eschweilera species were genetically differentiated according to soil chemistry. Our results suggest that varied topography represents a powerful driver of processes modulating tropical forest biodiversity with differential adaptations that stabilize local coexistence of closely related tree species.

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.

Isolation of microsatellite loci in the African tree species Staudtia kamerunensis (Myristicaceae) using high-throughput sequencing

Staudtia kamerunensis (Myristicaceae) or 'Niové' is an evergreen tree widespread in Central African moist forests. The bark and seeds are used in traditional medicine, yet the tree is mainly harvested for its high quality, multi-purpose timber. To facilitate sustainable harvesting and conservation of the species, we aim to develop microsatellite markers that can be used to study the mating system, gene flow, genetic diversity and population structure. Genomic DNA of S. kamerunensis was sequenced on an Illumina MiSeq platform, generating 195,720 paired-end reads with 3671 sequences containing microsatellites. Amplification tests resulted in the development of 16 highly polymorphic microsatellite loci of which 14 were tested in 183 individuals of S. kamerunensis from three populations. The number of detected alleles per locus ranged from 15 to 39 and the average observed and expected heterozygosity across loci and populations were Ho = 0.713 (0.14-0.97) and He = 0.879 (0.19-0.95) respectively. The high levels of polymorphism observed in the newly developed microsatellite markers demonstrate their usefulness to study gene flow, population structure and spatial distribution of genetic diversity in S. kamerunensis.

Altitudinal gradients, biogeographic history and microhabitat adaptation affect fine-scale spatial genetic structure in African and Neotropical populations of an ancient tropical tree species

The analysis of fine-scale spatial genetic structure (FSGS) within populations can provide insights into eco-evolutionary processes. Restricted dispersal and locally occurring genetic drift are the primary causes for FSGS at equilibrium, as described in the isolation by distance (IBD) model. Beyond IBD expectations, spatial, environmental or historical factors can affect FSGS. We examined FSGS in seven African and Neotropical populations of the late-successional rain forest tree Symphonia globulifera L. f. (Clusiaceae) to discriminate the influence of drift-dispersal vs. landscape/ecological features and historical processes on FSGS. We used spatial principal component analysis and Bayesian clustering to assess spatial genetic heterogeneity at SSRs and examined its association with plastid DNA and habitat features. African populations (from Cameroon and São Tomé) displayed a stronger FSGS than Neotropical populations at both marker types (mean Sp = 0.025 vs. Sp = 0.008 at SSRs) and had a stronger spatial genetic heterogeneity. All three African populations occurred in pronounced altitudinal gradients, possibly restricting animal-mediated seed dispersal. Cyto-nuclear disequilibria in Cameroonian populations also suggested a legacy of biogeographic history to explain these genetic patterns. Conversely, Neotropical populations exhibited a weaker FSGS, which may reflect more efficient wide-ranging seed dispersal by Neotropical bats and other dispersers. The population from French Guiana displayed an association of plastid haplotypes with two morphotypes characterized by differential habitat preferences. Our results highlight the importance of the microenvironment for eco-evolutionary processes within persistent tropical tree populations.

Advances in ecological genomics in forest trees and applications to genetic resources conservation and breeding

Forest trees are an unparalleled group of organisms in their combined ecological, economic and societal importance. With widespread distributions, predominantly random mating systems and large population sizes, most tree species harbour extensive genetic variation both within and among populations. At the same time, demographic processes associated with Pleistocene climate oscillations and land-use change have affected contemporary range-wide diversity and may impinge on the potential for future adaptation. Understanding how these adaptive and neutral processes have shaped the genomes of trees species is therefore central to their management and conservation. As for many other taxa, the advent of high-throughput sequencing methods is expected to yield an understanding of the interplay between the genome and environment at a level of detail and depth not possible only a few years ago. An international conference entitled 'Genomics and Forest Tree Genetics' was held in May 2016, in Arcachon (France), and brought together forest geneticists with a wide range of research interests to disseminate recent efforts that leverage contemporary genomic tools to probe the population, quantitative and evolutionary genomics of trees. An important goal of the conference was to discuss how such data can be applied to both genome-enabled breeding and the conservation of forest genetic resources under land use and climate change. Here, we report discoveries presented at the meeting and discuss how the ecological genomic toolkit can be used to address both basic and applied questions in tree biology.