Chao Li Huang

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.

Adaptive divergence with gene flow in incipient speciation of Miscanthus floridulus/sinensis complex (Poaceae).

Young incipient species provide ideal materials for untangling the process of ecological speciation in the presence of gene flow. The Miscanthus floridulus/sinensis complex exhibits diverse phenotypic and ecological differences despite recent divergence (approximately 1.59 million years ago). To elucidate the process of genetic differentiation during early stages of ecological speciation, we analyzed genomic divergence in the Miscanthus complex using 72 randomly selected genes from a newly assembled transcriptome. In this study, rampant gene flow was detected between species, estimated as M = 3.36 × 10(-9) to 1.20 × 10(-6) , resulting in contradicting phylogenies across loci. Nevertheless, beast analyses revealed the species identity and the effects of extrinsic cohesive forces that counteracted the non-stop introgression. As expected, early in speciation with gene flow, only 3-13 loci were highly diverged; two to five outliers (approximately 2.78-6.94% of the genome) were characterized by strong linkage disequilibrium, and asymmetrically distributed among ecotypes, indicating footprints of diversifying selection. In conclusion, ecological speciation of incipient species of Miscanthus probably followed the parapatric model, whereas allopatric speciation cannot be completely ruled out, especially between the geographically isolated northern and southern M. sinensis, for which no significant gene flow across oceanic barriers was detected. Divergence between local ecotypes in early-stage speciation began at a few genomic regions under the influence of natural selection and divergence hitchhiking that overcame gene flow.

Ecological genomics in Xanthomonas : the nature of genetic adaptation with homologous recombination and host shifts

BackgroundComparative genomics provides insights into the diversification of bacterial species. Bacterial speciation usually takes place with lasting homologous recombination, which not only acts as a cohering force between diverging lineages but brings advantageous alleles favored by natural selection, and results in ecologically distinct species, e.g., frequent host shift in Xanthomonas pathogenic to various plants.ResultsUsing whole-genome sequences, we examined the genetic divergence in Xanthomonas campestris that infected Brassicaceae, and X. citri, pathogenic to a wider host range. Genetic differentiation between two incipient races of X. citri pv. mangiferaeindicae was attributable to a DNA fragment introduced by phages. In contrast to most portions of the genome that had nearly equivalent levels of genetic divergence between subspecies as a result of the accumulation of point mutations, 10% of the core genome involving with homologous recombination contributed to the diversification in Xanthomonas, as revealed by the correlation between homologous recombination and genomic divergence. Interestingly, 179 genes were under positive selection; 98 (54.7%) of these genes were involved in homologous recombination, indicating that foreign genetic fragments may have caused the adaptive diversification, especially in lineages with nutritional transitions. Homologous recombination may have provided genetic materials for the natural selection, and host shifts likely triggered ecological adaptation in Xanthomonas. To a certain extent, we observed positive selection nevertheless contributed to ecological divergence beyond host shifting.ConclusionAltogether, mediated with lasting gene flow, species formation in Xanthomonas was likely governed by natural selection that played a key role in helping the deviating populations to explore novel niches (hosts) or respond to environmental cues, subsequently triggering species diversification.

Amyloid core formed of full-length recombinant mouse prion protein involves sequence 127-143 but not sequence 107-126.

The principal event underlying the development of prion disease is the conversion of soluble cellular prion protein (PrPC) into its disease-causing isoform, PrPSc. This conversion is associated with a marked change in secondary structure from predominantly α-helical to a high β-sheet content, ultimately leading to the formation of aggregates consisting of ordered fibrillar assemblies referred to as amyloid. In vitro, recombinant prion proteins and short prion peptides from various species have been shown to form amyloid under various conditions and it has been proposed that, theoretically, any protein and peptide could form amyloid under appropriate conditions. To identify the peptide segment involved in the amyloid core formed from recombinant full-length mouse prion protein mPrP(23–230), we carried out seed-induced amyloid formation from recombinant prion protein in the presence of seeds generated from the short prion peptides mPrP(107–143), mPrP(107–126), and mPrP(127–143). Our results showed that the amyloid fibrils formed from mPrP(107–143) and mPrP(127–143), but not those formed from mPrP(107–126), were able to seed the amyloidogenesis of mPrP(23–230), showing that the segment residing in sequence 127–143 was used to form the amyloid core in the fibrillization of mPrP(23–230).

Highly diverse endophytes in roots of Cycas bifida (Cycadaceae), an ancient but endangered gymnosperm

