Genetic dissection of the human leukocyte antigen region by use of haplotypes of Tasmanians with multiple sclerosis

Association of multiple sclerosis (MS) with the human leukocyte antigen (HLA) class II haplotype DRB1*1501-DQB1*0602 is the most consistently replicated finding of genetic studies of the disease. However, the high level of linkage disequilibrium (LD) in the HLA region has hindered the identification of other loci that single-marker tests for association are unlikely to resolve. In order to address this issue, we generated haplotypes spanning 14.754 Mb (5 cM) across the entire HLA region. The haplotypes, which were inferred by genotyping relatives of 152 patients with MS and 105 unaffected control subjects of Tasmanian ancestry, define a genomic segment from D6S276 to D6S291, including 13 microsatellite markers integrated with allele-typing data for DRB1 and DQB1. Association to the DRB1*1501-DQB1*0602 haplotype was replicated. In addition, we found that the class I/extended class I region, defined by a genomic segment of approximately 400 kb between MOGCA and D6S265, harbors genes that independently increase risk of, or provide protection from, MS. Log-linear modeling analysis of constituent haplotypes that represent genomic regions containing class I (MOGCA-D6S265), class III (TNFa-TNFd-D6S273), and class II (DRB1-DQB1) genes indicated that having class I and class II susceptibility variants on the same haplotype provides an additive effect on risk. Moreover, we found no evidence for a disease locus in the class III region defined by a 150-kb genomic segment containing the TNF locus and 14 other genes. A global overview of LD performed using GOLD identified two discrete blocks of LD in the HLA region that correspond well with previous findings. We propose that the analysis of haplotypes, by use of the types of approaches outlined in the present article, should make it possible to more accurately define the contribution of the HLA to MS.     

Mapping multiple QTL of different effects: comparison of a simple sequential testing strategy and multiple QTL mapping

