Testing allele homogeneity: the problem of nested hypotheses

Background

Coronary artery disease (CAD), and one of its intermediate risk factors, dyslipidemia, possess a demonstrable genetic component, although the genetic architecture is incompletely defined. We previously reported a linkage peak on chromosome 5q31-33 for early-onset CAD where the strength of evidence for linkage was increased in families with higher mean low density lipoprotein-cholesterol (LDL-C). Therefore, we sought to fine-map the peak using association mapping of LDL-C as an intermediate disease-related trait to further define the etiology of this linkage peak. The study populations consisted of 1908 individuals from the CATHGEN biorepository of patients undergoing cardiac catheterization; 254 families (N = 827 individuals) from the GENECARD familial study of early-onset CAD; and 162 aorta samples harvested from deceased donors. Linkage disequilibrium-tagged SNPs were selected with an average of one SNP per 20 kb for 126.6-160.2 MB (region of highest linkage) and less dense spacing (one SNP per 50 kb) for the flanking regions (117.7-126.6 and 160.2-167.5 MB) and genotyped on all samples using a custom Illumina array. Association analysis of each SNP with LDL-C was performed using multivariable linear regression (CATHGEN) and the quantitative trait transmission disequilibrium test (QTDT; GENECARD). SNPs associated with the intermediate quantitative trait, LDL-C, were then assessed for association with CAD (i.e., a qualitative phenotype) using linkage and association in the presence of linkage (APL; GENECARD) and logistic regression (CATHGEN and aortas).

Results

We identified four genes with SNPs that showed the strongest and most consistent associations with LDL-C and CAD: EBF1, PPP2R2B, SPOCK1, and PRELID2. The most significant results for association of SNPs with LDL-C were: EBF1, rs6865969, p = 0.01; PPP2R2B, rs2125443, p = 0.005; SPOCK1, rs17600115, p = 0.003; and PRELID2, rs10074645, p = 0.0002). The most significant results for CAD were EBF1, rs6865969, p = 0.007; PPP2R2B, rs7736604, p = 0.0003; SPOCK1, rs17170899, p = 0.004; and PRELID2, rs7713855, p = 0.003.

Conclusion

Using an intermediate disease-related quantitative trait of LDL-C we have identified four novel CAD genes, EBF1, PRELID2, SPOCK1, and PPP2R2B. These four genes should be further examined in future functional studies as candidate susceptibility loci for cardiovascular disease mediated through LDL-cholesterol pathways.     

Testing hybridization hypotheses and evaluating the evolutionary potential of hybrids in mangrove plant species

Natural hybridization is of marked importance from global to local biological diversity. In mangroves, species ranges overlap extensively with one another and species share a long overlap of flowering time. Although hybridization has been suggested, patterns of hybridization and the evolutionary potential of hybrids are not yet fully understood. This study provides molecular evidence for the parental origins and status of hybrids in the dominant mangrove genus Rhizophora based on comparisons of chloroplast and nuclear phylogenies and estimations of genetic relatedness and structure from inter‐simple sequence repeat (ISSR) markers. Phylogenetic analyses indicate that almost all species can act as maternal parents to hybrids and that hybridization can be bidirectional. Bayesian analyses indicate that hybrids are simple F₁s, and no trace of backcrossing was detected within populations. Hybridization, for the most part, occurs almost only locally and dispersal of hybrid individuals is limited beyond the hybrid sites.     

Testing hypotheses in order

In certain circumstances, one wishes to test one hypothesisonly if certain other hypotheses have been rejected. This orderingof hypotheses simplifies the task of controlling the probabilityof rejecting any true hypothesis. In an example from an observationalstudy, a treated group is shown to be further from both of twocontrol groups than the two control groups are from each other.     

Paleoecology and coalescence: phylogeographic analysis of hypotheses from the fossil record

The application of principles from coalescence theory to genealogical relationships within species can provide insights into the process of diversification and the influence of biogeography on distributional patterns. There are several features that make some organisms more suitable for detailed studies of historical processes; in particular, limited dispersal, which serves to conserve the patterns of genetic variation that developed during colonization. We describe the potential benefits of studies that integrate analyses of genetic variation with information from the fossil pollen record and present recent examples of the application of quantitative methods of phylogeographic analysis.     

Testing "species pair" hypotheses: evolutionary processes in the lichen-forming species complex Porpidia flavocoerulescens and Porpidia melinodes

