Genetic diversity and biogeography of Cunninghamia konishii (Cupressaceae), an island species in Taiwan: a comparison with Cunninghamia lanceolata, a mainland species in China

A taxonomic revision of the genus Angelica (Apiaceae) in Taiwan with a new species A. aliensis

BACKGROUND

Angelica L. sensu lato is a taxonomically complex genus, and many studies have utilized morphological and molecular features to resolve its classification issues. In Taiwan, there are six taxa within Angelica, and their taxonomic treatments have been a subject of controversy. In this study, we conducted a comprehensive analysis incorporating morphological and molecular (cpDNA and nrDNA) characteristics to revise the taxonomic treatments of Angelica in Taiwan.

RESULTS

As a result of our research, we have revised the classification between A. dahurica var. formosana and A. pubescens and merged two varieties of A. morrisonicola into a single taxon. A new taxon, A. aliensis, has been identified and found to share a close relationship with A. tarokoensis. Based on the morphological and molecular characteristics data, it has been determined that the former three taxa should be grouped into the Eurasian Angelica clade, while the remaining four taxa should belong to the littoral Angelica clade. Furthermore, Angelica species in Taiwan distributed at higher altitudes displayed higher genetic diversity, implying that the central mountain range of Taiwan serves as a significant reservoir of plant biodiversity. Genetic drift, such as bottlenecks, has been identified as a potential factor leading to the fixation or reduction of genetic diversity of populations in most Angelica species. We provide key to taxa, synopsis, phenology, and distribution for each taxon of Taiwan.

CONCLUSIONS

Our comprehensive analysis of morphological and molecular features has shed light on the taxonomic complexities within Angelica in Taiwan, resolving taxonomic issues and providing valuable insights into the phylogenetic relationships of Angelica in Taiwan.

Effect of geographic isolation on genetic variation and population structure of Euphrasia nankotaizanensis, a threatened endemic alpine herb in Taiwan

Euphrasia nankotaizanensis (Orobanchaceae) is a rare alpine herb that is endemic to Taiwan. Only four small populations remain in Xue, Nanhu, and Cilai Mountains of Taiwan. The distribution of alpine herbs is severely threatened by climate change, which influences genetic variation and population structure. In this study, we investigated the effects of the natural isolation of alpine habitats on the genetic diversity and geographic structure of populations of E. nankotaizanensis using chloroplast (cp) and nuclear DNA (nrDNA) markers. We found lower levels of genetic diversity in E. nankotaizanensis than in other alpine plants and little to no genetic variation within populations, which could be mainly attributed to the small population size and genetic drift. Only one nrDNA haplotype was present in each population. The lack of monophyly of the four populations in cpDNA probably resulted from lineage sorting or occasional long-distance seed dispersal. Phylogeographic analysis suggested that Nanhu Mountain was probably a refugium over the glacial maxima, agreeing with the potential refugia in central Taiwan. The STRUCTURE and AMOVA analyses revealed significant genetic differentiation in nrDNA among the mountains, which resulted from geographical isolation among these mountains. Estimates of the effective population size (Ne) and demography reflected lower Ne values and a recent population decline, probably implying a greater extinction risk for E. nankotaizanensis. We observed genetic depletion and considerable genetic differentiation among mountain populations, which should be considered in future conservation efforts for this species. In addition, this study provides important insights into the long-term potential of alpine herbs in Taiwan, which are useful for a better prediction of their responses to future climate change.

Genetic diversity and structure of the 4th cycle breeding population of Chinese fir (Cunninghamia lanceolata (lamb.) hook)

Studying population genetic structure and diversity is crucial for the marker-assisted selection and breeding of coniferous tree species. In this study, using RAD-seq technology, we developed 343,644 high-quality single nucleotide polymorphism (SNP) markers to resolve the genetic diversity and population genetic structure of 233 Chinese fir selected individuals from the 4^th cycle breeding program, representing different breeding generations and provenances. The genetic diversity of the 4^th cycle breeding population was high with nucleotide diversity (P_i ) of 0.003, and H_o and H_e of 0.215 and 0.233, respectively, indicating that the breeding population has a broad genetic base. The genetic differentiation level between the different breeding generations and different provenances was low (F_st < 0.05), with population structure analysis results dividing the 233 individuals into four subgroups. Each subgroup has a mixed branch with interpenetration and weak population structure, which might be related to breeding rather than provenance, with aggregation from the same source only being in the local branches. Our results provide a reference for further research on the marker-assisted selective breeding of Chinese fir and other coniferous trees.

Characterization and Prediction of Mechanical and Chemical Properties of Luanta Fir Wood with Vacuum Hydrothermal Treatment

