Tsai Wen Hsu

Contrasting phylogeographical patterns between mainland and island taxa of the Pinus luchuensis complex

Species whose geographical distribution encompasses both mainland and island populations provide an ideal system for examining isolation and genetic divergence. In this study, paternally transmitted chloroplast DNA (cpDNA) and maternally transmitted mitochondrial DNA (mtDNA) were used to estimate population structure and phylogeography of Pinus luchuensis, a species found in eastern China (ssp. hwangshanensis), Taiwan (ssp. taiwanensis), and the Ryukyu Archipelago (ssp. luchuensis). Gene genealogies of both mtDNA and cpDNA reveal two major lineages. Molecular dating indicates that these lineages diverged before the colonization of P. luchuensis subspecies in Taiwan and the Ryukyu Archipelago. Both mtDNA and cpDNA show a lack of correspondence between molecular phylogeny and subspecies designation. Phylogeographical analysis suggests that paraphyly of the subspecies is the result of recent divergence rather than secondary contacts. In spite of the short divergence history of P. luchuensis on islands, the island populations show the same degree of genetic divergence as mainland populations. Low levels of genetic diversity in the mainland ssp. hwangshanensis suggest demographic bottlenecks. In contrast, the high heterogeneity of genetic composition for island populations is likely to be associated with a history of multiple colonization from the mainland. The spatial apportionment of organelle DNA polymorphisms is consistent with a pattern of stepwise colonization on island populations.

Paraphyly of organelle DNAs in Cycas Sect. Asiorientales due to ancient ancestral polymorphisms

BackgroundThis study addresses the apportionment of genetic diversity between Cycas revoluta and C. taitungensis, species that constitute the section Asiorientales and represent a unique, basal lineage of the Laurasian genus Cycas. Fossil evidence indicates divergence of the section from the rest of Cycas at least 30 million years ago. Geographically, C. taitungensis is limited to Taiwan whereas C. revoluta is found in the Ryukyu Archipelago and on mainland China.ResultsThe phylogenies of ribosomal ITS region of mtDNA and the intergenic spacer between atpB and rbcL genes of cpDNA were reconstructed. Phylogenetic analyses revealed paraphyly of both loci in the two species and also in the section Asiorientales. The lack of reciprocal monophyly between these long isolated sections is likely due to persistent shared ancestral polymorphisms. Molecular dating estimated that mt- and cp DNA lineages coalesced to the most recent common ancestors (TMRCA) about 327 (mt) and 204 MYA (cp), corresponding with the divergence of cycad sections in the Mesozoic.ConclusionFates of newly derived mutations of cycads follow Klopfstein et al.s surfing model where the majority of new mutations do not spread geographically and remain at low frequencies or are eventually lost by genetic drift. Only successful surfing mutations reach very high frequencies and occupy a large portion of a species range. These mutations exist as dominant cytotypes across populations and species. Geographical subdivision is lacking in both species, even though recurrent gene flow by both pollen and seed is severely limited. In total, the contrasting levels between historical and ongoing gene flow, large population sizes, a long lifespan, and slow mutation rates in both organelle DNAs have all likely contributed to the unusually long duration of paraphyly in cycads.

Footprints of natural and artificial selection for photoperiod pathway genes in Oryza

Asian rice, Oryza sativa, consists of two major subspecies, indica and japonica, which are physiologically differentiated and adapted to different latitudes. Genes for photoperiod sensitivity are likely targets of selection along latitude. We examined the footprints of natural and artificial selections for four major genes of the photoperiod pathway, namely PHYTOCHROME B (PhyB), HEADING DATE 1 (Hd1), HEADING DATE 3a (Hd3a), and EARLY HEADING DATE 1 (Ehd1), by investigation of the patterns of nucleotide polymorphisms in cultivated and wild rice. Geographical subdivision between tropical and subtropical O.xa0rufipogon was found for all of the photoperiod genes in plants divided by the Tropic of Cancer (TOC). All of these genes, except for PhyB, were characterized by the existence of clades that split a long time ago and that corresponded to latitudinal subdivisions, and revealed a likely diversifying selection. Ssp. indica showed close affinity to tropical O.xa0rufipogon for all genes, while ssp. japonica, which has a much wider range of distribution, displayed complex patterns of differentiation from O.xa0rufipogon, which reflected various agricultural needs in relation to crop yield. In japonica, all genes, except Hd3a, were genetically differentiated at the TOC, while geographical subdivision occurred at 31°N in Hd3a, probably the result of varying photoperiods. Many other features of the photoperiod genes revealed domestication signatures, which included high linkage disequilibrium (LD) within genes, the occurrence of frequent and recurrent non-functional Hd1 mutants in cultivated rice, crossovers between subtropical and tropical alleles of Hd1, and significant LD between Hd1 and Hd3a in japonica and indica.

