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.     

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 cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping. I. Basic theory and an analysis of alcohol dehydrogenase activity in Drosophila

Because some genes have been cloned that have a known biochemical or physiological function, genetic variation can be measured in a population at loci that may directly influence a phenotype of interest. With this measured genotype approach, specific alleles or haplotypes in the probed DNA region can be assigned phenotypic effects. In this paper we address several problems encountered in implementing the measured genotype approach with restriction site data. A number of analytical problems arise in part as a consequence of the linkage disequilibrium that is commonly encountered when dealing with small DNA regions: 1) different restriction site polymorphisms are not statistically independent, 2) the sites being measured are not likely to be the direct cause of the associated phenotypic effects, 3) haplotype classes may be phenotypically heterogeneous, and 4) the sites that are most strongly associated with phenotypic effects are not necessarily the most closely linked to the actual genetic cause of the effects. When recombination and gene conversion are rare, the primary cause of linkage disequilibrium is history (mutational origin, genetic drift, hitchhiking, etc.). We deal with historical association directly by producing a cladogram that partially reconstructs the evolutionary history of the present-day haplotype variability. The cladogram defines a nested analysis of variance that simultaneously detects phenotypic effects, localizes the effects within the cladogram, and identifies haplotypes that are potentially heterogeneous in their phenotypic associations. The power of this approach is illustrated by an analysis of the associations between alcohol dehydrogenase (ADH) activity and restriction site variability in a 13-kb fragment surrounding the ADH locus in Drosophila melanogaster.     

Separating Population Structure from Population History: A Cladistic Analysis of the Geographical Distribution of Mitochondrial DNA Haplotypes in the Tiger Salamander, Ambystoma Tigrinum

Nonrandom associations of alleles or haplotypes with geographical location can arise from restricted gene flow, historical events (fragmentation, range expansion, colonization), or any mixture of these factors. In this paper, we show how a nested cladistic analysis of geographical distances can be used to test the null hypothesis of no geographical association of haplotypes, test the hypothesis that significant associations are due to restricted gene flow, and identify patterns of significant association that are due to historical events. In this last case, criteria are given to discriminate among contiguous range expansion, long-distance colonization, and population fragmentation. The ability to make these discriminations depends critically upon an adequate geographical sampling design. These points are illustrated with a worked example: mitochondrial DNA haplotypes in the salamander Ambystoma tigrinum. For this example, prior information exists about restricted gene flow and likely historical events, and the nested cladistic analyses were completely concordant with this prior information. This concordance establishes the plausibility of this nested cladistic approach, but much future work will be necessary to demonstrate robustness and to explore the power and accuracy of this procedure.     

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.     

A microbial population-species interface: nested cladistic and coalescent inference with multilocus data

Using sequence data from seven nuclear loci in 385 isolates of the haploid, plant parasitic, ascomycete fungus, Sclerotinia, divergence times of populations and of species were distinguished. The evolutionary history of haplotypes on both population and species scales was reconstructed using a combination of parsimony, maximum likelihood and coalescent methods, implemented in a specific order. Analysis of site compatibility revealed recombination blocks from which alternative (marginal) networks were inferred, reducing uncertainty in the network due to recombination. Our own modifications of Templeton and co-workers' cladistic inference method and a coalescent approach detected the same phylogeographic processes. Assuming neutrality and a molecular clock, the boundary between divergent populations and species is an interval of time between coalescence (to a common ancestor) of populations and coalescence of species.     

Tree scanning: a method for using haplotype trees in phenotype/genotype association studies