Since the chemical composition of wood is closely related to its mechanical properties, chemical analysis techniques such as near-infrared (NIR) spectroscopy provide a reasonable non-destructive method for predicting wood strength. In this study, we used NIR spectra with principal component analysis (PCA) to reveal that vacuum hydrothermal (VH) treatment causes degradation of hemicellulose as well as the amorphous region of cellulose, resulting in lower hydroxyl and acetyl group content. These processes increase the crystallinity of the luanta fir wood (Cunninghamia konishii Hayata), which, in turn, effectively increases its compressive strength (σc,max), hardness, and modulus of elasticity (MOE). The PCA results also revealed that the primary factors affecting these properties are the hemicellulose content, hydroxyl groups in the cellulose amorphous region, the wood moisture content, and the relative lignin content. Moreover, the ratios of performance deviation (RPDs) for the σc,max, shear strength (σs,max), hardness, and modulus of rupture (MOR) models were 1.49, 1.24, 1.13, and 2.39, indicating that these models can be used for wood grading (1.0 < RPD < 2.5). Accordingly, NIR can serve as a useful tool for predicting the mechanical properties of VH-treated wood.

Differentiation patterns of emperor moths (Lepidoptera: Saturniidae: Saturniinae) of a continental island: divergent evolutionary history driven by Pleistocene glaciations

Background

On the basis of molecular dating, Pleistocene glaciations have been proposed as the major driving force of biota speciation in the Palearctic and the pre-Quaternary origin of Amazonian taxa. However, the major driving factors in East Asia remain unclear. All 16 saturniine species inhabiting Taiwan with congeners of populations, subspecies, or species in East Asia constitute research objects for addressing the mode of speciation because of the repeated formation and disappearance of a landbridge from the Asian mainland to Taiwan during glacial cycles.

Methods

The genetic divergences of mitochondrial cytochrome c oxidase subunit I (COI) and 16S rDNA and the nuclear 28S rDNA of the saturniine species from Taiwan and the Asian mainland were assessed to determine the monophyly of each genus and species of Saturniinae. Moreover, 519 saturniine COI sequences of 114 taxa from adjacent East and Southeast Asian populations and closely related species were retrieved from GenBank and analyzed. The differentiation timing and possible origination of the insular saturniines were elucidated based on phylogenetic relationships, haplotype networks, and lineage calibrations.

Results

Approximately 90% of intraspecific COI divergence was <2%; all divergences exceeding 2% originated from comparisons between allopatric populations or subspecies. Relationship analyses revealed that multiple introductions likely occurred in insular saturniines and that some East Asian saturniines were paraphyletic as deduced by analyzing endemic insular species. Calibration dating revealed that Taiwanese endemic saturniines split from sibling Asian species 0.2-2.7 million years ago (Mya), whereas subspecific-level and population-level splitting events occurred 0.1-1.7 Mya and 0.2-1.2 Mya, respectively. Moreover, phylogenetic patterns combined with geographical distributions revealed that hill-distributed Taiwanese saturniines are closely related to those from southern China and Southeast Asia, whereas saturniines inhabiting altitudes higher than 1,500 m in Taiwan have siblings distributed in temperate Northeast Asia.

Discussion

The Global DNA Barcoding Initiative was successfully applied to study the population genetic structure in species. Most Formosan saturniines are distinct and monophyletic, reflecting the vicariant barrier of the Taiwan Strait; Pleistocene glacial cycles provided opportunities for insular saturniines to experience repeated isolation from and secondary contact with the continental mainland. Each insular saturniine may have evolved with a unique differentiation timing pattern that possibly emerged in the Early, Middle, or Late Pleistocene with these patterns differing from the consistent pattern that occurred in the temperate Palearctic and tropical Amazonian regions. Moreover, multiple migrations or artificial genetic admixtures may have also occurred, as suggested by the coexistence of two divergent lineages in a few Taiwanese saturniines.

Genome skimming and microsatellite analysis reveal contrasting patterns of genetic diversity in a rare sandhill endemic (Erysimum teretifolium, Brassicaceae)

Barriers between islands often inhibit gene flow creating patterns of isolation by distance. In island species, the majority of genetic diversity should be distributed among isolated populations. However, a self-incompatible mating system leads to higher genetic variation within populations and very little between-population subdivision. We examine these two contrasting predictions in Erysimum teretifolium, a rare self-incompatible plant endemic to island-like sandhill habitats in Santa Cruz County, California. We used genome skimming and nuclear microsatellites to assess the distribution of genetic diversity within and among eight of the 13 remaining populations. Phylogenetic analyses of the chloroplast genomes revealed a deep separation of three of the eight populations. The nuclear ribosomal DNA cistron showed no genetic subdivision. Nuclear microsatellites suggest 83% of genetic variation resides within populations. Despite this, 18 of 28 between-population comparisons exhibited significant population structure (mean F_ST = 0.153). No isolation by distance existed among all populations, however when one outlier population was removed from the analysis due to uncertain provenance, significant isolation by distance emerged (r² = 0.5611, p = 0.005). Population census size did not correlate with allelic richness as predicted on islands. Bayesian population assignment detected six genetic groupings with substantial admixture. Unique genetic clusters were concentrated at the periphery of the species’ range. Since the overall distribution of nuclear genetic diversity reflects E. tereifolium’s self-incompatible mating system, the vast majority of genetic variation could be sampled within any individual population. Yet, the chloroplast genome results suggest a deep split and some of the nuclear microsatellite analyses indicate some island-like patterns of genetic diversity. Restoration efforts intending to maximize genetic variation should include representatives from both lineages of the chloroplast genome and, for maximum nuclear genetic diversity, should include representatives of the smaller, peripheral populations.

