Inference of differentially expressed genes using generalized linear mixed models in a pairwise fashion.

Background: Technological advances involving RNA-Seq and Bioinformatics allow quantifying the transcriptional levels of genes in cells, tissues, and cell lines, permitting the identification of Differentially Expressed Genes (DEGs). DESeq2 and edgeR are well-established computational tools used for this purpose and they are based upon generalized linear models (GLMs) that consider only fixed effects in modeling. However, the inclusion of random effects reduces the risk of missing potential DEGs that may be essential in the context of the biological phenomenon under investigation. The generalized linear mixed models (GLMM) can be used to include both effects. Methods: We present DEGRE (Differentially Expressed Genes with Random Effects), a user-friendly tool capable of inferring DEGs where fixed and random effects on individuals are considered in the experimental design of RNA-Seq research. DEGRE preprocesses the raw matrices before fitting GLMMs on the genes and the derived regression coefficients are analyzed using the Wald statistical test. DEGRE offers the Benjamini-Hochberg or Bonferroni techniques for P-value adjustment. Results: The datasets used for DEGRE assessment were simulated with known identification of DEGs. These have fixed effects, and the random effects were estimated and inserted to measure the impact of experimental designs with high biological variability. For DEGs' inference, preprocessing effectively prepares the data and retains overdispersed genes. The biological coefficient of variation is inferred from the counting matrices to assess variability before and after the preprocessing. The DEGRE is computationally validated through its performance by the simulation of counting matrices, which have biological variability related to fixed and random effects. DEGRE also provides improved assessment measures for detecting DEGs in cases with higher biological variability. We show that the preprocessing established here effectively removes technical variation from those matrices. This tool also detects new potential candidate DEGs in the transcriptome data of patients with bipolar disorder, presenting a promising tool to detect more relevant genes. Conclusions: DEGRE provides data preprocessing and applies GLMMs for DEGs' inference. The preprocessing allows efficient remotion of genes that could impact the inference. Also, the computational and biological validation of DEGRE has shown to be promising in identifying possible DEGs in experiments derived from complex experimental designs. This tool may help handle random effects on individuals in the inference of DEGs and presents a potential for discovering new interesting DEGs for further biological investigation.

Zika virus disrupts gene expression in human myoblasts and myotubes: Relationship with susceptibility to infection.

The tropism of Zika virus (ZIKV) has been described in the nervous system, blood, placenta, thymus, and skeletal muscle. We investigated the mechanisms of skeletal muscle susceptibility to ZIKV using an in vitro model of human skeletal muscle myogenesis, in which myoblasts differentiate into myotubes. Myoblasts were permissive to ZIKV infection, generating productive viral particles, while myotubes controlled ZIKV replication. To investigate the underlying mechanisms, we used gene expression profiling. First, we assessed gene changes in myotubes compared with myoblasts in the model without infection. As expected, we observed an increase in genes and pathways related to the contractile muscle system in the myotubes, a reduction in processes linked to proliferation, migration and cytokine production, among others, confirming the myogenic capacity of our system in vitro. A comparison between non-infected and infected myoblasts revealed more than 500 differentially expressed genes (DEGs). In contrast, infected myotubes showed almost 2,000 DEGs, among which we detected genes and pathways highly or exclusively expressed in myotubes, including those related to antiviral and innate immune responses. Such gene modulation could explain our findings showing that ZIKV also invades myotubes but does not replicate in these differentiated cells. In conclusion, we showed that ZIKV largely (but differentially) disrupts gene expression in human myoblasts and myotubes. Identifying genes involved in myotube resistance can shed light on potential antiviral mechanisms against ZIKV infection.

The transposable element-rich genome of the cereal pest Sitophilus oryzae.

BACKGROUND: The rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions. RESULTS: We sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino acids required for insect development and cuticle biosynthesis. CONCLUSIONS: Here we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes.

Testing the genomic stability of the Brazilian yellow fever vaccine strain using next-generation sequencing data.

