A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared to its predecessor. Gene annotations are now more complete, significantly improving the mapping precision of genomic, transcriptomic, and proteomics data sets. We jointly analyzed 163 short-read whole genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ~20.0 million sequence variations, of which 18.7 thousand are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Evaluation and characterization of expression quantitative trait analysis methods in the Hybrid Rat Diversity Panel.

The Hybrid Rat Diversity Panel (HRDP) is a stable and well-characterized set of more than 90 inbred rat strains that can be leveraged for systems genetics approaches to understanding the genetic and genomic variation associated with complex disease. The HRDP exhibits substantial between-strain diversity while retaining substantial within-strain isogenicity, allowing for the precise mapping of genetic variation associated with complex phenotypes and providing statistical power to identify associated variants. In order to robustly identify associated genetic variants, it is important to account for the population structure induced by inbreeding. To this end, we investigate the performance of four plausible approaches towards modeling quantitative traits in the HRDP and quantify their operating characteristics. In particular, we investigate three approaches based on genome-wide mixed model analysis, and one approach based on ordinary least squares linear regression. Towards facilitating study planning and design, we conduct extensive simulations to investigate the power of genetic association analyses in the HRDP, and characterize the impressive attained power. In simulation of eQTL data in the HRDP, we find that a mixed model approach that leverages leave-one-chromosome-out kinship estimation attains the highest power while controlling type I error.

Beyond Genes: Inclusion of Alternative Splicing and Alternative Polyadenylation to Assess the Genetic Architecture of Predisposition to Voluntary Alcohol Consumption in Brain of the HXB/BXH Recombinant Inbred Rat Panel.

Post transcriptional modifications of RNA are powerful mechanisms by which eukaryotes expand their genetic diversity. For instance, researchers estimate that most transcripts in humans undergo alternative splicing and alternative polyadenylation. These splicing events produce distinct RNA molecules, which in turn yield distinct protein isoforms and/or influence RNA stability, translation, nuclear export, and RNA/protein cellular localization. Due to their pervasiveness and impact, we hypothesized that alternative splicing and alternative polyadenylation in brain can contribute to a predisposition for voluntary alcohol consumption. Using the HXB/BXH recombinant inbred rat panel (a subset of the Hybrid Rat Diversity Panel), we generated over one terabyte of brain RNA sequencing data (total RNA) and identified novel splice variants (via StringTie) and alternative polyadenylation sites (via aptardi) to determine the transcriptional landscape in the brains of these animals. After establishing an analysis pipeline to ascertain high quality transcripts, we quantitated transcripts and integrated genotype data to identify candidate transcript coexpression networks and individual candidate transcripts associated with predisposition to voluntary alcohol consumption in the two-bottle choice paradigm. For genes that were previously associated with this trait (e.g., Lrap, Ift81, and P2rx4) (Saba et al., Febs. J., 282, 3556-3578, Saba et al., Genes. Brain. Behav., 20, e12698), we were able to distinguish between transcript variants to provide further information about the specific isoforms related to the trait. We also identified additional candidate transcripts associated with the trait of voluntary alcohol consumption (i.e., isoforms of Mapkapk5, Aldh1a7, and Map3k7). Consistent with our previous work, our results indicate that transcripts and networks related to inflammation and the immune system in brain can be linked to voluntary alcohol consumption. Overall, we have established a pipeline for including the quantitation of alternative splicing and alternative polyadenylation variants in the transcriptome in the analysis of the relationship between the transcriptome and complex traits.

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.

A long non-coding RNA (Lrap) modulates brain gene expression and levels of alcohol consumption in rats.

LncRNAs are important regulators of quantitative and qualitative features of the transcriptome. We have used QTL and other statistical analyses to identify a gene coexpression module associated with alcohol consumption. The "hub gene" of this module, Lrap (Long non-coding RNA for alcohol preference), was an unannotated transcript resembling a lncRNA. We used partial correlation analyses to establish that Lrap is a major contributor to the integrity of the coexpression module. Using CRISPR/Cas9 technology, we disrupted an exon of Lrap in Wistar rats. Measures of alcohol consumption in wild type, heterozygous and knockout rats showed that disruption of Lrap produced increases in alcohol consumption/alcohol preference. The disruption of Lrap also produced changes in expression of over 700 other transcripts. Furthermore, it became apparent that Lrap may have a function in alternative splicing of the affected transcripts. The GO category of "Response to Ethanol" emerged as one of the top candidates in an enrichment analysis of the differentially expressed transcripts. We validate the role of Lrap as a mediator of alcohol consumption by rats, and also implicate Lrap as a modifier of the expression and splicing of a large number of brain transcripts. A defined subset of these transcripts significantly impacts alcohol consumption by rats (and possibly humans). Our work shows the pleiotropic nature of non-coding elements of the genome, the power of network analysis in identifying the critical elements influencing phenotypes, and the fact that not all changes produced by genetic editing are critical for the concomitant changes in phenotype.

