Examining Phylogenetic Relationships Among Gibbon Genera Using Whole Genome Sequence Data Using an Approximate Bayesian Computation Approach

Gibbons are believed to have diverged from the larger great apes ∼16.8 MYA and today reside in the rainforests of Southeast Asia. Based on their diploid chromosome number, the family Hylobatidae is divided into four genera, Nomascus, Symphalangus, Hoolock, and Hylobates. Genetic studies attempting to elucidate the phylogenetic relationships among gibbons using karyotypes, mitochondrial DNA (mtDNA), the Y chromosome, and short autosomal sequences have been inconclusive . To examine the relationships among gibbon genera in more depth, we performed second-generation whole genome sequencing (WGS) to a mean of ∼15× coverage in two individuals from each genus. We developed a coalescent-based approximate Bayesian computation (ABC) method incorporating a model of sequencing error generated by high coverage exome validation to infer the branching order, divergence times, and effective population sizes of gibbon taxa. Although Hoolock and Symphalangus are likely sister taxa, we could not confidently resolve a single bifurcating tree despite the large amount of data analyzed. Instead, our results support the hypothesis that all four gibbon genera diverged at approximately the same time. Assuming an autosomal mutation rate of 1 × 10⁻⁹/site/year this speciation process occurred ∼5 MYA during a period in the Early Pliocene characterized by climatic shifts and fragmentation of the Sunda shelf forests. Whole genome sequencing of additional individuals will be vital for inferring the extent of gene flow among species after the separation of the gibbon genera.     

Phylogenetic Relationships of Amaryllidaceae Based on matK Sequence Data

matK gene, which is located in the chloroplast genome and evolves more quickly than the rbcL gene. A total of 31 species representing 31 of the 59 genera in the family were examined in this study. We also used 21 species from another ten families of Asparagales, four species from three families of Liliales and Acorus as outgroups. We obtained partial sequences of matK with lengths of 1,109–1,148 bp, corresponding to positions 230 to 1,343 of the Oryza sativa matK gene. The pairwise percentage sequence divergence ranged from 0 to 19.1% for all the species examined except Acorus, and 0 to 4.6% within Amaryllidaceae. Two methods of phylogenetic analysis, the Maximum Parsimony and Neighbor-Joining methods, were used. The trees obtained from these two analyses were fundamentally consistent. In both trees, the Amaryllidaceae sensu Dahlgren et al. formed a well-supported monophyletic clade with 100% bootstrap support. Amaryllidaceae were included in the Asparagales; however, its phylogenetic position within the Asparagales was not clearly resolved. Judging from the NJ tree, Agapanthus might be a sister group of the Amaryllidaceae, although bootstrap support for this was low. Character-state mapping was used to infer a center of origin and the biogeographic history of Amaryllidaceae. The result supports the hypothesis that the family evolved in Africa and subsequently spread to other continents, further suggesting that South America is the center of secondary diversification. Received 6 January 1999/ Accepted in revised form 8 April 1999     

A Comparison Between Several Goodness of Fit Tests for Cox's Model

A comparison is made between two approaches to testing goodness of fit of Cox's regression model for survival data. The first approach is based on the inclusion of time dependent covariates, whereas the second one is based on the autocovariance of successive contributions to the derivative of the loglikelihood. It appears that the second test is most appropriate for testing in situations where the structure of the departure from proportional hazards is not known a priori. An approximate expression for the relative efficiency of the two test procedures is presented.     

A Power Law Model for Survival Data Using GLIM

This paper proposes a regression model for the Weibull survival distribution of which the scale parameter is a power function of covariates. The estimation of parameters for partially censored data is pursued by using a statistical package called GLIM. Two sets of carcinogenic data are used to illustrate this procedure.     

Analysis of Dichotomized Factorial Data Using Cox's Proportional Hazards Model and Related Tests

