Guilherme Loss de Morais

Evolutionary view of acyl-CoA diacylglycerol acyltransferase (DGAT), a key enzyme in neutral lipid biosynthesis

BackgroundTriacylglycerides (TAGs) are a class of neutral lipids that represent the most important storage form of energy for eukaryotic cells. DGAT (acyl-CoA: diacylglycerol acyltransferase; EC 2.3.1.20) is a transmembrane enzyme that acts in the final and committed step of TAG synthesis, and it has been proposed to be the rate-limiting enzyme in plant storage lipid accumulation. In fact, two different enzymes identified in several eukaryotic species, DGAT1 and DGAT2, are the main enzymes responsible for TAG synthesis. These enzymes do not share high DNA or protein sequence similarities, and it has been suggested that they play non-redundant roles in different tissues and in some species in TAG synthesis. Despite a number of previous studies on the DGAT1 and DGAT2 genes, which have emphasized their importance as potential obesity treatment targets to increase triacylglycerol accumulation, little is known about their evolutionary timeline in eukaryotes. The goal of this study was to examine the evolutionary relationship of the DGAT1 and DGAT2 genes across eukaryotic organisms in order to infer their origin.ResultsWe have conducted a broad survey of fully sequenced genomes, including representatives of Amoebozoa, yeasts, fungi, algae, musses, plants, vertebrate and invertebrate species, for the presence of DGAT1 and DGAT2 gene homologs. We found that the DGAT1 and DGAT2 genes are nearly ubiquitous in eukaryotes and are readily identifiable in all the major eukaryotic groups and genomes examined. Phylogenetic analyses of the DGAT1 and DGAT2 amino acid sequences revealed evolutionary partitioning of the DGAT protein family into two major DGAT1 and DGAT2 clades. Protein secondary structure and hydrophobic-transmembrane analysis also showed differences between these enzymes. The analysis also revealed that the MGAT2 and AWAT genes may have arisen from DGAT2 duplication events.ConclusionsIn this study, we identified several DGAT1 and DGAT2 homologs in eukaryote taxa. Overall, the data show that DGAT1 and DGAT2 are present in most eukaryotic organisms and belong to two different gene families. The phylogenetic and evolutionary analyses revealed that DGAT1 and DGAT2 evolved separately, with functional convergence, despite their wide molecular and structural divergence.

Chromosomal copy number variation reveals differential levels of genomic plasticity in distinct Trypanosoma cruzi strains

BackgroundTrypanosoma cruzi, the etiologic agent of Chagas disease, is currently divided into six discrete typing units (DTUs), named TcI–TcVI. CL Brener, the reference strain of the T. cruzi genome project, is a hybrid with a genome assembled into 41 putative chromosomes. Gene copy number variation (CNV) is well documented as an important mechanism to enhance gene expression and variability in T. cruzi. Chromosomal CNV (CCNV) is another level of gene CNV in which whole blocks of genes are expanded simultaneously. Although the T. cruzi karyotype is not well defined, several studies have demonstrated a significant variation in the size and content of chromosomes between different T. cruzi strains. Despite these studies, the extent of diversity in CCNV among T. cruzi strains based on a read depth coverage analysis has not been determined.ResultsWe identify the CCNV in T. cruzi strains from the TcI, TcII and TcIII DTUs, by analyzing the depth coverage of short reads from these strains using the 41 CL Brener chromosomes as reference. This study led to the identification of a broader extent of CCNV in T. cruzi than was previously speculated. The TcI DTU strains have very few aneuploidies, while the strains from TcII and TcIII DTUs present a high degree of chromosomal expansions. Chromosome 31, which is the only chromosome that is supernumerary in all six T. cruzi samples evaluated in this study, is enriched with genes related to glycosylation pathways, highlighting the importance of glycosylation to parasite survival.ConclusionsIncreased gene copy number due to chromosome amplification may contribute to alterations in gene expression, which represents a strategy that may be crucial for parasites that mainly depend on post-transcriptional mechanisms to control gene expression.

Metatranscriptomic analysis of small RNAs present in soybean deep sequencing libraries

A large number of small RNAs unrelated to the soybean genome were identified after deep sequencing of soybean small RNA libraries. A metatranscriptomic analysis was carried out to identify the origin of these sequences. Comparative analyses of small interference RNAs (siRNAs) present in samples collected in open areas corresponding to soybean field plantations and samples from soybean cultivated in greenhouses under a controlled environment were made. Different pathogenic, symbiotic and free-living organisms were identified from samples of both growth systems. They included viruses, bacteria and different groups of fungi. This approach can be useful not only to identify potentially unknown pathogens and pests, but also to understand the relations that soybean plants establish with microorganisms that may affect, directly or indirectly, plant health and crop production.

Unusual RNA plant virus integration in the soybean genome leads to the production of small RNAs.

