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.

Bioremediation with White-Rot Fungi at Fisherville Mill: Analyses of Gene Expression and Number 6 Fuel Oil Degradation

Extracellular enzymes that white-rot fungi secrete during lignin decay have been proposed as promising agents for oxidizing pollutants. We investigated the abilities of the white-rot fungi Punctularia strigosozonata, Irpex lacteus, Trichaptum biforme, Phlebia radiata, Trametes versicolor, and Pleurotus ostreatus to degrade Number 6 fuel oil in wood sawdust cultures. Our goals are to advise bioremediation efforts at a brownfield redevelopment site on the Blackstone River in Grafton, Massachusetts and to contribute to the understanding of decay mechanisms in white-rot fungi. All species tested degraded a C10 alkane. When cultivated for 6 months, Irpex lacteus, T. biforme, P. radiata, T. versicolor and P. ostreatus also degraded a C14 alkane and the polycyclic aromatic hydrocarbon phenanthrene. Gene expression analyses of P. strigosozonata indicate differential gene expression in the presence of Number 6 oil and on pine and aspen sawdust.

Gene number in an invertebrate chordate, Ciona intestinalis

Gene number can be considered a pragmatic measure of biological complexity, but reliable data is scarce. Estimates for vertebrates are 50-100,000 genes per haploid genome, whereas invertebrate estimates fall below 25,000. We wished to test the hypothesis that the origin of vertebrates coincided with extensive gene creation. A prediction is that gene number will differ sharply between invertebrate and vertebrate members of the chordate phylum. A gene number estimation method requiring limited sequence sampling of genomic DNA was developed and validated by using data for Caenorhabditis elegans. Using the method, we estimated that the invertebrate chordate Ciona intestinalis has 15,500 protein-coding genes (±3,700). This number is significantly lower than gene numbers of vertebrate chordates, but similar to those of invertebrates in distantly related phyla. The data indicate that evolution of vertebrates was accompanied by a dramatic increase in protein-coding capacity of the genome.

Cloning and characterization of the Flavobacterium johnsoniae (Cytophaga johnsonae) gliding motility gene, gldA

The mechanism of bacterial gliding motility (active movement over surfaces without the aid of flagella) is not known. A large number of nonmotile mutants of the gliding bacterium Flavobacterium johnsoniae (Cytophaga johnsonae) have been previously isolated, and genetic techniques to analyze these mutants have recently been developed. We complemented a nonmotile mutant of F. johnsoniae (UW102-09) with a library of wild-type DNA by using the shuttle cosmid pCP17. The complementing plasmid (pCP100) contained an insert of 13 kbp, and restored motility to 4 of 61 independently isolated nonmotile mutants. A 1.3-kbp fragment that encompassed a single ORF, gldA, complemented all four mutants. Disruption of the chromosomal copy of gldA in wild-type F. johnsoniae UW101 eliminated gliding motility. The predicted protein produced by gldA has strong sequence similarity to ATP binding cassette transport proteins.

Thermostable and site-specific DNA binding of the gene product ORF56 from the Sulfolobus islandicus plasmid pRN1, a putative archael plasmid copy control protein

There is still a lack of information on the specific characteristics of DNA-binding proteins from hyperthermophiles. Here we report on the product of the gene orf56 from plasmid pRN1 of the acidophilic and thermophilic archaeon Sulfolobus islandicus. orf56 has not been characterised yet but low sequence similarily to several eubacterial plasmid-encoded genes suggests that this 6.5 kDa protein is a sequence-specific DNA-binding protein. The DNA-binding properties of ORF56, expressed in Escherichia coli, have been investigated by EMSA experiments and by fluorescence anisotropy measurements. Recombinant ORF56 binds to double-stranded DNA, specifically to an inverted repeat located within the promoter of orf56. Binding to this site could down-regulate transcription of the orf56 gene and also of the overlapping orf904 gene, encoding the putative initiator protein of plasmid replication. By gel filtration and chemical crosslinking we have shown that ORF56 is a dimeric protein. Stoichiometric fluorescence anisotropy titrations further indicate that ORF56 binds as a tetramer to the inverted repeat of its target binding site. CD spectroscopy points to a significant increase in ordered secondary structure of ORF56 upon binding DNA. ORF56 binds without apparent cooperativity to its target DNA with a dissociation constant in the nanomolar range. Quantitative analysis of binding isotherms performed at various salt concentrations and at different temperatures indicates that approximately seven ions are released upon complex formation and that complex formation is accompanied by a change in heat capacity of –6.2 kJ/mol. Furthermore, recombinant ORF56 proved to be highly thermostable and is able to bind DNA up to 85°C.

