Genomewide association analysis in diverse inbred mice: power and population structure

The discovery of quantitative trait loci (QTL) in model organisms has relied heavily on the ability to perform controlled breeding to generate genotypic and phenotypic diversity. Recently, we and others have demonstrated the use of an existing set of diverse inbred mice (referred to here as the mouse diversity panel, MDP) as a QTL mapping population. The use of the MDP population has many advantages relative to traditional F(2) mapping populations, including increased phenotypic diversity, a higher recombination frequency, and the ability to collect genotype and phenotype data in community databases. However, these methods are complicated by population structure inherent in the MDP and the lack of an analytical framework to assess statistical power. To address these issues, we measured gene expression levels in hypothalamus across the MDP. We then mapped these phenotypes as quantitative traits with our association algorithm, resulting in a large set of expression QTL (eQTL). We utilized these eQTL, and specifically cis-eQTL, to develop a novel nonparametric method for association analysis in structured populations like the MDP. These eQTL data confirmed that the MDP is a suitable mapping population for QTL discovery and that eQTL results can serve as a gold standard for relative measures of statistical power.     

关联分析及其在植物遗传学研究中的应用

植物的很多重要经济性状均属于复杂性状。基于连锁分析的QTL作图是研究复杂性状的有效手段, 但其尚存在一定的局限性。随着现代生物学的发展, 一种基于连锁不平衡的新剖分复杂性状方法--关联分析法, 开始应用于植物遗传学研究。与QTL作图法相比, 应用关联分析法具有不需要构建特殊的群体, 可同时对多个等位基因进行分析, 定位QTL精度可达到单基因水平等优势。该文介绍了关联分析方法学的基础和特性, 简述了其在植物遗传学研究中的进展情况, 并对其未来发展和在植物遗传学研究中的应用进行了展望。     

Segregation of regulatory polymorphisms with effects on the gluteus medius transcriptome in a purebred pig population

Background

The main goal of the present study was to analyse the genetic architecture of mRNA expression in muscle, a tissue with an outmost economic importance for pig breeders. Previous studies have used F₂ crosses to detect porcine expression QTL (eQTL), so they contributed with data that mostly represents the between-breed component of eQTL variation. Herewith, we have analysed eQTL segregation in an outbred Duroc population using two groups of animals with divergent fatness profiles. This approach is particularly suitable to analyse the within-breed component of eQTL variation, with a special emphasis on loci involved in lipid metabolism.

Methodology/Principal Findings

GeneChip Porcine Genome arrays (Affymetrix) were used to determine the mRNA expression levels of gluteus medius samples from 105 Duroc barrows. A whole-genome eQTL scan was carried out with a panel of 116 microsatellites. Results allowed us to detect 613 genome-wide significant eQTL unevenly distributed across the pig genome. A clear predominance of trans- over cis-eQTL, was observed. Moreover, 11 trans-regulatory hotspots affecting the expression levels of four to 16 genes were identified. A Gene Ontology study showed that regulatory polymorphisms affected the expression of muscle development and lipid metabolism genes. A number of positional concordances between eQTL and lipid trait QTL were also found, whereas limited evidence of a linear relationship between muscle fat deposition and mRNA levels of eQTL regulated genes was obtained.

Conclusions/Significance

Our data provide substantial evidence that there is a remarkable amount of within-breed genetic variation affecting muscle mRNA expression. Most of this variation acts in trans and influences biological processes related with muscle development, lipid deposition and energy balance. The identification of the underlying causal mutations and the ascertainment of their effects on phenotypes would allow gaining a fundamental perspective about how complex traits are built at the molecular level.     

遗传与基因表达数据的整合—— eQTL的方法及应用

高通量的基因型分析和芯片技术的发展使人们能够进一步研究哪些遗传差异最终影响基因的表达。通过表达数量性状座位(eQTL)作图方法可对基因表达水平的遗传基础进行解析。与传统的QTL分析方法一样, eQTL的主要目标是鉴别表达性状座位所在的染色体区域。但由于表达谱数据成千上万, 而传统的QTL分析方法最多分析几十个性状, 因此需要考虑这类实验设计的特点以及统计分析方法。本文详细介绍了eQTL定位过程及其研究方法, 重点从个体选择、基因芯片实验设计、基因表达数据的获得与标准化、作图方法及结果分析等方面进行了综述, 指出了当前eQTL研究存在的问题和局限性。最后介绍了eQTL研究在估计基因表达遗传率、挖掘候选基因、构建基因调控网络、理解基因间及基因与环境的互作的应用进展。     

