Deep sequencing of New World screw-worm transcripts to discover genes involved in insecticide resistance

Background: The New World screw-worm (NWS), Cochliomyia hominivorax, is one of the most important myiasis-causing flies, causing severe losses to the livestock industry. In its current geographical distribution, this species has been controlled by the application of insecticides, mainly organophosphate (OP) compounds, but a number of lineages have been identified that are resistant to such chemicals. Despite its economic importance, only limited genetic information is available for the NWS. Here, as a part of an effort to characterize the C. hominivorax genome and identify putative genes involved in insecticide resistance, we sampled its transcriptome by deep sequencing of polyadenylated transcripts using the 454 sequencing technology. Results: Deep sequencing on the 454 platform of three normalized libraries (larval, adult male and adult female) generated a total of 548,940 reads. Eighteen candidate genes coding for three metabolic detoxification enzyme families, cytochrome P450 monooxygenases, glutathione S transferases and carboxyl/cholinesterases were selected and gene expression levels were measured using quantitative real-time polymerase chain reaction (qRT-PCR). Of the investigated candidates, only one gene was expressed differently between control and resistant larvae with, at least, a 10-fold down-regulation in the resistant larvae. The presence of mutations in the acetylcholinesterase (target site) and carboxylesterase E3 genes was investigated and all of the resistant flies presented E3 mutations previously associated with insecticide resistance. Conclusions: Here, we provided the largest database of NWS expressed sequence tags that is an important resource, not only for further studies on the molecular basis of the OP resistance in NWS fly, but also for functional and comparative studies among Calliphoridae flies. Among our candidates, only one gene was found differentially expressed in resistant individuals, and its role on insecticide resistance should be further investigated. Furthermore, the absence of mutations in the OP target site and the high frequency of mutant carboxylesterase E3 indicate that metabolic resistance mechanisms have evolved predominantly in this species.

Alternative splicing enriched cDNA libraries identify breast cancer-associated transcripts

Background: Alternative splicing (AS) is a central mechanism in the generation of genomic complexity and is a major contributor to transcriptome and proteome diversity. Alterations of the splicing process can lead to deregulation of crucial cellular processes and have been associated with a large spectrum of human diseases. Cancer-associated transcripts are potential molecular markers and may contribute to the development of more accurate diagnostic and prognostic methods and also serve as therapeutic targets. Alternative splicing-enriched cDNA libraries have been used to explore the variability generated by alternative splicing. In this study, by combining the use of trapping heteroduplexes and RNA amplification, we developed a powerful approach that enables transcriptome-wide exploration of the AS repertoire for identifying AS variants associated with breast tumor cells modulated by ERBB2 (HER-2/neu) oncogene expression. Results: The human breast cell line (C5.2) and a pool of 5 ERBB2 over-expressing breast tumor samples were used independently for the construction of two AS-enriched libraries. In total, 2,048 partial cDNA sequences were obtained, revealing 214 alternative splicing sequence-enriched tags (ASSETs). A subset with 79 multiple exon ASSETs was compared to public databases and reported 138 different AS events. A high success rate of RT-PCR validation (94.5%) was obtained, and 2 novel AS events were identified. The influence of ERBB2-mediated expression on AS regulation was evaluated by capillary electrophoresis and probe-ligation approaches in two mammary cell lines (Hb4a and C5.2) expressing different levels of ERBB2. The relative expression balance between AS variants from 3 genes was differentially modulated by ERBB2 in this model system. Conclusions: In this study, we presented a method for exploring AS from any RNA source in a transcriptome-wide format, which can be directly easily adapted to next generation sequencers. We identified AS transcripts that were differently modulated by ERBB2-mediated expression and that can be tested as molecular markers for breast cancer. Such a methodology will be useful for completely deciphering the cancer cell transcriptome diversity resulting from AS and for finding more precise molecular markers.

Var transcription profiling of Plasmodium falciparum 3D7: assignment of cytoadherent phenotypes to dominant transcripts