Genome survey of Chinese fir (Cunninghamia lanceolata): Identification of genomic SSRs and demonstration of their utility in genetic diversity analysis

Chinese fir (Cunninghamia lanceolata) is an important coniferous species that accounts for 20–30% of the total commercial timber production in China. Though traditional breeding of Chinese fir has achieved remarkable success, molecular-assisted breeding has made little progress due to limited availability of genomic information. In this study, a survey of Chinese fir genome was performed using the Illumina HiSeq Xten sequencing platform. K-mer analysis indicated that Chinese fir has a large genome of approximately 11.6 Gb with 74.89% repetitive elements and is highly heterozygous. Meanwhile, its genome size was estimated to be 13.2 Gb using flow cytometry. A total of 778.02 Gb clean reads were assembled into 10,982,272 scaffolds with an N50 of 1.57 kb. In total, 362,193 SSR loci were detected with a frequency of 13.18 kb. Dinucleotide repeats were the most abundant (up to 73.6% of the total SSRs), followed by trinucleotide and tetranucleotide repeats. Forty-six polymorphic pairs were developed, and 298 alleles were successfully amplified from 199 Chinese fir clones. The average PIC value was 0.53, indicating that the identified genomic SSR (gSSR) markers have a high degree of polymorphism. In addition, these breeding resources were divided into three groups, and a limited gene flow existed among these inferred groups.

Testing the Effect of Mountain Ranges as a Physical Barrier to Current Gene Flow and Environmentally Dependent Adaptive Divergence in Cunninghamia konishii (Cupressaceae)

Populations can be genetically isolated by differences in their ecology or environment that hampered efficient migration, or they may be isolated solely by geographic distance. Moreover, mountain ranges across a species’ distribution area might have acted as barriers to gene flow. Genetic variation was quantified using amplified fragment length polymorphism (AFLP) and 13 selective amplification primer combinations used generated a total of 482 fragments. Here, we tested the barrier effects of mountains on gene flow and environmentally dependent local adaptation of Cunninghamia konishii occur in Taiwan. A pattern of genetic isolation by distance was not found and variation partitioning revealed that environment explained a relatively larger proportion of genetic variation than geography. The effect of mountains as barriers to genetic exchange, despite low population differentiation indicating a high rate of gene flow, was found within the distribution range of C. konishii. Twelve AFLP loci were identified as potential selective outliers using genome-scan methods (BAYESCAN and DFDIST) and strongly associated with environmental variables using regression approaches (LFMM, Samβada, and rstanarm) demonstrating adaptive divergence underlying local adaptation. Annual mean temperature, annual precipitation, and slope could be the most important environmental factors causally associated with adaptive genetic variation in C. konishii. The study revealed the existence of physical barriers to current gene flow and environmentally dependent adaptive divergence, and a significant proportion of the rate of gene flow may represent a reflection of demographic history.

Nuclear and chloroplast DNA phylogeography suggests an Early Miocene southward expansion of Lithocarpus (Fagaceae) on the Asian continent and islands

BACKGROUND

Most genera of Fagaceae are thought to have originated in the temperate regions except for the genus Lithocarpus, the stone oaks. Lithocarpus is distributed in subtropical and tropical Asia, and its ancestral population is hypothesized to be distributed in tropical regions in Borneo and Indochina. Borneo and the nearby islands (the Greater Sunda Islands) were connected to the Malay Peninsula and Indochina prior to the Pliocene epoch and formed the former Sundaland continent. The Southeast Asian Lithocarpus, is thought to have dispersed between continental Asia and the present Sundaland. The drastic climate changes during the Pliocene and Pleistocene epochs which caused periodic sea-level changes is often used to explain the cause of its diversity. The aim of this study was to establish phylogenetic relationships by analyzing nuclear (nrDNA) and chloroplast (cpDNA) DNA in order to describe and analyze the origin, causes of diversification and historical biogeography of Lithocarpus.

RESULTS

Phylogeny reconstructed through the multiple-species coalescent method with nrDNA and cpDNA revealed that the continental-Asian taxa were clustered at the basal lineages. The derived lineages of tropical Lithocarpus, with the inference of a subtropical ancestral state, imply a southward migration in the Early Miocene period with subsequent in situ diversification in the Greater Sunda Islands. The gradual decrease in temperature since the Middle Miocene period is proposed as a cause of the increase in the net diversification rate.