Pairs of taxa are commonly found in lichen-forming ascomycetes that differ primarily in their reproductive modes: one taxon reproduces sexually, the other vegetatively. The evolutionary processes underlying such "species pairs" are unknown. The species pair formed by Porpidia flavocoerulescens (sexual) and Porpidia melinodes (vegetative) was chosen to investigate four previously proposed hypotheses. These hypotheses posit that species pairs are either two monophyletic, independently evolving species with contrasting reproductive mode; a single outcrossing species polymorphic with regard to its reproductive modes; a sexual mother lineage frequently giving rise to asexual spin-offs; or a complex of cryptic species. The phylogenetic patterns observed within the species pair in the present study were analyzed using stringent hypothesis testing and visualizations of relationships and conflict based on tree and network reconstructions. DNA sequences at the three analyzed loci revealed the same four to five deeply divergent lineages. A detailed analysis of DNA-sequence variability revealed closely linked gene loci, but high levels of conflict within each of the gene fragments, as well as between observed genetic lineages. The observed patterns of phylogenetic relationships, linkage, and conflict are not congruent with any of the previously proposed species pair hypotheses. Rather, it is proposed that the observed results can be explained by conflicting reproductive and nutritional requirements imposed by an obligate symbiotic lifestyle. These interacting constraints produce recurring selective sweeps within predominantly vegetatively reproducing lineages and are the main forces that shape the evolution within the investigated species pair.     

The causes of phylogenetic conflict in a classic Drosophila species group

Bifurcating phylogenies are frequently used to describe the evolutionary history of groups of related species. However, simple bifurcating models may poorly represent the evolutionary history of species that have been exchanging genes. Here, we show that the history of three well-known closely related species, Drosophila pseudoobscura, D. persimilis and D. p. bogotana, is not well represented by a bifurcating phylogenetic tree. The phylogenetic relationships among these species vary widely between different genomic regions. Much of this phylogenetic variation can be explained by the potential of different genomic regions to introgress between species, as measured in laboratory studies. We argue that the utility of multiple markers in species-level phylogenetic studies can be greatly enhanced by knowledge of genomic location and, in the case of hybridizing species, by knowledge of the functional or linkage relationships among the markers and regions of the genome that reduce hybrid fitness.     

Testing species boundaries in an ancient species complex with deep phylogeographic history: genus Xantusia (Squamata: Xantusiidae)

Identification of species in natural populations has recently received increased attention with a number of investigators proposing rigorous methods for species delimitation. Morphologically conservative species (or species complexes) with deep phylogenetic histories (and limited gene flow) are likely to pose particular problems when attempting to delimit species, yet this is crucial to comparative studies of the geography of speciation. We apply two methods of species delimitation to an ancient group of lizards (genus Xantusia) that occur throughout southwestern North America. Mitochondrial cytochrome b and nicotinamide adenine dinucleotide dehydrogenase subunit 4 gene sequences were generated from samples taken throughout the geographic range of Xantusia. Maximum likelihood, Bayesian, and nested cladogram analyses were used to estimate relationships among haplotypes and to infer evolutionary processes. We found multiple well-supported independent lineages within Xantusia, for which there is considerable discordance with the currently recognized taxonomy. High levels of sequence divergence (21.3%) suggest that the pattern in Xantusia may predate the vicariant events usually hypothesized for the fauna of the Baja California peninsula, and the existence of deeply divergent clades (18.8%-26.9%) elsewhere in the complex indicates the occurrence of ancient sundering events whose genetic signatures were not erased by the late Wisconsin vegetation changes. We present a revised taxonomic arrangement for this genus consistent with the distinct mtDNA lineages and discuss the phylogeographic history of this genus as a model system for studies of speciation in North American deserts.     

Integrating phylogenetic and population genetic analyses of multiple loci to test species divergence hypotheses in Passerina buntings

Phylogenetic and population genetic analyses of DNA sequence data from 10 nuclear loci were used to test species divergence hypotheses within Passerina buntings, with special focus on a strongly supported, but controversial, sister relationship between Passerina amoena and P. caerulea inferred from a previous mitochondrial study. Here, a maximum-likelihood analysis of a concatenated 10-locus data set, as well as minimize-deep-coalescences and maximum-likelihood analyses of the locus-specific gene trees, recovered the traditional sister relationship between P. amoena and P. cyanea. In addition, a more recent divergence time estimate between P. amoena and P. cyanea than between P. amoena and P. caerulea provided evidence for the traditional sister relationship. These results provide a compelling example of how lineage sorting stochasticity can lead to incongruence between gene trees and species trees, and illustrate how phylogenetic and population genetic analyses can be integrated to investigate evolutionary relationships between recently diverged taxa.     

