Inference of Microbial Recombination Rates from Metagenomic Data

Metagenomic sequencing projects from environments dominated by a small number of species produce genome-wide population samples. We present a two-site composite likelihood estimator of the scaled recombination rate, ρ = 2N_ec, that operates on metagenomic assemblies in which each sequenced fragment derives from a different individual. This new estimator properly accounts for sequencing error, as quantified by per-base quality scores, and missing data, as inferred from the placement of reads in a metagenomic assembly. We apply our estimator to data from a sludge metagenome project to demonstrate how this method will elucidate the rates of exchange of genetic material in natural microbial populations. Surprisingly, for a fixed amount of sequencing, this estimator has lower variance than similar methods that operate on more traditional population genetic samples of comparable size. In addition, we can infer variation in recombination rate across the genome because metagenomic projects sample genetic diversity genome-wide, not just at particular loci. The method itself makes no assumption specific to microbial populations, opening the door for application to any mixed population sample where the number of individuals sampled is much greater than the number of fragments sequenced.     

Phylogenetic Inference and the Misplaced Premise of Substitution Rates

Acta Biotheoretica - Three competing ‘methods’ have been endorsed for inferring phylogenetic hypotheses: parsimony, likelihood, and Bayesianism. The latter two have been claimed...     

PALM: A Paralleled and Integrated Framework for Phylogenetic Inference with Automatic Likelihood Model Selectors

Background

Selecting an appropriate substitution model and deriving a tree topology for a given sequence set are essential in phylogenetic analysis. However, such time consuming, computationally intensive tasks rely on knowledge of substitution model theories and related expertise to run through all possible combinations of several separate programs. To ensure a thorough and efficient analysis and avert tedious manipulations of various programs, this work presents an intuitive framework, the phylogenetic reconstruction with automatic likelihood model selectors (PALM), with convincing, updated algorithms and a best-fit model selection mechanism for seamless phylogenetic analysis.

Methodology

As an integrated framework of ClustalW, PhyML, MODELTEST, ProtTest, and several in-house programs, PALM evaluates the fitness of 56 substitution models for nucleotide sequences and 112 substitution models for protein sequences with scores in various criteria. The input for PALM can be either sequences in FASTA format or a sequence alignment file in PHYLIP format. To accelerate the computing of maximum likelihood and bootstrapping, this work integrates MPICH2/PhyML, PalmMonitor and Palm job controller across several machines with multiple processors and adopts the task parallelism approach. Moreover, an intuitive and interactive web component, PalmTree, is developed for displaying and operating the output tree with options of tree rooting, branches swapping, viewing the branch length values, and viewing bootstrapping score, as well as removing nodes to restart analysis iteratively.

Significance

The workflow of PALM is straightforward and coherent. Via a succinct, user-friendly interface, researchers unfamiliar with phylogenetic analysis can easily use this server to submit sequences, retrieve the output, and re-submit a job based on a previous result if some sequences are to be deleted or added for phylogenetic reconstruction. PALM results in an inference of phylogenetic relationship not only by vanquishing the computation difficulty of ML methods but also providing statistic methods for model selection and bootstrapping. The proposed approach can reduce calculation time, which is particularly relevant when querying a large data set. PALM can be accessed online at http://palm.iis.sinica.edu.tw.     

删失数据下非线性半参数回归模型中参数的经验似然推断

考察了响应变量在随机删失情形下的非线性半参数回归模型,构造了未知参数的经验对数似然比统计量和调整经验对数似然比统计量,证明在一定条件下,所构造的经验似然比统计量渐近于X~2分布,并由此构造出未知参数的置信域.此外,又构造了未知参数的最小二乘估计量,证明了它的渐近性质.通过模拟研究表明,经验似然方法在置信域的覆盖概率以及精度方面要优于最小二乘法.     

Assessing Variability by Joint Sampling of Alignments and Mutation Rates

When two sequences are aligned with a single set of alignment parameters, or when mutation parameters are estimated on the basis of a single ``optimal' sequence alignment, the variability of both the alignment and the estimated parameters can be seriously underestimated. To obtain a more realistic impression of the actual uncertainty, we propose sampling sequence alignments and mutation parameters simultaneously from their joint posterior distribution given the two original sequences. We illustrate our method with human and orangutan sequences from the hyper variable region I and with gene–pseudogene pairs. Received: 16 November 2000 / Accepted: 15 May 2001     

Accuracy and Power of the Likelihood Ratio Test for Comparing Evolutionary Rates Among Genes