We use evolutionary trees of haplotypes to study phenotypic associations by exhaustively examining all possible biallelic partitions of the tree, a technique we call tree scanning. If the first scan detects significant associations, additional rounds of tree scanning are used to partition the tree into three or more allelic classes. Two worked examples are presented. The first is a reanalysis of associations between haplotypes at the Alcohol Dehydrogenase locus in Drosophila melanogaster that was previously analyzed using a nested clade analysis, a more complicated technique for using haplotype trees to detect phenotypic associations. Tree scanning and the nested clade analysis yield the same inferences when permutation testing is used with both approaches. The second example is an analysis of associations between variation in various lipid traits and genetic variation at the Apolipoprotein E (APOE) gene in three human populations. Tree scanning successfully identified phenotypic associations expected from previous analyses. Tree scanning for the most part detected more associations and provided a better biological interpretative framework than single SNP analyses. We also show how prior information can be incorporated into the tree scan by starting with the traditional three electrophoretic alleles at APOE. Tree scanning detected genetically determined phenotypic heterogeneity within all three electrophoretic allelic classes. Overall, tree scanning is a simple, powerful, and flexible method for using haplotype trees to detect phenotype/genotype associations at candidate loci.     

A Review on Cladistics

The theoretical bases and approaches of cladistics and some specific problems that, directly or indirectly, rely on cladistic analysis for their revolution, are outlined and discussed. Seven sections comprise this paper: a ) the philosophical foundation of cladistics; b) the theoretical tenets of cladistics; c) the operational procedure of cladisties; d) three schools of classification; e) cladistics and biogeography; f) cladistics and hybrid recognition; and g) is cladistic systematics a scientific theory ? Considerations of scientific methodology involve philosophical questions. From this point, Popper'falsificationism serves a good foundation. Popper emphasizes that all scientific knowledge is hypothetical-deductive, consisting of general statements (theories) that can never be confirmed or verified but only falsified. The theories, that can be tested most effectively, are preferable. Cladistics, aiming at generating accurately expressed and strictly testable systematic hypotheses, is well compatible with this requirement. The principles central to the cladistic theory and methodology are: the Principle of Synapomorphy; the Principle of Strict Monophyly; and the Principle of Strict Parsimony. The first requires forming nested groups by nesting statements about shared evolutionary novelties (synapomorphy) postulated from observed similarities and is the primary one. The second is mainly methodological, subject to modification and compromise. The principle of strict parsimony specifies the most preferable hypothesis (namely the one exhibiting the most congruence in the synapomorphy pattern). The operational procedure that might be followed in formulating and testing hypotheses of the synapomorphy pattern (the cladogram itself) consists of five steps. The erections of monophyletic groups, to a greater or lesser extent, rely on the hypothesis of the previous systematic studies and is the starting point for cladistic analysis. Character analysis, which focuses on character distribution and determination of the polarities, decides the reconstructed phylogeny. A detailed discussion on the methodological principles for identifying transformation sequence is presented. Many algorithms have been designated to infer the cladogram, and are basically of parsimony techniques and Compatibility techiques. The thus yielded cladograms, with their expected pattern of congruent synapomorphies, are tests of a particular hypothesis of synapomorphy and reciprocally synapomorphies are tests of cladistic hypothesis (cladogram). Such reciprocity is a strong stimulus to profound understanding on phylogenetic process and phyletic relationships. The cladogram and the Linnaean classification have the identical logic structure and the set-membership of the two can be made isomorphic. There are three principal approaches to biological classification : cladistics, phenetics and evolutionary classification. Cladistics is the determination of the branching pattern of evolution, and in the context of classification, the development of nested sets based on cladograms. Phenetics is the classification by overall similarities, without regard to evolutionary considerations. Evolutionary classification attempts to consider all meaningful aspects of phylogeny and to use these for making a classification. The last approach has been done intuitively, without explicit methods. An enumeration of their differences and a discussion on their relative merits are presented. Three theoretical approaches have been proposed for interpreting biogeographical history: the phylogenetic theory of biogeography, classical evolutionary biogeography and vicariance biogeography. The former two show some similarities in that they usually look upon biogeography in terms of centers of origin and dispersal from the centers. But the first puts a strong emphasis on the construction of hypotheses about the phylogenetic relationships of the organisms in question and the subsequent inference of their geographic relationships; the second advocates a theory which does not have a precise deductive link with phylogenetic construction and often results in wildly narratative-type hypotheses. The vicariance approach de-emphasizes the concepts of centers of origin and dispersal and attempts to analyse distribution patterns in terms of subdivision (vicariance) of ancestral biotas. The development of the theory of plate tectonics and its universal acceptance enormously stimulate biogeographers to look at the world's continents and oceans from a mobilist point, which, along with the establishment of the rigorous tool of the phylogenetic analysis (cladistics), profoundly reshapes the above three theories. Hybridization and polyploidy are outstanding features of many plant groups. But hybridization, or reticulate evolution, is inconsistent with the basic concepts of cladistics which is an ever-branching pattern. Cladists have suggested several approaches. One of them analyses all the taxa by a standard cladistic procedure and closely examines the cladograms for polytomies and character conflicts that may indicate possible hybrids. Such generated hypothesis of hybridization can be corroborated or falsified by other forms of data, such as distribution, polyploidy, karyotype and pollen fertility. There are three criteria to justify a theory to be scientific: a) whether it is a theory composed of hypotheses strictly falsifiable; b) whether it has predictive effect; and c) whether it has a explanatory value. Cladistic systematics aims at generating cladograms, which are hypotheses of the nested pattern of synapomorphy, phylogenetic process and phyletic relationships, susceptible to testing by postulated synapomorphies. The predictive effect of systematics relies on the acceptance of hypotheses of congruence about the correlation of characters, which has been well founded. For non-systematic biologists, phylogenetic classification can be used as axiom to form a preliminary and fundamental explanation.     

