Gene expression profiling of the left ventricles in a rat model of intrinsic aerobic running capacity

Sex-related differences in gene expression by porcine aortic valvular interstitial cells

While many large-scale risk factors for calcific aortic valve disease (CAVD) have been identified, the molecular etiology and subsequent pathogenesis of CAVD have yet to be fully understood. Specifically, it is unclear what biological phenomena underlie the significantly higher occurrence of CAVD in the male population. We hypothesized the existence of intrinsic, cellular-scale differences between male and female valvular interstitial cells (VICs) that contribute to male sex being a risk factor for CAVD. Differences in gene expression profiles between healthy male and female porcine VICs were investigated via microarray analysis. Mean expression values of each probe set in the male samples were compared to the female samples, and biological processes were analyzed for overrepresentation using Gene Ontology term enrichment analysis. There were 183 genes identified as significantly (fold change>2; P<0.05) different in male versus female aortic valve leaflets. Within this significant gene list there were 298 overrepresented biological processes, several of which are relevant to pathways identified in CAVD pathogenesis. In particular, pathway analysis indicated that cellular proliferation, apoptosis, migration, ossification, angiogenesis, inflammation, and extracellular matrix reorganization were all significantly represented in the data set. These gene expression findings also translated into functional differences in VIC behavior in the in vitro environment, as sex-related differences in proliferation and apoptosis were confirmed in VIC populations cultured in vitro. These data suggest that a sex-related propensity for CAVD exists on the cellular level in healthy subjects, a phenomenon that could have significant clinical implications. These findings also strongly support discontinuing the use of mixed-sex VIC cultures, thereby changing the current standard in the field.

Systematic approaches towards the development of host-directed antiviral therapeutics

Since the onset of antiviral therapy, viral resistance has compromised the clinical value of small-molecule drugs targeting pathogen components. As intracellular parasites, viruses complete their life cycle by hijacking a multitude of host-factors. Aiming at the latter rather than the pathogen directly, host-directed antiviral therapy has emerged as a concept to counteract evolution of viral resistance and develop broad-spectrum drug classes. This approach is propelled by bioinformatics analysis of genome-wide screens that greatly enhance insights into the complex network of host-pathogen interactions and generate a shortlist of potential gene targets from a multitude of candidates, thus setting the stage for a new era of rational identification of drug targets for host-directed antiviral therapies. With particular emphasis on human immunodeficiency virus and influenza virus, two major human pathogens, we review screens employed to elucidate host-pathogen interactions and discuss the state of database ontology approaches applicable to defining a therapeutic endpoint. The value of this strategy for drug discovery is evaluated, and perspectives for bioinformatics-driven hit identification are outlined.

Serious limitations of the QTL/microarray approach for QTL gene discovery

BACKGROUND

It has been proposed that the use of gene expression microarrays in nonrecombinant parental or congenic strains can accelerate the process of isolating individual genes underlying quantitative trait loci (QTL). However, the effectiveness of this approach has not been assessed.

RESULTS

Thirty-seven studies that have implemented the QTL/microarray approach in rodents were reviewed. About 30% of studies showed enrichment for QTL candidates, mostly in comparisons between congenic and background strains. Three studies led to the identification of an underlying QTL gene. To complement the literature results, a microarray experiment was performed using three mouse congenic strains isolating the effects of at least 25 biometric QTL. Results show that genes in the congenic donor regions were preferentially selected. However, within donor regions, the distribution of differentially expressed genes was homogeneous once gene density was accounted for. Genes within identical-by-descent (IBD) regions were less likely to be differentially expressed in chromosome 2, but not in chromosomes 11 and 17. Furthermore, expression of QTL regulated in cis (cis eQTL) showed higher expression in the background genotype, which was partially explained by the presence of single nucleotide polymorphisms (SNP).

CONCLUSIONS

The literature shows limited successes from the QTL/microarray approach to identify QTL genes. Our own results from microarray profiling of three congenic strains revealed a strong tendency to select cis-eQTL over trans-eQTL. IBD regions had little effect on rate of differential expression, and we provide several reasons why IBD should not be used to discard eQTL candidates. In addition, mismatch probes produced false cis-eQTL that could not be completely removed with the current strains genotypes and low probe density microarrays. The reviewed studies did not account for lack of coverage from the platforms used and therefore removed genes that were not tested. Together, our results explain the tendency to report QTL candidates as differentially expressed and indicate that the utility of the QTL/microarray as currently implemented is limited. Alternatives are proposed that make use of microarray data from multiple experiments to overcome the outlined limitations.

Heart transcriptome of the bank vole (Myodes glareolus): towards understanding the evolutionary variation in metabolic rate

