Oviduct-Embryo Interactions in Cattle: Two-Way Traffic or a One-Way Street?

This study examined the effect of the presence of single or multiple embryos on the transcriptome of the bovine oviduct. In experiment 1, cyclic (nonbred, n = 6) and pregnant (artificially inseminated, n = 11) heifers were slaughtered on Day 3 after estrus, and the ampulla and isthmic regions of the oviduct ipsilateral to the corpus luteum were separately flushed. Oviductal epithelial cells from the isthmus region, in which all oocytes/embryos were located, were snap-frozen for microarray analysis. In experiment 2, heifers were divided into cyclic (nonbred, n = 6) or pregnant (multiple embryo transfer, n = 10) groups. In vitro-produced presumptive zygotes were transferred endoscopically to the ipsilateral oviduct on Day 1.5 postestrus (n = 50 zygotes/heifer). Heifers were slaughtered on Day 3, and oviductal isthmus epithelial cells were recovered for RNA sequencing. Microarray analysis in experiment 1 failed to detect any difference in the transcriptome of the oviductal isthmus induced by the presence of a single embryo. In experiment 2, following multiple embryo transfer, RNA sequencing revealed 278 differentially expressed genes, of which 123 were up-regulated and 155 were down-regulated in pregnant heifers. Most of the down-regulated genes were related to immune function. In conclusion, the presence of multiple embryos in the oviduct resulted in the detection of differentially expressed genes in the oviductal isthmus; failure to detect changes in the oviduct transcriptome in the presence of a single embryo may be due to the effect being local and undetectable under the conditions of this study.

Next-generation sequencing identifies TGF-ß1-associated gene expression profiles in renal epithelial cells reiterated in human diabetic nephropathy.

Transforming growth factor-beta (TGF-ß1) is implicated in the onset and progression of renal fibrosis and diabetic nephropathy (DN), leading to a loss of epithelial characteristics of tubular cells. The transcriptional profile of renal tubular epithelial cells stimulated with TGF-ß1 was assessed using RNA-Seq, with 2027 differentially expressed genes identified. Promoter analysis of transcription factor binding sites in the TGF-ß1 responsive gene set predicted activation of multiple transcriptional networks, including NFκB. Comparison of RNA-Seq with microarray data from identical experimental conditions identified low abundance transcripts exclusive to RNA-Seq data. We compared these findings to human disease by analyzing transcriptomic data from renal biopsies of patients with DN versus control groups, identifying a shared subset of 179 regulated genes. ARK5, encoding an AMP-related kinase, and TGFBI - encoding transforming growth factor, beta-induced protein were induced by TGF-ß1 and also upregulated in human DN. Suppression of ARK5 attenuated fibrotic responses of renal epithelia to TGF-ß1 exposure; and silencing of TGFBI induced expression of the epithelial cell marker - E-cadherin. We identified low abundance transcripts in sequence data and validated expression levels of several transcripts (ANKRD56, ENTPD8) in tubular enriched kidney biopsies of DN patients versus living donors. In conclusion, we have defined a TGF-ß1-driven pro-fibrotic signal in renal epithelial cells that is also evident in the DN renal transcriptome.

mab21l2 transgenics reveal novel expression patterns of mab21l1 and mab21l2, and conserved promoter regulation without sequence conservation.

mab21l1 and mab21l2 paralogs have widespread and dynamic expression patterns during vertebrate development. Both genes are expressed in the developing eye, midbrain, neural tube, and branchial arches. Our goal was to identify promoter regions with activity in mab21l2 expression domains. Assays of mab21l2 promoter-EGFP constructs in zebrafish embryos confirm that constructs containing 7.2 or 4.9 kb of mab21l2 promoter region are sufficient to drive expression in known (e.g., tectum, branchial arches) and unexpected domains (e.g., lens and retinal amacrine cells). A comparative analysis identifies complementary and novel expression domains of endogenous mab21l2 (e.g., lens and ventral iridocorneal canal) and mab21l1 (e.g., retinal amacrine and ganglion cells). Interestingly, therefore, despite the absence of conserved non-coding elements, a 4.9-kb mab21l2 promoter is sufficient to recapitulate expression in tissues unique to mab21l1 or mab21l2.

Chemostat-based transcriptional analysis of growth rate change and BCL-2 over-expression in NS0 cells.

