Physiogenomic resources for rat models of heart, lung and blood disorders.

Cardiovascular disorders are influenced by genetic and environmental factors. The TIGR rodent expression web-based resource (TREX) contains over 2,200 microarray hybridizations, involving over 800 animals from 18 different rat strains. These strains comprise genetically diverse parental animals and a panel of chromosomal substitution strains derived by introgressing individual chromosomes from normotensive Brown Norway (BN/NHsdMcwi) rats into the background of Dahl salt sensitive (SS/JrHsdMcwi) rats. The profiles document gene-expression changes in both genders, four tissues (heart, lung, liver, kidney) and two environmental conditions (normoxia, hypoxia). This translates into almost 400 high-quality direct comparisons (not including replicates) and over 100,000 pairwise comparisons. As each individual chromosomal substitution strain represents on average less than a 5% change from the parental genome, consomic strains provide a useful mechanism to dissect complex traits and identify causative genes. We performed a variety of data-mining manipulations on the profiles and used complementary physiological data from the PhysGen resource to demonstrate how TREX can be used by the cardiovascular community for hypothesis generation.

Comparison of three next-generation sequencing platforms for metagenomic sequencing and identification of pathogens in blood.

BACKGROUND: The introduction of benchtop sequencers has made adoption of whole genome sequencing possible for a broader community of researchers than ever before. Concurrently, metagenomic sequencing (MGS) is rapidly emerging as a tool for interrogating complex samples that defy conventional analyses. In addition, next-generation sequencers are increasingly being used in clinical or related settings, for instance to track outbreaks. However, information regarding the analytical sensitivity or limit of detection (LoD) of benchtop sequencers is currently lacking. Furthermore, the specificity of sequence information at or near the LoD is unknown. RESULTS: In the present study, we assess the ability of three next-generation sequencing platforms to identify a pathogen (viral or bacterial) present in low titers in a clinically relevant sample (blood). Our results indicate that the Roche-454 Titanium platform is capable of detecting Dengue virus at titers as low as 1X102.5 pfu/mL, corresponding to an estimated 5.4X104 genome copies/ml maximum. The increased throughput of the benchtop sequencers, the Ion Torrent PGM and Illumina MiSeq platforms, enabled detection of viral genomes at concentrations as low as 1X104 genome copies/mL. Platform-specific biases were evident in sequence read distributions as well as viral genome coverage. For bacterial samples, only the MiSeq platform was able to provide sequencing reads that could be unambiguously classified as originating from Bacillus anthracis. CONCLUSION: The analytical sensitivity of all three platforms approaches that of standard qPCR assays. Although all platforms were able to detect pathogens at the levels tested, there were several noteworthy differences. The Roche-454 Titanium platform produced consistently longer reads, even when compared with the latest chemistry updates for the PGM platform. The MiSeq platform produced consistently greater depth and breadth of coverage, while the Ion Torrent was unequaled for speed of sequencing. None of the platforms were able to verify a single nucleotide polymorphism responsible for antiviral resistance in an Influenza A strain isolated from the 2009 H1N1 pandemic. Overall, the benchtop platforms perform well for identification of pathogens from a representative clinical sample. However, unlike identification, characterization of pathogens is likely to require higher titers, multiple libraries and/or multiple sequencing runs.

Neuroplasticity, axonal guidance and micro-RNA genes are associated with morphine self-administration behavior.

Neuroadaptations in the ventral striatum (VS) and ventral midbrain (VMB) following chronic opioid administration are thought to contribute to the pathogenesis and persistence of opiate addiction. In order to identify candidate genes involved in these neuroadaptations, we utilized a behavior-genetics strategy designed to associate contingent intravenous drug self-administration with specific patterns of gene expression in inbred mice differentially predisposed to the rewarding effects of morphine. In a Yoked-control paradigm, C57BL/6J mice showed clear morphine-reinforced behavior, whereas DBA/2J mice did not. Moreover, the Yoked-control paradigm revealed the powerful consequences of self-administration versus passive administration at the level of gene expression. Morphine self-administration in the C57BL/6J mice uniquely up- or down-regulated 237 genes in the VS and 131 genes in the VMB. Interestingly, only a handful of the C57BL/6J self-administration genes (<3%) exhibited a similar expression pattern in the DBA/2J mice. Hence, specific sets of genes could be confidently assigned to regional effects of morphine in a contingent- and genotype-dependent manner. Bioinformatics analysis revealed that neuroplasticity, axonal guidance and micro-RNAs (miRNAs) were among the key themes associated with drug self-administration. Noteworthy were the primary miRNA genes H19 and micro-RNA containing gene (Mirg), processed, respectively, to mature miRNAs miR-675 and miR-154, because they are prime candidates to mediate network-like changes in responses to chronic drug administration. These miRNAs have postulated roles in dopaminergic neuron differentiation and mu-opioid receptor regulation. The strategic approach designed to focus on reinforcement-associated genes provides new insight into the role of neuroplasticity pathways and miRNAs in drug addiction.

