A Cladistic Analysis of Phenotypic Associations with Haplotypes Inferred from Restriction Endonuclease Mapping and DNA Sequence Data. III. Cladogram Estimation

We previously developed a cladistic approach to identify subsets of haplotypes defined by restriction endonuclease mapping or DNA sequencing that are associated with significant phenotypic deviations. Our approach was limited to segments of DNA in which little recombination occurs. In such cases, a cladogram can be constructed from the restriction site or sequence data that represents the evolutionary steps that interrelate the observed haplotypes. The cladogram is used to define a nested statistical design to identify mutational steps associated with significant phenotypic deviations. The central assumption behind this strategy is that any undetected mutation causing a phenotypic effect is embedded within the same evolutionary history that is represented by the cladogram. The power of this approach depends upon the confidence one has in the particular cladogram used to draw inferences. In this paper, we present a strategy for estimating the set of cladograms that are consistent with a particular sample of either restriction site or nucleotide sequence data and that includes the possibility of recombination. We first evaluate the limits of parsimony in constructing cladograms. Once these limits have been determined, we construct the set of parsimonious and nonparsimonious cladograms that is consistent with these limits. Our estimation procedure also identifies haplotypes that are candidates for being products of recombination. If recombination is extensive, our algorithm subdivides the DNA region into two or more subsections, each having little or no internal recombination. We apply this estimation procedure to three data sets to illustrate varying degrees of cladogram ambiguity and recombination.     

Phylogenetics of Freshwater Black Basses (Centrarchidae: Micropterus) Inferred from Restriction Endonuclease Analysis of Mitochondrial DNA

Geographic isolation and habitat specialization has aided in the evolution and genetic integrity of the micropterid bass species of North America. Members of the genus Micropterus form a close natural unit with little morphologic and meristic variation. Our goals were to measure the genetic characteristics of and distances between six black bass species by using mitochondrial DNA analysis. Mitochondrial DNA restriction fragment length polymorphisms were examined in Guadalupe bass (M. treculi), largemouth bass (M. salmoides), shoal bass (M. cataractae), smallmouth bass (M. dolomieu), spotted bass (M. punctulatus), and Suwannee bass (M. notius), using 15 restriction endonucleases. The bluegill (Lepomis macrochirus) was used as an outgroup. The phylogeny inferred from Dollo parsimony cladistic analysis concurred with published results from allozyme analyses, yet it was inconsistent with published meristic analyses. Genetic distances between species ranged from 0.0659 to 0.2145, with the largemouth and Suwannee basses showing the greatest divergence from the other black basses. The Guadalupe, smallmouth, and spotted basses were most diverged from the bluegill. The black basses diverged over a broad time frame, with estimated black bass speciation occurring during late Miocene-early Pliocene (3.30-10.73 MYA).     

Genetic Differentiation Inferred from Data on Restriction Analysis of Mitochondrial DNA in the Northern and Southern Forms of the Dolly Varden Char

Restriction fragment length polymorphism of the ATPase6/ND4L and cytochrome b regions of mitochondrial DNA (mtDNA) was studied in populations of the northern and southern forms of the Dolly Varden char Salvelinus malma from the Asian coast of the northern part of the Pacific Ocean. Seven genotypes of mtDNA were found, of which AAAA/ACCC and BBAC/ACEC characterize the northern and southern groups of populations, respectively. By the identity test, the studied populations of the Chukotka and Kamchatka peninsulas and the northern coast of the Sea of Okhotsk significantly differed from the studied populations of Sakhalin Island, Primorye, and Kuril Islands.     

Phylogenetic Relationships among Bufonoid Frogs (Anura: Neobatrachia) Inferred from Mitochondrial DNA Sequences

