Simulating and detecting autocorrelation of molecular evolutionary rates among lineages

Evolutionary timescales can be estimated from genetic data using phylogenetic methods based on the molecular clock. To account for molecular rate variation among lineages, a number of relaxed‐clock models have been developed. Some of these models assume that rates vary among lineages in an autocorrelated manner, so that closely related species share similar rates. In contrast, uncorrelated relaxed clocks allow all of the branch‐specific rates to be drawn from a single distribution, without assuming any correlation between rates along neighbouring branches. There is uncertainty about which of these two classes of relaxed‐clock models are more appropriate for biological data. We present an R package, NELSI, that allows the evolution of DNA sequences to be simulated according to a range of clock models. Using data generated by this package, we assessed the ability of two Bayesian phylogenetic methods to distinguish among different relaxed‐clock models and to quantify rate variation among lineages. The results of our analyses show that rate autocorrelation is typically difficult to detect, even when there is complete taxon sampling. This provides a potential explanation for past failures to detect rate autocorrelation in a range of data sets.     

Genetic analysis of early holocene skeletal remains from Alaska and its implications for the settlement of the Americas

Mitochondrial and Y-chromosome DNA were analyzed from 10,300-year-old human remains excavated from On Your Knees Cave on Prince of Wales Island, Alaska (Site 49-PET-408). This individual's mitochondrial DNA (mtDNA) represents the founder haplotype of an additional subhaplogroup of haplogroup D that was brought to the Americas, demonstrating that widely held assumptions about the genetic composition of the earliest Americans are incorrect. The amount of diversity that has accumulated in the subhaplogroup over the past 10,300 years suggests that previous calibrations of the mtDNA clock may have underestimated the rate of molecular evolution. If substantiated, the dates of events based on these previous estimates are too old, which may explain the discordance between inferences based on genetic and archaeological evidence regarding the timing of the settlement of the Americas. In addition, this individual's Y-chromosome belongs to haplogroup Q-M3*, placing a minimum date of 10,300 years ago for the emergence of this haplogroup.     

Phylogenetics of early branching eudicots:Comparing phylogenetic signal across plastid introns,spacers,and genes

<正>Recent phylogenetic analyses revealed a grade with Ranunculales,Sabiales,Proteales,Trochodendrales, and Buxales as first branching eudicots,with the respective positions of Proteales and Sabiales still lacking statistical confidence.As previous analyses of conserved plastid genes remain inconclusive,we aimed to use and evaluate a representative set of plastid introns(group I:trnL;group II:petD,rpll6,trnK) and intergenic spacers(trnL-F,petB-petD, atpB-rbcL,rps3-rpl16) in comparison to the rapidly evolving matK and slowly evolving atpB and rbcL genes. Overall patterns of micro structural mutations converged across genomic regions,underscoring the existence of a general mutational pattern throughout the plastid genome.Phylogenetic signal differed strongly between functionally and structurally different genomic regions and was highest in matK,followed by spacers,then group II and group I introns.The more conserved atpB and rbcL coding regions showed distinctly lower phylogenetic information content.Parsimony,maximum likelihood,and Bayesian phylogenetic analyses based on the combined dataset of non-coding and rapidly evolving regions(14 000 aligned characters) converged to a backbone topology of eudicots with Ranunculales branching first,a Proteales-Sabiales clade second,followed by Trochodendrales and Buxales. Gunnerales generally appeared as sister to all remaining core eudicots with maximum support.Our results show that a small number of intron and spacer sequences allow similar insights into phylogenetic relationships of eudicots compared to datasets of many combined genes.The non-coding proportion of the plastid genome thus can be considered an important information source for plastid phylogenomics.     

Molecular analysis of phylogenetic relationships among Myrmecophytic macaranga species (Euphorbiaceae)

Many species of the paleotropical pioneer tree genus Macaranga Thou. (Euphorbiaceae) live in association with ants. Various types of mutualistic interactions exist, ranging from the attraction of unspecific ant visitors to obligate myrmecophytism. In the latter, nesting space and food bodies are exchanged for protection by highly specific ant partners (mainly species of the myrmicine genus Crematogaster). As a first step toward elucidating the coevolution of ant-plant interactions in the Macaranga-Crematogaster system, we have initiated a molecular investigation of the plant partners' phylogeny. Nuclear ribosomal DNA internal transcribed spacer (ITS) sequences were analyzed for 73 accessions from 47 Macaranga species, representing 17 sections or informally described species groups. Three accessions from the putative sister taxon Mallotus Lour, were included as outgroups. Cladograms of the ITS data revealed Macaranga to be nested within Mallotus. ITS sequences are highly similar within section Pachystemon s.str., suggesting a relatively recent and rapid radiation of obligate myrmecophytes within this section. Forty-three accessions, mainly of ant-inhabited species, were additionally investigated by random amplified polymorphic DNA (RAPD) and microsatellite-primed PCR (MP-PCR) techniques. Phenetic analysis of RAPD and MP-PCR banding profiles generally confirmed the ITS results. Best resolutions for individual clades were obtained when ITS and RAPD/MP-PCR data were combined into a single matrix and analyzed phenetically. The combined analysis suggests multiple (four) rather than a single evolutionary origin of myrmecophytism, at least one reversal from obligate myrmecophytism to nonmyrmecophytism, and one loss of mutualistic specifity.     