The multiMiR R package and database: integration of microRNA-target interactions along with their disease and drug associations.

microRNAs (miRNAs) regulate expression by promoting degradation or repressing translation of target transcripts. miRNA target sites have been catalogued in databases based on experimental validation and computational prediction using various algorithms. Several online resources provide collections of multiple databases but need to be imported into other software, such as R, for processing, tabulation, graphing and computation. Currently available miRNA target site packages in R are limited in the number of databases, types of databases and flexibility. We present multiMiR, a new miRNA-target interaction R package and database, which includes several novel features not available in existing R packages: (i) compilation of nearly 50 million records in human and mouse from 14 different databases, more than any other collection; (ii) expansion of databases to those based on disease annotation and drug microRNAresponse, in addition to many experimental and computational databases; and (iii) user-defined cutoffs for predicted binding strength to provide the most confident selection. Case studies are reported on various biomedical applications including mouse models of alcohol consumption, studies of chronic obstructive pulmonary disease in human subjects, and human cell line models of bladder cancer metastasis. We also demonstrate how multiMiR was used to generate testable hypotheses that were pursued experimentally.

Multi-step preparation technique to recover multiple metabolite compound classes for in-depth and informative metabolomic analysis.

Metabolomics is an emerging field which enables profiling of samples from living organisms in order to obtain insight into biological processes. A vital aspect of metabolomics is sample preparation whereby inconsistent techniques generate unreliable results. This technique encompasses protein precipitation, liquid-liquid extraction, and solid-phase extraction as a means of fractionating metabolites into four distinct classes. Improved enrichment of low abundance molecules with a resulting increase in sensitivity is obtained, and ultimately results in more confident identification of molecules. This technique has been applied to plasma, bronchoalveolar lavage fluid, and cerebrospinal fluid samples with volumes as low as 50 µl.  Samples can be used for multiple downstream applications; for example, the pellet resulting from protein precipitation can be stored for later analysis. The supernatant from that step undergoes liquid-liquid extraction using water and strong organic solvent to separate the hydrophilic and hydrophobic compounds. Once fractionated, the hydrophilic layer can be processed for later analysis or discarded if not needed. The hydrophobic fraction is further treated with a series of solvents during three solid-phase extraction steps to separate it into fatty acids, neutral lipids, and phospholipids. This allows the technician the flexibility to choose which class of compounds is preferred for analysis. It also aids in more reliable metabolite identification since some knowledge of chemical class exists.

Transient and persistent metabolomic changes in plasma following chronic cigarette smoke exposure in a mouse model.

Cigarette smoke exposure is linked to the development of a variety of chronic lung and systemic diseases in susceptible individuals. Metabolomics approaches may aid in defining disease phenotypes, may help predict responses to treatment, and could identify biomarkers of risk for developing disease. Using a mouse model of chronic cigarette smoke exposure sufficient to cause mild emphysema, we investigated whether cigarette smoke induces distinct metabolic profiles and determined their persistence following smoking cessation. Metabolites were extracted from plasma and fractionated based on chemical class using liquid-liquid and solid-phase extraction prior to performing liquid chromatography mass spectrometry-based metabolomics. Metabolites were evaluated for statistically significant differences among group means (p-value≤0.05) and fold change ≥1.5). Cigarette smoke exposure was associated with significant differences in amino acid, purine, lipid, fatty acid, and steroid metabolite levels compared to air exposed animals. Whereas 60% of the metabolite changes were reversible, 40% of metabolites remained persistently altered even following 2 months of smoking cessation, including nicotine metabolites. Validation of metabolite species and translation of these findings to human plasma metabolite signatures induced by cigarette smoking may lead to the discovery of biomarkers or pathogenic pathways of smoking-induced disease.

Cytosolic phospholipase A(2)α and eicosanoids regulate expression of genes in macrophages involved in host defense and inflammation.

The role of Group IVA cytosolic phospholipase A2 (cPLA2α) activation in regulating macrophage transcriptional responses to Candida albicans infection was investigated. cPLA2α releases arachidonic acid for the production of eicosanoids. In mouse resident peritoneal macrophages, prostacyclin, prostaglandin E2 and leukotriene C4 were produced within minutes of C. albicans addition before cyclooxygenase 2 expression. The production of TNFα was lower in C. albicans-stimulated cPLA2α(+/+) than cPLA2α(-/-) macrophages due to an autocrine effect of prostaglandins that increased cAMP to a greater extent in cPLA2α(+/+) than cPLA2α(-/-) macrophages. For global insight, differential gene expression in C. albicans-stimulated cPLA2α(+/+) and cPLA2α(-/-) macrophages (3 h) was compared by microarray. cPLA2α(+/+) macrophages expressed 86 genes at lower levels and 181 genes at higher levels than cPLA2α(-/-) macrophages (≥2-fold, p<0.05). Several pro-inflammatory genes were expressed at lower levels (Tnfα, Cx3cl1, Cd40, Ccl5, Csf1, Edn1, CxCr7, Irf1, Irf4, Akna, Ifnγ, several IFNγ-inducible GTPases). Genes that dampen inflammation (Socs3, Il10, Crem, Stat3, Thbd, Thbs1, Abca1) and genes involved in host defense (Gja1, Csf3, Trem1, Hdc) were expressed at higher levels in cPLA2α(+/+) macrophages. Representative genes expressed lower in cPLA2α(+/+) macrophages (Tnfα, Csf1) were increased by treatment with a prostacyclin receptor antagonist and protein kinase A inhibitor, whereas genes expressed at higher levels (Crem, Nr4a2, Il10, Csf3) were suppressed. The results suggest that C. albicans stimulates an autocrine loop in macrophages involving cPLA2α, cyclooxygenase 1-derived prostaglandins and increased cAMP that globally effects expression of genes involved in host defense and inflammation.

New sample preparation approach for mass spectrometry-based profiling of plasma results in improved coverage of metabolome.

Sample preparation remains a challenge in untargeted metabolomics studies and no method currently results in complete extraction of all metabolite classes in human plasma. Because a large variety of molecules, with vast differences in dynamic range, could be involved in human disease, there is an urgent need to develop analytical techniques that result in comprehensive coverage of metabolites. Furthermore, analysis of more focused molecular classes could be necessary to more fully interrogate markers of human disease. However, such techniques, which generally involve multiple steps, often result in high variability. We have optimized a combined liquid-liquid and solid phase extraction method for plasma and have compared that to traditional methanol precipitation using spiked internal standards as controls. This method, based largely on previously published methods, results in 5 separate fractions enriched for aqueous species, phospholipids, fatty acids, neutral lipids, and hydrophobic lipids. Using liquid chromatography mass spectrometry as the analytical method, we detect over 3806 metabolites using the new method versus 1851 metabolites using methanol alone. Qualitative analysis and quantitative analysis of both internal standards (ISTDs) and endogenous metabolites demonstrate excellent reproducibility with CV's below 15% for the combined method compared to 30% using the methanol method. While both methods have applications for clinical metabolomics, fractionation resulted in greater overall coverage and can be used for initial classification of molecular species.

Galectin-9 functionally impairs natural killer cells in humans and mice.

Galectin-9 is a pleiotropic immune modulator affecting numerous cell types of innate and adaptive immunity. Patients with chronic infection with either hepatitis C virus (HCV) or HIV have elevated circulating levels. Limited data exist on the regulation of natural killer (NK) cell function through interaction with galectin-9. We found that galectin-9 ligation downregulates multiple immune-activating genes, including eight involved in the NK cell-mediated cytotoxicity pathway, impairs lymphokine-activated killing, and decreases the proportion of gamma interferon (IFN-γ)-producing NK cells that had been stimulated with interleukin-12 (IL-12)/IL-15. We demonstrate that the transcriptional and functional changes induced by galectin-9 are independent of Tim-3. Consistent with these results for humans, we find that the genetic absence of galectin-9 in mice is associated with greater IFN-γ production by NK cells and enhanced degranulation. We also show that in the setting of a short-term (4-day) murine cytomegalovirus infection, terminally differentiated NKs accumulate in the livers of galectin-9 knockout mice, and that hepatic NKs spontaneously produce significantly more IFN-γ in this setting. Taken together, our results indicate that galectin-9 engagement impairs the function of NK cells, including cytotoxicity and cytokine production.

Comparative proteomic profiling of patients with atopic dermatitis based on history of eczema herpeticum infection and Staphylococcus aureus colonization.

BACKGROUND: Atopic dermatitis (AD) is the most common inflammatory skin disorder in the general population worldwide, and the majority of patients are colonized with Staphylococcus aureus. Eczema herpeticum is a disseminated herpes simplex virus infection that occurs in a small subset of patients. OBJECTIVES: The goal was to conduct proteomic profiling of patients with AD based on S. aureus colonization status and history of eczema herpeticum. We hoped to identify new biomarkers for improved diagnosis and prediction of eczema herpeticum and S. aureus susceptibility and to generate new hypotheses regarding disease pathogenesis. METHODS: Skin taping was performed on nonlesional skin of nonatopic control subjects and on lesional and nonlesional skin of patients with AD. Subjects were classified according to the history of eczema herpeticum and S. aureus colonization. Proteins were analyzed by using mass spectrometry; diagnostic groups were compared for statistically significant differences in protein expression. RESULTS: Proteins related to the skin barrier (filaggrin-2, corneodesmosin, desmoglein-1, desmocollin-1, and transglutaminase-3) and generation of natural moisturizing factor (arginase-1, caspase-14, and gamma-glutamyl cyclotransferase) were expressed at significantly lower levels in lesional versus nonlesional sites of patients with AD with and without history of eczema herpeticum; epidermal fatty acid-binding protein was expressed at significantly higher levels in patients with methicillin-resistant S. aureus. CONCLUSION: This noninvasive, semiquantitative profiling method has revealed novel proteins likely involved in the pathogenesis of AD. The lower expression of skin barrier proteins and enzymes involved in the generation of the natural moisturizing factor could further exacerbate barrier defects and perpetuate water loss from the skin. The greater expression of epidermal fatty acid-binding protein, especially in patients colonized with methicillin-resistant S. aureus, might perpetuate the inflammatory response through eicosanoid signaling.

Proteomic definition of the cell wall of Mycobacterium tuberculosis.

The cell envelope of Mycobacterium tuberculosis (Mtb) is complex and diverse; composed of proteins intermingled in a matrix of peptidoglycan, mycolic acids, lipids, and carbohydrates. Proteomic studies of the Mtb cell wall have been limited; nonetheless, the characterization of resident and secreted proteins associated with the cell wall are critical to understanding bacterial survival and immune modulation in the host. In this study, the cell wall proteome was defined in order to better understand its unique biosynthetic and secretion processes. Mtb cell wall was subjected to extraction with organic solvents to remove noncovalently bound lipids and lipoglycans and remaining proteins were solubilized with either SDS, Guanidine-HCl, or TX-114. These extracts were analyzed by two-dimensional gel electrophoresis and mass-spectrometry and resulted in the identification of 234 total proteins. The lipoproteome of Mtb, enriched in the TX-114 extract, was further resolved by multidimensional chromatography and mass spectrometry to identify an additional 294 proteins. A query of the 528 total protein identifications against Neural Network or Hidden Markov model algorithms predicted secretion signals in 87 proteins. Classification of these 528 proteins also demonstrated that 35% are involved in small molecule metabolism and 25% are involved in macromolecule synthesis and degradation building upon evidence that the Mtb cell wall is actively engaged in mycobacterial survival and remodeling.

Isolation of mycobacterium species genomic DNA.

A myriad of methods has been reported for the isolation of genomic DNA from Mycobacterium spp.; some methods use mechanical disruption of the bacterial cells, whereas others use some form of chemical or enzymatic lysis. Regardless of the approach, the end points remain efficient breaking of the complex mycobacterial cell wall and release of high-quality DNA that is suitable for manipulation and analyses by molecular genetic techniques. This chapter providers detailed methods for the large and small isolation of mycobacterial genomic DNA.

Mycobacterium tuberculosis malate synthase is a laminin-binding adhesin.

Mycobacterium tuberculosis (M. tb) uses the glyoxalate bypass for intracellular survival in vivo. These studies provide evidence that the M. tb malate synthase (MS) has adapted to function as an adhesin which binds to laminin and fibronectin. This binding is achieved via the unique C-terminal region of the M. tb MS. The ability to function as an adhesin necessitates extracellular localization. We provide evidence that despite the absence of a Sec-translocation signal sequence the M. tb MS is secreted/excreted, and is anchored on the cell wall by an undefined mechanism. The MS of Mycobacterium smegmatis is cytoplasmic but the M. tb MS expressed in M. smegmatis localizes to the cell wall and enhances the adherence of the bacteria to lung epithelial A549 cells. Antibodies to the C-terminal laminin/fibronectin-binding domain interfere with the binding of the M. tb MS to laminin and fibronectin and reduce the adherence of M. tb to A549 cells. Coupled to the earlier evidence of in vivo expression of M. tb MS during active but not latent infection in humans, these studies show that a housekeeping enzyme of M. tb contributes to its armamentarium of virulence promoting factors.