A genetical genomics approach to genome scans increases power for QTL mapping

We describe a method for integrating gene expression information into genome scans and show that this can substantially increase the statistical power of QTL mapping. The method has three stages. First, standard clustering methods identify small (size 5-20) groups of genes with similar expression patterns. Second, each gene group is tested for a causative genetic locus shared with the clinical trait of interest. This is done using an EM algorithm approach that treats genotype at the putative causative locus as an unobserved variable and combines expression information from all of the genes in the group to infer genotype information at the locus. Finally, expression QTL (eQTL) are mapped for each gene group that shares a causative locus with the clinical trait. Such eQTL are candidates for the causative locus. Simulation results show that this method has far superior power to standard QTL mapping techniques in many circumstances. We applied this method to existing data on mouse obesity. Our method identified 27 putative body weight QTL, whereas standard QTL mapping produced only one. Furthermore, most gene groups with body weight QTL included cis genes, so candidate genes could be immediately identified. Eleven body weight QTL produced 16 candidate genes that have been previously associated with body weight or body weight-related traits, thus validating our method. In addition, 15 of the 16 other loci produced 32 candidate genes that have not been associated with body weight. Thus, this method shows great promise for finding new causative loci for complex traits.     

遗传基因组学(Genetical genomics)的研究进展

遗传基因组学(geneticalgenomics)是将微阵列技术和数量性状座位(QTL)分析结合起来,全基因组水平上定位基因表达的QTL(eQTL).它为研究复杂性状的分子机理和调控网络提供全新的手段.遗传基因组这个概念和研究策略在2001年由Janson和Nap首先提出,到目前为止,遗传基因组学已应用于酵母、老鼠、人以及玉米等植物.研究结果表明:基因表达水平的差异是可遗传的复杂性状;eQTL可以分为顺式作用eQTL和反式作用eQTL,顺式作用eQTL就是某个基因的eQTL定位到该基因所在的基因组区域,表明可能是该基因本身的差别引起mRNA水平的差别,反式作用就是eQTL定位到其他基因组区域,表明其他基因的差别控制该基因mRNA水平的差异.将eQTL结果、基因功能注解以及多种统计分析方法相结合,不仅能更准确地鉴别控制复杂性状及其相关基因表达的候选基因,而且能构建相应的基因调控网络.     

Architecture of energy balance traits in emerging lines of the Collaborative Cross

The potential utility of the Collaborative Cross (CC) mouse resource was evaluated to better understand complex traits related to energy balance. A primary focus was to examine if genetic diversity in emerging CC lines (pre-CC) would translate into equivalent phenotypic diversity. Second, we mapped quantitative trait loci (QTL) for 15 metabolism- and exercise-related phenotypes in this population. We evaluated metabolic and voluntary exercise traits in 176 pre-CC lines, revealing phenotypic variation often exceeding that seen across the eight founder strains from which the pre-CC was derived. Many phenotypic correlations existing within the founder strains were no longer significant in the pre-CC population, potentially representing reduced linkage disequilibrium (LD) of regions harboring multiple genes with effects on energy balance or disruption of genetic structure of extant inbred strains with substantial shared ancestry. QTL mapping revealed five significant and eight suggestive QTL for body weight (Chr 4, 7.54 Mb; CI 3.32-10.34 Mb; Bwq14), body composition, wheel running (Chr 16, 33.2 Mb; CI 32.5-38.3 Mb), body weight change in response to exercise (1: Chr 6, 77.7Mb; CI 72.2-83.4 Mb and 2: Chr 6, 42.8 Mb; CI 39.4-48.1 Mb), and food intake during exercise (Chr 12, 85.1 Mb; CI 82.9-89.0 Mb). Some QTL overlapped with previously mapped QTL for similar traits, whereas other QTL appear to represent novel loci. These results suggest that the CC will be a powerful, high-precision tool for examining the genetic architecture of complex traits such as those involved in regulation of energy balance.     

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.     

Genetic Regulation of Bone Metabolism in the Chicken: Similarities and Differences to Mammalian Systems