Origin and diversification of Thymelaea (Thymelaeaceae): inferences from a phylogenetic study based on ITS (rDNA) sequences

Thymelaea is a Mediterranean genus belonging to a primarily tropical and subtropical family. This genus is here presented as a particular case on which the hypothesis of an in situ evolution of the Mediterranean flora from a Tertiary subtropical stock can be phylogenetically tested. To better understand the evolutionary history of Thymelaea, molecular phylogenies based on ITS (rDNA) sequences are estimated. Phylogenetic analyses strongly support the monophyly of Thymelaea and the sister position of Daphne. Unexpected relationships between sequences of T. granatensis, T. pubescens and T. sanamunda are found and explained as a case of cryptic nuclear introgression. Phylogenetic analyses also reveal the polyphyletic character of Diarthron, sensu Tan, suggesting that Dendrostellera and Stelleropsis should be reinstated as independent genera. Based on the resulting phylogenetic hypothesis and using statistical tests it is suggested that Thymelaea experienced an initial rapid diversification. Time estimates for the origin and radiation of Thymelaea based on independent sources of evidence—paleontological, paleoclimatic and paleogeographic data, on one hand, and ITS substitution rates reported for other angiosperm groups, on the other—give strikingly similar results. The early radiation of Thymelaea is inferred to have occurred in the Upper Miocene, after this lineage diverged from that of Daphne in the Lower or Middle Miocene. Both evolutionary events are related to the overall climatic deterioration (declining temperature and increasing dryness) registered throughout the Miocene in Eurasia. The rapid and extensive diversification of Thymelaea is placed in the Iberian Peninsula or here and in the neighbouring North African territory.     

A time‐calibrated phylogenetic approach to assessing the phylogeography,colonization history and phenotypic evolution of snakes in the Japanese Izu Islands

Aim We infer the biogeography and colonization history of a dispersal‐limited terrestrial vertebrate, the Japanese four‐lined ratsnake (Elaphe quadrivirgata), to reveal the number of times mainland populations have invaded the Izu Archipelago of Japan, the mainland sources of these colonists, and the time‐scale of colonization. We compare these results with those of past studies in an attempt to uncover general biogeographical patterns. Moreover, we briefly examine the significance of colonization history when evaluating the evolution of body size and melanism of the Izu Island E. quadrivirgata populations. Location The Izu Islands (Oshima, Toshima, Niijima, Shikine, Kozu, Tadanae and Mikura), a volcanic archipelago off the Pacific coast of central Japan. Methods We obtained DNA sequences of the mitochondrial cytochrome b gene (1117 base pairs) from 373 individual snakes sampled from seven of the Izu Islands and 25 mainland localities. We employed partitioned Bayesian phylogenetic analyses assuming a relaxed molecular clock to estimate phylogenetic relationships among extant haplotypes and to give an explicit temporal scale to the timing of clade divergence, colonization history and tempo of body‐size evolution. Moreover, we employed model‐based biogeographical analysis to calculate the minimum number of times E. quadrivirgata colonized the Izu Islands. Results We found evidence that three separate regions of the Izu Archipelago have been colonized independently from mainland ancestors within the past 0.58–0.20 Ma. The Izu Peninsula plus Oshima and Mikura were both colonized independently from lineages inhabiting eastern mainland Japan. The Toshima, Niijima, Shikine, Kozu and Tadanae populations all derive from a single colonization from western mainland Japan. Oshima has been subject to three or possibly four colonizations. Main conclusions These results support the hypothesis that the extreme body‐size disparity among island populations of this ratsnake evolved in situ. Moreover, the fact that the dwarf, melanistic population inhabiting Oshima descends from multiple mainland colonization events is evidence of an extremely strong natural selection pressure resulting in the rapid evolution of this unique morphology. These results contrast with theoretical predictions that natural selection pressures should play a decreased role on islands close to the mainland and/or subject to frequent or recent immigration.     

