Modularity and evolutionary constraints in a baculovirus gene regulatory network

Abstract Background The structure of regulatory networks remains an open question in our understanding of complex biological systems. Interactions during complete viral life cycles present unique opportunities to understand how host-parasite network take shape and behave. The Anticarsia gemmatalis multiple nucleopolyhedrovirus (AgMNPV) is a large double-stranded DNA virus, whose genome may encode for 152 open reading frames (ORFs). Here we present the analysis of the ordered cascade of the AgMNPV gene expression. Results We observed an earlier onset of the expression than previously reported for other baculoviruses, especially for genes involved in DNA replication. Most ORFs were expressed at higher levels in a more permissive host cell line. Genes with more than one copy in the genome had distinct expression profiles, which could indicate the acquisition of new functionalities. The transcription gene regulatory network (GRN) for 149 ORFs had a modular topology comprising five communities of highly interconnected nodes that separated key genes that are functionally related on different communities, possibly maximizing redundancy and GRN robustness by compartmentalization of important functions. Core conserved functions showed expression synchronicity, distinct GRN features and significantly less genetic diversity, consistent with evolutionary constraints imposed in key elements of biological systems. This reduced genetic diversity also had a positive correlation with the importance of the gene in our estimated GRN, supporting a relationship between phylogenetic data of baculovirus genes and network features inferred from expression data. We also observed that gene arrangement in overlapping transcripts was conserved among related baculoviruses, suggesting a principle of genome organization. Conclusions Albeit with a reduced number of nodes (149), the AgMNPV GRN had a topology and key characteristics similar to those observed in complex cellular organisms, which indicates that modularity may be a general feature of biological gene regulatory networks.

Crystal structure of the SarR protein from Staphylococcus aureus

The expression of virulence determinants in Staphylococcus aureus is controlled by global regulatory loci (e.g., sarA and agr). The sar (Staphylococcus accessory regulator) locus is composed of three overlapping transcripts (sarA P1, P3, and P2, transcripts initiated from the P1, P3, and P2 promoters, respectively), all encoding the 124-aa SarA protein. The level of SarA, the major regulatory protein, is partially controlled by the differential activation of the sarA promoters. We previously partially purified a 13.6-kDa protein, designated SarR, that binds to the sarA promoter region to down-modulate sarA transcription from the P1 promoter and subsequently SarA expression. SarR shares sequence similarity to SarA, and another SarA homolog, SarS. Here we report the 2.3 Å-resolution x-ray crystal structure of the dimeric SarR-MBP (maltose binding protein) fusion protein. The structure reveals that the SarR protein not only has a classic helix–turn–helix module for DNA binding at the major grooves, but also has an additional loop region involved in DNA recognition at the minor grooves. This interaction mode could represent a new functional class of the “winged helix” family. The dimeric SarR structure could accommodate an unusually long stretch of ≈27 nucleotides with two or four bending points along the course, which could lead to the bending of DNA by 90° or more, similar to that seen in the catabolite activator protein (CAP)–DNA complex. The structure also demonstrates the molecular basis for the stable dimerization of the SarR monomers and possible motifs for interaction with other proteins.

Ribosomal and Immune Transcripts Associate with Relapse in Acquired ADAMTS13-Deficient Thrombotic Thrombocytopenic Purpura.

Approximately 40% of patients who survive acute episodes of thrombotic thrombocytopenic purpura (TTP) associated with severe acquired ADAMTS13 deficiency experience one or more relapses. Risk factors for relapse other than severe ADAMTS13 deficiency and ADAMTS13 autoantibodies are unknown. ADAMTS13 autoantibodies, TTP episodes following infection or type I interferon treatment and reported ensuing systemic lupus erythematosus in some patients suggest immune dysregulation. This cross-sectional study asked whether autoantibodies against RNA-binding proteins or peripheral blood gene expression profiles measured during remission are associated with history of prior relapse in acquired ADAMTS13-deficient TTP. Peripheral blood from 38 well-characterized patients with autoimmune ADAMTS13-deficient TTP in remission was examined for autoantibodies and global gene expression. A subset of TTP patients (9 patients, 24%) exhibited a peripheral blood gene signature composed of elevated ribosomal transcripts that associated with prior relapse. A non-overlapping subset of TTP patients (9 patients, 24%) displayed a peripheral blood type I interferon gene signature that associated with autoantibodies to RNA-binding proteins but not with history of relapse. Patients who had relapsed bimodally expressed higher HLA transcript levels independently of ribosomal transcripts. Presence of any one potential risk factor (ribosomal gene signature, elevated HLA-DRB1, elevated HLA-DRB5) associated with relapse (OR = 38.4; p = 0.0002) more closely than any factor alone or all factors together. Levels of immune transcripts typical of natural killer (NK) and T lymphocytes positively correlated with ribosomal gene expression and number of prior episodes but not with time since the most recent episode. Flow cytometry confirmed elevated expression of cell surface markers encoded by these transcripts on T and/or NK cell subsets of patients who had relapsed. These data associate elevated ribosomal and immune transcripts with relapse history in acquired, ADAMTS13-deficient TTP.

Rescue of embryonic lethality in mdmx null mice by loss of p53 suggests a non-overlapping pathway with MDM2 to regulate p53

