Biases in assessing the evolutionary history of the angiosperm flora of China

In their recent paper published in Nature (2018, 554, 234‐238), Lu et al. use phylogenetic approaches to determine the proportion of the Chinese angiosperm genera that originated during the Miocene or later, and contrast divergence times and phylogenetic dispersion between eastern and western China. One of their key conclusions is that 66% of the angiosperm genera in China originated in the Miocene or later. However, an analysis of 300 angiosperm genera shows that 139 (76.8%) of the 181 genera considered as originating in the Miocene or later by Lu et al. have fossil records before the Miocene. Thus, the evolutionary history of Chinese angiosperm flora has been substantially underestimated in Lu et al. In addition, the results of Lu et al. have been biased by using an incomplete phylogeny.     

Mixed-model reanalysis of primate data suggests tissue and species biases in oligonucleotide-based gene expression profiles

An emerging issue in evolutionary genetics is whether it is possible to use gene expression profiling to identify genes that are associated with morphological, physiological, or behavioral divergence between species and whether these genes have undergone positive selection. Some of these questions were addressed in a recent study (Enard et al. 2002) of the difference in gene expression among human, chimp, and orangutan, which suggested an accelerated rate of divergence in gene expression in the human brain relative to liver. Reanalysis of the Affymetrix data set using analysis of variance methods to quantify the contributions of individuals and species to variation in expression of 12,600 genes indicates that as much as one-quarter of the genome shows divergent expression between primate species at the 5% level. The magnitude of fold change ranges from 1.2-fold up to 8-fold. Similar conclusions apply to reanalysis of Enard et al. 2002 parallel murine data set. However, biases inherent to short oligonucleotide microarray technology may account for some of the tissue and species effects. At high significance levels, more differences were observed in the liver than in the brain in each of the pairwise species comparisons, so it is not clear that expression divergence is accelerated in the human brain. Further, there is an apparent bias toward upregulation of gene expression in the brain in both primates and mice, whereas genes are equally likely to be up- or downregulated in the liver when these species diverge. A small subset of genes that are candidates for adaptive divergence may be identified on the basis of a high ratio of interspecific to intraspecific divergence.     

Prediction of treatment response using gene expression profiles

This paper concerns prediction of clinical outcome from gene expression profiles using work in a different area, nonlinear system identification. In particular, the approach can predict long-term treatment response from data of a landmark article by Golub et al. (Golub, T. R.; Slonim, D. K.; Tamayo, P.; Huard, C.; Gaasenbeek, M.; Mesirov, J. P. et al. Science 1999, 286, 531-537) that has not previously been achieved with these data. The present paper shows that, for these data, gene expression profiles taken at time of diagnosis of acute myeloid leukemia contain information predictive of eventual response to chemotherapy. This was not evident in previous work; indeed, the Golub et al. article did not find a set of genes strongly correlated with clinical outcome. However, the present approach can accurately predict outcome class of gene expression profiles even when the genes do not have large differences in expression levels between the classes.     

Modeling endophilin-mediated Aβ disposal in glioma cells

Autophagy dysregulation has emerged in age-related neurological diseases (Ulland et al.; Matheoud et al.; Ashkenazi et al.). Alzheimer Disease (AD), the most common progressive neurodegenerative disorder, is characterized by the accumulation of amyloid-β (Aβ) plaques caused by aberrant Aβ metabolism (Qiang et al.; Sevigny et al.; Ittner et al.). Glia constitute the brain immune system and ingest extracellular Aβ for degradation via the autophagy-lysosome machinery (Ries and Sastre; Cho et al.). Here, we model the molecular rationale for this clearance process in glioma cells by showing that miR34a inhibits autophagy-mediated disposal of Aβ fibrils and identifying two novel direct targets of miR34a, endophilin-3 and cathepsin B (CTSB, a previously reported enzyme for Aβ degrading (Sun et al.)). Bioinformatics analyses revealed that endophilin-3 expresses at a significantly lower level in neurodegenerative diseases. Its gain-of-function substantially promotes both uptake and degradation of Aβ while small interfering RNA (siRNA)-mediated endophilin-3 knockdown slowed down Aβ clearance and blocked autolysosome formation. Mechanistically, gene ontology (GO) analysis of the endophilin-3 interactome identified by mass spectrometry uncovered enriched components involved in actin binding (with the highest score). Importantly, we validated that the actin-binding protein phostensin interacted with endophilin-3. Phostensin knockdown restored endophilin-3-mediated up-regulation of Aβ clearance. Thus, our findings indicate that miR34a inhibits Aβ clearance by targeting endophilin-3 and CTSB at multiple steps including uptake and autophagy-mediated degradation.     