Phylogenetic relationships and historical biogeography of tribes and genera in the subfamily Nymphalinae (Lepidoptera: Nymphalidae)

We infer for the first time the phylogenetic relationships of genera and tribes in the ecologically and evolutionarily well‐studied subfamily Nymphalinae using DNA sequence data from three genes: 1450 bp of cytochrome oxidase subunit I (COI) (in the mitochondrial genome), 1077 bp of elongation factor 1‐alpha (EF1‐α) and 400–403 bp of wingless (both in the nuclear genome). We explore the influence of each gene region on the support given to each node of the most parsimonious tree derived from a combined analysis of all three genes using Partitioned Bremer Support. We also explore the influence of assuming equal weights for all characters in the combined analysis by investigating the stability of clades to different transition/transversion weighting schemes. We find many strongly supported and stable clades in the Nymphalinae. We are also able to identify ‘rogue’ taxa whose positions are weakly supported (the different gene regions are in conflict with each other) and unstable. Our main conclusions are: (1) the tribe Coeini as currently constituted is untenable, and Smyrna, Colobura and Tigridia are part of Nymphalini; (2) ‘Kallimini’ is paraphyletic with regard to Melitaeini and should be split into three tribes: Kallimini s.s., Junoniini and Victorinini; (3) Junoniini, Victorinini, Melitaeini and the newly circumscribed Nymphalini are strongly supported monophyletic groups, and (4) Precis and Junonia are not synonymous or even sister groups. The species Junonia coenia, a model system in developmental biology, clearly belongs in the genus Junonia. A dispersal‐vicariance analysis suggests that dispersal has had a major effect on the distributions of extant species, and three biotic regions are identified as being centres of diversification of three major clades: the Palaearctic for the Nymphalis‐group, the Afrotropics for Junoniini and the Nearctic for Melitaeini. © 2005 The Linnean Society of London, Biological Journal of the Linnean Society, 2005, 86 , 227–251.     

Molecular phylogeny of an endemic radiation of Cuban toads (Bufonidae: Peltophryne) based on mitochondrial and nuclear genes

Aim In this study we present a molecular phylogenetic and phylogeographical analysis of Peltophryne (Anura: Bufonidae), an endemic genus of Antillean toads, to investigate the spatial and temporal origins of the genus, with particular focus on the eight Cuban species. Location Greater Antilles, with extensive sampling of the Cuban archipelago. Methods We obtained DNA sequence data from two mitochondrial genes, cytochrome c oxidase subunit I (COI) and ribosomal RNA (16S), for 124 toads representing all eight Cuban species, and combined this with published data from Hispaniola (one of three species) and Puerto Rico (one of one species) to establish a molecular phylogeny for Peltophryne. In addition, we explored the phylogeographical structure of widespread Cuban species. For a subset of 42 toads we also obtained DNA sequence data from two nuclear genes, recombination activator‐1 (RAG‐1) and chemokine receptor 4 (CXCR‐4). We combined our molecular data with published DNA sequences from a global sample of bufonid toads to place the spatial and temporal origins of Peltophryne in the Caribbean within a fuller geographical and phylogenetic context. Results All phylogenetic analyses supported the monophyly of West Indian toads. The ancestor of Peltophyrne diverged from its mainland source around the Eocene–Oligocene boundary, with a subsequent radiation across the Caribbean islands taking place during the Miocene. Cuban species are monophyletic with a basal split in the early–middle Miocene that separates extant small‐bodied from large‐bodied species. Extensive mitochondrial DNA (mtDNA) sampling within widespread Cuban species revealed contrasting phylogeographical patterns. Peltophryne taladai and P. empusa showed deeply divergent lineages, whereas no geographical structure was observed in the widespread P. peltocephala. Main conclusions Our timeline for Peltophryne diversification is consistent with a biogeographical model requiring no long‐distance overwater dispersal. Although confidence intervals on divergence time estimates are wide, the stem age of Peltophyrne coincides with the hypothesized GAARlandia landspan or archipelago, which may have connected South America briefly with the Antilles. The ages of Peltophryne for Puerto Rico, Hispaniola and Cuba are consistent with a recently proposed vicariance scenario for the region. Our molecular results support the recognition of all eight species in Cuba, and provide evidence of possible cryptic species.     

Historical biogeography, phylogenetic relationships and intraspecific diversity of agamid lizards in the Central Asian deserts of Kazakhstan and Uzbekistan

