eQTL Viewer: visualizing how sequence variation affects genome-wide transcription

BACKGROUND: Expression Quantitative Trait Locus (eQTL) mapping methods have been used to identify the genetic basis of gene expression variations. To map eQTL, thousands of expression profiles are related with sequence polymorphisms across the genome through their correlated variations. These eQTL distribute in many chromosomal regions, each of which can include many genes. The large number of mapping results produced makes it difficult to consider simultaneously the relationships between multiple genomic regions and multiple expressional profiles. There is a need for informative bioinformatics tools to assist the visualization and interpretation of these mapping results. RESULTS: We have developed a web-based tool, called eQTL Viewer, to visualize the relationships between the expression trait genes and the candidate genes in the eQTL regions using Scalable Vector Graphics. The plot generated by eQTL Viewer has the capacity to display mapping results with high resolutions at a variety of scales, and superimpose biological annotations onto the mapping results dynamically. CONCLUSION: Our tool provides an efficient and intuitive way for biologists to explore transcriptional regulation patterns, and to generate hypotheses on the genetic basis of transcriptional regulations

Age and Genetic Background Modify Hybrid Male Sterility in House Mice.

Hybrid male sterility (HMS) contributes to reproductive isolation commonly observed among house mouse (Mus musculus) subspecies, both in the wild and in laboratory crosses. Incompatibilities involving specific Prdm9 alleles and certain Chromosome (Chr) X genotypes are known determinants of fertility and HMS, and previous work in the field has demonstrated that genetic background modifies these two major loci. We constructed hybrids that have identical genotypes at Prdm9 and identical X chromosomes, but differ widely across the rest of the genome. In each case, we crossed female PWK/PhJ mice representative of the M. m. musculus subspecies to males from a classical inbred strain representative of M. m. domesticus: 129S1/SvImJ, A/J, C57BL/6J, or DBA/2J. We detected three distinct trajectories of fertility among the hybrids using breeding experiments. The PWK129S1 males were always infertile. PWKDBA2 males were fertile, despite their genotypes at the major HMS loci. We also observed age-dependent changes in fertility parameters across multiple genetic backgrounds. The PWKB6 and PWKAJ males were always infertile before 12 weeks and after 35 weeks. However, some PWKB6 and PWKAJ males were transiently fertile between 12 and 35 weeks. This observation could resolve previous contradictory reports about the fertility of PWKB6. Taken together, these results point to multiple segregating HMS modifier alleles, some of which have age-related modes of action. The ultimate identification of these alleles and their age-related mechanisms will advance understanding both of the genetic architecture of HMS and of how reproductive barriers are maintained between house mouse subspecies

Genetic Analysis of Hematological Parameters in Incipient Lines of the Collaborative Cross

Hematological parameters, including red and white blood cell counts and hemoglobin concentration, are widely used clinical indicators of health and disease. These traits are tightly regulated in healthy individuals and are under genetic control. Mutations in key genes that affect hematological parameters have important phenotypic consequences, including multiple variants that affect susceptibility to malarial disease. However, most variation in hematological traits is continuous and is presumably influenced by multiple loci and variants with small phenotypic effects. We used a newly developed mouse resource population, the Collaborative Cross (CC), to identify genetic determinants of hematological parameters. We surveyed the eight founder strains of the CC and performed a mapping study using 131 incipient lines of the CC. Genome scans identified quantitative trait loci for several hematological parameters, including mean red cell volume (Chr 7 and Chr 14), white blood cell count (Chr 18), percent neutrophils/lymphocytes (Chr 11), and monocyte number (Chr 1). We used evolutionary principles and unique bioinformatics resources to reduce the size of candidate intervals and to view functional variation in the context of phylogeny. Many quantitative trait loci regions could be narrowed sufficiently to identify a small number of promising candidate genes. This approach not only expands our knowledge about hematological traits but also demonstrates the unique ability of the CC to elucidate the genetic architecture of complex traits

Genetically Determined Severity of Anti-Myeloperoxidase Glomerulonephritis