BACKGROUND

Understanding the genetic basis of adaptive changes has been a major goal of evolutionary biology. In complex organisms without sequenced genomes, de novo transcriptome assembly using a longer read sequencing technology followed by expression profiling using short reads is likely to provide comprehensive identification of adaptive variation at the expression level and sequence polymorphisms in coding regions. We performed sequencing and de novo assembly of the bank vole heart transcriptome in lines selected for high metabolism and unselected controls.

RESULTS

A single 454 Titanium run produced over million reads, which were assembled into 63,581 contigs. Searches against the SwissProt protein database and the ENSEMBL collection of mouse transcripts detected similarity to 11,181 and 14,051 genes, respectively. As judged by the representation of genes from the heart-related Gene Ontology categories and UniGenes detected in the mouse heart, our detection of the genes expressed in the heart was nearly complete (> 95% and almost 90% respectively). On average, 38.7% of the transcript length was covered by our sequences, with notably higher (45.0%) coverage of coding regions than of untranslated regions (24.5% of 5' and 32.7% of 3'UTRs). Lower sequence conservation between mouse and bank vole in untranslated regions was found to be partially responsible for poorer UTR representation. Our data might suggest a widespread transcription from noncoding genomic regions, a finding not reported in previous studies regarding transcriptomes in non-model organisms. We also identified over 19 thousand putative single nucleotide polymorphisms (SNPs). A much higher fraction of the SNPs than expected by chance exhibited variant frequency differences between selection regimes.

CONCLUSION

Longer reads and higher sequence yield per run provided by the 454 Titanium technology in comparison to earlier generations of pyrosequencing proved beneficial for the quality of assembly. An almost full representation of genes known to be expressed in the mouse heart was identified. Usage of the extensive genomic resources available for the house mouse, a moderately (20-40 mln years) divergent relative of the voles, enabled a comprehensive assessment of the transcript completeness. Transcript sequences generated in the present study allowed the identification of candidate SNPs associated with divergence of selection lines and constitute a valuable permanent resource forming a foundation for RNAseq experiments aiming at detection of adaptive changes both at the level of gene expression and sequence variants, that would facilitate studies of the genetic basis of evolutionary divergence.

Selecting exercise regimens and strains to modify obesity and diabetes in rodents: an overview

Exercise is part of a healthy lifestyle and frequently is an important component in combating chronic diseases, such as obesity and diabetes. Understanding the molecular events initiated by regular exercise is best studied in laboratory animals, with mice and rats being favoured for a number of reasons. However, the wide variety of rodent strains available for biomedical research often makes it challenging to select an animal strain suitable for studying specific disease outcomes. In the present review we focus on exercise as a management strategy for obesity and diabetes and we discuss: (i) exercise paradigms in humans shown to ameliorate signs and symptoms of obesity and diabetes; (ii) different rodent strains in terms of their advantages, disadvantages and limitations when using specific forms of exercise; (iii) the strengths and weaknesses of commonly used laboratory methods for rodent exercise; and (iv) the unintended consequences of exercise that are often manifested by increased hormonal and oxidative stress responses.

Retroviral gene insertion in breast milk mediated lymphomagenesis

We have demonstrated breast milk transmitted MoMuLV-ts1 retrovirus infection and subsequent lymphoma development in offspring of uninfected mothers suckled by infected surrogate mothers. Additionally, we have shown that the lymphoma development occurs as a result of viral gene integration into host genome. A total of 146 pups from Balb/C mice were divided into 5 groups; one control and 4 experimental. All offspring suckled from surrogate infected or control mothers, except one group of infected pups left with their biological mothers. Thirteen of 91 infected pups developed lymphoma. Inverse-PCR, DNA cloning, and quantitative real-time PCR (qRT-PCR) were used to study the virus integration sites (VIS) and alterations in gene expression. VIS were randomly distributed throughout the genome. The majority of insertion sites were found in chromosomes 10, 12 and 13. A total of 209 proviral genomic insertion sites were located with 52 intragenic and 157 intergenic sites. We have identified 29 target genes. Four genes including Tacc3, Aurka, Gfi1 and Ahi1 showed the maximum upregulation of mRNA expression. These four genes can be considered as candidate genes based on their association with cancer. Upregulation of these genes may be involved in this type of lymphoma development. This model provides an important opportunity to gain insight into the relationship of viral gene insertion into host genome and development of lymphoma via natural transmission route such as breast milk.

Molecular networks discriminating mouse bladder responses to intravesical bacillus Calmette-Guerin (BCG), LPS, and TNF-alpha

Background

