Haplotype Clusters Research Articles

Exceptional cryptic diversity and multiple origins of parthenogenesis in a freshwater ostracod

The persistence of asexual reproduction in many taxa depends on a balance between the origin of new asexual lineages and the extinction of old ones. This turnover determines the diversity of extant asexual populations and so influences the interaction between sexual and asexual modes of reproduction. Species with mixed reproduction, like the freshwater ostracod (Crustacea) morphospecies Eucypris virens, are a good model to examine these dynamics. This species is also a geographic parthenogen, in which sexual females and males co-exist with asexual females in the circum-Mediterranean area only, whereas asexual females occur all over Europe. A molecular phylogeny of E. virens based on the mitochondrial COI and 16S fragments is presented. It is characterised by many distinct clusters of haplotypes which are either exclusively sexual or asexual, with only one exception, and are often separated by deep branches. Analysis of the phylogeny reveals an astonishing cryptic diversity, which indicates the existence of a species complex with more than 40 cryptic taxa. We therefore suggest a revision of the single species status of E. virens. The phylogeny indicates multiple transitions from diverse sexual ancestor populations to asexuality. Although many transitions appear to be ancient, we argue that this may be an artefact of the existence of unsampled or extinct sexual lineages.

A novel method for haplotype clustering and visualization

In population genetics, it is common to represent the haplotype diversity at a genomic region between multiple populations using well-constructed visual representations. This typically requires the chromosomes from these populations to be grouped according to some definition of haplotypic similarity. Here, we introduce a novel algorithm for clustering haplotypes with the specific aim of addressing haplotype diversity within or between populations. The algorithm allows for missing data in the haplotypes and appropriately downweighs single nucleotide polymorphisms with higher extent of missingness. By identifying the canonical haplotypes in a genomic region, defined as the haplotype forms, which most chromosomes are similar to, the algorithm maps each chromosome to either a unique canonical haplotype or as a mosaic of the identified canonical haplotypes. This mapping can subsequently be utilized for producing graphical visualizations of the haplotype clustering for understanding the extent of haplotype diversity in the region. The clustering application has been implemented in R for distribution as haplosim, and we also provide a visualization script hapvisual for graphical display of the clustering results. The outcome of such analysis can be informative in understanding the extent of haplotype diversity between populations, in addressing the reproducibility of established association signals across multiple populations, and also in the investigation of positive selection in the human genome.

HOST SPECIALIZATION DIFFERENTIATES CRYPTIC SPECIES OF FEATHER-FEEDING LICE

Parasite species with differentiated host-specific populations provide a natural opportunity to explore factors involved in parasite diversification. Columbicola macrourae is a species of ectoparasitic feather louse currently recognized from 15 species of New World pigeons and doves. Mitochondrial sequences reveal five divergent haplotype clusters within C. macrourae, suggesting cryptic species. Each cluster is relatively host specific, with only one or a few hosts. We conducted a reciprocal transfer experiment with two of these lineages to test whether host use has an adaptive component. Our results demonstrate that the fitness of each lineage is considerably higher on its native host than on the novel host suggesting that one or more selective agents favor host specialization by the different lineages. In addition, we were able to morphologically separate individual lice from the two experimental lineages using discriminant function analysis. Furthermore, differences in the size of these louse lineages match differences in the size of their respective hosts, paralleling the strong correlation between parasite and host body size across the genus Columbicola. Together, these results suggest that selection in this cryptic species complex reflects selection across the whole genus, and that this selection, in part, contributes to the maintenance of host specialization.

Genetic structure of Schrenck newt Salamandrella schrenckii populations by mitochondrial cytochrome b variation

The nucleotide sequence variation of the mitochondrial cytochrome b gene was studied in Schrenck newt Salamandrella schrenckii (Strauch, 1870) from populations of Primorye and the Khabarovsk region. Phylogenetic analysis revealed two haplotype clusters, southern cluster 1 and northern cluster 2, with a divergence of 3%. Analysis of the mtDNA and cytochrome b amino acid sequence variations made it possible to assume that the modern range of Schrenck newt was colonized from south Primorye northwards. In contrast to the southern cluster, the northern one demonstrated all the signs of demographic expansion (a unimodal distribution of pairwise nucleotide differences, specific results of tests for selective neutrality of mtDNA variation, and a good correspondence of genetic parameters to those expected from demographic expansion models).