As an ancient seed plant, cycads are one of the few gymnosperms that develop a root symbiosis with cyanobacteria, which has allowed cycads to cope with harsh geologic and climatic conditions during the evolutionary process. However, the endophytic microbes in cycad roots remain poorly identified. In this study, using next-generation sequencing techniques, we investigated the microbial diversity and composition of both the coralloid and regular roots of Cycas bifida (Dyer) K.D. Hill. Highly diverse endophytic communities were observed in both the coralloid and regular roots. Of the associated bacteria, the top five families were the Nostocaceae, Sinobacteraceae, Bradyrhizobiaceae, Bacillaceae, and Hyphomicrobiaceae. The Nectriaceae, Trichocomaceae, and Incertae sedis were the predominant fungal families in all root samples. A significant difference in the endophytic bacterial community was detected between coralloid roots and regular roots, but no difference was observed between the fungal communities in the two root types. Cyanobacteria were more dominant in coralloid roots than in regular roots. The divergence of cycad root structures and the modified physiological processes may have contributed to the abundance of cyanobionts in coralloid roots. Consequently, the colonization of cyanobacteria inhibits the assemblage of other endophytes. Our results contribute to an understanding of the species diversity and composition of the cycad-endophyte microbiome and provide an abbreviated list of potential ecological roles of the core microbes present.

Common stress transcriptome analysis reveals functional and genomic architecture differences between early and delayed response genes

To identify the similarities among responses to diverse environmental stresses, we analyzed the transcriptome response of rice roots to three rhizotoxic perturbations (chromium, ferulic acid and mercury) and identified common early-transient, early-constant and delayed gene inductions. Common early response genes were mostly associated with signal transduction and hormones, and delayed response genes with lipid metabolism. Network component analysis revealed complicated interactions among common genes, the most highly connected signaling hubs being PP2C68, MPK5, LRR-RLK and NPR1. Gene architecture studies revealed different conserved promoter motifs and a different ratio of CpG island distribution between early and delayed genes. In addition, early-transient genes had more exons and a shorter first exon. IMEter was used to calculate the transcription regulation effects of introns, with greater effects for the first introns of early-transient than delayed genes. The higher Ka/Ks (non-synonymous/synonymous mutation) ratio of early-constant genes than early-transient, delayed and the genome median demonstrates the rapid evolution of early-constant genes. Our results suggest that finely tuned transcriptional control in response to environmental stress in rice depends on genomic architecture and signal intensity and duration.

Frequent gene flow blurred taxonomic boundaries of sections in Lilium L. (Liliaceae)

Gene flow between species may last a long time in plants. Reticulation inevitably causes difficulties in phylogenetic reconstruction. In this study, we looked into the genetic divergence and phylogeny of 20 Lilium species based on multilocus analyses of 8 genes of chloroplast DNA (cpDNA), the internally transcribed nuclear ribosomal DNA (nrITS) spacer and 20 loci extracted from the expressed sequence tag (EST) libraries of L. longiflorum Thunb. and L. formosanum Wallace. The phylogeny based on the combined data of the maternally inherited cpDNA and nrITS was largely consistent with the taxonomy of Lilium sections. This phylogeny was deemed the hypothetical species tree and uncovered three groups, i.e., Cluster A consisting of 4 taxa from the sections Pseudolirium and Liriotypus, Cluster B consisting of the 4 taxa from the sections Leucolirion, Archelirion and Daurolirion, and Cluster C comprising 10 taxa mostly from the sections Martagon and Sinomartagon. In contrast, systematic inconsistency occurred across the EST loci, with up to 19 genes (95%) displaying tree topologies deviating from the hypothetical species tree. The phylogenetic incongruence was likely attributable to the frequent genetic exchanges between species/sections, as indicated by the high levels of genetic recombination and the IMa analyses with the EST loci. Nevertheless, multilocus analysis could provide complementary information among the loci on the species split and the extent of gene flow between the species. In conclusion, this study not only detected frequent gene flow among Lilium sections that resulted in phylogenetic incongruence but also reconstructed a hypothetical species tree that gave insights into the nature of the complex relationships among Lilium species.

Multilocus Analyses Reveal Postglacial Demographic Shrinkage of Juniperus morrisonicola (Cupressaceae), a Dominant Alpine Species in Taiwan

Postglacial climate changes alter geographical distributions and diversity of species. Such ongoing changes often force species to migrate along the latitude/altitude. Altitudinal gradients represent assemblage of environmental, especially climatic, variable factors that influence the plant distributions. Global warming that triggered upward migrations has therefore impacted the alpine plants on an island. In this study, we examined the genetic structure of Juniperus morrisonicola, a dominant alpine species in Taiwan, and inferred historical, demographic dynamics based on multilocus analyses. Lower levels of genetic diversity in north indicated that populations at higher latitudes were vulnerable to climate change, possibly related to historical alpine glaciers. Neither organellar DNA nor nuclear genes displayed geographical subdivisions, indicating that populations were likely interconnected before migrating upward to isolated mountain peaks, providing low possibilities of seed/pollen dispersal across mountain ranges. Bayesian skyline plots suggested steady population growth of J. morrisonicola followed by recent demographic contraction. In contrast, most lower-elevation plants experienced recent demographic expansion as a result of global warming. The endemic alpine conifer may have experienced dramatic climate changes over the alternation of glacial and interglacial periods, as indicated by a trend showing decreasing genetic diversity with the altitudinal gradient, plus a fact of upward migration.