Horizontal gene transfer (HGT) is known to be a major force in genome evolution. The acquisition of genes from viruses by eukaryotic genomes is a well-studied example of HGT, including rare cases of non-retroviral RNA virus integration. The present study describes the integration of cucumber mosaic virus RNA-1 into soybean genome. After an initial metatranscriptomic analysis of small RNAs derived from soybean, the de novo assembly resulted a 3029-nt contig homologous to RNA-1. The integration of this sequence in the soybean genome was confirmed by DNA deep sequencing. The locus where the integration occurred harbors the full RNA-1 sequence followed by the partial sequence of an endogenous mRNA and another sequence of RNA-1 as an inverted repeat and allowing the formation of a hairpin structure. This region recombined into a retrotransposon located inside an exon of a soybean gene. The nucleotide similarity of the integrated sequence compared to other Cucumber mosaic virus sequences indicates that the integration event occurred recently. We described a rare event of non-retroviral RNA virus integration in soybean that leads to the production of a double-stranded RNA in a similar fashion to virus resistance RNAi plants.

Mycoplasma non-coding RNA: identification of small RNAs and targets

BackgroundBacterial non-coding RNAs act by base-pairing as regulatory elements in crucial biological processes. We performed the identification of trans-encoded small RNAs (sRNA) from the genomes of Mycoplama hyopneumoniae, Mycoplasma flocculare and Mycoplasma hyorhinis, which are Mycoplasma species that have been identified in the porcine respiratory system.ResultsA total of 47, 15 and 11 putative sRNAs were predicted in M. hyopneumoniae, M. flocculare and M. hyorhinis, respectively. A comparative genomic analysis revealed the presence of species or lineage specific sRNA candidates. Furthermore, the expression profile of some M. hyopneumoniae sRNAs was determined by a reverse transcription amplification approach, in three different culture conditions. All tested sRNAs were transcribed in at least one condition. A detailed investigation revealed a differential expression profile for two M. hyopneumoniae sRNAs in response to oxidative and heat shock stress conditions, suggesting that their expression is influenced by environmental signals. Moreover, we analyzed sRNA-mRNA hybrids and accessed putative target genes for the novel sRNA candidates. The majority of the sRNAs showed interaction with multiple target genes, some of which could be linked to pathogenesis and cell homeostasis activity.ConclusionThis study contributes to our knowledge of Mycoplasma sRNAs and their response to environmental changes. Furthermore, the mRNA target prediction provides a perspective for the characterization and comprehension of the function of the sRNA regulatory mechanisms.

Intraclonal Genome Stability of the Metallo-β-lactamase SPM-1-producing Pseudomonas aeruginosa ST277, an Endemic Clone Disseminated in Brazilian Hospitals

Carbapenems represent the mainstay therapy for the treatment of serious P. aeruginosa infections. However, the emergence of carbapenem resistance has jeopardized the clinical use of this important class of compounds. The production of SPM-1 metallo-β-lactamase has been the most common mechanism of carbapenem resistance identified in P. aeruginosa isolated from Brazilian medical centers. Interestingly, a single SPM-1-producing P. aeruginosa clone belonging to the ST277 has been widely spread within the Brazilian territory. In the current study, we performed a next-generation sequencing of six SPM-1-producing P. aeruginosa ST277 isolates. The core genome contains 5899 coding genes relative to the reference strain P. aeruginosa PAO1. A total of 26 genomic islands were detected in these isolates. We identified remarkable elements inside these genomic islands, such as copies of the blaSPM−1 gene conferring resistance to carbapenems and a type I-C CRISPR-Cas system, which is involved in protection of the chromosome against foreign DNA. In addition, we identified single nucleotide polymorphisms causing amino acid changes in antimicrobial resistance and virulence-related genes. Together, these factors could contribute to the marked resistance and persistence of the SPM-1-producing P. aeruginosa ST277 clone. A comparison of the SPM-1-producing P. aeruginosa ST277 genomes showed that their core genome has a high level nucleotide similarity and synteny conservation. The variability observed was mainly due to acquisition of genomic islands carrying several antibiotic resistance genes.

Genome-wide diversity and differentiation in New World populations of the human malaria parasite Plasmodium vivax

Background The Americas were the last continent colonized by humans carrying malaria parasites. Plasmodium falciparum from the New World shows very little genetic diversity and greater linkage disequilibrium, compared with its African counterparts, and is clearly subdivided into local, highly divergent populations. However, limited available data have revealed extensive genetic diversity in American populations of another major human malaria parasite, P. vivax. Methods We used an improved sample preparation strategy and next-generation sequencing to characterize 9 high-quality P. vivax genome sequences from northwestern Brazil. These new data were compared with publicly available sequences from recently sampled clinical P. vivax isolates from Brazil (BRA, total n = 11 sequences), Peru (PER, n = 23), Colombia (COL, n = 31), and Mexico (MEX, n = 19). Principal findings/Conclusions We found that New World populations of P. vivax are as diverse (nucleotide diversity π between 5.2 × 10−4 and 6.2 × 10−4) as P. vivax populations from Southeast Asia, where malaria transmission is substantially more intense. They display several non-synonymous nucleotide substitutions (some of them previously undescribed) in genes known or suspected to be involved in antimalarial drug resistance, such as dhfr, dhps, mdr1, mrp1, and mrp-2, but not in the chloroquine resistance transporter ortholog (crt-o) gene. Moreover, P. vivax in the Americas is much less geographically substructured than local P. falciparum populations, with relatively little between-population genome-wide differentiation (pairwise FST values ranging between 0.025 and 0.092). Finally, P. vivax populations show a rapid decline in linkage disequilibrium with increasing distance between pairs of polymorphic sites, consistent with very frequent outcrossing. We hypothesize that the high diversity of present-day P. vivax lineages in the Americas originated from successive migratory waves and subsequent admixture between parasite lineages from geographically diverse sites. Further genome-wide analyses are required to test the demographic scenario suggested by our data.