This paper describes how Cox's Proportional Hazards model may be used to analyze dichotomized factorial data obtained from a right-censored epidemiological study where time to response is of interest. Exact maximum likelihood estimates of the relative mortality rates are derived for any number of prognostic factors, but for the sake of simplicity, the mathematical details are presented for the case of two factors. This method is not based on the life table procedure. Kaplan-Meier estimates are obtained for the survival function of the internal control population, Which are in turn used to determine the expected number of deaths in the study population. The asymptotic (large sample) joint sampling distribution of the relative mortality rates is derived and some relevant simultaneous and conditional statistical tests are discussed. The relative mortality rates of several prognostic factors may be jointly considered as the multivariate extension of the familiar standard mortality ratio (SMR) of epidemiological studies. A numerical example is discussed to illustrate the method.     

Relative Efficiencies of a Class of Goodness of Fit Tests for Randomly Censored Data

We investigate the relative performances of a class of goodness of fit procedures for randomly censored data. For purposes of planning experiments, we quantify the loss of information induced by censorship. We evaluate efficiencies against particular alternatives of interest in survival studies, as the amount of censorship increases. We caution against attributing various power and efficiency properties of the goodness of fit criteria that are obtained under no censorship to situations where the censorship is far from negligible.     

基于matR基因序列分析的山茶科系统关系

通过线粒体matR基因序列分析探讨了山茶科的分类学范围和系统演化关系。结果显示,传统山茶科的两个核心——山茶亚科（Theoideae或Camellioideae）和厚皮香亚科（Ternstroemioideae）不构成姐妹群关系,山茶亚科是一个支持率很高的单系类群,厚皮香亚科没有形成单系;山茶亚科下可区分出3个明显的分支,基部的分支由紫茎属（Stewartia）和舟柄茶属（Hartia）组成,木荷属（Schima）、美洲荷属（Franklirda）和美国大头茶属（Gordonia）构成第2个分支,该分支与由山茶属（Camellia）、核果茶属（Pyrenaria）、多瓣核果茶属（Parapyrenaria）、石笔木属（Tutcheria）、大头荣属（Polyspora）和圆籽荷属（Aptersperma）组成的第3个分支互为姐妹群。研究结果很好地支持了Prince和Parks等学者提出的的狭义山茶科（仅含山茶亚科）和狭义大头茶属的概念以及科下3个族（紫茎族Stewartieae、大头茶族Gordonieae和山茶族Theeae）的划分。但本研究更为清晰地揭示了科下3个族间的系统关系,即紫茎族是最基部的分支,山茶族与大头茶族间有更近的亲缘关系。同时,本文认为,厚皮香（亚）科是否为单系类群值得进一步研究。     

Using Genome-Wide Protein Sequence Data to Predict Amino Acid Conservation

For most proteins, multiple sequence alignments are a viable method to identify functionally and structurally important amino acids, but for most organisms, there is a subset of proteins that are unique or found in a few closely related organisms. For these proteins, it is not possible to produce sequence alignments that are useful in identifying functionally or structurally important amino acids. We have investigated the relationship between amino acid conservation and five factors (the amino acid’s identity, N-terminal neighbor, C-terminal neighbor, the local hydropathy of surrounding amino acids, and the local expected net charge of the surrounding amino acids based on the primary sequence) in Escherichia coli proteins. For four of the factors examined (all but the amino acid’s identity), there is a significant relationship with conservation for some of the standard 20 amino acids. Using the combination of all five factors, we show that it is possible to calculate a score based on the primary sequences of a subset of E. coli proteins that has statistically significant predictive value with respect to predicting conserved amino acids in other E. coli proteins and Saccharomyces cerevisiae proteins. As these five variables show significant relationships with conservation, we have termed them conservation factors. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Phylogenetic Diversity of New Zealand Gelidiales as Revealed by rbcL Sequence Data