A novel strategy for genetic dissection of complex traits: the population of specific chromosome substitution strains from laboratory and wild mice

The mouse is an irreplaceable model for understanding the precise genetic mechanisms of mammalian physiological pathways. Thousands of quantitative trait loci (QTLs) have been mapped onto the mouse genome during the last two decades. However, only a few genes’ underlying complex traits have been successfully identified, and rapid fine mapping of QTL genes still remains a challenge for mouse geneticists. Currently, the Collaborative Cross (CC) has proceeded to the goal of establishing more than 1,000 recombinant inbred strains for the sub-centimorgan mapping resolution of QTL loci. In this article, a novel complementary strategy, designated as population of specific chromosome substitution strains or PSCSS, is proposed for rapid fine mapping of QTLs on the substituted chromosome. One specific chromosome (Chr 1) of recipient mouse strain C57BL/6 J has been substituted by homologous counterparts from five different inbred strains (C3H/He, FVB/N, AKR, NOD/LtJ, NZW/LacJ), an outbred line Kunmin mouse in China, and 23 wild mice captured in different localities. The primary genetic studies on the structure of these wild donor chromosomes (Chr 1) show that these donor chromosomes harbor extensive genetic polymorphisms, with a high density of SNP distribution, abundant variations of STR alleles, and a high level of historical recombination accumulation. These specific chromosome substitution strains eventually form a special population that has the identical genetic background of the recipient strain and differs only in the donor chromosomes. With simple association studies, known QTLs on the donor chromosome can be rapidly mapped in high resolution without requirement of further crosses. This approach, taking advantage of the extensive genetic polymorphisms of wild resources and chromosome substitution strategy, brings a new outlook for genetic dissection of complex traits.     

A random model approach to mapping quantitative trait loci for complex binary traits in outbred populations.

Mapping quantitative trait loci (QTL) for complex binary traits is more challenging than for normally distributed traits due to the nonlinear relationship between the observed phenotype and unobservable genetic effects, especially when the mapping population contains multiple outbred families. Because the number of alleles of a QTL depends on the number of founders in an outbred population, it is more appropriate to treat the effect of each allele as a random variable so that a single variance rather than individual allelic effects is estimated and tested. Such a method is called the random model approach. In this study, we develop the random model approach of QTL mapping for binary traits in outbred populations. An EM-algorithm with a Fisher-scoring algorithm embedded in each E-step is adopted here to estimate the genetic variances. A simple Monte Carlo integration technique is used here to calculate the likelihood-ratio test statistic. For the first time we show that QTL of complex binary traits in an outbred population can be scanned along a chromosome for their positions, estimated for their explained variances, and tested for their statistical significance. Application of the method is illustrated using a set of simulated data.     

Gene expression differences in mice divergently selected for methamphetamine sensitivity

In an effort to identify genes that may be important for drug-abuse liability, we mapped behavioral quantitative trait loci (bQTL) for sensitivity to the locomotor stimulant effect of methamphetamine (MA) using two mouse lines that were selectively bred for high MA-induced activity (HMACT) or low MA-induced activity (LMACT). We then examined gene expression differences between these lines in the nucleus accumbens, using 20 U74Av2 Affymetrix microarrays and quantitative polymerase chain reaction (qPCR). Expression differences were detected for several genes, including Casein Kinase 1 Epsilon (Csnkle), glutamate receptor, ionotropic, AMPA1 (GluR1), GABA B1 receptor (Gabbr1), and dopamine- and cAMP-regulated phosphoprotein of 32 kDa (Darpp-32). We used the database to identify QTL that regulate the expression of the genes identified by the microarrays (expression QTL; eQTL). This approach identified an eQTL for Csnkle on Chromosome 15 (LOD=3.8) that comapped with a bQTL for the MA stimulation phenotype (LOD=4.5), suggesting that a single allele may cause both traits. The chromosomal region containing this QTL has previously been associated with sensitivity to the stimulant effects of cocaine. These results suggest that selection was associated with (and likely caused) altered gene expression that is partially attributable to different frequencies of gene expression polymorphisms. Combining classical genetics with analysis of whole-genome gene expression and bioinformatic resources provides a powerful method for provisionally identifying genes that influence complex traits. The identified genes provide excellent candidates for future hypothesis-driven studies, translational genetic studies, and pharmacological interventions.     

