Human QTL linkage mapping

Hunting for Genes Underlying Emotionality in the Laboratory Rat: Maps, Tools and Traps

Scientists have systematically investigated the hereditary bases of behaviors since the 19th century, moved by either evolutionary questions or clinically-motivated purposes. The pioneer studies on the genetic selection of laboratory animals had already indicated, one hundred years ago, the immense complexity of analyzing behaviors that were influenced by a large number of small-effect genes and an incalculable amount of environmental factors. Merging Mendelian, quantitative and molecular approaches in the 1990s made it possible to map specific rodent behaviors to known chromosome regions. From that point on, Quantitative Trait Locus (QTL) analyses coupled with behavioral and molecular techniques, which involved in vivo isolation of relevant blocks of genes, opened new avenues for gene mapping and characterization. This review examines the QTL strategy applied to the behavioral study of emotionality, with a focus on the laboratory rat. We discuss the challenges, advances and limitations of the search for Quantitative Trait Genes (QTG) playing a role in regulating emotionality. For the past 25 years, we have marched the long journey from emotionality-related behaviors to genes. In this context, our experiences are used to illustrate why and how one should move forward in the molecular understanding of complex psychiatric illnesses. The promise of exploring genetic links between immunological and emotional responses are also discussed. New strategies based on humans, rodents and other animals (such as zebrafish) are also acknowledged, as they are likely to allow substantial progress to be made in the near future.

Mapping of dwarfing QTL of Ari1327, a semi-dwarf mutant of upland cotton

BACKGROUND

Upland Cotton (Gossypium hirsutum L.) has few cotton varieties suitable for mechanical harvesting. The plant height of the cultivar is one of the key features that need to modify. Hence, this study was planned to locate the QTL for plant height in a 60Co γ treated upland cotton semi-dwarf mutant Ari1327.

RESULTS

Interestingly, bulk segregant analysis (BSA) and genotyping by sequencing (GBS) methods exhibited that candidate QTL was co-located in the region of 5.80-9.66 Mb at D01 chromosome in two F2 populations. Using three InDel markers to genotype a population of 1241 individuals confirmed that the offspring's phenotype is consistent with the genotype. Comparative analysis of RNA-seq between the mutant and wild variety exhibited that Gh_D01G0592 was identified as the source of dwarfness from 200 genes. In addition, it was also revealed that the appropriate use of partial separation markers in QTL mapping can escalate linkage information.

CONCLUSIONS

Overwhelmingly, the results will provide the basis to reveal the function of candidate genes and the utilization of excellent dwarf genetic resources in the future.

Opinion: Is gene mapping in wild populations useful for understanding and predicting adaptation to global change?

Changing environmental conditions will inevitably alter selection pressures. Over the long term, populations have to adapt to these altered conditions by evolutionary change to avoid extinction. Quantifying the 'evolutionary potential' of populations to predict whether they will be able to adapt fast enough to forecasted changes is crucial to fully assess the threat for biodiversity posed by climate change. Technological advances in sequencing and high-throughput genotyping have now made genomic studies possible in a wide range of species. Such studies, in theory, allow an unprecedented understanding of the genomics of ecologically relevant traits and thereby a detailed assessment of the population's evolutionary potential. Aimed at a wider audience than only evolutionary geneticists, this paper gives an overview of how gene-mapping studies have contributed to our understanding and prediction of evolutionary adaptations to climate change, identifies potential reasons why their contribution to understanding adaptation to climate change may remain limited, and highlights approaches to study and predict climate change adaptation that may be more promising, at least in the medium term.

Quantitative trait locus mapping reveals complex genetic architecture of quantitative virulence in the wheat pathogen Zymoseptoria tritici

We conducted a comprehensive analysis of virulence in the fungal wheat pathogen Zymoseptoria tritici using quantitative trait locus (QTL) mapping. High-throughput phenotyping based on automated image analysis allowed the measurement of pathogen virulence on a scale and with a precision that was not previously possible. Across two mapping populations encompassing more than 520 progeny, 540 710 pycnidia were counted and their sizes and grey values were measured. A significant correlation was found between pycnidia size and both spore size and number. Precise measurements of percentage leaf area covered by lesions provided a quantitative measure of host damage. Combining these large and accurate phenotypic datasets with a dense panel of restriction site-associated DNA sequencing (RADseq) genetic markers enabled us to genetically dissect pathogen virulence into components related to host damage and those related to pathogen reproduction. We showed that different components of virulence can be under separate genetic control. Large- and small-effect QTLs were identified for all traits, with some QTLs specific to mapping populations, cultivars and traits and other QTLs shared among traits within the same mapping population. We associated the presence of four accessory chromosomes with small, but significant, increases in several virulence traits, providing the first evidence for a meaningful function associated with accessory chromosomes in this organism. A large-effect QTL involved in host specialization was identified on chromosome 7, leading to the identification of candidate genes having a large effect on virulence.

