The Framingham Heart Study 100K SNP genome-wide association study resource: overview of 17 phenotype working group reports

BACKGROUND:The Framingham Heart Study (FHS), founded in 1948 to examine the epidemiology of cardiovascular disease, is among the most comprehensively characterized multi-generational studies in the world. Many collected phenotypes have substantial genetic contributors; yet most genetic determinants remain to be identified. Using single nucleotide polymorphisms (SNPs) from a 100K genome-wide scan, we examine the associations of common polymorphisms with phenotypic variation in this community-based cohort and provide a full-disclosure, web-based resource of results for future replication studies.METHODS:Adult participants (n = 1345) of the largest 310 pedigrees in the FHS, many biologically related, were genotyped with the 100K Affymetrix GeneChip. These genotypes were used to assess their contribution to 987 phenotypes collected in FHS over 56 years of follow up, including: cardiovascular risk factors and biomarkers; subclinical and clinical cardiovascular disease; cancer and longevity traits; and traits in pulmonary, sleep, neurology, renal, and bone domains. We conducted genome-wide variance components linkage and population-based and family-based association tests.RESULTS:The participants were white of European descent and from the FHS Original and Offspring Cohorts (examination 1 Offspring mean age 32 +/- 9 years, 54% women). This overview summarizes the methods, selected findings and limitations of the results presented in the accompanying series of 17 manuscripts. The presented association results are based on 70,897 autosomal SNPs meeting the following criteria: minor allele frequency [greater than or equal to] 10%, genotype call rate [greater than or equal to] 80%, Hardy-Weinberg equilibrium p-value [greater than or equal to] 0.001, and satisfying Mendelian consistency. Linkage analyses are based on 11,200 SNPs and short-tandem repeats. Results of phenotype-genotype linkages and associations for all autosomal SNPs are posted on the NCBI dbGaP website at http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?id=phs000007.CONCLUSION:We have created a full-disclosure resource of results, posted on the dbGaP website, from a genome-wide association study in the FHS. Because we used three analytical approaches to examine the association and linkage of 987 phenotypes with thousands of SNPs, our results must be considered hypothesis-generating and need to be replicated. Results from the FHS 100K project with NCBI web posting provides a resource for investigators to identify high priority findings for replication.

Genome-wide association with bone mass and geometry in the Framingham Heart Study

BACKGROUND:Osteoporosis is characterized by low bone mass and compromised bone structure, heritable traits that contribute to fracture risk. There have been no genome-wide association and linkage studies for these traits using high-density genotyping platforms.METHODS:We used the Affymetrix 100K SNP GeneChip marker set in the Framingham Heart Study (FHS) to examine genetic associations with ten primary quantitative traits: bone mineral density (BMD), calcaneal ultrasound, and geometric indices of the hip. To test associations with multivariable-adjusted residual trait values, we used additive generalized estimating equation (GEE) and family-based association tests (FBAT) models within each sex as well as sexes combined. We evaluated 70,987 autosomal SNPs with genotypic call rates [greater than or equal to]80%, HWE p [greater than or equal to] 0.001, and MAF [greater than or equal to]10% in up to 1141 phenotyped individuals (495 men and 646 women, mean age 62.5 yrs). Variance component linkage analysis was performed using 11,200 markers.RESULTS:Heritability estimates for all bone phenotypes were 30-66%. LOD scores [greater than or equal to]3.0 were found on chromosomes 15 (1.5 LOD confidence interval: 51,336,679-58,934,236 bp) and 22 (35,890,398-48,603,847 bp) for femoral shaft section modulus. The ten primary phenotypes had 12 associations with 100K SNPs in GEE models at p < 0.000001 and 2 associations in FBAT models at p < 0.000001. The 25 most significant p-values for GEE and FBAT were all less than 3.5 x 10-6 and 2.5 x 10-5, respectively. Of the 40 top SNPs with the greatest numbers of significantly associated BMD traits (including femoral neck, trochanter, and lumbar spine), one half to two-thirds were in or near genes that have not previously been studied for osteoporosis. Notably, pleiotropic associations between BMD and bone geometric traits were uncommon. Evidence for association (FBAT or GEE p < 0.05) was observed for several SNPs in candidate genes for osteoporosis, such as rs1801133 in MTHFR; rs1884052 and rs3778099 in ESR1; rs4988300 in LRP5; rs2189480 in VDR; rs2075555 in COLIA1; rs10519297 and rs2008691 in CYP19, as well as SNPs in PPARG (rs10510418 and rs2938392) and ANKH (rs2454873 and rs379016). All GEE, FBAT and linkage results are provided as an open-access results resource at http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?id=phs000007.CONCLUSION:The FHS 100K SNP project offers an unbiased genome-wide strategy to identify new candidate loci and to replicate previously suggested candidate genes for osteoporosis.