The use of parsimony in testing phylogenetic hypotheses

With the advance of cladistic theory differences in principle between it and other systematic techniques are few but of fundamental importance. In the mechanics of classification they are confined to ranking and the rejection of paraphyletic taxa. In cladistic analysis, leading to cladograms, trees and phylogeny reconstruction, inconsistencies in apparent synapomorphies are said to be resolved using Popper's hypothetico-deductive method together with the principle of parsi However, not only do cladists not use Popper's methodology, which is inconsistent with parsimony, but their use of parsimony is invalid as a test of phylo The only accepted extrinsic test of a classification is that enunciated by John Stuart Mill. It has been claimed that cladistic classifications yield the best results when judged by Mill's criteria, but this is only possibly the case with analytic classifications produced by numerical techniques. No satisfactory test exists in normal (synthetic) cladism for distinguishing synapomorphy from homoplasy. The effects of this are particularly dire in cladograms and classifications involving fossils in which a Stufenreihe arrangement is adopted.     

Testing morphology-based hypotheses of phylogenetic relationships in Parmeliaceae (Ascomycota) using three ribosomal markers and the nuclear RPB1 gene

Parmeliaceae is the largest family of lichen-forming fungi with more than 2000 species and includes taxa with different growth forms. Morphology was widely employed to distinguish groups within this large, cosmopolitan family. In this study we test these morphology-based groupings using DNA sequence data from three nuclear and one mitochondrial marker from 120 taxa that include 59 genera and represent the morphological and chemical diversity in this lineage. Parmeliaceae is strongly supported as monophyletic and six well-supported main clades can be distinguished within the family. The relationships among them remain unresolved. The clades largely agree with the morphology-based groupings and only the placement of four of the genera studied is rejected by molecular data, while four other genera belong to clades previously unrecognised. The classification of these previously misplaced genera, however, has already been questioned by some authors based on morphological evidence. These results support morphological characters as important for the identification of monophyletic clades within Parmeliaceae.     

Testing hypotheses about the sister group of the passeriformes using an independent 30-locus data set

Although many phylogenetic studies have focused on developing hypotheses about relationships, advances in data collection and computation have increased the feasibility of collecting large independent data sets to rigorously test controversial hypotheses or carefully assess artifacts that may be misleading. One such relationship in need of independent evaluation is the position of Passeriformes (perching birds) in avian phylogeny. This order comprises more than half of all extant birds, and it includes one of the most important avian model systems (the zebra finch). Recent large-scale studies using morphology, mitochondrial, and nuclear sequence data have generated very different hypotheses about the sister group of Passeriformes, and all conflict with an older hypothesis generated using DNA-DNA hybridization. We used novel data from 30 nuclear loci, primarily introns, for 28 taxa to evaluate five major a priori hypotheses regarding the phylogenetic position of Passeriformes. Although previous studies have suggested that nuclear introns are ideal for the resolution of ancient avian relationships, introns have also been criticized because of the potential for alignment ambiguities and the loss of signal due to saturation. To examine these issues, we generated multiple alignments using several alignment programs, varying alignment parameters, and using guide trees that reflected the different a priori hypotheses. Although different alignments and analyses yielded slightly different results, our analyses excluded all but one of the five a priori hypotheses. In many cases, the passerines were sister to the Psittaciformes (parrots), and taxa were members of a larger clade that includes Falconidae (falcons) and Cariamidae (seriemas). However, the position of Coliiformes (mousebirds) was highly unstable in our analyses of 30 loci, and this represented the primary source of incongruence among analyses. Mousebirds were united with passerines or parrots in some analyses, suggesting an additional hypothesis that needs to be considered in future studies. There was no clear evidence that base-compositional convergence, saturation, or long-branch attraction affected our conclusions. These results provide independent evidence excluding four major hypotheses about the position of passerines, allowing the extensive studies on this group to be placed in a more rigorous evolutionary framework.     

Deep and wide valleys drive nested phylogeographic patterns across a montane bird community