Nucleotide sequences of portions of the mitochondrial 12S and 16S ribosomal RNA genes were used to extend a recent study of anuran phylogeny (Hayet al., Mol. Biol. Evol.12: 928–937, 1995) and to further evaluate phylogenetic relationships within the Neobatrachia. An analysis of almost 900 nucleotides from each of 8 new representatives of the Dendrobatidae, Hylidae, Leptodactylidae, and Myobatrachidae, plus 14 available members of the Neobatrachia provides support for 2 major lineages (Bufonoidea and Ranoidea) within this anuran suborder. The neotropical Bufonoidea and their derivatives are monophyletic. There is an interesting association of the 2 Australian myobatrachids with the South African Heleophrynidae, and the Sooglossidae is one of the basal bufonoid lineages. Within the New World bufonoid frogs, a monophyletic Dendrobatidae is strongly supported. An Australian hylid (Pelodryadinae) shows close affinity with the South American hylid Phyllomedusinae. A group composed of Hylinae (HylaandSmilisca), Centrolenidae, Bufonidae, and the hylid Hemiphractinae, with the latter two clustered, was supported significantly. The addition of new taxa has more clearly defined some relationships within the suborder Neobatrachia and has indicated that the families Hylidae, Leptodactylidae, and Myobatrachidae may not be monophyletic.     

Phylogenetic Relations among Percid Fishes as Inferred from Mitochondrial CytochromebDNA Sequence Data

Hypotheses of relationship among genera of Percidae have been conflicting. Based on different phylogenetic premises, the evolution of small benthic forms in Percidae has been interpreted as resulting from either convergence or common ancestry. In order to assess various phylogenetic hypotheses of Percidae we collected complete sequences (1140 bp) of mitochondrially encoded cytochromebfor 21 species of percids. Seven species representing four additional families of Perciformes were used as outgroups. Maximum parsimony and minimum evolution analyses both recovered single shortest trees, and the results of these analyses were generally congruent with one another. All analyses consistently recovered three monophyletic groups in Percidae: Etheostomatinae (Ammocrypta, Crystallaria, Etheostoma,andPercina), Percinae (PercaandGymnocephalus), and Luciopercinae (Stizostedion, Zingel,andRomanichthys). As a result of this analysis we present a revised classification of Percidae and discuss the phylogenetic evidence for the independent evolution of small benthic species within Etheostomatinae and Luciopercinae.     

A Cladistic Analysis of Phenotype Associations with Haplotypes Inferred from Restriction Endonuclease Mapping. II. the Analysis of Natural Populations

Genes that code for products involved in the physiology of a phenotype are logical candidates for explaining interindividual variation in that phenotype. We present a methodology for discovering associations between genetic variation at such candidate loci (assayed through restriction endonuclease mapping) with phenotypic variation at the population level. We confine our analyses to DNA regions in which recombination is very rare. In this case, the genetic variation at the candiate locus can be organized into a cladogram that represents the evolutionary relationships between the observed haplotypes. Any mutation causing a significant phenotypic effect should be imbedded within the same historical structure defined by the cladogram. We showed, in the first paper of this series, how to use the cladogram to define a nested analysis of variance (NANOVA) that was very efficient at detecting and localizing phenotypically important mutations. However, the NANOVA of haplotype effects could only be applied to populations of homozygous genotypes. In this paper, we apply the quantitative genetic concept of average excess to evaluate the phenotypic effect of a haplotype or group of haplotypes stratified and contrasted according to the nested design defined by the cladogram. We also show how a permutational procedure can be used to make statistical inferences about the nested average excess values in populations containing heterozygous as well as homozygous genotypes. We provide two worked examples that investigate associations between genetic variation at or near the Alcohol dehydrogenase (Adh) locus and Adh activity in Drosophila melanogaster, and associations between genetic variation at or near some apolipoprotein loci and various lipid phenotypes in a human population.     

A Molecular Phylogeny of Hemiptera Inferred from Mitochondrial Genome Sequences