A genome-wide association study for blood lipid phenotypes in the Framingham Heart Study

BACKGROUND:Blood lipid levels including low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and triglycerides (TG) are highly heritable. Genome-wide association is a promising approach to map genetic loci related to these heritable phenotypes.METHODS:In 1087 Framingham Heart Study Offspring cohort participants (mean age 47 years, 52% women), we conducted genome-wide analyses (Affymetrix 100K GeneChip) for fasting blood lipid traits. Total cholesterol, HDL-C, and TG were measured by standard enzymatic methods and LDL-C was calculated using the Friedewald formula. The long-term averages of up to seven measurements of LDL-C, HDL-C, and TG over a ~30 year span were the primary phenotypes. We used generalized estimating equations (GEE), family-based association tests (FBAT) and variance components linkage to investigate the relationships between SNPs (on autosomes, with minor allele frequency [greater than or equal to]10%, genotypic call rate [greater than or equal to]80%, and Hardy-Weinberg equilibrium p [greater than or equal to] 0.001) and multivariable-adjusted residuals. We pursued a three-stage replication strategy of the GEE association results with 287 SNPs (P < 0.001 in Stage I) tested in Stage II (n ~1450 individuals) and 40 SNPs (P < 0.001 in joint analysis of Stages I and II) tested in Stage III (n~6650 individuals).RESULTS:Long-term averages of LDL-C, HDL-C, and TG were highly heritable (h2 = 0.66, 0.69, 0.58, respectively; each P < 0.0001). Of 70,987 tests for each of the phenotypes, two SNPs had p < 10-5 in GEE results for LDL-C, four for HDL-C, and one for TG. For each multivariable-adjusted phenotype, the number of SNPs with association p < 10-4 ranged from 13 to 18 and with p < 10-3, from 94 to 149. Some results confirmed previously reported associations with candidate genes including variation in the lipoprotein lipase gene (LPL) and HDL-C and TG (rs7007797; P = 0.0005 for HDL-C and 0.002 for TG). The full set of GEE, FBAT and linkage results are posted at the database of Genotype and Phenotype (dbGaP). After three stages of replication, there was no convincing statistical evidence for association (i.e., combined P < 10-5 across all three stages) between any of the tested SNPs and lipid phenotypes.CONCLUSION:Using a 100K genome-wide scan, we have generated a set of putative associations for common sequence variants and lipid phenotypes. Validation of selected hypotheses in additional samples did not identify any new loci underlying variability in blood lipids. Lack of replication may be due to inadequate statistical power to detect modest quantitative trait locus effects (i.e., < 1% of trait variance explained) or reduced genomic coverage of the 100K array. GWAS in FHS using a denser genome-wide genotyping platform and a better-powered replication strategy may identify novel loci underlying blood lipids.

Genetic correlates of longevity and selected age-related phenotypes: a genome-wide association study in the Framingham Study

