Phylogenomic analyses of Camellia support reticulate evolution among major clades.

Camellia (Theaceae) is a morphologically highly diverse genus of flowering plants and includes many famous species with high economic value, and the phylogeny of this genus is not fully resolved. We used 95 transcriptomes from 87 Camellia species and identified 1481 low-copy genes to conduct a detailed analysis of the phylogeny of this genus according to various data-screening criteria. The results show that, very different from the two existing classification systems of Camellia, 87 species are grouped into 8 main clades and two independent species, and that all 8 clades except Clade 8 were strongly supported by almost all the coalescent or concatenated trees using different gene subsets. However, the relationships among these clades were weakly supported and different from analyses using different gene subsets; furthermore, they do not agree with the phylogeny from chloroplast genomes of Camellia. Additional analyses support reticulate evolution (probably resulting from introgression or hybridization) among some major Camellia lineages, providing explanation for extensive gene tree conflicts. Furthermore, we inferred that together with the formation of East Asian subtropical evergreen broad-leaved forests, Camellia underwent a radiative divergence of major clades at 23 âˆ¼ 19 Ma in the late Miocene then had a subsequent species burst at 10 âˆ¼ 5 Ma. Principal component and cluster analyses provides new insights into morphological changes underlying the evolution of Camellia and a reference to further clarify subgenus and sections of this genus. The comprehensive study here including a nuclear phylogeny and other analyses reveal the rapid evolutionary history of Camellia.

Comparative morphological study of skeletal muscle weight among the red jungle fowl (Gallus gallus) and various fowl breeds (Gallus domesticus).

We examined the weight distribution of skeletal muscles of the red jungle fowl, then compared these values with those of domesticated populations to determine how muscle distribution has changed by selecting breeding. Sonia, Fayoumi, and Rhode Island Red were selected for comparison from livestock breeds, while Japanese Shamo and Thai fighting cocks were selected from cockfighting groups. Principal component analysis was applied using body size-free data. The mass distribution of muscles clearly differed between the wild, livestock, and cockfighting groups, demonstrating that muscle distribution has changed after selecting breeding, coupled with functional demands of each group. The red jungle fowl, which has the ability to fly, could be clearly distinguished from the flightless domesticated populations due to differences in flight pectoral muscle size. The cervical muscles in the wild population were smaller than in the domesticated groups; these do not contribute to flight. The gluteal muscles were larger in the fighting cock group, functionally coupled to their traditionally preferred upright posture. Wild bird populations typically exhibit reduced weight of their hind limbs, associated with flight, but as the red jungle fowl displays largely terrestrial behavior, these muscles are similar in arrangement and relative size to those of the livestock groups. We showed that the mass distribution pattern of skeletal muscles expresses selecting breeding strategy and clearly reflects the specific traits for each group.

Fine-scale genetic diversity and putative ecotypes of oxymonad protists coinhabiting the hindgut of Reticulitermes speratus.

The hindgut of lower termites is generally coinhabited by multiple morphologically identifiable protist species. However, it is unclear how many protist species truly coexist in this miniaturized environment, and moreover, it is difficult to define the fundamental unit of protist diversity. Species delineation of termite gut protists has therefore been guided without a theory-based concept of species. Here, we focused on the hindgut of the termite Reticulitermes speratus, where 10 or 11 morphologically distinct oxymonad cell types, that is, morphospecies, coexist. We elucidated the phylogenetic structure of all co-occurring oxymonads and addressed whether their diversity can be explained by the "ecotype" hypothesis. Oxymonad-specific 18S rRNA gene amplicon sequencing analyses of whole-gut samples, combined with single-cell 18S rRNA sequencing of the oxymonad morphospecies, identified 210 one-nucleotide-level variants. The phylogenetic analysis of these variants revealed the presence of microdiverse clusters typically within 1% sequence divergence. Each known oxymonad morphospecies comprised one to several monophyletic or paraphyletic microdiverse clusters. Using these sequence data sets, we conducted computational simulation to predict the rates of ecotype formation and periodic selection, and to demarcate putative ecotypes. Our simulations suggested that the oxymonad genetic divergence is constrained primarily by strong selection, in spite of limited population size and possible bottlenecks during intergenerational transmission. A total of 33 oxymonad ecotypes were predicted, and most of the putative ecotypes were consistently detected among different colonies and host individuals. These findings provide a possible theoretical basis for species diversity and underlying mechanisms of coexistence of termite gut protists.

