Discovery of antimicrobial peptides in the global microbiome with machine learning.

Novel antibiotics are urgently needed to combat the antibiotic-resistance crisis. We present a machine-learning-based approach to predict antimicrobial peptides (AMPs) within the global microbiome and leverage a vast dataset of 63,410 metagenomes and 87,920 prokaryotic genomes from environmental and host-associated habitats to create the AMPSphere, a comprehensive catalog comprising 863,498 non-redundant peptides, few of which match existing databases. AMPSphere provides insights into the evolutionary origins of peptides, including by duplication or gene truncation of longer sequences, and we observed that AMP production varies by habitat. To validate our predictions, we synthesized and tested 100 AMPs against clinically relevant drug-resistant pathogens and human gut commensals both in vitro and in vivo. A total of 79 peptides were active, with 63 targeting pathogens. These active AMPs exhibited antibacterial activity by disrupting bacterial membranes. In conclusion, our approach identified nearly one million prokaryotic AMP sequences, an open-access resource for antibiotic discovery.

Challenges in computational discovery of bioactive peptides in 'omics data.

Peptides have a plethora of activities in biological systems that can potentially be exploited biotechnologically. Several peptides are used clinically, as well as in industry and agriculture. The increase in available 'omics data has recently provided a large opportunity for mining novel enzymes, biosynthetic gene clusters, and molecules. While these data primarily consist of DNA sequences, other types of data provide important complementary information. Due to their size, the approaches proven successful at discovering novel proteins of canonical size cannot be naïvely applied to the discovery of peptides. Peptides can be encoded directly in the genome as short open reading frames (smORFs), or they can be derived from larger proteins by proteolysis. Both of these peptide classes pose challenges as simple methods for their prediction result in large numbers of false positives. Similarly, functional annotation of larger proteins, traditionally based on sequence similarity to infer orthology and then transferring functions between characterized proteins and uncharacterized ones, cannot be applied for short sequences. The use of these techniques is much more limited and alternative approaches based on machine learning are used instead. Here, we review the limitations of traditional methods as well as the alternative methods that have recently been developed for discovering novel bioactive peptides with a focus on prokaryotic genomes and metagenomes.

Seasonality of biogeochemically relevant microbial genes in a coastal ocean microbiome.

Microbes drive the biogeochemical cycles of marine ecosystems through their vast metabolic diversity. While we have a fairly good understanding of the spatial distribution of these metabolic processes in various ecosystems, less is known about their seasonal dynamics. We investigated the annual patterns of 21 biogeochemical relevant functions in an oligotrophic coastal ocean site by analysing the presence of key genes, analysing high-rank gene taxonomy and the dynamics of nucleotide variants. Most genes presented seasonality: photoheterotrophic processes were enriched during spring, phosphorous-related genes were dominant during summer, coinciding with potential phosphate limitation, and assimilatory nitrate reductases appeared mostly during summer and autumn, correlating negatively with nitrate availability. Additionally, we identified the main taxa driving each function at each season and described the role of underrecognized taxa such as Litoricolaceae in carbon fixation (rbcL), urea degradation (ureC), and CO oxidation (coxL). Finally, the seasonality of single variants of some families presented a decoupling between the taxonomic abundance patterns and the functional gene patterns, implying functional specialization of the different genera. Our study unveils the seasonality of key biogeochemical functions and the main taxonomic groups that harbour these relevant functions in a coastal ocean ecosystem.

Microbial population genomes from the Amazon River reveal possible modulation of the organic matter degradation process in tropical freshwaters.

Rivers connect the carbon cycle in land with that in aquatic ecosystems by transporting and transforming terrestrial organic matter (TeOM). The Amazon River receives huge loads of TeOM from the surrounding rainforest, promoting a substantial microbial heterotrophic activity and consequently, CO2 outgassing. In the Amazon River, microbes degrade up to 55% of the lignin present in the TeOM. Yet, the main microbial genomes involved in TeOM degradation were unknown. Here, we characterize 51 population genomes (PGs) representing some of the most abundant microbes in the Amazon River deriving from 106 metagenomes. The 51 reconstructed PGs are among the most abundant microbes in the Amazon River, and 53% of them are not able to degrade TeOM. Among the PGs capable of degrading TeOM, 20% were exclusively cellulolytic, while the others could also oxidize lignin. The transport and consumption of lignin oxidation byproducts seemed to be decoupled from the oxidation process, being apparently performed by different groups of microorganisms. By connecting the genomic features of abundant microbes in the Amazon River with the degradation machinery of TeOM, we suggest that a complex microbial consortium could explain the quick turnover of TeOM previously observed in this ecosystem.

