Atila van Nas

Mapping the Genetic Architecture of Gene Expression in Human Liver

Genetic variants that are associated with common human diseases do not lead directly to disease, but instead act on intermediate, molecular phenotypes that in turn induce changes in higher-order disease traits. Therefore, identifying the molecular phenotypes that vary in response to changes in DNA and that also associate with changes in disease traits has the potential to provide the functional information required to not only identify and validate the susceptibility genes that are directly affected by changes in DNA, but also to understand the molecular networks in which such genes operate and how changes in these networks lead to changes in disease traits. Toward that end, we profiled more than 39,000 transcripts and we genotyped 782,476 unique single nucleotide polymorphisms (SNPs) in more than 400 human liver samples to characterize the genetic architecture of gene expression in the human liver, a metabolically active tissue that is important in a number of common human diseases, including obesity, diabetes, and atherosclerosis. This genome-wide association study of gene expression resulted in the detection of more than 6,000 associations between SNP genotypes and liver gene expression traits, where many of the corresponding genes identified have already been implicated in a number of human diseases. The utility of these data for elucidating the causes of common human diseases is demonstrated by integrating them with genotypic and expression data from other human and mouse populations. This provides much-needed functional support for the candidate susceptibility genes being identified at a growing number of genetic loci that have been identified as key drivers of disease from genome-wide association studies of disease. By using an integrative genomics approach, we highlight how the gene RPS26 and not ERBB3 is supported by our data as the most likely susceptibility gene for a novel type 1 diabetes locus recently identified in a large-scale, genome-wide association study. We also identify SORT1 and CELSR2 as candidate susceptibility genes for a locus recently associated with coronary artery disease and plasma low-density lipoprotein cholesterol levels in the process.

Dosage compensation is less effective in birds than in mammals

Background In animals with heteromorphic sex chromosomes, dosage compensation of sex-chromosome genes is thought to be critical for species survival. Diverse molecular mechanisms have evolved to effectively balance the expressed dose of X-linked genes between XX and XY animals, and to balance expression of X and autosomal genes. Dosage compensation is not understood in birds, in which females (ZW) and males (ZZ) differ in the number of Z chromosomes. Results Using microarray analysis, we compared the male:female ratio of expression of sets of Z-linked and autosomal genes in two bird species, zebra finch and chicken, and in two mammalian species, mouse and human. Male:female ratios of expression were significantly higher for Z genes than for autosomal genes in several finch and chicken tissues. In contrast, in mouse and human the male:female ratio of expression of X-linked genes is quite similar to that of autosomal genes, indicating effective dosage compensation even in humans, in which a significant percentage of genes escape X-inactivation. Conclusion Birds represent an unprecedented case in which genes on one sex chromosome are expressed on average at constitutively higher levels in one sex compared with the other. Sex-chromosome dosage compensation is surprisingly ineffective in birds, suggesting that some genomes can do without effective sex-specific sex-chromosome dosage compensation mechanisms.

Validation of candidate causal genes for obesity that affect shared metabolic pathways and networks

A principal task in dissecting the genetics of complex traits is to identify causal genes for disease phenotypes. We previously developed a method to infer causal relationships among genes through the integration of DNA variation, gene transcription and phenotypic information. Here we have validated our method through the characterization of transgenic and knockout mouse models of genes predicted to be causal for abdominal obesity. Perturbation of eight out of the nine genes, with Gas7, Me1 and Gpx3 being newly confirmed, resulted in significant changes in obesity-related traits. Liver expression signatures revealed alterations in common metabolic pathways and networks contributing to abdominal obesity and overlapped with a macrophage-enriched metabolic network module that is highly associated with metabolic traits in mice and humans. Integration of gene expression in the design and analysis of traditional F2 intercross studies allows high-confidence prediction of causal genes and identification of pathways and networks involved.

An Integration of Genome-Wide Association Study and Gene Expression Profiling to Prioritize the Discovery of Novel Susceptibility Loci for Osteoporosis-Related Traits