Functional genomics: Probing plant gene function and expression with transposons

Transposable elements provide a convenient and flexible means to disrupt plant genes, so allowing their function to be assessed. By engineering transposons to carry reporter genes and regulatory signals, the expression of target genes can be monitored and to some extent manipulated. Two strategies for using transposons to assess gene function are outlined here: First, the PCR can be used to identify plants that carry insertions into specific genes from among pools of heavily mutagenized individuals (site-selected transposon mutagenesis). This method requires that high copy transposons be used and that a relatively large number of reactions be performed to identify insertions into genes of interest. Second, a large library of plants, each carrying a unique insertion, can be generated. Each insertion site then can be amplified and sequenced systematically. These two methods have been demonstrated in maize, Arabidopsis, and other plant species, and the relative merits of each are discussed in the context of plant genome research.

Disruption of the Cbfa2 gene causes necrosis and hemorrhaging in the central nervous system and blocks definitive hematopoiesis.

The CBFA2 (AML1) gene encodes a DNA-binding subunit of the heterodimeric core-binding factor. The CBFA2 gene is disrupted by the (8;21), (3;21), and (12;21) chromosomal translocations associated with leukemias and myelodysplasias in humans. Mice lacking a CBF alpha 2 protein capable of binding DNA die between embryonic days 11.5 and 12.5 due to hemorrhaging in the central nervous system (CNS), at the nerve/CNS interfaces of cranial and spinal nerves, and in somitic/intersomitic regions along the presumptive spinal cord. Hemorrhaging is preceded by symmetric, bilateral necrosis in these regions. Definitive erythropoiesis and myelopoiesis do not occur in Cbfa2-deficient embryos, and disruption of one copy of the Cbfa2 gene significantly reduces the number of progenitors for erythroid and myeloid cells.

Multiplex genotype determination at a large number of gene loci.

To facilitate large-scale genotype analysis, an efficient PCR-based multiplex approach has been developed. For simultaneously amplifying the target sequences at a large number of genetic loci, locus-specific primers containing 5' universal tails are used. Attaching the universal tails to the target sequences in the initial PCR steps allows replacement of all specific primers with a pair of primers identical to the universal tails and converts the multiplex amplification into "uniplex." Simultaneous amplification of 26 genetic loci with this approach is described. The multiplex amplification can be coupled with genotype determination. By incorporating a single-base mismatch between a primer and the template into the target sequences, a polymorphic site can be converted into a desirable restriction fragment length polymorphism when it is necessary. In this way, the allelic PCR products for the polymorphic loci can be discriminated by gel electrophoresis after restriction enzyme digestion. In this study, 32 loci were typed in such a multiplex way.

The penicillin gene cluster is amplified in tandem repeats linked by conserved hexanucleotide sequences.