We investigated the transcriptional response of NS0 cells undergoing controlled nutrient growth change in continuous chemostat culture using mouse microarrays. A 50% reduction in growth rate resulted in detectable alterations in the expression of 29 genes in NS0 cells. Notably, expression of genes in three major biological processes, namely transcriptional, translational, and protein processing functions, were modified. To further elucidate the advantage of the chemostat environment for establishment of "omic" data sets, an expression profile of the over-expressed gene bcl-2 in NS0 cells was probed. Functional analysis revealed that the underlying altered molecular mechanism was particularly associated with G1 cell cycle progression, protein synthesis, and apoptosis. Importantly, these findings agreed with the physical function of the cells. Despite an increase in survival rate, bcl-2 over-expression resulted in a decrease of specific productivity, glucose consumption, oxygen uptake rate and intracellular protein content, indicating a lower energy generating metabolism. Further, a prolongation of G1 cell cycle phase was evident on lowering the growth rate. Overall, the application of microarray analysis to chemostat-grown cultures offers an excellent combination for the interpretation of transcriptomic profiles to elucidate the molecular mechanisms during nutrient growth change and bcl-2 over-expression.

Bi-directional gene set enrichment and canonical correlation analysis identify key diet-sensitive pathways and biomarkers of metabolic syndrome.

BACKGROUND: Currently, a number of bioinformatics methods are available to generate appropriate lists of genes from a microarray experiment. While these lists represent an accurate primary analysis of the data, fewer options exist to contextualise those lists. The development and validation of such methods is crucial to the wider application of microarray technology in the clinical setting. Two key challenges in clinical bioinformatics involve appropriate statistical modelling of dynamic transcriptomic changes, and extraction of clinically relevant meaning from very large datasets. RESULTS: Here, we apply an approach to gene set enrichment analysis that allows for detection of bi-directional enrichment within a gene set. Furthermore, we apply canonical correlation analysis and Fisher's exact test, using plasma marker data with known clinical relevance to aid identification of the most important gene and pathway changes in our transcriptomic dataset. After a 28-day dietary intervention with high-CLA beef, a range of plasma markers indicated a marked improvement in the metabolic health of genetically obese mice. Tissue transcriptomic profiles indicated that the effects were most dramatic in liver (1270 genes significantly changed; p < 0.05), followed by muscle (601 genes) and adipose (16 genes). Results from modified GSEA showed that the high-CLA beef diet affected diverse biological processes across the three tissues, and that the majority of pathway changes reached significance only with the bi-directional test. Combining the liver tissue microarray results with plasma marker data revealed 110 CLA-sensitive genes showing strong canonical correlation with one or more plasma markers of metabolic health, and 9 significantly overrepresented pathways among this set; each of these pathways was also significantly changed by the high-CLA diet. Closer inspection of two of these pathways--selenoamino acid metabolism and steroid biosynthesis--illustrated clear diet-sensitive changes in constituent genes, as well as strong correlations between gene expression and plasma markers of metabolic syndrome independent of the dietary effect. CONCLUSION: Bi-directional gene set enrichment analysis more accurately reflects dynamic regulatory behaviour in biochemical pathways, and as such highlighted biologically relevant changes that were not detected using a traditional approach. In such cases where transcriptomic response to treatment is exceptionally large, canonical correlation analysis in conjunction with Fisher's exact test highlights the subset of pathways showing strongest correlation with the clinical markers of interest. In this case, we have identified selenoamino acid metabolism and steroid biosynthesis as key pathways mediating the observed relationship between metabolic health and high-CLA beef. These results indicate that this type of analysis has the potential to generate novel transcriptome-based biomarkers of disease.

Inter-laboratory comparison of human renal proximal tubule (HK-2) transcriptome alterations due to Cyclosporine A exposure and medium exhaustion.