On Estimating the Relationship between Longitudinal Measurements and Time‐to‐Event Data Using a Simple Two‐Stage Procedure

Summary Ye, Lin, and Taylor (2008, Biometrics 64 , 1238–1246) proposed a joint model for longitudinal measurements and time‐to‐event data in which the longitudinal measurements are modeled with a semiparametric mixed model to allow for the complex patterns in longitudinal biomarker data. They proposed a two‐stage regression calibration approach that is simpler to implement than a joint modeling approach. In the first stage of their approach, the mixed model is fit without regard to the time‐to‐event data. In the second stage, the posterior expectation of an individual's random effects from the mixed‐model are included as covariates in a Cox model. Although Ye et al. (2008) acknowledged that their regression calibration approach may cause a bias due to the problem of informative dropout and measurement error, they argued that the bias is small relative to alternative methods. In this article, we show that this bias may be substantial. We show how to alleviate much of this bias with an alternative regression calibration approach that can be applied for both discrete and continuous time‐to‐event data. Through simulations, the proposed approach is shown to have substantially less bias than the regression calibration approach proposed by Ye et al. (2008) . In agreement with the methodology proposed by Ye et al. (2008) , an advantage of our proposed approach over joint modeling is that it can be implemented with standard statistical software and does not require complex estimation techniques.     

Bayesian analysis of combined chloroplast loci, using multiple calibrations, supports the recent arrival of Melastomataceae in Africa and Madagascar

A new biogeographic scenario for Melastomataceae (Morley and Dick, American Journal of Botany 90(11) pp. 1638-1645, 2003) accepts an ndhF-based phylogeny for the family by Renner et al. (American Journal of Botany 88(7): 1290-1300, 2001), but rejects those authors' divergence time estimates. Morley and Dick concluded that Gondwanan vicariance, rather than the more recent long dispersal proposed by Renner et al. explains the presence of the family in Africa and Madagascar. To assess the strength of this conclusion, a Bayesian analysis was conducted on three times the amount of sequence data used before (ndhF, rbcL, rpl16; 3100 base pairs [bp], excluding all gaps). The Bayesian approach to divergence time estimation does not rely on a strict molecular clock and employs multiple simultaneous minimal or maximal bounds on node ages. Reliance on northern mid-latitude fossils of Melastomataceae for calibrations was avoided or reduced by using alternative fossil and tectonic calibrations, including all those suggested by Morley and Dick. Results reaffirm the relatively recent spread of melastome lineages among the southern continents and refute the breakup of Gondwana as a plausible explanation for the presence of Dissochaeteae/Sonerileae in Madagascar and Africa and the presence of Melastomeae in Africa and Southeast Asia. Melastomeae appear to have reached Africa around 17-15 million years (my) ago, while Dissochaeteae and Sonerileae apparently reached Madagascar at 17-15 and 20-18 my ago. I also explored the effects of constraining Melastomeae to minimally 76 my old (to have reached Africa by island hopping as postulated by Morley and Dick). This resulted in an estimate for their arrival in Africa of 35 my ago and for Dissochaeteae and Sonerileae in Madagascar of 28 and 33 my ago, still implying long-distance dispersal. The Bayesian 95% credibility ranges around these dates, however, are large. Regardless of the increasing sophistication of molecular estimates of divergence time, Gondwanan scenarios will remain untestable as long as biases in the fossil record can justifiably be invoked to explain away the absence of fossils.     

Using n‐dimensional hypervolumes for species distribution modelling: A response to Qiao et al. ()