Transcriptional and Linkage Analyses Identify Loci that Mediate the Differential Macrophage Response to Inflammatory Stimuli and Infection

Macrophages display flexible activation states that range between pro-inflammatory (classical activation) and anti-inflammatory (alternative activation). These macrophage polarization states contribute to a variety of organismal phenotypes such as tissue remodeling and susceptibility to infectious and inflammatory diseases. Several macrophage- or immune-related genes have been shown to modulate infectious and inflammatory disease pathogenesis. However, the potential role that differences in macrophage activation phenotypes play in modulating differences in susceptibility to infectious and inflammatory disease is just emerging. We integrated transcriptional profiling and linkage analyses to determine the genetic basis for the differential murine macrophage response to inflammatory stimuli and to infection with the obligate intracellular parasite Toxoplasma gondii. We show that specific transcriptional programs, defined by distinct genomic loci, modulate macrophage activation phenotypes. In addition, we show that the difference between AJ and C57BL/6J macrophages in controlling Toxoplasma growth after stimulation with interferon gamma and tumor necrosis factor alpha mapped to chromosome 3, proximal to the Guanylate binding protein (Gbp) locus that is known to modulate the murine macrophage response to Toxoplasma. Using an shRNA-knockdown strategy, we show that the transcript levels of an RNA helicase, Ddx1, regulates strain differences in the amount of nitric oxide produced by macrophage after stimulation with interferon gamma and tumor necrosis factor. Our results provide a template for discovering candidate genes that modulate macrophage-mediated complex traits.

Macrophages provide a first line of defense against invading pathogens and play an important role in the initiation and resolution of immune responses. When in contact with pathogens or immune factors, such as cytokines, macrophages assume activation states that range between pro-inflammatory (classical activation) and anti-inflammatory (alternative activation). Even though it is known that macrophages from different individuals are biased towards one of the various activation states, the genetic factors that define individual differences in macrophage activation are not fully understood. Additionally, although macrophages are important in infectious disease pathogenesis, how individual differences in macrophage activation contribute to individual differences in susceptibility to infectious disease is just emerging. We used macrophages from genetically segregating mice to show that discrete transcriptional programs, which are modulated by specific genomic regions, modulate differences in macrophage activation. Murine macrophages differences in controlling Toxoplasma growth mapped to chromosome 3, proximal to the Guanylate binding protein (Gbp) locus that is known to modulate the murine macrophage response to Toxoplasma. Using a shRNA-mediated knockdown approach, we show that the DEAD box polypeptide 1 (Ddx1) modulates nitric oxide production in macrophages stimulated with interferon gamma and tumor necrosis factor. These findings are a step towards the identification of genes that regulate macrophage phenotypes and disease outcome.

Discrimination of korean native chicken lines using fifteen selected microsatellite markers

In order to evaluate the genetic diversity and discrimination among five Korean native chicken lines, a total of 86 individuals were genotyped using 150 microsatellite (MS) markers, and 15 highly polymorphic MS markers were selected. Based on the highest value of the number of alleles, the expected heterozygosity (He) and polymorphic information content (PIC) for the selected markers ranged from 6 to 12, 0.466 to 0.852, 0.709 to 0.882 and 0.648 to 0.865, respectively. Using these markers, the calculated genetic distance (Fst), the heterozygote deficit among chicken lines (Fit) and the heterozygote deficit within chicken line (Fis) values ranged from 0.0309 to 0.2473, 0.0013 to 0.4513 and -0.1002 to 0.271, respectively. The expected probability of identity values in random individuals (PI), random half-sib (PI half-sibs ) and random sibs (PI sibs ) were estimated at 7.98×10(-29), 2.88×10(-20) and 1.25×10(-08), respectively, indicating that these markers can be used for traceability systems in Korean native chickens. The unrooted phylogenetic neighbor-joining (NJ) tree was constructed using 15 MS markers that clearly differentiated among the five native chicken lines. Also, the structure was estimated by the individual clustering with the K value of 5. The selected 15 MS markers were found to be useful for the conservation, breeding plan, and traceability system in Korean native chickens.

