Genetic association studies under the population stratification, family pedigree and application to genome-wide association studies

This dissertation has three separate parts: the first part deals with the general pedigree association testing incorporating continuous covariates; the second part deals with the association tests under population stratification using the conditional likelihood tests; the third part deals with the genome-wide association studies based on the real rheumatoid arthritis (RA) disease data sets from Genetic Analysis Workshop 16 (GAW16) problem 1. Many statistical tests are developed to test the linkage and association using either case-control status or phenotype covariates for family data structure, separately. Those univariate analyses might not use all the information coming from the family members in practical studies. On the other hand, the human complex disease do not have a clear inheritance pattern, there might exist the gene interactions or act independently. In part I, the new proposed approach MPDT is focused on how to use both the case control information as well as the phenotype covariates. This approach can be applied to detect multiple marker effects. Based on the two existing popular statistics in family studies for case-control and quantitative traits respectively, the new approach could be used in the simple family structure data set as well as general pedigree structure. The combined statistics are calculated using the two statistics; A permutation procedure is applied for assessing the p-value with adjustment from the Bonferroni for the multiple markers. We use simulation studies to evaluate the type I error rates and the powers of the proposed approach. Our results show that the combined test using both case-control information and phenotype covariates not only has the correct type I error rates but also is more powerful than the other existing methods. For multiple marker interactions, our proposed method is also very powerful. Selective genotyping is an economical strategy in detecting and mapping quantitative trait loci in the genetic dissection of complex disease. When the samples arise from different ethnic groups or an admixture population, all the existing selective genotyping methods may result in spurious association due to different ancestry distributions. The problem can be more serious when the sample size is large, a general requirement to obtain sufficient power to detect modest genetic effects for most complex traits. In part II, I describe a useful strategy in selective genotyping while population stratification is present. Our procedure used a principal component based approach to eliminate any effect of population stratification. The paper evaluates the performance of our procedure using both simulated data from an early study data sets and also the HapMap data sets in a variety of population admixture models generated from empirical data. There are one binary trait and two continuous traits in the rheumatoid arthritis dataset of Problem 1 in the Genetic Analysis Workshop 16 (GAW16): RA status, AntiCCP and IgM. To allow multiple traits, we suggest a set of SNP-level F statistics by the concept of multiple-correlation to measure the genetic association between multiple trait values and SNP-specific genotypic scores and obtain their null distributions. Hereby, we perform 6 genome-wide association analyses using the novel one- and two-stage approaches which are based on single, double and triple traits. Incorporating all these 6 analyses, we successfully validate the SNPs which have been identified to be responsible for rheumatoid arthritis in the literature and detect more disease susceptibility SNPs for follow-up studies in the future. Except for chromosome 13 and 18, each of the others is found to harbour susceptible genetic regions for rheumatoid arthritis or related diseases, i.e., lupus erythematosus. This topic is discussed in part III.

Genetic Distance between Species Predicts Novel Trait Expression in Their Hybrids

Interspecific hybridization can generate transgressive hybrid phenotypes with extreme trait values exceeding the combined range of the parental species. Such variation can enlarge the working surface for natural selection, and may facilitate the evolution of novel adaptations where ecological opportunity exists. The number of quantitative trait loci fixed for different alleles in different species should increase with time since speciation. If transgression is caused by complementary gene action or epistasis, hybrids between more distant species should be more likely to display transgressive phenotypes. To test this prediction we collected data on transgression frequency from the literature, estimated genetic distances between the hybridizing species from gene sequences, and calculated the relationship between the two using phylogenetically controlled methods. We also tested if parental phenotypic divergence affected the occurrence of transgression. We found a highly significant positive correlation between transgression frequency and genetic distance in eudicot plants explaining 43% of the variance in transgression frequency. In total, 36% of the measured traits were transgressive. The predicted effect of time since speciation on transgressive segregation was unconfounded by the potentially conflicting effects of phenotypic differentiation between species. Our analysis demonstrates that the potential impact hybridization may have on phenotypic evolution is predictable from the genetic distance between species.