Myeloperoxidase (MPO) is a target antigen for antineutrophil cytoplasmic autoantibodies (ANCA). There is evidence that MPO-ANCA cause necrotizing and crescentic glomerulonephritis (NCGN) and vasculitis. NCGN severity varies among patients with ANCA disease, and genetic factors influence disease severity. The role of genetics in MPO-ANCA NCGN severity was investigated using 13 inbred mouse strains, F1 and F2 hybrids, bone marrow chimeras, and neutrophil function assays. Mouse strains include founders of the Collaborative Cross. Intravenous injection of anti-MPO IgG induced glomerular crescents in >60% of glomeruli in 129S6/SvEv and CAST/EiJ mice, but <1% in A/J, DBA/1J, DBA/2J, NOD/LtJ, and PWK/PhJ mice. C57BL6J, 129S1/SvImJ, LP/J, WSB/EiJ, NZO/HILtJ, and C3H mice had intermediate severity. High-density genotypes at 542,190 single nucleotide polymorphisms were used to identify candidate loci for disease severity by identifying genomic regions that are different between 129S6/SvEv and 129S1/SvImJ mice, which are genetically similar but phenotypically distinct. C57BL/6 × 129S6 F2 mice were genotyped at 76 SNPs to capture quantitative trait loci for disease severity. The absence of a dominant quantitative trait locus suggests that differences in severity are the result of multiple gene interactions. In vivo studies using bone marrow chimeric mice and in vitro studies of neutrophil activation by anti-MPO IgG indicated that severity of NCGN is mediated by genetically determined differences in the function of neutrophils

From Classical Genetics to Quantitative Genetics to Systems Biology: Modeling Epistasis

Gene expression data has been used in lieu of phenotype in both classical and quantitative genetic settings. These two disciplines have separate approaches to measuring and interpreting epistasis, which is the interaction between alleles at different loci. We propose a framework for estimating and interpreting epistasis from a classical experiment that combines the strengths of each approach. A regression analysis step accommodates the quantitative nature of expression measurements by estimating the effect of gene deletions plus any interaction. Effects are selected by significance such that a reduced model describes each expression trait. We show how the resulting models correspond to specific hierarchical relationships between two regulator genes and a target gene. These relationships are the basic units of genetic pathways and genomic system diagrams. Our approach can be extended to analyze data from a variety of experiments, multiple loci, and multiple environments

Integrative Genetic Analysis of Allergic Inflammation in the Murine Lung

Airway allergen exposure induces inflammation among individuals with atopy that is characterized by altered airway gene expression, elevated levels of T helper type 2 cytokines, mucus hypersecretion, and airflow obstruction. To identify the genetic determinants of the airway allergen response, we employed a systems genetics approach. We applied a house dust mite mouse model of allergic airway disease to 151 incipient lines of the Collaborative Cross, a new mouse genetic reference population, and measured serum IgE, airway eosinophilia, and gene expression in the lung. Allergen-induced serum IgE and airway eosinophilia were not correlated. We detected quantitative trait loci (QTL) for airway eosinophilia on chromosome (Chr) 11 (71.802–87.098 megabases [Mb]) and allergen-induced IgE on Chr 4 (13.950–31.660 Mb). More than 4,500 genes expressed in the lung had gene expression QTL (eQTL), the majority of which were located near the gene itself. However, we also detected approximately 1,700 trans-eQTL, and many of these trans-eQTL clustered into two regions on Chr 2. We show that one of these loci (at 147.6 Mb) is associated with the expression of more than 100 genes, and, using bioinformatics resources, fine-map this locus to a 53 kb-long interval. We also use the gene expression and eQTL data to identify a candidate gene, Tlcd2, for the eosinophil QTL. Our results demonstrate that hallmark allergic airway disease phenotypes are associated with distinct genetic loci on Chrs 4 and 11, and that gene expression in the allergically inflamed lung is controlled by both cis and trans regulatory factors

The Founder Strains of the Collaborative Cross Express a Complex Combination of Advantageous and Deleterious Traits for Male Reproduction

