Integrating fossils in a molecular‐based phylogeny and testing them as calibration points for divergence time estimates in Menispermaceae

Abstract The phylogeny of extant Menispermaceae (Ranunculales) is reconstructed based on DNA sequences of two chloroplast genes (rbcL and atpB) from 94 species belonging to 56 genera. Fossilized endocarps represent 34 genera. The positions of these are inferred using 30 morphological characters and the molecular phylogeny as a backbone constraint. Nine of the thirteen nodes that are each dated by a fossil are used as calibration points for the estimates of molecular divergence times. BEAST is used to estimate stem age (121.2 Myr) and crown age (105.4 Myr) for Menispermaceae. This method does not require an input tree topology and can also account for rate heterogeneity among lineages. The sensitivity of these estimates to fossil constraints is then evaluated by a cross‐validation procedure. The estimated origin for Menispermaceae is dated to the mid‐Jurassic if the customary maximum age of 125 Myr for eudicots is not implemented. All constraints when used alone failed to estimate node ages in some parts of the tree. Fossils from the Palaeocene and Eocene impose strict constraints. Likewise, the use of Prototinomiscium as a dating constraint for Menispermaceae appears to be a conservative approach.     

Integrating fossils in a molecular-based phylogeny and testing them as calibration points for divergence time estimates in Menispermaceae

The phylogeny of extant Menispermaceae (Ranunculales) is reconstructed based on DNA sequences of two chloroplast genes (rbcL and atpB) from 94 species belonging to 56 genera. Fossilized endocarps represent 34 genera. The positions of these are inferred using 30 morphological characters and the molecular phylogeny as a backbone constraint. Nine of the thirteen nodes that are each dated by a fossil are used as calibration points for the estimates of molecular divergence times. BEAST is used to estimate stem a...     

Revisiting the phylogeny of Ranunculeae: Implications for divergence time estimation and historical biogeography

Previous studies based on different molecular datasets have generated conflicting topologies for Ranunculeae. Here, we revisit the phylogeny of Ranunculeae by analyzing the individual matK/trnK, psbJ-petA, and internal transcribed spacer (ITS) data, the combined matK/trnK, psbJ-petA, and ITS dataset, and the combinedrbcL, trnL-F, matK/trnK, psbJ-petA, and ITS dataset. Based on the tree-based comparisons, blast searches against NCBI of the sequences, and close examination of the alignment, we found that 10 psbJ-petA sequences previously used were questionable (erroneous or problematic) and responsible for previous conflicting topologies. After omitting these questionable sequences, we provide a new phylogeny for Ranunculeae, in which Beckwithia–Cyrtorhyncha, Kumlienia, andPeltocalathos were replaced. These new replacements are supported by corresponding morphological characters. Moreover, three previously proposed intercontinental disjunct distributions within Ranunculus were also refuted. In our framework, our divergence time and biogeographic analyses indicate that divergence time estimates and the ancestral areas reconstructed for 10 of the 15 nodes in the genus-level phylogeny were influenced by elimination of the questionable sequences. The most recent common ancestor of Ranunculeae was inferred to be present in Europe and North America during the late Eocene. Clades I and II began to diversify in Europe and North America, respectively, and subsequently migrated to other continents. This study shows that it is necessary to analyze individual chloroplast DNA region datasets separately to detect questionable sequences early in the study. The combined dataset including the questionable sequences resulted in an erroneous phylogenetic tree, and the use of this tree subsequently affected age estimates and biogeographic analyses.     

Genomic divergence between the migratory and stationary ecotypes of Atlantic cod

Atlantic cod displays a range of phenotypic and genotypic variations, which includes the differentiation into coastal stationary and offshore migratory types of cod that co‐occur in several parts of its distribution range and are often sympatric on the spawning grounds. Differentiation of these ecotypes may involve both historical separation and adaptation to ecologically distinct environments, the genetic basis of which is now beginning to be unravelled. Genomic analyses based on recent sequencing advances are able to document genomic divergence in more detail and may facilitate the exploration of causes and consequences of genome‐wide patterns. We examined genomic divergence between the stationary and migratory types of cod in the Northeast Atlantic, using next‐generation sequencing of pooled DNA from each of two population samples. Sequence data was mapped to the published cod genome sequence, arranged in more than 6000 scaffolds (611 Mb). We identified 25 divergent scaffolds (26 Mb) with a higher than average gene density, against a backdrop of overall moderate genomic differentiation. Previous findings of localized genomic divergence in three linkage groups were confirmed, including a large (15 Mb) genomic region, which seems to be uniquely involved in the divergence of migratory and stationary cod. The results of the pooled sequencing approach support and extend recent findings based on single‐nucleotide polymorphism markers and suggest a high degree of reproductive isolation between stationary and migratory cod in the North‐east Atlantic.     

The impact of fossil stratigraphic ranges on tip-calibration,and the accuracy and precision of divergence time estimates