Effects of the number of markers per haplotype and clustering of haplotypes on the accuracy of QTL mapping and prediction of genomic breeding values

The aim of this paper was to compare the effect of haplotype definition on the precision of QTL-mapping and on the accuracy of predicted genomic breeding values. In a multiple QTL model using identity-by-descent (IBD) probabilities between haplotypes, various haplotype definitions were tested i.e. including 2, 6, 12 or 20 marker alleles and clustering base haplotypes related with an IBD probability of > 0.55, 0.75 or 0.95. Simulated data contained 1100 animals with known genotypes and phenotypes and 1000 animals with known genotypes and unknown phenotypes. Genomes comprising 3 Morgan were simulated and contained 74 polymorphic QTL and 383 polymorphic SNP markers with an average r2 value of 0.14 between adjacent markers. The total number of haplotypes decreased up to 50% when the window size was increased from two to 20 markers and decreased by at least 50% when haplotypes related with an IBD probability of > 0.55 instead of > 0.95 were clustered. An intermediate window size led to more precise QTL mapping. Window size and clustering had a limited effect on the accuracy of predicted total breeding values, ranging from 0.79 to 0.81. Our conclusion is that different optimal window sizes should be used in QTL-mapping versus genome-wide breeding value prediction.

Towards a comprehensive barcode library for arctic life - Ephemeroptera, Plecoptera, and Trichoptera of Churchill, Manitoba, Canada

BackgroundThis study reports progress in assembling a DNA barcode reference library for Ephemeroptera, Plecoptera, and Trichoptera ("EPTs") from a Canadian subarctic site, which is the focus of a comprehensive biodiversity inventory using DNA barcoding. These three groups of aquatic insects exhibit a moderate level of species diversity, making them ideal for testing the feasibility of DNA barcoding for routine biotic surveys. We explore the correlation between the morphological species delineations, DNA barcode-based haplotype clusters delimited by a sequence threshold (2%), and a threshold-free approach to biodiversity quantification--phylogenetic diversity.ResultsA DNA barcode reference library is built for 112 EPT species for the focal region, consisting of 2277 COI sequences. Close correspondence was found between EPT morphospecies and haplotype clusters as designated using a standard threshold value. Similarly, the shapes of taxon accumulation curves based upon haplotype clusters were very similar to those generated using phylogenetic diversity accumulation curves, but were much more computationally efficient.ConclusionThe results of this study will facilitate other lines of research on northern EPTs and also bode well for rapidly conducting initial biodiversity assessments in unknown EPT faunas.

Analysis of RYR1 Haplotype Profile in Patients with Malignant Hyperthermia

This study represents a new approach to characterising patients at risk of malignant hyperthermia (MH) through the use of a recently published method for identifying high-risk haplotypes in candidate genes. We present analysis based upon the largest standardised and genotyped database of MH patients worldwide. We used unphased RYR1 SNP data directly to (1) assess RYR1 haplotype frequency differences between susceptible cases and control groups and (2) analyse population-based association via clustering of RYR1 haplotypes based on disease risk. Our results show a significant difference in RYR1 haplotype frequency between susceptible cases and UK Caucasian population controls. Furthermore we identify a high-risk cluster of haplotypes that is associated with the commonest UK MH mutation p.G2434R/c.7300G>A. These results demonstrate the applicability of this new and practical method for population based association analysis.

Phylogeography of native ploidy levels and invasive tetraploids of Solidago gigantea