Diversity and phylogenetic relationships of New Zealand representatives of the red algal order Gelidiales have been examined using rbcL sequence data. Extensive field collections have been made from throughout the New Zealand region. Six genera have been reported previously from New Zealand (Capreolia, Gelidium, Pterocladia, Pterocladiella, Pterocladiastrum, Ptilophora). This research has revealed species with very restricted local distributions, as well as the discovery of several undescribed, cryptic taxa. The common and widespread Gelidium caulacantheum is confirmed to be more closely related to Capreolia than to other species of Gelidium. The generic concept of Capreolia, based on life history characters, will need to be modified to accommodate additional species possessing “Gelidium” life histories. A species endemic to New Zealand, Gelidium ceramoides, has been found to differ significantly from all other members of the Gelidiales and requires reclassification in another genus and order. Examination of field collections and herbarium specimens in addition to molecular sequence data have led us to conclude that specimens previously placed in the genera Ptilophora and Pterocladiastrum belong within Pterocladia lucida.     

The Power of Relative Rates Tests Depends on the Data

One of the most useful features of molecular phylogenetic analyses is the potential for estimating dates of divergence of evolutionary lineages from the DNA of extant species. But lineage-specific variation in rate of molecular evolution complicates molecular dating, because a calibration rate estimated from one lineage may not be an accurate representation of the rate in other lineages. Many molecular dating studies use a ``clock test' to identify and exclude sequences that vary in rate between lineages. However, these clock tests should not be relied upon without a critical examination of their effectiveness at removing rate variable sequences from any given data set, particularly with regard to the sequence length and number of variable sites. As an illustration of this problem we present a power test of a frequently employed triplet relative rates test. We conclude that (1) relative rates tests are unlikely to detect moderate levels of lineage-specific rate variation (where one lineage has a rate of molecular evolution 1.5 to 4.0 times the other) for most commonly used sequences in molecular dating analyses, and (2) this lack of power is likely to result in substantial error in the estimation of dates of divergence. As an example, we show that the well-studied rate difference between murid rodents and great apes will not be detected for many of the sequences used to date the divergence between these two lineages and that this failure to detect rate variation is likely to result in consistent overestimation the date of the rodent–primate split. Received: 9 June 1999 / Accepted: 22 October 1999     

Some Limitations of Public Sequence Data for Phylogenetic Inference (in Plants)

The GenBank database contains essentially all of the nucleotide sequence data generated for published molecular systematic studies, but for the majority of taxa these data remain sparse. GenBank has value for phylogenetic methods that leverage data–mining and rapidly improving computational methods, but the limits imposed by the sparse structure of the data are not well understood. Here we present a tree representing 13,093 land plant genera—an estimated 80% of extant plant diversity—to illustrate the potential of public sequence data for broad phylogenetic inference in plants, and we explore the limits to inference imposed by the structure of these data using theoretical foundations from phylogenetic data decisiveness. We find that despite very high levels of missing data (over 96%), the present data retain the potential to inform over 86.3% of all possible phylogenetic relationships. Most of these relationships, however, are informed by small amounts of data—approximately half are informed by fewer than four loci, and more than 99% are informed by fewer than fifteen. We also apply an information theoretic measure of branch support to assess the strength of phylogenetic signal in the data, revealing many poorly supported branches concentrated near the tips of the tree, where data are sparse and the limiting effects of this sparseness are stronger. We argue that limits to phylogenetic inference and signal imposed by low data coverage may pose significant challenges for comprehensive phylogenetic inference at the species level. Computational requirements provide additional limits for large reconstructions, but these may be overcome by methodological advances, whereas insufficient data coverage can only be remedied by additional sampling effort. We conclude that public databases have exceptional value for modern systematics and evolutionary biology, and that a continued emphasis on expanding taxonomic and genomic coverage will play a critical role in developing these resources to their full potential.     

利用肠道菌群的分布和16S DNA序列研究鲤科肠道菌系统演化关系

