Tzen Yuh Chiang

Phylogeography of the Qinghai-Tibetan Plateau endemic Juniperus przewalskii (Cupressaceae) inferred from chloroplast DNA sequence variation

The vegetation of the northeast Qinghai‐Tibetan Plateau is dominated by alpine meadow and desert‐steppe with sparse forests scattered within it. To obtain a better understanding of the phylogeography of one constituent species of the forests in this region, we examined chloroplast trnT‐trnF and trnS‐trnG sequence variation within Juniperus przewalskii, a key endemic tree species. Sequence data were obtained from 392 trees in 20 populations covering the entire distribution range of the species. Six cpDNA haplotypes were identified. Significant population subdivision was detected (GST = 0.772, NST = 0.834), suggesting low levels of recurrent gene flow among populations and significant phylogeographic structure (NST > GST, P < 0.05). Eight of the nine disjunct populations surveyed on the high‐elevation northeast plateau were fixed for a single haplotype (A), while the remaining, more westerly population, contained the same haplotype at high frequency together with two low frequency haplotypes (C and F). In contrast, most populations that occurred at lower altitudes at the plateau edge were fixed or nearly fixed for one of two haplotypes, A or E. However, two plateau edge populations had haplotype compositions different from the rest. In one, four haplotypes (A, B, D and E) were present at approximately equivalent frequencies, which might reflect a larger refugium in the area of this population during the last glacial period. Phylogenetic analysis indicated that the most widely distributed haplotype A is not ancestral to other haplotypes. The contrasting phylogeographic structures of the haplotype‐rich plateau edge area and the almost haplotype‐uniform plateau platform region indicate that the plateau platform was recolonized by J. przewalskii during the most recent postglacial period. This is supported by the findings of a nested clade analysis, which inferred that postglacial range expansion from the plateau edge followed by recent fragmentation is largely responsible for the present‐day spatial distribution of cpDNA haplotypes within the species.

Organelle DNA phylogeography of Cycas taitungensis, a relict species in Taiwan.

The phylogegraphic pattern of Cycas taitungensis, an endemic species with two remaining populations in Taiwan, was investigated based on genetic variability and phylogeny of the atpB‐rbcL noncoding spacer of chloroplast DNA (cpDNA) and the ribosomal DNA (rDNA) internal transcribed spacer (ITS) of mitochondrial DNA (mtDNA). High levels of genetic variation at both organelle loci, due to frequent intramolecular recombination, and low levels of genetic differentiation were detected in the relict gymnosperm. The apportionment of genetic variation within and between populations agreed with a migrant‐pool model, which describes a migratory pattern with colonists recruited from a random sample of earlier existing populations. Phylogenies obtained from cpDNA and mtDNA were discordant according to neighbour‐joining analyses. In total four chlorotypes (clades I–IV) and five mitotypes (clades A–E) were identified based on minimum spanning networks of each locus. Significant linkage disequilibrium in mitotype–chlorotype associations excluded the possibility of the recurrent homoplasious mutations as the major force causing phylogenetic inconsistency. The most abundant chlorotype I was associated with all mitotypes and the most abundant mitotype C with all chlorotypes; no combinations of rare mitotypes with rare chlorotypes were found. According to nested clade analyses, such nonrandom associations may be ascribed to relative ages among alleles associated with the geological history through which cycads evolved. Nested in networks as interior nodes coupled with wide geographical distribution, the most dominant cytotypes of CI and EI may represent ancestral haplotypes of C. taitungensis with a possible long existence prior to the Pleistocene glacial maximum. In contrast, rare chlorotypes and mitotypes with restricted and patchy distribution may have relatively recent origins. Newly evolved genetic elements of mtDNA, with a low frequency, were likely to be associated with the dominant chlorotype, and vice versa, resulting in the nonrandom mitotype–chlorotype associations. Paraphyly of CI and EI cytotypes, leading to the low level of genetic differentiation between cycad populations, indicated a short period for isolation, which allowed low possibilities of the attainment of coalescence at polymorphic ancestral alleles.