Inferring multiple refugia and phylogeographical patterns in Pinus massoniana based on nucleotide sequence variation and DNA fingerprinting.

Background Pinus massoniana, an ecologically and economically important conifer, is widespread across central and southern mainland China and Taiwan. In this study, we tested the central–marginal paradigm that predicts that the marginal populations tend to be less polymorphic than the central ones in their genetic composition, and examined a founders effect in the island population. Methodology/Principal Findings We examined the phylogeography and population structuring of the P. massoniana based on nucleotide sequences of cpDNA atpB-rbcL intergenic spacer, intron regions of the AdhC2 locus, and microsatellite fingerprints. SAMOVA analysis of nucleotide sequences indicated that most genetic variants resided among geographical regions. High levels of genetic diversity in the marginal populations in the south region, a pattern seemingly contradicting the central–marginal paradigm, and the fixation of private haplotypes in most populations indicate that multiple refugia may have existed over the glacial maxima. STRUCTURE analyses on microsatellites revealed that genetic structure of mainland populations was mediated with recent genetic exchanges mostly via pollen flow, and that the genetic composition in east region was intermixed between south and west regions, a pattern likely shaped by gene introgression and maintenance of ancestral polymorphisms. As expected, the small island population in Taiwan was genetically differentiated from mainland populations. Conclusions/Significance The marginal populations in south region possessed divergent gene pools, suggesting that the past glaciations might have low impacts on these populations at low latitudes. Estimates of ancestral population sizes interestingly reflect a recent expansion in mainland from a rather smaller population, a pattern that seemingly agrees with the pollen record.

Development of 12 genic microsatellite loci for a biofuel grass, Miscanthus sinensis (Poaceae)

PREMISE OF THE STUDYnMiscanthus, a nonfood plant with high potential as a biofuel, has been used in Europe and the United States. The selection of a cultivar with high biomass, photosynthetic efficiency, and stress resistance from wild populations has become an important issue. New genic microsatellite markers will aid the assessment of genetic diversity for different strains.nnnMETHODS AND RESULTSnTwelve polymorphic microsatellite markers derived from the transcriptome of Miscanthus sinensis fo. glaber were identified and screened on 80 individuals of M. sinensis. The number of alleles per locus ranged from 6 to 12, and the mean expected heterozygosity was 0.75. Cross-taxa transferability revealed that all loci can be applied to all varieties of M. sinensis, as well as the closely related species M. floridulus.nnnCONCLUSIONSnThese new genic microsatellite markers are useful for characterizing different traits in breeding programs or to select genes useful for biofuel.

Multilocus analysis of genetic divergence between outcrossing Arabidopsis species: evidence of genome‐wide admixture

• Outcrossing Arabidopsis species that diverged from their inbreeding relative Arabidopsis thaliana 5 million yr ago and display a biogeographical pattern of interspecific sympatry vs intraspecific allopatry provides an ideal model for studying impacts of gene introgression and polyploidization on species diversification. • Flow cytometry analyses detected ploidy polymorphisms of 2× and 4× in Arabidopsis lyrata ssp. kamchatica of Taiwan. Genomic divergence between species/subspecies was estimated based on 98 randomly chosen nuclear genes. Multilocus analyses revealed a mosaic genome in diploid A. l. kamchatica composed of Arabidopsis halleri-like and A. lyrata-like alleles. • Coalescent analyses suggest that the segregation of ancestral polymorphisms alone cannot explain the high inconsistency between gene trees across loci, and that gene introgression via diploid A. l. kamchatica likely distorts the molecular phylogenies of Arabidopsis species. However, not all genes migrated across species freely. Gene ontology analyses suggested that some nonmigrating genes were constrained by natural selection. • High levels of estimated ancestral polymorphisms between A. halleri and A. lyrata suggest that gene flow between these species has not completely ceased since their initial isolation. Polymorphism data of extant populations also imply recent gene flow between the species. Our study reveals that interspecific gene flow affects the genome evolution in Arabidopsis.

Evolutionary rates of commonly used nuclear and organelle markers of Arabidopsis relatives (Brassicaceae)