CONCLUSIONS

The historical ancestral origin of Lithocarpus has been suggested to be mainland Asia. Southward migration in the Early Miocene period with subsequent in situ diversification could explain the current diversity of stone oaks in Southeast Asia. This study also considered the multiple origins of stone oaks currently indigenous to the subtropical islands offshore and near mainland China. Our results provide phylogenetic evidence for a subtropical origin of Asian stone oaks and reveal the process of diversification and how it fits into the timeline of major geologic and climatic events rather than local, episodic, rate-shifting events.

The Contribution of Neutral and Environmentally Dependent Processes in Driving Population and Lineage Divergence in Taiwania (Taiwania cryptomerioides)

The question of what determines divergence both between and within species has been the central topic in evolutionary biology. Neutral drift and environmentally dependent divergence are predicted to play roles in driving population and lineage divergence. However, neutral drift may preclude adaptation if the rate of gene flow between populations is high. Here, we sampled populations of three Taiwania (Taiwania cryptomerioides) lineages occurring in Taiwan, the mainland of China (Yunnan-Myanmar border), and northern Vietnam, and tested the relative strength of neutral drift and divergent selection in shaping divergence of those populations and lineages. We quantified genetic and epigenetic variation, respectively, using amplified fragment length polymorphism (AFLP) and methylation-sensitive amplification polymorphism (MSAP). Analysis of 1413 AFLP and 462 MSAP loci using frequency-based genome scan methods and generalized linear models (GLMs) found no potential selective outliers when only Taiwanese populations were examined, suggesting that neutral drift was the predominant evolutionary process driving differentiation between those populations. However, environmentally associated divergence was found when lineages were compared. Thirty-two potential selective outliers were identified based on genome scans and their associations with environmental variables were tested with GLMs, generalized linear mixed effect models (GLMMs), and model selection with a model averaging approach. Ten loci (six AFLP and four MSAP) were found to be strongly associated with environmental variables, particularly monthly temperature variation and normalized difference vegetation index (NDVI) using model selection and a model averaging approach. Because only a small portion of genetic and epigenetic loci were found to be potential selective outliers, neutral evolutionary process might also have played crucial roles in driving lineage divergence, particularly between geographically and genetically isolated island and mainland Asia lineages. Nevertheless, the vast amount of neutral drift causing genetic and epigenetic variations might have the potential for adaptation to future climate changes. These could be important for the survival of Taiwania in different geographic areas.

Adaptive Genetic Divergence Despite Significant Isolation-by-Distance in Populations of Taiwan Cow-Tail Fir (Keteleeria davidiana var. formosana)

Double digest restriction site-associated DNA sequencing (ddRADseq) is a tool for delivering genome-wide single nucleotide polymorphism (SNP) markers for non-model organisms useful in resolving fine-scale population structure and detecting signatures of selection. This study performs population genetic analysis, based on ddRADseq data, of a coniferous species, Keteleeria davidiana var. formosana, disjunctly distributed in northern and southern Taiwan, for investigation of population adaptive divergence in response to environmental heterogeneity. A total of 13,914 SNPs were detected and used to assess genetic diversity, FST outlier detection, population genetic structure, and individual assignments of five populations (62 individuals) of K. davidiana var. formosana. Principal component analysis (PCA), individual assignments, and the neighbor-joining tree were successful in differentiating individuals between northern and southern populations of K. davidiana var. formosana, but apparent gene flow between the southern DW30 population and northern populations was also revealed. Fifteen of 23 highly differentiated SNPs identified were found to be strongly associated with environmental variables, suggesting isolation-by-environment (IBE). However, multiple matrix regression with randomization analysis revealed strong IBE as well as significant isolation-by-distance. Environmental impacts on divergence were found between populations of the North and South regions and also between the two southern neighboring populations. BLASTN annotation of the sequences flanking outlier SNPs gave significant hits for three of 23 markers that might have biological relevance to mitochondrial homeostasis involved in the survival of locally adapted lineages. Species delimitation between K. davidiana var. formosana and its ancestor, K. davidiana, was also examined (72 individuals). This study has produced highly informative population genomic data for the understanding of population attributes, such as diversity, connectivity, and adaptive divergence associated with large- and small-scale environmental heterogeneity in K. davidiana var. formosana.

Postglacial range expansion and the role of ecological factors in driving adaptive evolution of Musa basjoo var. formosana