Accumulating mutations in series of haplotypes at the KIT and MITF loci are major determinants of white markings in Franches-Montagnes horses.

Coat color and pattern variations in domestic animals are frequently inherited as simple monogenic traits, but a number are known to have a complex genetic basis. While the analysis of complex trait data remains a challenge in all species, we can use the reduced haplotypic diversity in domestic animal populations to gain insight into the genomic interactions underlying complex phenotypes. White face and leg markings are examples of complex traits in horses where little is known of the underlying genetics. In this study, Franches-Montagnes (FM) horses were scored for the occurrence of white facial and leg markings using a standardized scoring system. A genome-wide association study (GWAS) was performed for several white patterning traits in 1,077 FM horses. Seven quantitative trait loci (QTL) affecting the white marking score with p-values p≤10(-4) were identified. Three loci, MC1R and the known white spotting genes, KIT and MITF, were identified as the major loci underlying the extent of white patterning in this breed. Together, the seven loci explain 54% of the genetic variance in total white marking score, while MITF and KIT alone account for 26%. Although MITF and KIT are the major loci controlling white patterning, their influence varies according to the basic coat color of the horse and the specific body location of the white patterning. Fine mapping across the MITF and KIT loci was used to characterize haplotypes present. Phylogenetic relationships among haplotypes were calculated to assess their selective and evolutionary influences on the extent of white patterning. This novel approach shows that KIT and MITF act in an additive manner and that accumulating mutations at these loci progressively increase the extent of white markings.

Evolution of dominance under frequency-dependent intraspecific competition

A population-genetic analysis is performed of a two-locus two-allele model, in which the primary locus has a major effect on a quantitative trait that is under frequency-dependent disruptive selection caused by intraspecific competition for a continuum of resources. The modifier locus determines the degree of dominance at the trait level. We establish the conditions when a modifier allele can invade and when it becomes fixed if sufficiently frequent. In general, these are not equivalent because an unstable internal equilibrium may exist and the condition for successful invasion of the modifier is more restrictive than that for eventual fixation from already high frequency. However, successful invasion implies global fixation, i.e., fixation from any initial condition. Modifiers of large effect can become fixed, and also invade, in a wider parameter range than modifiers of small effect. We also study modifiers with a direct, frequency-independent deleterious fitness effect. We show that they can invade if they induce a sufficiently high level of dominance and if disruptive selection on the ecological trait is strong enough. For deleterious modifiers, successful invasion no longer implies global fixation because they can become stuck at an intermediate frequency due to a stable internal equilibrium. Although the conditions for invasion and for fixation if sufficiently frequent are independent of the linkage relation between the two loci, the rate of spread depends strongly on it. The present study provides further support to the view that evolution of dominance may be an efficient mechanism to remove unfit heterozygotes that are maintained by balancing selection. It also demonstrates that an invasion analysis of mutants of very small effect is insufficient to obtain a full understanding of the evolutionary dynamics under frequency-dependent selection.

EVOLUTION OF DOMINANCE UNDER FREQUENCY-DEPENDENT INTRASPECIFIC COMPETITION IN AN ASSORTATIVELY MATING POPULATION

We study the evolution of higher levels of dominance as a response to negative frequency-dependent selection. In contrast to previous studies, we focus on the effect of assortative mating on the evolution of dominance under frequency-dependent intraspecific competition. We analyze a two-locus two-allele model, in which the primary locus has a major effect on a quantitative trait that is under a mixture of frequency-independent stabilizing selection, density-dependent selection, and frequency-dependent selection caused by intraspecific competition for a continuum of resources. The second (modifier) locus determines the degree of dominance at the trait level. Additionally, the population mates assortatively with respect to similarities in the ecological trait. Our analysis shows that the parameter region in which dominance can be established decreases if small levels of assortment are introduced. In addition, the degree of dominance that can be established also decreases. In contrast, if assortment is intermediate, sexual selection for extreme types can be established, which leads to evolution of higher levels of dominance than under random mating. For modifiers with large effects, intermediate levels of assortative mating are most favorable for the evolution of dominance. For large modifiers, the speed of fixation can even be higher for intermediate levels of assortative mating than for random mating.