The molecular clock provides the only viable means of establishing realistic evolutionary timescales but it remains unclear how best to calibrate divergence time analyses. Calibrations can be applied to the tips and/or to the nodes of a phylogeny. Tip-calibration is an attractive approach since it allows fossil species to be included alongside extant relatives in molecular clock analyses. However, most fossil species are known from multiple stratigraphical horizons and it remains unclear how such age ranges should be interpreted to codify tip-calibrations. We use simulations and empirical data to explore the impact on precision and accuracy of different approaches to informing tip-calibrations. In particular, we focus on the effect of using tip-calibrations defined using the oldest vs youngest stratigraphic occurrences, the full stratigraphical range, as well as confidence intervals on these data points. The results of our simulations show that using different calibration approaches leads to different divergence-time estimates and demonstrate that concentrating tip-calibrations near the root of the dated phylogeny improves both precision and accuracy of estimated divergence times. Finally, our results indicate that the highest levels of accuracy and precision are achieved when fossil tips are calibrated based on the fossil occurrence from which the morphological data were derived. These trends were corroborated by analysis of an empirical dataset for Ursidae. Overall, we conclude that tip-dating analyses should, in particular, employ tip calibrations close to the root of the tree and they should be calibrated based on the age of the fossil used to inform the morphological data used in Total Evidence Dating.     

Rapid,economical single‐nucleotide polymorphism and microsatellite discovery based on de novo assembly of a reduced representation genome in a non‐model organism: a case study of Atlantic cod Gadus morhua

By combining next‐generation sequencing technology (454) and reduced representation library (RRL) construction, the rapid and economical isolation of over 25 000 potential single‐nucleotide polymorphisms (SNP) and >6000 putative microsatellite loci from c. 2% of the genome of the non‐model teleost, Atlantic cod Gadus morhua from the Celtic Sea, south of Ireland, was demonstrated. A small‐scale validation of markers indicated that 80% (11 of 14) of SNP loci and 40% (6 of 15) of the microsatellite loci could be amplified and showed variability. The results clearly show that small‐scale next‐generation sequencing of RRL genomes is an economical and rapid approach for simultaneous SNP and microsatellite discovery that is applicable to any species. The low cost and relatively small investment in time allows for positive exploitation of ascertainment bias to design markers applicable to specific populations and study questions.     

Resequencing and comparison of whole mitochondrial genome to gain insight into the evolutionary status of the Shennongjia golden snub‐nosed monkey (SNJ R. roxellana)

Shennongjia Rhinopithecus roxellana (SNJ R. roxellana) is the smallest geographical population of R. roxellana. The phylogenetic relationships among its genera and species and the biogeographic processes leading to their current distribution are largely unclear. To address these issues, we resequenced and obtained a new, complete mitochondrial genome of SNJ R. roxellana by next‐generation sequencing and standard Sanger sequencing. We analyzed the gene composition, constructed a phylogenetic tree, inferred the divergence ages based on complete mitochondrial genome sequences, and analyzed the genetic divergence of 13 functional mtDNA genes. The phylogenetic tree and divergence ages showed that R. avunculus (the Tonkin snub‐nosed monkey) was the first to diverge from the Rhinopithecus genus ca. 2.47 million years ago (Ma). Rhinopithecus bieti and Rhinopithecus strykeri formed sister groups, and the second divergence from the Rhinopithecus genus occurred ca. 1.90 Ma. R. roxellana and R. brelichi diverged from the Rhinopithecus genus third, ca. 1.57 Ma. SNJ R. roxellana was the last to diverge within R. roxellana species in 0.08 Ma, and the most recent common ancestor of R. roxellana is 0.10 Ma. The analyses on gene composition showed SNJ R. roxellana was the newest geographic population of R. roxellana. The work will help to develop a more accurate protection policy for SNJ R. roxellana and facilitate further research on selection and adaptation of R. roxellana.     

The distribution of the coalescence time and the number of pairwise nucleotide differences in a model of population divergence or speciation with an initial period of gene flow

This paper is concerned with a model of “isolation with an initial period of migration”, where a panmictic ancestral population split into n descendant populations which exchanged migrants symmetrically at a constant rate for a period of time and subsequently became completely isolated. In the limit as the population split occurred an infinitely long time ago, the model becomes an “isolation after migration” model, describing completely isolated descendant populations which arose from a subdivided ancestral population. The probability density function of the coalescence time of a pair of genes and the probability distribution of the number of pairwise nucleotide differences are derived for both models. Whilst these are theoretical results of interest in their own right, they also give an exact analytical expression for the likelihood, for data consisting of the numbers of nucleotide differences between pairs of DNA sequences where each pair is at a different, independent locus. The behaviour of the distribution of the number of pairwise nucleotide differences under these models is illustrated and compared to the corresponding distributions under the “isolation with migration” and “complete isolation” models. It is shown that the distribution of the number of nucleotide differences between a pair of DNA sequences from different descendant populations in the model of “isolation with an initial period of migration” can be quite different from that under the “isolation with migration model”, even if the average migration rate over time (and hence the total number of migrants) is the same in both scenarios. It is also illustrated how the results can be extended to other demographic scenarios that can be described by a combination of isolated panmictic populations and “symmetric island” models.     