Expression quantitative trait analyses to identify causal genetic variants for type 2 diabetes susceptibility

Type 2 diabetes (T2D) is a common metabolic disorder which is caused by multiple genetic perturbations affecting different biological pathways. Identifying genetic factors modulating the susceptibility of this complex heterogeneous metabolic phenotype in different ethnic and racial groups remains challenging. Despite recent success, the functional role of the T2D susceptibility variants implicated by genome-wide association studies (GWAS) remains largely unknown. Genetic dissection of transcript abundance or expression quantitative trait (eQTL) analysis unravels the genomic architecture of regulatory variants. Availability of eQTL information from tissues relevant for glucose homeostasis in humans opens a new avenue to prioritize GWAS-implicated variants that may be involved in triggering a causal chain of events leading to T2D. In this article, we review the progress made in the field of eQTL research and knowledge gained from those studies in understanding transcription regulatory mechanisms in human subjects. We highlight several novel approaches that can integrate eQTL analysis with multiple layers of biological information to identify ethnic-specific causal variants and gene-environment interactions relevant to T2D pathogenesis. Finally, we discuss how the eQTL analysis mediated search for “missing heritability” may lead us to novel biological and molecular mechanisms involved in susceptibility to T2D.

Reading and language disorders: the importance of both quantity and quality

Reading and language disorders are common childhood conditions that often co-occur with each other and with other neurodevelopmental impairments. There is strong evidence that disorders, such as dyslexia and Specific Language Impairment (SLI), have a genetic basis, but we expect the contributing genetic factors to be complex in nature. To date, only a few genes have been implicated in these traits. Their functional characterization has provided novel insight into the biology of neurodevelopmental disorders. However, the lack of biological markers and clear diagnostic criteria have prevented the collection of the large sample sizes required for well-powered genome-wide screens. One of the main challenges of the field will be to combine careful clinical assessment with high throughput genetic technologies within multidisciplinary collaborations.

The endocrine stress response is linked to one specific locus on chromosome 3 in a mouse model based on extremes in trait anxiety

Background

The hypothalamic-pituitary-adrenal (HPA) axis is essential to control physiological stress responses in mammals. Its dysfunction is related to several mental disorders, including anxiety and depression. The aim of this study was to identify genetic loci underlying the endocrine regulation of the HPA axis.

Method

High (HAB) and low (LAB) anxiety-related behaviour mice were established by selective inbreeding of outbred CD-1 mice to model extremes in trait anxiety. Additionally, HAB vs. LAB mice exhibit comorbid characteristics including a differential corticosterone response upon stress exposure. We crossbred HAB and LAB lines to create F1 and F2 offspring. To identify the contribution of the endocrine phenotypes to the total phenotypic variance, we examined multiple behavioural paradigms together with corticosterone secretion-based phenotypes in F2 mice by principal component analysis. Further, to pinpoint the genomic loci of the quantitative trait of the HPA axis stress response, we conducted genome-wide multipoint oligogenic linkage analyses based on Bayesian Markov chain Monte Carlo approach as well as parametric linkage in three-generation pedigrees, followed by a two-dimensional scan for epistasis and association analysis in freely segregating F2 mice using 267 single-nucleotide polymorphisms (SNPs), which were identified to consistently differ between HAB and LAB mice as genetic markers.

Results

HPA axis reactivity measurements and behavioural phenotypes were represented by independent principal components and demonstrated no correlation. Based on this finding, we identified one single quantitative trait locus (QTL) on chromosome 3 showing a very strong evidence for linkage (2ln (L-score) > 10, LOD > 23) and significant association (lowest Bonferroni adjusted p < 10^-28) to the neuroendocrine stress response. The location of the linkage peak was estimated at 42.3 cM (95% confidence interval: 41.3 - 43.3 cM) and was shown to be in epistasis (p-adjusted < 0.004) with the locus at 35.3 cM on the same chromosome. The QTL harbours genes involved in steroid synthesis and cardiovascular effects.

Conclusion

