Replication and narrowing of gene expression quantitative trait loci using inbred mice

Population-based dose-response analysis of liver transcriptional response to trichloroethylene in mouse

Studies of gene expression are common in toxicology and provide important clues to mechanistic understanding of adverse effects of chemicals. Most prior studies have been performed in a single strain or cell line; however, gene expression is heavily influenced by the genetic background, and these genotype-expression differences may be key drivers of inter-individual variation in response to chemical toxicity. In this study, we hypothesized that the genetically diverse Collaborative Cross mouse population can be used to gain insight and suggest mechanistic hypotheses for the dose- and genetic background-dependent effects of chemical exposure. This hypothesis was tested using a model liver toxicant trichloroethylene (TCE). Liver transcriptional responses to TCE exposure were evaluated 24 h after dosing. Transcriptomic dose-responses were examined for both TCE and its major oxidative metabolite trichloroacetic acid (TCA). As expected, peroxisome- and fatty acid metabolism-related pathways were among the most dose-responsive enriched pathways in all strains. However, nearly half of the TCE-induced liver transcriptional perturbation was strain-dependent, with abundant evidence of strain/dose interaction, including in the peroxisomal signaling-associated pathways. These effects were highly concordant between the administered TCE dose and liver levels of TCA. Dose-response analysis of gene expression at the pathway level yielded points of departure similar to those derived from the traditional toxicology studies for both non-cancer and cancer effects. Mapping of expression-genotype-dose relationships revealed some significant associations; however, the effects of TCE on gene expression in liver appear to be highly polygenic traits that are challenging to positionally map. This study highlights the usefulness of mouse population-based studies in assessing inter-individual variation in toxicological responses, but cautions that genetic mapping may be challenging because of the complexity in gene exposure-dose relationships.

Genetic and epigenetic determinants of inter-individual variability in responses to toxicants

It is well established that genetic variability has a major impact on susceptibility to common diseases, responses to drugs and toxicants, and influences disease-related outcomes. The appreciation that epigenetic marks also vary across the population is growing with more data becoming available from studies in humans and model organisms. In addition, the links between genetic variability, toxicity outcomes and epigenetics are being actively explored. Recent studies demonstrate that gene-by-environment interactions involve both chromatin states and transcriptional regulation, and that epigenetics provides important mechanistic clues to connect expression-related quantitative trait loci (QTL) and disease outcomes. However, studies of Gene×Environment×Epigenetics further extend the complexity of the experimental designs and create a challenge for selecting the most informative epigenetic readouts that can be feasibly performed to interrogate multiple individuals, exposures, tissue types and toxicity phenotypes. We propose that among the many possible epigenetic experimental methodologies, assessment of chromatin accessibility coupled with total RNA levels provides a cost-effective and comprehensive option to sufficiently characterize the complexity of epigenetic and regulatory activity in the context of understanding the inter-individual variability in responses to toxicants.

Variable cardiac α-actin (Actc1) expression in early adult skeletal muscle correlates with promoter methylation

Different genes encode the α-actin isoforms that are predominantly expressed in heart and skeletal muscle. Mutations in the skeletal muscle α-actin gene (ACTA1) cause muscle diseases that are mostly lethal in the early postnatal period. We previously demonstrated that the disease phenotype of ACTA1 mouse models could be rescued by transgenic over-expression of cardiac α-actin (ACTC1). ACTC1 is the predominant striated α-actin isoform in the heart but is also expressed in developing skeletal muscle. To develop a translatable therapy, we investigated the genetic regulation of Actc1 expression. Using strains from The Collaborative Cross (CC) genetic resource, we found that Actc1 varies in expression by up to 24-fold in skeletal muscle. We defined significant expression quantitative trait loci (eQTL) associated with early adult Actc1 expression in soleus and heart. eQTL in both heart and soleus mapped to the Actc1 locus and replicate an eQTL mapped for Actc1 in BXD heart and quadriceps. We built on this previous work by analysing genes within the eQTL peak regions to prioritise likely candidates for modifying Actc1 expression. Additionally we interrogated the CC founder haplotype contributions to enable prioritisation of genetic variants for functional analyses. Methylation around the Actc1 transcriptional start site in early adult skeletal muscle negatively correlated with Actc1 expression in a strain-dependent manner, while other marks of regulatory potential (histone modification and chromatin accessibility) were unaltered. This study provides novel insights into the complex genetic regulation of Actc1 expression in early adult skeletal muscles.