Positional identification of variants of Adamts16 linked to inherited hypertension.

A previously reported blood pressure (BP) quantitative trait locus on rat Chromosome 1 was isolated in a short congenic segment spanning 804.6 kb. The 804.6 kb region contained only two genes, LOC306664 and LOC306665. LOC306664 is predicted to translate into A Disintegrin-like and Metalloproteinase with Thrombospondin Motifs-16 (Adamts16). LOC306665 is a novel gene. All predicted exons of both LOC306664 and LOC306665 were sequenced. Non-synonymous variants were identified in only one of these genes, LOC306664. These variants were naturally existing polymorphisms among inbred, outbred and wild rats. The full-length rat transcript of Adamts16 was detected in multiple tissues. Similar to ADAMTS16 in humans, expression of Adamts16 was prominent in the kidney. Renal transcriptome analysis suggested that a network of genes related to BP was differential between congenic and S rats. These genes were also differentially expressed between kidney cell lines with or without knock-down of Adamts16. Adamts16 is conserved between rats and humans. It is a candidate gene within the homologous region on human Chromosome 5, which is linked to systolic and diastolic BP in the Quebec Family Study. Multiple variants, including an Ala to Pro variant in codon 90 (rs2086310) of human ADAMTS16, were associated with human resting systolic BP (SBP). Replication study in GenNet confirmed the association of two variants of ADAMTS16 with SBP, including rs2086310. Overall, our report represents a high resolution positional cloning and translational study for Adamts16 as a candidate gene controlling BP.

TLX1 and NOTCH coregulate transcription in T cell acute lymphoblastic leukemia cells.

BACKGROUND: The homeobox gene TLX1 (for T-cell leukemia homeobox 1, previously known as HOX11) is inappropriately expressed in a major subgroup of T cell acute lymphoblastic leukemia (T-ALL) where it is strongly associated with activating NOTCH1 mutations. Despite the recognition that these genetic lesions cooperate in leukemogenesis, there have been no mechanistic studies addressing how TLX1 and NOTCH1 functionally interact to promote the leukemic phenotype. RESULTS: Global gene expression profiling after downregulation of TLX1 and inhibition of the NOTCH pathway in ALL-SIL cells revealed that TLX1 synergistically regulated more than 60% of the NOTCH-responsive genes. Structure-function analysis demonstrated that TLX1 binding to Groucho-related TLE corepressors was necessary for maximal transcriptional regulation of the NOTCH-responsive genes tested, implicating TLX1 modulation of the NOTCH-TLE regulatory network. Comparison of the dataset to publicly available biological databases indicated that the TLX1/NOTCH-coregulated genes are frequently targeted by MYC. Gain- and loss-of-function experiments confirmed that MYC was an essential mediator of TLX1/NOTCH transcriptional output and growth promotion in ALL-SIL cells, with TLX1 contributing to the NOTCH-MYC regulatory axis by posttranscriptional enhancement of MYC protein levels. Functional classification of the TLX1/NOTCH-coregulated targets also showed enrichment for genes associated with other human cancers as well as those involved in developmental processes. In particular, we found that TLX1, NOTCH and MYC coregulate CD1B and RAG1, characteristic markers of early cortical thymocytes, and that concerted downregulation of the TLX1 and NOTCH pathways resulted in their irreversible repression. CONCLUSIONS: We found that TLX1 and NOTCH synergistically regulate transcription in T-ALL, at least in part via the sharing of a TLE corepressor and by augmenting expression of MYC. We conclude that the TLX1/NOTCH/MYC network is a central determinant promoting the growth and survival of TLX1+ T-ALL cells. In addition, the TLX1/NOTCH/MYC transcriptional network coregulates genes involved in T cell development, such as CD1 and RAG family members, and therefore may prescribe the early cortical stage of differentiation arrest characteristic of the TLX1 subgroup of T-ALL.