The very prominent effect on stress-induced corticosterone secretion of the genomic locus on chromosome 3 and its involvement in epistasis highlights the critical role of this specific locus in the regulation of the HPA axis.

Genetic determinants of metabolism in health and disease: from biochemical genetics to genome-wide associations

Increasingly sophisticated measurement technologies have allowed the fields of metabolomics and genomics to identify, in parallel, risk factors of disease; predict drug metabolism; and study metabolic and genetic diversity in large human populations. Yet the complementarity of these fields and the utility of studying genes and metabolites together is belied by the frequent separate, parallel applications of genomic and metabolomic analysis. Early attempts at identifying co-variation and interaction between genetic variants and downstream metabolic changes, including metabolic profiling of human Mendelian diseases and quantitative trait locus mapping of individual metabolite concentrations, have recently been extended by new experimental designs that search for a large number of gene-metabolite associations. These approaches, including metabolomic quantitiative trait locus mapping and metabolomic genome-wide association studies, involve the concurrent collection of both genomic and metabolomic data and a subsequent search for statistical associations between genetic polymorphisms and metabolite concentrations across a broad range of genes and metabolites. These new data-fusion techniques will have important consequences in functional genomics, microbial metagenomics and disease modeling, the early results and implications of which are reviewed.

Design considerations for genetic linkage and association studies

This chapter describes the main issues that genetic epidemiologists usually consider in the design of linkage and association studies. For linkage, we briefly consider the situation of rare, highly penetrant alleles showing a disease pattern consistent with Mendelian inheritance investigated through parametric methods in large pedigrees or with autozygosity mapping in inbred families, and we then turn our focus to the most common design, affected sibling pairs, of more relevance for common, complex diseases. Theoretical and more practical power and sample size calculations are provided as a function of the strength of the genetic effect being investigated. We also discuss the impact of other determinants of statistical power such as disease heterogeneity, pedigree, and genotyping errors, as well as the effect of the type and density of genetic markers. Linkage studies should be as large as possible to have sufficient power in relation to the expected genetic effect size. Segregation analysis, a formal statistical technique to describe the underlying genetic susceptibility, may assist in the estimation of the relevant parameters to apply, for instance. However, segregation analyses estimate the total genetic component rather than a single-locus effect. Locus heterogeneity should be considered when power is estimated and at the analysis stage, i.e. assuming smaller locus effect than the total the genetic component from segregation studies. Disease heterogeneity should be minimised by considering subtypes if they are well defined or by otherwise collecting known sources of heterogeneity and adjusting for them as covariates; the power will depend upon the relationship between the disease subtype and the underlying genotypes. Ultimately, identifying susceptibility alleles of modest effects (e.g. RR≤1.5) requires a number of families that seem unfeasible in a single study. Meta-analysis and data pooling between different research groups can provide a sizeable study, but both approaches require even a higher level of vigilance about locus and disease heterogeneity when data come from different populations. All necessary steps should be taken to minimise pedigree and genotyping errors at the study design stage as they are, for the most part, due to human factors. A two-stage design is more cost-effective than one stage when using short tandem repeats (STRs). However, dense single-nucleotide polymorphism (SNP) arrays offer a more robust alternative, and due to their lower cost per unit, the total cost of studies using SNPs may in the future become comparable to that of studies using STRs in one or two stages. For association studies, we consider the popular case-control design for dichotomous phenotypes, and we provide power and sample size calculations for one-stage and multistage designs. For candidate genes, guidelines are given on the prioritisation of genetic variants, and for genome-wide association studies (GWAS), the issue of choosing an appropriate SNP array is discussed. A warning is issued regarding the danger of designing an underpowered replication study following an initial GWAS. The risk of finding spurious association due to population stratification, cryptic relatedness, and differential bias is underlined. GWAS have a high power to detect common variants of high or moderate effect. For weaker effects (e.g. relative risk<1.2), the power is greatly reduced, particularly for recessive loci. While sample sizes of 10,000 or 20,000 cases are not beyond reach for most common diseases, only meta-analyses and data pooling can allow attaining a study size of this magnitude for many other diseases. It is acknowledged that detecting the effects from rare alleles (i.e. frequency<5%) is not feasible in GWAS, and it is expected that novel methods and technology, such as next-generation resequencing, will fill this gap. At the current stage, the choice of which GWAS SNP array to use does not influence the power in populations of European ancestry. A multistage design reduces the study cost but has less power than the standard one-stage design. If one opts for a multistage design, the power can be improved by jointly analysing the data from different stages for the SNPs they share. The estimates of locus contribution to disease risk from genome-wide scans are often biased, and relying on them might result in an underpowered replication study. Population structure has so far caused less spurious associations than initially feared, thanks to systematic ethnicity matching and application of standard quality control measures. Differential bias could be a more serious threat and must be minimised by strictly controlling all the aspects of DNA acquisition, storage, and processing.