Epigenetic modifications and their relation to caste and sex determination and adult division of labor in the stingless bee Melipona scutellaris.

Stingless bees of the genus Melipona, have long been considered an enigmatic case among social insects for their mode of caste determination, where in addition to larval food type and quantity, the genotype also has a saying, as proposed over 50 years ago by Warwick E. Kerr. Several attempts have since tried to test his Mendelian two-loci/two-alleles segregation hypothesis, but only recently a single gene crucial for sex determination in bees was evidenced to be sex-specifically spliced and also caste-specifically expressed in a Melipona species. Since alternative splicing is frequently associated with epigenetic marks, and the epigenetic status plays a major role in setting the caste phenotype in the honey bee, we investigated here epigenetic chromatin modification in the stingless bee Melipona scutellaris. We used an ELISA-based methodology to quantify global methylation status and western blot assays to reveal histone modifications. The results evidenced DNA methylation/demethylation events in larvae and pupae, and significant differences in histone methylation and phosphorylation between newly emerged adult queens and workers. The epigenetic dynamics seen in this stingless bee species represent a new facet in the caste determination process in Melipona bees and suggest a possible mechanism that is likely to link a genotype component to the larval diet and adult social behavior of these bees.

A catalog of small proteins from the global microbiome.

Small open reading frames (smORFs) shorter than 100 codons are widespread and perform essential roles in microorganisms, where they encode proteins active in several cell functions, including signal pathways, stress response, and antibacterial activities. However, the ecology, distribution and role of small proteins in the global microbiome remain unknown. Here, we construct a global microbial smORFs catalog (GMSC) derived from 63,410 publicly available metagenomes across 75 distinct habitats and 87,920 high-quality isolate genomes. GMSC contains 965 million non-redundant smORFs with comprehensive annotations. We find that archaea harbor more smORFs proportionally than bacteria. We moreover provide a tool called GMSC-mapper to identify and annotate small proteins from microbial (meta)genomes. Overall, this publicly-available resource demonstrates the immense and underexplored diversity of small proteins.

Computational exploration of the global microbiome for antibiotic discovery.

Novel antibiotics are urgently needed to combat the antibiotic-resistance crisis. We present a machine learning-based approach to predict prokaryotic antimicrobial peptides (AMPs) by leveraging a vast dataset of 63,410 metagenomes and 87,920 microbial genomes. This led to the creation of AMPSphere, a comprehensive catalog comprising 863,498 non-redundant peptides, the majority of which were previously unknown. We observed that AMP production varies by habitat, with animal-associated samples displaying the highest proportion of AMPs compared to other habitats. Furthermore, within different human-associated microbiota, strain-level differences were evident. To validate our predictions, we synthesized and experimentally tested 50 AMPs, demonstrating their efficacy against clinically relevant drug-resistant pathogens both in vitro and in vivo. These AMPs exhibited antibacterial activity by targeting the bacterial membrane. Additionally, AMPSphere provides valuable insights into the evolutionary origins of peptides. In conclusion, our approach identified AMP sequences within prokaryotic microbiomes, opening up new avenues for the discovery of antibiotics.

Sugarcane cystatins: From discovery to biotechnological applications.

Phytocystatins are tight-binding cysteine protease inhibitors produced by plants. The first phytocystatin described was isolated from Oryza sativa and, since then, cystatins from several plant species were reported, including from sugarcane. Sugarcane cystatins were unraveled in Sugarcane EST project database, after sequencing of cDNA libraries from various sugarcane tissues at different developmental stages and six sugarcane cystatins were cloned, expressed and characterized (CaneCPI-1 to CaneCPI-6). These recombinant proteins were produced in different expression systems and inhibited several cysteine proteases, including human cathepsins B and L, which can be involved in pathologies, such as cancer. In this review, we summarize a comprehensive history of all sugarcane cystatins, presenting an updated phylogenetic analysis; chromosomal localization, and genomic organization. We also present protein docking of CaneCPI-5 in the active site of human cathepsin B, insights about canecystatins structures; recombinant expression in different systems, comparison of their inhibitory activities against human cysteine cathepsins B, K, L, S, V, falcipains from Plasmodium falciparum and a cathepsin L-like from the sugarcane weevil Sphenophorus levis; and enlighten their potential and current applications in agriculture and health.