Sequences for multiple protein-coding genes are now commonly available from several, often closely related species. These data sets offer intriguing opportunities to test hypotheses regarding whether different types of genes evolve under different selective pressures. Although maximum likelihood (ML) models of codon substitution that are suitable for such analyses have been developed, little is known about the statistical properties of these tests. We use a previously developed fixed-sites model and computer simulations to examine the accuracy and power of the likelihood ratio test (LRT) in comparing the nonsynonymous-to-synonymous substitution rate ratio (=dN/dS) between two genes. Our results show that the LRT applied to fixed-sites models may be inaccurate in some cases when setting significance thresholds using a ² approximation. Instead, we use a parametric bootstrap to describe the distribution of the LRT statistic for fixed-sites models and examine the power of the test as a function of sampling variables and properties of the genes under study. We find that the power of the test is high (>80%) even when sampling few taxa (e.g., six species) if sequences are sufficiently diverged and the test is largely unaffected by the tree topology used to simulate data. Our simulations show fixed-sites models are suitable for comparing substitution parameters among genes evolving under even strong evolutionary constraint ( 0.05), although relative rate differences of 25% or less may be difficult to detect.Reviewing Editor: Dr. Rosmus Nielsen     

Relationships of Cetacea (Artiodactyla) Among Mammals: Increased Taxon Sampling Alters Interpretations of Key Fossils and Character Evolution

Background

Integration of diverse data (molecules, fossils) provides the most robust test of the phylogeny of cetaceans. Positioning key fossils is critical for reconstructing the character change from life on land to life in the water.

Methodology/Principal Findings

We reexamine relationships of critical extinct taxa that impact our understanding of the origin of Cetacea. We do this in the context of the largest total evidence analysis of morphological and molecular information for Artiodactyla (661 phenotypic characters and 46,587 molecular characters, coded for 33 extant and 48 extinct taxa). We score morphological data for Carnivoramorpha, †Creodonta, Lipotyphla, and the †raoellid artiodactylan †Indohyus and concentrate on determining which fossils are positioned along stem lineages to major artiodactylan crown clades. Shortest trees place Cetacea within Artiodactyla and close to †Indohyus, with †Mesonychia outside of Artiodactyla. The relationships of †Mesonychia and †Indohyus are highly unstable, however - in trees only two steps longer than minimum length, †Mesonychia falls inside Artiodactyla and displaces †Indohyus from a position close to Cetacea. Trees based only on data that fossilize continue to show the classic arrangement of relationships within Artiodactyla with Cetacea grouping outside the clade, a signal incongruent with the molecular data that dominate the total evidence result.

Conclusions/Significance

Integration of new fossil material of †Indohyus impacts placement of another extinct clade †Mesonychia, pushing it much farther down the tree. The phylogenetic position of †Indohyus suggests that the cetacean stem lineage included herbivorous and carnivorous aquatic species. We also conclude that extinct members of Cetancodonta (whales + hippopotamids) shared a derived ability to hear underwater sounds, even though several cetancodontans lack a pachyostotic auditory bulla. We revise the taxonomy of living and extinct artiodactylans and propose explicit node and stem-based definitions for the ingroup.     

Ecological and Ecomorphological Specialization Are Not Associated with Diversification Rates in Muroid Rodents (Rodentia: Muroidea)

Multiple diversification rate shifts explain uneven clade richness in muroid rodents. Previous muroid studies have shown that extrinsic factors, notwithstanding ecological opportunity, are poor predictors of clade diversity. Here, we use a 297-muroid species chronogram that is sampled proportional to total clade diversity, along with various trait-dependent diversification approaches to investigate the association between diversification rates with intrinsic attributes—diet, habitat, body mass, and relative tail length. We found some association between both dietary specialization and body mass, as well as between habitat specialization with relative tail lengths using phylogenetic analyses of variance. However, there was no significant association between diversification rates with the evolution of these traits in muroid rodents. We also show that several of the state-dependent diversification approaches are highly susceptible to Type I error—a result that is in accordance with recent criticisms of these methods. Finally, we discuss several potential causes for the lack of association between the examined trait data with diversification rates, ranging from methodological biases (e.g. method conservativism) to biology (e.g. behavioral plasticity and ecological opportunism of muroid rodents).     

A General HIV Incidence Inference Scheme Based on Likelihood of Individual Level Data and a Population Renewal Equation

We derive a new method to estimate the age specific incidence of an infection with a differential mortality, using individual level infection status data from successive surveys. The method consists of a) an SI-type model to express the incidence rate in terms of the prevalence and its derivatives as well as the difference in mortality rate, and b) a maximum likelihood approach to estimate the prevalence and its derivatives. Estimates can in principle be obtained for any chosen age and time, and no particular assumptions are made about the epidemiological or demographic context. This is in contrast with earlier methods for estimating incidence from prevalence data, which work with aggregated data, and the aggregated effect of demographic and epidemiological rates over the time interval between prevalence surveys. Numerical simulation of HIV epidemics, under the presumption of known excess mortality due to infection, shows improved control of bias and variance, compared to previous methods. Our analysis motivates for a) effort to be applied to obtain accurate estimates of excess mortality rates as a function of age and time among HIV infected individuals and b) use of individual level rather than aggregated data in order to estimate HIV incidence rates at times between two prevalence surveys.     