Hybridisation‐based target enrichment of phenology genes to dissect the genetic basis of yield and adaptation in barley

Barley (Hordeum vulgare L.) is a major cereal grain widely used for livestock feed, brewing malts and human food. Grain yield is the most important breeding target for genetic improvement and largely depends on optimal timing of flowering. Little is known about the allelic diversity of genes that underlie flowering time in domesticated barley, the genetic changes that have occurred during breeding, and their impact on yield and adaptation. Here, we report a comprehensive genomic assessment of a worldwide collection of 895 barley accessions based on the targeted resequencing of phenology genes. A versatile target‐capture method was used to detect genome‐wide polymorphisms in a panel of 174 flowering time‐related genes, chosen based on prior knowledge from barley, rice and Arabidopsis thaliana. Association studies identified novel polymorphisms that accounted for observed phenotypic variation in phenology and grain yield, and explained improvements in adaptation as a result of historical breeding of Australian barley cultivars. We found that 50% of genetic variants associated with grain yield, and 67% of the plant height variation was also associated with phenology. The precise identification of favourable alleles provides a genomic basis to improve barley yield traits and to enhance adaptation for specific production areas.     

A Resurgence in Field Research is Essential to Better Understand the Diversity,Ecology, and Evolution of Microbial Eukaryotes

The discovery and characterization of protist communities from diverse environments are crucial for understanding the overall evolutionary history of life on earth. However, major questions about the diversity, ecology, and evolutionary history of protists remain unanswered, notably because data obtained from natural protist communities, especially of heterotrophic species, remain limited. In this review, we discuss the challenges associated with “field protistology”, defined here as the exploration, characterization, and interpretation of microbial eukaryotic diversity within the context of natural environments or field experiments, and provide suggestions to help fill this important gap in knowledge. We also argue that increased efforts in field studies that combine molecular and microscopical methods offer the most promising path toward (1) the discovery of new lineages that expand the tree of eukaryotes; (2) the recognition of novel evolutionary patterns and processes; (3) the untangling of ecological interactions and functions, and their roles in larger ecosystem processes; and (4) the evaluation of protist adaptations to a changing climate.     

Evaluation of two approaches to genotyping major histocompatibility complex class I in a passerine-CE-SSCP and 454 pyrosequencing

Genes of the highly dynamic major histocompatibility complex (MHC) are directly linked to individual fitness and are of high interest in evolutionary ecology and conservation genetics. Gene duplication and positive selection usually lead to high levels of polymorphism in the MHC region, making genotyping of MHC a challenging task. Here, we compare the performance of two methods for MHC class I genotyping in a passerine with highly duplicated MHC class I genes: capillary electrophoresis-single-strand conformation polymorphism (CE-SSCP) analysis and 454 GS FLX Titanium pyrosequencing. According to our findings, the number of MHC variants (called alleles for simplicity) detected by CE-SSCP is significantly lower than detected by 454. To resolve discrepancies between the two methods, we cloned and Sanger sequenced a MHC class I amplicon for an individual with high number of alleles. We found a perfect congruence between cloning/Sanger sequencing results and 454. Thus, in case of multi-locus amplification, CE-SSCP considerably underestimates individual MHC diversity. However, numbers of alleles detected by both methods are significantly correlated, although the correlation is weak (r = 0.32). Thus, in systems with highly duplicated MHC, 454 provides more reliable information on individual diversity than CE-SSCP.     

RNA‐sequencing profiles hippocampal gene expression in a validated model of cancer‐induced depression

To investigate the pathophysiology of cancer‐induced depression (CID), we have recently developed a validated CID mouse model. Given that the efficacy of antidepressants in cancer patients is controversial, it remains unclear whether CID is a biologically distinct form of depression. We used RNA‐sequencing (RNA‐seq) to investigate differentially expressed genes (DEGs) in hippocampi of animals from our CID model relative a positive control model of depressive‐like behavior induced with chronic corticosterone (CORT). To validate RNA‐seq results, we performed quantitative real‐time RT‐PCR (qRT‐PCR) on a subset of DEGs. Enrichment analysis using DAVID was performed on DEGs to identify enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and biological process gene ontologies (GO:BP). qRT‐PCR results significantly predicted RNA‐seq results. RNA‐seq revealed that most DEGs identified in the CORT model overlapped with the CID model. Enrichment analyses identified KEGG pathways and GO:BP terms associated with ion homeostasis and neuronal communication for both the CORT and CID model. In addition, CID DEGs were enriched in pathways and terms relating to neuronal development, intracellular signaling, learning and memory. This study is the first to investigate CID at the mRNA level. We have shown that most hippocampal mRNA changes that are associated with a depressive‐like state are also associated with cancer. Several other changes occur at the mRNA level in cancer, suggesting that the CID model may represent a biologically distinct form of a depressive‐like state.