Conservation genetics and phylogeography of endangered and endemic shrub Tetraena mongolica (Zygophyllaceae) in Inner Mongolia, China

BackgroundTetraena mongolica (Zygophyllaceae), an endangered endemic species in western Inner Mongolia, China. For endemic species with a limited geographical range and declining populations, historical patterns of demography and hierarchical genetic structure are important for determining population structure, and also provide information for developing effective and sustainable management plans. In this study, we assess genetic variation, population structure, and phylogeography of T. mongolica from eight populations. Furthermore, we evaluate the conservation and management units to provide the information for conservation.ResultsSequence variation and spatial apportionment of the atp B-rbc L noncoding spacer region of the chloroplast DNA were used to reconstruct the phylogeography of T. mongolica. A total of 880 bp was sequenced from eight extant populations throughout the whole range of its distribution. At the cpDNA locus, high levels of genetic differentiation among populations and low levels of genetic variation within populations were detected, indicating that most seed dispersal was restricted within populations.ConclusionsDemographic fluctuations, which led to random losses of genetic polymorphisms from populations, due to frequent flooding of the Yellow River and human disturbance were indicated by the analysis of BEAST skyline plot. Nested clade analysis revealed that restricted gene flow with isolation by distance plus occasional long distance dispersal is the main evolutionary factor affecting the phylogeography and population structure of T. mongolica. For setting a conservation management plan, each population of T. mongolica should be recognized as a conservation unit.

Phylogeography of Kandelia candel in East Asiatic mangroves based on nucleotide variation of chloroplast and mitochondrial DNAs.

Vivipary with precocious seedlings in mangrove plants was thought to be a hindrance to long‐range dispersal. To examine the extent of seedling dispersal across oceans, we investigated the phylogeny and genetic structure among East Asiatic populations of Kandelia candel based on organelle DNAs. In total, three, 28 and seven haplotypes of the chloroplast DNA (cpDNA) atpB‐rbcL spacer, cpDNA trnL‐trnF spacer, and mitochondrial DNA (mtDNA) internal transcribed spacer (ITS) were identified, respectively, from 202 individuals. Three data sets suggested consistent phylogenies recovering two differentiated lineages corresponding to geographical regions, i.e. northern South‐China‐Sea + East‐China‐Sea region and southern South‐China‐Sea region (Sarawak). Phylogenetically, the Sarawak population was closely related to the Ranong population of western Peninsula Malaysia instead of other South‐China‐Sea populations, indicating its possible origin from the Indian Ocean Rim. No geographical subdivision was detected within the northern geographical region. An analysis of molecular variance (amova) revealed low levels of genetic differentiation between and within mainland and island populations (ΦCT = 0.015, ΦSC = 0.037), indicating conspicuous long‐distance seedling dispersal across oceans. Significant linkage disequilibrium excluded the possibility of recurrent homoplasious mutations as the major force causing phylogenetic discrepancy between mtDNA and the trnL‐trnF spacer within the northern region. Instead, relative ages of alleles contributed to nonrandom chlorotype–mitotype associations and tree inconsistency. Widespread distribution and random associations (χ2 = 0.822, P = 0.189) of eight hypothetical ancestral cytotypes indicated the panmixis of populations of the northern geographical region as a whole. In contrast, rare and recently evolved alleles were restricted to marginal populations, revealing some preferential directional migration.

Effects of GC Bias in Next-Generation-Sequencing Data on De Novo Genome Assembly