Isohemigossypolone: Antiophidic properties of a naphthoquinone isolated from Pachira aquatica Aubl.

We investigated the antiophidic properties of isohemigossypolone (ISO), a naphthoquinone isolated from the outer bark of the Pachira aquatic Aubl. The inhibition of phospholipase A2, coagulant, fibrinogenolytic, hemorrhagic and myotoxic activities induced by Bothrops pauloensis venom (Pb) was investigated. For this, we use samples resulting from the incubation of Pb with ISO in different concentrations (1:1, 1:5 and 1:10 w/w), we also evaluated the condition of treatment using ISO after 15 min of venom inoculation. The activities of phospholipase A2, coagulant, fibrinogenolytic, hemorrhagic and myotoxic induced by the B. pauloensis venom were significantly inhibited when the ISO was pre-incubated with the crude venom. For in vivo neutralization tests, the results were observed even when the ISO was applied after 15 min of inoculation of the venom or metalloprotease (BthMP). Also, to identify the inhibition mechanism, we performed in silico assays, across simulations of molecular coupling and molecular dynamics, it was possible to identify the modes of interaction between ISO and bothropic toxins BmooMPα-I, Jararacussin-I and BNSP-7. The present study shows that naphthoquinone isohemigossypolone isolated from the P. aquatica plant inhibited part of the local and systemic damage caused by venom proteins, demonstrating the pharmacological potential of this compound in neutralizing the harmful effects caused by snakebites.

Macrel: antimicrobial peptide screening in genomes and metagenomes.

MOTIVATION: Antimicrobial peptides (AMPs) have the potential to tackle multidrug-resistant pathogens in both clinical and non-clinical contexts. The recent growth in the availability of genomes and metagenomes provides an opportunity for in silico prediction of novel AMP molecules. However, due to the small size of these peptides, standard gene prospection methods cannot be applied in this domain and alternative approaches are necessary. In particular, standard gene prediction methods have low precision for short peptides, and functional classification by homology results in low recall. RESULTS: Here, we present Macrel (for metagenomic AMP classification and retrieval), which is an end-to-end pipeline for the prospection of high-quality AMP candidates from (meta)genomes. For this, we introduce a novel set of 22 peptide features. These were used to build classifiers which perform similarly to the state-of-the-art in the prediction of both antimicrobial and hemolytic activity of peptides, but with enhanced precision (using standard benchmarks as well as a stricter testing regime). We demonstrate that Macrel recovers high-quality AMP candidates using realistic simulations and real data. AVAILABILITY: Macrel is implemented in Python 3. It is available as open source at https://github.com/BigDataBiology/macrel and through bioconda. Classification of peptides or prediction of AMPs in contigs can also be performed on the webserver: https://big-data-biology.org/software/macrel.

Uncovering the genomic potential of the Amazon River microbiome to degrade rainforest organic matter.

BACKGROUND: The Amazon River is one of the largest in the world and receives huge amounts of terrestrial organic matter (TeOM) from the surrounding rainforest. Despite this TeOM is typically recalcitrant (i.e. resistant to degradation), only a small fraction of it reaches the ocean, pointing to a substantial TeOM degradation by the river microbiome. Yet, microbial genes involved in TeOM degradation in the Amazon River were barely known. Here, we examined the Amazon River microbiome by analysing 106 metagenomes from 30 sampling points distributed along the river. RESULTS: We constructed the Amazon River basin Microbial non-redundant Gene Catalogue (AMnrGC) that includes ~ 3.7 million non-redundant genes, affiliating mostly to bacteria. We found that the Amazon River microbiome contains a substantial gene-novelty compared to other relevant known environments (rivers and rainforest soil). Genes encoding for proteins potentially involved in lignin degradation pathways were correlated to tripartite tricarboxylates transporters and hemicellulose degradation machinery, pointing to a possible priming effect. Based on this, we propose a model on how the degradation of recalcitrant TeOM could be modulated by labile compounds in the Amazon River waters. Our results also suggest changes of the microbial community and its genomic potential along the river course. CONCLUSIONS: Our work contributes to expand significantly our comprehension of the world's largest river microbiome and its potential metabolism related to TeOM degradation. Furthermore, the produced gene catalogue (AMnrGC) represents an important resource for future research in tropical rivers. Video abstract.