Identification of candidate genes in the type 2 diabetes modifier locus using expression QTL

To identify new genetic determinants relevant to type 2 diabetes (T2D), diabetic F2 progeny were generated by intercrossing F1 mice obtained from a cross of BKS.Cg-Lepr(db)+/+m and DBA/2, and T2D-related phenotypes were measured. In the F2 population, increased susceptibility to diabetes and obesity was observed. We also detected the major quantitative trait loci (QTL) modifying the severity of diabetes on chromosome 9, where peaks of logarithm of odds (LOD) overlapped for three traits. To identify candidate genes in the QTL intervals, we combined "expression QTL" (eQTL), taking mRNA levels as quantitative traits, and "interstrain sequence variations, including cSNPs." As a result, four genes were identified from cosegregation of clinical QTL with eQTL and 13 genes were found from interstrain cSNPs as candidates in the T2D modifier QTL. Our combined approach shows the acceleration of the discovery of candidate genes in the QTL of interest, spanning several megabases.     

Impact of Natural Genetic Variation on Gene Expression Dynamics

DNA sequence variation causes changes in gene expression, which in turn has profound effects on cellular states. These variations affect tissue development and may ultimately lead to pathological phenotypes. A genetic locus containing a sequence variation that affects gene expression is called an “expression quantitative trait locus” (eQTL). Whereas the impact of cellular context on expression levels in general is well established, a lot less is known about the cell-state specificity of eQTL. Previous studies differed with respect to how “dynamic eQTL” were defined. Here, we propose a unified framework distinguishing static, conditional and dynamic eQTL and suggest strategies for mapping these eQTL classes. Further, we introduce a new approach to simultaneously infer eQTL from different cell types. By using murine mRNA expression data from four stages of hematopoiesis and 14 related cellular traits, we demonstrate that static, conditional and dynamic eQTL, although derived from the same expression data, represent functionally distinct types of eQTL. While static eQTL affect generic cellular processes, non-static eQTL are more often involved in hematopoiesis and immune response. Our analysis revealed substantial effects of individual genetic variation on cell type-specific expression regulation. Among a total number of 3,941 eQTL we detected 2,729 static eQTL, 1,187 eQTL were conditionally active in one or several cell types, and 70 eQTL affected expression changes during cell type transitions. We also found evidence for feedback control mechanisms reverting the effect of an eQTL specifically in certain cell types. Loci correlated with hematological traits were enriched for conditional eQTL, thus, demonstrating the importance of conditional eQTL for understanding molecular mechanisms underlying physiological trait variation. The classification proposed here has the potential to streamline and unify future analysis of conditional and dynamic eQTL as well as many other kinds of QTL data.     

Molecular mechanisms of fiber differential development between G. barbadense and G. hirsutum revealed by genetical genomics

Cotton fiber qualities including length, strength and fineness are known to be controlled by genes affecting cell elongation and secondary cell wall (SCW) biosynthesis, but the molecular mechanisms that govern development of fiber traits are largely unknown. Here, we evaluated an interspecific backcrossed population from G. barbadense cv. Hai7124 and G. hirsutum acc. TM-1 for fiber characteristics in four-year environments under field conditions, and detected 12 quantitative trait loci (QTL) and QTL-by-environment interactions by multi-QTL joint analysis. Further analysis of fiber growth and gene expression between TM-1 and Hai7124 showed greater differences at 10 and 25 days post-anthesis (DPA). In this two period important for fiber performances, we integrated genome-wide expression profiling with linkage analysis using the same genetic materials and identified in total 916 expression QTL (eQTL) significantly (P<0.05) affecting the expression of 394 differential genes. Many positional cis-/trans-acting eQTL and eQTL hotspots were detected across the genome. By comparative mapping of eQTL and fiber QTL, a dataset of candidate genes affecting fiber qualities was generated. Real-time quantitative RT-PCR (qRT-PCR) analysis confirmed the major differential genes regulating fiber cell elongation or SCW synthesis. These data collectively support molecular mechanism for G. hirsutum and G. barbadense through differential gene regulation causing difference of fiber qualities. The down-regulated expression of abscisic acid (ABA) and ethylene signaling pathway genes and high-level and long-term expression of positive regulators including auxin and cell wall enzyme genes for fiber cell elongation at the fiber developmental transition stage may account for superior fiber qualities.