Osteoporosis is a complex disorder and commonly leads to fractures in elderly persons. Genome-wide association studies (GWAS) have become an unbiased approach to identify variations in the genome that potentially affect health. However, the genetic variants identified so far only explain a small proportion of the heritability for complex traits. Due to the modest genetic effect size and inadequate power, true association signals may not be revealed based on a stringent genome-wide significance threshold. Here, we take advantage of SNP and transcript arrays and integrate GWAS and expression signature profiling relevant to the skeletal system in cellular and animal models to prioritize the discovery of novel candidate genes for osteoporosis-related traits, including bone mineral density (BMD) at the lumbar spine (LS) and femoral neck (FN), as well as geometric indices of the hip (femoral neck-shaft angle, NSA; femoral neck length, NL; and narrow-neck width, NW). A two-stage meta-analysis of GWAS from 7,633 Caucasian women and 3,657 men, revealed three novel loci associated with osteoporosis-related traits, including chromosome 1p13.2 (RAP1A, p = 3.6×10−8), 2q11.2 (TBC1D8), and 18q11.2 (OSBPL1A), and confirmed a previously reported region near TNFRSF11B/OPG gene. We also prioritized 16 suggestive genome-wide significant candidate genes based on their potential involvement in skeletal metabolism. Among them, 3 candidate genes were associated with BMD in women. Notably, 2 out of these 3 genes (GPR177, p = 2.6×10−13; SOX6, p = 6.4×10−10) associated with BMD in women have been successfully replicated in a large-scale meta-analysis of BMD, but none of the non-prioritized candidates (associated with BMD) did. Our results support the concept of our prioritization strategy. In the absence of direct biological support for identified genes, we highlighted the efficiency of subsequent functional characterization using publicly available expression profiling relevant to the skeletal system in cellular or whole animal models to prioritize candidate genes for further functional validation.

Elucidating the role of gonadal hormones in sexually dimorphic gene coexpression networks.

We previously used high-density expression arrays to interrogate a genetic cross between strains C3H/HeJ and C57BL/6J and observed thousands of differences in gene expression between sexes. We now report analyses of the molecular basis of these sex differences and of the effects of sex on gene expression networks. We analyzed liver gene expression of hormone-treated gonadectomized mice as well as XX male and XY female mice. Differences in gene expression resulted in large part from acute effects of gonadal hormones acting in adulthood, and the effects of sex chromosomes, apart from hormones, were modest. We also determined whether there are sex differences in the organization of gene expression networks in adipose, liver, skeletal muscle, and brain tissue. Although coexpression networks of highly correlated genes were largely conserved between sexes, some exhibited striking sex dependence. We observed strong body fat and lipid correlations with sex-specific modules in adipose and liver as well as a sexually dimorphic network enriched for genes affected by gonadal hormones. Finally, our analyses identified chromosomal loci regulating sexually dimorphic networks. This study indicates that gonadal hormones play a strong role in sex differences in gene expression. In addition, it results in the identification of sex-specific gene coexpression networks related to genetic and metabolic traits.

Expression Quantitative Trait Loci: Replication, Tissue- and Sex-Specificity in Mice

By treating the transcript abundance as a quantitative trait, gene expression can be mapped to local or distant genomic regions relative to the gene encoding the transcript. Local expression quantitative trait loci (eQTL) generally act in cis (that is, control the expression of only the contiguous structural gene), whereas distal eQTL act in trans. Distal eQTL are more difficult to identify with certainty due to the fact that significant thresholds are very high since all regions of the genome must be tested, and confounding factors such as batch effects can produce false positives. Here, we compare findings from two large genetic crosses between mouse strains C3H/HeJ and C57BL/6J to evaluate the reliability of distal eQTL detection, including “hotspots” influencing the expression of multiple genes in trans. We found that >63% of local eQTL and >18% of distal eQTL were replicable at a threshold of LOD > 4.3 between crosses and 76% of local and >24% of distal eQTL at a threshold of LOD > 6. Additionally, at LOD > 4.3 four tissues studied (adipose, brain, liver, and muscle) exhibited >50% preservation of local eQTL and >17% preservation of distal eQTL. We observed replicated distal eQTL hotspots between the crosses on chromosomes 9 and 17. Finally, >69% of local eQTL and >10% of distal eQTL were preserved in most tissues between sexes. We conclude that most local eQTL are highly replicable between mouse crosses, tissues, and sex as compared to distal eQTL, which exhibited modest replicability.

Copy number variation influences gene expression and metabolic traits in mice

Copy number variants (CNVs) are genomic segments which are duplicated or deleted among different individuals. CNVs have been implicated in both Mendelian and complex traits, including immune and behavioral disorders, but the study of the mechanisms by which CNVs influence gene expression and clinical phenotypes in humans is complicated by the limited access to tissues and by population heterogeneity. We now report studies of the effect of 19 CNVs on gene expression and metabolic traits in a mouse intercross between strains C57BL/6J and C3H/HeJ. We found that 83% of genes predicted to occur within CNVs were differentially expressed. The expression of most CNV genes was correlated with copy number, but we also observed evidence that gene expression was altered in genes flanking CNVs, suggesting that CNVs may contain regulatory elements for these genes. Several CNVs mapped to hotspots, genomic regions influencing expression of tens or hundreds of genes. Several metabolic traits including cholesterol, triglycerides, glucose and body weight mapped to three CNVs in the genome, in mouse chromosomes 1, 4 and 17. Predicted CNV genes, such as Itlna, Defcr-1, Trim12 and Trim34 were highly correlated with these traits. Our results suggest that CNVs have a significant impact on gene expression and that CNVs may be playing a role in the mechanisms underlying metabolic traits in mice.