Expression quantitative trait loci infer the regulation of isoflavone accumulation in soybean (Glycine max L. Merr.) seed

BACKGROUND

Mapping expression quantitative trait loci (eQTL) of targeted genes represents a powerful and widely adopted approach to identify putative regulatory variants. Linking regulation differences to specific genes might assist in the identification of networks and interactions. The objective of this study is to identify eQTL underlying expression of four gene families encoding isoflavone synthetic enzymes involved in the phenylpropanoid pathway, which are phenylalanine ammonia-lyase (PAL; EC 4.3.1.5), chalcone synthase (CHS; EC 2.3.1.74), 2-hydroxyisoflavanone synthase (IFS; EC1.14.13.136) and flavanone 3-hydroxylase (F3H; EC 1.14.11.9). A population of 130 recombinant inbred lines (F5:11), derived from a cross between soybean cultivar 'Zhongdou 27' (high isoflavone) and 'Jiunong 20' (low isoflavone), and a total of 194 simple sequence repeat (SSR) markers were used in this study. Overlapped loci of eQTLs and phenotypic QTLs (pQTLs) were analyzed to identify the potential candidate genes underlying the accumulation of isoflavone in soybean seed.

RESULTS

Thirty three eQTLs (thirteen cis-eQTLs and twenty trans-eQTLs) underlying the transcript abundance of the four gene families were identified on fifteen chromosomes. The eQTLs between Satt278-Sat_134, Sat_134-Sct_010 and Satt149-Sat_234 underlie the expression of both IFS and CHS genes. Five eQTL intervals were overlapped with pQTLs. A total of eleven candidate genes within the overlapped eQTL and pQTL were identified.

CONCLUSIONS

These results will be useful for the development of marker-assisted selection to breed soybean cultivars with high or low isoflavone contents and for map-based cloning of new isoflavone related genes.

Variation in Rubisco activase (RCAβ) gene promoters and expression in soybean [Glycine max (L.) Merr]

Understanding the genetic basis of Rubisco activase (RCA) gene regulation and altering its expression levels to optimize Rubisco activation may provide an approach to enhance plant productivity. However, the genetic mechanisms and the effect of RCA expression on phenotype are still unknown in soybean. This work analysed the expression of RCA genes and demonstrated that two RCA isoforms presented different expression patterns. Compared with GmRCAα, GmRCAβ was expressed at higher mRNA and protein levels. In addition, GmRCAα and GmRCAβ were positively correlated with chlorophyll fluorescence parameters and seed yield, suggesting that changes in expression of RCA has a potential applicability in breeding for enhanced soybean productivity. To identify the genetic factors that cause expression level variation of GmRCAβ, expression quantitative trait loci (eQTL) mapping was combined with allele mining in a natural population including 219 landraces. The eQTL mapping showed that a combination of both cis- and trans-acting eQTLs might control GmRCAβ expression. As promoters can affect both cis- and trans-acting eQTLs by altering cis-acting regulatory elements or transcription factor binding sites, this work subsequently focused on the promoter region of GmRCAβ. Single-nucleotide polymorphisms in the GmRCAβ promoter were identified and shown to correlate with expression level diversity. These SNPs were classified into two groups, A and B. Further transient expression showed that GUS expression driven by the group A promoter was stronger than that by the group B promoter, suggesting that promoter sequence types could influence gene expression levels. These results would improve understanding how variation within promoters affects gene expression and, ultimately, phenotypic diversity in natural populations.