Background: Cytoadherence of Plasmodium falciparum-infected red blood cells is mediated by var gene-encoded P. falciparum erythrocyte membrane protein-1 and host receptor preference depends in most cases on which of the 50-60 var genes per genome is expressed. Enrichment of phenotypically homogenous parasites by panning on receptor expressing cells is fundamental for the identification of the corresponding var transcript. Methods: P. falciparum 3D7 parasites were panned on several transfected CHO-cell lines and their var transcripts analysed by i) reverse transcription/PCR/cloning/sequencing using a universal DBL alpha specific oligonucleotide pair and ii) by reverse transcription followed by quantitative PCR using 57 different oligonucleotide pairs. Results: Each cytoadherence selected parasite line also adhered to untransfected CHO-745 cells and upregulation of the var gene PFD995/PFD1000c was consistently associated with cytoadherence to all but one CHO cell line. In addition, parasites panned on different CHO cell lines revealed candidate var genes which reproducibly associated to the respective cytoadherent phenotype. The transcription profile obtained by RT-PCR/cloning/sequencing differed significantly from that of RT-quantitative PCR. Conclusion: Transfected CHO cell lines are of limited use for the creation of monophenotypic cytoadherent parasite lines. Nevertheless, 3D7 parasites can be reproducibly selected for the transcription of different determined var genes without genetic manipulation. Most importantly, var transcription analysis by RT-PCR/cloning/sequencing may lead to erroneous interpretation of var transcription profiles.

Cloning and characterization of transcripts differentially expressed in the pulp of ripening papaya

Papaya (Carica papaya) is a relevant tropical crop and physico-chemical changes take place very quickly, as a consequence of activation of biochemical pathways by de nova synthesis of several proteins. Thus, in order to have information on the changes in gene expression in ripening papaya, transcripts from the pulp of unripe and ripe fruit were profiled by differential-display RT-PCR (DDRT-PCR). Seventy transcript derived fragments (TDFs) isolated from gels were re-amplified by PCR and differential expression of 40 papaya genes was confirmed by reverse northern blotting. Twenty-nine positively cloned TDFs were sequenced, and 17 were putatively identified by homology search. Ten of these genes were downregulated during ripening and UDP-glucose glucosyltransferase, alpha-2 importin, RNase L inhibitor-like protein, and a syntaxin protein were identified. Among the up-regulated genes there was a carboxylesterase, an integral membrane Yip1 family protein, a glycosyl hydrolase family-like protein and an endopolygalacturonase. Considering their relatedness to papaya quality, the fragments of genes potentially implicated in carbohydrate metabolism and pulp softening may be considered of interest for further studies. According to the results, differential display was a feasible approach to investigate differences in gene expression during fruit ripening, and can provide interesting information about those fruits whose genomic data is scarce, as is the case of papayas. (c) 2009 Elsevier B.V. All rights reserved.

Wolbachia neither induces nor suppresses transcripts encoding antimicrobial peptides

Wolbachia are intracellular maternally inherited microorganisms that are associated with reproductive abnormalities such as cytoplasmic incompatibility (CI), feminization and parthenogenesis in the various arthropod species they infect. Surveys indicate that these bacteria infect more than 16% of all insect species as well as isopods, mites and nematodes, making Wolbachia one of the most ubiquitous parasites yet described. However, nothing is known about the interactions of this bacterium with the host's immune system. We studied the expression of inducible antimicrobial markers in the adults of two Wolbachia infected insect species, Drosophila simulans and Aedes albopictus. The lack of available immune markers in the mosquito species led us to clone part of the defensin gene from this species, which was found to be very similar to the other mosquito defensins cloned from Anopheles gambiae and Aedes aegypti. Comparisons of the expression pattern of the antibacterial markers between Wolbachia-infected and cured lines, and also between bacteria-challenged and unchallenged adults indicated that Wolbachia does not either constitutively induce or suppress the transcription of these antibacterial genes. In addition, no difference in the transcription of these genes was found between double and single Wolbachia-infected strains or between strains in which Wolbachia has different tissue tropisms.

Gene expression profiling of papillary thyroid carcinoma identifies transcripts correlated with BRAF mutational status and lymph node metastasis