Concerted and birth-and-death evolution of 26S ribosomal DNA in Camellia L.

BACKGROUND AND AIMS: The ribosomal DNA (rDNA) gene family, encoding ribosomal RNA (rRNA), has long been regarded as an archetypal example illustrating the model of concerted evolution. However, controversy is arising, as rDNA in many eukaryotic species has been proved to be polymorphic. Here, a metagenomic strategy was applied to detect the intragenomic polymorphism as well as the evolutionary patterns of 26S rDNA across the genus Camellia. METHODS: Degenerate primer pairs were designed to amplify the 26S rDNA fragments from different Camellia species. The amplicons were then paired-end sequenced on the Illumina MiSeq platform. KEY RESULTS: An extremely high level of rDNA polymorphism existed universally in Camellia. However, functional rDNA was still the major component of the family, and was relatively conserved among different Camellia species. Sequence variations mainly came from rRNA pseudogenes and favoured regions that are rich in GC. Specifically, some rRNA pseudogenes have existed in the genome for a long time, and have even experienced several expansion events, which has greatly enriched the abundance of rDNA polymorphism. CONCLUSIONS: Camellia represents a group in which rDNA is subjected to a mixture of concerted and birth-and-death evolution. Some rRNA pseudogenes may still have potential functions. Conversely, when released from selection constraint, they can evolve in the direction of decreasing GC content and structural stability through a methylation-induced process, and finally be eliminated from the genome.

Evaluation of reported sediment samples from 20 Ma using a molecular phylogenetic approach: comment on Liu et al. (2017).

Liu et al. reported the cultivation and DNA sequencing of 69 fungal isolates (Ascomycota and Basidiomycota) from ancient subseafloor sediments, suggesting that they represent living fungal populations that have persisted for over 20 million years. Because these findings could bring about a paradigm shift in our understanding of the spatial breadth of the deep subsurface biosphere as well as the longevity of ancient DNA, it is extremely important to verify that their samples represent pure ancient fungi from 20 million years ago without contamination by modern species. For this purpose, we estimated the divergence times of Dikarya (Ascomycota + Basidiomycota) and Mucoromycota fungi assuming that the fungal isolates were actually sampled from 20 Ma (mega-annum) sediments and evaluated the validity of the sample ages. Using this approach, we estimate that the age of the last common ancestor of Dikarya and Mucoromycota fungi greatly exceeds the age of the Earth. Our finding emphasizes the importance of using reliable approaches to confirm the dating of ancient samples.

Phylogeographic Analysis of Madagascan Goats Using mtDNA Control Region and SRY Gene Sequences.

In the present study, we estimated the genetic diversity and relationships, as well as the propagation routes, of Madagascan goats using mtDNA control region and SRY gene sequences. The mtDNA sequences of 40 Madagascan goats revealed 10 haplotypes and a quite low nucleotide diversity (0.0014), suggesting a founder and/or bottleneck effect resulting from goat migration to Madagascar island. The analysis of sequences identical to Madagascan haplotypes indicated close genetic relationships between goats from Madagascar and Africa. Sequence analysis of the SRY gene in 40 male Madagascan goats revealed two haplotypes: Y1A (45%) and Y2A (55%). The paternal result indicated genetic influences from Africa, South Asia, and the Near East proximal to Madagascar. The analyses of the mtDNA control region and SRY gene sequences suggested a genetic relationship between Africa and Madagascar. Moreover, SRY sequences indicated influences from South Asia and the Near East. These phylogenetic results provide important genetic information for elucidating the propagation routes of Madagascan goats.

Demographic analysis of cyanobacteria based on the mutation rates estimated from an ancient ice core.