Microscale mapping of alteration conditions and potential biosignatures in basaltic-ultramafic rocks on early Earth and beyond

Subseafloor environments preserved in Archean greenstone belts provide an analogue for investigating potential subsurface habitats on Mars. The c. 3.5-3.4 Ga pillow lava metabasalts of the mid-Archean Barberton greenstone belt, South Africa, have been argued to contain the earliest evidence for microbial subseafloor life. This includes candidate trace fossils in the form of titanite microtextures, and sulfur isotopic signatures of pyrite preserved in metabasaltic glass of the c. 3.472 Ga Hooggenoeg Formation. It has been contended that similar microtextures in altered martian basalts may represent potential extraterrestrial biosignatures of microbe-fluid-rock interaction. But despite numerous studies describing these putative early traces of life, a detailed metamorphic characterization of the microtextures and their host alteration conditions in the ancient pillow lava metabasites is lacking. Here, we present a new nondestructive technique with which to study the in situ metamorphic alteration conditions associated with potential biosignatures in mafic-ultramafic rocks of the Hooggenoeg Formation. Our approach combines quantitative microscale compositional mapping by electron microprobe with inverse thermodynamic modeling to derive low-temperature chlorite crystallization conditions. We found that the titanite microtextures formed under subgreenschist to greenschist facies conditions. Two chlorite temperature groups were identified in the maps surrounding the titanite microtextures and record peak metamorphic conditions at 315 ± 40°C (XFe3+(chlorite) = 25-34%) and lower-temperature chlorite veins/microdomains at T = 210 ± 40°C (lower XFe3+(chlorite) = 40-45%). These results provide the first metamorphic constraints in textural context on the Barberton titanite microtextures and thereby improve our understanding of the local preservation conditions of these potential biosignatures. We suggest that this approach may prove to be an important tool in future studies to assess the biogenicity of these earliest candidate traces of life on Earth. Furthermore, we propose that this mapping approach could also be used to investigate altered mafic-ultramafic extraterrestrial samples containing candidate biosignatures.

A Non-Synonymous HMGA2 Variant Decreases Height in Shetland Ponies and Other Small Horses

The identification of quantitative trait loci (QTL) such as height and their underlying causative variants is still challenging and often requires large sample sizes. In humans hundreds of loci with small effects control the heritable portion of height variability. In domestic animals, typically only a few loci with comparatively large effects explain a major fraction of the heritability. We investigated height at withers in Shetland ponies and mapped a QTL to ECA 6 by genome-wide association (GWAS) using a small cohort of only 48 animals and the Illumina equine SNP70 BeadChip. Fine-mapping revealed a shared haplotype block of 793 kb in small Shetland ponies. The HMGA2 gene, known to be associated with height in horses and many other species, was located in the associated haplotype. After closing a gap in the equine reference genome we identified a non-synonymous variant in the first exon of HMGA2 in small Shetland ponies. The variant was predicted to affect the functionally important first AT-hook DNA binding domain of the HMGA2 protein (c.83G>A; p.G28E). We assessed the functional impact and found impaired DNA binding of a peptide with the mutant sequence in an electrophoretic mobility shift assay. This suggests that the HMGA2 variant also affects DNA binding in vivo and thus leads to reduced growth and a smaller stature in Shetland ponies. The identified HMGA2 variant also segregates in several other pony breeds but was not found in regular-sized horse breeds. We therefore conclude that we identified a quantitative trait nucleotide for height in horses.

Single nucleotide polymorphism (SNP) selection for genotype-phenotype association studies

