Methods for computing the standard errors of branching points in an evolutionary tree and their application to molecular data from humans and apes

Statistical methods for computing the standard errors of the branching points of an evolutionary tree are developed. These methods are for the unweighted pair-group method-determined (UPGMA) trees reconstructed from molecular data such as amino acid sequences, nucleotide sequences, restriction-sites data, and electrophoretic distances. They were applied to data for the human, chimpanzee, gorilla, orangutan, and gibbon species. Among the four different sets of data used, DNA sequences for an 895-nucleotide segment of mitochondrial DNA (Brown et al. 1982) gave the most reliable tree, whereas electrophoretic data (Bruce and Ayala 1979) gave the least reliable one. The DNA sequence data suggested that the chimpanzee is the closest and that the gorilla is the next closest to the human species. The orangutan and gibbon are more distantly related to man than is the gorilla. This topology of the tree is in agreement with that for the tree obtained from chromosomal studies and DNA-hybridization experiments. However, the difference between the branching point for the human and the chimpanzee species and that for the gorilla species and the human-chimpanzee group is not statistically significant. In addition to this analysis, various factors that affect the accuracy of an estimated tree are discussed.      

A highly conserved nuclear gene for low-level phylogenetics: elongation factor-1 alpha recovers morphology-based tree for heliothine moths

Molecular systematists need increased access to nuclear genes. Highly conserved, low copy number protein-encoding nuclear genes have attractive features for phylogenetic inference but have heretofore been applied mostly to very ancient divergences. By virtue of their synonymous substitutions, such genes should contain a wealth of information about lower-level taxonomic relationships as well, with the advantage that amino acid conservatism makes both alignment and primer definition straightforward. We tested this postulate for the elongation factor-1 alpha (EF-1 alpha) gene in the noctuid moth subfamily Heliothinae, which has probably diversified since the middle Tertiary. We sequenced 1,240 bp in 18 taxa representing heliothine groupings strongly supported by previous morphological and allozyme studies. The single most parsimonious gene tree and the neighbor-joining tree for all nucleotides show almost complete concordance with the morphological tree. Homoplasy and pairwise divergence levels are low, transition/transversion ratios are high, and phylogenetic information is spread evenly across gene regions. The EF-1 alpha gene and presumably other highly conserved genes hold much promise for phylogenetics of Tertiary age eukaryote groups.      

Mesodern-inducing factors and mesodermal patterning

Identification of the signalling molecules involved in mesoderm formation in amphibian embryos still presents problems. None of the original candidates, such as activin, have been definitively ruled out, and the new factors, such as the nodal-related genes, have come on to the scene. Of the original candidates, activin has been definitively shown to act as a morphogen, whereas bone morphogenetic protein (BMP)-4 has emerged as a ventral inducer and an inhibitor of neural differentiation. The effects of BMP-4 are antagonized by chordin, a molecule related to the product of the Drosophila gene short gastrulation.     

Ant genomics sheds light on the molecular regulation of social organization

Ants are powerful model systems for the study of cooperation and sociality. In this review, we discuss how recent advances in ant genomics have contributed to our understanding of the evolution and organization of insect societies at the molecular level.     

A comparative description of mitochondrial DNA differentiation in selected avian and other vertebrate genera

Levels of mitochondrial DNA (mtDNA) sequence divergence between species within each of several avian (Anas, Aythya, Dendroica, Melospiza, and Zonotrichia) and nonavian (Lepomis and Hyla) vertebrate genera were compared. An analysis of digestion profiles generated by 13-18 restriction endonucleases indicates little overlap in magnitude of mtDNA divergence for the avian versus nonavian taxa examined. In 55 interspecific comparisons among the avian congeners, the fraction of identical fragment lengths (F) ranged from 0.26 to 0.96 (F = 0.46), and, given certain assumptions, these translate into estimates of nucleotide sequence divergence (p) ranging from 0.007 to 0.088; in 46 comparisons among the fish and amphibian congeners, F values ranged from 0.00 to 0.36 (F = 0.09), yielding estimates of P greater than 0.070. The small mtDNA distances among avian congeners are associated with protein-electrophoretic distances (D values) less than approximately 0.2, while the mtDNA distances among assayed fish and amphibian congeners are associated with D values usually greater than 0.4. Since the conservative pattern of protein differentiation previously reported for many avian versus nonavian taxa now appears to be paralleled by a conservative pattern of mtDNA divergence, it seems increasingly likely that many avian species have shared more recent common ancestors than have their nonavian taxonomic counterparts. However, estimates of avian divergence times derived from mtDNA- and protein-calibrated clocks cannot readily be reconciled with some published dates based on limited fossil remains. If the earlier paleontological interpretations are valid, then protein and mtDNA evolution must be somewhat decelerated in birds. The empirical and conceptual issues raised by these findings are highly analogous to those in the long-standing debate about rates of molecular evolution and times of separation of ancestral hominids from African apes.      

Equivalence of multibreed animal models and hierarchical Bayes analysis for maternally influenced traits

Background

It has been argued that multibreed animal models should include a heterogeneous covariance structure. However, the estimation of the (co)variance components is not an easy task, because these parameters can not be factored out from the inverse of the additive genetic covariance matrix. An alternative model, based on the decomposition of the genetic covariance matrix by source of variability, provides a much simpler formulation. In this study, we formalize the equivalence between this alternative model and the one derived from the quantitative genetic theory. Further, we extend the model to include maternal effects and, in order to estimate the (co)variance components, we describe a hierarchical Bayes implementation. Finally, we implement the model to weaning weight data from an Angus × Hereford crossbred experiment.

Methods

Our argument is based on redefining the vectors of breeding values by breed origin such that they do not include individuals with null contributions. Next, we define matrices that retrieve the null-row and the null-column pattern and, by means of appropriate algebraic operations, we demonstrate the equivalence. The extension to include maternal effects and the estimation of the (co)variance components through the hierarchical Bayes analysis are then straightforward. A FORTRAN 90 Gibbs sampler was specifically programmed and executed to estimate the (co)variance components of the Angus × Hereford population.

Results

In general, genetic (co)variance components showed marginal posterior densities with a high degree of symmetry, except for the segregation components. Angus and Hereford breeds contributed with 50.26% and 41.73% of the total direct additive variance, and with 23.59% and 59.65% of the total maternal additive variance. In turn, the contribution of the segregation variance was not significant in either case, which suggests that the allelic frequencies in the two parental breeds were similar.

Conclusion

The multibreed maternal animal model introduced in this study simplifies the problem of estimating (co)variance components in the framework of a hierarchical Bayes analysis. Using this approach, we obtained for the first time estimates of the full set of genetic (co)variance components. It would be interesting to assess the performance of the procedure with field data, especially when interbreed information is limited.