肠道微生物与寄主具有复杂的、多方面的相互依存效应,这种依存效应所产生的共生关系或协同进化关系既可反映寄主间的系统演化关系,也可显示肠道微生物间的系统演化关系,共生关系或协同进化关系是由于寄主与肠道微生物两者之间存在着相互自然选择作用所形成的,在长期的进化历程中逐步发生的共生关系信息很可能被记录在DNA序列中。本文通过检测鱼鲤鱼科8种鱼中9种肠道菌群的分布含量对这9种菌群进行分析,且利用从GenBank调取这9种肠道细菌菌属的43个种或亚种的16S DNA序列的构建NJ树和MP树,将这6个科9个属43个种或亚种分为革兰氏阴性和革兰氏阳性两大类群（一级分枝）。在这两类群中,又以科为单位分为6个亚类群（二级分枝）,而肠杆菌科中则以属为单位分为4个小类群（三级分枝）,此外球状菌与杆状菌也能截然分开。将16S DNA的NJ树隐去所有的种,以属为单位所得到的以分枝形式的无根树在拓扑结构上与菌群分布含量（寄主范围）所构建的无根树相近,但芽孢杆菌在两种无根树的位置中有较大的差异。如果提高检测水平,扩大所检测的寄主对象,这种差异有可能消除。     

The Dawn of Open Access to Phylogenetic Data

The scientific enterprise depends critically on the preservation of and open access to published data. This basic tenet applies acutely to phylogenies (estimates of evolutionary relationships among species). Increasingly, phylogenies are estimated from increasingly large, genome-scale datasets using increasingly complex statistical methods that require increasing levels of expertise and computational investment. Moreover, the resulting phylogenetic data provide an explicit historical perspective that critically informs research in a vast and growing number of scientific disciplines. One such use is the study of changes in rates of lineage diversification (speciation – extinction) through time. As part of a meta-analysis in this area, we sought to collect phylogenetic data (comprising nucleotide sequence alignment and tree files) from 217 studies published in 46 journals over a 13-year period. We document our attempts to procure those data (from online archives and by direct request to corresponding authors), and report results of analyses (using Bayesian logistic regression) to assess the impact of various factors on the success of our efforts. Overall, complete phylogenetic data for of these studies are effectively lost to science. Our study indicates that phylogenetic data are more likely to be deposited in online archives and/or shared upon request when: (1) the publishing journal has a strong data-sharing policy; (2) the publishing journal has a higher impact factor, and; (3) the data are requested from faculty rather than students. Importantly, our survey spans recent policy initiatives and infrastructural changes; our analyses indicate that the positive impact of these community initiatives has been both dramatic and immediate. Although the results of our study indicate that the situation is dire, our findings also reveal tremendous recent progress in the sharing and preservation of phylogenetic data.     

Analysis of Comparative Sequence and Genomic Data to Verify Phylogenetic Relationship and Explore a New Subfamily of Bacterial Lipases

Thermostable and organic solvent-tolerant enzymes have significant potential in a wide range of synthetic reactions in industry due to their inherent stability at high temperatures and their ability to endure harsh organic solvents. In this study, a novel gene encoding a true lipase was isolated by construction of a genomic DNA library of thermophilic Aneurinibacillus thermoaerophilus strain HZ into Escherichia coli plasmid vector. Sequence analysis revealed that HZ lipase had 62% identity to putative lipase from Bacillus pseudomycoides. The closely characterized lipases to the HZ lipase gene are from thermostable Bacillus and Geobacillus lipases belonging to the subfamily I.5 with ≤ 57% identity. The amino acid sequence analysis of HZ lipase determined a conserved pentapeptide containing the active serine, GHSMG and a Ca²⁺-binding motif, GCYGSD in the enzyme. Protein structure modeling showed that HZ lipase consisted of an α/β hydrolase fold and a lid domain. Protein sequence alignment, conserved regions analysis, clustal distance matrix and amino acid composition illustrated differences between HZ lipase and other thermostable lipases. Phylogenetic analysis revealed that this lipase represented a new subfamily of family I of bacterial true lipases, classified as family I.9. The HZ lipase was expressed under promoter P_lac using IPTG and was characterized. The recombinant enzyme showed optimal activity at 65°C and retained ≥ 97% activity after incubation at 50°C for 1h. The HZ lipase was stable in various polar and non-polar organic solvents.     

Gene Conversion at the Gray Locus of SORDARIA FIMICOLA: Fit of the Experimental Data to a Hybrid DNA Model of Recombination

