Phylotranscriptomic and ecological analyses reveal the evolution and morphological adaptation of Abies

Coniferous forests are under severe threat of the rapid anthropogenic climate warming. Abies (firs), the fourth-largest conifer genus, is a keystone component of the boreal and temperate dark-coniferous forests and harbors a remarkably large number of relict taxa. However, the uncertainty of the phylogenetic and biogeographic history of Abies significantly impedes our prediction of future dynamics and efficient conservation of firs. In this study, using 1,533 nuclear genes generated from transcriptome sequencing and a complete sampling of all widely recognized species, we have successfully reconstructed a robust phylogeny of global firs, in which four clades are strongly supported and all intersectional relationships are resolved, although phylogenetic discordance caused mainly by incomplete lineage sorting and hybridization was detected. Molecular dating and ancestral area reconstruction suggest a Northern Hemisphere high-latitude origin of Abies during the Late Cretaceous, but all extant firs diversified during the Miocene to the Pleistocene, and multiple continental and intercontinental dispersals took place in response to the late Neogene climate cooling and orogenic movements. Notably, four critically endangered firs endemic to subtropical mountains of China, including A. beshanzuensis, A. ziyuanensis, A. fanjingshanensis and A. yuanbaoshanensis from east to west, have different origins and evolutionary histories. Moreover, three hotspots of species richness, including western North America, central Japan, and the Hengduan Mountains, were identified in Abies. Elevation and precipitation, particularly precipitation of the coldest quarter, are the most significant environmental factors driving the global distribution pattern of fir species diversity. Some morphological traits are evolutionarily constrained, and those linked to elevational variation (e.g., purple cone) and cold resistance (e.g., pubescent branch and resinous bud) may have contributed to the diversification of global firs. Our study sheds new light on the spatiotemporal evolution of global firs, which will be of great help to forest management and species conservation in a warming world.

Comparisons of Chloroplast Genome Mutations among 13 Samples of Oil-Tea Camellia from South China

The differences in cpDNA SNPs and InDels of 13 samples from single trees of different species or populations of oil-tea camellia in South China were examined in this study, and phylogenetic trees were reconstructed based on CDSs and non-CDSs of cpDNAs to research the evolutionary relationships among all samples. The SNPs of all samples included all kinds of substitutions, and the frequency of the transition from AT to GC was highest; meanwhile, the frequencies of all kinds of transversions differed among the samples, and the SNPs exhibited polymorphism. The SNPs were distributed in all the different functional regions of cpDNAs, and approximately half of all SNPs in exons led to missense mutations and the gain or loss of termination codons. There were no InDels in the exons of any cpDNA samples, except those retrieved from Camellia gigantocarpa, although this InDel did not lead to a frame shift. The InDels of all cpDNA samples were unevenly distributed in the intergenic region and upstream and downstream of genes. The genes, regions of the same gene, sites and mutation types in the same region related to the distributions of SNPs, and InDels were inconsistent among samples. The 13 samples were divided into 2 clades and 7 or 6 subclades, and the samples of species from the same sections of the Camellia genus did not belong to the same subclades. Meanwhile, the genetic relationship between the samples of Camellia vietnamensis and the undetermined species from Hainan Province or the population of C. gauchowensis in Xuwen was closer than that between C. vietnamensis and the population of C. gauchowensis in Luchuan, and the genetic relationship among C. osmantha, C. vietnamensis and C. gauchowensis was very close. In sum, SNPs and InDels in the different cpDNAs resulted in variable phenotypes among the different species or populations, and they could be developed into molecular markers for studies on species and population identification and phylogenetic relationships. The conclusion from the identification of undetermined species from Hainan Province and the phylogenetic relationships among 13 oil-tea camellia samples based on cpCDS and cpnon-CDS sequences were the same as those from the former report.

Comparison of the Chloroplast Genome Sequences of 13 Oil-Tea Camellia Samples and Identification of an Undetermined Oil-Tea Camellia Species From Hainan Province