Linkage analysis in the next-generation sequencing era

Linkage analysis was developed to detect excess co-segregation of the putative alleles underlying a phenotype with the alleles at a marker locus in family data. Many different variations of this analysis and corresponding study design have been developed to detect this co-segregation. Linkage studies have been shown to have high power to detect loci that have alleles (or variants) with a large effect size, i.e. alleles that make large contributions to the risk of a disease or to the variation of a quantitative trait. However, alleles with a large effect size tend to be rare in the population. In contrast, association studies are designed to have high power to detect common alleles which tend to have a small effect size for most diseases or traits. Although genome-wide association studies have been successful in detecting many new loci with common alleles of small effect for many complex traits, these common variants often do not explain a large proportion of disease risk or variation of the trait. In the past, linkage studies were successful in detecting regions of the genome that were likely to harbor rare variants with large effect for many simple Mendelian diseases and for many complex traits. However, identifying the actual sequence variant(s) responsible for these linkage signals was challenging because of difficulties in sequencing the large regions implicated by each linkage peak. Current 'next-generation' DNA sequencing techniques have made it economically feasible to sequence all exons or the whole genomes of a reasonably large number of individuals. Studies have shown that rare variants are quite common in the general population, and it is now possible to combine these new DNA sequencing methods with linkage studies to identify rare causal variants with a large effect size. A brief review of linkage methods is presented here with examples of their relevance and usefulness for the interpretation of whole-exome and whole-genome sequence data.

Estimation and visualization of identity-by-descent within pedigrees simplifies interpretation of complex trait analysis

Linkage analysis identifies markers that appear to be co-inherited with a trait within pedigrees. The inheritance of a chromosomal segment may be probabilistically reconstructed, with missing data complicating inference. Inheritance patterns are further obscured in the analysis of complex traits, where variants in one or more genes may contribute to phenotypic variation within a pedigree. In this case, determining which relatives share a trait variant is not simple. We describe how to represent these patterns of inheritance for marker loci. We summarize how to sample patterns of inheritance consistent with genotypic and pedigree data using gl_auto, available in MORGAN v3.0. We describe identification of classes of equivalent inheritance patterns with the program IBDgraph. We finally provide an example of how these programs may be used to simplify interpretation of linkage analysis of complex traits in general pedigrees. We borrow information across loci in a parametric linkage analysis of a large pedigree. We explore the contribution of each equivalence class to a linkage signal, illustrate estimated patterns of identity-by-descent sharing, and identify a haplotype tagging the chromosomal segment driving the linkage signal. Haplotype carriers are more likely to share the linked trait variant, and can be prioritized for subsequent DNA sequencing.

The Maximum-Likelihood-Binomial method revisited: a robust approach for model-free linkage analysis of quantitative traits in large sibships

Model-free linkage analysis methods, based on identity-by-descent allele sharing, are commonly used for complex trait analysis. The Maximum-Likelihood-Binomial (MLB) approach, which is based on the hypothesis that parental alleles are binomially distributed among affected sibs, is particularly popular. An extension of this method to quantitative traits (QT) has been proposed (MLB-QTL), based on the introduction of a latent binary variable capturing information about the linkage between the QT and the marker. Interestingly, the MLB-QTL method does not require the decomposition of sibships into constituent sibpairs and requires no prior assumption about the distribution of the QT. We propose a new formulation of the MLB method for quantitative traits (nMLB-QTL) that explicitly takes advantage of the independence of paternal and maternal allele transmission under the null hypothesis of no linkage. Simulation studies under H₀ showed that the nMLB-QTL method generated very consistent type I errors. Furthermore, simulations under the alternative hypothesis showed that the nMLB-QTL method was slightly, but systematically more powerful than the MLB-QTL method, whatever the genetic model, residual correlation, ascertainment strategy and sibship size considered. Finally, the power of the nMLB-QTL method is illustrated by a chromosome-wide linkage scan for a quantitative endophenotype of leprosy infection. Overall, the nMLB-QTL method is a robust, powerful, and flexible approach for detecting linkage with quantitative phenotypes, particularly in studies of non Gaussian phenotypes in large sibships.