Classically, Hemiptera is comprised of two suborders: Homoptera and Heteroptera. Homoptera includes Cicadomorpha, Fulgoromorpha and Sternorrhyncha. However, according to previous molecular phylogenetic studies based on 18S rDNA, Fulgoromorpha has a closer relationship to Heteroptera than to other hemipterans, leaving Homoptera as paraphyletic. Therefore, the position of Fulgoromorpha is important for studying phylogenetic structure of Hemiptera. We inferred the evolutionary affiliations of twenty-five superfamilies of Hemiptera using mitochondrial protein-coding genes and rRNAs. We sequenced three mitogenomes, from Pyrops candelaria, Lycorma delicatula and Ricania marginalis, representing two additional families in Fulgoromorpha. Pyrops and Lycorma are representatives of an additional major family Fulgoridae in Fulgoromorpha, whereas Ricania is a second representative of the highly derived clade Ricaniidae. The organization and size of these mitogenomes are similar to those of the sequenced fulgoroid species. Our consensus phylogeny of Hemiptera largely supported the relationships (((Fulgoromorpha,Sternorrhyncha),Cicadomorpha),Heteroptera), and thus supported the classic phylogeny of Hemiptera. Selection of optimal evolutionary models (exclusion and inclusion of two rRNA genes or of third codon positions of protein-coding genes) demonstrated that rapidly evolving and saturated sites should be removed from the analyses.     

基于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个族间的系统关系,即紫茎族是最基部的分支,山茶族与大头茶族间有更近的亲缘关系。同时,本文认为,厚皮香（亚）科是否为单系类群值得进一步研究。     

Molecular Phylogenetics of Honey Bee Subspecies (Apis melliferaL.) Inferred from Mitochondrial DNA Sequence

A mitochondrial DNA region encompassing part of the NADH dehydrogenase subunit 2 and isoleucine transfer RNA genes was PCR amplified, cloned, and sequenced for 14 morphometrically identifiedApis melliferasubspecies and the New World “Africanized” honeybee. Twenty different haplotypes were detected and phylogenetic analyses supported the existence of 3 or 4 major subspecies groups similar to those based on morphometric measurements. However, some discrepancies are reported concerning the subspecies composition of each group. Based on the sequence divergence ofDrosophila(2% per Myr) we found that the four lineages may have diverged around 0.67 Myr. The variability found in this region enables us to infer phylogenetic relationships and test hypotheses concerning subspecies origin, dispersion, and biogeography.     

Differentiation and Genetic Position of Slavs among Eurasian Ethnic Groups as Inferred from Variation in Mitochondrial DNA

The distribution of identical and similar (phylogenetically related) types of hypervariable segment 1 (HVS1) of the mitochondrial DNA (mtDNA) was studied in human populations belonging to three Slavonic groups and nine ethnogeographic groups of Eurasia (total sample size 2772 people). The results testified to a common origin of West, South, and East Slavs and revealed a central place of West Slavs among all Slavonic ethnic groups. Mixing was shown to play a substantial role in the formation of specific features of all three Slavonic gene pools. The mitochondrial gene pools of the Slavonic ethnic groups proved to preserve features suggesting a common ancestor for these and South European populations (especially those of the Balkan Peninsula).     

A Cladistic Analysis of Phenotypic Associations with Haplotypes Inferred from Restriction Endonuclease Mapping. IV. Nested Analyses with Cladogram Uncertainty and Recombination