An Algorithm for Cladistic Taxonomy—Method of Maximal Same Step Length

This paper deals with the numerical cladistic taxonomy. A method for constructing evolutionary tree (method of maximal same step length) is proposed in the applications and practice of cladistic taxonomy. Its algorithm runs as follows: 1) According to the order of evolution, characters are coded with nonnegative integers, producing the original data matrix. 2) Calculate the same step coefficients Sij (i≠j) by the formula (3) and form the coefficient matrix. 3. Find the maximal value S_pq of the same step coefficients in the coefficient matrix. 4) According to the maximal same step length S_pq, the most recent common ancestor CTU, of CTUp and CTUq can be determinated by (2). 5) draw the cladistic edges of cladogram representing the evolutionary relationship from OTUT to OTUp and OTUq. If the number of CTUs in the data matrix≤2, go to (2), otherwise stop. An example of 6 species from the family Campanulaceae is given for illustration (See Table 1). In general case, the evolutionary length of the cladogram obtained by this method is shorter than that by monothetic and other methods. Its algorithm is easily performed and is especially suitable for computerizing.相似文献   

Information obtained in cladistic analysis.

The application of elementary equations from information theory to the elements involved in cladistic analysis is formalized mathematically. An equation is derived that quantifies the amount of information obtained by constructing a cladogram from a cladistic data matrix. Given particular conditions, the amount of information obtained increases monotonically with increases of the number of taxa involved and, so, may be used directly as a comparative measure of species richness for sister groups; in general, however, the amount of information obtained is related to the distribution of character states on the cladogram(s) deduced. An example is presented in which clades representing 11 phyla in the animal kingdom are compared in terms of information yielded. The amount of information obtained is consistent for different numbers of taxa and characters used in classifications. Speculative evolutionary explanations are presented for differences of information yielded among the phyla analyzed.     

The Role of Functional Analysis in Phylogenetic Inference: Examples from the History of the Xiphosura

SYNOPSIS. Conventional cladistic analyses of phylogeny can beinterpreted as operating at the level of phylogenetic trees.They assume that all "evolutionary steps" (transitions fromone character state to the next, along a morphocline) are independentand equal, and, on that basis, select the cladogram which isconsistent with the most parsimonious trees. Evaluation of theassumptions of independence and equality requires considerationof hypotheses at the levelof scenarios. In some cases, argumentsbased on functional analysis can suggest revised interpretationsof either homology or polarity. If properly formulated, thesearguments can alter the evaluation of parsimony for trees tothe extent that even the choice of cladogram is affected. Thestructure of scenario level arguments is identical to that ofarguments operating at tree level. Examples of phylogeneticinference in the context of xiphosurans (horseshoe crabs), usingboth comparative morphological and functional analysis, illustratethis approach. In different cases, orthodox interpretationsof relationship are either challenged or corroborated. Althoughthe introduction of functional analysis into the process ofphylogenetic inference may appear to compromise the usefulnessof the reconstructed phylogeny for testing hypotheses concerningthe role of natural selection in evolution, it actually increasesthe strength of such tests.     