Genetic variation evolves during postglacial range expansion of a species and is important for adapting to varied environmental conditions. It is crucial for the future survival of a species. We investigate the nuclear DNA sequence variation to provide evidence of postglacial range expansion of Musa basjoo var. formosana, a wild banana species, and test for adaptive evolution of amplified fragment length polymorphic (AFLP) loci underlying local adaptation in association with environmental variables. Postglacial range expansion was suggested by phylogeographical analyses based on sequence variation of the second intron of copper zinc superoxide dismutase 2 gene. Two glacial refugia were inferred by the average F_ST parameter (mean F_ST of a population against the remaining populations). Using variation partitioning by redundancy analysis, we found a significant amount of explained AFLP variation attributed to environmental and spatially-structured environmental effects. By combining genome scan methods and multiple univariate logistic regression, four AFLP loci were found to be strongly associated with environmental variables, including temperature, precipitation, soil moisture, wet days, and surface coverage activity representing vegetation greenness. These environmental variables may have played various roles as ecological drivers for adaptive evolution of M. basjoo var. formosana during range expansion after the last glacial maximum.

Genetic diversity and variation of Chinese fir from Fujian province and Taiwan, China, based on ISSR markers

Genetic diversity and variation among 11 populations of Chinese fir from Fujian province and Taiwan were assessed using inter-simple sequence repeat (ISSR) markers to reveal the evolutionary relationship in their distribution range in this report. Analysis of genetic parameters of the different populations showed that populations in Fujian province exhibited a greater level of genetic diversity than did the populations in Taiwan. Compared to Taiwan populations, significant limited gene flow were observed among Fujian populations. An UPGMA cluster analysis showed that the most individuals of Taiwan populations formed a single cluster, whereas 6 discrete clusters were formed by each population from Fujian. All populations were divided into 3 main groups and that all 5 populations from Taiwan were gathered into a subgroup combined with 2 populations, Dehua and Liancheng, formed one of the 3 main groups, which indicated relative stronger relatedness. It is supported by a genetic structure analysis. All those results are suggesting different levels of genetic diversity and variation of Chinese fir between Fujian and Taiwan, and indicating different patterns of evolutionary process and local environmental adaption.

New insights into the phylogeny of Burasaieae (Menispermaceae) with the recognition of a new genus and emphasis on the southern Taiwanese and mainland Chinese disjunction

Taiwan is a continental island lying at the boundary between the Eurasian and the Philippine tectonic plates and possesses high biodiversity. Southern Taiwan, viz. Hengchun Peninsula, is notably floristically different from northern Taiwan. The floristic origin and relationships of the Hengchun Peninsula have been rarely investigated in a phylogenetic context. In this study, data from six plastid and nuclear sequences were used to reconstruct phylogenetic relationships within Burasaieae (Menispermaceae), which mainly inhabits tropical rainforests. The tree-based comparisons indicate that the position of Tinospora sensu stricto conflicts significantly between the cpDNA and ITS trees. However, alternative hypothesis tests from the ITS data did not reject the result of the cpDNA data, which suggests that tree-based comparisons might sometimes generate an artificial incongruence, especially when markers with high homoplasy are used. Based on the combined cpDNA and ITS data, we present an inter-generic phylogenetic framework for Burasaieae. Sampled species of Tinospora are placed in three different clades, including Tinospora dentata from southern Taiwan and T. sagittata from mainland China in an unresolved position alongside six lineages of Burasaieae. By integrating lines of evidence from molecular phylogeny, divergence times, and morphology, we recognize the three Tinospora clades as three different genera, including Tinospora sensu stricto, a new genus (Paratinospora) for T. dentata and T. sagittata, and Hyalosepalum resurrected. Tinospora dentata, now endemic to the Hengchun Peninsula, originated from the Late Eocene (ca. 39Ma), greatly predating the formation of Taiwan. Our study suggests that the flora of the Hengchun Peninsula contains some ancient components that might have migrated from mainland China.

Genomic scanning using AFLP to detect loci under selection in the moss Funaria hygrometrica along a climate gradient in the Sierra Nevada Mountains, Spain

The common cord moss Funaria hygrometrica has a worldwide distribution and thrives in a wide variety of environments. Here, we studied the genetic diversity in F. hygrometrica along an abiotic gradient in the Mediterranean high mountain of Sierra Nevada (Spain) using a genome scan method. Eighty-four samples from 17 locations from 24 to 2700 m were fingerprinted based on their amplified fragment length polymorphism (AFLP) banding pattern. Using PCA and Bayesian inference we found that the genetic diversity was structured in three or four clusters, respectively. Using a genome scan method we identified 13 outlier loci, which showed a signature of positive selection. Partial Mantel tests were performed between the Euclidean distance matrices of geographic and climatic variables, versus the pair-wise genetic distance of the AFLP dataset and AFLP-positive outliers dataset. AFLP-positive outlier data were significantly correlated with the gradient of the climatic variables, suggesting adaptive variation among populations of F. hygrometrica along the Sierra Nevada Mountains. We highlight the additional analyses necessary to identify the nature of these loci, and their biological role in the adaptation process.

Differentiation in stag beetles, Neolucanus swinhoei complex (Coleoptera: Lucanidae): four major lineages caused by periodical Pleistocene glaciations and separation by a mountain range