Next-generation-sequencing (NGS) has revolutionized the field of genome assembly because of its much higher data throughput and much lower cost compared with traditional Sanger sequencing. However, NGS poses new computational challenges to de novo genome assembly. Among the challenges, GC bias in NGS data is known to aggravate genome assembly. However, it is not clear to what extent GC bias affects genome assembly in general. In this work, we conduct a systematic analysis on the effects of GC bias on genome assembly. Our analyses reveal that GC bias only lowers assembly completeness when the degree of GC bias is above a threshold. At a strong GC bias, the assembly fragmentation due to GC bias can be explained by the low coverage of reads in the GC-poor or GC-rich regions of a genome. This effect is observed for all the assemblers under study. Increasing the total amount of NGS data thus rescues the assembly fragmentation because of GC bias. However, the amount of data needed for a full rescue depends on the distribution of GC contents. Both low and high coverage depths due to GC bias lower the accuracy of assembly. These pieces of information provide guidance toward a better de novo genome assembly in the presence of GC bias.

Phylogeography of North American populations of the moss species Hylocomium splendens based on the nucleotide sequence of internal transcribed spacer 2 of nuclear ribosomal DNA

Sequence variation of the internal transcribed sequence 2 (ITS2) region of nuclear ribosomal DNA (nrDNA) was investigated in 10 North American populations of Hylocomium splendens. Cladistic analyses supported the monophyly of this moss species, rooted at Hylocomiastrum and Neodolichomitra. Three geographically based groups (Great Lakes Forest, Appalachian Mountains, and Pacific Northwest) were identified by a minimum spanning network. Significant genetic differentiation was detected (FST = 0.197 − 0.390) among three geographical regions in North America. Although high genetic divergence exists within H. splendens, these results do not suggest sufficient divergence for designating sibling species.

Mitochondrial DNA phylogeography of Acrossocheilus paradoxus (Cyprinidae) in Taiwan

Nucleotide sequences of 3′ end of the cytochrome b gene, tRNA genes, D‐loop control region, and the 5′ end of the 12S rRNA of mitochondrial DNA (mtDNA) were used to assess the genetic and phylogeographic structure of Acrossocheilus paradoxus populations, a Cyprinidae fish of Taiwan. A hierarchical examination of populations in 12 major streams from three geographical regions using an analysis of molecular variance ( amova) indicates high genetic differentiation both among populations (ΦST = 0.511, P < 0.001) and among regions (ΦCT = 0.368, P < 0.001). Limited migration largely contributed to the genetic differentiation. High nucleotide diversity (1.13%) and haplotype diversity (0.80%) were detected among populations. The degree of genetic differentiation was correlated with geographical distance between populations, a result consistent with the one‐dimensional stepping stone models. A neighbour‐joining tree recovered by ( dambe) supports the pattern of isolation by distance and reveals a closer relationship between populations of the central and southern regions. A minimum spanning network based on nucleotide substitutions reflected migration routes from populations of the central region to the northern and southern regions, respectively. Postglacial colonization and expansion can explain the phylogeographical pattern. Single and ancient migration events may have allowed the northern region to attain the monophyly of mtDNA alleles. In contrast, most populations within geographical regions are either paraphyletic or polyphyletic due to the relatively shorter time period for coalescence. Both low haplotype number and genetic variability suggest a bottleneck event in the Chingmei population of northern Taiwan. Based on coalescence theory, the monophyly of the Tungkang population of the southern region may be associated with a founder event.

Chloroplast DNA phylogeography of Cunninghamia konishii (Cupressaceae), an endemic conifer of Taiwan.