Birds have a unique bone physiology, due to the demands placed on them through egg production. In particular their medullary bone serves as a source of calcium for eggshell production during lay and undergoes continuous and rapid remodelling. We take advantage of the fact that bone traits have diverged massively during chicken domestication to map the genetic basis of bone metabolism in the chicken. We performed a quantitative trait locus (QTL) and expression QTL (eQTL) mapping study in an advanced intercross based on Red Junglefowl (the wild progenitor of the modern domestic chicken) and White Leghorn chickens. We measured femoral bone traits in 456 chickens by peripheral computerised tomography and femoral gene expression in a subset of 125 females from the cross with microarrays. This resulted in 25 loci for female bone traits, 26 loci for male bone traits and 6318 local eQTL loci. We then overlapped bone and gene expression loci, before checking for an association between gene expression and trait values to identify candidate quantitative trait genes for bone traits. A handful of our candidates have been previously associated with bone traits in mice, but our results also implicate unexpected and largely unknown genes in bone metabolism. In summary, by utilising the unique bone metabolism of an avian species, we have identified a number of candidate genes affecting bone allocation and metabolism. These findings can have ramifications not only for the understanding of bone metabolism genetics in general, but could also be used as a potential model for osteoporosis as well as revealing new aspects of vertebrate bone regulation or features that distinguish avian and mammalian bone.     

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.     

A review of statistical methods for expression quantitative trait loci mapping

With high-throughput technologies now widely available, investigators can easily measure thousands of phenotypes for quantitative trait loci (QTL) mapping. Microarray measurements are particularly amenable to QTL mapping, as evidenced by a number of recent studies demonstrating utility across a broad range of biological endeavors. The early success stories have impelled a rapid increase in both the number and complexity of expression QTL (eQTL) experiments. Consequently, there is a need to consider the statistical principles involved in the design and analysis of these experiments and the methods currently being used. In this article we review these principles and methods and discuss the open questions most likely to yield significant progress toward increasing the amount of meaningful information obtained from eQTL mapping experiments.     

Quantitative trait loci for chemical body composition traits in pigs and their positional associations with body tissues, growth and feed intake

In this study, quantitative trait loci (QTL) for chemical and physical body composition, growth and feed intake in pigs were identified in a three-generation full-sib population, developed by crossing Pietrain sires with a commercial dam line. Phenotypic data from 315 F₂ animals were available for protein and lipid deposition measured in live animals by the deuterium dilution technique at 30-, 60-, 90-, 120- and 140-kg body weight. At 140-kg body weight, carcass characteristics were measured by the AutoFOM grading system and after dissection. Three hundred and eighty-six animals from 49 families were genotyped for 51 molecular markers covering chromosomes SSC2, SSC4, SSC8, SSC9, SSC10 and SSC14. Novel QTL for protein (lipid) content at 60-kg body weight and protein (lipid) accretion from 120 to 140 kg were detected on SSC9 near several previously detected QTL for lean and fat tissue in neck, shoulder and ham cuts. Another QTL for lipid accretion was found on SSC8, closely associated with a QTL for intramuscular fat content. QTL for daily feed intake were detected on SSC2 and SSC10. The favourable allele of a QTL for food conversion ratio (FCR) on SSC2 was associated with alleles for increased lean tissue and decreased fat tissue. Because no QTL for growth rate were found in the region, the QTL for FCR is most likely due to a change in body composition. These QTL provide insights into the genomic regulation of chemical or physical body composition and its association with feed intake, feed efficiency and growth.     

Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast- and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.     

Quantitative trait loci for body growth and sex determination in the hermaphrodite teleost fish Sparus aurata L.

Gilthead sea bream (Sparus aurata L.) is an important marine fish in Mediterranean aquaculture. Sex determination by age and/or body weight is a critical life‐history trait, the genetic basis for which is largely unknown in this sequential hermaphrodite species. Herein, we performed a partial genome scan to map quantitative trait loci (QTL) affecting body weight and sex using 74 informative microsatellite markers from 10 paternal half‐sib families to construct nine linkage groups (LG). In total, four growth‐related QTL (two chromosome‐wide and two genome‐wide) and six QTL related to sex determination (three pairs in three different LGs) were detected (two chromosome‐wide and one genome‐wide). The proportion of phenotypic variation explained by the body‐weight QTL ranged from 9.3% to 17.2%, showing their potential for use in marker‐assisted selection. The results obtained offer solid ground to investigate the structure and function of the genomic regions involved in the mechanisms of sex reversal.     

猪重要胴体性状的遗传定位