The tumor suppressor p53 is mutated in over 50% of human sporadic tumors originating from diverse tissues. p53 responds to DNA damage and cell stress by activating the transcription of a variety of target genes, the protein products of which then initiate either growth arrest or apoptosis. ^ A p53 target with a particularly intriguing function is the oncogene MDM2. MDM2 functions, in part, by binding to and inhibiting p53's activity. Overexpression of MDM2, by gene amplification, has been found in 30% of human sarcomas harboring a wild type p53, indicating that an increase in MDM2 levels is sufficient for p53 inactivation. Mice carrying a homozygous null allele for mdm2 exhibit an early embryonic lethality that is completely rescued in a p53-null background. These data indicate that MDM2's only critical function in early mouse embryogenesis is the negative regulation of p53. ^ The mdmx gene is the first additional member of the mdm2 gene family to be isolated. MDMX, like MDM2, contains a RING-finger domain, ATP binding domain and a p53 binding domain, which retains the ability to bind and inhibit p53 transactivation in vitro. However, mdmx does not appear to be transcriptionally regulated by p53. We have cloned and characterized the murine mdmx genomic locus from a mouse 129 genomic library. The mdmx gene contains 11 exons, spans approximately 37 Kb of DNA, and is located on mouse chromosome 1. The genomic organization of the mdmx gene is identical to that of mdm2 except at the 5′ end of the gene near the p53 responsive element. Northern expression analysis of mdmx transcripts during mouse embryogenesis and in adult tissues revealed constitutive and ubiquitous expression throughout adult tissues and embryonic development. To determine the in vivo function of MDMX, mice carrying a null allele of mdmx have been generated. Mdmx homozygous null mice are early embryonic lethal. Mdmx null mice do not develop beyond 9.5 dpc and can be discerned by gross dissection as early as 7.5 dpc. Utilizing TUNEL and BrdU assays on 7.5 dpc histological sections we have determined that the mutant embryos are dying due to increased levels of growth arrest, but not apoptosis. Surprisingly, Mdmx homozygous null mice are viable in a p53 null background, indicating that MDMX is also very important in the negative regulation of p53. ^

Multiple imprinted sense and antisense transcripts, differential methylation and tandem repeats in a putative imprinting control region upstream of mouse Igf2

The mouse insulin-like growth factor 2 (Igf2) locus is a complex genomic region that produces multiple transcripts from alternative promoters. Expression at this locus is regulated by parental imprinting. However, despite the existence of putative imprinting control elements in the Igf2 upstream region, imprinted transcriptional repression is abolished by null mutations at the linked H19 locus. To clarify the extent to which the Igf2 upstream region contains autonomous imprinting control elements we have performed functional and comparative analyses of the region in the mouse and human. Here we report the existence of multiple, overlapping imprinted (maternally repressed) sense and antisense transcripts that are associated with a tandem repeat in the mouse Igf2 upstream region. Regions flanking the repeat exhibit tissue-specific parental allelic methylation patterns, suggesting the existence of tissue-specific control elements in the upstream region. Studies in H19 null mice indicate that both parental allelic methylation and monoallelic expression of the upstream transcripts depends on an intact H19 gene acting in cis. The homologous region in human IGF2 is structurally conserved, with the significant exception that it does not contain a tandem repeat. Our results support the proposal that tandem repeats act to target methylation to imprinted genetic loci.

From hematopoiesis to neuropoiesis: Evidence of overlapping genetic programs

It is reasonable to propose that gene expression profiles of purified stem cells could give clues for the molecular mechanisms of stem cell behavior. We took advantage of cDNA subtraction to identify a set of genes selectively expressed in mouse adult hematopoietic stem cells (HSC) as opposed to bone marrow (BM). Analysis of HSC-enriched genes revealed several key regulatory gene candidates, including two novel seven transmembrane (7TM) receptors. Furthermore, by using cDNA microarray techniques we found a large set of HSC-enriched genes that are expressed in mouse neurospheres (a population greatly enriched for neural progenitor cells), but not present in terminally differentiated neural cells. In situ hybridization demonstrated that many of them, including one HSC-enriched 7TM receptor, were selectively expressed in the germinal zones of fetal and adult brain, the regions harboring mouse neural stem cells. We propose that at least some of the transcripts that are selectively and commonly expressed in two or more types of stem cells define a functionally conserved group of genes evolved to participate in basic stem cell functions, including stem cell self-renewal.

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.

Investigating Giant (Gt) Repression in the Formation of Partially Overlapping Pair-Rule Stripes

Drosophila pair-rule genes are expressed in striped patterns with a precise order of overlap between stripes of different genes. We investigated the role of Giant (Gt) in the regulation of even-skipped, hairy, runt, and fushi tarazu stripes formed in the vicinity of Gt expression domains. In gt null embryos, specific stripes of eve, h, run, and ftz are disrupted. With an ectopic expression system, we verified that stripes affected in the mutant are also repressed. Simultaneously hybridizing gt misxpressing embryos with two pair-rule gene probes, we were able to distinguish differences in the repression of pairs of stripes that overlap extensively. Together, our results showed Gt repression roles in the regulation of two groups of partially overlapping stripes and that Gt morphogen activity is part of the mechanism responsible for the differential positioning of these stripes borders. We discuss the possibility that other factors regulate Gt stripe targets as well. Developmental Dynamics 239:2989-2999, 2010. (C) 2010 Wiley-Liss, Inc.

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.

Overlapping motifs in the membrane-proximal region of cytokine receptor accessory and signaling subunits

The membrane-proximal cytoplasmic region of cytokine receptors (CRs) is highly conserved and essential for receptor activation. In particular this region is essential for the activation of members of the Janus family of protein kinases (JAK) which results in initiation of receptor signaling. We have examined the sequence of this region in a number of CR signaling and accessory subunits with a view to better delineating motifs that play an important role in initiating receptor activity. Here, we have delineated two distinct proline-rich motifs in the membrane-proximal domains of cytokine receptors. Their configuration and distribution among CR subunits strongly suggest a model in which the two motifs act in a concerted manner to induce full receptor and JAK activation. (C) 2004 Elsevier Ltd. All rights reserved.

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.