Molecular networks in Dahl salt-sensitive hypertension based on transcriptome analysis of a panel of consomic rats.

The Dahl salt-sensitive (SS) rat is a widely used model of human salt-sensitive hypertension and renal injury. We studied the molecular networks that underlie the complex disease phenotypes in the SS model, using a design that involved two consomic rat strains that were protected from salt-induced hypertension and one that was not protected. Substitution of Brown Norway (BN) chromosome 13 or 18, but not 20, into the SS genome was found to significantly attenuate salt-induced hypertension and albuminuria. Gene expression profiles were examined in the kidneys of SS and consomic SS-13(BN), SS-18(BN), and SS-20(BN) rats with a total of 240 cDNA microarrays. The substituted chromosome was overrepresented in genes differentially expressed between a consomic strain and SS rats on a 0.4% salt diet. F5, Serpinc1, Slc19a2, and genes represented by three other expressed sequence tags (ESTs), which are located on chromosome 13, were found to be differentially expressed between SS-13(BN) and all other strains examined. Likewise, Acaa2, B4galt6, Colec12, Hsd17b4, and five other ESTs located on chromosome 18 exhibited expression patterns unique to SS-18(BN). On exposure to a 4% salt diet, there were 184 ESTs in the renal cortex and 346 in the renal medulla for which SS-13(BN) and SS-18(BN) shared one expression pattern, while SS and SS-20(BN) shared another, mirroring the phenotypic segregation among the four strains. Molecular networks that might contribute to the development of Dahl salt-sensitive hypertension and albuminuria were constructed with an approach that merged biological knowledge-driven analysis and data-driven Bayesian probabilistic analysis.

Differential MicroRNA Analyses of Burkholderia pseudomallei- and Francisella tularensis-Exposed hPBMCs Reveal Potential Biomarkers.

Increasing evidence that microRNAs (miRNAs) play important roles in the immune response against infectious agents suggests that miRNA might be exploitable as signatures of exposure to specific infectious agents. In order to identify potential early miRNA biomarkers of bacterial infections, human peripheral blood mononuclear cells (hPBMCs) were exposed to two select agents, Burkholderia pseudomallei K96243 and Francisella tularensis SHU S4, as well as to the nonpathogenic control Escherichia coli DH5α. RNA samples were harvested at three early time points, 30, 60, and 120 minutes postexposure, then sequenced. RNAseq analyses identified 87 miRNAs to be differentially expressed (DE) in a linear fashion. Of these, 31 miRNAs were tested using the miScript miRNA qPCR assay. Through RNAseq identification and qPCR validation, we identified differentially expressed miRNA species that may be involved in the early response to bacterial infections. Based upon its upregulation at early time points postexposure in two different individuals, hsa-mir-30c-5p is a miRNA species that could be studied further as a potential biomarker for exposure to these gram-negative intracellular pathogens. Gene ontology functional analyses demonstrated that programmed cell death is the first ranking biological process associated with miRNAs that are upregulated in F. tularensis-exposed hPBMCs.

MicroRNAs Are Involved in the Development of Morphine-Induced Analgesic Tolerance and Regulate Functionally Relevant Changes in Serpini1.