BACKGROUND: Family studies and heritability estimates provide evidence for a genetic contribution to variation in the human life span. METHODS:We conducted a genome wide association study (Affymetrix 100K SNP GeneChip) for longevity-related traits in a community-based sample. We report on 5 longevity and aging traits in up to 1345 Framingham Study participants from 330 families. Multivariable-adjusted residuals were computed using appropriate models (Cox proportional hazards, logistic, or linear regression) and the residuals from these models were used to test for association with qualifying SNPs (70, 987 autosomal SNPs with genotypic call rate [greater than or equal to]80%, minor allele frequency [greater than or equal to]10%, Hardy-Weinberg test p [greater than or equal to] 0.001).RESULTS:In family-based association test (FBAT) models, 8 SNPs in two regions approximately 500 kb apart on chromosome 1 (physical positions 73,091,610 and 73, 527,652) were associated with age at death (p-value < 10-5). The two sets of SNPs were in high linkage disequilibrium (minimum r2 = 0.58). The top 30 SNPs for generalized estimating equation (GEE) tests of association with age at death included rs10507486 (p = 0.0001) and rs4943794 (p = 0.0002), SNPs intronic to FOXO1A, a gene implicated in lifespan extension in animal models. FBAT models identified 7 SNPs and GEE models identified 9 SNPs associated with both age at death and morbidity-free survival at age 65 including rs2374983 near PON1. In the analysis of selected candidate genes, SNP associations (FBAT or GEE p-value < 0.01) were identified for age at death in or near the following genes: FOXO1A, GAPDH, KL, LEPR, PON1, PSEN1, SOD2, and WRN. Top ranked SNP associations in the GEE model for age at natural menopause included rs6910534 (p = 0.00003) near FOXO3a and rs3751591 (p = 0.00006) in CYP19A1. Results of all longevity phenotype-genotype associations for all autosomal SNPs are web posted at http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?id=phs000007. CONCLUSION: Longevity and aging traits are associated with SNPs on the Affymetrix 100K GeneChip. None of the associations achieved genome-wide significance. These data generate hypotheses and serve as a resource for replication as more genes and biologic pathways are proposed as contributing to longevity and healthy aging.

Rule Extraction, Fuzzy ARTMAP, and Medical Databases

This paper shows how knowledge, in the form of fuzzy rules, can be derived from a self-organizing supervised learning neural network called fuzzy ARTMAP. Rule extraction proceeds in two stages: pruning removes those recognition nodes whose confidence index falls below a selected threshold; and quantization of continuous learned weights allows the final system state to be translated into a usable set of rules. Simulations on a medical prediction problem, the Pima Indian Diabetes (PID) database, illustrate the method. In the simulations, pruned networks about 1/3 the size of the original actually show improved performance. Quantization yields comprehensible rules with only slight degradation in test set prediction performance.

A medical study on wireless inertial measurement technology as a tool for identifying patients at risk of death or imminent clinical deterioration

This paper provides a system description and preliminary results for an ongoing clinical study currently being carried out at the Mid-Western Regional Hospital, Nenagh, Ireland. The goal of the trial is to determine if wireless inertial measurement technology can be employed to identify elderly patients at risk of death or imminent clinical deterioration. The system measures cumulative movement and provides a score that will help provide a robust early warning to clinical staff of clinical deterioration. In addition the study examines some of the logistical barriers to the adoption of wearable wireless technology in front-line medical care.

Molecular genetics of the imprinted gene - SPROUTY3 (SPRY3)

Sprouty proteins are key regulators of cell growth and branching morphogenesis during development. Human SPRY3 which maps to the pseudoautosomal region 2, undergoes random X-inactivation in females and preferential Y-inactivation in males, behaving as though genetically X-linked. Spry3 is widely expressed in neuronal tissues, being found at high levels in the cerebellum and particularly in the Purkinje cells which, notably, are deficient in the autistic brain. Spry3 is also highly expressed in other ganglia in adults including retinal ganglion cells, dorsal root ganglion and superior cervical ganglion. SPRY3 enhancer can drive SPRY3 expression in the lung airway, which is consistent with a role in branching morphogenesis and the function of the original Drosophila Spry gene, which is critical for lung morphogenesis, providing a possible explanation for an observed anatomic abnormality in the autistic lung airway. In the human and mouse, the SPRY3 core promoter contains an AG-rich repeat and we found evidence of coexpression, promoter binding and regulation of SPRY3 expression by transcription factors EGR1, ZNF263 and PAX6. Spry3 over-expression in mouse superior cervical ganglion cells inhibits axon branching and Spry3 knockdown in those cells increases axon branching, consistent with known functions of other Sprouty proteins. Novel SPRY3 upstream transcripts that I characterised originate from three start sites in the X-linked F8A3 – TMLHE gene region, which is recently implicated in autism causation. Arising from these findings, I propose that the lung airway abnormality and low levels of blood carnitine found in autism suggest that deregulation of SPRY3 may underpin a subset of autism cases.

The evolution of the medical professions in eighteenth-century Ireland: An institutional perspective