The aim of this study was to explore, by computer simulation, the mapping of QTLs in a realistic but complex situation of many (linked) QTLs with different effects, and to compare two QTL mapping methods. A novel method to dissect genetic variation on multiple chromosomes using molecular markers in backcross and F₂ populations derived from inbred lines was suggested, and its properties tested using simulations. The rationale for this sequential testing method was to explicitly test for alternative genetic models. The method consists of a series of four basic statistical tests to decide whether variance was due to a single QTL, two QTLs, multiple QTLs, or polygenes, starting with a test to detect genetic variance associated with a particular chromosome. The method was able to distinguish between different QTL configurations, in that the probability to `detect' the correct model was high, varying from 0.75 to 1. For example, for a backcross population of 200 and an overall heritability of 50%, in 78% of replicates a polygenic model was detected when that was the underlying true model. To test the method for multiple chromosomes, QTLs were simulated on 10 chromosomes, following a geometric series of allele effects, assuming positive alleles were in coupling in the founder lines For these simulations, the sequential testing method was compared to the established Multiple QTL Mapping (MQM) method. For a backcross population of 400 individuals, power to detect genetic variance was low with both methods when the heritability was 0.40. For example, the power to detect genetic variation on a chromosome on which 6 QTLs explained 12.6% of the genetic variance, was less than 60% for both methods. For a large heritability (0.90), the power of MQM to detect genetic variance and to dissect QTL configurations was generally better, due to the simultaneous fitting of markers on all chromosomes. It is concluded that when testing different QTL configurations on a single chromosome using the sequential testing procedure, regions of other chromosomes which explain a significant amount of variation should be fitted in the model of analysis. This study reinforces the need for large experiments in plants and other species if the aim of a genome scan is to dissect quantitative genetic variation.     

An algorithm for inferring complex haplotypes in a region of copy-number variation

Recent studies have extensively examined the large-scale genetic variants in the human genome known as copy-number variations (CNVs), and the universality of CNVs in normal individuals, along with their functional importance, has been increasingly recognized. However, the absence of a method to accurately infer alleles or haplotypes within a CNV region from high-throughput experimental data hampers the finer analyses of CNV properties and applications to disease-association studies. Here we developed an algorithm to infer complex haplotypes within a CNV region by using data obtained from high-throughput experimental platforms. We applied this algorithm to experimental data and estimated the population frequencies of haplotypes that can yield information on both sequences and numbers of DNA copies. These results suggested that the analysis of such complex haplotypes is essential for accurately detecting genetic differences within a CNV region between population groups.     

Limitations of a two-step moment method for mapping linked multiple QTL

As an alternative to multiple-interval mapping a two-step moment method was recently proposed to map linked multiple quantitative trait loci (QTLs). The advantage of this moment method was supposed to be its simplicity and computational efficiency, especially in detecting closely linked QTLs within a marker bracket, but also in mapping QTLs in different marker intervals. Using simulations it is shown that the two-step moment method may give poor results compared with multiple-interval mapping, irrespective of whether the QTLs are in the same or in different marker intervals, especially if linked QTLs are in repulsion. The criteria of comparison are number of identified QTLs, likelihood ratio test statistics, means and empirical standard errors of the QTL position and QTL effects estimates, and the accuracy of the residual variance estimates. Further, the joint conditional probabilities of QTL genotypes for two putative QTLs within a marker interval were derived and compared with the unmodified approach ignoring the non-independence of the conditional probabilities.     

Feline tetherin is characterized by a short N-terminal region and is counteracted by the feline immunodeficiency virus envelope glycoprotein

Tetherin (BST2) is the host cell factor that blocks the particle release of some enveloped viruses. Two putative feline tetherin proteins differing at the level of the N-terminal coding region have recently been described and tested for their antiviral activity. By cloning and comparing the two reported feline tetherins (called here cBST2(504) and cBST2*) and generating specific derivative mutants, this study provides evidence that feline tetherin has a shorter intracytoplasmic domain than those of other known homologues. The minimal tetherin promoter was identified and assayed for its ability to drive tetherin expression in an alpha interferon-inducible manner. We also demonstrated that cBST2(504) is able to dimerize, is localized at the cellular membrane, and impairs human immunodeficiency virus type 1 (HIV-1) particle release, regardless of the presence of the Vpu antagonist accessory protein. While cBST2(504) failed to restrict wild-type feline immunodeficiency virus (FIV) egress, FIV mutants, bearing a frameshift at the level of the envelope-encoding region, were potently blocked. The transient expression of the FIV envelope glycoprotein was able to rescue mutant particle release from feline tetherin-positive cells but did not antagonize human BST2 activity. Moreover, cBST2(504) was capable of specifically immunoprecipitating the FIV envelope glycoprotein. Finally, cBST2(504) also exerted its function on HIV-2 ROD10 and on the simian immunodeficiency virus SIVmac239. Taken together, these results show that feline tetherin does indeed have a short N-terminal region and that the FIV envelope glycoprotein is the predominant factor counteracting tetherin restriction.     

Genetic structure of a QTL hotspot on chromosome 2 in sweet cherry indicates positive selection for favorable haplotypes

A genomic region of particular interest for sweet cherry (Prunus avium L.) breeding is a quantitative trait locus (QTL) “hotspot” on chromosome 2. QTLs for fruit size, firmness, sweetness, and flowering time are reported to map to this region. An understanding of genetic diversity, allele sources, linkage relationships, and historical recombinations is critical to enable the combining of favorable alleles at multiple loci. The objectives of this study were to characterize, visualize, and interpret the genetic structure of this previously identified QTL hotspot within North American sweet cherry breeding germplasm, using a pedigree-based haploblocking approach. Across the 29.4 cM (6.3 Mbp) region defined by single nucleotide polymorphism (SNP) information from the RosBREED cherry 6K SNP array v1, a total of 12 recombination events falling into six inter-marker regions were traced within the pedigree of elite and wild germplasm (n = 55). These recombinations defined five haploblocks containing 5–15 markers and exhibiting 7–11 haplotypes each. Over the entire QTL hotspot, 30 extended haplotypes were identified for which parental gametes could be determined. When the haploblocks and their haplotypes were used to explore genetic diversity, ancestry, and recombination patterns, and then integrated with previous QTL results for fruit size, the results indicated that favorable alleles at this QTL hotspot are under positive selection in breeding. The genetic framework provided by a haploblock approach and knowledge of haplotype-level diversity sets the stage for assigning breeding utility to these haplotypes.     

Identification of a 1.1-Mb region for a carcass weight QTL onbovine Chromosome 14

We previously mapped a quantitative trait locus for carcass weight, designated Carcass Weight-1 (CW-1), to bovine Chromosome 14 using a purebred Wagyu pedigree based on progeny design analysis. To refine the critical region within 8.1 cM flanked by microsatellites BMS1941 and INRA094, we constructed a bacterial artificial chromosome (BAC) contig composed of 60 tiled BAC clones and prepared a high-density physical map including 80 microsatellites, of which 55 were developed in this study. We conducted linkage disequilibrium (LD) mapping in the CW-1 region with 47 microsatellites using paternal half-sib pedigrees whose sires exhibited homozygous CW-1 Q alleles in the region. The LD mapping study significantly narrowed the CW-1 locus to the 1.1-Mb region between microsatellites DIK7012 and DIK7020. Finally, we surveyed the 1.1-Mb-region genotypes of 1700 steers from 11 bulls having −/−, Q/−, or Q/Q alleles in the region, and we examined the effect of the CW-1 Q allele on carcass weight. The presence of the first Q increased carcass weight by 23.6 kg (95% confidence interval [CI], 17.6–29.5 kg), and the second Q increased carcass weight an additional 15.2 kg (95% CI, 10.7–19.7 kg). These results indicate the presence of a gene responsible for carcass weight within the 1.1-Mb region.Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank databases under the accession number AB204952-AB205006.Kazunori Mizoshita, and Atsushi Takano contributed equally to this work.     

DeFries-Fulker multiple regression analysis of sibship QTL data: a SAS macro

A SAS macro package for performing multipoint QTL mapping using the DeFries-Fulker multiple regression method is presented.     

R/qtl: high-throughput multiple QTL mapping

MOTIVATION: R/qtl is free and powerful software for mapping and exploring quantitative trait loci (QTL). R/qtl provides a fully comprehensive range of methods for a wide range of experimental cross types. We recently added multiple QTL mapping (MQM) to R/qtl. MQM adds higher statistical power to detect and disentangle the effects of multiple linked and unlinked QTL compared with many other methods. MQM for R/qtl adds many new features including improved handling of missing data, analysis of 10,000 s of molecular traits, permutation for determining significance thresholds for QTL and QTL hot spots, and visualizations for cis-trans and QTL interaction effects. MQM for R/qtl is the first free and open source implementation of MQM that is multi-platform, scalable and suitable for automated procedures and large genetical genomics datasets. AVAILABILITY: R/qtl is free and open source multi-platform software for the statistical language R, and is made available under the GPLv3 license. R/qtl can be installed from http://www.rqtl.org/. R/qtl queries should be directed at the mailing list, see http://www.rqtl.org/list/. CONTACT: kbroman@biostat.wisc.edu.     

The familial adenomatous polyposis region exhibits many different haplotypes

In the present study, we used five different polymorphic markers to construct the haplotype at the adenomatous polyposis coli (APC) locus in families with familial adenomatous polyposis (FAP) and in the normal Italian population. Non-ambiguous haplotypes were reconstructed from 246 normal chromosomes and 65 FAP chromosomes. In the control population, the four polymorphisms intragenic to APC gave rise to 16 haplotypes, the most common of which (II and XV) accounted for over 50% of all chromosomes. In FAP patients, 13 haplotypes were found but their distribution was not statistically different from normal subjects. Eighty complete chromosomal haplotypes (many fewer than the theoretical maximum of 208) for the five polymorphic sites assayed were observed in the control population, 35 being found in the FAP patients. We compared the distribution of these haplotypes within the two groups; no statistically significant differences between normal and FAP chromosomes were found. The elevated heterogeneity of FAP chromosomes was clearly confirmed by the observation that 19 patients who carried one or other of the two most common APC mutations (nt 3183 and nt 3927) showed 18 different haplotypes. On the basis of these results, we were not able to identify a founder FAP chromosome. Various mechanisms are presented to explain this observation. Received: 5 November 1997 / Accepted: 3 February 1998     

Characterization of the class III region in different MHC haplotypes by pulsed-field gel electrophoresis

The class III region of the human major histocompatibility complex (MHC) in sevenHLA haplotypes has been analyzed using pulsed-field gel electrophoresis (PFGE), restriction enzymes that cut genomic DNA infrequently, and Southern blotting. In particular, extensive mapping with the enzymeBss HII, which generates fragments in the size range 8–280 kilobases (kb), has revealed that in the haplotypes studied the DNA content of the class III region does not appear to vary other than as previously observed at theC4 andCYP21 loci.     

Generation of a sequence characterized amplified region probe for authentication of crocodilian species

A 209-base pair (bp) crocodilian-specific sequence characterized amplified region (SCAR) was identified from a 425-bp randomly amplified polymorphic DNA (RAPD) fragment. The 209-bp SCAR was produced from amplifications of DNA extracted from fresh and/or dry meat samples from at least three species of Crocodylus, Caiman crocodylus, and Alligator mississippiensis. No amplification was observed from DNA of other common animal species. The use of SCAR opens the way for quick authentication of crocodilian samples for conservation biology and trade regulation.     

Identification of a 3.7-Mb region for a marbling QTL on bovine chromosome 4 by identical-by-descent and association analysis

QTL mapping for growth and carcass traits was performed using a paternal half-sib family composed of 325 Japanese Black cattle offspring. Nine QTL were detected at the 1% chromosome-wise significance level at a false discovery rate of less than 0.1. These included two QTL for marbling on BTA 4 and 18, two QTL for carcass weight on BTA 14 and 24, two QTL for longissimus muscle area on BTA 1 and 4, two QTL for subcutaneous fat thickness on BTA 1 and 15 and one QTL for rib thickness on BTA 6. Although the marbling QTL on BTA 4 has been replicated with significant linkages in two Japanese Black cattle sires, the three Q (more marbling) haplotypes, each inherited maternally, were apparently different. To compare the three Q haplotypes in more detail, high-density microsatellite markers for the overlapping regions were developed within the 95% CIs (65 markers in 44–78 cM). A detailed haplotype comparison indicated that a small region (<3.7 Mb) around 46 cM was shared between the Qs of the two sires, whose dams were related. An association of this region with marbling was shown by a regression analysis using the local population, in which the two sires were produced and this was confirmed by an association study using a population collected throughout Japan. These results strongly suggest that the marbling QTL on BTA 4 is located in the 3.7-Mb region at around 46 cM.     

Analysis of extended HLA haplotypes in multiple sclerosis and narcolepsy families confirms a predisposing effect for the class I region in Tasmanian MS patients

Human leucocyte antigen (HLA)-DRB1*15 is associated with predisposition to multiple sclerosis (MS), although conjecture surrounds the possible involvement of an alternate risk locus in the class I region of the HLA complex. We have shown previously that an alternate MS risk allele(s) may be encompassed by the telomerically extended DRB1*15 haplotype, and here, we have attempted to map the putative variant. Thirteen microsatellite markers encompassing a 6.79-megabase (D6S2236-G51152) region, and the DRB1 and DQB1 genes, were genotyped in 166 MS simplex families and 104 control families from the Australian State of Tasmania and 153 narcolepsy simplex families (trios) from the USA. Complementary approaches were used to investigate residual predisposing effects of microsatellite alleles comprising the extended DRB1*15 haplotype taking into account the strong predisposing effect of DRB1*15: (1) Disease association of the extended DRB1*15 haplotype was compared for MS and narcolepsy families--predisposing effects were observed for extended class I microsatellite marker alleles in MS families, but not narcolepsy families; (2) disease association of the extended DRB1*15 haplotype was investigated after conditioning MS and control haplotypes on the absence of DRB1*15--a significant predisposing effect was observed for a 627-kb haplotype (D6S258 allele 8-MOGCA allele 4; MOG, myelin oligodendrocyte glycoprotein) spanning the extended class I region. MOGCA allele 4 displayed the strongest predisposing effect in DRB1*15-conditioned haplotypes (p = 0.0006; OR 2.83 [1.54-5.19]). Together, these data confirm that an alternate MS risk locus exists in the extended class I region in Tasmanian MS patients independent of DRB1*15.     

Development of sequence characterized amplified region (SCAR) primers for the detection of Phyto.5.2, a major QTL for resistance to Phytophthora capsici Leon. in pepper

Phytophthora capsici causes devastating disease on many crop species, including Capsicum. Resistance in Capsicum annuum is genetically and physiologically complex. A panel of Capsicum germplasm that included genotypes from both C. annuum and C. chinense showing highly resistant, highly susceptible and intermediate or tolerant responses to the pathogen, respectively, was screened with a series of randomly amplified polymorphic sequence primers to determine which genomic regions contribute to the highest level of resistance. One primer, OpD04, amplified a single band only in those C. annuum and C. chinense genotypes showing the highest level of resistance. The amplified product was cloned, sequenced and used to design longer primers in order to generate a sequence characterized amplified region marker which was then mapped in a reference mapping population and a screened population segregating for resistance to P. capsici. These primers were observed to define a locus on pepper chromosome 5 tightly linked to Phyto.5.2, one of six quantitative trait loci (QTL) previously reported to contribute to P. capsici resistance. These results indicate that the Phyto.5.2 QTL may be widely distributed in highly resistant germplasm and provide improved resolution for this QTL. This work also defines the first breeding tools for this system, allowing for the rapid selection of genotypes likely to be highly resistant to P. capsici.     

Using haplotypes for the prediction of allelic identity to fine-map QTL: characterization and properties

Background

Numerous methods have been developed over the last decade to predict allelic identity at unobserved loci between pairs of chromosome segments along the genome. These loci are often unobserved positions tested for the presence of quantitative trait loci (QTL). The main objective of this study was to understand from a theoretical standpoint the relation between linkage disequilibrium (LD) and allelic identity prediction when using haplotypes for fine mapping of QTL. In addition, six allelic identity predictors (AIP) were also compared in this study to determine which one performed best in theory and application.

Results

A criterion based on a simple measure of matrix distance was used to study the relation between LD and allelic identity prediction when using haplotypes. The consistency of this criterion with the accuracy of QTL localization, another criterion commonly used to compare AIP, was evaluated on a set of real chromosomes. For this set of chromosomes, the criterion was consistent with the mapping accuracy of a simulated QTL with either low or high effect. As measured by the matrix distance, the best AIP for QTL mapping were those that best captured LD between a tested position and a QTL. Moreover the matrix distance between a tested position and a QTL was shown to decrease for some AIP when LD increased. However, the matrix distance for AIP with continuous predictions in the [0,1] interval was algebraically proven to decrease less rapidly up to a lower bound with increasing LD in the simplest situations, than the discrete predictor based on identity by state between haplotypes (IBS _hap), for which there was no lower bound. The expected LD between haplotypes at a tested position and alleles at a QTL is a quantity that increases naturally when the tested position gets closer to the QTL. This behavior was demonstrated with pig and unrelated human chromosomes.

Conclusions

When the density of markers is high, and therefore LD between adjacent loci can be assumed to be high, the discrete predictor IBS _hap is recommended since it predicts allele identity correctly when taking LD into account.