Genetic-genomic replication to identify candidate mouse atherosclerosis modifier genes

Objective

Genetics plays a large role in atherosclerosis susceptibility in humans and mice. We attempted to confirm previously determined mouse atherosclerosis‐associated loci and use bioinformatics and transcriptomics to create a catalog of candidate atherosclerosis modifier genes at these loci.

Methods and Results

A strain intercross was performed between AKR and DBA/2 mice on the apoE^−/− background generating 166 F₂ progeny. Using the phenotype log10 of the aortic root lesion area, we identified 3 suggestive atherosclerosis quantitative trait loci (Ath QTLs). When combined with our prior strain intercross, we confirmed 3 significant Ath QTLs on chromosomes 2, 15, and 17, with combined logarithm of odds scores of 5.9, 5.3, and 5.6, respectively, which each met the genome‐wide 5% false discovery rate threshold. We identified all of the protein coding differences between these 2 mouse strains within the Ath QTL intervals. Microarray gene expression profiling was performed on macrophages and endothelial cells from this intercross to identify expression QTLs (eQTLs), the loci that are associated with variation in the expression levels of specific transcripts. Cross tissue eQTLs and macrophage eQTLs that replicated from a prior strain intercross were identified. These bioinformatic and eQTL analyses produced a comprehensive list of candidate genes that may be responsible for the Ath QTLs.

Conclusions

Replication studies for clinical traits as well as gene expression traits are worthwhile in identifying true versus false genetic associations. We have replicated 3 loci on mouse chromosomes 2, 15, and 17 that are associated with atherosclerosis. We have also identified protein coding differences and multiple replicated eQTLs, which may be useful in the identification of atherosclerosis modifier genes.

An interval of the obesity QTL Nob3.38 within a QTL hotspot on chromosome 1 modulates behavioral phenotypes

A region on mouse distal chromosome 1 (Chr. 1) that is highly enriched in quantitative trait loci (QTLs) controlling neural and behavioral phenotypes overlaps with the peak region of a major obesity QTL (Nob3.38), which we identified in an intercross of New Zealand Obese (NZO) mice with C57BL/6J (B6). By positional cloning we recently identified a microdeletion within this locus causing the disruption of Ifi202b that protects from adiposity by suppressing expression of 11β-Hsd1. Here we show that the Nob3.38 segment also corresponds with the QTL rich region (Qrr1) on Chr. 1 and associates with increased voluntary running wheel activity, Rota-rod performance, decreased grip strength, and anxiety-related traits. The characterization of a subcongenic line carrying 14.2 Mbp of Nob3.38 with a polymorphic region of 4.4 Mbp indicates that the microdeletion and/or other polymorphisms in its proximity alter body weight, voluntary activity, and exploration. Since 27 out of 32 QTL were identified in crosses with B6, we hypothesized that the microdeletion and or adjacent SNPs are unique for B6 mice and responsible for some of the complex Qrr1-mediated effects. Indeed, a phylogenic study of 28 mouse strains revealed a NZO-like genotype for 22 and a B6-like genotype for NZW/LacJ and 4 other C57BL strains. Thus, we suggest that a Nob3.38 interval (173.0-177.4 Mbp) does not only modify adiposity but also neurobehavioral traits by a haplotype segregating with C57BL strains.

Computational tools for discovery and interpretation of expression quantitative trait loci

Expression quantitative trait locus (eQTL) analysis is rapidly moving from a cutting-edge concept in genomics to a mature area of investigation, with important connections to genome-wide association studies for human disease, pharmacogenomics and toxicogenomics. Despite the importance of the topic, many investigators must develop their own code or use tools not specifically suited for eQTL analysis. Convenient computational tools are becoming available, but they are not widely publicized, and investigators who are interested in discovery or eQTL, or in using them to interpret genome-wide association study results may have difficulty navigating the available resources. The purpose of this review is to help investigators find appropriate programs for eQTL analysis and interpretation.