Ploidy level is an important aspect of the genetic makeup of a plant, and can strongly influence ecological characteristics such as invasiveness. We used a phylogeographical approach to elucidate the history of polyploidization and colonization success of diploid and tetraploid Solidago gigantea Aiton (Asteraceae) within its native range in North America. We were also able to identify the probable source material of the haplotype lineages invasive in Europe and Asia, where only tetraploid plants occur. To do this, we sequenced 1275 bp of chloroplast intergenic spacer DNA in 268 individuals from 57 populations. In addition, we performed a crossing experiment, which supported the hypothesis that chloroplast inheritance in this species is maternal. The phylogeographical analysis showed a complex pattern of 20 haplotypes of diploid and tetraploid plants. In North America, we found significant differentiation among regions, private haplotypes, and isolation by distance. Ploidy levels were more differentiated in the northern regions than in the South. The haplotype network was shallow and included one tetraploid-only, star-shaped cluster of haplotypes that were particularly successful colonizers. Post-glacial migration of diploid S. gigantea occurred mainly northwards east of the Appalachian Mountains, and to a lesser degree also southward. Our data suggest that tetraploids have formed several times in North America. Haplotype number and diversity were lower in European populations than in the native range, and we found evidence that four haplotypes were introduced to Europe from two source areas, New England and the Southern Appalachian Mountains.

Search for nucleotide diversity patterns of local adaptation in dehydrins and other cold-related candidate genes in Scots pine (Pinus sylvestris L.)

Nucleotide variation at several cold candidate genes including seven members of the dehydrin gene family was surveyed in haplotypes of Scots pine (Pinus sylvestris) sampled in populations showing divergence for cold tolerance in Europe. Patterns of nucleotide diversity, linkage disequilibrium, and frequency spectrum of alleles were compared between north and south populations to search for signs of directional selection potentially underlying adaptation to cold. Significant differentiation between populations in allelic frequency or haplotype structure was detected at dhn1, dhn3, and abaH loci. Allelic dimorphism with no evidence of haplotype clustering by geographical distribution was found at dhn9. An excess of fixed non-synonymous mutations as compared to the outgroup P. pinaster pine species was found at dhn1. Differences in nucleotide polymorphisms were found between the members of the Kn class of dehydrin upregulated during cold acclimation (average πsil = 0.004) as compared to the SKn class (average πsil = 0.024). The multilocus nucleotide diversity at silent sites (θW = 0.009) was moderate compared to other conifer species, but higher than previous estimates for Scots pine. There was an excess of rare and high frequency derived variants as revealed by significantly negative multilocus value of Tajima’s D (D = −0.72, P < 0.01) and negative mean value of Fay and Wu H statistics (H = −0.50). The level of linkage disequilibrium decayed rapidly with an average expected r2 of 0.2 at about 200 bp. Overall, there was a positive correlation between polymorphism and divergence at ten loci when outgroup sequence was available. The discovered polymorphism will be used for further evaluation of the adaptive role of genes through association mapping studies.

PXR Pharmacogenetics: Association of Haplotypes with Hepatic CYP3A4 and ABCB1 Messenger RNA Expression and Doxorubicin Clearance in Asian Breast Cancer Patients

To characterize pregnane X receptor (PXR) polymorphic variants in healthy Asian populations [Chinese, Malay and Indian (n=100 each)], and to investigate the association between PXR haplotypes and hepatic mRNA expression of PXR and its downstream target genes, CYP3A4 and ABCB1, as well as their influence on the clearance of doxorubicin in Asian breast cancer patients. PXR genotyping was done by direct DNA sequencing, and PXR haplotypes and haplotype clusters were derived by expectation-maximization algorithm. Genotype-phenotype correlations were done using Mann-Whitney U test and Kruskal-Wallis test. Significant interethnic variations were observed in PXR pharmacogenetics among the three Asian ethnic groups. The expression of PXR mRNA in liver tissues harboring the PXR*1B haplotype clusters was 4-fold lower compared with the non-PXR*1B (*1A + *1C) haplotype clusters [PXR*1B versus PXR*1A; P=0.015; PXR*1B versus PXR*1C; P=0.023]. PXR*1B-bearing liver tissues were associated with significantly lower expression of CYP3A4 (PXR*1B versus non-PXR*1B, P=0.030) and ABCB1 (PXR*1B versus non-PXR*1B, P=0.060) compared with non-PXR*1B-bearing liver tissues. Doxorubicin clearance in breast cancer patients harboring the PXR*1B haplotypes was significantly lower compared with patients carrying the non-PXR*1B haplotypes [PXR*1B versus non-PXR*1B, CL/BSA (L h(-1) m(-2)): 20.84 (range, 8.68-29.24) versus 24.85 (range, 13.80-55.66), P=0.022]. This study showed that PXR*1B was associated with reduced hepatic mRNA expression of PXR and its downstream targets, CYP3A4 and ABCB1. Genotype-phenotype correlates in breast cancer patients showed PXR*1B to be significantly associated with lower doxorubicin clearance, suggesting that PXR haplotype constitution could be important in influencing interindividual and interethnic variations in disposition of its putative drug substrates.