The Central Asian agamid lizards are ecologically and morphologically diverse, occurring across a broad range of desert environments in this biogeographically important region. It is probable that past climatic shifts have significantly influenced the diversification patterns and distributions of the agamid lizards of this region. To assess this within a phylogenetic framework we sequenced a 1200 bp region of mitochondrial DNA and a 1200 bp nuclear gene (RAG-1), incorporating both inter- and intraspecific sampling across Central Asian agamids. Our topology and divergence time estimates support an Eocene origin of the Agaminae subfamily on the Indian subcontinent, coinciding with the collision of India into Eurasia. The onset of aridification in Central Asia during the Late Oligocene, resulting from the retreat of the Paratethys Sea and the intensified uplift of the Tibetan–Himalayan complex, probably played an important role in the diversification of Phrynocephalus, one of the three genera studied. Intensification of aridity and geologic events in the Plio-Pleistocene and Quaternary glacial cycling probably had a significant influence on intraspecific diversification patterns within Phrynocephalus.     

Genetic variation and phylogenetic relationships of a pantropical species group in Polystachya (Orchidaceae)

Amplified fragment length polymorphism (AFLP) markers were used to investigate the relationships among Polystachya accessions from a group of closely related pantropical tetraploids. Before starting with the fingerprinting analyses, the polyploid accessions were first included in a phylogenetic analysis using low‐copy nuclear DNA data to establish their relationships, which confirmed that they belonged to a species group of closely related allotetraploids. Neo‐ and Palaeotropical polyploid accessions formed two hybrid clades with apparently independent origins. Sampling for the AFLP analyses included single accessions from much of the range of the genus and populations from Costa Rica (CR) and Sri Lanka (SL) to compare population structure and genetic diversity in these two areas in more detail. A splits graph of the complete AFLP data showed three major clusters corresponding to three sources of population sampling (P. concreta, SL; P. foliosa, CR; P. masayensis, CR), with individual accessions from Africa and Indian Ocean islands showing a closer relationship to P. concreta from SL than to the two CR species. Individual accessions from the Neotropics occurred in more isolated positions in the splits network, with little resolution. Some P. foliosa accessions clustered with P. masayensis, suggesting some hybridization between the two species, and this was confirmed by Bayesian structure analysis. However, the splits network, structure and analyses of molecular variance indicated a generally high level of genetic divergence between the two CR species, despite their recent hybrid origin, occurrence in largely the same localities and occasional hybridization. Polystachya foliosa from CR had a higher degree of population‐level genetic structure (Φ_ST = 0.291) than P. masayensis from CR (Φ_ST = 0.161) and P. concreta from SL (Φ_ST = 0.138), possibly because of its occurrence within a larger and more environmentally diverse continuous range than the other two species. Genetic divergence between Neo‐ and Palaeotropical members of the pantropical tetraploid group of Polystachya and the nonmonophyly of P. concreta suggested that P. concreta s.l. should be split and the use of this epithet should be confined to the Neotropics (the type is from Martinique). Other names should be used in Africa and the Asian tropics. © 2011 The Linnean Society of London, Botanical Journal of the Linnean Society, 2011, 165 , 235–250.     

Species limits and phylogenetic relationships of red‐finned cryptic species of the seasonal killifish genus Hypsolebias from the Brazilian semi‐arid Caatinga (Teleostei: Cyprinodontiformes: Rivulidae)

Members of the Hypsolebias antenori species group comprise a diverse clade of morphologically similar seasonal killifishes occurring in a vast region of the semi‐arid savannah of northeastern Brazil. The present paper focuses on an assemblage of three allopatric cryptic species (H. antenori from isolated coastal river drainages, Hypsolebias igneus from the São Francisco River basin and Hypsolebias coamazonicus sp. nov. from the Parnaíba River basin) sharing almost identical colour patterns, including the presence of an orangish red anal fin in males, thus herein named as the red‐finned assemblage. A tree‐based approach using mt‐DNA (cytochrome b) supports delimitation of all three species, but indicates that the red‐finned assemblage is paraphyletic – H. igneus and H. coamazonicus are closely related to Hypsolebias nudiorbitatus, whereas H. antenori is the sister group to a clade comprising all 13 species of the H. antenori group included in the analysis. Morphological characters are useful to diagnose species, but are not informative for most clades highly supported by molecular data. H. coamazonicus is distinguished from all other congeners by the possession of a dark grey or black stripe on the dorsal fin in males. The basal position of H. antenori is related to uplift episodes involving the Araripe‐Borborema plateau during the Miocene, which isolated the coastal area inhabited by H. antenori from the remaining areas of the Caatinga. The sister group relationship between H. igneus and H. coamazonicus is attributed to a past connection between the São Francisco and Paranaíba River until the Tertiary.     