With hundreds of single nucleotide polymorphisms (SNPs) in a candidate gene and millions of SNPs across the genome, selecting an informative subset of SNPs to maximize the ability to detect genotype-phenotype association is of great interest and importance. In addition, with a large number of SNPs, analytic methods are needed that allow investigators to control the false positive rate resulting from large numbers of SNP genotype-phenotype analyses. This dissertation uses simulated data to explore methods for selecting SNPs for genotype-phenotype association studies. I examined the pattern of linkage disequilibrium (LD) across a candidate gene region and used this pattern to aid in localizing a disease-influencing mutation. The results indicate that the r2 measure of linkage disequilibrium is preferred over the common D′ measure for use in genotype-phenotype association studies. Using step-wise linear regression, the best predictor of the quantitative trait was not usually the single functional mutation. Rather it was a SNP that was in high linkage disequilibrium with the functional mutation. Next, I compared three strategies for selecting SNPs for application to phenotype association studies: based on measures of linkage disequilibrium, based on a measure of haplotype diversity, and random selection. The results demonstrate that SNPs selected based on maximum haplotype diversity are more informative and yield higher power than randomly selected SNPs or SNPs selected based on low pair-wise LD. The data also indicate that for genes with small contribution to the phenotype, it is more prudent for investigators to increase their sample size than to continuously increase the number of SNPs in order to improve statistical power. When typing large numbers of SNPs, researchers are faced with the challenge of utilizing an appropriate statistical method that controls the type I error rate while maintaining adequate power. We show that an empirical genotype based multi-locus global test that uses permutation testing to investigate the null distribution of the maximum test statistic maintains a desired overall type I error rate while not overly sacrificing statistical power. The results also show that when the penetrance model is simple the multi-locus global test does as well or better than the haplotype analysis. However, for more complex models, haplotype analyses offer advantages. The results of this dissertation will be of utility to human geneticists designing large-scale multi-locus genotype-phenotype association studies. ^

The genetic contribution of the major histocompatibility complex to bleomycin-induced lung injury

Pulmonary fibrosis (PF) is the result of a variety of environmental and cancer treatment related insults and is characterized by excessive deposition of collagen. Gas exchange in the alveoli is impaired as the normal lung becomes dense and collapsed leading to a loss of lung volume. It is now accepted that lung injury and fibrosis are in part genetically regulated. ^ Bleomycin is a chemotherapeutic agent used for testicular cancer and lymphomas that induces significant pulmonary toxicity. We delivered bleomycin to mice subcutaneously via a miniosmotic pump in order to elicit lung injury (LI) and quantified the %LI morphometrically using video imaging software. We previously identified a quantitative trait loci, Blmpf-1(LOD=17.4), in the Major Histocompatibility Complex (MHC), but the exact genetic components involved have remained unknown. ^ In the current studies, Blmpf-1 was narrowed to an interval spanning 31.9-32.9Mb on Chromosome 17 using MHC Congenic mice. This region includes the MHC Class II and III genes, and is flanked by the TNF-alpha super locus and MHC Class I genes. Knockout mice of MHC Class I genes (B2mko), MHC Class II genes (Cl2ko), and TNF-alpha (TNF-/-) and its receptors (p55-/-, p75-/-, and p55/p75-/-) were treated with bleomycin in order to ascertain the role of these genes in the pathogenesis of lung injury. ^ Cl2ko mice had significantly better survival and %LI when compared to treated background BL/6 (B6, P<.05). In contrast, B2mko showed no differences in survival or %LI compared to B6. This suggests that the MHC Class II locus contains susceptibility genes for bleomycin-induced lung injury. ^ TNF-alpha, a Class III gene, was examined and it was found that TNF-/- and p55-/- mice had higher %LI and lower survival when compared to B6 (P<.05). In contrast, p75-/- mice had significantly reduced %LI when compared to TNF-/-, p55-/-, and B6 mice as well as higher survival (P<.01). These data contradict the current paradigm that TNF-alpha is a profibrotic mediator of lung injury and suggest a novel and distinct role for the p55 and p75 receptors in mediating lung injury. ^

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

Importance of epistasis as the genetic basis of heterosis in an elite rice hybrid

