贝叶斯支端定年法推断分异时间和演化速率

贝叶斯支端定年法是近些年开发的推断类群分异时间和演化速率的方法。它克服了传统分步计算的缺陷,但涉及的统计学知识也更多。本文从贝叶斯统计计算的角度分层剖析了支端定年法的原理和计算过程,按照分异时间的先验分布、演化速率的先验分布、特征状态变化的模型和马氏链蒙特卡罗算法几个部分,叙述并讨论了定年计算中的主要模型和算法。旨在一定程度上为古生物学家分析实际数据提供参考。     

Using zebrafish to investigate cypriniform evolutionary novelties: functional development and evolutionary diversification of the kinethmoid

Although the zebrafish has become a popular model organism for biomedical studies, we propose that the wealth of morphological novelties that characterize this cypriniform fish makes it well suited for investigating the development of evolutionary innovations. Morphological novelties associated with feeding in cypriniform fishes include: a unique structure of the pharyngeal jaws in which the lower pharyngeal jaws are enlarged and opposed to a pad on the basioccipital process; a palatal organ found on the roof of the buccal chamber that is thought to help process detrital food within the buccal chamber; and, the kinethmoid, a novel ossification that effects a unique means of premaxillary protrusion. We present new morphological and developmental data and review functional data regarding the role of the kinethmoid in premaxillary protrusion in the zebrafish. Premaxillary protrusion plays an important role in effective prey acquisition in teleosts and the evolution of a unique means of premaxillary protrusion within Cypriniformes may have led to a number of trophic radiations within this clade. Ontogenetic data from zebrafish show that substantial premaxillary protrusion is not seen until these fish have undergone metamorphosis at which point the adductor mandibulae musculature becomes divided and all ligamentous attachments become established. A comparative study of families within Cypriniformes shows diverse morphologies of the kinethmoid. The morphological diversification that characterizes the kinethmoid suggests that this feeding structure has played a role in trophic radiations within Cypriniformes, since the morphology of this feature is correlated with feeding habits.     

Nonlinear relationships among evolutionary rates identify regions of functional divergence in heat-shock protein 70 genes

The neutral theory predicts that, in comparisons among related genes, the number of amino acid replacements per site in a given gene region should be a linear function of that in another region of the same gene, unless the genes have diverged functionally in one region. Therefore, nonlinearity of this relationship can be used to identify regions of possible functional divergence among members of a multigene family. This method of analysis was applied to members of the heat-shock protein 70 (HSP70) gene family, which encode highly conserved ATP- dependent chaperone proteins found in all organisms. A nonlinear relationship was found between the rate of amino acid replacement in the conserved IA domain of the ATPase portion of the molecule and that in other ATPase domains and the peptide-binding domain. These results suggest that genes in the HSP70 subfamily C (dnaK of bacteria and SSC1 of yeast) may have diverged functionally from other subfamilies in the ATPase domains, especially IIB, whereas SSB1 of yeast has diverged markedly in the peptide-binding domain. Functional divergence within these regions is consistent with what is known about functional differences between the HSP70 subfamilies in yeast.      

Magnetic resonance study of exchangeable protons in human carbonic anhydrases.

A titratable exchangeable proton resonance assignable to a histidine imidazole ring N--H proton is observed approximately minus 15 ppm downfield from tetramethylsilane. The chemical shift of this resonance is affected by sulfonamide and anion inhibitors, and by removal of zinc or replacement of zinc by cobalt, indicating that the proton is located at or near the active site. The pH dependence of the chemical shift of this resonance, which is abolished by inhibitors, reflects the titration of a group with a pK-a of 7.3 in human carbonic anhydrase B and smaller than or equal to 7.1 in human carbonic anhydrase C. These pK-a values are interpreted to be due to the ionization of a neutral imidazole to form the imidazolate anion coordinated to zinc. A mechanism for enzymatic catalysis involving reversible deprotonation and coordination of a histidine to the metal is consistent with these studies.     

Xenopus laevis integrins. Structural conservation and evolutionary divergence of integrin beta subunits

We report the sequences of cDNA clones for two different integrin beta subunits isolated from a Xenopus laevis neurula cDNA library. mRNAs corresponding to both genes are first detected at gastrulation. We show that these two beta subunits are very highly related (98% identity in amino acid sequence) and probably arose at the time of tetraploidization of the X. laevis genome around 50 million years ago. Comparison of these sequences with those of various other vertebrate integrin beta subunit establishes that all species analyzed to date contain a highly conserved integrin beta subunit (beta 1). The interspecies homologies within this class of integrin beta subunits (82-86% identity in amino acid sequence) are much greater than those among the three different beta subunits which are known in humans (40-48% identity in amino acid sequence). Analysis of the homologies clearly indicates duplication and divergence of this multigene family more than 500 million years ago prior to the appearance of the vertebrates. We also observe cross-hybridization between cDNA probes for chicken integrin beta subunits and genomic DNAs of several invertebrate species. Despite the divergence in sequence among different integrin beta subunits, certain features of their structure are remarkably conserved.     