INTERRELATIONSHIPS OF THE SPANISH MACKERELS (PISCES: SCOMBRIDAE: SCOMBEROMORUS) AND THEIR COPEPOD PARASITES

Abstract— A cladistic analysis of the eighteen species of Spanish mackerels (Scomberomorus) was conducted using the double-lined mackerel (Grammatorcynus) as the outgroup. Based on fifty-eight osteological and morphological characters, six species groups are recognized. Comparisons of the cladogram for the regalis species group were made with a cladogram of the nine species of parasitic copepods that infest the six members of the regalis group. Incongruences were analyzed. The parasite tree was then "forced" onto the host cladogram to determine which evolutionary events of the parasites could be explained by evolutionary events of the hosts. Hypotheses of coevolution are supported in some cases, hypotheses of dispersal are proposed in others.     

Using Haplotype Trees for Phylogeographic and Species Inference in Fish Populations

Genetic variation is now routinely screened at the DNA sequence level in many studies. If the DNA region being screened has not experienced excessive amounts of recombination, it is often possible to reconstruct the evolutionary history of the genetic variation in the form of a haplotype tree. This tree estimates the evolutionary pathway that interconnects all the different haplotypes (sequence variants) observed in the sample. This haplotype tree can be used to define a series of nested branches (clades) that reflects the relative temporal history of the haplotypes and groups of haplotypes. Geographical information can then be overlaid upon this temporal series to test for significant associations between geography and temporal position in the haplotype tree. This allows a reconstruction of how the genetic variation arose and spread in both space and time. Such reconstructions can yield many insights into the joint roles of recurrent events such as gene flow and of historical events such as fragmentation or range expansion. These points are illustrated with studies on the chub, Leuciscus cephalus. There is also a need to extend such nested phylogeographic analyses to a phylo/reticulate geographic analysis that incorporates both assortment and recombination between and within DNA regions. A preliminary phylo/reticulate geographic analysis is presented of the transferrin locus in the brown trout, Salmo trutta, species complex that reveals the importance of hybridization in the recent evolutionary history of this group. This example shows the inadequacy of a strictly phylogenetic approach and illustrates the need to incorporate reticulate evolution. The results of nested clade phylogeographic analysis and the new phylo/reticulate geographic analysis are then used for inferring species status of the marbled trout. The results indicate that an old hybridization event may have played a role in the origin of the marbled trout. Currently the marbled trout is primarily endangered by hybridization with introduced brown trout. These results show both the positive and negative impacts of hybridization upon biodiversity. Such phylo/reticulate geographic studies will challenge both our concepts of species and our conservation management strategies.     

Statistical power analysis of neutrality tests under demographic expansions, contractions and bottlenecks with recombination

Several tests have been proposed to detect departures of nucleotide variability patterns from neutral expectations. However, very different kinds of evolutionary processes, such as selective events or demographic changes, can produce similar deviations from these tests, thus making interpretation difficult when a significant departure of neutrality is detected. Here we study the effects of demography and recombination upon neutrality tests by analyzing their power under sudden population expansions, sudden contractions, and bottlenecks. We evaluate tests based on the frequency spectrum of mutations and the distribution of haplotypes and explore the consequences of using incorrect estimates of the rates of recombination when testing for neutrality. We show that tests that rely on haplotype frequencies-especially Fs and ZnS, which are based, respectively, on the number of different haplotypes and on the r2 values between all pairs of polymorphic sites-are the most powerful for detecting expansions on nonrecombining genomic regions. Nevertheless, they are strongly affected by misestimations of recombination, so they should not be used when recombination levels are unknown. Instead, class I tests, particularly Tajima's D or R2, are recommended.     