Montane species distributions interrupted by valleys can lead to range fragmentation, differentiation and ultimately speciation. Paleoclimatic fluctuations may accentuate or reduce such diversification by temporally altering the extent of montane habitat and may affect species differentially. We examined how an entire montane bird community of the Western Ghats—a linear, coastal tropical mountain range—responds to topographic valleys that host different habitats. Using genetic data from 23 species (356 individuals) collected across nine locations, we examined if different species in the community reveal spatial concordance in population differentiation, and whether the timing of these divergences correlate with climatic events. Our results reveal a nested effect of valleys, with several species (10 of 23) demonstrating the oldest divergences associated with the widest and deepest valley in the mountain range, the Palghat Gap. Further, a subset of these 10 species revealed younger divergences across shallower, narrower valleys. We recovered discordant divergence times for all valley-affected montane birds, mostly in the Pleistocene, supporting the Pliestocene-pump hypotheses and highlighting the role of climatic fluctuations during this period in driving species evolution. A majority of species remain unaffected by valleys, perhaps owing to geneflow or extinction–recolonization dynamics. Studying almost the entire community allowed us to uncover a range of species’ responses, including some generalizable and other unpredicted patterns.     

Testing hypotheses regarding the genetics of adaptation

Many of the hypotheses regarding the genetics of adaptation require that one know specific details about the genetic basis of complex traits, such as the number and effects of the loci involved. Developments in molecular biology have made it possible to create relatively dense maps of markers that can potentially be used to map genes underlying specific traits. However, there are a number of reasons to doubt that such mapping will provide the level of resolution necessary to specifically address many evolutionary questions. Moreover, evolutionary change is built upon the substitution of individual mutations, many of which may now be cosegregating in the same allele. In order for this developing area not to become a mirage that traps the efforts of an entire field, the genetic dissection of adaptive traits should be conducted within a strict hypothesis-testing framework and within systems that promise a reasonable chance of identifying the specific genetic changes of interest. Continuing advances in molecular technology may lead the way here, but some form of genetic testing is likely to be forever required.     

An ITS-based phylogenetic analysis of theEuryptera species group inLomatium (Apiaceae)

Lomatium, the largest genus of Apiaceae in western North America, includes many narrow endemics whose relationships are uncertain. Although no infrageneric classification exists forLomatium, several informal groups have been recognized. TheEuryptera group comprises seven narrowly endemic species distributed primarily in California. We conducted parsimony and maximum likelihood (ML) analyses using sequences of the internal transcribed spacers (ITS) of the nuclear ribosomal DNA from species of theEuryptera group and several other species ofLomatium. When considered with distribution, morphological, and cpDNA data, the ITS analyses are consistent with the monophyly of theEuryptera group and suggests that speciation in this group has occurred through geographical divergence. Inferences from ITS data also identify putative progenitors of the polyploidEuryptera species.     

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.     

Extensive gene flow blurs phylogeographic but not phylogenetic signal in Olea europaea L.

Genetic structure and evolutionary patterns of the wild olive tree (Olea europaea L.) were investigated with AFLP fingerprinting data at three geographic levels: (a) phylogenetic relationships of the six currently recognized subspecies in Eurasia and Africa; (b) lineage identification in subsp. europaea of the Mediterranean basin; and (c) phylogeography in the western Mediterranean. Two statistical approaches (Bayesian inference and analysis of molecular variance) were used to analyse the AFLP fingerprints. To determine the congruency and transferability of results across studies previous RAPD and ISSR data were analysed in a similar manner. Comparisons proved that qualitative results were mostly congruent but quantitative values differed, depending on the method of analysis. Neighbour-Joining analysis of AFLP phenotypes supported current classification of subspecies. At a Mediterranean scale no clear cut phylogeographic pattern was recovered, likely due to extensive gene flow between populations of subsp. europaea. Gene flow estimates calculated with conventional F-statistics showed that reproductive barriers separated neither populations nor lineages of O. europaea. Genetic divergence between eastern and western parts of the Mediterranean basin was observed only when geographical and population information were incorporated into the analyses through hierarchical analysis of molecular variance (AMOVA). Within the western Mediterranean, the highest genetic diversity was found in two regions: on both sides of the Strait of Gibraltar and in the Balearic archipelago. Additionally, long-lasting isolation of the northern-most populations of the Iberian Peninsula appeared to be responsible for a significant divergence.Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .