Facilitating Complex Trait Analysis via Reduced Complexity Crosses.

Genetically diverse inbred strains are frequently used in quantitative trait mapping to identify sequence variants underlying trait variation. Poor locus resolution and high genetic complexity impede variant discovery. As a solution, we explore reduced complexity crosses (RCCs) between phenotypically divergent, yet genetically similar, rodent substrains. RCCs accelerate functional variant discovery via decreasing the number of segregating variants by orders of magnitude. The simplified genetic architecture of RCCs often permit immediate identification of causal variants or rapid fine-mapping of broad loci to smaller intervals. Whole-genome sequences of substrains make RCCs possible by supporting the development of array- and targeted sequencing-based genotyping platforms, coupled with rapid genome editing for variant validation. In summary, RCCs enhance discovery-based genetics of complex traits.

Pooled analysis of radiation hybrids identifies loci for growth and drug action in mammalian cells.

Genetic screens in mammalian cells commonly focus on loss-of-function approaches. To evaluate the phenotypic consequences of extra gene copies, we used bulk segregant analysis (BSA) of radiation hybrid (RH) cells. We constructed six pools of RH cells, each consisting of ∼2500 independent clones, and placed the pools under selection in media with or without paclitaxel. Low pass sequencing identified 859 growth loci, 38 paclitaxel loci, 62 interaction loci, and three loci for mitochondrial abundance at genome-wide significance. Resolution was measured as ∼30 kb, close to single-gene. Divergent properties were displayed by the RH-BSA growth genes compared to those from loss-of-function screens, refuting the balance hypothesis. In addition, enhanced retention of human centromeres in the RH pools suggests a new approach to functional dissection of these chromosomal elements. Pooled analysis of RH cells showed high power and resolution and should be a useful addition to the mammalian genetic toolkit.

VHL L169P Variant Does Not Alter Cellular Hypoxia Tension in Clear Cell Renal Cell Carcinoma.

In the current era of tumor genome sequencing, single amino acid missense variants in the von Hippel-Lindau (VHL) tumor suppressor gene are frequently identified in clear cell renal carcinoma (ccRCC). Due to the incomplete knowledge of the structural architecture of VHL protein, the functional significance of many missense mutations cannot be assigned. L169P is one such missense mutation identified in the case of aggressive, metastatic ccRCC. Here, we characterized the biochemical activity, transcriptomic hypoxia signature and biological functions of the L169P variant. Lentiviral vector expressing either wildtype (WT) or L169P VHL were used to transduce two VHL-deficient human ccRCC cell lines, 786-O and RCC4. The stability of the VHL protein and the expression level of VHL, HIF1α and HIF2α were analyzed. The impact of restoring L169P or WT VHL on the hypoxia gene expression program in 786-O cells was assessed by mRNA sequencing (RNAseq) and computed hypoxic scores. The impact of restoring VHL expression on the growth of ccRCC models was assessed in cell cultures and in chorioallantoic membrane (CAM) xenografts. In the 786-O cells, the protein stability of L169P VHL was comparable to WT VHL. No obvious difference in the capability of degrading HIF1α and HIF2α was observed between WT and L169P VHL in the 786-O or RCC4 cells. The hypoxic scores were not significantly different in the 786-O cells expressing either wildtype or L169P VHL. From the cellular function perspective, both WT and L169P VHL slowed cell proliferation in vitro and in vivo. The L169P VHL variant is comparable to WT VHL in terms of protein stability, ability to degrade HIF1α factors and ability to regulate hypoxia gene expression, as well as in the suppression of ccRCC tumor cell growth. Taken together, our data indicate that the L169P VHL variant alone is unlikely to drive the oncogenesis of sporadic ccRCC.

Extreme phenotypic diversity in operant response to intravenous cocaine or saline infusion in the hybrid mouse diversity panel.

Cocaine self-administration is a complexly determined trait, with a substantial proportion of individual differences being determined by genetic variation. However, the relevant genetic variants that drive heritable differences in cocaine use remain undiscovered. Cocaine intravenous self-administration (IVSA) procedures in laboratory animals provide opportunities to prospectively investigate neurogenetic influences on the acquisition of voluntary cocaine use. Here, we provide information on cocaine (or saline-as a control) IVSA in 84 members of the hybrid mouse diversity panel (HMDP), an array of genetically distinct classical or recombinant inbred strains. We found cocaine IVSA to be substantially heritable in this population, with strain-level intake ranging for near 0 to >25 mg/kg/session. Though saline IVSA was also found to be heritable, a modest genetic correlation between cocaine and saline IVSA indicates that operant responding for the cocaine reinforcer was influenced, at least in part, by unique genetic variants. Genome-wide association studies (GWAS) of infusions earned in cocaine and saline groups revealed significant quantitative trait loci (QTL) on Chromosomes 3 and 14 for cocaine, but not saline, IVSA. Positional candidates were further prioritized through use of bulk RNA sequencing data that revealed genes with cis-eQTL and genetic correlation to number of infusions. Additionally, these data identify reference strains with extreme cocaine IVSA phenotypes, revealing them as polygenic models of risk and resilience to cocaine reinforcement. This work is part of an ongoing effort to characterize genetic variation that moderates cocaine IVSA that may, in turn, provide a more comprehensive understanding of cocaine risk genetics and neurobiology.

Reduced Folate Carrier: an Entry Receptor for a Novel Feline Leukemia Virus Variant.

Feline leukemia virus (FeLV) is horizontally transmitted among cats and causes a variety of hematopoietic disorders. Five subgroups of FeLV, A to D and T, each with distinct receptor usages, have been described. Recently, we identified a new FeLV Env (TG35-2) gene from a pseudotyped virus that does not belong to any known subgroup. FeLV-A is the primary virus from which other subgroups have emerged via mutation or recombination of the subgroup A env gene. Retrovirus entry into cells is mediated by the interaction of envelope protein (Env) with specific cell surface receptors. Here, phenotypic screening of a human/hamster radiation hybrid panel identified SLC19A1, a feline reduced folate carrier (RFC) and potential receptor for TG35-2-phenotypic virus. RFC is a multipass transmembrane protein. Feline and human RFC cDNAs conferred susceptibility to TG35-2-pseudotyped virus when introduced into nonpermissive cells but did not render these cells permissive to other FeLV subgroups or feline endogenous retrovirus. Moreover, human cells with genomic deletion of RFC were nonpermissive for TG35-2-pseudotyped virus infection, but the introduction of feline and human cDNAs rendered them permissive. Mutation analysis of FeLV Env demonstrated that amino acid substitutions within variable region A altered the specificity of the Env-receptor interaction. We isolated and reconstructed the full-length infectious TG35-2-phenotypic provirus from a naturally FeLV-infected cat, from which the FeLV Env (TG35-2) gene was previously isolated, and compared the replication of the virus in hematopoietic cell lines with that of FeLV-A 61E by measuring the viral RNA copy numbers. These results provide a tool for further investigation of FeLV infectious disease.IMPORTANCE Feline leukemia virus (FeLV) is a member of the genus Gammaretrovirus, which causes malignant diseases in cats. The most prevalent FeLV among cats is FeLV subgroup A (FeLV-A), and specific binding of FeLV-A Env to its viral receptor, thiamine transporter feTHTR1, is the first step of infection. In infected cats, novel variants of FeLV with altered receptor specificity for viral entry have emerged by mutation or recombination of the env gene. A novel FeLV variant arose from a subtle mutation of FeLV-A Env, which altered the specific interaction of the virus with its receptor. RFC, a folate transporter, is a potential receptor for the novel FeLV variant. The perturbation of specific retrovirus-receptor interactions under selective pressure by the host results in the emergence of novel viruses.

The Hybrid Mouse Diversity Panel: a resource for systems genetics analyses of metabolic and cardiovascular traits.

The Hybrid Mouse Diversity Panel (HMDP) is a collection of approximately 100 well-characterized inbred strains of mice that can be used to analyze the genetic and environmental factors underlying complex traits. While not nearly as powerful for mapping genetic loci contributing to the traits as human genome-wide association studies, it has some important advantages. First, environmental factors can be controlled. Second, relevant tissues are accessible for global molecular phenotyping. Finally, because inbred strains are renewable, results from separate studies can be integrated. Thus far, the HMDP has been studied for traits relevant to obesity, diabetes, atherosclerosis, osteoporosis, heart failure, immune regulation, fatty liver disease, and host-gut microbiota interactions. High-throughput technologies have been used to examine the genomes, epigenomes, transcriptomes, proteomes, metabolomes, and microbiomes of the mice under various environmental conditions. All of the published data are available and can be readily used to formulate hypotheses about genes, pathways and interactions.

Comparative analysis of proteome and transcriptome variation in mouse.

The relationships between the levels of transcripts and the levels of the proteins they encode have not been examined comprehensively in mammals, although previous work in plants and yeast suggest a surprisingly modest correlation. We have examined this issue using a genetic approach in which natural variations were used to perturb both transcript levels and protein levels among inbred strains of mice. We quantified over 5,000 peptides and over 22,000 transcripts in livers of 97 inbred and recombinant inbred strains and focused on the 7,185 most heritable transcripts and 486 most reliable proteins. The transcript levels were quantified by microarray analysis in three replicates and the proteins were quantified by Liquid Chromatography-Mass Spectrometry using O(18)-reference-based isotope labeling approach. We show that the levels of transcripts and proteins correlate significantly for only about half of the genes tested, with an average correlation of 0.27, and the correlations of transcripts and proteins varied depending on the cellular location and biological function of the gene. We examined technical and biological factors that could contribute to the modest correlation. For example, differential splicing clearly affects the analyses for certain genes; but, based on deep sequencing, this does not substantially contribute to the overall estimate of the correlation. We also employed genome-wide association analyses to map loci controlling both transcript and protein levels. Surprisingly, little overlap was observed between the protein- and transcript-mapped loci. We have typed numerous clinically relevant traits among the strains, including adiposity, lipoprotein levels, and tissue parameters. Using correlation analysis, we found that a low number of clinical trait relationships are preserved between the protein and mRNA gene products and that the majority of such relationships are specific to either the protein levels or transcript levels. Surprisingly, transcript levels were more strongly correlated with clinical traits than protein levels. In light of the widespread use of high-throughput technologies in both clinical and basic research, the results presented have practical as well as basic implications.

A genome-wide map of human genetic interactions inferred from radiation hybrid genotypes.

Using radiation hybrid genotyping data, 99% of all possible gene pairs across the mammalian genome were tested for interactions based on co-retention frequencies higher (attraction) or lower (repulsion) than chance. Gene interaction networks constructed from six independent data sets overlapped strongly. Combining the data sets resulted in a network of more than seven million interactions, almost all attractive. This network overlapped with protein-protein interaction networks on multiple measures and also confirmed the relationship between essentiality and centrality. In contrast to other biological networks, the radiation hybrid network did not show a scale-free distribution of connectivity but was Gaussian-like, suggesting a closer approach to saturation. The radiation hybrid (RH) network constitutes a platform for understanding the systems biology of the mammalian cell.

Genetic pathways regulating the longitudinal acquisition of cocaine self-administration in a panel of inbred and recombinant inbred mice.

To identify addiction genes, we evaluate intravenous self-administration of cocaine or saline in 84 inbred and recombinant inbred mouse strains over 10 days. We integrate the behavior data with brain RNA-seq data from 41 strains. The self-administration of cocaine and that of saline are genetically distinct. We maximize power to map loci for cocaine intake by using a linear mixed model to account for this longitudinal phenotype while correcting for population structure. A total of 15 unique significant loci are identified in the genome-wide association study. A transcriptome-wide association study highlights the Trpv2 ion channel as a key locus for cocaine self-administration as well as identifying 17 additional genes, including Arhgef26, Slc18b1, and Slco5a1. We find numerous instances where alternate splice site selection or RNA editing altered transcript abundance. Our work emphasizes the importance of Trpv2, an ionotropic cannabinoid receptor, for the response to cocaine.

Learning pair-wise gene functional similarity by multiplex gene expression maps.

BACKGROUND: The relationships between the gene functional similarity and gene expression profile, and between gene function annotation and gene sequence have been studied extensively. However, not much work has considered the connection between gene functions and location of a gene's expression in the mammalian tissues. On the other hand, although unsupervised learning methods have been commonly used in functional genomics, supervised learning cannot be directly applied to a set of normal genes without having a target (class) attribute. RESULTS: Here, we propose a supervised learning methodology to predict pair-wise gene functional similarity from multiplex gene expression maps that provide information about the location of gene expression. The features are extracted from expression maps and the labels denote the functional similarities of pairs of genes. We make use of wavelet features, original expression values, difference and average values of neighboring voxels and other features to perform boosting analysis. The experimental results show that with increasing similarities of gene expression maps, the functional similarities are increased too. The model predicts the functional similarities between genes to a certain degree. The weights of the features in the model indicate the features that are more significant for this prediction. CONCLUSIONS: By considering pairs of genes, we propose a supervised learning methodology to predict pair-wise gene functional similarity from multiplex gene expression maps. We also explore the relationship between similarities of gene maps and gene functions. By using AdaBoost coupled with our proposed weak classifier we analyze a large-scale gene expression dataset and predict gene functional similarities. We also detect the most significant single voxels and pairs of neighboring voxels and visualize them in the expression map image of a mouse brain. This work is very important for predicting functions of unknown genes. It also has broader applicability since the methodology can be applied to analyze any large-scale dataset without a target attribute and is not restricted to gene expressions.

Hybrid mouse diversity panel: a panel of inbred mouse strains suitable for analysis of complex genetic traits.

We have developed an association-based approach using classical inbred strains of mice in which we correct for population structure, which is very extensive in mice, using an efficient mixed-model algorithm. Our approach includes inbred parental strains as well as recombinant inbred strains in order to capture loci with effect sizes typical of complex traits in mice (in the range of 5% of total trait variance). Over the last few years, we have typed the hybrid mouse diversity panel (HMDP) strains for a variety of clinical traits as well as intermediate phenotypes and have shown that the HMDP has sufficient power to map genes for highly complex traits with resolution that is in most cases less than a megabase. In this essay, we review our experience with the HMDP, describe various ongoing projects, and discuss how the HMDP may fit into the larger picture of common diseases and different approaches.

Endogenous 3,4-dihydroxyphenylalanine and dopaquinone modifications on protein tyrosine: links to mitochondrially derived oxidative stress via hydroxyl radical.

Oxidative modifications of protein tyrosines have been implicated in multiple human diseases. Among these modifications, elevations in levels of 3,4-dihydroxyphenylalanine (DOPA), a major product of hydroxyl radical addition to tyrosine, has been observed in a number of pathologies. Here we report the first proteome survey of endogenous site-specific modifications, i.e. DOPA and its further oxidation product dopaquinone in mouse brain and heart tissues. Results from LC-MS/MS analyses included 50 and 14 DOPA-modified tyrosine sites identified from brain and heart, respectively, whereas only a few nitrotyrosine-containing peptides, a more commonly studied marker of oxidative stress, were detectable, suggesting the much higher abundance for DOPA modification as compared with tyrosine nitration. Moreover, 20 and 12 dopaquinone-modified peptides were observed from brain and heart, respectively; nearly one-fourth of these peptides were also observed with DOPA modification on the same sites. For both tissues, these modifications are preferentially found in mitochondrial proteins with metal binding properties, consistent with metal-catalyzed hydroxyl radical formation from mitochondrial superoxide and hydrogen peroxide. These modifications also link to a number of mitochondrially associated and other signaling pathways. Furthermore, many of the modification sites were common sites of previously reported tyrosine phosphorylation, suggesting potential disruption of signaling pathways. Collectively, the results suggest that these modifications are linked with mitochondrially derived oxidative stress and may serve as sensitive markers for disease pathologies.

Single-cell analysis of gene expression in the substantia nigra pars compacta of a pesticide-induced mouse model of Parkinson's disease.

Exposure to pesticides in humans increases the risk of Parkinson's disease (PD), but the mechanisms remain poorly understood. To elucidate these pathways, we dosed C57BL/6J mice with a combination of the pesticides maneb and paraquat. Behavioral analysis revealed motor deficits consistent with PD. Single-cell RNA sequencing of substantia nigra pars compacta revealed both cell-type-specific genes and genes expressed differentially between pesticide and control, including Fam241b, Emx2os, Bivm, Gm1439, Prdm15, and Rai2. Neurons had the largest number of significant differentially expressed genes, but comparable numbers were found in astrocytes and less so in oligodendrocytes. In addition, network analysis revealed enrichment in functions related to the extracellular matrix. These findings emphasize the importance of support cells in pesticide-induced PD and refocus our attention away from neurons as the sole agent of this disorder.

Effects of genome-wide copy number variation on expression in mammalian cells.

BACKGROUND: There is only a limited understanding of the relation between copy number and expression for mammalian genes. We fine mapped cis and trans regulatory loci due to copy number change for essentially all genes using a human-hamster radiation hybrid (RH) panel. These loci are called copy number expression quantitative trait loci (ceQTLs). RESULTS: Unexpected findings from a previous study of a mouse-hamster RH panel were replicated. These findings included decreased expression as a result of increased copy number for 30% of genes and an attenuated relationship between expression and copy number on the X chromosome suggesting an Xist independent form of dosage compensation. In a separate glioblastoma dataset, we found conservation of genes in which dosage was negatively correlated with gene expression. These genes were enriched in signaling and receptor activities. The observation of attenuated X-linked gene expression in response to increased gene number was also replicated in the glioblastoma dataset. Of 523 gene deserts of size > 600 kb in the human RH panel, 325 contained trans ceQTLs with -log10 P > 4.1. Recently discovered genes, ultra conserved regions, noncoding RNAs and microRNAs explained only a small fraction of the results, suggesting a substantial portion of gene deserts harbor as yet unidentified functional elements. CONCLUSION: Radiation hybrids are a useful tool for high resolution mapping of cis and trans loci capable of affecting gene expression due to copy number change. Analysis of two independent radiation hybrid panels show agreement in their findings and may serve as a discovery source for novel regulatory loci in noncoding regions of the genome.

A Poisson model for random multigraphs.

MOTIVATION: Biological networks are often modeled by random graphs. A better modeling vehicle is a multigraph where each pair of nodes is connected by a Poisson number of edges. In the current model, the mean number of edges equals the product of two propensities, one for each node. In this context it is possible to construct a simple and effective algorithm for rapid maximum likelihood estimation of all propensities. Given estimated propensities, it is then possible to test statistically for functionally connected nodes that show an excess of observed edges over expected edges. The model extends readily to directed multigraphs. Here, propensities are replaced by outgoing and incoming propensities. RESULTS: The theory is applied to real data on neuronal connections, interacting genes in radiation hybrids, interacting proteins in a literature curated database, and letter and word pairs in seven Shaskespearean plays. AVAILABILITY: All data used are fully available online from their respective sites. Source code and software is available from http://code.google.com/p/poisson-multigraph/.

Mapping protein abundance patterns in the brain using voxelation combined with liquid chromatography and mass spectrometry.

Voxelation creates expression atlases by high-throughput analysis of spatially registered cubes or voxels harvested from the brain. The modality independence of voxelation allows a variety of bioanalytical techniques to be used to map abundance. Protein expression patterns in the brain can be obtained using liquid chromatography (LC) combined with mass spectrometry (MS). Here we describe the methodology of voxelation as it pertains particularly to LC-MS proteomic analysis: sample preparation, instrumental set up and analysis, peptide identification and protein relative abundance quantitation. We also briefly describe some of the advantages, limitations and insights into the brain that can be obtained using combined proteomic and transcriptomic maps.

High-resolution mapping of gene expression using association in an outbred mouse stock.

Quantitative trait locus (QTL) analysis is a powerful tool for mapping genes for complex traits in mice, but its utility is limited by poor resolution. A promising mapping approach is association analysis in outbred stocks or different inbred strains. As a proof of concept for the association approach, we applied whole-genome association analysis to hepatic gene expression traits in an outbred mouse population, the MF1 stock, and replicated expression QTL (eQTL) identified in previous studies of F2 intercross mice. We found that the mapping resolution of these eQTL was significantly greater in the outbred population. Through an example, we also showed how this precise mapping can be used to resolve previously identified loci (in intercross studies), which affect many different transcript levels (known as eQTL "hotspots"), into distinct regions. Our results also highlight the importance of correcting for population structure in whole-genome association studies in the outbred stock.

Cost-Effective Mapping of Genetic Interactions in Mammalian Cells.

Comprehensive maps of genetic interactions in mammalian cells are daunting to construct because of the large number of potential interactions, ~ 2 × 108 for protein coding genes. We previously used co-inheritance of distant genes from published radiation hybrid (RH) datasets to identify genetic interactions. However, it was necessary to combine six legacy datasets from four species to obtain adequate statistical power. Mapping resolution was also limited by the low density PCR genotyping. Here, we employ shallow sequencing of nascent human RH clones as an economical approach to constructing interaction maps. In this initial study, 15 clones were analyzed, enabling construction of a network with 225 genes and 2,359 interactions (FDR < 0.05). Despite its small size, the network showed significant overlap with the previous RH network and with a protein-protein interaction network. Consumables were ≲$50 per clone, showing that affordable, high quality genetic interaction maps are feasible in mammalian cells.

Integration of evidence across human and model organism studies: A meeting report.

The National Institute on Drug Abuse and Joint Institute for Biological Sciences at the Oak Ridge National Laboratory hosted a meeting attended by a diverse group of scientists with expertise in substance use disorders (SUDs), computational biology, and FAIR (Findability, Accessibility, Interoperability, and Reusability) data sharing. The meeting's objective was to discuss and evaluate better strategies to integrate genetic, epigenetic, and 'omics data across human and model organisms to achieve deeper mechanistic insight into SUDs. Specific topics were to (a) evaluate the current state of substance use genetics and genomics research and fundamental gaps, (b) identify opportunities and challenges of integration and sharing across species and data types, (c) identify current tools and resources for integration of genetic, epigenetic, and phenotypic data, (d) discuss steps and impediment related to data integration, and (e) outline future steps to support more effective collaboration-particularly between animal model research communities and human genetics and clinical research teams. This review summarizes key facets of this catalytic discussion with a focus on new opportunities and gaps in resources and knowledge on SUDs.

Differential expression and function of ABCG1 and ABCG4 during development and aging.

ABCG1 and ABCG4 are highly homologous members of the ATP binding cassette (ABC) transporter family that regulate cellular cholesterol homeostasis. In adult mice, ABCG1 is known to be expressed in numerous cell types and tissues, whereas ABCG4 expression is limited to the central nervous system (CNS). Here, we show significant differences in expression of these two transporters during development. Examination of beta-galactosidase-stained tissue sections from Abcg1(-/-)LacZ and Abcg4(-/-)LacZ knockin mice shows that ABCG4 is highly but transiently expressed both in hematopoietic cells and in enterocytes during development. In contrast, ABCG1 is expressed in macrophages and in endothelial cells of both embryonic and adult liver. We also show that ABCG1 and ABCG4 are both expressed as early as E12.5 in the embryonic eye and developing CNS. Loss of both ABCG1 and ABCG4 results in accumulation in the retina and/or brain of oxysterols, in altered expression of liver X receptor and sterol-regulatory element binding protein-2 target genes, and in a stress response gene. Finally, behavioral tests show that Abcg4(-/-) mice have a general deficit in associative fear memory. Together, these data indicate that loss of ABCG1 and/or ABCG4 from the CNS results in changes in metabolic pathways and in behavior.