Demographic History of Shorea curtisii (Dipterocarpaceae) Inferred from Chloroplast DNA Sequence Variations

The genome of Shorea leprosula (Dipterocarpaceae) highlights the ecological relevance of drought in aseasonal tropical rainforests

Hyperdiverse tropical rainforests, such as the aseasonal forests in Southeast Asia, are supported by high annual rainfall. Its canopy is dominated by the species-rich tree family of Dipterocarpaceae (Asian dipterocarps), which has both ecological (e.g., supports flora and fauna) and economical (e.g., timber production) importance. Recent ecological studies suggested that rare irregular drought events may be an environmental stress and signal for the tropical trees. We assembled the genome of a widespread but near threatened dipterocarp, Shorea leprosula, and analyzed the transcriptome sequences of ten dipterocarp species representing seven genera. Comparative genomic and molecular dating analyses suggested a whole-genome duplication close to the Cretaceous-Paleogene extinction event followed by the diversification of major dipterocarp lineages (i.e. Dipterocarpoideae). Interestingly, the retained duplicated genes were enriched for genes upregulated by no-irrigation treatment. These findings provide molecular support for the relevance of drought for tropical trees despite the lack of an annual dry season.

Genetic Diversity and Demographic History of an Upper Hill Dipterocarp (Shorea platyclados): Implications for Conservation

Southeast Asian rainforests at upper hill elevations are increasingly vulnerable to degradation because most lowland forest areas have been converted to different land uses. As such, understanding the genetics of upper hill species is becoming more crucial for their future management and conservation. Shorea platyclados is an important, widespread upper hill dipterocarp in Malaysia. To elucidate the genetic structure of S. platyclados and ultimately provide guidelines for a conservation strategy for this species, we carried out a comprehensive study of the genetic diversity and demographic history of S. platyclados. Twenty-seven populations of S. platyclados across its range in Malaysia were genotyped at 15 polymorphic microsatellite loci and sequenced at seven noncoding chloroplast DNA (cpDNA) regions. A total of 303 alleles were derived from the microsatellite loci, and 29 haplotypes were identified based on 2892 bp of concatenated cpDNA sequences. The populations showed moderately high genetic diversity (mean HE = 0.680 for microsatellite gene diversity and HT = 0.650 for total haplotype diversity) and low genetic differentiation (FST = 0.060). Bayesian clustering divided the studied populations into two groups corresponding to western and eastern Malaysia. Bottleneck analysis did not detect any recent bottleneck events. Extended Bayesian skyline analyses showed a model of constant size for the past population history of this species. Based on our findings, priority areas for in situ and ex situ conservation and a minimum population size are recommended for the sustainable utilization of S. platyclados.

Potential of Genome-Wide Association Studies and Genomic Selection to Improve Productivity and Quality of Commercial Timber Species in Tropical Rainforest, a Case Study of Shorea platyclados

Shorea platyclados (Dark Red Meranti) is a commercially important timber tree species in Southeast Asia. However, its stocks have dramatically declined due, inter alia, to excessive logging, insufficient natural regeneration and a slow recovery rate. Thus, there is a need to promote enrichment planting and develop effective technique to support its rehabilitation and improve timber production through implementation of Genome-Wide Association Studies (GWAS) and Genomic Selection (GS). To assist such efforts, plant materials were collected from a half-sib progeny population in Sari Bumi Kusuma forest concession, Kalimantan, Indonesia. Using 5900 markers in sequences obtained from 356 individuals, we detected high linkage disequilibrium (LD) extending up to >145 kb, suggesting that associations between phenotypic traits and markers in LD can be more easily and feasibly detected with GWAS than with analysis of quantitative trait loci (QTLs). However, the detection power of GWAS seems low, since few single nucleotide polymorphisms linked to any focal traits were detected with a stringent false discovery rate, indicating that the species’ phenotypic traits are mostly under polygenic quantitative control. Furthermore, Machine Learning provided higher prediction accuracies than Bayesian methods. We also found that stem diameter, branch diameter ratio and wood density were more predictable than height, clear bole, branch angle and wood stiffness traits. Our study suggests that GS has potential for improving the productivity and quality of S. platyclados, and our genomic heritability estimates may improve the selection of traits to target in future breeding of this species.

Nuclear and chloroplast DNA phylogeography reveals Pleistocene divergence and subsequent secondary contact of two genetic lineages of the tropical rainforest tree species Shorea leprosula (Dipterocarpaceae) in South-East Asia

Tropical rainforests in South-East Asia have been affected by climatic fluctuations during past glacial eras. To examine how the accompanying changes in land areas and temperature have affected the genetic properties of rainforest trees in the region, we investigated the phylogeographic patterns of a widespread dipterocarp species, Shorea leprosula. Two types of DNA markers were used: expressed sequence tag-based simple sequence repeats and chloroplast DNA (cpDNA) sequence variations. Both sets of markers revealed clear genetic differentiation between populations in Borneo and those in the Malay Peninsula and Sumatra (Malay/Sumatra). However, in the south-western part of Borneo, genetic admixture of the lineages was observed in the two marker types. Coalescent simulation based on cpDNA sequence variation suggested that the two lineages arose 0.28-0.09 million years before present and that following their divergence migration from Malay/Sumatra to Borneo strongly exceeded migration in the opposite direction. We conclude that the genetic structure of S. leprosula was largely formed during the middle Pleistocene and was subsequently modified by eastward migration across the subaerially exposed Sunda Shelf.

Contrasting cpDNA variation in two Indonesian endemic lowland dipterocarp species and implications for their conservation

Shorea javanica (Dipterocarpaceae) is an economically important dammar-producing tree, endemic to the tropical lowland forests of Sumatra and Java, Indonesia. However, its total population size is limited and this species is endangered. Shorea selanica is one of the very limited numbers of species in genus Shorea (Dipterocarpaceae) that grows in Wallacean. This species can be found only in the central part of the Moluccas, eastern Indonesia. Six populations (77 individuals in total) were sampled for S. javanica while three populations (27 individuals in total) were sampled for S. selanica. To determine genetic variation and population structure, three non-coding chloroplast DNA regions of trnL-trnF, psbC-trnS, trnS-trnfM and two non-coding chloroplast DNA regions of trnT-trnL, trnL-trnF were sequenced from S. javanica and S. selanica, respectively. There was no variation in the chloroplast DNA regions from S. javanica, except for one unique indel polymorphism. Nucleotide diversity within S. selanica populations ranged from 0 (Seram) to 0.00044 (Buru), with a pooled value of 0.00041. S. javanica was determined as having no population structure while high levels of genetic differentiation was found among populations of S. selanica (F(ST) = 0.702). Different pattern of population structure among the two species in this study suggested the need for distinct management and conservation strategies for each species. For S. javanica, connectivity within and among populations, including augmentation of population size across the species range, should be promoted. For S. selanica, an in situ conservation plan that defines core areas completely free from perturbation within each population is necessary.