Recombinant expression, characterization and phylogenetic studies of novels cystatins-like proteins of sweet orange (Citrus sinensis) and clementine (Citrus clementina).

Phytocystatins are plant cystatins that are related to several physiological processes regulating endogenous cysteine proteases involved in seed development and germination, programmed cell death and response to stress conditions. In addition, phytocystatins can act in plant defense against exogenous peptidases from herbivorous insects, pathogens and nematodes. Considering that Citrus fruits are important to human nutrition and represent a high value crop in worldwide agriculture, in the present work, we performed the identification of putative cystatins from Citrus sinensis and from Citrus clementine and submitted them to phylogenetic analysis. Six cystatins from each species were identified as orthologous and classified into three well supported phylogenetic groups. Five cystatins representative of the phylogenetic groups were recombinantly expressed and the in vitro studies revealed them to be potent inhibitors against the cysteine peptidases papain, legumain, human cathepsins (B, L, S, K) and a cathepsin B-like from Diaphorina citri (the Asian Citrus psyllid). Our findings provide the C. clementina and C. sinensis cystatins classification and an enzyme-inhibitor interactions profile, which may reflect an evolutionary process of Citrus cystatins related to gene functions as initial germination rates and seedlings development as well associated to plant defense against pathogens, as insects and nematodes.

Flood Season Microbiota from the Amazon Basin Lakes: Analysis with Metagenome Sequencing.

Despite an apparent geographic separation of the Amazon water bodies, they are an interconnected system. During floods, the microbiota of rivers, lakes, and soil combines. This study used metagenomics sequencing to survey the microbiota of the Amazon Basin lakes during flood season, showing important patterns in microbial communities.

Metagenome Sequencing of Prokaryotic Microbiota Collected from Rivers in the Upper Amazon Basin.

Tropical freshwater environments, like rivers, are important reservoirs of microbial life. This study employed metagenomic sequencing to survey prokaryotic microbiota in the Solimões, Purus, and Urucu Rivers of the Amazon Basin in Brazil. We report a rich and diverse microbial community.

A transcriptomic survey of Migdolus fryanus (sugarcane rhizome borer) larvae.

Sugarcane, a major crop grown in the tropical and subtropical areas of the world, is produced mainly for sucrose, which is used as a sweetener or for the production of bioethanol. Among the numerous pests that significantly affect the yield of sugarcane, the sugarcane rhizome borer (Migdolus fryanus, a cerambycidae beetle) is known to cause severe damage to the crops in Brazil. The absence of molecular information about this insect reinforces the need for studies and an effective method to control this pest. In this study, RNA-Seq technology was employed to study different parts of M. fryanus larvae. The generated data will help in further investigations about the taxonomy, development, and adaptation of this insect. RNA was extracted from six different parts (head, fat body, integument, hindgut, midgut, and foregut) using Trizol methodology. Using Illumina paired-end sequencing technology and the Trinity platform, trimming and de novo assembly was performed, resulting in 44,567 contigs longer than 200 nt for a reunion of data from all transcriptomes, with a mean length of 1,095.27 nt. Transcripts were annotated using BLAST against different protein databanks (Uniprot/Swissprot, PFAM, KEEG, SignalP 4.1, Gene Ontology, and CAZY) and were compared for similarity using a Venn diagram. Differential expression patterns were studied for select genes through qPCR and FPKM comprising important protein families (digestive peptidases, glucosyl hydrolases, serine protease inhibitors and otopetrin), which allowed a better understanding of the insect's digestion, immunity and gravity sensorial mechanisms.

Metagenomics Analysis of Microorganisms in Freshwater Lakes of the Amazon Basin.

The Amazon Basin is the largest hydrographic basin on the planet, and the dynamics of its aquatic microorganisms strongly impact global biogeochemical cycles. However, it remains poorly studied. This metagenome project was performed to obtain a snapshot of prokaryotic microbiota from four important lakes in the Amazon Basin.