Hypervolume approaches are used to quantify functional diversity and quantify environmental niches for species distribution modelling. Recently, Qiao et al. ( 2016 ) criticized our geometrical kernel density estimation (KDE) method for measuring hypervolumes. They used a simulation analysis to argue that the method yields high error rates and makes biased estimates of fundamental niches. Here, we show that (a) KDE output depends in useful ways on dataset size and bias, (b) other species distribution modelling methods make equally stringent but different assumptions about dataset bias, (c) simulation results presented by Qiao et al. ( 2016 ) were incorrect, with revised analyses showing performance comparable to other methods, and (d) hypervolume methods are more general than KDE and have other benefits for niche modelling. As a result, our KDE method remains a promising tool for species distribution modelling.     

Comment on "Unbiased statistical analysis for multi-stage proteomic search strategies"

Everett et al. recently reported on a statistical bias that arises in the target-decoy approach to false discovery rate estimation in two-pass proteomics search strategies as exemplified by X!Tandem. This bias can cause serious underestimation of the false discovery rate. We argue here that the "unbiased" solution proposed by Everett et al., however, is also biased and under certain circumstances can also result in a serious underestimate of the FDR, especially at the protein level.     

Implications of uniformly distributed,empirically informed priors for phylogeographical model selection: A reply to Hickerson et al.

Establishing that a set of population‐splitting events occurred at the same time can be a potentially persuasive argument that a common process affected the populations. Recently, Oaks et al. ( 2013 ) assessed the ability of an approximate‐Bayesian model‐choice method (msBayes ) to estimate such a pattern of simultaneous divergence across taxa, to which Hickerson et al. ( 2014 ) responded. Both papers agree that the primary inference enabled by the method is very sensitive to prior assumptions and often erroneously supports shared divergences across taxa when prior uncertainty about divergence times is represented by a uniform distribution. However, the papers differ about the best explanation and solution for this problem. Oaks et al. ( 2013 ) suggested the method's behavior was caused by the strong weight of uniformly distributed priors on divergence times leading to smaller marginal likelihoods (and thus smaller posterior probabilities) of models with more divergence‐time parameters (Hypothesis 1); they proposed alternative prior probability distributions to avoid such strongly weighted posteriors. Hickerson et al. ( 2014 ) suggested numerical‐approximation error causes msBayes analyses to be biased toward models of clustered divergences because the method's rejection algorithm is unable to adequately sample the parameter space of richer models within reasonable computational limits when using broad uniform priors on divergence times (Hypothesis 2). As a potential solution, they proposed a model‐averaging approach that uses narrow, empirically informed uniform priors. Here, we use analyses of simulated and empirical data to demonstrate that the approach of Hickerson et al. ( 2014 ) does not mitigate the method's tendency to erroneously support models of highly clustered divergences, and is dangerous in the sense that the empirically derived uniform priors often exclude from consideration the true values of the divergence‐time parameters. Our results also show that the tendency of msBayes analyses to support models of shared divergences is primarily due to Hypothesis 1, whereas Hypothesis 2 is an untenable explanation for the bias. Overall, this series of papers demonstrates that if our prior assumptions place too much weight in unlikely regions of parameter space such that the exact posterior supports the wrong model of evolutionary history, no amount of computation can rescue our inference. Fortunately, as predicted by fundamental principles of Bayesian model choice, more flexible distributions that accommodate prior uncertainty about parameters without placing excessive weight in vast regions of parameter space with low likelihood increase the method's robustness and power to detect temporal variation in divergences.     

The age of the common ancestor of eukaryotes and prokaryotes: statistical inferences

In this paper, a simple distance measure was used to estimate the age (T) of the common ancestor of eukaryotes and prokaryotes which takes the rate variation among sites and the pattern of amino acid substitutions into account. Our new estimate of T based on Doolittle et al.'s data is about 2.5 billion years ago (Ga), with 95% confidence interval from 2.1 to 2.9 Ga. This result indicates (1) that Doolittle et al.'s estimate (approximately 2.0 Ga) seems too recent, and (2) that the traditional view about the divergence time between eukaryotes and prokaryotes (T0 = 3.5 Ga) can be rejected at the 0.1% significance level.      

A comparison of two popular statistical methods for estimating the time to most recent common ancestor (TMRCA) from a sample of DNA sequences

We have compared two statistical methods of estimating the time to most recent common ancestor (TMRCA) from a sample of DNA sequences, which have been proposed by Templeton (1993) and Bandeltet al. (1995). Monte-Carlo simulations were used for generating DNA sequence data. Different evolutionary scenarios were simulated and the estimation procedures were evaluated. We have found that for both methods (i) the estimates are insensitive to demographic parameters and (ii) the standard deviations of the estimates are too high for these methods to be reliably used in practice.     