Nested clade analyses of phylogeographic data: testing hypotheses about gene flow and population history

Since the 1920s, population geneticists have had measures that describe how genetic variation is distributed spatially within a species' geographical range. Modern genetic survey techniques frequently yield information on the evolutionary relationships among the alleles or haplotypes as well as information on allele frequencies and their spatial distributions. This evolutionary information is often expressed in the form of an estimated haplotype or allele tree. Traditional statistics of population structure, such as F statistics, do not make use of evolutionary genealogical information, so it is necessary to develop new statistical estimators and tests that explicitly incorporate information from the haplotype tree. One such technique is to use the haplotype tree to define a nested series of branches (clades), thereby allowing an evolutionary nested analysis of the spatial distribution of genetic variation. Such a nested analysis can be performed regarding the geographical sampling locations either as categorical or continuous variables (i.e. some measure of spatial distance). It is shown that such nested phylogeographical analyses have more power to detect geographical associations than traditional, nonhistorical analyses and, as a consequence, allow a broader range of gene-flow parameters to be estimated in a precise fashion. More importantly, such nested analyses can discriminate between phylogeographical associations due to recurrent but restricted gene flow vs. historical events operating at the population level (e.g. past fragmentation, colonization, or range expansion events). Restricted gene flow and historical events can be intertwined, and the cladistic analyses can reconstruct their temporal juxtapositions, thereby yielding great insight into both the evolutionary history and population structure of the species. Examples are given that illustrate these properties, concentrating on the detection of range expansion events.     

Nested cladistic analysis indicates population fragmentation shapes genetic diversity in a freshwater mussel

Recently developed phylogeographic analyses that incorporate genealogical relationships of alleles offer the exciting prospect of disentangling historical from contemporary events. However, the relative advantages and shortfalls of this approach remain to be studied. We compared the nested cladistic method to the more traditional analysis of variance approach in a study of intraspecific genetic variation in the freshwater mussel, Lampsilis hydiana. We surveyed 257 specimens for nucleotide sequence level variation in a fragment of the mitochondrial 16S rRNA gene. When compared side by side, nested cladistic analysis and analysis of molecular variance (AMOVA) identified fragmentation of Arkansas river populations from remaining populations to the southwest. Nested cladistic analysis identified a second, more recent separation of Ouachita and Upper Saline river populations that was not detected by AMOVA. Differences among analytical methods probably arise from treatment of spatial hierarchical information: hierarchical groups emerge via a parsimony criterion in nested cladistic analysis but must be specified a priori in AMOVA. Both methods identified significant genetic structure among localities within hierarchical groups. Results from AMOVA suggested little gene flow among local populations with an island model. However, inferences about process that gave rise to patterns at this level were not possible in nested cladistic analysis, because an ancestral (interior) haplotype was not observed for a key one-step clade in the parsimony network. Our results suggest that, under some circumstances, nested cladistic analysis has lower power than more traditional analysis of variance to infer processes at the local population level.     

Cladistic structure within the human Lipoprotein lipase gene and its implications for phenotypic association studies

Haplotype variation in 9.7 kb of genomic DNA sequence from the human lipoprotein lipase (LPL) gene was scored in three populations: African-Americans from Jackson, Mississippi (24 individuals), Finns from North Karelia, Finland (24), and non-Hispanic whites from Rochester, Minnesota (23). Earlier analyses had indicated that recombination was common but concentrated into a hotspot and that recurrent mutations at multiple sites may have occurred. We show that much evolutionary structure exists in the haplotype variation on either side of the recombinational hotspot. By peeling off significant recombination events from a tree estimated under the null hypothesis of no recombination, we also reveal some cladistic structure not disrupted by recombination during the time to coalescence of this variation. Additional cladistic structure is estimated to have emerged after recombination. Many apparent multiple mutational events at sites still remain after removing the effects of the detected recombination/gene conversion events. These apparent multiple events are found primarily at sites identified as highly mutable by previous studies, strengthening the conclusion that they are true multiple events. This analysis portrays the complexity of the interplay among many recombinational and mutational events that would be needed to explain the patterns of haplotype diversity in this gene. The cladistic structure in this region is used to identify four to six single-nucleotide polymorphisms (SNPs) that would provide disequilibrium coverage over much of this region. These sites may be useful in identifying phenotypic associations with variable sites in this gene. Evolutionary considerations also imply that the SNPs in the 3' region should have general utility in most human populations, but the 5' SNPs may be more population specific. Choosing SNPs at random would generally not provide adequate disequilibrium coverage of the sequenced region.     