Purpose: To identify papillary thyroid carcinoma (PTC)-associated transcripts, we compared the gene expression profiles of three Serial Analysis of Gene Expression libraries generated from thyroid tumors and a normal thyroid tissue. Experimental Design: Selected transcripts were validated in a panel of 57 thyroid tumors using quantitative PCR (qPCR). An independent set of 71 paraffin-embedded sections was used for validation using immunohistochemical analysis. To determine if PTC-associated gene expression could predict lymph node involvement, a separate cohort of 130 primary PTC (54 metastatic and 76 nonmetastatic) was investigated. The BRAF(V600E) mutational status was compared with qPCR data to identify genes that might be regulated by abnormal BRAF/MEK/extracellular signal-regulated kinase signaling. Results: We identified and validated new PTC-associated transcripts. Three genes (CST6, CXCL14, and DHRS3) are strongly associated with PTC. Immunohistochemical analysis of CXCL14 confirmed the qPCR data and showed protein expression in PTC epithelial cells. We also observed that CST6, CXCL14, DHRS3, and SPP1 were associated with PTC lymph node metastasis, with CST6, CXCL14, and SPP1 being positively correlated with metastasis and DHRS3 being negatively correlated. Finally, we found a strong correlation between CST6 and CXCL14 expression and BRAF(V600E) mutational status, suggesting that these genes may be induced subsequently to BRAF activation and therefore may be downstream in the BRAF/MEK/extracellular signal-regulated kinase signaling pathway. Conclusion: CST6, CXCL14, DHRS3, and SPP1 may play a role in PTC pathogenesis and progression and are possible molecular targets for FTC therapy.

Characterization of bovine transcripts preferentially expressed in testis and with a putative role in spermatogenesis

Although the number of genes known to be associated with bovine spermatogenesis has increased in the past few years, regulation of this biological process remains poorly understood. Therefore, discovery of new male fertility genetic markers is of great value for assisted selection in commercially important cattle breeds, e.g., Nelore, that have delayed reproductive maturation and low fertility rates. The objective of the present study was to identify sequences associated with spermatogenesis that could be used as fertility markers. With RT-PCR, the following five transcripts preferentially expressed in adult testis were detected: TET(656) detected only in adult testis; TET(868) and TET(515) expressed preferentially in adult testis but also detected in fetal gonads of both sexes; and TET(456) and TET(262). expressed primarily in the testis, but also present in very low amounts in somatic tissues. Based on their homologies and expression profiles, we inferred that they had putative roles in spermatogenesis. Detection of sequences differentially expressed in testis, ovary, or both, was a useful approach for identifying new genes related to bovine spermatogenesis. The data reported here contributed to discovery of gene pathways involved in bovine spermatogenesis, with potential for prediction of fertility. (C) 2011 Elsevier Inc. All rights reserved.

Prognostic Impact of MYCN, DDX1, TrkA, and TrkC Gene Transcripts Expression in Neuroblastoma

Background. The aims of this study were to define the mRNA expression profiles of MYCN, DDX1, TrkA, and TrkC in biopsy tumor samples from 64 Brazilian patients with neuroblastomas of different risk stages and to correlate altered expression with prognostic values. Procedure. Patients were retrospectively classified into low- (n = 11), intermediate- (n = 18), and high-risk (n = 35) groups using standard criteria. The mRNA levels of the above genes were measured by quantitative real-time polymerase chain reaction. Univariate analyses were performed and survival curves were plotted by the Kaplan-Meier method. Results. Of the 64 patients, 53% were female and 62.5% were older than 18 months. The 5-year overall survival (OS) for the entire cohort was 40.3%, with inferior median OS in patients identified in the intermediate- and high-risk groups. A significant difference in OS with respect to TrkA mRNA expression was found for the high-risk group vs. either the low- or intermediate-risk groups (P < 0.01, log rank test). Within the intermediate-risk group, neuroblastoma patients with positive TrkA mRNA expression had better clinical outcomes than patients with no TrkA transcript expression (P = 0.004). Another difference in OS was only found between the intermediate- and high-risk groups (P < 0.027, log rank test). No significant correlation of mRNA expression and survival outcome could be detected for the MYCN, DDX1. Conclusions. Positive expression of TrkA mRNA may be a clinically useful addition to the current risk classification system, allowing the identification of NB tumors with favorable prognosis. Pediatr Blood Cancer 2011; 56: 749-756. (c) 2010 Wiley-Liss, Inc.