Recovering the genetic divergence between species is one of the major interests in the evolutionary biology. It requires accurate estimation of the neutral substitution rates. Arabidopsis thaliana, the first whole-genome sequenced plant, and its out-crossing relatives provide an ideal model for examining the split between sister species. In the study, rates of molecular evolution at markers frequently used for systematics and population genetics, including 14 nuclear genes spanning most chromosomes, three noncoding regions of chloroplast genome, and one intron of mitochondrial genome, between A. thaliana and four relatives were estimated. No deviation from neutrality was detected in the genes examined. Based on the known divergence between A. thaliana and its sisters about 8.0-17.6 MYA, evolutionary rates of the eighteen genes were estimated. Accordingly, the ratio of rates of synonymous substitutions among mitochondrial, chloroplast and nuclear genes was calculated with an average and 95% confidence interval of 1 (0.25-1.75): 15.77 (7.48-114.09): 74.79 (36.27-534.61). Molecular evolutionary rates of nuclear genes varied, with a range of 0.383-0.856×10(-8) for synonymous substitutions per site per year and 0.036-0.081×10(-9) for nonsynonymous substitutions per site per year. Compared with orthologs in Populus, a long life-span tree, genes in Arabidopsis evolved faster in an order of magnitude at the gene level, agreeing with a generation time hypothesis. The estimated substitution rates of these genes can be used as a reference for molecular dating.

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

BackgroundA 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.ResultsSystematic 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 (FST = 0.56-0.72) between populations in the Yushan Mountain Range and populations of the other mountain ranges.ConclusionBoth 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/founders events.

Characterization of 42 Microsatellite Markers from Poison Ivy, Toxicodendron radicans (Anacardiaceae)

Poison ivy, Toxicodendron radicans, and poison oaks, T. diversilobum and T. pubescens, are perennial woody species of the Anacardiaceae and are poisonous, containing strong allergens named urushiols that cause allergic contact dermatitis. Poison ivy is a species distributed from North America to East Asia, while T. diversilobum and T. pubescens are distributed in western and eastern North America, respectively. Phylogreography and population structure of these species remain unclear. Here, we developed microsatellite markers, via constructing a magnetic enriched microsatellite library, from poison ivy. We designed 51 primer pairs, 42 of which successfully yielded products that were subsequently tested for polymorphism in poison oak, and three subspecies of poison ivy. Among the 42 loci, 38 are polymorphic, while 4 are monomorphic. The number of alleles and the expected heterozygosity ranged from 1 to 12 and from 0.10 to 0.87, respectively, in poison ivy, while varied from 2 to 8 and, from 0.26 to 0.83, respectively in poison oak. Genetic analysis revealed distinct differentiation between poison ivy and poison oak, whereas slight genetic differentiation was detected among three subspecies of poison ivy. These highly polymorphic microsatellite fingerprints enable biologists to explore the population genetics, phylogeography, and speciation in Toxicodendron.

Genetic Divergence and Biogeographical Patterns in Amentotaxus argotaenia Species Complex

The Amentotaxus argotaenia (Taxaceae) species complex is comprised of four relict and endangered gymnosperms. Amentotaxus argotaenia and A. yunnanensis are scattered throughout South and Central China, and A. formosana, and A. poilanei are restricted to Taiwan and Vietnam. In the present study, we examined genetic divergence and biogeographical patterns in the species complex based on chloroplast DNA atpI-atpH and rpS16-trnK intergenic spacers, the mitochondrial DNA nad1 intron, and microsatellite markers. Genealogical analyses revealed systematic inconsistencies between organellar DNA markers, with paraphyletic species inferred based on cpDNA, versus a monophyletic A. formosana, and clustering of A. poilanei within A. yunnanensis at mtDNA. AMOVA revealed that most organellar DNA variants resided among species and populations, a pattern also supported by clustering with STRUCTURE analyses on microsatellites. Higher levels of genetic diversity in the Yongxin population of A. argotaenia reflected features of a refugium. Most speciation events were rather recent, with divergence no earlier than 2.50 MYA, although many ancient lineages derived 4.12–13.84 MYA were maintained within species. Significant genetic structuring within A. argotaenia was detected among three geographical regions—a phylogeographical pattern likely attributable to past fragmentations. In contrast to paraphyletic A. argotaenia based on the maintenance of ancestral polymorphisms, interspecific hybridization likely contributed to polyphyly in A. yunnanensis. Furthermore, despite the low likelihood of recurrent gene flow due to geographical isolation, IM analyses revealed continued interspecific gene flow after species divergence, reflecting a parapatric speciation model.