Long-term opioid treatment results in reduced therapeutic efficacy and in turn leads to an increase in the dose required to produce equivalent pain relief and alleviate break-through or insurmountable pain. Altered gene expression is a likely means for inducing long-term neuroadaptations responsible for tolerance. Studies conducted by our laboratory (Tapocik et al., 2009) revealed a network of gene expression changes occurring in canonical pathways involved in neuroplasticity, and uncovered miRNA processing as a potential mechanism. In particular, the mRNA coding the protein responsible for processing miRNAs, Dicer1, was positively correlated with the development of analgesic tolerance. The purpose of the present study was to test the hypothesis that miRNAs play a significant role in the development of analgesic tolerance as measured by thermal nociception. Dicer1 knockdown, miRNA profiling, bioinformatics, and confirmation of high value targets were used to test the proposition. Regionally targeted Dicer1 knockdown (via shRNA) had the anticipated consequence of eliminating the development of tolerance in C57BL/6J (B6) mice, thus supporting the involvement of miRNAs in the development of tolerance. MiRNA expression profiling identified a core set of chronic morphine-regulated miRNAs (miR's 27a, 9, 483, 505, 146b, 202). Bioinformatics approaches were implemented to identify and prioritize their predicted target mRNAs. We focused our attention on miR27a and its predicted target serpin peptidase inhibitor clade I (Serpini1) mRNA, a transcript known to be intricately involved in dendritic spine density regulation in a manner consistent with chronic morphine's consequences and previously found to be correlated with the development of analgesic tolerance. In vitro reporter assay confirmed the targeting of the Serpini1 3'-untranslated region by miR27a. Interestingly miR27a was found to positively regulate Serpini1 mRNA and protein levels in multiple neuronal cell lines. Lastly, Serpini1 knockout mice developed analgesic tolerance at a slower rate than wild-type mice thus confirming a role for the protein in analgesic tolerance. Overall, these results provide evidence to support a specific role for miR27a and Serpini1 in the behavioral response to chronic opioid administration (COA) and suggest that miRNA expression and mRNA targeting may underlie the neuroadaptations that mediate tolerance to the analgesic effects of morphine.

Cross-talk of expression quantitative trait loci within 2 interacting blood pressure quantitative trait loci.

Genetic dissection of the S rat genome has provided strong evidence for the presence of 2 interacting blood pressure quantitative trait loci (QTLs), termed QTL1 and QTL2, on rat chromosome 5. However, the identities of the underlying interacting genetic factors remain unknown. Further experiments targeted to identify the interacting genetic factors by the substitution mapping approach alone are difficult because of the interdependency of natural recombinations to occur at the 2 QTLs. We hypothesized that the interacting genetic factors underlying these 2 QTLs may interact at the level of gene transcription and thereby represent expression QTLs or eQTLs. To detect these interacting expression QTLs, a custom QTL chip containing the annotated genes within QTL1 and QTL2 was developed and used to conduct a transcriptional profiling study of S and 2 congenic strains that retain either 1 or both of the QTLs. The results uncovered an interaction between 2 transcription factor genes, Dmrta2 and Nfia. Furthermore, the "biological signature" elicited by these 2 transcription factors was differential between the congenic strain that retained Lewis alleles at both QTL1 and QTL2 compared with the congenic strain that retained Lewis alleles at QTL1 alone. A network of transcription factors potentially affecting blood pressure could be traced, lending support to our hypothesis.

Assessing unmodified 70-mer oligonucleotide probe performance on glass-slide microarrays.

BACKGROUND: Long oligonucleotide microarrays are potentially more cost- and management-efficient than cDNA microarrays, but there is little information on the relative performance of these two probe types. The feasibility of using unmodified oligonucleotides to accurately measure changes in gene expression is also unclear. RESULTS: Unmodified sense and antisense 70-mer oligonucleotides representing 75 known rat genes and 10 Arabidopsis control genes were synthesized, printed and UV cross-linked onto glass slides. Printed alongside were PCR-amplified cDNA clones corresponding to the same genes, enabling us to compare the two probe types simultaneously. Our study was designed to evaluate the mRNA profiles of heart and brain, along with Arabidopsis cRNA spiked into the labeling reaction at different relative copy number. Hybridization signal intensity did not correlate with probe type but depended on the extent of UV irradiation. To determine the effect of oligonucleotide concentration on hybridization signal, 70-mers were serially diluted. No significant change in gene-expression ratio or loss in hybridization signal was detected, even at the lowest concentration tested (6.25 microm). In many instances, signal intensity actually increased with decreasing concentration. The correlation coefficient between oligonucleotide and cDNA probes for identifying differentially expressed genes was 0.80, with an average coefficient of variation of 13.4%. Approximately 8% of the genes showed discordant results with the two probe types, and in each case the cDNA results were more accurate, as determined by real-time PCR. CONCLUSIONS: Microarrays of UV cross-linked unmodified oligonucleotides provided sensitive and specific measurements for most of the genes studied.