Tollip or not Tollip: what are the evolving questions behind it?

Tollip plays an important role in the interleukin-1 receptor IL-1R and Toll pathways. As a modulator of the immune pathway, it indirectly controls the amount of antimicrobial peptides. This could indicate a vital step in maintaining animal immune systems and preventing infection. Evolutionary questions are crucial to understanding the conservation and functioning of the biochemical pathways like the Tollip-mediated one. Through an analysis of 36 sequences of the Tollip protein from different animal taxa, downloaded from Kyoto Encyclopedia of Genes and Genomes (KEGG) databank, we inferred diverse evolutionary parameters, such as molecular selection and structure conservation, by analyzing residue by residue, beyond the canonical parameters to this type of study, as maximum likelihood trees. We found that Tollip presented different trends in its evolving history. In primates, the protein is becoming more unstable, just the opposite is observed in the arthropod group. The most interesting finding was the concentration of positively selected residues at amino terminal ends. Some observed topological incongruences in maximum likelihood trees of complete and curated Tollip data sets could be explained through horizontal transfers, evidenced by recombination detection. These results suggest that there is more to be researched and understood about this protein.

A novel ß-fructofuranosidase in Coleoptera: Characterization of a ß-fructofuranosidase from the sugarcane weevil, Sphenophorus levis.

ß-fructofuranosidases or invertases (EC 3.2.1.26) catalyze the hydrolysis of sucrose into fructose and glucose. ß-fructofuranosidases have been widely described in microorganisms, but were not known in the animal kingdom until very recently. There are studies reporting lepidopteran ß-fructofuranosidases, but no ß-fructofuranosidase gene sequence or encoding transcript has previously been identified in beetles. Considering the scarcity of functional studies on insect ß-fructofuranosidases and their apparent non-occurrence among coleopterans, the aim of the present study was to investigate the occurrence and characterize a ß-fructofuranosidase transcript identified in a cDNA library from the sugarcane weevil, Sphenophorus levis (Curculionidae). To validate that the ß-fructofuranosidase sequence (herein denominated Sl-ß-fruct) is indeed encoded by the S. levis genome, PCRs were performed using genomic DNA extracted from the larval fat body as well as DNA from the midgut with microbial content. Amplification of Sl-ß-fruct gene using larval fat body DNA indicated its presence in the insect's genomic DNA. The Sl-ß-fruct gene was cloned in Pichia pastoris to produce the recombinant enzyme (rSl-ß-fruct). Molecular weight of the recombinant protein was about 64 kDa, indicating possible glycosylation, since the theoretical weight was 54.8 kDa. The substrate specificity test revealed that rSl-ß-fruct hydrolyzes sucrose and raffinose, but not melibiose or maltose, thereby confirming invertase activity. The pH curve revealed greatest activity at pH 5.0, demonstrating rSl-ß-fruct to be an acidic ß-fructofuranosidase. Quantitative PCR (qRT-PCR) analyses indicated that the production of mRNA only occurs in the midgut and reaches the greatest expression level in 30-day-old larvae, which is the expected pattern for digestive enzymes. Chromatography of glycosidases from S. levis midguts showed two enzymes acting as ß-fructofuranosidase, indicating the presence of a Sl-ß-fruct isoform or a ß-fructofuranosidase from insect intestinal microbiota. Moreover, it was found that α-glucosidases do not act on sucrose hydrolysis. Phylogenetic analyses indicated this enzyme to be similar to enzymes found in other coleopteran and lepidopteran ß-fructofuranosidases, but also closely similar to bacterial enzymes, suggesting potential horizontal gene transfer. Despite this, the enzyme seems to be restricted to different groups of bacteria, which suggests distinct origin events. The present study expands the concept of the occurrence of ß-fructofuranosidase in insects. Despite the few descriptions of this gene in the animal kingdom, it is possible to state that ß-fructofuranosidase is crucial to the establishment of some insects throughout their evolutionary history, especially members of the Lepidoptera and Coleoptera clades.

Epigenetic modifications and their relation to caste and sex determination and adult division of labor in the stingless bee Melipona scutellaris