Despite the crucial role of cyanobacteria in various ecosystems, little is known about their evolutionary histories, especially microevolutionary dynamics, because of the lack of knowledge regarding their mutation rates. Here we directly estimated cyanobacterial mutation rates based on ancient DNA analyses of ice core samples collected from Kyrgyz Republic that dates back to ~12,500 cal years before present. We successfully sequenced the 16S rRNA and 16S-23S internal transcribed spacer (ITS) region. Two cyanobacterial operational taxonomic units (OTUs) were detected from the ancient ice core samples, and these OTUs are shared with those from the modern glacier surface. The mutation rate of ITS region was estimated by comparing ancient and modern populations, and were at the magnitude of 10-7substitutions/sites/year. By using a model selection framework, we also demonstrated that the ancient sequences from the ice sample were not contaminated from modern samples. Bayesian demographic analysis based on coalescent theory revealed that cyanobacterial population sizes increased over Asia regions during the Holocene. Thus, our results enhance our understanding of the enigmatic timescale of cyanobacterial microevolution, which has the potential to elucidate the environmental responses of cyanobacteria to the drastic climatic change events of the Quaternary.

Phylogeography of Sophora moorcroftiana Supports Wu's Hypothesis on the Origin of Tibetan Alpine Flora.

Wu hypothesized that the Tibetan flora originated mostly from the paleotropical Tertiary flora in the Hengduan Mountains by adapting to the cold and arid environments associated with the strong uplift of the Qinghai-Tibet Plateau (QTP). Here, we combine the phylogeographic history of Sophora moorcroftiana with that of Sophora davidii to explore the speciation of S. moorcroftiana to test this hypothesis. We collected 151 individuals from 17 populations and sequenced 2 chloroplast fragments and the internal transcribed spacer of rDNA. Five chlorotypes and 9 ribotypes were detected but no significant phylogeographic structure was revealed. The integrated results of phylogeographic studies of these 2 species clearly support the progenitor-derivative relationship between them. We infer that the western peripheral population of S. davidii migrated westwards from the Hengduan Mountains to the middle reaches of the Yarlung Zangbo River and differentiated from its ancestor in the process of adaptation to increasingly cold and arid environments with the uplift of the QTP and finally evolved into S. moorcroftiana during the Late Pliocene. In addition, our findings shed light on the idea that natural selection, as imposed by climate differentiation (especially mean diurnal range and precipitation seasonality), directly drove this peripatric speciation event after geographic isolation. The speciation of S. moorcroftiana is a strong case supporting Wu's hypothesis about the origin of Tibet's flora.

Was the universal common ancestry proved?

The question of whether or not all life on Earth shares a single common ancestor has been a central problem of evolutionary biology since Darwin. Although the theory of universal common ancestry (UCA) has gathered a compelling list of circumstantial evidence, as given in ref. 2, there has been no attempt to test statistically the UCA hypothesis among the three domains of life (eubacteria, archaebacteria and eukaryotes) by using molecular sequences. Theobald recently challenged this problem with a formal statistical test, and concluded that the UCA hypothesis holds. Although his attempt is the first step towards establishing the UCA theory with a solid statistical basis, we think that the test of Theobald is not sufficient enough to reject the alternative hypothesis of the separate origins of life, despite the Akaike information criterion (AIC) of model selection giving a clear distinction between the competing hypotheses.

Compression-only CPR training in elementary schools and student attitude toward CPR.

BACKGROUND: Little is known about the effectiveness of systematic cardiopulmonary resuscitation (CPR) training for elementary school children. METHODS: We introduced systematic training of chest compression-only CPR and automated external defibrillator (AED) use to elementary school students aged 10-12 years at 17 schools. The questionnaire compared student attitudes towards CPR and their knowledge about it before and after CPR training. We also evaluated parent and teacher views about CPR training in school education. The primary outcome was positive attitude, defined as "yes" and "maybe yes" on a 5 point Likert-type scale of student attitudes towards CPR.1 RESULTS: A total of 2047 elementary school students received CPR training. Of them, 1899 (92.8%) responded to the questionnaire regarding their attitude towards CPR before and after the training. Before training, 50.2% answered "yes" and 30.3% answered "maybe yes", to the question: "If someone suddenly collapses in front of you, can you do something such as check response or call emergency?" After training, their answers changed to 75.6% and 18.3% for "yes" and "maybe yes", respectively. Many of the students (72.3%, 271/370) who did not have a positive attitude before CPR training had a positive attitude after the training (P < 0.001). Most students understood how to perform CPR (97.7%) and use an AED (98.5%). Parents (96.2%, 1173/1220) and teachers (98.3%, 56/57) answered that it was "good" and "maybe good" for children to receive the training at elementary schools. CONCLUSION: Systematic chest compression-only CPR training helped elementary school students to improve their attitude towards CPR.

