Profiling of Luteal Transcriptome during Prostaglandin F2-Alpha Treatment in Buffalo Cows: Analysis of Signaling Pathways Associated with Luteolysis

In several species including the buffalo cow, prostaglandin (PG) F-2 alpha is the key molecule responsible for regression of corpus luteum (CL). Experiments were carried out to characterize gene expression changes in the CL tissue at various time points after administration of luteolytic dose of PGF(2 alpha) in buffalo cows. Circulating progesterone levels decreased within 1 h of PGF(2 alpha) treatment and evidence of apoptosis was demonstrable at 18 h post treatment. Microarray analysis indicated expression changes in several of immediate early genes and transcription factors within 3 h of treatment. Also, changes in expression of genes associated with cell to cell signaling, cytokine signaling, steroidogenesis, PG synthesis and apoptosis were observed. Analysis of various components of LH/CGR signaling in CL tissues indicated decreased LH/CGR protein expression, pCREB levels and PKA activity post PGF(2 alpha) treatment. The novel finding of this study is the down regulation of CYP19A1 gene expression accompanied by decrease in expression of E-2 receptors and circulating and intra luteal E-2 post PGF(2 alpha) treatment. Mining of microarray data revealed several differentially expressed E-2 responsive genes. Since CYP19A1 gene expression is low in the bovine CL, mining of microarray data of PGF(2 alpha)-treated macaques, the species with high luteal CYP19A1 expression, showed good correlation between differentially expressed E-2 responsive genes between both the species. Taken together, the results of this study suggest that PGF(2 alpha) interferes with luteotrophic signaling, impairs intraluteal E-2 levels and regulates various signaling pathways before the effects on structural luteolysis are manifest.

Quantitative RNA-seq analysis of the campylobacter jejuni transcriptome

Campylobacter jejuni is the most common bacterial cause of foodborne disease in the developed world. Its general physiology and biochemistry, as well as the mechanisms enabling it to colonize and cause disease in various hosts, are not well understood, and new approaches are required to understand its basic biology. High-throughput sequencing technologies provide unprecedented opportunities for functional genomic research. Recent studies have shown that direct Illumina sequencing of cDNA (RNA-seq) is a useful technique for the quantitative and qualitative examination of transcriptomes. In this study we report RNA-seq analyses of the transcriptomes of C. jejuni (NCTC11168) and its rpoN mutant. This has allowed the identification of hitherto unknown transcriptional units, and further defines the regulon that is dependent on rpoN for expression. The analysis of the NCTC11168 transcriptome was supplemented by additional proteomic analysis using liquid chromatography-MS. The transcriptomic and proteomic datasets represent an important resource for the Campylobacter research community. © 2011 SGM.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Profiling of the BRCA1 transcriptome through microarray and ChIP-chip analysis

A role for BRCA1 in the direct and indirect regulation of transcription is well established. However, a comprehensive view of the degree to which BRCA1 impacts transcriptional regulation on a genome-wide level has not been defined. We performed genome-wide expression profiling and ChIP-chip analysis, comparison of which revealed that although BRCA1 depletion results in transcriptional changes in 1294 genes, only 44 of these are promoter bound by BRCA1. However, 27 of these transcripts were linked to transcriptional regulation possibly explaining the large number of indirect transcriptional changes observed by microarray analysis. We show that no specific consensus sequence exists for BRCA1 DNA binding but rather demonstrate the presence of a number of known and novel transcription factor (TF)- binding sites commonly found on BRCA1 bound promoters. Co-immunoprecipitations confirmed that BRCA1 interacts with a number of these TFs including AP2-a, PAX2 and ZF5. Finally, we show that BRCA1 is bound to a subset of promoters of genes that are not altered by BRCA1 loss, but are transcriptionally regulated in a BRCA1-dependent manner upon DNA damage. These data suggest a model, whereby BRCA1 is present on defined promoters as part of an inactive complex poised to respond to various genotoxic stimuli. © The Author(s) 2011. Published by Oxford University Press.

Evolutionary analysis of novel serine proteases in the venom gland transcriptome of Bitis gabonica rhinoceros

Background: Serine proteases are major components of viper venom and target various stages of the blood coagulation system in victims and prey. A better understanding of the diversity of serine proteases and other enzymes present in snake venom will help to understand how the complexity of snake venom has evolved and will aid the development of novel therapeutics for treating snake bites. Methodology and Principal Findings: Four serine protease-encoding genes from the venom gland transcriptome of Bitis gabonica rhinoceros were amplified and sequenced. Mass spectrometry suggests the four enzymes corresponding to these genes are present in the venom of B. g. rhinoceros. Two of the enzymes, rhinocerases 2 and 3 have substitutions to two of the serine protease catalytic triad residues and are thus unlikely to be catalytically active, though they may have evolved other toxic functions. The other two enzymes, rhinocerases 4 and 5, have classical serine protease catalytic triad residues and thus are likely to be catalytically active, however they have glycine rather than the more typical aspartic acid at the base of the primary specificity pocket (position 189). Based on a detailed analysis of these sequences we suggest that alternative splicing together with individual amino acid mutations may have been involved in their evolution. Changes within amino acid segments which were previously proposed to undergo accelerated change in venom serine proteases have also been observed. Conclusions and Significance: Our study provides further insight into the diversity of serine protease isoforms present within snake venom and discusses their possible functions and how they may have evolved. These multiple serine protease isoforms with different substrate specificities may enhance the envenomation effects and help the snake to adapt to new habitats and diets. Our findings have potential for helping the future development of improved therapeutics for snake bites.

Analysis of Bothrops jararacussu venomous gland transcriptome focusing on structural and functional aspects: I - gene expression profile of highly expressed phospholipases A(2)

Snake venom glands are a rich source of bioactive molecules such as peptides, proteins and enzymes that show important pharmacological activity leading to in local and systemic effects as pain, edema, bleeding and muscle necrosis. Most studies on pharmacologically active peptides and proteins from snake venoms have been concerned with isolation and structure elucidation through methods of classical biochemistry. As an attempt to examine the transcripts expressed in the venom gland of Bothrops jararacussu and to unveil the toxicological and pharmacological potential of its products at the molecular level, we generated 549 expressed sequence tags (ESTs) from a directional cDNA library. Sequences obtained from single-pass sequencing of randomly selected cDNA clones could be identified by similarities searches on existing databases, resulting in 197 sequences with significant similarity to phospholipase A(2) (PLA(2)), of which 83.2% were Lys49-PLA(2) homologs (BOJU-1), 0.1% were basic Asp49-PLA(2)s (BOJU-II) and 0.6% were acidic Asp49-PLA(2)s (BOJU-III). Adjoining this very abundant class of proteins we found 88 transcripts codifying for putative sequences of metalloproteases, which after clustering and assembling resulted in three full-length sequences: BOJUMET-I, BOJUMET-II and BOJUMET-III; as well as 25 transcripts related to C-type lectin like protein including a full-length cDNA of a putative galactose binding C-type lectin and a cluster of eight serine-proteases transcripts including a full-length cDNA of a putative serine protease. Among the full-length sequenced clones we identified a nerve growth factor (Bj-NGF) with 92% identity with a human NGF (NGHUBM) and an acidic phospholipase A2 (BthA-I-PLA(2)) displaying 85-93% identity with other snake venom toxins. Genetic distance among PLA(2)s from Bothrops species were evaluated by phylogenetic analysis. Furthermore, analysis of full-length putative Lys49-PLA(2) through molecular modeling showed conserved structural domains, allowing the characterization of those proteins as group II PLA(2)s. The constructed cDNA library provides molecular clones harboring sequences that can be used to probe directly the genetic material from gland venom of other snake species. Expression of complete cDNAs or their modified derivatives will be useful for elucidation of the structure-function relationships of these toxins and peptides of biotechnological interest. (C) 2004 Elsevier SAS. All rights reserved.

Analysis of the Nicotiana tabacum Stigma/Style Transcriptome Reveals Gene Expression Differences between Wet and Dry Stigma Species

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Impact of breed and sex on porcine endocrine transcriptome: a bayesian biometrical analysis

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Analysis of the workers head transcriptome of the Asian subterranean termite, Coptotermes gestroi

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The Vampirome: Transcriptome and proteome analysis of the principal and accessory submaxillary glands of the vampire bat Desmodus rotundus, a vector of human rabies

Vampire bats are notorious for being the sole mammals that strictly feed on fresh blood for their survival. While their saliva has been historically associated with anticoagulants, only one antihemostatic (plasminogen activator) has been molecularly and functionally characterized. Here, RNAs from both principal and accessory submaxillary (submandibular) salivary glands of Desmodus rotundus were extracted, and ~. 200. million reads were sequenced by Illumina. The principal gland was enriched with plasminogen activators with fibrinolytic properties, members of lipocalin and secretoglobin families, which bind prohemostatic prostaglandins, and endonucleases, which cleave neutrophil-derived procoagulant NETs. Anticoagulant (tissue factor pathway inhibitor, TFPI), vasodilators (PACAP and C-natriuretic peptide), and metalloproteases (ADAMTS-1) were also abundantly expressed. Members of the TSG-6 (anti-inflammatory), antigen 5/CRISP, and CCL28-like (antimicrobial) protein families were also sequenced. Apyrases (which remove platelet agonist ADP), phosphatases (which degrade procoagulant polyphosphates), and sphingomyelinase were found at lower transcriptional levels. Accessory glands were enriched with antimicrobials (lysozyme, defensin, lactotransferrin) and protease inhibitors (TIL-domain, cystatin, Kazal). Mucins, heme-oxygenase, and IgG chains were present in both glands. Proteome analysis by nano LC-MS/MS confirmed that several transcripts are expressed in the glands. The database presented herein is accessible online at http://exon.niaid.nih.gov/transcriptome/D_rotundus/Supplemental-web.xlsx. These results reveal that bat saliva emerges as a novel source of modulators of vascular biology. Biological significance: Vampire bat saliva emerges as a novel source of antihemostatics which modulate several aspects of vascular biology. © 2013.