Taxonomic debates on Neolucanus swinhoei complex consisting of N. swinhoei, N. doro doro, N. doro horaguchii, and N. euganiae, distributed exclusively in Taiwan, have been ongoing for several decades because of their overlapping morphological characters. To clarify their taxonomic status and phylogeographical history, we analyzed nine morphological characteristics and four molecular amplicons. Phylogenetic inferences based on COI+16S rDNA+wingless showed one eastern and three western lineages, with the latter consisting of one low-hill and two montane lineages. Intermingled DNA sequences from different populations within each lineage, many low FST values, and a high variance component between lineages indicate the possibility of gene flow among populations. However, positive relationships were observed between the genetic divergences of 16S rDNA and its FST values with geographic distance. A divergence estimation based on COI+16S revealed that these beetles might have originated from Asian mainland and differentiated into western and eastern lineages ca. 1Mya, with the differentiation of the western lineages occurring approximately 0.50-0.75Mya. Isolation by mountain ranges and limited flying capability of these beetles as well as populations retreat to and expansion from refugia in response to glaciation cycles have resulted in the current distribution of N. swinhoei complex. Although most morphological characters are variable and undistinguishable, multi-dimensional scaling analysis based on measurable characteristics could recognize hill N. swinhoei as a cluster distinct from the others. However, based on the realities of genetic admixture, shared phylogeographical history and overlapping characteristics, all of these stag beetles should be regarded as Neolucanus swinhoei Bates, 1866.

Divergent selection and local adaptation in disjunct populations of an endangered conifer, Keteleeria davidiana var. formosana (Pinaceae)

The present study investigated the genetic diversity, population structure, F_ST outliers, and extent and pattern of linkage disequilibrium in five populations of Keteleeria davidiana var. formosana, which is listed as a critically endangered species by the Council of Agriculture, Taiwan. Twelve amplified fragment length polymorphism primer pairs generated a total of 465 markers, of which 83.74% on average were polymorphic across populations, with a mean Nei’s genetic diversity of 0.233 and a low level of genetic differentiation (approximately 6%) based on the total dataset. Linkage disequilibrium and HICKORY analyses suggested recent population bottlenecks and inbreeding in K. davidiana var. formosana. Both STRUCTURE and BAPS observed extensive admixture of individual genotypes among populations based on the total dataset in various clustering scenarios, which probably resulted from incomplete lineage sorting of ancestral variation rather than a high rate of recent gene flow. Our results based on outlier analysis revealed generally high levels of genetic differentiation and suggest that divergent selection arising from environmental variation has been driven by differences in temperature, precipitation, and humidity. Identification of ecologically associated outliers among environmentally disparate populations further support divergent selection and potential local adaptation.

Genetic population structure of the alpine species Rhododendron pseudochrysanthum sensu lato (Ericaceae) inferred from chloroplast and nuclear DNA

Background

A complex of incipient species with different degrees of morphological or ecological differentiation provides an ideal model for studying species divergence. We examined the phylogeography and the evolutionary history of the Rhododendron pseudochrysanthum s. l.

Results

Systematic inconsistency was detected between gene genealogies of the cpDNA and nrDNA. Rooted at R. hyperythrum and R. formosana, both trees lacked reciprocal monophyly for all members of the complex. For R. pseudochrysanthum s.l., the spatial distribution of the cpDNA had a noteworthy pattern showing high genetic differentiation (F_ST = 0.56-0.72) between populations in the Yushan Mountain Range and populations of the other mountain ranges.

Conclusion

Both incomplete lineage sorting and interspecific hybridization/introgression may have contributed to the lack of monophyly among R. hyperythrum, R. formosana and R. pseudochrysanthum s.l. Independent colonizations, plus low capabilities of seed dispersal in current environments, may have resulted in the genetic differentiation between populations of different mountain ranges. At the population level, the populations of Central, and Sheishan Mountains may have undergone postglacial demographic expansion, while populations of the Yushan Mountain Range are likely to have remained stable ever since the colonization. In contrast, the single population of the Alishan Mountain Range with a fixed cpDNA haplotype may have experienced bottleneck/founder's events.

Historical spatial range expansion and a very recent bottleneck of Cinnamomum kanehirae Hay. (Lauraceae) in Taiwan inferred from nuclear genes

Background

Species in the varied geographic topology of Taiwan underwent obvious demographic changes during glacial periods. Cinnamomum kanehirae has been exploited for timber and to obtain medicinal fungi for the past 100 years. Understanding anthropogenic factors influencing the demography of this species after the last glacial maximum (LGM) is critically important for the conservation of this species.

Results