A hybrid DNA (hDNA) model of recombination has been algebraically formulated, which allows the prediction of frequencies of postmeiotic segregation and conversion of a given allele and their probability of being associated with a crossing over. The model considered is essentially the "Aviemore model." In contrast to some other interpretations of recombination, it states that gene conversion can only result from the repair of heteroduplex hDNA, with postmeiotic segregation resulting from unrepaired heteroduplexes. The model also postulates that crossing over always occurs distally to the initiation site of the hDNA. Eleven types of conversion and postmeiotic segregation with or without associated crossover were considered. Their theoretical frequencies are given by 11 linear equations with ten variables, four describing heteroduplex repair, four giving the probability of hDNA formation and its topological properties and two giving the probability that crossing over occurs at the left or right of the converting allele. Using the experimental data of Kitani and coworkers on conversion at the six best studied gray alleles of Sordaria fimicola, we found that the model considered fit the data at a P level above or very close (allele h4) to the 5% level of sampling error provided that the hDNA is partly asymmetric. The best fitting solutions are such that the hDNA has an equal probability of being formed on either chromatid or, alternatively, that both DNA strands have the same probability of acting as the invading strand during hDNA formation. The two mismatches corresponding to a given allele are repaired with different efficiencies. Optimal solutions are found if one allows for repair to be more efficient on the asymmetric hDNA than on the symmetric one. In the case of allele g1, our data imply that the direction of repair is nonrandom with respect to the strand on which it occurs.     

The Evolutionary Gain of Spliceosomal Introns: Sequence and Phase Preferences

doi:10.1093/molbev/msh201 Wei-Gang Qiu, Nick     

An Arsenic Exposure Model: Probabilistic Validation Using Empirical Data

This paper describes development of a multi-pathway arsenic exposure model. The model uses information on arsenic concentrations in food, water, soil, and dust, combined with estimates of intake and medium-specific absorption. Urinary arsenic is predicted assuming that 60% of absorbed arsenic is excreted in urine under steady state conditions. Fecal arsenic is predicted assuming all unabsorbed arsenic is excreted in feces. We applied this model at a former copper smelter site. Site specific distributions were available for the following parameters: soil and dust arsenic concentration (geometric mean approximately 100 to 200?ppm and 50 to 100?ppm, respectively); the combined childhood soil and dust ingestion rate (geometric mean of 20?mg/d); soil and dust arsenic relative bioavailability (geometric mean 0.20 and 0.28, respectively); exposure duration; water arsenic concentration; air arsenic concentration; and total arsenic in food. Monte Carlo simulation was used to predict daily arsenic uptake and excretion in urine and feces for children. Predicted urine arsenic levels were less than measured levels (73% to 88% of measured values, depending on region of site). On the other hand, predicted fecal arsenic levels exceeded measured levels by a factor of 1.7 to 4.6. We were able to improve the correspondence between predicted and measured arsenic excretion rates by decreasing the assumed value of the combined soil and dust ingestion rate, and increasing the assumed bioavailability of arsenic in soil and dust.     

A Phylogenetic Analysis of Heterorhabditis (Nemata: Rhabditidae) Based on Internal Transcribed Spacer 1 DNA Sequence Data

Internal transcribed spacer 1 sequences were used to infer phylogenetic relationships among 8 of the 9 described species and one putative species of the entomopathogenic nematode genus Heterorhabditis. Sequences were aligned and optimized based on pairwise genetic distance and parsimony criteria and subjected to a variety of sequence alignment parameters. Phylogenetic trees were constructed with maximum parsimony, cladistic, distance, and maximum likelihood algorithms. Our results gave strong support for four pairs of sister species, while relationships between these pairs also were resolved but less well supported. The ITS1 region of the nuclear ribosomal repeat was a reliable source of homologous characters for resolving relationships between closely related taxa but provided more tenuous resolution among more divergent lineages. A high degree of sequence identity and lack of autapomorphic characters suggest that sister species pairs within three distinct lineages may be mutually conspecific. Application of these molecular data and current morphological knowledge to the delimitation of species is hindered by an incomplete understanding of their variability in natural populations.