Statistical performance of cladistic strategies for haplotype grouping in pharmacogenetics

Haplotypes comprising multiple single nucleotide polymorphisms (SNPs) are popular covariates for capturing the key genetic variation present over a region of interest in the DNA sequence. Although haplotypes can provide a clearer assessment of genetic variation in a region than their component SNPs considered individually, the multi-allelic nature of haplotypes increases the complexity of the statistical models intended to discover association with outcomes of interest. Cladistic methods cluster haplotypes according to the estimates of their genealogical closeness and have been proposed recently as strategies for reducing model complexity and increasing power. Two examples are methods based on a haplotype nesting algorithm described by Templeton et al. (Genetics 1987; 117:343-351) and hierarchical clustering of haplotypes as described by Durrant et al. (Am. J. Hum. Genet. 2004; 75:35-43). In the context of assessing the pharmacogenetic effects of candidate genes, for which high-density SNP data have been gathered, we have conducted a simulation-based case study of the testing and estimation properties of two strategies based on Templeton's algorithm (TA), one being that described by Seltman et al. (Am. J. Hum. Genet. 2001; 68:1250-1263; Genet. Epidemiol. 2003; 25:48-58), as well as the method of Durrant et al. using data from a diabetes clinical trial. Even after adjusting for multiplicity, improvements in power can be realized using cladistic approaches with treatment group sizes in the range expected for standard trials, although these gains may be sensitive to the cladistic structure used. Differences in the relative performance of the cladistic approaches examined were observed with the clustering approach of Durrant et al. showing statistical properties superior to the methods based on TA.

CRH Haplotype as a Factor Influencing Cerebrospinal Fluid Levels of Corticotropin-Releasing Hormone, Hypothalamic-Pituitary-Adrenal Axis Activity, Temperament, and Alcohol Consumption in Rhesus Macaques

Both highly stress-reactive and novelty-seeking individuals are susceptible to alcohol use disorders. Variation in stress reactivity, exploration, and response to novelty have been attributed to differences in corticotropin-releasing hormone (CRH) system function. As such, CRH gene variation may influence risk for alcohol use and dependence. To determine whether CRH variation influences relevant intermediate phenotypes, behavior, and alcohol consumption in rhesus macaques. We sequenced the rhesus macaque CRH locus (rhCRH) and performed cladistic clustering of haplotypes. In silico analysis, gel shift, and in vitro reporter assays were performed to identify functional variants. Cerebrospinal fluid (CSF) and blood samples were obtained, and levels of CRH and corticotropin (ACTH) were measured by radioimmunoassay. Behavioral data were collected from macaques during infancy. Among adolescent/adult animals, we recorded responses to an unfamiliar conspecific and measured levels of ethanol consumption. National Institutes of Health Animal Center. Rhesus macaques. Animals were genotyped for a single-nucleotide polymorphism disrupting a glucocorticoid response element, rhCRH -2232 C>G, and the effects of this allele on CSF levels of CRH, plasma levels of ACTH, behavior, and ethanol consumption were assessed by analysis of variance. We show that -2232C>G alters DNA x protein interactions and confers decreased sensitivity of the CRH promoter to glucocorticoids in vitro. Consistent with the known effects of glucocorticoids on CRH expression in the brain, carriers of the G allele had lower CSF levels of CRH but higher levels of ACTH. Infants carrying the G allele were more exploratory and bold, and among adolescent and adult male macaques, the G allele was associated with exploratory/bold responding to an unfamiliar male. Adults with the C/G genotype also exhibited increased alcohol consumption in the social group, a model for high-risk alcohol-seeking behavior. Haplotypes that differ in terms of corticosteroid sensitivity have been identified in humans. Our data may suggest that functionally similar CRH variants could influence risk for externalizing disorders in human subjects.