Study of high- and low-current-configuration homes from the 1988 Denver Childhood Cancer Study

An epidemiological study conducted by Savitz et al. reported that residential wire codes were more strongly associated with childhood cancer than were measured magnetic fields, a peculiar result because wire codes were originally developed to be a surrogate for residential magnetic fields. The primary purpose of the study reported here, known as the Back to Denver (BTD) study, was to obtain data to help in the interpretation of the original results of Savitz et al. The BTD study included 81 homes that had been occupied by case and control subjects of Savitz et al., stratified by wire code as follows: 18 high current configuration (HCC) case homes; 20 HCC control homes; 20 low current configuration (LCC) case homes; and 23 LCC control homes. Analysis of new data acquired in these homes led to the following previously unpublished conclusions. The home-averaged (i.e., mean of fields measured in subjects' bedrooms, family/living rooms, and rooms where meals normally eaten) spot 60 Hz, 180 Hz, and harmonic (i.e., 60-420 Hz) magnetic fields were associated with wire codes. The 180 Hz and harmonic components, but not the 60 Hz component, were associated with case/control status. Measured static magnetic fields were only weakly correlated (rapproximately 0.2) between rooms in homes. The BTD data provide little support for, but are too sparse to definitively test, the 1995 resonance hypothesis proposed by Bowman et al. Case and control homes had similar concentrations of copper in their tap water. Copper concentration was not associated with wire codes nor with the level of electric current carried by a home's water pipe. These results of the BTD study suggest that future case/control studies investigating power frequency magnetic fields might wish to include measurements of 180 Hz or harmonic magnetic fields in order to examine their associations (if any) with disease status.     

Speciation in the Central American Seaway: the importance of taxon sampling in the identification of trans-isthmian geminate pairs

Aim To create a molecular phylogenetic hypothesis for the closely related serranid genera Alphestes Bloch and Schneider and Dermatolepis Gill and assess the role of the Panamanian Isthmus in speciation within these reef fishes. Location Tropical eastern Pacific, Caribbean, and Indian Oceans. Methods Sequence data from one nuclear (TMO‐4C4) and three mitochondrial genes (16S, 12S, and cytochrome b) were used in maximum parsimony and maximum likelihood analyses. Results Here we show that previously hypothesized trans‐isthmian geminate species are not each other's closest living relatives. Species of Alphestes Bloch and Schneider in the eastern Pacific are sister taxa indicating post‐closure speciation. Within Dermatolepis Gill, we identify a sister group relationship between the Caribbean and western Indian Ocean species, a rarely reported biogeographic pattern. Based on sequence divergence, speciation among the three species of Dermatolepis was, however, nearly simultaneous around the time of the isthmian closure event. Main conclusions Our molecular phylogenetic analysis of two closely related genera of reef fishes, each with presumed trans‐isthmian geminates, cautions against the uncritical use of morphological similarity in identification of geminates, as well as the assumption that trans‐isthmian sister groups date to the isthmian closure event. These findings suggest that in some instances incomplete sampling of species within a clade including putative geminates may lead to improper conclusions regarding the pattern and timing of speciation, as well as incorrect estimation of the rate at which evolution has proceeded.     

Relaxed clocks and inferences of heterogeneous patterns of nucleotide substitution and divergence time estimates across whales and dolphins (Mammalia: Cetacea)

Various nucleotide substitution models have been developed to accommodate among lineage rate heterogeneity, thereby relaxing the assumptions of the strict molecular clock. Recently developed "uncorrelated relaxed clock" and "random local clock" (RLC) models allow decoupling of nucleotide substitution rates between descendant lineages and are thus predicted to perform better in the presence of lineage-specific rate heterogeneity. However, it is uncertain how these models perform in the presence of punctuated shifts in substitution rate, especially between closely related clades. Using cetaceans (whales and dolphins) as a case study, we test the performance of these two substitution models in estimating both molecular rates and divergence times in the presence of substantial lineage-specific rate heterogeneity. Our RLC analyses of whole mitochondrial genome alignments find evidence for up to ten clade-specific nucleotide substitution rate shifts in cetaceans. We provide evidence that in the uncorrelated relaxed clock framework, a punctuated shift in the rate of molecular evolution within a subclade results in posterior rate estimates that are either misled or intermediate between the disparate rate classes present in baleen and toothed whales. Using simulations, we demonstrate abrupt changes in rate isolated to one or a few lineages in the phylogeny can mislead rate and age estimation, even when the node of interest is calibrated. We further demonstrate how increasing prior age uncertainty can bias rate and age estimates, even while the 95% highest posterior density around age estimates decreases; in other words, increased precision for an inaccurate estimate. We interpret the use of external calibrations in divergence time studies in light of these results, suggesting that rate shifts at deep time scales may mislead inferences of absolute molecular rates and ages.     

Do NMDA receptor kinetics regulate the end of critical periods of plasticity?

Increasing NR2A subunit expression and the associated shortening of the NMDA-EPSC are thought to underlie the loss of diverse types of sensory cortical plasticity. Lu and colleagues (Lu et al., 2001 [this issue of Neuron]) now report that mice lacking the NR2A subunit display normal duration of critical periods of barrel cortex plasticity. Shortening of the NMDA-EPSC is therefore not responsible for the end of these critical periods.     

On the Use of Gene Ontology Annotations to Assess Functional Similarity among Orthologs and Paralogs: A Short Report

A recent paper (Nehrt et al., PLoS Comput. Biol. 7:e1002073, 2011) has proposed a metric for the "functional similarity" between two genes that uses only the Gene Ontology (GO) annotations directly derived from published experimental results. Applying this metric, the authors concluded that paralogous genes within the mouse genome or the human genome are more functionally similar on average than orthologous genes between these genomes, an unexpected result with broad implications if true. We suggest, based on both theoretical and empirical considerations, that this proposed metric should not be interpreted as a functional similarity, and therefore cannot be used to support any conclusions about the "ortholog conjecture" (or, more properly, the "ortholog functional conservation hypothesis"). First, we reexamine the case studies presented by Nehrt et al. as examples of orthologs with divergent functions, and come to a very different conclusion: they actually exemplify how GO annotations for orthologous genes provide complementary information about conserved biological functions. We then show that there is a global ascertainment bias in the experiment-based GO annotations for human and mouse genes: particular types of experiments tend to be performed in different model organisms. We conclude that the reported statistical differences in annotations between pairs of orthologous genes do not reflect differences in biological function, but rather complementarity in experimental approaches. Our results underscore two general considerations for researchers proposing novel types of analysis based on the GO: 1) that GO annotations are often incomplete, potentially in a biased manner, and subject to an "open world assumption" (absence of an annotation does not imply absence of a function), and 2) that conclusions drawn from a novel, large-scale GO analysis should whenever possible be supported by careful, in-depth examination of examples, to help ensure the conclusions have a justifiable biological basis.     