There is an acknowledged need to promote and further develop in vitro techniques in order to achieve the goal of improved risk assessment of chemicals and pharmaceuticals to humans. The EU 6th framework project "PREDICTOMICS" was established in order to contribute to the further development of in vitro toxicology, with a particular focus on emerging techniques including toxicogenomics. DNA microarray technology is being used more frequently in the in vitro field, however, only very few studies have assessed the reproducibility of this technique with respect to in vitro toxicology. To this end we conducted an interlaboratory comparison to test the reproducibility of transcriptomic changes induced by the immunosuppressive agent, Cyclosporine A (CsA) on the human renal proximal tubular cell line, HK-2 cell. Four European laboratories took part in this study. Under standardised conditions, each laboratory treated HK-2 cells with 5microM CsA for 12 and 48h. RNA was isolated and hybridised to Affymetrix HGU-133 plus two arrays at three different sites. Analysis of the transcription profiles demonstrated that one laboratory clustered away from the other laboratories, potentially due to an inclusion of a trypsinisation step by this laboratory. Once the genes responsible for this separate clustering were removed all laboratories showed similar expression profiles. There was a major impact of time since feed, due to medium exhaustion in the 48h arrays compared to the 12h arrays, regardless of CsA treatment. Biological processes including general vesicle transport, amino acid metabolism, amino acid transport and amino acid biosynthesis were over-represented due to time since feed, while cell cycle, DNA replication, mitosis and DNA metabolism were under-represented. CsA responsive genes were involved in cell cycle, the p53 pathway and Wnt signaling. Additionally there was an overlap of differentially expressed genes due to CsA and medium exhaustion which is most likely due to CsA induced glycolysis. The glucose deprivation dependent genes HspA5 and GP96 and the Hsp70 chaperones DNAJ/Hsp40, DNAJ/HspB9, DNAJ/HspC3 DNAJ/HspC10 were induced by both CsA and medium exhaustion. We conclude that under standardised conditions the application of Affymetrix DNA microarrays to in vitro toxiciological studies are satisfactorily reproducible. However, confounding factors such as medium exhaustion must also be considered in such analyses.

Complete Genome Sequence of Clostridium estertheticum DSM 8809, a Microbe Identified in Spoiled Vacuum Packed Beef.

Blown pack spoilage (BPS) is a major issue for the beef industry. Etiological agents of BPS involve members of a group of Clostridium species, including Clostridium estertheticum which has the ability to produce gas, mostly carbon dioxide, under anaerobic psychotrophic growth conditions. This spore-forming bacterium grows slowly under laboratory conditions, and it can take up to 3 months to produce a workable culture. These characteristics have limited the study of this commercially challenging bacterium. Consequently information on this bacterium is limited and no effective controls are currently available to confidently detect and manage this production risk. In this study the complete genome of C. estertheticum DSM 8809 was determined by SMRT® sequencing. The genome consists of a circular chromosome of 4.7 Mbp along with a single plasmid carrying a potential tellurite resistance gene tehB and a Tn3-like resolvase-encoding gene tnpR. The genome sequence was searched for central metabolic pathways that would support its biochemical profile and several enzymes contributing to this phenotype were identified. Several putative antibiotic/biocide/metal resistance-encoding genes and virulence factors were also identified in the genome, a feature that requires further research. The availability of the genome sequence will provide a basic blueprint from which to develop valuable biomarkers that could support and improve the detection and control of this bacterium along the beef production chain.

Transcriptomic Profiling of Sequential Tumors from Breast Cancer Patients Provides a Global View of Metastatic Expression Changes Following Endocrine Therapy.

PURPOSE: Disease recurrence is a common problem in breast cancer and yet the mechanisms enabling tumor cells to evade therapy and colonize distant organs remain unclear. We sought to characterize global expression changes occurring with metastatic disease progression in the endocrine-resistant setting. EXPERIMENTAL DESIGN: Here, for the first time, RNAsequencing has been performed on matched primary, nodal, and liver metastatic tumors from tamoxifen-treated patients following disease progression. Expression of genes commonly elevated in the metastases of sequenced patients was subsequently examined in an extended matched patient cohort with metastatic disease from multiple sites. The impact of tamoxifen treatment on endocrine-resistant tumors in vivo was investigated in a xenograft model. RESULTS: The extent of patient heterogeneity at the gene level was striking. Less than 3% of the genes differentially expressed between sequential tumors were common to all patients. Larger divergence was observed between primary and liver tumors than between primary and nodal tumors, reflecting both the latency to disease progression and the genetic impact of intervening therapy. Furthermore, an endocrine-resistant in vivo mouse model demonstrated that tamoxifen treatment has the potential to drive disease progression and establish distant metastatic disease. Common functional pathways altered during metastatic, endocrine-resistant progression included extracellular matrix receptor interactions and focal adhesions. CONCLUSIONS: This novel global analysis highlights the influence of primary tumor biology in determining the transcriptomic profile of metastatic tumors, as well as the need for adaptations in cell-cell communications to facilitate successful tumor cell colonization of distant host organs.

