Conservation genomic investigation of an endangered conifer, Thuja sutchuenensis, reveals low genetic diversity but also low genetic load.

Endangered species generally have small populations with low genetic diversity and a high genetic load. Thuja sutchuenensis is an endangered conifer endemic to southwestern China. It was once considered extinct in the wild, but in 1999 was rediscovered. However, little is known about its genetic load. We collected 67 individuals from five wild, isolated T. sutchuenensis populations, and used 636,151 SNPs to analyze the level of genetic diversity and genetic load in T. sutchuenensis to delineate the conservation units of T. sutchuenensis, based on whole transcriptome sequencing data, as well as target capture sequencing data. We found that populations of T. sutchuenensis could be divided into three groups. These groups had low levels genetic diversity and were moderately genetically differentiated. Our findings also indicate that T. sutchuenensis suffered two severe bottlenecks around the Last Glaciation Period and Last Glacial Maximum. Among Thuja species, T. sutchuenensis presented the lowest genetic load and hence might have purged deleterious mutations efficiently through purifying selection. However, distribution of fitness effects analysis indicated a high extinction risk for T. sutchuenensis. Multiple lines of evidence identified three management units for T. sutchuenensis. Although T. sutchuenensis possesses a low genetic load, low genetic diversity, suboptimal fitness, and anthropogenic pressures all present an extinction risk for this rare conifer. This might also hold true for many endangered plant species in the mountains all over the world.

Incomplete lineage sorting and local extinction shaped the complex evolutionary history of the Paleogene relict conifer genus, Chamaecyparis (Cupressaceae).

Inferring accurate biogeographic history of plant taxa with an East Asia (EA)-North America (NA) is usually hindered by conflicting phylogenies and a poor fossil record. The current distribution of Chamaecyparis (false cypress; Cupressaceae) with four species in EA, and one each in western and eastern NA, and its relatively rich fossil record, make it an excellent model for studying the EA-NA disjunction. Here we reconstruct phylogenomic relationships within Chamaecyparis using > 1400 homologous nuclear and 61 plastid genes. Our phylogenomic analyses using concatenated and coalescent approaches revealed strong cytonuclear discordance and conflicting topologies between nuclear gene trees. Incomplete lineage sorting (ILS) and hybridization are possible explanations of conflict; however, our coalescent analyses and simulations suggest that ILS is the major contributor to the observed phylogenetic discrepancies. Based on a well-resolved species tree and four fossil calibrations, the crown lineage of Chamaecyparis is estimated to have originated in the upper Cretaceous, followed by diversification events in the early and middle Paleogene. Ancestral area reconstructions suggest that Chamaecyparis had an ancestral range spanning both EA and NA. Fossil records further indicate that this genus is a relict of the "boreotropical" flora, and that local extinctions of European species were caused by global cooling. Overall, our results unravel a complex evolutionary history of a Paleogene relict conifer genus, which may have involved ILS, hybridization and the extinction of local species.

Seeing through the hedge: Phylogenomics of Thuja (Cupressaceae) reveals prominent incomplete lineage sorting and ancient introgression for Tertiary relict flora.

The Eastern Asia (EA) - North America (NA) disjunction is a well-known biogeographic pattern of the Tertiary relict flora; however, few studies have investigated the evolutionary history of this disjunction using a phylogenomic approach. Here, we used 2369 single copy nuclear genes and nearly full plastomes to reconstruct the evolutionary history of the small Tertiary relict genus Thuja, which consists of five disjunctly distributed species. The nuclear species tree strongly supported an EA clade Thuja standishii-Thuja sutchuenensis and a "disjunct clade", where western NA species T. plicata is sister to an EA-eastern NA disjunct Thuja occidentalis-Thuja koraiensis group. Our results suggested that the observed topological discordance among the gene trees as well as the cytonuclear discordance is mainly due to incomplete lineage sorting, probably facilitated by the fast diversification of Thuja around the Early Miocene and the large effective population sizes of ancestral lineages. Furthermore, approximately 20% of the T. sutchuenensis nuclear genome is derived from an unknown ancestral lineage of Thuja, which might explain the close resemblance of its cone morphology to that of an ancient fossil species. Overall, our study demonstrates that single genes may not resolve interspecific relationships for disjunct taxa, and that more reliable results will come from hundreds or thousands of loci, revealing a more complex evolutionary history. This will steadily improve our understanding of their origin and evolution.