Accuracy of rate estimation using relaxed-clock models with a critical focus on the early metazoan radiation

In recent years, a number of phylogenetic methods have been developed for estimating molecular rates and divergence dates under models that relax the molecular clock constraint by allowing rate change throughout the tree. These methods are being used with increasing frequency, but there have been few studies into their accuracy. We tested the accuracy of several relaxed-clock methods (penalized likelihood and Bayesian inference using various models of rate change) using nucleotide sequences simulated on a nine-taxon tree. When the sequences evolved with a constant rate, the methods were able to infer rates accurately, but estimates were more precise when a molecular clock was assumed. When the sequences evolved under a model of auto-correlated rate change, rates were accurately estimated using penalized likelihood and by Bayesian inference using lognormal and exponential models of rate change, while other models did not perform as well. When the sequences evolved under a model of uncorrelated rate change, only Bayesian inference using an exponential rate model performed well. Collectively, the results provide a strong recommendation for using the exponential model of rate change if a conservative approach to divergence time estimation is required. A case study is presented in which we use a simulation-based approach to examine the hypothesis of elevated rates in the Cambrian period, and it is found that these high rate estimates might be an artifact of the rate estimation method. If this bias is present, then the ages of metazoan divergences would be systematically underestimated. The results of this study have implications for studies of molecular rates and divergence dates.     

Screening,Identification and Distribution of Endophytic Associative Diazotrophs Isolated from Rice Plants