Genesis of phenotypic and genotypic diversity in land plants: The present as the key to the past

The evolutionary history of vascular plants is reviewed by extrapolation back through time from a wide range of data recently derived from the present flora, using as the central theme evolutionary inferences gained from phylogenies reconstructed as cladograms. Any region of the genome can be used to infer relationships, but only a combination of knowledge of morphology and the developmental genes that underpin morphology can allow evolutionary interpretation of macroevolutionary transitions; this in turn is necessary to identify bona fide evolutionary radiations and any putative causal key innovations. Such studies require clades to be delimited not by the inclusion of particular extant ‘crown’ species but rather by specific apo‐morphies, thereby giving important phylogenetic roles to extinct as well as extant species. Dating phylogenetic divergences via molecular clocks remains seriously inaccurate, and ultimately relies primarily on fossil benchmarks. First principles suggest that evolution of most regions of the genome is fundamentally gradual, whereas evolution of regions especially prone to strong selection pressure, and of the many facets of the phenotype, is punctuational, being characterized through time dominantly by stasis. Sequence data have proved valuable for inferring monophyletic groups, but within the now widely accepted context of monophyly the taxonomic hierarchy should primarily reflect degrees of morphological rather than molecular divergence. Incongruence among contrasting data sets is best explained by understanding the biological constraints operating on each type of phylogenetic information. The conventional ‘uniformitarian’ view of evolution has only limited applicability as one traces the history of land plants through time. Diversity increased in stepwise fashion, reflecting either attainments of complexity and/or fitness thresholds by the lineage (intrinsic) or the availability of unusually permissive environments, often following major perturbations (extrinsic). The Quaternary period demonstrates especially well the resilience, and ease of migration, of the Earth's vegetation. A higher frequency of generation of novel phenotypes in the deep past is possible, but a far higher frequency of their establishment is certain; together, these factors generate an evolutionary pattern of nested radiations that is fractal, as saturation of the resource space rendered the environment decreasingly permissive through time. In the immediate future, evolutionary‐developmental genetics will have increasing value for testing homology, interpreting homoplasy and elucidating evolutionary constraints, and will become easier to pursue as whole‐genome sequences of additional ‘model’ species further invigorate comparative genomics. Complexity of gene regulation, both by other genes and by the cellular and extra‐cellular environment, appears a particularly fruitful area for further research. Nonetheless, environmental filtering of evolutionary novelties (whether instantaneously isolated mutant ‘prospecies’ or classic neoDarwinian ‘selfish genes’ selectively spreading through panmictic populations) can only be effectively understood by longer term monitoring of populations in the wild, to better capture rare evolutionary and ecological events and to better assess the efficacy of traditional microevolutionary processes. We believe that the resulting renaissance in macroevolutionary studies will encourage a broader systematic perspective ‐ one that better encompasses the remarkable diversity of evolutionary processes that together generated the present diversity of life.     

Identification and development of long non‐coding RNA‐associated regulatory network in colorectal cancer

Colorectal cancer (CRC) is one of the leading causes of cancer‐associated death globally. Long non‐coding RNAs (lncRNAs) have been identified as micro RNA (miRNA) sponges in a competing endogenous RNA (ceRNA) network and are involved in the regulation of mRNA expression. This study aims to construct a lncRNA‐associated ceRNA network and investigate the prognostic biomarkers in CRC. A total of 38 differentially expressed (DE) lncRNAs, 23 DEmiRNAs and 27 DEmRNAs were identified by analysing the expression profiles of CRC obtained from The Cancer Genome Atlas (TCGA). These RNAs were chosen to develop a ceRNA regulatory network of CRC, which comprised 125 edges. Survival analysis showed that four lncRNAs, six miRNAs and five mRNAs were significantly associated with overall survival. A potential regulatory axis of ADAMTS9‐AS2/miR‐32/PHLPP2 was identified from the network. Experimental validation was performed using clinical samples by quantitative real‐time PCR (qRT‐PCR), which showed that expression of the genes in the axis was associated with clinicopathological features and the correlation among them perfectly conformed to the ‘ceRNA theory’. Overexpression of ADAMTS9‐AS2 in colon cancer cell lines significantly inhibited the miR‐32 expression and promoted PHLPP2 expression, while ADAMTS9‐AS2 knockdown had the opposite effects. The constructed novel ceRNA network may provide a comprehensive understanding of the mechanisms of CRC carcinogenesis. The ADAMTS9‐AS2/miR‐32/PHLPP2 regulatory axis may serve as a potential therapeutic target for CRC.