We previously developed an analytical strategy based on cladistic theory to identify subsets of haplotypes that are associated with significant phenotypic deviations. Our initial approach was limited to segments of DNA in which little recombination occurs. In such cases, a cladogram can be constructed from the restriction site data to estimate the evolutionary steps that interrelate the observed haplotypes to one another. The cladogram is then used to define a nested statistical design for identifying mutational steps associated with significant phenotypic deviations. The central assumption behind this strategy is that a mutation responsible for a particular phenotypic effect is embedded within the evolutionary history that is represented by the cladogram. The power of this approach depends on the accuracy of the cladogram in portraying the evolutionary history of the DNA region. This accuracy can be diminished both by recombination and by uncertainty in the estimated cladogram topology. In a previous paper, we presented an algorithm for estimating the set of likely cladograms and recombination events. In this paper we present an algorithm for defining a nested statistical design under cladogram uncertainty and recombination. Given the nested design, phenotypic associations can be examined using either a nested analysis of variance (for haploids or homozygous strains) or permutation testing (for outcrossed, diploid gene regions). In this paper we also extend this analytical strategy to include categorical phenotypes in addition to quantitative phenotypes. Some worked examples are presented using Drosophila data sets. These examples illustrate that having some recombination may actually enhance the biological inferences that may derived from a cladistic analysis. In particular, recombination can be used to assign a physical localization to a given subregion for mutations responsible for significant phenotypic effects.     

A Cladistic Analysis of Phenotypic Associations with Haplotypes Inferred from Restriction Endonuclease Mapping or DNA Sequencing. V. Analysis of Case/Control Sampling Designs: Alzheimer''s Disease and the Apoprotein E Locus

Present-day associations between haplotypes at a candidate locus and phenotypes exist when phenotypically important mutations occurred at some point during the evolution of the current array of genetic variation. A cladistic statistical design can be defined that focuses power by using the evolutionary history of the candidate DNA region. This paper shows how cladistic methodology is used for the analysis of case/control data, a common sampling design in genetic/disease association studies. A worked example is presented of the associations for sporadic early and late-onset forms of Alzheimer's disease with the 19q13.2 chromosomal region that includes the loci for apoproteins E, CI, and CII. This analysis confirms earlier reports of a strong association of the ApoE &4 allele with Alzheimer's disease but indicates that it is premature to condsider this association causal, particularly for early onset cases. Associations were also found with the &2 allele, as previously reported, and with the 1 allele at the ApoCI locus. However, this analysis indicates that it is inappropriate both statistically and medically to use single markers as risk predictors when haplotype data are available, even when the mutation leading to the marker is identified as having a strong phenotypic association.     

Phylogenetic Relationships of the Macaques (Cercopithecidae: Macaca), Inferred from Mitochondrial DNA Sequences

To study the phylogenetic relationships of the macaques, five gene fragments were sequenced from 40 individuals of eight species: Macaca mulatta, M. cyclopis, M. fascicularis, M. arctoides, M. assamensis, M. thibetana, M. silenus, and M. leonina. In addition, sequences of M. sylvanus were obtained from Genbank. A baboon was used as the outgroup. The phylogenetic trees were constructed using maximum-parsimony and Bayesian methods. Because five gene fragments were from the mitochondrial genome and were inherited as a single entity without recombination, we combined the five genes into a single analysis. The parsimony bootstrap proportions we obtained were higher than those from earlier studies based on the combined mtDNA dataset. Excluding M. arctoides, our results are generally consistent with the classification of Delson (1980). Our phylogenetic analyses agree with earlier studies suggesting that the mitochondrial lineages of M. arctoides share a close evolutionary relationship with the mitochondrial lineages of the fascicularis group of macaques (and M. fascicularis, specifically). M. mulatta (with respect to M. cyclopis), M. assamensis assamensis (with respect to M. thibetana), and M. leonina (with respect to M. silenus) are paraphyletic based on our analysis of mitochondrial genes.     

Phylogenetic Relationships of Bears (the Ursidae) Inferred from Mitochondrial DNA Sequences

The phylogenetic relationships among some bear species are still open questions. We present here mitochondrial DNA sequences of D-loop region, cytochrome b, 12S rRNA, tRNA^Pro, and tRNA^Thr genes from all bear species and the giant panda. A series of evolutionary trees with concordant topology has been derived based on the combined data set of all of the mitochondrial DNA sequences, which may have resolved the evolutionary relationships of all bear species: the ancestor of the spectacled bear diverged first, followed by the sloth bear; the brown bear and polar bear are sister taxa relative to the Asiatic black bear; the closest relative of the American black bear is the sun bear. Primers for forensic identification of the giant panda and bears are proposed. Analysis of these data, in combination with data from primates and antelopes, suggests that relative substitutional rates between different mitochondrial DNA regions may vary greatly among different taxa of the vertebrates.     