Populations of C. kanehirae were classified into four geographic regions: northwestern (NW), west-central (WC), southwestern (SW), and southeastern (SE). In total, 113 individuals from 19 localities were sampled, and variations in the chalcone synthase gene (Chs) intron and leafy (Lfy) intron-2 sequences of nuclear DNA were examined in order to assess phylogeographic patterns, the timescales of demographic and evolutionary events, and recent anthropogenic effects. In total, 210 Chs and 170 Lfy sequences, which respectively constituted 36 and 35 haplotypes, were used for the analyses. Estimates of the migration rate (M) through time revealed a pattern of frequent gene flow during previous and the present interglacials. The isolation-by-distance test showed that there generally was no significant correlation between genetic and geographic distances. The level of among-region genetic differentiation was significant when comparing eastern to western populations. However, no significant among-region genetic differentiation was found in comparisons among the four geographic regions. Moreover, essentially no genetic structuring was found for the three regions west of the CMR. A fit of spatial range expansion was found for pooled and regional samples according to the non-significant values of the sum of squared deviations. Using the Bayesian skyline plot (BSP) method, a recent bottleneck after the LGM expansion was detected in both regional and pooled samples.

Conclusions

Common haplotype distributions among geographic regions and the relatively shallow genetic structuring displayed are the result of historical gene flows. Southward dispersals in an earlier time frame from the NW region and in a later time frame from the SE region were inferred. The BSP analysis suggested a postglacial expansion event. Recent trends, however, refer to a bottleneck due to human interventions observed for both pooled and regional C. kanehirae samples.

Population genetics of Astragalus bibullatus (Fabaceae) using AFLPs

Astragalus bibullatus (Fabaceae) is an endangered plant species endemic to limestone cedar glades in Tennessee. Data from 134 amplified fragment length polymorphism (AFLP) fragments indicate that populations of this species are genetically very similar, with genetic identity values ranging from 0.976 to 0.991. Approximately 10% of the species' total genetic variation is due to differences among populations (F(ST) = 0.102), and a principle coordinate analysis based on genetic distance among individuals revealed considerable overlap between populations. Averaging across populations, 38.4% of the AFLP markers were polymorphic, and the mean expected heterozygosity was 0.120. These estimates are higher than what has previously been reported for this species based on a survey of allozyme variation. Despite the overall similarity of populations of this species, a model-based clustering approach revealed the presence of 2 (possibly 3) genetically distinct subgroups. The results of this study highlight the utility of DNA-based markers for conservation genetic studies in genetically depauperate species and reveal that detectable levels of genetic substructuring may be present even in relatively undifferentiated species.

Molecular differentiation and diversity of Forcipomyia taiwana (Diptera: Ceratopogonidae) based on the mitochondrial cytochrome oxidase II sequence

Forcipomyia taiwana (Shiraki), a biting midge, is one of the most annoying blood-sucking pests in Taiwan. In this study, partial DNA sequences of cytochrome c oxidase II from 113 individuals collected from 11 locations around the island were analyzed to delineate the differentiation pattern and possible dispersal processes of F. taiwana in Taiwan. The uncorrected nucleotide divergences, composed of mostly transition substitutions, were high (up to 2.7%) among the samples. Average comparable variations (approximately equal to 0.7%) were found within and between populations. Phylogenetic analysis suggested that several distinct lineages exist and some can be found simultaneously in some populations. A relationship between sequence divergences among populations and their relative geographical distances was observed. Moreover, haplotype diversity was high in all populations, and low to middle levels (Fst = 0.004-0.288) of genetic differentiation were found among populations. Linearized calibration from sequence divergences and phylogenetic analysis showed that different ancestral lineages of F. taiwana possibly emerged as early as 0.6 million years ago. Taken together, genetic exchanges among these divergently ancestral lineages, likely caused by recent artificial events, have possibly led to the similarly diversified compositions of F. taiwana populations all around Taiwan nowadays.

Uniform genetic diversity, low differentiation, and neutral evolution characterize contemporary refuge populations of Taiwan fir (Abies kawakamii, Pinaceae)

Based on fossil pollen, the distribution range of Taiwan fir [Abies kawakamii (Hay.) Ito] (Pinaceae) is smaller than it was 50 000 years ago. To characterize the present refuge populations of A. kawakamii, which survive only in subalpine forests in Taiwan, we surveyed nuclear genes and chloroplast intergenic spacers to assess the genetic diversity of Taiwan fir. Populations maintain high genetic diversity and contain similar numbers of haplotypes for the GapC (cytosolic glyceraldehyde 3-phosphate dehydrogenase) fragment. Haplotypes for GapC are generally widespread, and population-specific haplotypes accounted for 2.5% of the total. Differentiation among populations is very low (G(ST) = 0.01). Only three haplotypes were detected for the cpDNA marker, and every population had one or two haplotypes. In a neutrality test, the variation in nucleotides did not deviate from that expected with neutral evolution for either marker. A retreat route to higher elevations was not evident from either the GapC or cpDNA markers. Hsuehshan was the site of the most divergent population in Taiwan. We concluded that uniform genetic diversity, low differentiation, low numbers of population-specific haplotypes, and neutral evolution characterize contemporary refuge populations of Taiwan fir.