Abstract Stingless bees of the genus Melipona, have long been considered an enigmatic case among social insects for their mode of caste determination, where in addition to larval food type and quantity, the genotype also has a saying, as proposed over 50 years ago by Warwick E. Kerr. Several attempts have since tried to test his Mendelian two-loci/two-alleles segregation hypothesis, but only recently a single gene crucial for sex determination in bees was evidenced to be sex-specifically spliced and also caste-specifically expressed in a Melipona species. Since alternative splicing is frequently associated with epigenetic marks, and the epigenetic status plays a major role in setting the caste phenotype in the honey bee, we investigated here epigenetic chromatin modification in the stingless bee Melipona scutellaris. We used an ELISA-based methodology to quantify global methylation status and western blot assays to reveal histone modifications. The results evidenced DNA methylation/demethylation events in larvae and pupae, and significant differences in histone methylation and phosphorylation between newly emerged adult queens and workers. The epigenetic dynamics seen in this stingless bee species represent a new facet in the caste determination process in Melipona bees and suggest a possible mechanism that is likely to link a genotype component to the larval diet and adult social behavior of these bees.

Correlation between yield components in F6 soybean progenies derived from seven biparental crosses / Correlação entre caracteres em progênies F6 de soja provenientes de sete cruzamentos biparentais

Breeding programs strive to obtain cultivars with superior traits. The association between these characters, such as those estimated by correlation coeficients, are important to the breeding selection process. In this study we analyzed phenotypic and genotypic correlations between yield and other agronomic traits in soybean progenies derived from various crosses to support breeding line selection in the soybean breading program of the Federal University of Uberlandia (UFU). The experiment was carried out in the 2011/2012 growing season at the Capim Branco research station in Uberlandia, Minas Gerais, Brazil. Seventy-one soybean lines were evaluated using a randomized complete block design with three replications. In these plots plants were evaluated for the number of pods, grain yield per plant, number of days for flowering, plant height at flowering, number of days for maturity, plant height at maturity, height of the first pod insertion, number of grains per pod and grain yield per plant. The estimates for genotypic correlations were equal or greater than those for the phenotypic correlations for all the evaluated traits. This suggests a small effect of environmental factors on the trait expression. An exception of this behavior was the correlation between plant height at flowering and plant height at maturity. In this case, although the signs for the genotypic and phenotypic correlations were the same, the value for the last was greater than that for the former. The late flowering and late-maturing plants were taller. Also, selecting plants for higher number of pods ultimately contributed for the selection of plants with higher grain yield. The selection of late-flowering and late-maturing plants have yielded taller plantsin the same way that the selection for a higher pod number could promote the selection and improvement of grains yield trait.

O objetivo dos programas de melhoramento genético de qualquer cultura é a obtenção de cultivares que superem com vantagens as pré-existentes. Desse modo, o conhecimento das relações entre caracteres, tais como estimadas pelas correlações, tem apresentado grande importância para o sucesso dos processos de seleção no melhoramento. Os objetivos deste estudo foram avaliar as correlações fenotípicas e genotípicas entre caracteres agronômicos de linhagens de soja oriundas de diferentes cruzamentos, visando fornecer subsídios para fins de seleção a ser adotado no programa de melhoramento genético de soja da Universidade Federal de Uberlândia. O experimento foi conduzido em uma área experimental localizada na Fazenda Capim Branco, pertencente à Universidade Federal de Uberlândia, no município de Uberlândia, Minas Gerais, na safra 2011/2012. Foram avaliados 71 linhagens de soja em delineamento de blocos completos casualisados com três repetições. Avaliaram-se os caracteres: número total de vagens, produção de grãos por planta, número de dias para o florescimento, altura da planta no florescimento, número de dias para maturidade, altura da planta na maturidade, altura de inserção de primeira vagem, número de grãos por vagem e produção de grãos por planta. Para todos os caracteres, com exceção da correlação entre altura da planta no florescimento e maturidade em que a apesar de apresentarem mesma direção, a magnitude da correlação fenotípica superou a correlação genotípica, as correlações genotípicas foram superiores as fenotípicas e de igual sinal, indicando menor influência do ambiente na expressão dos caracteres. A seleção de plantas mais tardias no florescimento e maturidade, resulta em plantas mais altas nesses estádios de desenvolvimento, bem como a seleção de plantas com maior número de vagens favorece a seleção e melhoramento para o caractere produção de grãos por planta.