Global characterization of the SRC-1 transcriptome identifies ADAM22 as an ER-independent mediator of endocrine-resistant breast cancer.

The development of breast cancer resistance to endocrine therapy results from an increase in cellular plasticity that permits the emergence of a hormone-independent tumor. The steroid coactivator protein SRC-1, through interactions with developmental proteins and other nonsteroidal transcription factors, drives this tumor adaptability. In this discovery study, we identified ADAM22, a non-protease member of the ADAM family of disintegrins, as a direct estrogen receptor (ER)-independent target of SRC-1. We confirmed SRC-1 as a regulator of ADAM22 by molecular, cellular, and in vivo studies. ADAM22 functioned in cellular migration and differentiation, and its levels were increased in endocrine resistant-tumors compared with endocrine-sensitive tumors in mouse xenograft models of human breast cancer. Clinically, ADAM22 was found to serve as an independent predictor of poor disease-free survival. Taken together, our findings suggest that SRC-1 switches steroid-responsive tumors to a steroid-resistant state in which the SRC-1 target gene ADAM22 has a critical role, suggesting this molecule as a prognostic and therapeutic drug target that could help improve the treatment of endocrine-resistant breast cancer.

AIB1:ERα transcriptional activity is selectively enhanced in aromatase inhibitor-resistant breast cancer cells.

PURPOSE: The use of aromatase inhibitors (AI) in the treatment of estrogen receptor (ER)-positive, postmenopausal breast cancer has proven efficacy. However, inappropriate activation of ER target genes has been implicated in the development of resistant tumors. The ER coactivator protein AIB1 has previously been associated with initiation of breast cancer and resistance to endocrine therapy. EXPERIMENTAL DESIGN: Here, we investigated the role of AIB1 in the deregulation of ER target genes occurring as a consequence of AI resistance using tissue microarrays of patients with breast cancer and cell line models of resistance to the AI letrozole. RESULTS: Expression of AIB1 associated with disease recurrence (P = 0.025) and reduced disease-free survival time (P = 0.0471) in patients treated with an AI as first-line therapy. In a cell line model of resistance to letrozole (LetR), we found ERα/AIB1 promoter recruitment and subsequent expression of the classic ER target genes pS2 and Myc to be constitutively upregulated in the presence of both androstenedione and letrozole. In contrast, the recruitment of the ERα/AIB1 transcriptional complex to the nonclassic ER target cyclin D1 and its subsequent expression remained sensitive to steroid treatment and could be inhibited by treatment with letrozole. Molecular studies revealed that this may be due in part to direct steroid regulation of c-jun-NH(2)-kinase (JNK), signaling to Jun and Fos at the cyclin D1 promoter. CONCLUSION: This study establishes a role for AIB1 in AI-resistant breast cancer and describes a new mechanism of ERα/AIB1 gene regulation which could contribute to the development of an aggressive tumor phenotype.

Early proteomic analysis may allow noninvasive identification of hepatitis C response to treatment with pegylated interferon α-2b and ribavirin.

BACKGROUND AND AIM: Chronic hepatitis C virus (HCV) infection represents a significant disease burden worldwide. Approximately 170 million people are chronically infected. HCV can lead to liver fibrosis, cirrhosis and hepatocellular carcinoma. Current standard treatment with pegylated interferon and ribavirin is suboptimal and up to 60% of patients fail to respond. Week 4 and 12 HCV RNA is used as a marker of response with nonresponders at 12 weeks discontinuing treatment. Earlier identification of nonresponders using novel biomarkers would be beneficial in preventing unnecessary toxicities and cost. This study profiled the proteomic response to treatment in HCV patients within the first 24 h using surface-enhanced laser desorption-ionization time-of-flight mass spectrometry (SELDI-TOF MS). METHODS: Serum from 25 HCV infected individuals during the initial 24 h of treatment was profiled using SELDI-TOF MS. Arrays were analyzed on the ProteinChip Reader and time-of-flight spectra were generated. Peak detection was performed by Biomarker Wizard software and analyzed using BioConductor packages. RESULTS: Significant differences were seen between the proteomic profiles of responders and nonresponders to treatment. Overall 70 peaks differentiated responders from nonresponders. A random forest classifier identified a panel of 20 peaks, which differentiated responders from nonresponders with 87.4% accuracy. The CM10 chip revealed 16 peaks identifying genotype 1 responders from nonresponders. CONCLUSION: This study identifies early proteomic spectra as potential predictors of HCV treatment response using SELDI-TOF MS. This illustrates the importance of early biomarkers in the prediction of response within the first 24 h, which may aid in tailoring HCV treatment regimens.