A comparative assessment of mandible shape in a consomic strain panel of the house mouse (Mus musculus)--implications for epistasis and evolvability of quantitative traits

Background

Expectations of repeatedly finding associations between given genes and phenotypes have been borne out by studies of parallel evolution, especially for traits involving absence or presence of characters. However, it has rarely been asked whether the genetic basis of quantitative trait variation is conserved at the intra- or even at the interspecific level. This question is especially relevant for shape, where the high dimensionality of variation seems to require a highly complex genetic architecture involving many genes.

Results

We analyse here the genetic effects of chromosome substitution strains carrying M. m. musculus chromosomes in a largely M. m. domesticus background on mandible shape and compare them to the results of previously published QTL mapping data between M. m. domesticus strains. We find that the distribution of genetic effects and effect sizes across the genome is consistent between the studies, while the specific shape changes associated with the chromosomes are different. We find also that the sum of the effects from the different M. m. musculus chromosomes is very different from the shape of the strain from which they were derived, as well as all known wild type shapes.

Conclusions

Our results suggest that the relative chromosome-wide effect sizes are comparable between the long separated subspecies M. m. domesticus and M. m. musculus, hinting at a relative stability of genes involved in this complex trait. However, the absolute effect sizes and the effect directions may be allele-dependent, or are context dependent, i.e. epistatic interactions appear to play an important role in controlling shape.

MicroRNA expression in the livers of inbred mice

MicroRNAs are short, non-coding RNA sequences that regulate genes at the post-transcriptional level and have been shown to be important in development, tissue differentiation, and disease. Limited attention has been given to the natural variation in miRNA expression across genetically diverse populations even though it is well established that genetic polymorphisms can have a profound effect on mRNA levels. Expression level of 577 miRNAs in the livers of 70 strains of inbred mice was assessed, and we found that miRNA expression is highly stable across different strains. Globally, the expression of miRNA target transcripts does not correlate with miRNA expression, primarily due to the low variance of miRNA but high variance of mRNA expression across strains. Our results show that there is little genetic effect on the baseline miRNA levels in murine liver. The stability of mouse liver miRNA expression in a genetically diverse population suggests that treatment-induced disruptions in liver miRNA expression, a phenomenon established for a large number of toxicants, may indicate an important mechanism for the disturbance of normal liver function, and may prove to be a useful genetic background-independent biomarker of toxicant effect.

Hippocampal gene expression analysis highlights Ly6a/Sca-1 as candidate gene for previously mapped novelty induced behaviors in mice

In this study, we show that the covariance between behavior and gene expression in the brain can help further unravel the determinants of neurobehavioral traits. Previously, a QTL for novelty induced motor activity levels was identified on murine chromosome 15 using consomic strains. With the goal of narrowing down the linked region and possibly identifying the gene underlying the quantitative trait, gene expression data from this F₂-population was collected and used for expression QTL analysis. While genetic variation in these mice was limited to chromosome 15, eQTL analysis of gene expression showed strong cis-effects as well as trans-effects elsewhere in the genome. Using weighted gene co-expression network analysis, we were able to identify modules of co-expressed genes related to novelty induced motor activity levels. In eQTL analyses, the expression of Ly6a (a.k.a. Sca-1) was found to be cis-regulated by chromosome 15. Ly6a also surfaced in a group of genes resulting from the network analysis that was correlated with behavior. Behavioral analysis of Ly6a knock-out mice revealed reduced novelty induced motor activity levels when compared to wild type controls, confirming functional importance of Ly6a in this behavior, possibly through regulating other genes in a pathway. This study shows that gene expression profiling can be used to narrow down a previously identified behavioral QTL in mice, providing support for Ly6a as a candidate gene for functional involvement in novelty responsiveness.

Interstrain differences in the liver effects of trichloroethylene in a multistrain panel of inbred mice