The use of Open Reading frame ESTs (ORESTES) for analysis of the honey bee transcriptome

Abstract Background The ongoing efforts to sequence the honey bee genome require additional initiatives to define its transcriptome. Towards this end, we employed the Open Reading frame ESTs (ORESTES) strategy to generate profiles for the life cycle of Apis mellifera workers. Results Of the 5,021 ORESTES, 35.2% matched with previously deposited Apis ESTs. The analysis of the remaining sequences defined a set of putative orthologs whose majority had their best-match hits with Anopheles and Drosophila genes. CAP3 assembly of the Apis ORESTES with the already existing 15,500 Apis ESTs generated 3,408 contigs. BLASTX comparison of these contigs with protein sets of organisms representing distinct phylogenetic clades revealed a total of 1,629 contigs that Apis mellifera shares with different taxa. Most (41%) represent genes that are in common to all taxa, another 21% are shared between metazoans (Bilateria), and 16% are shared only within the Insecta clade. A set of 23 putative genes presented a best match with human genes, many of which encode factors related to cell signaling/signal transduction. 1,779 contigs (52%) did not match any known sequence. Applying a correction factor deduced from a parallel analysis performed with Drosophila melanogaster ORESTES, we estimate that approximately half of these no-match ESTs contigs (22%) should represent Apis-specific genes. Conclusions The versatile and cost-efficient ORESTES approach produced minilibraries for honey bee life cycle stages. Such information on central gene regions contributes to genome annotation and also lends itself to cross-transcriptome comparisons to reveal evolutionary trends in insect genomes.

Combined transcriptome and proteome analysis as a powerful approach to study genes under glucose repression in Bacillus subtilis

We used 2D protein gel electrophoresis and DNA microarray technologies to systematically analyze genes under glucose repression in Bacillus subtilis. In particular, we focused on genes expressed after the shift from glycolytic to gluconeogenic at the middle logarithmic phase of growth in a nutrient sporulation medium, which remained repressed by the addition of glucose. We also examined whether or not glucose repression of these genes was mediated by CcpA, the catabolite control protein of this bacterium. The wild-type and ccpA1 cells were grown with and without glucose, and their proteomes and transcriptomes were compared. 2D gel electrophoresis allowed us to identify 11 proteins, the synthesis of which was under glucose repression. Of these proteins, the synthesis of four (IolA, I, S and PckA) was under CcpA-independent control. Microarray analysis enabled us to detect 66 glucose-repressive genes, 22 of which (glmS, acoA, C, yisS, speD, gapB, pckA, yvdR, yxeF, iolA, B, C, D, E, F, G, H, I, J, R, S and yxbF ) were at least partially under CcpA-independent control. Furthermore, we found that CcpA and IolR, a repressor of the iol divergon, were involved in the glucose repression of the synthesis of inositol dehydrogenase encoded by iolG included in the above list. The CcpA-independent glucose repression of the iol genes appeared to be explained by inducer exclusion.

Identification and analysis of chromodomain-containing proteins encoded in the mouse transcriptome

The chromodomain is 40-50 amino acids in length and is conserved in a wide range of chromatic and regulatory proteins involved in chromatin remodeling. Chromodomain-containing proteins can be classified into families based on their broader characteristics, in particular the presence of other types of domains, and which correlate with different subclasses of the chromodomains themselves. Hidden Markov model (HMM)-generated profiles of different subclasses of chromodomains were used here to identify sequences encoding chromodomain-containing proteins in the mouse transcriptome and genome. A total of 36 different loci encoding proteins containing chromodomains, including 17 novel loci, were identified. Six of these loci (including three apparent pseudogenes, a novel HP1 ortholog, and two novel Msl-3 transcription factor-like proteins) are not present in the human genome, whereas the human genome contains four loci (two CDY orthologs and two apparent CDY pseuclogenes) that are not present in mouse. A number of these loci exhibit alternative splicing to produce different isoforms, including 43 novel variants, some of which lack the chromodomain. The likely functions of these proteins are discussed in relation to the known functions of other chromodomain-containing proteins within the same family.

Construction of representative transcript and protein sets of human, mouse, and rat as a platform for their transcriptome and proteome analysis

The number of mammalian transcripts identified by full-length cDNA projects and genome sequencing projects is increasing remarkably. Clustering them into a strictly nonredundant and comprehensive set provides a platform for functional analysis of the transcriptome and proteome, but the quality of the clustering and predictive usefulness have previously required manual curation to identify truncated transcripts and inappropriate clustering of closely related sequences. A Representative Transcript and Protein Sets (RTPS) pipeline was previously designed to identify the nonredundant and comprehensive set of mouse transcripts based on clustering of a large mouse full-length cDNA set (FANTOM2). Here we propose an alternative method that is more robust, requires less manual curation, and is applicable to other organisms in addition to mouse. RTPSs of human, mouse, and rat have been produced by this method and used for validation. Their comprehensiveness and quality are discussed by comparison with other clustering approaches. The RTPSs are available at ftp://fantom2.gsc.riken.go.jp/RTPS/. (C). 2004 Elsevier Inc. All rights reserved.