Ireland, in the eighteenth century, followed the classic tripartite division of regular medical practitioners into physicians, surgeons and apothecaries. At the beginning of the century surgeons and apothecaries were regarded as mere tradesmen, but by the end of the century both were regarded as professionals and had the right to regulate their respective professions. Practitioners in different regions of Europe developed in a different manner, and eighteenth-century practitioners in Ireland developed independently from their English counterparts. In common with Britain and Europe in the eighteenth century, the total number of practitioners increased in Ireland, and by the end of the century, apothecaries were the largest group in Dublin, closely followed by the surgeons. Surgeons and apothecaries at the start of the eighteenth century belonged to the same guild. However in mid-century, St Luke's guild of apothecaries was established and this provided the apothecaries with a new identity that allowed them to pursue auto regulation, rather than hitherto, when they had been regulated by the physicians. This was vital to the apothecaries as they were in direct commercial competition with both the physicians and the surgeons and faced increasing pressure from both druggists and the disparate group of practitioners known as the irregulars. The 1765 County Infirmaries Act established a hospital in virtually every county in Ireland, and cast the surgeon as the primary medical officer in the countrywide network of hospitals. This legislation, which was unique in Europe, had the unintended consequence of elevating the status of the surgeons, as prior to this physicians were always in the ascendancy in the voluntary hospitals in Ireland and Britain, in contrast to France. The status of the surgeons was further enhanced by the establishment of the College of Surgeons in Ireland in 1784, which provided them with a new corporate identity, the authority to regulate the profession countrywide, and, also, the ability to educate surgeons in Ireland. The establishment of the College of Surgeons placed further pressure on the apothecaries to demonstrate that they also had a recognisable identity, and the authority to regulate their own profession. This was achieved with the 1791 Apothecaries Act which established the Apothecaries Hall and give the apothecaries the right to regulate themselves. This innovative legislation deemed the apothecaries a profession, and was enacted twenty-four years prior to similar legislation in Britain. Commercial pressure from druggists and, probably, irregulars expedited the requirement of the apothecaries to establish a new corporate identity, in order to distance themselves from these groups. The changing status of both apothecaries and surgeons had little effect on the physicians as a group, and, despite being the beneficiaries of a generous bequest from Sir Patrick Dun in 1711 to provide medical chairs in Dublin, the physicians displayed an inertia during the eighteenth century that was not in keeping with the developments that occurred in the contemporary Dublin medical world. The fact that it took ninety-five years, and that five acts of parliament, two House of Commons enquiries and a House of Lords enquiry were required to ensure that Dun's wishes were brought to fruition demonstrates that the physicians did not develop at the same pace as the other medical groups in the city. Had Dun’s bequest been implemented as he desired, Dublin, with a number of voluntary hospitals, would have been well placed to provide comprehensive tuition for medical students in the eighteenth century. It was not until the nineteenth century that the city, and the populace, benefited from this legacy. This thesis will trace these developments in the context of changes that occurred in contemporary medical education and diagnosis in Ireland, Britain and France. It will demonstrate that Irish practitioners developed independently, influenced mainly by local issues, but also by those who had travelled abroad and returned to Ireland with new concepts and ideas, ensuring that Irish medical practitioners had the institutional structure that could encompass the diagnostic and regulatory changes that would become accepted in the nineteenth century.

Researcher practices on returning genetic research results.

BACKGROUND/AIMS: as genetic and genomic research proliferates, debate has ensued about returning results to participants. In addition to consideration of the benefits and harms to participants, researchers must also consider the logistical and financial feasibility of returning research results. However, little data exist of actual researcher practices. METHODS: we conducted an online survey of 446 corresponding authors of genetic/genomic studies conducted in the United States and published in 2006-2007 to assess the frequency with which they considered, offered to, or actually returned research results, what factors influenced these decisions, and the method of communicating results. RESULTS: the response rate was 24% (105/446). Fifty-four percent of respondents considered the issue of returning research results to participants, 28% offered to return individual research results, and 24% actually returned individual research results. Of those who considered the issue of returning research results during the study planning phase, the most common factors considered were whether research results were deemed clinically useful (18%) and respect for participants (13%). Researchers who had a medical degree and conducted studies on children were significantly more likely to offer to return or actually return individual results compared to those with a Ph.D. only. CONCLUSIONS: we speculate that issues associated with clinical validity and respect for participants dominated concerns of time and expense given the prominent and continuing ethical debates surrounding genetics and genomics research. The substantial number of researchers who did not consider returning research results suggests that researchers and institutional review boards need to devote more attention to a topic about which research participants are interested.