The penicillin biosynthetic genes (pcbAB, pcbC, penDE) of Penicillium chrysogenum AS-P-78 were located in a 106.5-kb DNA region that is amplified in tandem repeats (five or six copies) linked by conserved TTTACA sequences. The wild-type strains P. chrysogenum NRRL 1951 and Penicillium notatum ATCC 9478 (Fleming's isolate) contain a single copy of the 106.5-kb region. This region was bordered by the same TTTACA hexanucleotide found between tandem repeats in strain AS-P-78. A penicillin overproducer strain, P. chrysogenum E1, contains a large number of copies in tandem of a 57.9-kb DNA fragment, linked by the same hexanucleotide or its reverse complementary TGTAAA sequence. The deletion mutant P. chrysogenum npe10 showed a deletion of 57.9 kb that corresponds exactly to the DNA fragment that is amplified in E1. The conserved hexanucleotide sequence was reconstituted at the deletion site. The amplification has occurred within a single chromosome (chromosome I). The tandem reiteration and deletion appear to arise by mutation-induced site-specific recombination at the conserved hexanucleotide sequences.

Extraordinary number of gene rearrangements in the mitochondrial genomes of lice (Phthiraptera : Insecta)

The arrangement of genes in the mitochondrial (mt) genomes of most insects is the same, or near-identical, to that inferred to be ancestral for insects. We sequenced the entire mt genome of the small pigeon louse, Campanulotes bidentatus compar, and part of the mt genomes of nine other species of lice. These species were from six families and the three main suborders of the order Phthiraptera. There was no variation in gene arrangement among species within a family but there was much variation in gene arrangement among the three suborders of lice. There has been an extraordinary number of gene rearrangements in the mitochondrial genomes of lice!

The upstream Variable Number Tandem Repeat polymorphism of the monoamine oxidase type A gene influences trigeminal pain-related evoked responses

Monoamines have an important role in neural plasticity, a key factor in cortical pain processing that promotes changes in neuronal network connectivity. Monoamine oxidase type A (MAOA) is an enzyme that, due to its modulating role in monoaminergic activity, could play a role in cortical pain processing. The X-linked MAOA gene is characterized by an allelic variant of length, the MAOA upstream Variable Number Tandem Repeat (MAOA-uVNTR) region polymorphism. Two allelic variants of this gene are known, the high-activity MAOA (HAM) and low-activity MAOA (LAM). We investigated the role of MAOA-uVNTR in cortical pain processing in a group of healthy individuals measured by the trigeminal electric pain-related evoked potential (tPREP) elicited by repeated painful stimulation. A group of healthy volunteers was genotyped to detect MAOA-uVNTR polymorphism. Electrical tPREPs were recorded by stimulating the right supraorbital nerve with a concentric electrode. The N2 and P2 component amplitude and latency as well as the N2-P2 inter-peak amplitude were measured. The recording was divided into three blocks, each containing 10 consecutive stimuli and the N2-P2 amplitude was compared between blocks. Of the 67 volunteers, 37 were HAM and 30 were LAM. HAM subjects differed from LAM subjects in terms of amplitude of the grand-averaged and first-block N2-P2 responses (HAM>LAM). The N2-P2 amplitude decreased between the first and third block in HAM subjects but not LAM subjects. The MAOA-uVNTR polymorphism seemed to influence the brain response in a repeated tPREP paradigm and suggested a role of the MAOA as a modulator of neural plasticity related to cortical pain processing. Monoamines have an important role in neural plasticity, a key factor in cortical pain processing that promotes changes in neuronal network connectivity. Monoamine oxidase type A (MAOA) is an enzyme that, due to its modulating role in monoaminergic activity, could play a role in cortical pain processing. © 2014 Federation of European Neuroscience Societies and John Wiley & Sons Ltd.

The mitochondrial genome of the stingless bee Melipona bicolor (Hymenoptera, Apidae, Meliponini): sequence, gene organization and a unique tRNA translocation event conserved across the tribe Meliponini

At present a complete mtDNA sequence has been reported for only two hymenopterans, the Old World honey bee, Apis mellifera and the sawfly Perga condei. Among the bee group, the tribe Meliponini (stingless bees) has some distinction due to its Pantropical distribution, great number of species and large importance as main pollinators in several ecosystems, including the Brazilian rain forest. However few molecular studies have been conducted on this group of bees and few sequence data from mitochondrial genomes have been described. In this project, we PCR amplified and sequenced 78% of the mitochondrial genome of the stingless bee Melipona bicolor (Apidae, Meliponini). The sequenced region contains all of the 13 mitochondrial protein-coding genes, 18 of 22 tRNA genes, and both rRNA genes (one of them was partially sequenced). We also report the genome organization (gene content and order), gene translation, genetic code, and other molecular features, such as base frequencies, codon usage, gene initiation and termination. We compare these characteristics of M. bicolor to those of the mitochondrial genome of A. mellifera and other insects. A highly biased A+T content is a typical characteristic of the A. mellifera mitochondrial genome and it was even more extreme in that of M. bicolor. Length and compositional differences between M. bicolor and A. mellifera genes were detected and the gene order was compared. Eleven tRNA gene translocations were observed between these two species. This latter finding was surprising, considering the taxonomic proximity of these two bee tribes. The tRNA Lys gene translocation was investigated within Meliponini and showed high conservation across the Pantropical range of the tribe.

O papel do polimorfismo funcional VNTR da região promotora do gene MAOA nos transtornos psiquiátricos

INTRODUÇÃO: Muitos estudos têm investigado a associação do polimorfismo VNTR (número variável de repetições em série) localizado na região promotora do gene da enzima monoamina oxidase A (MAOA) com alterações no comportamento humano e em diversos transtornos psiquiátricos. OBJETIVO: O objetivo do presente trabalho foi revisar a literatura sobre a participação desse polimorfismo funcional na modulação do comportamento humano para o desenvolvimento dos transtornos psiquiátricos. MÉTODO: A pesquisa foi realizada na literatura em inglês, de janeiro de 1998 a junho de 2009, disponível no Medline, Embase, Web of Science e na base de dados PsycInfo, utilizando os seguintes termos: "MAOA e comportamento humano" e "MAOA e psiquiatria". RESULTADOS: Foram encontrados 3.873 estudos. Desses, 109 foram selecionados e incluídos na revisão. Encontrou-se associação de alelos de baixa atividade do VNTR com transtorno de personalidade antissocial, transtorno de conduta, transtorno de déficit de atenção e hiperatividade, jogo patológico e dependência de substâncias. Alelos da alta atividade da MAOA foram associados a depressão, ansiedade, neuroticismo e anorexia nervosa. Não se encontrou associação entre polimorfismos da MAOA e esquizofrenia e transtorno bipolar. CONCLUSÃO: Os principais achados dão suporte ao papel do polimorfismo VNTR da região promotora do gene da MAOA em alguns transtornos psiquiátricos, apesar das divergências encontradas devidas às dificuldades metodológicas de estudos em genética. De modo geral, os estudos associam os alelos de baixa atividade da MAOA com comportamentos impulsivos e agressivos ("comportamentos hiperativos"), enquanto os alelos de alta atividade do gene são mais associados a "comportamentos hipoativos".

The full Bayesian significance test for mixture models: results in gene expression clustering

Gene clustering is a useful exploratory technique to group together genes with similar expression levels under distinct cell cycle phases or distinct conditions. It helps the biologist to identify potentially meaningful relationships between genes. In this study, we propose a clustering method based on multivariate normal mixture models, where the number of clusters is predicted via sequential hypothesis tests: at each step, the method considers a mixture model of m components (m = 2 in the first step) and tests if in fact it should be m - 1. If the hypothesis is rejected, m is increased and a new test is carried out. The method continues (increasing m) until the hypothesis is accepted. The theoretical core of the method is the full Bayesian significance test, an intuitive Bayesian approach, which needs no model complexity penalization nor positive probabilities for sharp hypotheses. Numerical experiments were based on a cDNA microarray dataset consisting of expression levels of 205 genes belonging to four functional categories, for 10 distinct strains of Saccharomyces cerevisiae. To analyze the method's sensitivity to data dimension, we performed principal components analysis on the original dataset and predicted the number of classes using 2 to 10 principal components. Compared to Mclust (model-based clustering), our method shows more consistent results.

Chicken skeletal muscle-associated macroarray for gene discovery

Macro- and microarrays are well-established technologies to determine gene functions through repeated measurements of transcript abundance. We constructed a chicken skeletal muscle-associated array based on a muscle-specific EST database, which was used to generate a tissue expression dataset of similar to 4500 chicken genes across 5 adult tissues (skeletal muscle, heart, liver, brain, and skin). Only a small number of ESTs were sufficiently well characterized by BLAST searches to determine their probable cellular functions. Evidence of a particular tissue-characteristic expression can be considered an indication that the transcript is likely to be functionally significant. The skeletal muscle macroarray platform was first used to search for evidence of tissue-specific expression, focusing on the biological function of genes/transcripts, since gene expression profiles generated across tissues were found to be reliable and consistent. Hierarchical clustering analysis revealed consistent clustering among genes assigned to 'developmental growth', such as the ontology genes and germ layers. Accuracy of the expression data was supported by comparing information from known transcripts and tissue from which the transcript was derived with macroarray data. Hybridization assays resulted in consistent tissue expression profile, which will be useful to dissect tissue-regulatory networks and to predict functions of novel genes identified after extensive sequencing of the genomes of model organisms. Screening our skeletal-muscle platform using 5 chicken adult tissues allowed us identifying 43 'tissue-specific' transcripts, and 112 co-expressed uncharacterized transcripts with 62 putative motifs. This platform also represents an important tool for functional investigation of novel genes; to determine expression pattern according to developmental stages; to evaluate differences in muscular growth potential between chicken lines, and to identify tissue-specific genes.