POLYMORPHIC TAXA, MISSING VALUES AND CLADISTIC ANALYSIS

Abstract Missing values have been used in cladistic analyses when data are unavailable, inapplicable or sometimes when character states are variable within terminal taxa. The practice of scoring taxa as having "missing values" for polymorphic characters introduces errors into the calculation of cladogram lengths and consistency indices because some character change is hidden within terminals. Because these hidden character steps are not counted, the set of most parsimonious cladograms may differ from those that would be found if polymorphic taxa had been broken into monomorphic subunits. In some cases, the trees found when polymorphisms are scored as missing values may not include any of the most parsimonious trees found when the data are scored properly. Additionally, in some cases, polymorphic taxa may be found to be polyphyletic when broken into monomorphic subunits; this is undetected when polymorphisms are treated as missing. Because of these problems, terminal units in cladistic analysis should be based on unique, fixed combinations of characters. Polymorphic taxa should be subdivided into subunits that are monomorphic for each character used in the analysis. Disregarding errors in topology, the additional hidden steps in a cladogram in which polymorphisms are scored as missing can be calculated by a simple formula, based on the observation that if it is assumed that polymorphic terminals include all combinations of character states, 2 ^p − 1 additional steps are required for each taxon in which p polymorphic binary characters are scored as missing values. Thus, when several polymorphisms are scored as missing in the same taxon, very large errors can be introduced into the calculation of tree length.     

Localization of a small genomic region associated with elevated ACE

Defining the relationship between multiple polymorphisms in a small genomic region and an underlying quantitative trait locus (QTL) represents a major challenge in human genetics. Pedigree analyses have shown that angiotensin I-converting enzyme (ACE) levels are influenced by a QTL located within or close to the ACE gene and most likely resides in the 3' region of this locus. We genotyped seven polymorphisms spanning 13 kb in the 3' end of ACE in 159 Afro-Caribbean subjects to evaluate the linkage disequilibrium between these sites and to narrow the genomic region associated with an elevated ACE level using a cladistic analysis. The linkage disequilibrium measurement D' and a haplotype tree revealed three distinct haplotype segments, presumably because of recombination. The value of the linkage disequilibrium parameter p(excess) was highest for site 22982, which is located in the middle segment. A series of nested, cladistic analyses confirmed that the other two regions are unlikely to be the ACE-linked QTL and that the variant resides in the middle region. Analyses of the same polymorphisms in 98 unrelated Europeans in the Monitoring Trends and Determinants in Cardiovascular Diseases (MONICA) study resulted in fewer haplotypes than were observed among the Afro-Caribbean subjects, suggesting that populations with greater genetic diversity may be especially informative for fine-scale mapping.     

Cladistic analysis of the alderflies of America north of Mexico (Megaloptera: Sialidae)

Abstract. A cladistic analysis of thirty-five morphological characters is used to produce a cladogram for the twenty-four American species of Sialis found north of Mexico, two European species, and three species of Protosialis . Three species groups ( americana, califomica and infumata ) proposed by Ross 1937 are shown to be monophyletic. Sialis, Protosialis , and the aequalis species group as currently constituted, are paraphyletic; S.nevadensis is nested within Protosialis and Protosialis is nested within Sialis . Character evolution is evaluated in the light of a phylogenetic hypothesis, and structural modifications are found to correlate with presumed functional constraints.