Impact of gene variants on sex-specific regulation of human Scavenger receptor class B type 1 (SR-BI) expression in liver and association with lipid levels in a population-based study.

BACKGROUND: Several studies have noted that genetic variants of SCARB1, a lipoprotein receptor involved in reverse cholesterol transport, are associated with serum lipid levels in a sex-dependent fashion. However, the mechanism underlying this gene by sex interaction has not been explored. METHODS: We utilized both epidemiological and molecular methods to study how estrogen and gene variants interact to influence SCARB1 expression and lipid levels. Interaction between 35 SCARB1 haplotype-tagged polymorphisms and endogenous estradiol levels was assessed in 498 postmenopausal Caucasian women from the population-based Rancho Bernardo Study. We further examined associated variants with overall and SCARB1 splice variant (SR-BI and SR-BII) expression in 91 human liver tissues using quantitative real-time PCR. RESULTS: Several variants on a haplotype block spanning intron 11 to intron 12 of SCARB1 showed significant gene by estradiol interaction affecting serum lipid levels, the strongest for rs838895 with HDL-cholesterol (p=9.2x10(-4)) and triglycerides (p=1.3x10(-3)) and the triglyceride:HDL cholesterol ratio (p=2.7x10(-4)). These same variants were associated with expression of the SR-BI isoform in a sex-specific fashion, with the strongest association found among liver tissue from 52 young women<45 years old (p=0.002). CONCLUSIONS: Estrogen and SCARB1 genotype may act synergistically to regulate expression of SCARB1 isoforms and impact serum levels of HDL cholesterol and triglycerides. This work highlights the importance of considering sex-dependent effects of gene variants on serum lipid levels.

Medical education as moral formation: an Aristotelian account of medical professionalsim.

The medical professionalism movement, bolstered by many influential medical organizations and institutions, has in the last decade produced a number of conceptual definitions of professionalism and a number of concrete proposals for its measurement and teaching. These projects, however laudable, are misguided when they treat professionalism as a unitary descriptive concept rather than as a contested and therefore primarily evaluative one; when they conceive professionalism as a domain of medical practice separable in principle from other domains; and when they treat professionalism as, in principle, a specifiable goal or product of sufficiently well designed educational curricula. The logic of professionalism-as-product corresponds to the logic of techne (art or practical skill) in Aristotle's Nicomachean Ethics. Aristotle provides a cogent argument, however, that the moral excellences denoted by "professionalism" cannot be "produced" or even prespecified in the concrete; rather, they must be acquired through long practice under the careful concrete guidance of teachers who themselves embody these moral excellences. Phronesis (practical wisdom) rather than techne must therefore be the guiding logic of educational initiatives in medical professional formation, with particular emphasis on close mentorship and on the moral character both of students and of those who teach them.

Host genetics and HIV-1: the final phase?

This is a crucial transition time for human genetics in general, and for HIV host genetics in particular. After years of equivocal results from candidate gene analyses, several genome-wide association studies have been published that looked at plasma viral load or disease progression. Results from other studies that used various large-scale approaches (siRNA screens, transcriptome or proteome analysis, comparative genomics) have also shed new light on retroviral pathogenesis. However, most of the inter-individual variability in response to HIV-1 infection remains to be explained: genome resequencing and systems biology approaches are now required to progress toward a better understanding of the complex interactions between HIV-1 and its human host.

Geography and Genetics of Ecological Speciation in Pea Aphids