High-throughput proteomics detection of novel splice isoforms in human platelets.

Alternative splicing (AS) is an intrinsic regulatory mechanism of all metazoans. Recent findings suggest that 100% of multiexonic human genes give rise to splice isoforms. AS can be specific to tissue type, environment or developmentally regulated. Splice variants have also been implicated in various diseases including cancer. Detection of these variants will enhance our understanding of the complexity of the human genome and provide disease-specific and prognostic biomarkers. We adopted a proteomics approach to identify exon skip events - the most common form of AS. We constructed a database harboring the peptide sequences derived from all hypothetical exon skip junctions in the human genome. Searching tandem mass spectrometry (MS/MS) data against the database allows the detection of exon skip events, directly at the protein level. Here we describe the application of this approach to human platelets, including the mRNA-based verification of novel splice isoforms of ITGA2, NPEPPS and FH. This methodology is applicable to all new or existing MS/MS datasets.

Analysis of the hypervariable region of the Salmonella enterica genome associated with tRNA(leuX).

The divergence of Salmonella enterica and Escherichia coli is estimated to have occurred approximately 140 million years ago. Despite this evolutionary distance, the genomes of these two species still share extensive synteny and homology. However, there are significant differences between the two species in terms of genes putatively acquired via various horizontal transfer events. Here we report on the composition and distribution across the Salmonella genus of a chromosomal region designated SPI-10 in Salmonella enterica serovar Typhi and located adjacent to tRNA(leuX). We find that across the Salmonella genus the tRNA(leuX) region is a hypervariable hot spot for horizontal gene transfer; different isolates from the same S. enterica serovar can exhibit significant variation in this region. Many P4 phage, plasmid, and transposable element-associated genes are found adjacent to tRNA(leuX) in both Salmonella and E. coli, suggesting that these mobile genetic elements have played a major role in driving the variability of this region.

Composition, acquisition, and distribution of the Vi exopolysaccharide-encoding Salmonella enterica pathogenicity island SPI-7.

Vi capsular polysaccharide production is encoded by the viaB locus, which has a limited distribution in Salmonella enterica serovars. In S. enterica serovar Typhi, viaB is encoded on a 134-kb pathogenicity island known as SPI-7 that is located between partially duplicated tRNA(pheU) sites. Functional and bioinformatic analysis suggests that SPI-7 has a mosaic structure and may have evolved as a consequence of several independent insertion events. Analysis of viaB-associated DNA in Vi-positive S. enterica serovar Paratyphi C and S. enterica serovar Dublin isolates revealed the presence of similar SPI-7 islands. In S. enterica serovars Paratyphi C and Dublin, the SopE bacteriophage and a 15-kb fragment adjacent to the intact tRNA(pheU) site were absent. In S. enterica serovar Paratyphi C only, a region encoding a type IV pilus involved in the adherence of S. enterica serovar Typhi to host cells was missing. The remainder of the SPI-7 islands investigated exhibited over 99% DNA sequence identity in the three serovars. Of 30 other Salmonella serovars examined, 24 contained no insertions at the equivalent tRNA(pheU) site, 2 had a 3.7-kb insertion, and 4 showed sequence variation at the tRNA(pheU)-phoN junction, which was not analyzed further. Sequence analysis of the SPI-7 region from S. enterica serovar Typhi strain CT18 revealed significant synteny with clusters of genes from a variety of saprophytic bacteria and phytobacteria, including Pseudomonas aeruginosa and Xanthomonas axonopodis pv. citri. This analysis suggested that SPI-7 may be a mobile element, such as a conjugative transposon or an integrated plasmid remnant.