The comparison of chloroplast genome (cpDNA) sequences among different plant species is an important source of plant molecular phylogenetic data. In this paper, the cpDNA sequences of 13 different oil-tea camellia samples were compared to identify an undetermined oil-tea camellia species from Hainan Province. The cpDNA of the samples was sequenced and resequenced, and divergence hotspots and simple sequence repeat (SSR) variations were analyzed. Bayesian inference (BI) and maximum-likelihood (ML) phylogenetic trees were constructed based on the full cpDNA sequences. The cpDNA sequences were 156512∼157089 bp in length and had the circular tetrad structure typical of angiosperms. The inverted repeats (IRs) of different species included varying contractions and expansions. The cpDNA sequences of the samples of the undetermined species of oil-tea camellia from Hainan Province and Camellia gauchowensis from Xuwen County were identical. In total, 136 genes were annotated, including 91 protein-coding genes (PCGs), 37 tRNA genes and 8 rRNA genes. The GC content of the cpDNA was 37.3%. The small single-copy (SSC)/IR boundary was rich in variation. Divergence hotspots were mainly located in the intergenic space (IGS) and coding sequences (CDSs), and there were obvious differences in divergence hotspots among species. The same divergence hotspots were found in Camellia vietnamensis, Camellia gauchowensis and the undetermined species of oil-tea camellia from Hainan Province. A total of 191∼198 SSR loci were detected. Most of the SSRs included A or T, and the distribution of SSRs in the cpDNA was uneven. Different species shared common SSRs and exhibited unique SSRs. Based on the full cpDNA sequences, the evolutionary relationships of different species of Camellia were well identified. The thirteen samples were classified into 2 clades and 6 subclades, and the different sections of Camellia clustered on the same branch in 2 clades and 2 subclades. Camellia vietnamensis was more closely related to the undetermined species of oil-tea camellia from Hainan Province and the sample of Camellia gauchowensis from Xuwen County than to the sample of Camellia gauchowensis from Luchuan County. Camellia osmantha was closely related to Camellia gauchowensis and Camellia vietnamensis. In conclusion, the cpDNA of different oil-tea camellia species has a conserved tetrad structure with certain length polymorphisms. SSRs are expected to be developed as "barcodes" or "identity cards" for species identification. SSR variations and other factors result in abundant divergence hotspots in the CDSs and IGS (one non-CDS region), indicating that full cpDNA sequences can be used for the species identification and phylogenetic analysis of Camellia. Accordingly, the undetermined species of oil-tea camellia from Hainan Province is likely Camellia vietnamensis, Camellia vietnamensis and Camellia gauchowensis may be the same species, and additional genetic evidence is needed to determine whether Camellia osmantha is a new independent species. The previous division of related sections of Camellia may need readjustment based on full cpDNA sequences.

Ten Years of Provenance Trials and Application of Multivariate Random Forests Predicted the Most Preferable Seed Source for Silviculture of Abies sachalinensis in Hokkaido, Japan

Research highlights: Using 10-year tree height data obtained after planting from the range-wide provenance trials of Abies sachalinensis, we constructed multivariate random forests (MRF), a machine learning algorithm, with climatic variables. The constructed MRF enabled prediction of the optimum seed source to achieve good performance in terms of height growth at every planting site on a fine scale. Background and objectives: Because forest tree species are adapted to the local environment, local seeds are empirically considered as the best sources for planting. However, in some cases, local seed sources show lower performance in height growth than that showed by non-local seed sources. Tree improvement programs aim to identify seed sources for obtaining high-quality timber products by performing provenance trials. Materials and methods: Range-wide provenance trials for one of the most important silvicultural species, Abies sachalinensis, were established in 1980 at nine transplanting experimental sites. We constructed an MRF to estimate the responses of tree height at 10 years after planting at eight climatic variables at 1 km × 1 km resolution. The model was applied for prediction of tree height throughout Hokkaido Island. Results: Our model showed that four environmental variables were major factors affecting height growth—winter solar radiation, warmth index, maximum snow depth, and spring solar radiation. A tree height prediction map revealed that local seeds showed the best performance except in the southernmost region and several parts of northern regions. Moreover, the map of optimum seed provenance suggested that deployment of distant seed sources can outperform local sources in the southernmost and northern regions. Conclusions: We predicted that local seeds showed optimum growth, whereas non-local seeds had the potential to outperform local seeds in some regions. Several deployment options were proposed to improve tree growth.