Isolation of transcripts overexpressed in the human pathogen Trichophyton rubrum grown in lipid as carbon source

Trichophyton rubrum is the most common etiological agent of human dermatophytosis. Despite the incidence and medical importance of this dermatophyte, little is known about the mechanisms of host invasion and pathogenicity. Host invasion depends on the adaptive cellular responses of the pathogen that allow it to penetrate the skin layers, which are mainly composed of proteins and lipids. In this study, we used suppression subtractive hybridization to identify transcripts over-expressed in T rubrum cultured in lipid as carbon source. Among the subtractive cDNA clones isolated, 85 clones were positively screened by cDNA array dot blotting and were sequenced. The putative proteins encoded by the isolated transcripts showed similarities to fungal proteins involved in metabolism, signaling, defense, and virulence, such as the MDR/ABC transporter, glucan 1,3-beta-glucosidase, chitin synthase B, copper-sulfate-regulated protein, and serine/threonine phosphatase (calcineurin A). These results provide the first molecular insight into the genes differentially expressed during the adaptation of T. rubrum to a lipidic carbon source.

Transcripts of the MHM region on the chicken Z chromosome accumulate as non-coding RNA in the nucleus of female cells adjacent to the DMRT1 locus

The male hypermethylated (MHM) region, located near the middle of the short arm of the Z chromosome of chickens, consists of approximately 210 tandem repeats of a BamHI 2.2-kb sequence unit. Cytosines of the CpG dinucleotides of this region are extensively methylated on the two Z chromosomes in the male but much less methylated on the single Z chromosome in the female. The state of methylation of the MHM region is established after fertilization by about the 1-day embryonic stage. The MHM region is transcribed only in the female from the particular strand into heterogeneous, high molecular-mass, non-coding RNA, which is accumulated at the site of transcription, adjacent to the DMRT1 locus, in the nucleus. The transcriptional silence of the MHM region in the male is most likely caused by the CpG methylation, since treatment of the male embryonic fibroblasts with 5-azacytidine results in hypo-methylation and active transcription of this region. In ZZW triploid chickens, MHM regions are hypomethylated and transcribed on the two Z chromosomes, whereas MHM regions are hypermethylated and transcriptionally inactive on the three Z chromosomes in ZZZ triploid chickens, suggesting a possible role of the W chromosome on the state of the MHM region.

Quantitation of alternatively spliced NMDA receptor NR1 isoform mRNA transcripts in human brain by competitive RT-PCR

The N-methyl-D-aspartate (NMDA)-selective subtype of ionotropic glutamate receptor is of importance in neuronal differentiation and synapse consolidation, activity-dependent forms of synaptic plasticity, and excitatory amino acid-mediated neuronal toxicity [Neurosci. Res. Program, Bull. 19 (1981) 1; Lab. Invest. 68 (1993) 372]. NMDA receptors exist in vivo as tetrameric or pentameric complexes comprising proteins from two families of homologous subunits, designated NR1 and NR2(A-D) [Biochem. Biophys. Res. Commun. 185 (1992) 826]. The gene coding for the human NR1 subunit (hNR1) is composed of 21 exons, three of which (4, 20 and 21) can be differentially spliced to generate a total of eight distinct subunit variants. We detail here a competitive RT-PCR (cRT-PCR) protocol to quantify endogenous levels of hNR1 splice variants in autopsied human brain. Quantitation of each hNR1 splice variant is performed using standard curve methodology in which a known amount of synthetic ribonucleic acid competitor (internal standard) is co-amplified against total RNA. This method can be used for the quantitation of hNR1 mRNA levels in response to acute or chronic disease states, in particular in the glutamatergic-associated neuronal loss observed in Alzheimer's disease [J. Neurochem. 78 (2001) 175]. Furthermore, alterations in hNR1 mRNA expression may be reflected at the translational level, resulting in functional changes in the NMDA receptor. (C) 2003 Elsevier Science B.V. All rights reserved.