Allozyme variation of populations of Castanopsis carlesii (Fagaceae) revealing the diversity centres and areas of the greatest divergence in Taiwan

BACKGROUND AND AIMS

The genetic variation and divergence estimated by allozyme analysis were used to reveal the evolutionary history of Castanopsis carlesii in Taiwan. Two major questions were discussed concerning evolutionary issues: where are the diversity centres, and where are the most genetically divergent sites in Taiwan?

METHODS

Twenty-two populations of C. carlesii were sampled throughout Taiwan. Starch gel electrophoresis was used to assay allozyme variation. Genetic parameters and mean FST values of each population were analysed using the BIOSYS-2 program. Mean F(ST) values of each population against the remaining populations, considered as genetic divergence, were estimated using the FSTAT program.

KEY RESULTS

Average values of genetic parameters describing the within-population variation, the average number of alleles per locus (A=2.5), the effective number of alleles per locus (Ae=1.38), the allelic richness (Ar=2.38), the percentage of polymorphic loci (P=69%), and the expected heterozygosity (He=0.270) were estimated. High levels of genetic diversity were found for C. carlesii compared with other local plant species. Genetic differentiation between populations was generally low.

CONCLUSIONS

From the data of expected heterozygosity, one major diversity centre was situated in central Taiwan corroborating previous reports for other plant species. According to the mean FST value of each population, the most divergent populations were situated in two places. One includes populations located in north central Taiwan between 24.80 degrees N and 24.20 degrees N. The other is located in south-eastern Taiwan between 22.40 degrees N and 23.10 degrees N. These two regions are approximately convergent with the most divergent locations determined for several other plant species using chloroplast DNA markers published previously. An important finding obtained from this study is that unordered markers like allozymes can be used to infer past population histories as well as chloroplast DNA markers do.

Extensive gene flow blurs phylogeographic but not phylogenetic signal in Olea europaea L

Genetic structure and evolutionary patterns of the wild olive tree (Olea europaea L.) were investigated with AFLP fingerprinting data at three geographic levels: (a) phylogenetic relationships of the six currently recognized subspecies in Eurasia and Africa; (b) lineage identification in subsp. europaea of the Mediterranean basin; and (c) phylogeography in the western Mediterranean. Two statistical approaches (Bayesian inference and analysis of molecular variance) were used to analyse the AFLP fingerprints. To determine the congruency and transferability of results across studies previous RAPD and ISSR data were analysed in a similar manner. Comparisons proved that qualitative results were mostly congruent but quantitative values differed, depending on the method of analysis. Neighbour-Joining analysis of AFLP phenotypes supported current classification of subspecies. At a Mediterranean scale no clear cut phylogeographic pattern was recovered, likely due to extensive gene flow between populations of subsp. europaea. Gene flow estimates calculated with conventional F-statistics showed that reproductive barriers separated neither populations nor lineages of O. europaea. Genetic divergence between eastern and western parts of the Mediterranean basin was observed only when geographical and population information were incorporated into the analyses through hierarchical analysis of molecular variance (AMOVA). Within the western Mediterranean, the highest genetic diversity was found in two regions: on both sides of the Strait of Gibraltar and in the Balearic archipelago. Additionally, long-lasting isolation of the northern-most populations of the Iberian Peninsula appeared to be responsible for a significant divergence.

Potential refugia in Taiwan revealed by the phylogeographical study of Castanopsis carlesii Hayata (Fagaceae)

In this study, we examined spatial patterns of chloroplast DNA (cpDNA) variation in a total of 30 populations of Castanopsis carlesii Hayata (Fagaceae), a subtropical and temperate tree species, including 201 individuals sampled throughout Taiwan. By sequencing two cpDNA fragments using universal primers (the trnL intron and the trnV-trnM intergenic spacer), we found a total of 1663 bp and 21 polymorphic sites. These gave rise to a total of 28 cpDNA haplotypes. The level of differentiation among the populations studied was relatively high (GST = 0.723). Two ancestral haplotypes are widely distributed. The Central Mountain Ridge (CMR) of Taiwan represents an insurmountable barrier to the east-west gene flow of C. carlesii. Among the populations studied, three separated populations, at Lienhuachih, Fushan and Lichia, have high nucleotide diversity. Estimates of NST-GST for populations on both sides of the CMR indicate that no phylogeographical structure exists. According to the genealogical tree, number of rare haplotype and population genetic divergence, this study suggests that two potential refugia existed during the last glaciation: the first refugium was located in a region to the north of Hsuehshan Range (HR) and west of the CMR; the second refugium was located in south, especially southeastern Taiwan. In fact, the second refugium happens to be the same as that reported for Quercus glauca. A 'star-like' genealogy is characteristic when all haplotypes rapidly coalesce and is a general outcome of population expansion. The neutrality test and mismatch distribution also suggest demographic expansion recovering from a bottleneck.