Early diversification and permeable species boundaries in the Mediterranean firs

BACKGROUND AND AIMS

Inferring the evolutionary relationships of species and their boundaries is critical in order to understand patterns of diversification and their historical drivers. Despite Abies (Pinaceae) being the second most diverse group of conifers, the evolutionary history of Circum-Mediterranean firs (CMFs) remains under debate.

METHODS

We used restriction site-associated DNA sequencing (RAD-seq) on all proposed CMF taxa to investigate their phylogenetic relationships and taxonomic status.

KEY RESULTS

Based on thousands of genome-wide single nucleotide polymorphisms (SNPs), we present here the first formal test of species delimitation, and the first fully resolved, complete species tree for CMFs. We discovered that all previously recognized taxa in the Mediterranean should be treated as independent species, with the exception of Abies tazaotana and Abies marocana. An unexpectedly early pulse of speciation in the Oligocene-Miocene boundary is here documented for the group, pre-dating previous hypotheses by millions of years, revealing a complex evolutionary history encompassing both ancient and recent gene flow between distant lineages.

CONCLUSIONS

Our phylogenomic results contribute to shed light on conifers' diversification. Our efforts to resolve the CMF phylogenetic relationships help refine their taxonomy and our knowledge of their evolution.

Genetic variation in Tertiary relics: The case of eastern-Mediterranean Abies (Pinaceae)

The eastern‐Mediterranean Abies taxa, which include both widely distributed species and taxa with minuscule ranges, represent a good model to study the impacts of range size and fragmentation on the levels of genetic diversity and differentiation. To assess the patterns of genetic diversity and phylogenetic relationships among eastern‐Mediterranean Abies taxa, genetic variation was assessed by eight nuclear microsatellite loci in 52 populations of Abies taxa with a focus on those distributed in Turkey and the Caucasus. Both at the population and the taxon level, the subspecies or regional populations of Abies nordmanniana s.l. exhibited generally higher allelic richness, private allelic richness, and expected heterozygosity compared with Abies cilicica s.l. Results of both the structure analysis and distance‐based approaches showed a strong differentiation of the two A. cilicica subspecies from the rest as well as from each other, whereas the subspecies of A. nordmanniana were distinct but less differentiated. ABC simulations were run for a set of scenarios of phylogeny and past demographic changes. For A. ×olcayana, the simulation gave a poor support for the hypothesis of being a taxon resulting from a past hybridization, the same is true for Abies equi‐trojani: both they represent evolutionary branches of Abies bornmuelleriana.

The complete plastid genome sequence of Abies koreana (Pinaceae: Abietoideae)

The nucleotide sequence of the chloroplast genome from Abies koreana is the first to have complete genome sequence from genus Abies of family Pinaceae. The circular double-stranded DNA, which consists of 121,373 base pairs (bp), contains a pair of very short inverted repeat regions (IRa and IRb) of 264 bp each, which are separated by a small and large single-copy regions (SSC and LSC) of 54,197 and 66,648 bp, respectively. The genome contents of 114 genes (68 peptide-encoding genes, 35 tRNA genes, four rRNA genes, six open reading frames and one pseudogene) are similar to the chloroplast DNA of other species of Abietoideae. Loss of ndh genes was also identified in the genome of A. koreana like other genomes in the family Pinaceae. Thirteen genes contain one (11 genes) or two (rps12 and ycf3 genes) introns. In phylogenetic analysis, the tree confirms that Abies, Keteleeria and Cedrus are strongly supported as monophyletic. Other inverted repeat sequences located in 42-kb inversion points (1186 bp) include trnS-psaM-ycf12- ψtrnG genes.