Polymorphism and evolution of ribosomal DNA in tea (Camellia sinensis, Theaceae).

Ribosomal DNA (rDNA) repeats often exhibit a high level of within-species homogeneity, but intra-individual polymorphism of rDNA has been found in Camellia species. In order to reveal the evolutionary pattern of rDNA repeats in tea [Camellia sinensis (L.) O. Kuntze], we identified the 45S rDNA loci, estimated their copy number, and cloned partial regions of them from different PCR products and from digested genomic DNA. The results show that there are 3 loci (6 sites) containing 45S rDNA in the tea genome; they are located at the ends of the short arms of 6 chromosomes and consist of a maximum of ∼6500 repeat units. On one hand, 164 sequences of the 26S rDNA cloned from PCR products contained 147 haplotypes, consisting of 62% pseudogenes, 24% putative functional genes, and 14% PCR-mediated recombinants, suggesting that the 45S rDNA of tea maintains an extremely high level of polymorphism and divergence; on the other hand, rDNA fragments cloned directly from genomic DNA exhibited a very high level of homogeneity: only one of 33 rDNA fragments was from a pseudogene. These results demonstrate that although over 60% of the 26S rDNA sequences identified in PCR products belong to pseudogenes, most 45S rDNA repeats are functional genes and have undergone concerted evolution. In this study, strong PCR bias and PCR-mediated recombination greatly increased the apparent proportion of pseudogenes in PCR products. Phylogenetic analysis and genetic divergence values for 26S rDNA sequences obtained in this study show that many pseudogenes have originated independently from functional genes at different times, and despite thus escaping from concerted evolution, they have failed to be eliminated from the tea genome over a long period, some of them having even produced addition copies by rapid expansion. Importantly, our study suggests that in order to determine the true pattern of evolution of rDNA it is necessary to combine data from more than one method rather than relying only on sequences from PCR products.

Transcriptome profiling of the UV-B stress response in the desert shrub Lycium ruthenicum.

Ultraviolet-B (UV-B) is a natural component of the solar radiation. Due to its high energy, low dosages of UV-B can bring huge potential damage effect to organisms. Despite much research that has analyzed the gene expression changes of plants that under UV-B radiation, the transcriptome response of Lycium ruthenicum under the UV-B induction is still un-available. The aim of our study was to identify UV-B responsive genes and gain an insight into the underlying genetic basis of the pathobiology of UV-B related damage. We collected leaf samples from L. ruthenicum with and without UV-B exposure, and then performed a transcriptome profiling to comprehensively investigate their expression signatures. By employing the high throughput RNA-sequencing analysis of samples with and without UV-B exposure, we identified 1,913 up-regulated and 536 down-regulated genes at least by twofold changes. The activity of antioxidant enzyme related genes, including the superoxide dismutase, was decreased, genes related to the synthesis of secondary metabolites and defense responses, such as cinnamyl alcohol dehydrogenase, chalcone-flavanone isomerase and dihydroflavonol reductase were also downregulated. The expression patterns of 14 randomly selected genes resulted from quantitative real-time PCR were basically consistent with their transcript abundance changes identified by RNA-sequencing. We found that several biological pathways related to biotic and abiotic stresses, including cell defense, photosynthesis processes, energy metabolism, were involved in the process of UV-B stress response. A genome-wide screening of gene deregulation under UV-B induction would provide an insight into the understanding of the molecular bases and pathogenesis of UV-B responses.

Origin and spatial population structure of Malagasy native chickens based on mitochondrial DNA.

Since Malagasy human culture became established in a multi-layered way by genetic admixture of Austronesian (Indonesia), Bantu (East Africa) and West Asian populations, the Malagasy native livestock should also have originated from these regions. While recent genetic studies revealed that Malagasy native dogs and goats were propagated from Africa, the origin of Malagasy native chickens is still controversial. Here, we conducted a phylogeographic analysis of the native chickens, focusing on the historical relationships among the Indian Ocean rim countries and based on mitochondrial D-loop sequences. Although previous work suggested that the rare Haplogroup D occurs with high frequencies in Island Southeast Asia-Pacific, East Africa and Madagascar, the major mitochondrial lineage in Malagasy populations is actually not Haplogroup D but the Sub-haplogroup C2, which is also observed in East Africa, North Africa, India and West Asia. We demonstrate that the Malagasy native chickens were propagated directly from West Asia (including India and North Africa), and not via East Africa. Furthermore, they display clear genetic differentiation within Madagascar, separated into the Highland and Lowland regions as seen in the human genomic landscape on this island. Our findings provide new insights for better understanding the intercommunion of material/non-material cultures within and around Madagascar.