Role of arginine 59 in the gamma-class carbonic anhydrases.

The functional role of the highly conserved active site Arg 59 in the prototype of the gamma-class carbonic anhydrase Cam (carbonic anhydrase from Methanosarcina thermophila) was investigated. Variants (R59A, -C, -E, -H, -K, -M, and -Q) were prepared by site-directed mutagenesis and characterized by size exclusion chromatography (SEC), circular dichroism (CD) spectroscopy, and stopped-flow kinetic analyses. CD spectra indicated similar secondary structures for the wild type and the R59A and -K variants, independent of nondenaturing concentrations of guanidine hydrochloride (GdnHCl). SEC indicated that all variants purified as homotrimers like the wild type. SEC also revealed that the R59A and -K variants unfolded at > or = 1.5 M GdnHCl, compared to 3.0 M GdnHCl for the wild type. These results indicate that Arg 59 contributes to the thermodynamic stability of the Cam trimer. The R59K variant had k(cat) and k(cat)/K(m) values that were 8 and 5% of the wild-type values, respectively, while all other variants had k(cat) and k(cat)/K(m) values 10-100-fold lower than those of the wild type. The R59A, -C, -E, -M, and -Q variants exhibited 4-63-fold increases in k(cat) and 9-120-fold increases in k(cat)/K(m) upon addition of 100 mM GdnHCl, with the largest increases observed for the R59A variant, which was comparable to the R59K variant. The kinetic results indicate that a positive charge at position 59 is essential for the CO(2) hydration step of the overall catalytic mechanism.     

Characterization and function of carbonic anhydrases in the zooxanthellae-giant clam symbiosis.

Carbonic anhydrase (CA) has been purified from the host tissue of Tridacna gigas, a clam that lives in symbiosis with the dinoflagellate alga, Symbiodinium. At least two isoforms of CA were identified in both gill and mantle tissue. The larger (70 kDa) isoform is a glycoprotein with both N- and O-glycans attached and has highest homology to CAII. It is associated with the membrane fraction while the smaller (32 kDa) is present in the aqueous phase in both tissues. The 32 kDa CA has high homology with mammalian CAI at the N-terminus. Both isoforms cross-reacted with antibodies to CAII from chicken. Immunohistology demonstrated that the 70 kDa CA is present within the ciliated branchial filaments and cells lining the tertiary water channels in the gills of T. gigas. This is consistent with a role in the transport of inorganic carbon (Ci) to the haemolymph and therefore supply of Ci to the zooxanthellae. CA was also detected in mantle epithelial cells where it may also contribute to Ci supply to the zooxanthellae. The hyaline body and nerve tissue in the mantle express the 70 kDa CA where it may be involved in light sensing and nervous transmission.     

Using functional groups to investigate community response to environmental changes: two grassland case studies

Functional groups (FG) are an useful generalization to investigate environmental change effects on biotic communities. Assigning species to FGs is a contextual task and carries an arbitrary element, regardless of whether the grouping is obtained a priori or by sophisticated numerical methods. Using two grassland community case studies, we show that even simple FG allocation based on growth form, architecture and longevity (plants and mosses), or foraging characteristics (above-ground invertebrates) can be useful to increase our understanding of community processes. For example, the sensitivity of organisms to climate change increases with trophic rank and is higher in disturbed than in undisturbed communities. Complexity of interaction webs (in terms of web connectance), however, is larger in undisturbed than in disturbed communities. A significant and important relationship is likely to exist between anthropogenic disturbances, community complexity and the ecosystem effects of climate change. Trophic interactions may be disrupted much easier by climate changes in disturbed than in undisturbed communities where complexity may be buffering these effects.     

Structural study of interaction between brinzolamide and dorzolamide inhibition of human carbonic anhydrases

Carbonic anhydrases (CAs, EC 4.2.1.1) are metalloenzymes that catalyze the reversible hydration of carbon dioxide and bicarbonate. Their pivotal role in metabolism, ubiquitous nature, and multiple isoforms (CA I–XIV) has made CAs an attractive drug target in clinical applications. The usefulness of CA inhibitors (CAIs) in the treatment of glaucoma and epilepsy are well documented. In addition several isoforms of CAs (namely, CA IX) also serve as biological markers for certain tumors, and therefore they have the potential for useful applications in the treatment of cancer. This is a structural study on the binding interactions of the widely used CA inhibitory drugs brinzolamide (marketed as Azopt®) and dorzolamide (marketed as Trusopt®) with CA II and a CA IX-mimic, which was created via site-directed mutagenesis of CA II cDNA such that the active site resembles that of CA IX. Also the inhibition of CA II and CA IX and molecular docking reveal brinzolamide to be a more potent inhibitor among the other catalytically active CA isoforms compared to dorzolamide. The structures show that the tail end of the sulfonamide inhibitor is critical in forming stabilizing interactions that influence tight binding; therefore, for future drug design it is the tail moiety that will ultimately determine isoform specificity.     

Inhibition of carbonic anhydrases I and II by N-unsubstituted carbamate esters.