The live attenuated yellow fever (YF) vaccine was developed in the 1930s. Currently, the 17D and 17DD attenuated substrains are used for vaccine production. The 17D strain is used for vaccine production by several countries, while the 17DD strain is used exclusively in Brazil. The cell passages carried out through the seed-lot system of vaccine production influence the presence of quasispecies causing changes in the stability and immunogenicity of attenuated genotypes by increasing attenuation or virulence. Using next-generation sequencing, we carried out genomic characterization and genetic diversity analysis between vaccine lots of the Brazilian YF vaccine, produced by BioManguinhos-Fiocruz, and used during 11 years of vaccination in Brazil. We present 20 assembled and annotated genomes from the Brazilian 17DD vaccine strain, eight single nucleotide polymorphisms and the quasispecies spectrum reconstruction for the 17DD vaccine, through a pipeline here introduced. The V2IDA pipeline provided a relationship between low genetic diversity, maintained through the seed lot system, and the confirmation of genetic stability of lots of the Brazilian vaccine against YF. Our study sets precedents for use of V2IDA in genetic diversity analysis and in silico stability investigation of attenuated viral vaccines, facilitating genetic surveillance during the vaccine production process.

Complete Genome Sequence of Strain SS-5, a Magnetotactic Gammaproteobacterium Isolated from the Salton Sea, a Shallow, Saline, Endorheic Rift Lake Located on the San Andreas Fault in California.

We report the 3.7-Mb genome sequence of strain SS-5, a magnetotactic, sulfur-oxidizing rod and member of the family Chromatiaceae of the class Gammaproteobacteria, which biomineralizes membrane-bound, elongated, prismatic octahedral, magnetite nanocrystals. This genome sequence brings further diversity for understanding the origin and evolution of magnetotaxis and magnetosome biomineralization.

Structural analysis of a novel N-carbamoyl-d-amino acid amidohydrolase from a Brazilian Bradyrhizobium japonicum strain: In silico insights by molecular modelling, docking and molecular dynamics.

In this work we performed several in silico analyses to describe the relevant structural aspects of an enzyme N-Carbamoyl-d-amino acid amidohydrolase (d-NCAase) encoded on the genome of the Brazilian strain CPAC 15 (=SEMIA 5079) of Bradyrhizobium japonicum, a nonpathogenic species belonging to the order Rhizobiales. d-NCAase has wide applications particularly in the pharmaceutical industry, since it catalyzes the production of d-amino acids such as D-p-hydroxyphenylglycine (D-HPG), an intermediate in the synthesis of ß-lactam antibiotics. We applied a homology modelling approach and 50 ns of molecular dynamics simulations to predict the structure and the intersubunit interactions of this novel d-NCAase. Also, in order to evaluate the substrate binding site, the model was subjected to 50 ns of molecular dynamics simulations in the presence of N-Carbamoyl-d-p-hydroxyphenylglycine (Cp-HPG) (a d-NCAase canonical substrate) and water-protein/water-substrate interactions analyses were performed. Overall, the structural analysis and the molecular dynamics simulations suggest that d-NCAase of B. japonicum CPAC-15 has a homodimeric structure in solution. Here, we also examined the substrate specificity of the catalytic site of our model and the interactions with water molecules into the active binding site were comprehensively discussed. Also, these simulations showed that the amino acids Lys123, His125, Pro127, Cys172, Asp174 and Arg176 are responsible for recognition of ligand in the active binding site through several chemical associations, such as hydrogen bonds and hydrophobic interactions. Our results show a favourable environment for a reaction of hydrolysis that transforms N-Carbamoyl-d-p-hydroxyphenylglycine (Cp-HPG) into the active compound D-p-hydroxyphenylglycine (D-HPG). This work envisage the use of d-NCAase from the Brazilian Bradyrhizobium japonicum strain CPAC-15 (=SEMIA 5079) for the industrial production of D-HPG, an important intermediate for semi-synthesis of ß-lactam antibiotics such as penicillins, cephalosporins and amoxicillin.

Combined genomic and structural analyses of a cultured magnetotactic bacterium reveals its niche adaptation to a dynamic environment.