In this study, we investigated the genetic structure and phylogeographic pattern of the genus Cunninghamia, a member of the Cupressaceae restricted to mainland China and Taiwan, based on sequences of the trnD-trnT noncoding spacer of the chloroplast DNA. Maternal inheritance of chloroplasts was determined experimentally. No paternal leakage was detected. Both parsimony and neighbor-joining analyses revealed the polyphyly of Cunninghamia konishii, populations of which were nested in clades of C. lanceolata from mainland China. The nucleotide diversity of chloroplast DNA sequences within C. konishii (0.0118) was higher than that between species (0.0104), which agrees with a previous allozyme investigation. Based on mutational differences between sequences, a minimum spanning network consisting of five clades was constructed. Significant genetic differentiation (phiST = 0.130, P < 0.001) was detected between the clades based on AMOVA analyses. We infer several possible refugia in the Yunnan, Zhejiang, and Guangdong provinces of south China, all located in the minimum network as interior nodes. We also infer possible migration routes of Cunninghamia populations. The phylogeographic pattern shown in the reconstructed network suggests that the present-day Cunninghamia populations in Taiwan were derived from six different sources in continental Asia via long-distance seed dispersal. A migrant-pool model explains the heterogeneous composition of the organelle DNA in Taiwans populations and the low differentiation between populations of Taiwan and China (phiCT = 0.012, P = 0.454). In contrast with the genetic heterogeneity within geographic populations, many local populations have attained coalescence at the trnD-trnT alleles, which has led to significant differentiation at the population level.

Contrasting selection modes at the Adh1 locus in outcrossing Miscanthus sinensis vs. inbreeding Miscanthus condensatus (Poaceae)

We estimated DNA sequence variation of the Adh1 locus in the outcrossing Miscanthus sinensis (Poaceae) and its close selfing relative, M. condensatus. Tajimas test of selection is significantly negative for both overall exons and replacement sites in M. sinensis. Among its entire sample, nucleotide diversity of nonsynonymous sites is higher than that of synonymous sites. A McDonald and Kreitman test of neutrality indicates an excess of intraspecific replacement polymorphisms, suggesting possible directional selection toward advantageous mutants. However, frequent intragenic recombination suggests both purifying and positive selection is unlikely. Recent demographic expansions coupled with relaxation of purifying selection may have resulted in elevated genetic diversity at the Adh1 locus as well as the trnL-trnF intergenic spacer of cpDNA in this outcrossing species. In contrast, low levels of genetic diversity were detected at both the Adh1 locus and the cpDNA spacer in M. condensatus, consistent with bottlenecks associated with selfing in all populations. While Tajimas D and Fu and Lis F statistics did not reveal deviation from neutrality at the Adh1 locus in M. condensatus, 12 replacements vs. 10 synonymous changes were detected. Based on pairwise comparisons of the d(N)/d(S) ratio, lineages of closely related populations of the species distributed along saline habitats appeared to be under directional selection.

Low genetic diversity and significant population structuring in the relict Amentotaxus argotaenia complex (Taxaceae) based on ISSR fingerprinting

Amentotaxus, a genus of the Taxaceae, represents an ancient lineage that has long existed in Eurasia. All Amentotaxus species experienced frequent population expansion and contraction over periodical glaciations in Tertiary and Quaternary. Among them, Amentotaxus argotaenia complex consists of three morphologically alike species, A. argotaenia, Amentotaxus yunnanensis, and Amentotaxus formosana. This complex is distributed in the subtropical region of mainland China and Taiwan where many Pleistocene refugia have been documented. In this study, genetic diversity and population structuring within and between species were investigated based on the inter-simple sequence repeats (ISSR) fingerprinting. Mean genetic diversity within populations was estimated in three ways: (1) the percentage of polymorphic loci out of all loci (P) (2) Nei’s unbiased expected heterozygosity (He), and (3) Shannon’s index of phenotypic diversity. For a total of 310 individuals of 15 populations sampled from the three species, low levels of ISSR genetic variation within populations were detected, with P=4.66–16.58%, He=0.0176–0.0645 and Hpop=0.0263–0.0939, agreeing with their seriously threatened status. AMOVA analyses revealed that the differences between species only accounted for 27.38% of the total variation, whereas differences among populations and within populations were 57.70 and 14.92%, respectively, indicating substantial isolation between the patch-like populations. A neighbor-joining tree identified a close affinity between A. yunnanensis and A. formosana. Genetic drift due to small population size, plus limited current gene flow, resulted in significant genetic structuring. Low levels of intrapopulational genetic variation and considerable interpopulational divergence were also attributable to demographic bottlenecks during and/or after the Pleistocene glaciations.