Despite being a mainstay for treating superficial bladder carcinoma and a promising agent for interstitial cystitis, the precise mechanism of Bacillus Calmette-Guerin (BCG) remains poorly understood. It is particularly unclear whether BCG is capable of altering gene expression in the bladder target organ beyond its well-recognized pro-inflammatory effects and how this relates to its therapeutic efficacy. The objective of this study was to determine differentially expressed genes in the mouse bladder following chronic intravesical BCG therapy and to compare the results to non-specific pro inflammatory stimuli (LPS and TNF-α). For this purpose, C57BL/6 female mice received four weekly instillations of BCG, LPS, or TNF-α. Seven days after the last instillation, the urothelium along with the submucosa was removed from detrusor muscle and the RNA was extracted from both layers for cDNA array experiments. Microarray results were normalized by a robust regression analysis and only genes with an expression above a conditional threshold of 0.001 (3SD above background) were selected for analysis. Next, genes presenting a 3-fold ratio in regard to the control group were entered in Ingenuity Pathway Analysis (IPA) for a comparative analysis in order to determine genes specifically regulated by BCG, TNF-α, and LPS. In addition, the transcriptome was precipitated with an antibody against RNA polymerase II and real-time polymerase chain reaction assay (Q-PCR) was used to confirm some of the BCG-specific transcripts.

Results

Molecular networks of treatment-specific genes generated several hypotheses regarding the mode of action of BCG. BCG-specific genes involved small GTPases and BCG-specific networks overlapped with the following canonical signaling pathways: axonal guidance, B cell receptor, aryl hydrocarbon receptor, IL-6, PPAR, Wnt/β-catenin, and cAMP. In addition, a specific detrusor network expressed a high degree of overlap with the development of the lymphatic system. Interestingly, TNF-α-specific networks overlapped with the following canonical signaling pathways: PPAR, death receptor, and apoptosis. Finally, LPS-specific networks overlapped with the LPS/IL-1 mediated inhibition of RXR. Because NF-kappaB occupied a central position in several networks, we further determined whether this transcription factor was part of the responses to BCG. Electrophoretic mobility shift assays confirmed the participation of NF-kappaB in the mouse bladder responses to BCG. In addition, BCG treatment of a human urothelial cancer cell line (J82) also increased the binding activity of NF-kappaB, as determined by precipitation of the chromatin by a NF-kappaB-p65 antibody and Q-PCR of genes bearing a NF-kappaB consensus sequence. Next, we tested the hypothesis of whether small GTPases such as LRG-47 are involved in the uptake of BCG by the bladder urothelium.

Conclusion

As expected, BCG treatment induces the transcription of genes belonging to common pro-inflammatory networks. However, BCG also induces unique genes belonging to molecular networks involved in axonal guidance and lymphatic system development within the bladder target organ. In addition, NF-kappaB seems to play a predominant role in the bladder responses to BCG therapy. Finally, in intact urothelium, BCG-GFP internalizes in LRG-47-positive vesicles.

These results provide a molecular framework for the further study of the involvement of immune and nervous systems in the bladder responses to BCG therapy.

Repeated BCG treatment of mouse bladder selectively stimulates small GTPases and HLA antigens and inhibits single-spanning uroplakins

BACKGROUND

Despite being a mainstay for treating superficial bladder carcinoma and a promising agent for interstitial cystitis, the precise mechanism of Bacillus Calmette-Guerin (BCG) remains poorly understood. It is particularly unclear whether BCG is capable of altering gene expression beyond its well-recognized pro-inflammatory effects and how this relates to its therapeutic efficacy. The objective of this study was to determine differentially expressed genes in the mouse bladder following repeated intravesical BCG therapy.

METHODS

Mice were transurethrally instilled with BCG or pyrogen-free on days 1, 7, 14, and 21. Seven days after the last instillation, urothelia along with the submucosa was removed and amplified ds-DNA was prepared from control- and BCG-treated bladder mucosa and used to generate suppression subtractive hybridization (SSH). Plasmids from control- and BCG-specific differentially expressed clones and confirmed by Virtual Northern were then purified and the inserts were sequenced and annotated. Finally, chromatin immune precipitation combined with real-time polymerase chain reaction assay (ChIP/Q-PCR) was used to validate SSH-selected transcripts.

RESULTS

Repeated intravesical BCG treatment induced an up regulation of genes associated with antigen presentation (B2M, HLA-A, HLA-DQA1, HLA-DQB2, HLA-E, HLA-G, IGHG, and IGH) and representatives of two IFNgamma-induced small GTPase families: the GBPs (GBP1, GBP2, and GBP5) and the p47GTPases (IIGTP1, IIGTP2, and TGTP). Genes expressed in saline-treated bladders but down-regulated by BCG included: the single-spanning uroplakins (UPK3a and UPK2), SPRR2G, GSTM5, and RSP 19.