Rapid turnover of Plasmodium falciparum var gene transcripts and genotypes during natural non-symptomatic infections

The var genes of Plasmodium falciparum code for the antigenically variant erythrocyte membrane proteins 1 (PfEMP1), a major factor for cytoadherence and immune escape of the parasite. Herein, we analyzed the var gene transcript turnover in two ongoing, non-symptomatic infections at sequential time points during two weeks. The number of different circulating genomes was estimated by microsatellite analyses. In both infections, we observed a rapid turnover of plasmodial genotypes and var transcripts. The rapidly changing repertoire of var transcripts could have been caused either by swift elimination of circulating var-transcribing parasites stemming from different or identical genetic backgrounds, or by accelerated switching of var gene transcription itself.

Analysis of the intranuclear life of nonsense transcripts

Dissertação para obtenção do Grau de Doutor em Biologia, Especialidade de Biologia Molecular

Nineteen additional unpredicted transcripts from human chromosome 21.

The identification of all human chromosome 21 (HC21) genes is a necessary step in understanding the molecular pathogenesis of trisomy 21 (Down syndrome). The first analysis of the sequence of 21q included 127 previously characterized genes and predicted an additional 98 novel anonymous genes. Recently we evaluated the quality of this annotation by characterizing a set of HC21 open reading frames (C21orfs) identified by mapping spliced expressed sequence tags (ESTs) and predicted genes (PREDs), identified only in silico. This study underscored the limitations of in silico-only gene prediction, as many PREDs were incorrectly predicted. To refine the HC21 annotation, we have developed a reliable algorithm to extract and stringently map sequences that contain bona fide 3' transcript ends to the genome. We then created a specific 21q graphical display allowing an integrated view of the data that incorporates new ESTs as well as features such as CpG islands, repeats, and gene predictions. Using these tools we identified 27 new putative genes. To validate these, we sequenced previously cloned cDNAs and carried out RT-PCR, 5'- and 3'-RACE procedures, and comparative mapping. These approaches substantiated 19 new transcripts, thus increasing the HC21 gene count by 9.5%. These transcripts were likely not previously identified because they are small and encode small proteins. We also identified four transcriptional units that are spliced but contain no obvious open reading frame. The HC21 data presented here further emphasize that current gene prediction algorithms miss a substantial number of transcripts that nevertheless can be identified using a combination of experimental approaches and multiple refined algorithms.

Differential expression of notch signaling-related transcripts accompanies Pro-thymocyte proliferation and phenotype transition induced by epidermal growth factor plus insulin in fetal thymus organ cultures

Thymus regression upon stressing stimuli, such as infectious diseases, is followed by organ reconstitution, paralleling its development in ontogeny. A narrow window of thymus development was here studied, encompassing the pro-T lymphoid precursor expansion during specification stages, by the use of epidermal growth factor plus insulin (INS) in murine fetal thymus organ cultures. Aiming to disclose signaling pathways related to these stages, cultured thymus lobes had their RNA extracted, for the search of transcripts differentially expressed using RNAse protection assays and reverse transcriptase-polymerase chain reactions. We found no difference that could explain INS-driven thymocyte growth, in the pattern of transcripts for death/proliferation mediators, or for a series of growth factor receptors and transcriptional regulators known as essential for thymus development. Thymocyte suspensions from cultured lobes, stained for phenotype analysis by fluorescence activated cell sorting, showed a decreased staining for Notch1 protein at cell surfaces upon INS addition. We analyzed the expression of Notch-related elements, and observed the recruitment of a specific set of transcripts simultaneous and compatible with INS-driven thymocyte growth, namely, transcripts for Notch3, for its ligand Jagged2, and for Deltex1, a mediator of a poorly characterized alternative pathway downstream of the Notch receptor.