Phylogenetic position and evolutionary history of the turtle and whale barnacles (Cirripedia: Balanomorpha: Coronuloidea).

Barnacles of the superfamily Coronuloidea are obligate epibionts of various marine mammals, marine reptiles and large crustaceans. We used five molecular markers: 12S rDNA, 16S rDNA, 18S rDNA, 28S rDNA and Histone 3 to infer phylogenetic relationships among sixteen coronuloids, representing most of the recent genera of barnacles of this superfamily. Our analyses confirm the monophyly of Coronuloidea and that this superfamily and Tetraclitoidea are sister groups. The six-plated Austrobalanus clusters with these two superfamilies. Based on BEAST and ML trees, Austrobalanus is basal and sister to the Coronuloidea, but the NJ tree places Austrobalanus within the Tetraclitoidae, and in the MP tree it is sister to both Coronuloidea and Tetraclitoidae. Hence the position of Austrobalanus remains unresolved. Within the Coronuloidea we identified four clades. Chelonibia occupies a basal position within the Coronuloidea which is in agreement with previous studies. The grouping of the other clades does not conform to previous studies. Divergence time analyses show that some of the time estimates are congruent with the fossil record while some others are older, suggesting the possibility of gaps in the fossil record.

Evolution of snow algae, from cosmopolitans to endemics, revealed by DNA analysis of ancient ice.

Recent studies of microbial biogeography have revealed the global distribution of cosmopolitans and dispersal of regional endemics, but little is known about how these processes are affected by microbial evolution. Here, we compared DNA sequences from snow/glacier algae found in an 8000-year-old ice from a glacier in central Asia with those from modern snow samples collected at 34 snow samples from globally distributed sites at the poles and mid-latitudes, to determine the evolutionary relationship between cosmopolitan and endemic phylotypes of snow algae. We further applied a coalescent theory-based demographic model to the DNA sequences. We found that the genus Raphidonema (Trebouxiophyceae) was distributed over both poles and mid-latitude regions and was detected in different ice core layers, corresponding to distinct time periods. Our results indicate that the modern cosmopolitan phylotypes belonging to Raphidonema were persistently present long before the last glacial period. Furthermore, endemic phylotypes originated from ancestral cosmopolitan phylotypes, suggesting that modern regional diversity of snow algae in the cryosphere is a product of microevolution. These findings suggest that the cosmopolitans dispersed across the world and then derived new localized endemics, which thus improves our understanding of microbial community formation by microevolution in natural environments.

Climate heterogeneity shapes phylogeographic pattern of Hippophae gyantsensis (Elaeagnaceae) in the east Himalaya-Hengduan Mountains.

The interaction of recent orographic uplift and climate heterogeneity acted as a key role in the East Himalaya-Hengduan Mountains (EHHM) has been reported in many studies. However, how exactly the interaction promotes clade diversification remains poorly understood. In this study, we both used the chloroplast trnT-trnF region and 11 nuclear microsatellite loci to investigate the phylogeographic structure and population dynamics of Hippophae gyantsensis and estimate what role geological barriers or ecological factors play in the spatial genetic structure. The results showed that this species had a strong east-west phylogeographic structure, with several mixed populations identified from microsatellite data in central location. The intraspecies divergence time was estimated to be about 3.59 Ma, corresponding well with the recent uplift of the Tibetan Plateau. Between the two lineages, there was significant climatic differentiation without geographic barriers. High consistency between lineage divergence, climatic heterogeneity, and Qingzang Movement demonstrated that climatic heterogeneity but not geographic isolation drives the divergence of H. gyantsensis, and the recent regional uplift of the QTP, as the Himalayas, creates heterogeneous climates by affecting the flow of the Indian monsoon. The east group of H. gyantsensis experienced population expansion c. 0.12 Ma, closely associated with the last interglacial interval. Subsequently, a genetic admixture event between east and west groups happened at 26.90 ka, a period corresponding to the warm inter-glaciation again. These findings highlight the importance of the Quaternary climatic fluctuations in the recent evolutionary history of H. gyantsensis. Our study will improve the understanding of the history and mechanisms of biodiversity accumulation in the EHHM region.