CONCLUSION

Here we introduced a hypothesis-generator approach to determine key genes involved in the urothelium/sumbmucosa responses to BCG therapy. Urinary bladder responds to repeated BCG treatment by up-regulating not only antigen presentation-related genes, but also GBP and p47 small GTPases, both potentially serving to mount a resistance to the replication of the Mycobacterium. It will be of tremendous future interest to determine whether these immune response cascades play a role in the anti-cancer effects exerted by BCG.

Regulatory network of inflammation downstream of proteinase-activated receptors

Background

Protease-activated receptors (PAR) are present in the urinary bladder, and their expression is altered in response to inflammation. PARs are a unique class of G protein-coupled that carry their own ligands, which remain cryptic until unmasked by proteolytic cleavage. Although the canonical signal transduction pathway downstream of PAR activation and coupling with various G proteins is known and leads to the rapid transcription of genes involved in inflammation, the effect of PAR activation on the downstream transcriptome is unknown.

We have shown that intravesical administration of PAR-activating peptides leads to an inflammatory reaction characterized by edema and granulocyte infiltration. Moreover, the inflammatory response to intravesical instillation of known pro-inflammatory stimuli such as E. coli lipopolysaccharide (LPS), substance P (SP), and antigen was strongly attenuated by PAR1- and to a lesser extent by PAR2-deficiency.

Results

Here, cDNA array experiments determined inflammatory genes whose expression is dependent on PAR1 activation. For this purpose, we compared the alteration in gene expression in wild type and PAR1^-/- mice induced by classical pro-inflammatory stimuli (LPS, SP, and antigen). 75 transcripts were considered to be dependent on PAR-1 activation and further annotated in silico by Ingenuity Pathways Analysis (IPA) and gene ontology (GO). Selected transcripts were target validated by quantitative PCR (Q-PCR). Among PAR1-dependent transcripts, the following have been implicated in the inflammatory process: b2m, ccl7, cd200, cd63, cdbpd, cfl1, dusp1, fkbp1a, fth1, hspb1, marcksl1, mmp2, myo5a, nfkbia, pax1, plaur, ppia, ptpn1, ptprcap, s100a10, sim2, and tnfaip2. However, a balanced response to signals of injury requires a transient cellular activation of a panel of genes together with inhibitory systems that temper the overwhelming inflammation. In this context, the activation of genes such as dusp1 and nfkbia seems to counter-balance the inflammatory response to PAR activation by limiting prolonged activation of p38 MAPK and increased cytokine production. In contrast, transcripts such as arf6 and dcnt1 that are involved in the mechanism of PAR re-sensitization would tend to perpetuate the inflammatory reaction in response to common pro-inflammatory stimuli.

Conclusion

The combination of cDNA array results and genomic networks reveals an overriding participation of PAR1 in bladder inflammation, provides a working model for the involvement of downstream signaling, and evokes testable hypotheses regarding the transcriptome downstream of PAR1 activation.

It remains to be determined whether or not mechanisms targeting PAR1 gene silencing or PAR1 blockade will ameliorate the clinical manifestation of cystitis.

Congenic strains confirm aerobic running capacity quantitative trait loci on rat chromosome 16 and identify possible intermediate phenotypes

We previously identified two inbred rat strains divergent for treadmill aerobic running capacity (ARC), the low-performing Copenhagen (COP) and the high-performing DA rats, and used an F(2)(COPxDA) population to identify ARC quantitative trait loci (QTLs) on rat chromosome 16 (RNO16) and the proximal portion of rat chromosome 3 (RNO3). Two congenic rat strains were bred to further investigate these ARC QTLs by introgressing RNO16 and the proximal portion of RNO3 from DA rats into the genetic background of COP rats and were named COP.DA(chr 16) and COP.DA(chr 3), respectively. COP.DA(chr 16) rats had significantly greater ARC compared with COP rats (696.7 +/- 38.2 m vs. 571.9 +/- 27.5 m, P = 0.03). COP.DA(chr 3) rats had increased, although not significant, ARC compared with COP rats (643.6 +/- 40.9 m vs. 571.9 +/- 27.5 m). COP.DA(chr 16) rats had significantly greater subcutaneous abdominal fat, as well as decreased fasting triglyceride levels, compared with COP rats (P < 0.05), indicating that genes responsible for strain differences in fat metabolism are also located on RNO16. While this colocalization of QTLs may be coincidental, it is also possible that these differences in energy balance may be associated with the superior running performance of COP.DA(chr 16) consomic rats.