We previously showed that the zinc metalloenzyme carbonic anhydrases (CA I and II isozymes) bind "neutral" amides and related compounds as anions through coordination of their deprotonated amide nitrogen to the active site zinc (Rogers, J. I., Mukherjee, J., and Khalifah, R. G. (1987) Biochemistry 26, 5672-5679). Urethan, the ethyl carbamate ester, was among such compounds. The present study was designed to test whether other N-unsubstituted carbamate esters of pharmacological interest (as sedatives, hypnotics, anxiolytics, and skeletal muscle relaxants) were capable of binding to CA in the same manner. We studied the interaction of human CA I and II with urethan, phenyl carbamate, ethinamate, meprobamate, and methocarbamol. Phenyl carbamate studies were greatly complicated by its uncatalyzed hydrolysis via an elimination mechanism to form cyanate, a powerful CA inhibitor. In general, the compounds display: 1) slow on-off inhibition binding kinetics in the seconds range, 2) maximal inhibitor affinity at alkaline pH, and 3) characteristic three-band visible spectra of their complexes with cobalt-substituted CA I. These properties are shared with the previously studied amide inhibitors and are taken as evidence that the deprotonated carbamate nitrogen coordinates to the active site metal ion. CA I appeared to bind carbamate esters more tightly than CA II, an unusual 1000-fold selectivity being seen in the case of methocarbamol. The inhibition by these drugs is not sufficiently strong to implicate CA I and II in their mechanism of action. However, it does suggest the possible existence of previously unsuspected similarities between binding to CA and to their physiological receptors or targets, particularly the involvement of zinc.     

Functional divergence prediction from evolutionary analysis: a case study of vertebrate hemoglobin

It is a central assumption of evolution that gene duplications provide the genetic raw material from which to create proteins with new functions. The increasing availability in multigene family sequences that has resulted from genome projects has inspired the creation of novel in silico approaches to predict details of protein function. The underlying principle of all such approaches is to compare the evolutionary properties of homologous sequence positions in paralogous proteins. It has been proposed that the positions that show switches in substitution rate over time-i.e., "heterotachous sites," are good indicators of functional divergence. Here, we analyzed the alpha and beta paralogous subunits of hemoglobin in search for such signatures. We found as many heterotachous sites in comparisons between groups of paralogous subunits (alpha/beta) as between orthologous ones (alpha/alpha, beta/beta). Thus, the importance of substitution rate shifts as predictors of specialization between protein subfamilies might be reconsidered. Instead, such shifts may reflect a more general process of protein evolution, consistent with the fact that they can be compatible with function conservation. As an alternative, we focused on those residues showing highly constrained states in two sequence groups, but different in each group, and we named them CBD (for "constant but different"). As opposed to heterotachous positions, CBD sites were markedly overrepresented in paralogous (alpha/beta) comparisons, as opposed to orthologous ones (alpha/alpha, beta/beta), identifying them as likely signatures of functional specialization between the two subunits. When superimposed onto the three-dimensional structure of hemoglobin, CBD positions consistently appeared to cluster preferentially on inter-subunit surfaces, two contact areas crucial to function in vertebrate tetrameric hemoglobin. The identification and analysis of CBD sites by complementing structural information with evolutionary data may represent a promising direction for future studies dealing with the functional characterization of a growing number of multigene families identified by complete genome analyses.     

分子绝对进化速率与物种分歧时间之间的定量关系

通过对美国国家生物技术信息中心数据库Gen Bank提供的一些蛋白质和核苷酸序列进行比对和分析,发现生物分子绝对进化速率k与进化时间或物种分歧时间t之间存在下列定量关系:lnk=-Ea/Rt+lnK0,式中Ea为位点突变活化能,k0为分子极限绝对进化速率,R为常数,并对其生物学意义进行了初步的探讨;数据分析还揭示出物种的分子极限绝对进化速率与进化时间或物种分歧时间之间也服从相似的定量公式,也就是说生物分子进化过程可能同时受到序列位点突变和控制物种分子极限绝对进化速率进化的两个"分子钟"作用,即存在"双重分子钟"现象。     

A viral sampling design for testing the molecular clock and for estimating evolutionary rates and divergence times

MOTIVATION: The high pace of viral sequence change means that variation in the times at which sequences are sampled can have a profound effect both on the ability to detect trends over time in evolutionary rates and on the power to reject the Molecular Clock Hypothesis (MCH). Trends in viral evolutionary rates are of particular interest because their detection may allow connections to be established between a patient's treatment or condition and the process of evolution. Variation in sequence isolation times also impacts the uncertainty associated with estimates of divergence times and evolutionary rates. Variation in isolation times can be intentionally adjusted to increase the power of hypothesis tests and to reduce the uncertainty of evolutionary parameter estimates, but this fact has received little previous attention. RESULTS: We provide approximations for the power to reject the MCH when the alternative is that rates change in a linear fashion over time and when the alternative is that rates differ randomly among branches. In addition, we approximate the standard deviation of estimated evolutionary rates and divergence times. We illustrate how these approximations can be exploited to determine which viral sample to sequence when samples representing different dates are available.