Inference of Type-Specific HPV Transmissibility,Progression and Clearance Rates: A Mathematical Modelling Approach

Quantifying rates governing the clearance of Human Papillomavirus (HPV) and its progression to clinical disease, together with viral transmissibility and the duration of naturally-acquired immunity, is essential in estimating the impact of vaccination programmes and screening or testing regimes. However, the complex natural history of HPV makes this difficult. We infer the viral transmissibility, rate of waning natural immunity and rates of progression and clearance of infection of 13 high-risk and 2 non-oncogenic HPV types, making use of a number of rich datasets from Sweden. Estimates of viral transmissibility, clearance of initial infection and waning immunity were derived in a Bayesian framework by fitting a susceptible-infectious-recovered-susceptible (SIRS) transmission model to age- and type-specific HPV prevalence data from both a cross-sectional study and a randomised controlled trial (RCT) of primary HPV screening. The models fitted well, but over-estimated the prevalence of four high-risk types with respect to the data. Three of these types (HPV-33, -35 and -58) are among the most closely related phylogenetically to the most prevalent HPV-16. The fourth (HPV-45) is the most closely related to HPV-18; the second most prevalent type. We suggest that this may be an indicator of cross-immunity. Rates of progression and clearance of clinical lesions were additionally estimated from longitudinal data gathered as part of the same RCT. Our estimates of progression and clearance rates are consistent with the findings of survival analysis studies and we extend the literature by estimating progression and clearance rates for non-16 and non-18 high-risk types. We anticipate that such type-specific estimates will be useful in the parameterisation of further models and in developing our understanding of HPV natural history.     

Delineating Species with DNA Barcodes: A Case of Taxon Dependent Method Performance in Moths

The accelerating loss of biodiversity has created a need for more effective ways to discover species. Novel algorithmic approaches for analyzing sequence data combined with rapidly expanding DNA barcode libraries provide a potential solution. While several analytical methods are available for the delineation of operational taxonomic units (OTUs), few studies have compared their performance. This study compares the performance of one morphology-based and four DNA-based (BIN, parsimony networks, ABGD, GMYC) methods on two groups of gelechioid moths. It examines 92 species of Finnish Gelechiinae and 103 species of Australian Elachistinae which were delineated by traditional taxonomy. The results reveal a striking difference in performance between the two taxa with all four DNA-based methods. OTU counts in the Elachistinae showed a wider range and a relatively low (ca. 65%) OTU match with reference species while OTU counts were more congruent and performance was higher (ca. 90%) in the Gelechiinae. Performance rose when only monophyletic species were compared, but the taxon-dependence remained. None of the DNA-based methods produced a correct match with non-monophyletic species, but singletons were handled well. A simulated test of morphospecies-grouping performed very poorly in revealing taxon diversity in these small, dull-colored moths. Despite the strong performance of analyses based on DNA barcodes, species delineated using single-locus mtDNA data are best viewed as OTUs that require validation by subsequent integrative taxonomic work.     

Rates of Species Accumulation and Taxonomic Diversification during Phototrophic Biofilm Development Are Controlled by both Nutrient Supply and Current Velocity

The accumulation of new and taxonomically diverse species is a marked feature of community development, but the role of the environment in this process is not well understood. To address this problem, we subjected periphyton in laboratory streams to low (10-cm · s⁻¹), high (30-cm · s⁻¹), and variable (9- to 32-cm · s⁻¹) current velocity and low- versus high-nutrient inputs. We examined how current velocity and resource supply constrained (i) the rates of species accumulation, a measure of temporal beta-diversity, and (ii) the rates of diversification of higher taxonomic categories, defined here as the rate of higher taxon richness increase with the increase of species richness. Temporal biofilm dynamics were controlled by a strong nutrient-current interaction. Nutrients accelerated the rates of accumulation of new species, when flow velocity was not too stressful. Species were more taxonomically diverse under variable than under low-flow conditions, indicating that flow heterogeneity increased the niche diversity in the high-nutrient treatments. Conversely, the lower diversification rates under high- than under low-nutrient conditions at low velocity are explained with finer resource partitioning among species, belonging to a limited number of related genera. The overall low rates of diversification in high-current treatments suggest that the ability to withstand current stress was conserved within closely related species. Temporal heterogeneity of disturbance has been shown to promote species richness, but here we further demonstrate that it also affects two other components of biodiversity, i.e., temporal beta-diversity and diversification rate. Therefore, management efforts for preserving the inherent temporal heterogeneity of natural ecosystems will have detectable positive effects on biodiversity.     