During ecological speciation, divergent natural selection drives evolution of ecological specialization and genetic differentiation of populations on alternate environments. Populations diverging onto the same alternate environments may be geographically widespread, so that divergence may occur at an array of locations simultaneously. Spatial variation in the process of divergence may produce a pattern of differences in divergence among locations called the Geographic Mosaic of Divergence. Diverging populations may vary in their degree of genetic differentiation and ecological specialization among locations. My dissertation examines the pattern and evolutionary processes of divergence in pea aphids (Acyrthosiphon pisum) on alfalfa (Medicago sativa) and clover (Trifolium pretense). In Chapter One, I examined differences among North American aphid populations in genetic differentiation at nuclear, sequence-based markers and in ecological specialization, measured as aphid fecundity on each host plant. In the East, aphids showed high host-plant associated ecological specialization and high genetic differentiation. In the West, aphids from clover were genetically indistinguishable from aphids on alfalfa, and aphids from clover were less specialized. Thus, the pattern of divergence differed among locations, suggesting a Geographic Mosaic of Divergence. In Chapter Two, I examined genomic heterogeneity in divergence in aphids on alfalfa and clover across North America using amplified fragment length polymorphisms (AFLPs). The degree of genetic differentiation varied greatly among markers, suggesting that divergent natural selection drives aphid divergence in all geographic locations. Three of the same genetic markers were identified as evolving under divergent selection in the eastern and western regions, and additional divergent markers were identified in the East. In Chapter Three, I investigated population structure of aphids in North America, France, and Sweden using AFLPs. Aphids on the same host plant were genetically similar across many parts of their range, so the evolution of host plant specialization does not appear to have occurred independently in every location. While aphids on alfalfa and clover were genetically differentiated in most locations, aphids from alfalfa and clover were genetically similar in both western North America and Sweden. High gene flow from alfalfa onto clover may constrain divergence in these locations.

An integrated alcohol abuse and medical treatment model for patients with hepatitis C.

BACKGROUND: Patients with chronic hepatitis C virus (HCV) infection have high rates of alcohol consumption, which is associated with progression of fibrosis and lower response rates to HCV treatment. AIMS: This prospective cohort study examined the feasibility of a 24-week integrated alcohol and medical treatment to HCV-infected patients. METHODS: Patients were recruited from a hepatology clinic if they had an Alcohol Use Disorders Identification Test score >4 for women and >8 for men, suggesting hazardous alcohol consumption. The integrated model included patients receiving medical care and alcohol treatment within the same clinic. Alcohol treatment consisted of 6 months of group and individual therapy from an addictions specialist and consultation from a study team psychiatrist as needed. RESULTS: Sixty patients were initially enrolled, and 53 patients participated in treatment. The primary endpoint was the Addiction Severity Index (ASI) alcohol composite scores, which significantly decreased by 0.105 (41.7% reduction) between 0 and 3 months (P < 0.01) and by 0.128 (50.6% reduction) between 0 and 6 months (P < 0.01) after adjusting for covariates. Alcohol abstinence was reported by 40% of patients at 3 months and 44% at 6 months. Patients who did not become alcohol abstinent had reductions in their ASI alcohol composite scores from 0.298 at baseline to 0.219 (26.8% reduction) at 6 months (P = 0.08). CONCLUSION: This study demonstrated that an integrated model of alcohol treatment and medical care could be successfully implemented in a hepatology clinic with significant favorable impact on alcohol use and abstinence among patients with chronic HCV.

Working at the interface of phylogenetics and population genetics: a biogeographical analysis of Triaenops spp. (Chiroptera: Hipposideridae).

New applications of genetic data to questions of historical biogeography have revolutionized our understanding of how organisms have come to occupy their present distributions. Phylogenetic methods in combination with divergence time estimation can reveal biogeographical centres of origin, differentiate between hypotheses of vicariance and dispersal, and reveal the directionality of dispersal events. Despite their power, however, phylogenetic methods can sometimes yield patterns that are compatible with multiple, equally well-supported biogeographical hypotheses. In such cases, additional approaches must be integrated to differentiate among conflicting dispersal hypotheses. Here, we use a synthetic approach that draws upon the analytical strengths of coalescent and population genetic methods to augment phylogenetic analyses in order to assess the biogeographical history of Madagascar's Triaenops bats (Chiroptera: Hipposideridae). Phylogenetic analyses of mitochondrial DNA sequence data for Malagasy and east African Triaenops reveal a pattern that equally supports two competing hypotheses. While the phylogeny cannot determine whether Africa or Madagascar was the centre of origin for the species investigated, it serves as the essential backbone for the application of coalescent and population genetic methods. From the application of these methods, we conclude that a hypothesis of two independent but unidirectional dispersal events from Africa to Madagascar is best supported by the data.