The genetic basis of heterosis was investigated in an elite rice hybrid by using a molecular linkage map with 150 segregating loci covering the entire rice genome. Data for yield and three traits that were components of yield were collected over 2 years from replicated field trials of 250 F2:3 families. Genotypic variations explained from about 50% to more than 80% of the total variation. Interactions between genotypes and years were small compared with the main effects. A total of 32 quantitative trait loci (QTLs) were detected for the four traits; 12 were observed in both years and the remaining 20 were detected in only one year. Overdominance was observed for most of the QTLs for yield and also for a few QTLs for the component traits. Correlations between marker heterozygosity and trait expression were low, indicating that the overall heterozygosity made little contribution to heterosis. Digenic interactions, including additive by additive, additive by dominance, and dominance by dominance, were frequent and widespread in this population. The interactions involved large numbers of marker loci, most of which individually were not detectable on single-locus basis; many interactions among loci were detected in both years. The results provide strong evidence that epistasis plays a major role as the genetic basis of heterosis.

Interleukin 9: A candidate gene for asthma

Asthma is a complex heritable inflammatory disorder of the airways associated with clinical signs of atopy and bronchial hyperresponsiveness. Recent studies localized a major gene for asthma to chromosome 5q31-q33 in humans. Thus, this segment of the genome represents a candidate region for genes that determine susceptibility to bronchial hyperresponsiveness and atopy in animal models. Homologs of candidate genes on human chromosome 5q31-q33 are found in four regions in the mouse genome, two on chromosome 18, and one each on chromosomes 11 and 13. We assessed bronchial responsiveness as a quantitative trait in mice and found it linked to chromosome 13. Interleukin 9 (IL-9) is located in the linked region and was analyzed as a gene candidate. The expression of IL-9 was markedly reduced in bronchial hyporesponsive mice, and the level of expression was determined by sequences within the qualitative trait locus (QTL). These data suggest a role for IL-9 in the complex pathogenesis of bronchial hyperresponsiveness as a risk factor for asthma.

The weediness of wild plants: molecular analysis of genes influencing dispersal and persistence of johnsongrass, Sorghum halepense (L.) Pers.

Many major weeds rely upon vegetative dispersal by rhizomes and seed dispersal by "shattering" of the mature inflorescence. We report molecular analysis of these traits in a cross between cultivated and wild species of Sorghum that are the probable progenitors of the major weed "johnsongrass." By restriction fragment length polymorphism mapping, variation in the number of rhizomes producing above-ground shoots was associated with three quantitative trait loci (QTLs). Variation in regrowth (ratooning) after overwintering was associated with QTLs accounting for additional rhizomatous growth and with QTLs influencing tillering. Vegetative buds that become rhizomes are similar to those that become tillers--one QTL appears to influence the number of such vegetative buds available, and additional independent genes determine whether individual buds differentiate into tillers or rhizomes. DNA markers described herein facilitate cloning of genes associated with weediness, comparative study of rhizomatousness in other Poaceae, and assessment of gene flow between cultivated and weedy sorghums--a risk that constrains improvement of sorghum through biotechnology. Cloning of "weediness" genes may create opportunities for plant growth regulation, in suppressing propagation of weeds and enhancing productivity of major forage, turf, and "ratoon" crops.

Modularisation épistatique des loci à trait quantitatif associés à la pression artérielle et identification de gènes candidats pour l’hypertension