Association mapping by generalized linear regression with density‐based haplotype clustering

Haplotypes of closely linked single-nucleotide polymorphisms (SNPs) potentially offer greater power than individual SNPs to detect association between genetic variants and disease. We present a novel approach for association mapping in which density-based clustering of haplotypes reduces the dimensionality of the general linear model (GLM)-based score test of association implemented in the HaploStats software (Schaid et al. [2002] Am. J. Hum. Genet. 70:425-434). A flexible haplotype similarity score, a generalization of previously used measures, forms the basis, for grouping haplotypes of probable recent common ancestry. All haplotypes within a cluster are assigned the same regression coefficient within the GLM, and evidence for association is assessed with a score statistic. The approach is applicable to both binary and continuous trait data, and does not require prior phase information. Results of simulation studies demonstrated that clustering enhanced the power of the score test to detect association, under a variety of conditions, while preserving valid Type-I error. Improvement in performance was most dramatic in the presence of extreme haplotype diversity, while a slight improvement was observed even at low diversity. Our method also offers, for binary traits, a slight advantage in power over a similar approach based on an evolutionary model (Tzeng et al. [2006] Am. J. Hum. Genet. 78:231-242).

Population structure and phylogeography of Aphyocypris kikuchii (Oshima) based on mitochondrial DNA variation

Aphyocypris kikuchii is a cyprinid species endemic to northern and eastern Taiwan and is the only primary freshwater fish native east of the Coastal Mountain Range. In total, 92 individuals of A. kikuchii from seven populations in three regions of the island were surveyed for mitochondrial DNA (mtDNA) variation. High haplotype diversity (h= 0·989) and low nucleotide diversity (π= 0·009) of mtDNA were detected. Negative values of Tajima’s D and unimodal mismatch distributions probably reflect a history of recent demographic expansions from small populations. Three major haplotype clusters displayed geographically non‐overlapping distributions, indicating a long‐term isolation between regions. Hierarchical analysis of molecular variance showed significant genetic structuring among populations (ΦST= 0·66). Significant haplotype heterogeneity was also detected among populations within regions (ΦSC= 0·41, P &lt; 0·001) and among regions (ΦCT= 0·43, P &lt; 0·05). Molecular clock estimates of coalescence in the three major mtDNA lineages indicated coalescence in the most recent common ancestor c. 0·11–0·39 million years ago. Haplotypes of cluster B nested as interior nodes in the haplotype network, indicating that migrations from Shueilian (SL) populations to the northern region (cluster A) and to the eastern region (cluster C) may have occurred independently. Lineages A and B + C should be managed as two distinct evolutionarily significant units, while the northern, SL and southern groups should be managed as separate management units.

Detection of QTL on Bovine X Chromosome by Exploiting Linkage Disequilibrium

A fine-mapping method exploiting linkage disequilibrium was used to detect quantitative trait loci (QTL) on the X chromosome affecting milk production, body conformation and productivity traits. The pedigree comprised 22 paternal half-sib families of Black-and-White Holstein bulls in the Netherlands in a grand-daughter design for a total of 955 sons. Twenty-five microsatellite markers were genotyped to construct a linkage map on the chromosome X spanning 170 Haldane cM with an average inter-marker distance of 7.1 cM. A covariance matrix including elements about identical-by-descent probabilities between haplotypes regarding QTL allele effects was incorporated into the animal model, and a restricted maximum-likelihood method was applied for the presence of QTL using the LDVCM program. Significance thresholds were obtained by permuting haplotypes to phenotypes and by using a false discovery rate procedure. Seven QTL responsible for conformation types (teat length, rump width, rear leg set, angularity and fore udder attachment), behavior (temperament) and a mixture of production and health (durable prestation) were detected at the suggestive level. Some QTL affecting teat length, rump width, durable prestation and rear leg set had small numbers of haplotype clusters, which may indicate good classification of alleles for causal genes or markers that are tightly associated with the causal mutation. However, higher maker density is required to better refine the QTL position and to better characterize functionally distinct haplotypes which will provide information to find causal genes for the traits.