Surveys of inbred strains of mice are standard approaches to determine the heritability and range of phenotypic variation for biomedical traits. In addition, they may lead to the identification of novel phenotypes and models of human disease. Surprisingly, male reproductive phenotypes are among the least-represented traits in the Mouse Phenome Database. Here we report the results of a broad survey of the eight founder inbred strains of both the Collaborative Cross (CC) and the Diversity Outbred populations, two new mouse resources that are being used as platforms for systems genetics and sources of mouse models of human diseases. Our survey includes representatives of the three main subspecies of the house mice and a mix of classical and wild-derived inbred strains. In addition to standard staples of male reproductive phenotyping such as reproductive organ weights, sperm counts, and sperm morphology, our survey includes sperm motility and the first detailed survey of testis histology. As expected for such a broad survey, heritability varies widely among traits. We conclude that although all eight inbred strains are fertile, most display a mix of advantageous and deleterious male reproductive traits. The CAST/EiJ strain is an outlier, with an unusual combination of deleterious male reproductive traits including low sperm counts, high levels of morphologically abnormal sperm, and poor motility. In contrast, sperm from the PWK/PhJ and WSB/EiJ strains had the greatest percentages of normal morphology and vigorous motility. Finally, we report an abnormal testis phenotype that is highly heritable and restricted to the WSB/EiJ strain. This phenotype is characterized by the presence of a large, but variable, number of vacuoles in at least 10% of the seminiferous tubules. The onset of the phenotype between 2 and 3 wk of age is temporally correlated with the formation of the blood-testis barrier. We speculate that this phenotype may play a role in high rates of extinction in the CC project and in the phenotypes associated with speciation in genetic crosses that use the WSB/EiJ strain as representative of the Mus muculus domesticus subspecies

Expression Quantitative Trait Loci for Extreme Host Response to Influenza A in Pre-Collaborative Cross Mice

Outbreaks of influenza occur on a yearly basis, causing a wide range of symptoms across the human population. Although evidence exists that the host response to influenza infection is influenced by genetic differences in the host, this has not been studied in a system with genetic diversity mirroring that of the human population. Here we used mice from 44 influenza-infected pre-Collaborative Cross lines determined to have extreme phenotypes with regard to the host response to influenza A virus infection. Global transcriptome profiling identified 2671 transcripts that were significantly differentially expressed between mice that showed a severe (“high”) and mild (“low”) response to infection. Expression quantitative trait loci mapping was performed on those transcripts that were differentially expressed because of differences in host response phenotype to identify putative regulatory regions potentially controlling their expression. Twenty-one significant expression quantitative trait loci were identified, which allowed direct examination of genes associated with regulation of host response to infection. To perform initial validation of our findings, quantitative polymerase chain reaction was performed in the infected founder strains, and we were able to confirm or partially confirm more than 70% of those tested. In addition, we explored putative causal and reactive (downstream) relationships between the significantly regulated genes and others in the high or low response groups using structural equation modeling. By using systems approaches and a genetically diverse population, we were able to develop a novel framework for identifying the underlying biological subnetworks under host genetic control during influenza virus infection

Modeling Host Genetic Regulation of Influenza Pathogenesis in the Collaborative Cross

Genetic variation contributes to host responses and outcomes following infection by influenza A virus or other viral infections. Yet narrow windows of disease symptoms and confounding environmental factors have made it difficult to identify polymorphic genes that contribute to differential disease outcomes in human populations. Therefore, to control for these confounding environmental variables in a system that models the levels of genetic diversity found in outbred populations such as humans, we used incipient lines of the highly genetically diverse Collaborative Cross (CC) recombinant inbred (RI) panel (the pre-CC population) to study how genetic variation impacts influenza associated disease across a genetically diverse population. A wide range of variation in influenza disease related phenotypes including virus replication, virus-induced inflammation, and weight loss was observed. Many of the disease associated phenotypes were correlated, with viral replication and virus-induced inflammation being predictors of virus-induced weight loss. Despite these correlations, pre-CC mice with unique and novel disease phenotype combinations were observed. We also identified sets of transcripts (modules) that were correlated with aspects of disease. In order to identify how host genetic polymorphisms contribute to the observed variation in disease, we conducted quantitative trait loci (QTL) mapping. We identified several QTL contributing to specific aspects of the host response including virus-induced weight loss, titer, pulmonary edema, neutrophil recruitment to the airways, and transcriptional expression. Existing whole-genome sequence data was applied to identify high priority candidate genes within QTL regions. A key host response QTL was located at the site of the known anti-influenza Mx1 gene. We sequenced the coding regions of Mx1 in the eight CC founder strains, and identified a novel Mx1 allele that showed reduced ability to inhibit viral replication, while maintaining protection from weight loss