BACKGROUND: Magnetotactic bacteria (MTB) are a unique group of prokaryotes that have a potentially high impact on global geochemical cycling of significant primary elements because of their metabolic plasticity and the ability to biomineralize iron-rich magnetic particles called magnetosomes. Understanding the genetic composition of the few cultivated MTB along with the unique morphological features of this group of bacteria may provide an important framework for discerning their potential biogeochemical roles in natural environments. RESULTS: Genomic and ultrastructural analyses were combined to characterize the cultivated magnetotactic coccus Magnetofaba australis strain IT-1. Cells of this species synthesize a single chain of elongated, cuboctahedral magnetite (Fe3O4) magnetosomes that cause them to align along magnetic field lines while they swim being propelled by two bundles of flagella at velocities up to 300 µm s-1. High-speed microscopy imaging showed the cells move in a straight line rather than in the helical trajectory described for other magnetotactic cocci. Specific genes within the genome of Mf. australis strain IT-1 suggest the strain is capable of nitrogen fixation, sulfur reduction and oxidation, synthesis of intracellular polyphosphate granules and transporting iron with low and high affinity. Mf. australis strain IT-1 and Magnetococcus marinus strain MC-1 are closely related phylogenetically although similarity values between their homologous proteins are not very high. CONCLUSION: Mf. australis strain IT-1 inhabits a constantly changing environment and its complete genome sequence reveals a great metabolic plasticity to deal with these changes. Aside from its chemoautotrophic and chemoheterotrophic metabolism, genomic data indicate the cells are capable of nitrogen fixation, possess high and low affinity iron transporters, and might be capable of reducing and oxidizing a number of sulfur compounds. The relatively large number of genes encoding transporters as well as chemotaxis receptors in the genome of Mf. australis strain IT-1 combined with its rapid swimming velocities, indicate that cells respond rapidly to environmental changes.

Adaptive processes of Staphylococcus aureus isolates during the progression from acute to chronic bone and joint infections in patients.

Staphylococcus aureus bone and joint infection (BJI) is associated with significant rates of chronicity and relapse. In this study, we investigated how S. aureus is able to adapt to the human environment by comparing isolates from single patients with persisting or relapsing BJIs that were recovered during the initial and recurrent BJI episodes. In vitro and in vivo assays and whole-genome sequencing analyses revealed that the recurrent isolates induced a reduced inflammatory response, formed more biofilms, persisted longer in the intracellular compartments of host bone cells, were less cytotoxic and induced less mortality in a mouse infection model compared with the initial isolates despite the lack of significant changes at the genomic level. These findings suggest that S. aureus BJI chronicization is associated with an in vivo bacterial phenotypical adaptation that leads to decreased virulence and host immune escape, which is linked to increased intraosteoblastic persistence and biofilm formation.

Complete genome sequence of Burkholderia phenoliruptrix BR3459a (CLA1), a heat-tolerant, nitrogen-fixing symbiont of Mimosa flocculosa.

The genus Burkholderia represents a challenge to the fields of taxonomy and phylogeny and, especially, to the understanding of the contrasting roles as either opportunistic pathogens or bacteria with biotechnological potential. Few genomes of nonpathogenic strains, especially of diazotrophic symbiotic bacteria, have been sequenced to improve understanding of the genus. Here, we contribute with the complete genome sequence of Burkholderia phenoliruptrix strain BR3459a (CLA1), an effective diazotrophic symbiont of the leguminous tree Mimosa flocculosa Burkart, which is endemic to South America.

Exploration of the core metabolism of symbiotic bacteria.

BACKGROUND: A large number of genome-scale metabolic networks is now available for many organisms, mostly bacteria. Previous works on minimal gene sets, when analysing host-dependent bacteria, found small common sets of metabolic genes. When such analyses are restricted to bacteria with similar lifestyles, larger portions of metabolism are expected to be shared and their composition is worth investigating. Here we report a comparative analysis of the small molecule metabolism of symbiotic bacteria, exploring common and variable portions as well as the contribution of different lifestyle groups to the reduction of a common set of metabolic capabilities. RESULTS: We found no reaction shared by all the bacteria analysed. Disregarding those with the smallest genomes, we still do not find a reaction core, however we did find a core of biochemical capabilities. While obligate intracellular symbionts have no core of reactions within their group, extracellular and cell-associated symbionts do have a small core composed of disconnected fragments. In agreement with previous findings in Escherichia coli, their cores are enriched in biosynthetic processes whereas the variable metabolisms have similar ratios of biosynthetic and degradation reactions. Conversely, the variable metabolism of obligate intracellular symbionts is enriched in anabolism. CONCLUSION: Even when removing the symbionts with the most reduced genomes, there is no core of reactions common to the analysed symbiotic bacteria. The main reason is the very high specialisation of obligate intracellular symbionts, however, host-dependence alone is not an explanation for such absence. The composition of the metabolism of cell-associated and extracellular bacteria shows that while they have similar needs in terms of the building blocks of their cells, they have to adapt to very distinct environments. On the other hand, in obligate intracellular bacteria, catabolism has largely disappeared, whereas synthetic routes appear to have been selected for depending on the nature of the symbiosis. As more genomes are added, we expect, based on our simulations, that the core of cell-associated and extracellular bacteria continues to diminish, converging to approximately 60 reactions.