CLUMPHAP: a simple tool for performing haplotype‐based association analysis

The completion of the HapMap Project and the development of high-throughput single nucleotide polymorphism genotyping technologies have greatly enhanced the prospects of identifying and characterizing the genetic variants that influence complex traits. In principle, association analysis of haplotypes rather than single nucleotide polymorphisms may better capture an underlying causal variant, but the multiple haplotypes can lead to reduced statistical power due to the testing of (and need to correct for) a large number of haplotypes. This paper presents a novel method based on clustering similar haplotypes to address this issue. The method, implemented in the CLUMPHAP program, is an extension of the CLUMP program designed for the analysis of multi-allelic markers (Sham and Curtis [1995] Ann. Hum. Genet. 59(Pt1):97-105). CLUMPHAP performs a hierarchical clustering of the haplotypes and then computes the chi(2) statistic between each haplotype cluster and disease; the statistical significance of the largest of the chi(2) statistics is obtained by permutation testing. A significant result suggests that the presence of a disease-causing variant in the haplotype cluster is over-represented in cases. Using simulation studies, we have compared CLUMPHAP and more widely used approaches in terms of their statistical power to identify an untyped susceptibility locus. Our results show that CLUMPHAP tends to have greater power than the omnibus haplotype test and is comparable in power to multiple regression locus-coding approaches.

Fine Mapping of the Psoriasis Susceptibility Locus PSORS1 Supports HLA-C as the Susceptibility Gene in the Han Chinese Population

PSORS1 (psoriasis susceptibility gene 1) is a major susceptibility locus for psoriasis. Several fine-mapping studies have highlighted a 300-kb candidate region of PSORS1 where multiple biologically plausible candidate genes were suggested. The most recent study has indicated HLA-Cw6 as the primary PSORS1 risk allele within the candidate region in a Caucasian population. In this study, a family-based association analysis of the PSORS1 locus was performed by analyzing 10 polymorphic microsatellite markers from the PSORS1 region as well as HLA-B, HLA-C and CDSN loci in 163 Chinese families of psoriasis. Five marker loci show strong evidence (P<10−3), and one marker locus shows weak evidence (P = 0.04) for association. The haplotype cluster analysis showed that all the risk haplotypes are Cw6 positive and share a 369-kb region of homologous marker alleles which carries all the risk alleles, including HLA-Cw6 and CDSN*TTC, identified in this study. The recombinant haplotype analysis of the HLA-Cw6 and CDSN*TTC alleles in 228 Chinese families showed that the HLA-Cw6 −/CDSN*TTC + recombinant haplotype is clearly not associated with risk for psoriasis (T∶NT = 29:57, p = 0.0025) in a Chinese population, suggesting that the CDSN*TTC allele itself does not confer risk without the presence of the HLA-Cw6 allele. The further exclusion analysis of the non-risk HLA-Cw6 −/CDSN*TTC + recombinant haplotypes with common recombination breakpoints has allowed us to refine the location of PSORS1 to a small candidate region. Finally, we performed a conditional linkage analysis and showed that the HLA-Cw6 is a major risk allele but does not explain the full linkage evidence of the PSORS1 locus in a Chinese population. By performing a series of family-based association analyses of haplotypes as well as an exclusion analysis of recombinant haplotypes, we were able to refine the PSORS1 gene to a small critical region where HLA-C is a strong candidate to be the PSORS1 susceptibility gene.

Patterns of nucleotide variation and reproductive isolation between aMimulusallotetraploid and its progenitor species