Domestication relaxed selective constraints on the yak mitochondrial genome.

Wild yaks roaming in the high-altitude Qinghai-Tibetan plateau have to maintain high metabolic efficiency. However, domestic yaks do not require such high efficiency because of their limited activity. Hence, domestication may have caused the relaxation of selective constraints on the yak mitochondrial genome because mitochondrial mutations are extremely sensitive to energy-related selective pressures. We have tested this hypothesis by analyzing the mitochondrial genomes of 51 domestic yaks and 21 wild yaks. The results show that the ratio of nonsynonymous/synonymous substitutions in mitochondrial protein-coding genes is significantly higher in domestic yak lineages than those of wild yaks. This genetic difference suggests that the relaxation of selective constraints following the domestication in addition to bottlenecks has allowed faster accumulation of nonsilent substitutions in the yak mitochondrial genome, despite its short domestication history.

Some problems in proving the existence of the universal common ancestor of life on Earth.

Although overwhelming circumstantial evidence supports the existence of the universal common ancestor of all extant life on Earth, it is still an open question whether the universal common ancestor existed or not. Theobald (Nature 465, 219-222 (2010)) recently challenged this problem with a formal statistical test applied to aligned sequences of conservative proteins sampled from all domains of life and concluded that the universal common ancestor hypothesis holds. However, we point out that there is a fundamental flaw in Theobald's method which used aligned sequences. We show that the alignment gives a strong bias for the common ancestor hypothesis, and we provide an example that Theobald's method supports a common ancestor hypothesis for two apparently unrelated families of protein-encoding sequences (cytb and nd2 of mitochondria). This arouses suspicion about the effectiveness of the "formal" test.

Molecular Evolutionary Rate Predicts Intraspecific Genetic Polymorphism and Species-Specific Selection.

It is unknown what determines genetic diversity and how genetic diversity is associated with various biological traits. In this work, we provide insight into these issues. By comparing genetic variation of 14,671 mammalian gene trees with thousands of individual human, chimpanzee, gorilla, mouse, and dog/wolf genomes, we found that intraspecific genetic diversity can be predicted by long-term molecular evolutionary rates rather than de novo mutation rates. This relationship was established during the early stage of mammalian evolution. Moreover, we developed a method to detect fluctuations of species-specific selection on genes based on the deviations of intraspecific genetic diversity predicted from long-term rates. We showed that the evolution of epithelial cells, rather than connective tissue, mainly contributed to morphological evolution of different species. For humans, evolution of the immune system and selective sweeps caused by infectious diseases are the most representative examples of adaptive evolution.

Complete Mitochondrial Genome Analysis Clarifies the Enigmatic Origin of Haplogroup D in Japanese Native Chickens.

Japanese native chickens (JNCs) comprise approximately 50 breeds, making Japan a diversity hotspot for native chicken breeds. JNCs were established through the repeated introduction of chickens from foreign countries. Jidori, which is the generic name of JNC breeds whose ancestral morphology resembles that of their wild progenitor (red junglefowls), is generally thought to have propagated from north East Asia (Korea and north China) to ancient Japan. However, mitochondrial haplogroup D, which is abundant in Island Southeast Asia (ISEA) as well as the Pacific but relatively rare in other regions, can be observed in some Jidori breeds (e.g., Tosa-Jidori, Tokuji-Jidori) with high frequency, leading to speculation that chickens from ISEA or the Pacific also contributed genetically to JNCs. To test this hypothesis, we sequenced the mitochondrial genomes of Jidori breeds and conducted phylogeographic analysis. Our results indicate that the JNC Haplogroup D belongs to Sub-haplogroup D2, which is currently only observed in Xinjiang, northwest China, and not to Sub-haplogroup D1, which is widely distributed in the ISEA-Pacific region. The other mitochondrial haplogroups of Jidori examined in this study also showed affinity to those of chickens native to north East Asia. Therefore, our findings support the north East Asian origin hypothesis for Jidori.