The Genetic Architecture of Quantitative Traits Cannot Be Inferred from Variance Component Analysis

Classical quantitative genetic analyses estimate additive and non-additive genetic and environmental components of variance from phenotypes of related individuals without knowing the identities of quantitative trait loci (QTLs). Many studies have found a large proportion of quantitative trait variation can be attributed to the additive genetic variance (V_A), providing the basis for claims that non-additive gene actions are unimportant. In this study, we show that arbitrarily defined parameterizations of genetic effects seemingly consistent with non-additive gene actions can also capture the majority of genetic variation. This reveals a logical flaw in using the relative magnitudes of variance components to indicate the relative importance of additive and non-additive gene actions. We discuss the implications and propose that variance component analyses should not be used to infer the genetic architecture of quantitative traits.     

Analysis of Mutation Mechanisms in Human Mitochondrial DNA

The cause of the high variability of human mitochondrial DNA (mtDNA) remains largely unknown. Three mechanisms of mutagenesis that might account for the generation of nucleotide substitutions in mtDNA have been analyzed: deamination of DNA nitrous bases caused by deamination agents, tautomeric proton migration in nitrous bases, and the hydrolysis of the glycoside bond between the nitrous base and carbohydrate residue in nucleotides against the background of the free-radical damage of DNA polymerase γ. Quantum chemical calculations demonstrated that the hydrolysis of the N-glycoside bond is the most probable mechanism; it is especially prominent in the H strand, which remains free during mtDNA replication for a relatively long time. It has also been found that hydrolytic deamination of adenine in single-stranded regions of the H strand is a possible cause of the high frequency of T → C transitions in the mutation spectra of the L-chain of the major mtDNA noncoding region.     

FMMU白化豚鼠线粒体DNA RFLP分析研究

目的研究FMMU白化豚鼠的mtDNA,并与花色豚鼠mtDNA进行多态性分析比较,以确定其独特的生物学特性是否与mtDNA相关。方法用碱变性法提取FMMU白化豚鼠以及花色豚鼠的mtDNA,并用AvaⅠ、BalⅠ等12种限制性内切酶进行酶切和限制性片段长度多态性分析。结果与结论FMMU白化豚鼠mtDNA和花色豚鼠mtDNA的相对分子质量相同,约为16.7×103;FMMU白化豚鼠与花色豚鼠两品系的mtDNA经AvaⅠ、BalⅠ等内切酶酶切后有3-8个酶切位点,酶切图谱完全相同,经RFLP分析FMMU白化豚鼠与花色豚鼠的mtDNA之间缺乏多态性。本实验没有发现FMMU白化豚鼠的独特的生物学特性与mtDNA相关。     

犬Ⅰ型腺病毒DNA的酶切分析及分子克隆

犬Ⅰ型腺病毒(CAV-1)弱毒用限制性内切酶EcoR Ⅰ,BamH Ⅰ,Pst Ⅰ,Sph Ⅰ和Hind Ⅲ消化分析后其图谱与强毒株相比没有差异。将弱毒DNA用Pst Ⅰ完全消化后以鸟枪法克隆到载体质粒pBluescrip'SK中,经用光生物素标记的CAV-1 DNA杂交筛选以及Pst Ⅰ分析重组质粒证明已将分子量为5.5,3.5,2.85,1.2,0.32和0.28Kb的CAV-1DNA片段克隆到质粒中。克隆到的这些片段将可考虑进一步研究作为探针检测犬及狐狸等野生动物的腺病毒感染。