Quantifying the Likelihood of Co-existence for Communities with Asymmetric Competition

Trade-offs in performance of different ecological functions within a species are commonly offered as an explanation for co-existence in natural communities. Single trade-offs between competitive ability and other life history traits have been shown to support a large number of species, as a result of strong competitive asymmetry. We consider a single competition-fecundity trade-off in a homogeneous environment, and examine the effect of the form of asymmetry on the likelihood of species co-existing. We find conditions that allow co-existence of two species for a general competition function, and show that (1)?two species can only co-exist if the competition function is sufficiently steep when the species are similar; (2)?when competition is determined by a linear function, no more than two species can co-exist; (3)?when the competition between two individuals is determined by a discontinuous step function, this single trade-off can support an arbitrarily large number of species. Further, we show analytically that as the degree of asymmetry in competition increases, the probability of a given number of species co-existing also increases, but note that even in the most favourable conditions, large numbers of species co-existing along a single trade-off is highly unlikely. On this basis, we suggest it is unlikely that single trade-offs are able to support high levels of bio-diversity without interacting other processes.     

Analysis of Partially Balanced Incomplete Block Designs with Standards

The paper presents a detailed analysis of variance of PBIB designs supplemented by one or more standard treatments. Comparisons were carried out of treatment effects within the group of standard and test treatments and between the two groups. The analysis is supplemented by an example.     

Incomplete division of labor: error-prone multitaskers coexist with specialists

Division of labor (DoL) occurs when individual members of a group specialize by performing particular tasks toward some common goal. Under complete DoL, every individual acts as a specialist and so performs only one kind of task. But under incomplete DoL, some individuals may act as generalists and so have the capacity to perform more than one kind of task. This persistence of generalists in the presence of specialists presents a theoretical challenge, particularly if generalists must pay an extra cost, an inefficiency penalty, for their capacity to perform more than one type of task. Prior work focused on how such costs tend to drive evolution toward complete DoL, with only specialists typically remaining at equilibrium [Wahl, L.M., 2002a. Evolving the division of labor: generalists, specialists and task allocation. J. Theoret. Biol. 219, 371-388; Wahl, L.M., 2002b. The division of labor: genotypic versus phenotypic specialization. Am. Nat. 160, 135-145]. Relaxing this key assumption, we show that generalists, despite paying some extra costs, can coexist with specialists. Relaxing another assumption, we also show that this coexistence can hold even when generalists often perform the wrong task. How can stable multitasking emerge despite this flawed decision-making? From the perspective that cognitive errors matter only when they translate into fitness decrements, we observe that error-prone generalists may persist most commonly in situations in which their mistakes do little to jeopardize group success. Our findings show that incomplete DoL can emerge even when generalists often err and must pay extra costs for their multitasking capacity.     

湖南参属的归并及亚洲楤木属(五加科)的一个新种和一些名称的变更

分子与形态证据表明湖南参 (HunaniopanaxhypoglaucusC .J.Qi&T .R .Cao)起源于木属 (AraliaL .) ,湖南参的近缘种为寄生五叶参 (Araliaparasitica)和轮伞五叶参 (Araliaverticillata) ,形态证据支持湖南参为轮伞五叶参的姐妹种 ,此 2种具轮伞花序这一共衍征。为了保证木属的单系性 ,现将湖南参并入木属 ,并作如下新组合 :Araliahypoglauca (C .J .Qi&T .R .Cao)J.Wen&Y .F .Deng .本文报道了湖南参在广西的新分布。目前的形态学和分子数据均不支持早先提出的湖南参与树参属 (Dendropanax)的近缘关系。本文亦报道了木属的 1新种、 2新种组合和 1个新名称 :AraliashangianaJ .Wen .sp .nov .(向氏五叶参 ,新种 ) ;Araliaglabrifoliolata (C .B .Shang)J .Wen .comb .nov . (光叶五叶参 ,新组合 ) ;Araliastellata (King)J.Wen ,comb .nov . (星毛羽叶参 ,新组合 )和AraliadelavayiJ .Wen .nom .nov .(云南五叶参 ,新名称 )。     

On Testing Equality of Means of Correlated Variates with Incomplete Data

A statistic is proposed for testing the hypothesis of equality of the means of a bivariate normal distribution with unknown common variance and correlation coefficient when observations are missing on one of the variates. Expressions for the second and fourth central moments of the statistic are obtained. These moments are used to approximate the distribution of the statistic by a Student's t distribution under the null hypothesis. The powers of the test are computed and compared with those of the conventional paired t and the other known statistics.