Complete genome sequence of the MRSA isolate HC1335 from ST239 lineage displaying a truncated AgrC histidine kinase receptor

Methicillin-resistant Staphylococcus aureus (MRSA) is still one of the most important hospital pathogen globally. The multiresistant isolates of the ST239-SCCmecIII lineage are spread over large geographic regions, colonizing and infecting hospital patients in virtually all continents. The balance between fitness (adaptability) and virulence potential is likely to represent an important issue in the clonal shift dynamics leading the success of some specific MRSA clones over another. The accessory gene regulator (agr) is the master quorum sensing system of staphylococci playing a role in the global regulation of key virulence factors. Consequently, agr inactivation in S. aureus may represent a significant mechanism of genetic variability in the adaptation of this healthcare-associated pathogen. We report here the complete genome sequence of the methicillin-resistant S. aureus, isolate HC1335, a variant of the ST239 lineage, which presents a natural insertion of an IS256 transposase element in the agrC gene encoding AgrC histidine kinase receptor.

RFMirTarget: A Random Forest Classifier for Human miRNA Target Gene Prediction

MicroRNAs (miRNAs) are key regulators of eukaryotic gene expression whose fundamental role has been already identified in many cell pathways. The correct identification of miRNAs targets is a major challenge in bioinformatics. So far, machine learning-based methods for miRNA-target prediction have shown the best results in terms of specificity and sensitivity. However, despite its well-known efficiency in other classifying tasks, the random forest algorithm has not been employed in this problem. Therefore, in this work we present RFMirTarget, an efficient random forest miRNA-target prediction system. Our tool analyzes the alignment between a candidate miRNA-target pair and extracts a set of structural, thermodynamics, alignment and position-based features. Experiments have shown that RFMirTarget achieves a Matthew’s correlation coefficient nearly 48% greater than the performance reported for the MultiMiTar, which was trained upon the same data set. In addition, tests performed with RFMirTarget reinforce the importance of the seed region for target prediction accuracy.

Comparative Genomic Analysis of a Clinical Isolate of Klebsiella quasipneumoniae subsp. similipneumoniae, a KPC-2 and OKP-B-6 Beta-Lactamases Producer Harboring Two Drug-Resistance Plasmids from Southeast Brazil

The aim of this study was to unravel the genetic determinants responsible for multidrug (including carbapenems) resistance and virulence in a clinical isolate of Klebsiella quasipneumoniae subsp. similipneumoniae by whole-genome sequencing and comparative analyses. Eighty-three clinical isolates initially identified as carbapenem-resistant K. pneumoniae were collected from nosocomial infections in southeast Brazil. After RAPD screening, the KPC-142 isolate, showing the most divergent DNA pattern, was selected for complete genome sequencing in an Illumina HiSeq 2500 instrument. Reads were assembled into scaffolds, gaps between scaffolds were resolved by in silico gap filling and extensive bioinformatics analyses were performed, using multiple comparative analysis tools and databases. Genome sequencing allowed to correct the classification of the KPC-142 isolate as K. quasipneumoniae subsp. similipneumoniae. To the best of our knowledge this is the first complete genome reported to date of a clinical isolate of this subspecies harboring both class A beta-lactamases KPC-2 and OKP-B-6 from South America. KPC-142 has one 5.2 Mbp chromosome (57.8% G+C) and two plasmids: 190 Kbp pKQPS142a (50.7% G+C) and 11 Kbp pKQPS142b (57.3% G+C). The 3 Kbp region in pKQPS142b containing the blaKPC−2 was found highly similar to that of pKp13d of K. pneumoniae Kp13 isolated in Southern Brazil in 2009, suggesting the horizontal transfer of this resistance gene between different species of Klebsiella. KPC-142 additionally harbors an integrative conjugative element ICEPm1 that could be involved in the mobilization of pKQPS142b and determinants of resistance to other classes of antimicrobials, including aminoglycoside and silver. We present the completely assembled genome sequence of a clinical isolate of K. quasipneumoniae subsp. similipneumoniae, a KPC-2 and OKP-B-6 beta-lactamases producer and discuss the most relevant genomic features of this important resistant pathogen in comparison to several strains belonging to K. quasipneumoniae subsp. similipneumoniae (phylogroup II-B), K. quasipneumoniae subsp. quasipneumoniae (phylogroup II-A), K. pneumoniae (phylogroup I), and K. variicola (phylogroup III). Our study contributes to the description of the characteristics of a novel K. quasipneumoniae subsp. similipneumoniae strain circulating in South America that currently represent a serious potential risk for nosocomial settings.