29 strains of endophyte associated diazotrophs isolated from rice ( Oryza sutiva L. ) "Yuefu" plant were selected by in vitro acetylene reduced activity and 15 N2-flxing activity determination. They were identificated into 14 species of 9 genera: Agrobacterium tumefaciens (Smith et Townsend) Conn, A. radiobacter (Beijerinck et van Delden) Conn; Alcaligenes piechaudii Kiredjian et al., Al. denitrificans (Leifson et Hagh) Ruger et Tan; Bacillus sphaericus Meyer et Neide, B. licheniforrnis Weigmann Chester, B. cereus Frankland et Frankland; Chryseomonas luteola (Kodama, Kinnura et Komagata) Holmes et al.; Enterobacter cloacae ( Jordan ) Honnaeche et Edwards, E. sakazakii Famer et al.; E. agglomeraus (Beijerinck) Ewing et Fife; Pseudomonas pseudoalcaligenes Stanier, P. alcaligeaes Monias, P. putida Biotype A Evans and Aeromonas Kluyver et van Niel, Serratia Bizio, Staphyloccocus Rosenbach, Xanthomonas Dowson. Among them C. luteola, E. sakazakii, E. agglomeraus, P. pseudoalcaligenes have not been previously reported as diazotrophs. Studying on the distribution of endophytic dizoatrophs in rice seeds and rice plants has demonstrated that the diversity of endophytic associated diazotrophs in roots was more than that in other organs of rice plant. Endophytic associative diazotroph E. cloacae not only occurred on the surface but also occurred intercellularly in R48 rice roots as observed with scanning electron microscopy and transmission microscopy.     

Integrating fossil preservation biases in the selection of calibrations for molecular divergence time estimation

The selection of fossil data to use as calibration age priors in molecular divergence time estimates inherently links neontological methods with paleontological theory. However, few neontological studies have taken into account the possibility of a taphonomic bias in the fossil record when developing approaches to fossil calibration selection. The Sppil-Rongis effect may bias the first appearance of a lineage toward the recent causing most objective calibration selection approaches to erroneously exclude appropriate calibrations or to incorporate multiple calibrations that are too young to accurately represent the divergence times of target lineages. Using turtles as a case study, we develop a Bayesian extension to the fossil selection approach developed by Marshall (2008. A simple method for bracketing absolute divergence times on molecular phylogenies using multiple fossil calibrations points. Am. Nat. 171:726-742) that takes into account this taphonomic bias. Our method has the advantage of identifying calibrations that may bias age estimates to be too recent while incorporating uncertainty in phylogenetic parameter estimates such as tree topology and branch lengths. Additionally, this method is easily adapted to assess the consistency of potential calibrations to any one calibration in the candidate pool.     

Taxonomic status of Chinese bahaba (Bahaba taipingensis) and its phylogenetic relationship with other species in the family Sciaenidae

Sciaenidae is one of the largest fish families, but the phylogeny and taxonomy of these fishes are still being disputed. Furthermore, the taxonomic status of the Chinese bahaba (Bahaba taipingensis), which is an endemic species to China, had never been studied through molecular method. In this study, phylogenetic relationships among sciaenid species were reconstructed using DNA sequence data from the mitochondrial cytochrome b (cyt b) gene and the 16S ribosomal RNA (16S rRNA) gene through Bayesian inference (BI) analyses. The phylogenetic trees indicated that the Chinese bahaba closely related to Collichthys and Pseudosciaena. Previous studies (Meng et al. 2004. Prog Nat Sci 14(5):514-521, [in Chinese]) showed that the subfamily Pseudosciaeninae represented the latest evolutionary sort, which was more suitable with the current environment. Based on our 16S rRNA, the Chinese bahaba showed close relationship with Pseudosciaena, thus the divergence of the Chinese bahaba maybe also very late within the family Sciaenidae. It was probably that the oceanographic or ecological discontinuities of these species vary considerably causing particularly strong breaks. Furthermore, combined with previous studies, we suggested that there was only one genus, Otolithes, within the subfamily Otolithinae. Nevertheless, the taxonomic status of the Chinese bahaba and the phylogeny of the whole family were not completely solved. Additional samples of more species are required to develop a clearer picture of the evolutionary history of Sciaenidae.