Fluoride's effects on the formation of teeth and bones, and the influence of genetics

Fluorides are present in the environment. Excessive systemic exposure to fluorides can lead to disturbances of bone homeostasis (skeletal fluorosis) and enamel development (dental/enamel fluorosis). The severity of dental fluorosis is also dependent upon fluoride dose and the timing and duration of fluoride exposure. Fluoride's actions on bone cells predominate as anabolic effects both in vitro and in vivo. More recently, fluoride has been shown to induce osteoclastogenesis in mice. Fluorides appear to mediate their actions through the MAPK signaling pathway and can lead to changes in gene expression, cell stress, and cell death. Different strains of inbred mice demonstrate differential physiological responses to ingested fluoride. Genetic studies in mice are capable of identifying and characterizing fluoride-responsive genetic variations. Ultimately, this can lead to the identification of at-risk human populations who are susceptible to the unwanted or potentially adverse effects of fluoride action and to the elucidation of fundamental mechanisms by which fluoride affects biomineralization.

Mapping of quantitative trait loci for mycoplasma and tetanus antibodies and interferon-gamma in a porcine F(2) Duroc x Pietrain resource population

The aim of the present study was to detect quantitative trait loci (QTL) for innate and adaptive immunity in pigs. For this purpose, a Duroc x Pietrain F(2) resource population (DUPI) with 319 offspring was used to map QTL for the immune traits blood antibodies and interferon-gamma using 122 microsatellites covering all autosomes. Antibodies response to Mycoplasma hyopneumoniae and tetanus toxoid vaccine and the interferon-gamma (IFNG) serum concentration were measured at three different time points and were used as phenotypes. The differences of antibodies and interferon concentration between different time points were also used for the linkage mapping. Line-cross and imprinting QTL analysis, including two-QTL, were performed using QTL Express. A total of 30 QTL (12, 6, and 12 for mycoplasma, tetanus antibody, and IFNG, respectively) were identified at the 5% chromosome-wide-level significant, of which 28 were detected by line-cross and 2 by imprinting model. In addition, two QTL were identified on chromosome 5 using the two-QTL approach where both loci were in repulsion phase. Most QTL were detected on pig chromosomes 2, 5, 11, and 18. Antibodies were increased over time and immune traits were found to be affected by sex, litter size, parity, and month of birth. The results demonstrated that antibody and IFNG concentration are influenced by multiple chromosomal areas. The flanking markers of the QTL identified for IFNG on SSC5 did incorporate the position of the porcine IFNG gene. The detected QTL will allow further research in these QTL regions for candidate genes and their utilization in selection to improve the immune response and disease resistance in pig.

Interspecific resources: a major tool for quantitative trait locus cloning and speciation research

Positional cloning of the quantitative trait locus (QTL) still encounters numerous difficulties, which explains why thousands of QTL have been mapped, while only a few have been identified at the molecular level. Here, we focus on a specific mapping tool that exists in plant and animal model species: interspecific recombinant congenic strains (IRCSs) or interspecific nearly isogenic lines (NILs). Such panels exhibit a much higher sequence diversity than intraspecific sets, thus enhancing the contrasts between phenotypes. In animals, it allows statistical significance to be reached even when using a limited number of individuals. Therefore, we argue that interspecific resources may constitute a major genetic tool for positional cloning and for understanding some bases of speciation mechanisms.

Identifying quantitative trait loci (QTLs) and genes (QTGs) for alcohol-related phenotypes in mice

Alcoholism is a complex clinical disorder with genetic and environmental contributions. Although no animal model duplicates alcoholism, models for specific factors, such as the withdrawal syndrome, are useful to identify potential genetic determinants of liability in humans. Murine models have been invaluable to identify quantitative trait loci (QTLs) that influence a variety of alcohol responses. However, the QTL regions are typically large, at least initially, and contain numerous genes, making identification of the causal quantitative trait gene(s) (QTGs) challenging. Here, we present QTG identification strategies currently used in the field of alcohol genetics and discuss relevance to alcoholic human populations.