为了寻找影响猪重要胴体性状主基因在染色体的位置,我们以大白猪和梅山猪为父母本建立了F2资源家系。随机选留81头F2代个体,经屠宰获得猪胴体性状数据。结合家系个体的48个微卫星标记基因型,用线性模型最小二乘法对各胴体性状进行数量性状基因座（QTL）的区间定位。定位结果表明位于猪染色体（SSC）4号的瘦肉率和瘦肉量QTL达到基因组极显著水平;SSC1、2和4上眼肌面积QTL达到染色体显著水平;位于SSC1和4上的眼肌高度QTL与眼肌面积QTL在同一染色体区域;而眼肌宽度QTL位于SSC6;位于SSC7同一标记区间的皮重、皮率、骨重和骨率QTL表现出很好的一致性,均达到染色体显著水平。SSC6和7的体长QTL达到染色体显著水平。 Abstract: To detect quantitative trait loci (QTL) for body composition traits in pigs, a resource family with three-generation was developed by using Large White grand sires and Meishan grand dams. A total of 81 F2 progenies were phenotyped for body composition. All animals were genotyped for microsatellite markers. The main results are as follows:, the strongest linkages at genome-wise level of lean meat percentage and total meat content were detected on SSC1 and 4. QTLs for loin eye area were located on SSC1, 2 and 4, QTLs for loin eye height on SSC 1 and 4, and QTLs for loin eye width on SSC 6. The best positions estimated for QTLs of skin percentage and of skin weight were in the same marker interval. Two QTLs significant at genome-wise level or highly significant at chromosome-wide level for carcass length were located on SSC6 and 7.     

Genetic architecture modulates diet-induced hepatic mRNA and miRNA expression profiles in Diversity Outbred mice

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Finally, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility     

Integration of QTL and bioinformatic tools to identify candidate genes for triglycerides in mice

To identify genetic loci influencing lipid levels, we performed quantitative trait loci (QTL) analysis between inbred mouse strains MRL/MpJ and SM/J, measuring triglyceride levels at 8 weeks of age in F2 mice fed a chow diet. We identified one significant QTL on chromosome (Chr) 15 and three suggestive QTL on Chrs 2, 7, and 17. We also carried out microarray analysis on the livers of parental strains of 282 F2 mice and used these data to find cis-regulated expression QTL. We then narrowed the list of candidate genes under significant QTL using a "toolbox" of bioinformatic resources, including haplotype analysis; parental strain comparison for gene expression differences and nonsynonymous coding single nucleotide polymorphisms (SNP); cis-regulated eQTL in livers of F2 mice; correlation between gene expression and phenotype; and conditioning of expression on the phenotype. We suggest Slc25a7 as a candidate gene for the Chr 7 QTL and, based on expression differences, five genes (Polr3 h, Cyp2d22, Cyp2d26, Tspo, and Ttll12) as candidate genes for Chr 15 QTL. This study shows how bioinformatics can be used effectively to reduce candidate gene lists for QTL related to complex traits.     

QTL analysis of body weight and conformation score in commercial broiler chickens using variance component and half-sib analyses

The aim of the study was to investigate quantitative trait loci (QTL) in previously identified regions of chicken chromosomes 1, 4 and 5 relating to 40-day body weights and conformation scores. Half-sib (HS) and variance component analyses were implemented and compared using QTL Express software. Data were from a two-generation design and consisted of 100 dam families nested in 46 sire families with trait values for 2,708 offspring. Chicken chromosome 4 showed nominal significance for QTL affecting body weight and conformation, and linkage was confirmed for both traits on chromosome 5. Results varied according to method of analysis and with common parent in the HS method.     

Quantitative trait locus analysis for obesity reveals multiple networks of interacting loci

Obesity is a highly heritable and genetically complex trait with hundreds of potential loci identified. An intercross of 513 F₂ progeny between the SM/J × NZB/BINJ inbred mouse strains was generated to identify quantitative trait loci (QTL) that are involved in the weight of four fat pads: mesenteric, inguinal, gonadal, and retroperitoneal. Sex and lean body weight were treated as covariates in the analysis of these fat pads. This analysis uncoupled genetic effects related to overall body size from those influencing the adiposity of a mouse. We identified multiple significant QTL. QTL alleles associated with increased lean body weight and individual fat pad weights are contributed by the NZB background. Adiposity loci are distinct from these body size QTLs and high-adiposity alleles are contributed by the SM background. An extended network of epistatic QTL is also observed. A QTL on Chr 19 is the center of a network of eight interacting QTL, Chr 4 is the center of six, and Chr 17 the center of four interacting QTL. We conclude that interacting networks of multiple genes characterize the regulation of fat pad depots and body weight. Haplotype patterns and a literature-driven approach were used to generate hypotheses regarding the identity of the genes and pathways underlying the QTL.