Maternal Low-Protein Diet or Hypercholesterolemia Reduces Circulating Essential Amino Acids and Leads to Intrauterine Growth Restriction

OBJECTIVE—We have examined maternal mechanisms for adult-onset glucose intolerance, increased adiposity, and atherosclerosis using two mouse models for intrauterine growth restriction (IUGR): maternal protein restriction and hypercholesterolemia. RESEARCH DESIGN AND METHODS—For these studies, we measured the amino acid levels in dams from two mouse models for IUGR: 1) feeding C57BL/6J dams a protein-restricted diet and 2) feeding C57BL/6J LDL receptor–null (LDLR−/−) dams a high-fat (Western) diet. RESULTS—Both protein-restricted and hypercholesterolemic dams exhibited significantly decreased concentrations of the essential amino acid phenylalanine and the essential branched chain amino acids leucine, isoleucine, and valine. The protein-restricted diet for pregnant dams resulted in litters with significant IUGR. Protein-restricted male offspring exhibited catch-up growth by 8 weeks of age and developed increased adiposity and glucose intolerance by 32 weeks of age. LDLR−/− pregnant dams on a Western diet also had litters with significant IUGR. Male and female LDLR−/− Western-diet offspring developed significantly larger atherosclerotic lesions by 90 days compared with chow-diet offspring. CONCLUSIONS—In two mouse models of IUGR, we found reduced concentrations of essential amino acids in the experimental dams. This indicated that shared mechanisms may underlie the phenotypic effects of maternal hypercholesterolemia and maternal protein restriction on the offspring.

An integrative genetics approach to identify candidate genes regulating BMD: combining linkage, gene expression, and association.

Numerous quantitative trait loci (QTLs) affecting bone traits have been identified in the mouse; however, few of the underlying genes have been discovered. To improve the process of transitioning from QTL to gene, we describe an integrative genetics approach, which combines linkage analysis, expression QTL (eQTL) mapping, causality modeling, and genetic association in outbred mice. In C57BL/6J × C3H/HeJ (BXH) F2 mice, nine QTLs regulating femoral BMD were identified. To select candidate genes from within each QTL region, microarray gene expression profiles from individual F2 mice were used to identify 148 genes whose expression was correlated with BMD and regulated by local eQTLs. Many of the genes that were the most highly correlated with BMD have been previously shown to modulate bone mass or skeletal development. Candidates were further prioritized by determining whether their expression was predicted to underlie variation in BMD. Using network edge orienting (NEO), a causality modeling algorithm, 18 of the 148 candidates were predicted to be causally related to differences in BMD. To fine‐map QTLs, markers in outbred MF1 mice were tested for association with BMD. Three chromosome 11 SNPs were identified that were associated with BMD within the Bmd11 QTL. Finally, our approach provides strong support for Wnt9a, Rasd1, or both underlying Bmd11. Integration of multiple genetic and genomic data sets can substantially improve the efficiency of QTL fine‐mapping and candidate gene identification.

Systems genetics of susceptibility to obesity-induced diabetes in mice

Inbred strains of mice are strikingly different in susceptibility to obesity-driven diabetes. For instance, deficiency in leptin receptor (db/db) leads to hyperphagia and obesity in both C57BL/6 and DBA/2 mice, but only on the DBA/2 background do the mice develop beta-cell loss leading to severe diabetes, while C57BL/6 mice are relatively resistant. To further investigate the genetic factors predisposing to diabetes, we have studied leptin receptor-deficient offspring of an F2 cross between C57BL/6J (db/+) males and DBA/2J females. The results show that the genetics of diabetes susceptibility are enormously complex and a number of quantitative trait loci (QTL) contributing to diabetes-related traits were identified, notably on chromosomes 4, 6, 7, 9, 10, 11, 12, and 19. The Chr. 4 locus is likely due to a disruption of the Zfp69 gene in C57BL/6J mice. To identify candidate genes and to model coexpression networks, we performed global expression array analysis in livers of the F2 mice. Expression QTL (eQTL) were identified and used to prioritize candidate genes at clinical trait QTL. In several cases, clusters of eQTLs colocalized with clinical trait QTLs, suggesting a common genetic basis. We constructed coexpression networks for both 5 and 12 wk old mice and identified several modules significantly associated with clinical traits. One module in 12 wk old mice was associated with several measures of hepatic fat content as well as with other lipid- and diabetes-related traits. These results add to the understanding of the complex genetic interactions contributing to obesity-induced diabetes.