Trichloroethylene (TCE) is a widely used industrial chemical and a common environmental contaminant. It is a well-known carcinogen in rodents and a probable carcinogen in humans. Studies utilizing panels of mouse inbred strains afford a unique opportunity to understand both metabolic and genetic basis for differences in responses to TCE. We tested the hypothesis that strain- and liver-specific toxic effects of TCE are genetically controlled and that the mechanisms of toxicity and susceptibility can be uncovered by exploring responses to TCE using a diverse panel of inbred mouse strains. TCE (2100 mg/kg) or corn oil vehicle was administered by gavage to 6- to 8-week-old male mice of 15 mouse strains. Serum and liver were collected at 2, 8, and 24 h postdosing and were analyzed for TCE metabolites, hepatocellular injury, and gene expression of liver. TCE metabolism, as evident from the levels of individual oxidative and conjugative metabolites, varied considerably between strains. TCE treatment-specific effect on the liver transcriptome was strongly dependent on genetic background. Peroxisome proliferator-activated receptor-mediated molecular networks, consisting of the metabolism genes known to be induced by TCE, represent some of the most pronounced molecular effects of TCE treatment in mouse liver that are dependent on genetic background. Conversely, cell death, liver necrosis, and immune-mediated response pathways, which are altered by TCE treatment in liver, are largely genetic background independent. These studies provide better understanding of the mechanisms of TCE-induced toxicity anchored on metabolism and genotype-phenotype correlations that may define susceptibility or resistance.

Integrative analysis of low- and high-resolution eQTL

The study of expression quantitative trait loci (eQTL) is a powerful way of detecting transcriptional regulators at a genomic scale and for elucidating how natural genetic variation impacts gene expression. Power and genetic resolution are heavily affected by the study population: whereas recombinant inbred (RI) strains yield greater statistical power with low genetic resolution, using diverse inbred or outbred strains improves genetic resolution at the cost of lower power. In order to overcome the limitations of both individual approaches, we combine data from RI strains with genetically more diverse strains and analyze hippocampus eQTL data obtained from mouse RI strains (BXD) and from a panel of diverse inbred strains (Mouse Diversity Panel, MDP). We perform a systematic analysis of the consistency of eQTL independently obtained from these two populations and demonstrate that a significant fraction of eQTL can be replicated. Based on existing knowledge from pathway databases we assess different approaches for using the high-resolution MDP data for fine mapping BXD eQTL. Finally, we apply this framework to an eQTL hotspot on chromosome 1 (Qrr1), which has been implicated in a range of neurological traits. Here we present the first systematic examination of the consistency between eQTL obtained independently from the BXD and MDP populations. Our analysis of fine-mapping approaches is based on 'real life' data as opposed to simulated data and it allows us to propose a strategy for using MDP data to fine map BXD eQTL. Application of this framework to Qrr1 reveals that this eQTL hotspot is not caused by just one (or few) 'master regulators', but actually by a set of polymorphic genes specific to the central nervous system.

Sex-specific gene expression in the BXD mouse liver

Differences in clinical phenotypes between the sexes are well documented and have their roots in differential gene expression. While sex has a major effect on gene expression, transcription is also influenced by complex interactions between individual genetic variation and environmental stimuli. In this study, we sought to understand how genetic variation affects sex-related differences in liver gene expression by performing genetic mapping of genomewide liver mRNA expression data in a genetically defined population of naive male and female mice from C57BL/6J, DBA/2J, B6D2F1, and 37 C57BL/6J x DBA/2J (BXD) recombinant inbred strains. As expected, we found that many genes important to xenobiotic metabolism and other important pathways exhibit sexually dimorphic expression. We also performed gene expression quantitative trait locus mapping in this panel and report that the most significant loci that appear to regulate a larger number of genes than expected by chance are largely sex independent. Importantly, we found that the degree of correlation within gene expression networks differs substantially between the sexes. Finally, we compare our results to a recently released human liver gene expression data set and report on important similarities in sexually dimorphic liver gene expression between mouse and human. This study enhances our understanding of sex differences at the genome level and between species, as well as increasing our knowledge of the molecular underpinnings of sex differences in responses to xenobiotics.