Problématique: L’hypertension artérielle essentielle, facteur de risque majeur dans le développement des maladies cardiovasculaires, est un trait multigénique complexe dont les connaissances sur le déterminisme génétique nécessitent d’être approfondies. De nombreux loci à trait quantitatif (QTLs); soit des gènes responsables de faire varier la pression artérielle (PA), ont été identifiés chez l’humain et le modèle animal. Cependant, le mystère plane encore sur la façon dont ces gènes fonctionnent ensemble pour réguler la PA. Hypothèse et objectif: Plutôt qu’une addition de QTLs ayant chacun une action infinitésimale sur la PA, une interaction épistatique entre les gènes serait responsable du phénotype hypertendu. Ainsi, l’étude de cette épistasie entre les gènes impliqués, directement ou indirectement, dans l’homéostasie de la PA nous permettrait d’explorer de nouvelles voies de régulation moléculaire en cause dans cette maladie. Méthodes: Via la réalisation de souches congéniques de rats, où un segment chromosomique provenant d’une souche receveuse hypertendue (Dahl Salt Sensitive, SS/Jr) est remplacé par son homologue provenant d’une souche donneuse normotendue (Lewis, LEW), des QTLs peuvent être mis en évidence. Dans ce contexte, la combinaison de QTLs via la création de doubles ou multiples congéniques constitue la première démonstration fonctionnelle des interactions intergéniques. Résultats: Vingt-sept combinaisons au total nous ont menés à l’appréciation d’une modularisation des QTLs. Ces derniers ont été catégorisés selon deux principaux modules épistatiques (EMs) où les QTLs appartenant à un même EM sont épistatiques entre eux et participent à une même voie régulatrice. Les EMs/cascades agissent alors en parallèle pour réguler la PA. Grâce à l’existence de QTLs ayant des effets opposés sur la PA, nous avons pu établir l’ordre hiérarchique entre trois paires de QTLs. Cependant, lorsque cette suite régulatrice ne peut être déterminée, d’autres approches sont nécessaires. Nos travaux nous ont mené à l’identification d’un QTL situé sur le chromosome 16 du rat (C16QTL), appartenant au EM1 et qui révélerait une nouvelle voie de l’homéostasie de la PA. Le gène retinoblastoma-associated protein 140 (Rap140)/family with sequence similarity 208 member A (Fam208a), présentant une mutation non synonyme entre SS/Jr et LEW est le gène candidat le plus plausible pour représenter C16QTL. Celui-ci code pour un facteur de transcription et semblerait influencer l’expression de Solute carrier family 7 (cationic amino acid transporter, y+ system) member 12 (Slc7a12), spécifiquement et significativement sous exprimé dans les reins de la souche congénique portant C16QTL par rapport à la souche SS/Jr. Rap140/Fam208a agirait comme un inhibiteur de la transcription de Slc7a12 menant à une diminution de la pression chez Lewis. Conclusions: L’architecture complexe de la régulation de la PA se dévoile mettant en scène de nouveaux acteurs, pour la plupart inconnus pour leur implication dans la PA. L’étude de la nouvelle voie de signalisation Rap140/Fam208a - Slc7a12 nous permettra d’approfondir nos connaissances quant à l’homéostasie de la pression artérielle et de l’hypertension chez SS/Jr. À long terme, de nouveaux traitements anti-hypertenseurs, ciblant plus d’une voie de régulation à la fois, pourraient voir le jour.

A genetic linkage map in autotetraploid lucerne adapted to northern Australia, and use of the map to identify DNA markers linked to resistance to Phytophthora medicaginis

Phytophthora root rot, caused by Phytophthora medicaginis, is a major limitation to lucerne ( Medicago sativa L.) production in Australia and North America. Quantitative trait loci (QTLs) involved in resistance to P. medicaginis were identified in a lucerne backcross population of 120 individuals. A genetic linkage map was constructed for tetraploid lucerne using 50 RAPD ( randomly amplified polymorphic DNA), 104 AFLP (amplified fragment length polymorphism) markers, and one SSR ( simple sequence repeat or microsatellite) marker, which originated from the resistant parent (W116); 13 markers remain unlinked. The linkage map contains 18 linkage groups covering 2136.5 cM, with an average distance of 15.0 cM between markers. Four of the linkage groups contained only either 2 or 3 markers. Using duplex markers and repulsion phase linkages the map condensed to 7 homology groups and 2 unassigned linkage groups. Three regions located on linkage groups 2, 14, and 18, were identified as associated with root reaction and the QTLs explained 6 - 15% of the phenotypic variation. The research also indicates that different resistance QTLs are involved in conferring resistance in different organs. Two QTLs were identified as associated with disease resistance expressed after inoculation of detached leaves. The marker, W11-2 on group 18, identified as associated with root reaction, contributed 7% of the phenotypic variation in leaf response in our population. This marker appears to be linked to a QTL encoding a resistance factor contributing to both root and leaf reaction. One other QTL, not identified as associated with root reaction, was positioned on group 1 and contributed to 6% of the variation. This genetic linkage map provides an entry point for future molecular-based improvement of lucerne in Australia, and markers linked to the QTLs we have reported should be useful for marker-assisted selection for partial resistance to P. medicaginis in lucerne.