Here we report our characterization of a widespread, highly selfing Mimulus allotetraploid formed by interspecific hybridization between M. nasutus and M. guttatus. Nucleotide variation at two nuclear loci (mCYCA and mAP3) within and among tetraploid populations resolves two haplotype clusters for each locus: one shares near identity with sequences from M. nasutus and the other group shares substantial variation with M. guttatus. With respect to the two loci studied, each allotetraploid individual is a 'fixed heterozygote' carrying sequences from both clusters. Moreover, mCYCA variation is consistent with at least two evolutionary origins for the Mimulus allotetraploid. We show that the allotetraploid is strongly reproductively isolated from M. nasutus and M. guttatus; interploidy crosses produce almost no viable seeds. By extension, we infer strong triploid block and argue that Mimulus allotetraploid formation might proceed in one step via the union of unreduced gametes in an M. nasutus-M. guttatus F(1) hybrid. We also discuss the potential roles of mating system and flowering asynchrony in allotetraploid establishment.

Reticulate evolution in the natural range of the invasive wetland tree species Melaleuca quinquenervia

The Melaleuca leucadendra complex (broad-leaf paperbarks; Myrtaceae) is a dominant component of the tropical and sub-tropical biota of Australia, particularly in wetlands of high conservation significance. In Florida and other parts of the Americas, however, one member of the group (Melaleuca quinquenervia) is a serious ecological and economic weed. Understanding the relationships and evolution of the group is integral to both conservation and biocontrol efforts. Although the complex is currently considered to include up to 14 species, there has been some concern over taxonomic boundaries within the complex because most species are circumscribed only by combinations of characters, each of which also occurs in other species. Here, DNA sequence data derived from the chloroplast and two nuclear regions are used to explore the relationships of M. quinquenervia. We find little evidence for clear species boundaries within the M. leucadendra complex in general, with regional sharing of chloroplast haplotypes across morphologically defined taxa, indicating asymmetrical introgression or retention of ancestral haplotypes (lineage sorting). Phylogenies were further confounded by the recovery of multiple copies of both nuclear regions sequenced (ITS and rpb2) from many individuals. There was no clear evidence of polyploidy or pseudogenes, but multiple duplications of rpb2 could not be ruled out. Parsimony networks of the nuclear ITS region show some clustering of haplotypes by morphospecies but there is also evidence of both hybridisation and recombination. Signals of introgression were also evident in rpb2, supporting an hypothesis of recent or ongoing gene flow between M. quinquenervia and other members of the M. leucadendra complex. Both relaxed and fixed molecular-clock dating estimate the introgression to have occurred sometime within the past seven million years (95% CI: 0.7–18). The New Caledonian population of M. quinquenervia appears to have been established by dispersal from Australia during this period. M. quinquenervia is found to have alleles closely related to multiple different morphotaxa within the M. leucadendra complex, suggesting considerable past introgression into this taxon from some other members of the M. leucadendra complex, and this has implications for biocontrol efforts. The M. leucadendra complex appears to reflect early to intermediate stages of speciation, possibly driven by different ecologies.

Disease association tests by inferring ancestral haplotypes using a hidden markov model

Most genome-wide association studies rely on single nucleotide polymorphism (SNP) analyses to identify causal loci. The increased stringency required for genome-wide analyses (with per-SNP significance threshold typically approximately 10(-7)) means that many real signals will be missed. Thus it is still highly relevant to develop methods with improved power at low type I error. Haplotype-based methods provide a promising approach; however, they suffer from statistical problems such as abundance of rare haplotypes and ambiguity in defining haplotype block boundaries. We have developed an ancestral haplotype clustering (AncesHC) association method which addresses many of these problems. It can be applied to biallelic or multiallelic markers typed in haploid, diploid or multiploid organisms, and also handles missing genotypes. Our model is free from the assumption of a rigid block structure but recognizes a block-like structure if it exists in the data. We employ a Hidden Markov Model (HMM) to cluster the haplotypes into groups of predicted common ancestral origin. We then test each cluster for association with disease by comparing the numbers of cases and controls with 0, 1 and 2 chromosomes in the cluster. We demonstrate the power of this approach by simulation of case-control status under a range of disease models for 1500 outcrossed mice originating from eight inbred lines. Our results suggest that AncesHC has substantially more power than single-SNP analyses to detect disease association, and is also more powerful than the cladistic haplotype clustering method CLADHC. The software can be downloaded from http://www.imperial.ac.uk/medicine/people/l.coin.