GO FEAT: a rapid web-based functional annotation tool for genomic and transcriptomic data

RNA sequencing and gene co-expression network of in vitro matured oocytes and blastocysts of buffalo

In reproductive technologies, uncovering the molecular aspects of oocyte and embryo competence under different conditions is crucial for refining protocols and enhancing efficiency. RNA-seq generates high-throughput data and provides transcriptomes that can undergo additional computational analyses. This study presented the transcriptomic profiles of in vitro matured oocytes and blastocysts produced in vitro from buffalo crossbred (Bubalus bubalis), coupled with gene co-expression and module preservation analysis. Cumulus Oophorus Complexes, obtained from slaughterhouse-derived ovaries, were subjected to in vitro maturation to yield metaphase II oocytes (616) or followed in vitro fertilization and culture to yield blastocysts for sequencing (526). Oocyte maturation (72%, ±3.34 sd) and embryo development (21.3%, ±4.18 sd) rates were obtained from three in vitro embryo production routines following standard protocols. Sequencing of 410 metaphase II oocytes and 70 hatched blastocysts (grade 1 and 2) identified a total of 13,976 genes, with 62% being ubiquitously expressed (8,649). Among them, the differentially expressed genes (4,153) and the strongly variable genes with the higher expression (fold-change above 11) were highlighted in oocytes (BMP15, UCHL1, WEE1, NLRPs, KPNA7, ZP2, and ZP4) and blastocysts (APOA1, KRT18, ANXA2, S100A14, SLC34A2, PRSS8 and ANXA2) as representative indicators of molecular quality. Additionally, genes exclusively found in oocytes (224) and blastocysts (2,200) with specific biological functions were identified. Gene co-expression network and module preservation analysis revealed strong preservation of functional modules related to exosome components, steroid metabolism, cell proliferation, and morphogenesis. However, cell cycle and amino acid transport modules exhibited weak preservation, which may reflect differences in embryo development kinetics and the activation of cell signaling pathways between buffalo and bovine. This comprehensive transcriptomic profile serves as a valuable resource for assessing the molecular quality of buffalo oocytes and embryos in future in vitro embryo production assays.

Genetic linkage mapping and quantitative trait locus (QTL) analysis of sweet basil (Ocimum basilicum L.) to identify genomic regions associated with cold tolerance and major volatiles

Chilling sensitivity is one of the greatest challenges affecting the marketability and profitability of sweet basil (Ocimum basilicum L.) in the US and worldwide. Currently, there are no sweet basils commercially available with significant chilling tolerance and traditional aroma profiles. This study was conducted to identify quantitative trait loci (QTLs) responsible for chilling tolerance and aroma compounds in a biparental mapping population, including the Rutgers advanced breeding line that served as a chilling tolerant parent, 'CB15', the chilling sensitive parent, 'Rutgers Obsession DMR' and 200 F2 individuals. Chilling tolerance was assessed by percent necrosis using machine learning and aroma profiling was evaluated using gas chromatography (GC) mass spectrometry (MS). Single nucleotide polymorphism (SNP) markers were generated from genomic sequences derived from double digestion restriction-site associated DNA sequencing (ddRADseq) and converted to genotype data using a reference genome alignment. A genetic linkage map was constructed and five statistically significant QTLs were identified in response to chilling temperatures with possible interactions between QTLs. The QTL on LG24 (qCH24) demonstrated the largest effect for chilling response and was significant in all three replicates. No QTLs were identified for linalool, as the population did not segregate sufficiently to detect this trait. Two significant QTLs were identified for estragole (also known as methyl chavicol) with only qEST1 on LG1 being significant in the multiple-QTL model (MQM). QEUC26 was identified as a significant QTL for eucalyptol (also known as 1,8-cineole) on LG26. These QTLs may represent key mechanisms for chilling tolerance and aroma in basil, providing critical knowledge for future investigation of these phenotypic traits and molecular breeding.

Recent advances in genome annotation and synthetic biology for the development of microbial chassis

This article provides an overview of microbial host selection, synthetic biology, genome annotation, metabolic modeling, and computational methods for predicting gene essentiality for developing a microbial chassis. This article focuses on lactic acid bacteria (LAB) as a microbial chassis and strategies for genome annotation of the LAB genome. As a case study, Lactococcus lactis is chosen based on its well-established therapeutic applications such as probiotics and oral vaccine development. In this article, we have delineated the strategies for genome annotations of lactic acid bacteria. These strategies also provide insights into streamlining genome reduction without compromising the functionality of the chassis and the potential for minimal genome chassis development. These insights underscore the potential for the development of efficient and sustainable synthetic biology systems using streamlined microbial chassis with minimal genomes.

Genomic and metabolomic insights into the antimicrobial compounds and plant growth-promoting potential of Bacillus velezensis Q-426

BACKGROUND

The Q-426 strain isolated from compost samples has excellent antifungal activities against a variety of plant pathogens. However, the complete genome of Q-426 is still unclear, which limits the potential application of Q-426.

RESULTS

Genome sequencing revealed that Q-426 contains a single circular chromosome 4,086,827 bp in length, with 4691 coding sequences and an average GC content of 46.3%. The Q-426 strain has a high degree of collinearity with B. velezensis FZB42, B. velezensis SQR9, and B. amyloliquefaciens DSM7, and the strain was reidentified as B. velezensis Q-426 based on the homology analysis results. Many genes in the Q-426 genome have plant growth-promoting activity, including the secondary metabolites of lipopeptides. Genome mining revealed 14 clusters and 732 genes encoding secondary metabolites with predicted functions, including the surfactin, iturin, and fengycin families. In addition, twelve lipopeptides (surfactin, iturin and fengycin) were successfully detected from the fermentation broth of B. velezensis Q-426 by ultra-high performance liquid chromatography-quadrupole time-of-flight mass spectrometry (UHPLC-QTOF-MS/MS), which is consistent with the genome analysis results. We found that Q-426 produced indole-3-acetic acid (IAA) at 1.56 mg/l on the third day of incubation, which might promote the growth of plants. Moreover, we identified eighteen volatile compounds (VOCs, including 2-heptanone, 6-methylheptan-2-one, 5-methylheptan-2-one, 2-nonanone, 2-decanone, 2-undecanone, 2-dodecanone, 2-tridecanone, 2-tetradecanone, 2-nonadecanone, pentadecanoic acid, oleic acid, dethyl phthalate, dibutyl phthalate, methyl (9E,12E)-octadeca-9,12-dienoate), pentadecane, (6E,10E)-1,2,3,4,4a,5,8,9,12,12a-decahydro-1,4-methanobenzo[10]annulene, and nonanal) based on gas chromatograph-mass spectrometer (GC/MS) results.

CONCLUSIONS

We mined secondary metabolite-related genes from the genome based on whole-genome sequence results. Our study laid the theoretical foundation for the development of secondary metabolites and the application of B. velezensis Q-426. Our findings provide insights into the genetic characteristics responsible for the bioactivities and potential application of B. velezensis Q-426 as a plant growth-promoting strain in ecological agriculture.

Targeting Essential Hypothetical Proteins of Pseudomonas aeruginosa PAO1 for Mining of Novel Therapeutics: An In Silico Approach

As an omnipresent opportunistic bacterium, Pseudomonas aeruginosa PAO1 is responsible for acute and chronic infection in immunocompromised individuals. Currently, this bacterium is on WHO's red list where new antibiotics are urgently required for the treatment. Finding essential genes and essential hypothetical proteins (EHP) can be crucial in identifying novel druggable targets and therapeutics. This study is aimed at characterizing these EHPs and analyzing subcellular and physiochemical properties, PPI network, nonhomologous analysis against humans, virulence factor and novel drug target prediction, and finally structural analysis of the identified target employing around 42 robust bioinformatics tools/databases, the output of which was evaluated using the ROC analysis. The study discovered 18 EHPs from 336 essential genes, with domain and functional annotation revealing that 50% of these proteins belong to the enzyme category. The majority are cytoplasmic and cytoplasmic membrane proteins, with half being stable proteins subjected to PPIs network analysis. The network contains 261 nodes and 269 edges for 9 proteins of interest, with 11 hubs containing at least three nodes each. Finally, a pipeline builder predicts 7 proteins with novel drug targets, 5 nonhomologous proteins against human proteome, human antitargets, and human gut flora, and 3 virulent proteins. Among these, homology modeling of NP_249450 and NP_251676 was done, and the Ramachandran plot analysis revealed that more than 94% of the residues were in the preferred region. By analyzing functional attributes and virulence characteristics, the findings of this study may facilitate the development of innovative antibacterial drug targets and drugs of Pseudomonas aeruginosa PAO1.

Unraveling the Secrets of a Double-Life Fungus by Genomics: Ophiocordyceps australis CCMB661 Displays Molecular Machinery for Both Parasitic and Endophytic Lifestyles

Ophiocordyceps australis (Ascomycota, Hypocreales, Ophiocordycipitaceae) is a classic entomopathogenic fungus that parasitizes ants (Hymenoptera, Ponerinae, Ponerini). Nonetheless, according to our results, this fungal species also exhibits a complete set of genes coding for plant cell wall degrading Carbohydrate-Active enZymes (CAZymes), enabling a full endophytic stage and, consequently, its dual ability to both parasitize insects and live inside plant tissue. The main objective of our study was the sequencing and full characterization of the genome of the fungal strain of O. australis (CCMB661) and its predicted secretome. The assembled genome had a total length of 30.31 Mb, N50 of 92.624 bp, GC content of 46.36%, and 8,043 protein-coding genes, 175 of which encoded CAZymes. In addition, the primary genes encoding proteins and critical enzymes during the infection process and those responsible for the host-pathogen interaction have been identified, including proteases (Pr1, Pr4), aminopeptidases, chitinases (Cht2), adhesins, lectins, lipases, and behavioral manipulators, such as enterotoxins, Protein Tyrosine Phosphatases (PTPs), and Glycoside Hydrolases (GHs). Our findings indicate that the presence of genes coding for Mad2 and GHs in O. australis may facilitate the infection process in plants, suggesting interkingdom colonization. Furthermore, our study elucidated the pathogenicity mechanisms for this Ophiocordyceps species, which still is scarcely studied.

Bacteriobiota of the Cave Church of Sts. Peter and Paul in Serbia-Culturable and Non-Culturable Communities' Assessment in the Bioconservation Potential of a Peculiar Fresco Painting

The principal aim of this study was to determine bacterial diversity within the Cave Church of Sts. Peter and Paul, via culturable and non-culturable approaches, and elucidate the antifungal potential of autochthonous antagonistic bacterial isolates against biodeteriogenic fungi. Furthermore, whole-genome sequencing of selected bacterial antagonists and the analysis of genes included in the synthesis of secondary metabolites were performed. With the highest RA values, determined in metabarcoding analysis, phyla Actinobacteriota (12.08-54.00%) and Proteobacteria (25.34-44.97%) dominated most of the samples. A total of 44 different species, out of 96 obtained isolates, were determined as part of the culturable bacteriobiota, with the predominance of species from the genus Bacillus. Bacillus simplex was the only isolated species simultaneously present in all investigated substrata within the church. The best antagonistic activity against 10 biodeteriogenic fungi was documented for Streptomyces anulatus, followed by Bacillus altitudinis, Chryseobacterium viscerum, and Streptomyces sp. with their highest PGI% values ranging of from 55.9% to 80.9%. These promising results indicate that characterized bacteria are excellent candidates for developing biocontrol strategies for suppressing deteriogenic fungi responsible for the deterioration of investigated fresco painting. Finally, isolate 11-11MM, characterized as Streptomyces sp., represents a new species for science prompting the need for further study.

Phylogenomics and gene selection in Aspergillus welwitschiae: Possible implications in the pathogenicity in Agave sisalana

Aspergillus welwitschiae causes bole rot disease in sisal (Agave sisalana and related species) which affects the production of natural fibers in Brazil, the main worldwide producer of sisal fibers. This fungus is a saprotroph with a broad host range. Previous research established A. welwitschiae as the only causative agent of bole rot in the field, but little is known about the evolution of this species and its strains. In this work, we performed a comparative genomics analysis of 40 Aspergillus strains. We show the conflicting molecular identity of this species, with one sisal-infecting strain sharing its last common ancestor with Aspergillus niger, having diverged only 833 thousand years ago. Furthermore, our analysis of positive selection reveals sites under selection in genes coding for siderophore transporters, Sodium‑calcium exchangers, and Phosphatidylethanolamine-binding proteins (PEBPs). Herein, we discuss the possible impacts of these gene functions on the pathogenicity in sisal.

Aeromonas allosaccharophila Strain AE59-TE2 Is Highly Antagonistic towards Multidrug-Resistant Human Pathogens, What Does Its Genome Tell Us?

Multidrug-resistant bacteria are of critical importance and a problem for human health and food preservation; the discovery of new antimicrobial substances to control their proliferation is part of the solution. This work reports on 57 antagonistic Aeromonas strains, of which 38 strains were antagonistic towards problematic human pathogens. The genome of the most antagonistic strain was sequenced and identified as Aeromonas allosaccharophila. Its genome was fully annotated and mined for genes that might explain that activity. Strain AE59-TE was antagonistic toward clinically relevant gram-negative and gram-positive multidrug-resistant bacteria, including Klebsiella pneumoniae KPC, Escherichia coli ESBL, Salmonella typhimurium, and Staphylococcus aureus MRSA. Strain AE59-TE2 was identified by multilocus sequence analysis. Genome mining identified four genes homologous to the bacteriocin, zoocin A from Streptococcus equi and a gene 98% similar to cvpA linked to colicin V production. A. allosaccharophila strain AE59-TE2 produced antimicrobial activity against a broad range of bacteria, including important gram-negative bacteria, not typically targeted by bacteriocins. Herewere described novel zoocin genes that are promising for industrial applications in the food and health sectors. Interesting and important antagonistic activity is described combined with the first detailed genomic analysis of the species Aeromonas allosaccharophila.

Integrative transcriptomic and metabolomic analyses unveil tanshinone biosynthesis in Salvia miltiorrhiza root under N starvation stress

Salvia miltiorrhiza is a model plant for Chinese herbal medicine with significant pharmacologic effects due to its tanshinone components. Our previous study indicated that nitrogen starvation stress increased its tanshinone content. However, the molecular mechanism of this low nitrogen-induced tanshinone biosynthesis is still unclear. Thus, this study aimed to elucidate the molecular mechanism of tanshinone biosynthesis in S. miltiorrhiza under different N conditions [N-free (N0), low-N (Nl), and full-N (Nf, as control) conditions] by using transcriptome and metabolome analyses. Our results showed 3,437 and 2,274 differentially expressed unigenes between N0 and Nf as well as Nl and Nf root samples, respectively. N starvation (N0 and Nl) promoted the expression of the genes involved in the MVA and MEP pathway of tanshinone and terpenoid backbone biosynthesis. Gene ontology and KEGG analyses revealed that terpenoid backbone biosynthesis, hormone signal transduction, and phenylpropanoid biosynthesis were promoted under N starvation conditions, whereas starch and sucrose metabolisms, nitrogen and phosphorus metabolisms, as well as membrane development were inhibited. Furthermore, metabolome analysis showed that metabolite compounds and biosynthesis of secondary metabolites were upregulated. This study provided a novel insight into the molecular mechanisms of tanshinone production in S. miltiorrhiza in response to nitrogen stress.

Unraveling the Genomic Potential of the Thermophilic Bacterium Anoxybacillus flavithermus from an Antarctic Geothermal Environment

Antarctica is a mosaic of extremes. It harbors active polar volcanoes, such as Deception Island, a marine stratovolcano having notable temperature gradients over very short distances, with the temperature reaching up to 100 °C near the fumaroles and subzero temperatures being noted in the glaciers. From the sediments of Deception Island, we isolated representatives of the genus Anoxybacillus, a widely spread genus that is mainly encountered in thermophilic environments. However, the phylogeny of this genus and its adaptive mechanisms in the geothermal sites of cold environments remain unknown. To the best of our knowledge, this is the first study to unravel the genomic features and provide insights into the phylogenomics and metabolic potential of members of the genus Anoxybacillus inhabiting the Antarctic thermophilic ecosystem. Here, we report the genome sequencing data of seven A. flavithermus strains isolated from two geothermal sites on Deception Island, Antarctic Peninsula. Their genomes were approximately 3.0 Mb in size, had a G + C ratio of 42%, and were predicted to encode 3500 proteins on average. We observed that the strains were phylogenomically closest to each other (Average Nucleotide Identity (ANI) > 98%) and to A. flavithermus (ANI 95%). In silico genomic analysis revealed 15 resistance and metabolic islands, as well as genes related to genome stabilization, DNA repair systems against UV radiation threats, temperature adaptation, heat- and cold-shock proteins (Csps), and resistance to alkaline conditions. Remarkably, glycosyl hydrolase enzyme-encoding genes, secondary metabolites, and prophage sequences were predicted, revealing metabolic and cellular capabilities for potential biotechnological applications.

Exploration of Streptococcus core genome to reveal druggable targets and novel therapeutics against S. pneumoniae

Streptococcus pneumoniae (S. pneumoniae), the major etiological agent of community-acquired pneumonia (CAP) contributes significantly to the global burden of infectious diseases which is getting resistant day by day. Nearly 30% of the S. pneumoniae genomes encode hypothetical proteins (HPs), and better understandings of these HPs in virulence and pathogenicity plausibly decipher new treatments. Some of the HPs are present across many Streptococcus species, systematic assessment of these unexplored HPs will disclose prospective drug targets. In this study, through a stringent bioinformatics analysis of the core genome and proteome of S. pneumoniae PCS8235, we identified and analyzed 28 HPs that are common in many Streptococcus species and might have a potential role in the virulence or pathogenesis of the bacteria. Functional annotations of the proteins were conducted based on the physicochemical properties, subcellular localization, virulence prediction, protein-protein interactions, and identification of essential genes, to find potentially druggable proteins among 28 HPs. The majority of the HPs are involved in bacterial transcription and translation. Besides, some of them were homologs of enzymes, binding proteins, transporters, and regulators. Protein-protein interactions revealed HP PCS8235_RS05845 made the highest interactions with other HPs and also has TRP structural motif along with virulent and pathogenic properties indicating it has critical cellular functions and might go under unconventional protein secretions. The second highest interacting protein HP PCS8235_RS02595 interacts with the Regulator of chromosomal segregation (RocS) which participates in chromosome segregation and nucleoid protection in S. pneumoniae. In this interacting network, 54% of protein members have virulent properties and 40% contain pathogenic properties. Among them, most of these proteins circulate in the cytoplasmic area and have hydrophilic properties. Finally, molecular docking and dynamics simulation demonstrated that the antimalarial drug Artenimol can act as a drug repurposing candidate against HP PCS8235_RS 04650 of S. pneumoniae. Hence, the present study could aid in drugs against S. pneumoniae.

The Complete Genome of the “Flavescence Dorée” Phytoplasma Reveals Characteristics of Low Genome Plasticity

Members of the genus ‘Candidatus Phytoplasma’ are obligate intracellular bacteria restricted to phloem sieve elements and are able to colonize several tissues and the hemolymph in their insect vectors. The current unfeasibility of axenic culture and the low complexity of genomic sequences are obstacles in assembling complete chromosomes. Here, a method combining pathogen DNA enrichment from infected insects and dual deep-sequencing technologies was used to obtain the complete genome of a phytoplasma causing Grapevine Flavescence dorée. The de novo assembly generated a circular chromosome of 654,223 bp containing 506 protein-coding genes. Quality assessment of the draft showed a high degree of completeness. Comparative analysis with other phytoplasmas revealed the absence of potential mobile units and a reduced amount of putative phage-derived segments, suggesting a low genome plasticity. Phylogenetic analyses identified Candidatus Phytoplasma ziziphi as the closest fully sequenced relative. The “Flavescence dorée” phytoplasma strain CH genome also encoded for several putative effector proteins potentially playing a role in pathogen virulence. The availability of this genome provides the basis for the study of the pathogenicity mechanisms and evolution of the Flavescence dorée phytoplasma.

Polyphasic Analysis Reveals Potential Petroleum Hydrocarbon Degradation and Biosurfactant Production by Rare Biosphere Thermophilic Bacteria From Deception Island, an Active Antarctic Volcano

Extreme temperature gradients in polar volcanoes are capable of selecting different types of extremophiles. Deception Island is a marine stratovolcano located in maritime Antarctica. The volcano has pronounced temperature gradients over very short distances, from as high as 100°C in the fumaroles to subzero next to the glaciers. These characteristics make Deception a promising source of a variety of bioproducts for use in different biotechnological areas. In this study, we isolated thermophilic bacteria from sediments in fumaroles at two geothermal sites on Deception Island with temperatures between 50 and 100°C, to evaluate the potential capacity of these bacteria to degrade petroleum hydrocarbons and produce biosurfactants under thermophilic conditions. We isolated 126 thermophilic bacterial strains and identified them molecularly as members of genera Geobacillus, Anoxybacillus, and Brevibacillus (all in phylum Firmicutes). Seventy-six strains grew in a culture medium supplemented with crude oil as the only carbon source, and 30 of them showed particularly good results for oil degradation. Of 50 strains tested for biosurfactant production, 13 showed good results, with an emulsification index of 50% or higher of a petroleum hydrocarbon source (crude oil and diesel), emulsification stability at 100°C, and positive results in drop-collapse, oil spreading, and hemolytic activity tests. Four of these isolates showed great capability of degrade crude oil: FB2_38 (Geobacillus), FB3_54 (Geobacillus), FB4_88 (Anoxybacillus), and WB1_122 (Geobacillus). Genomic analysis of the oil-degrading and biosurfactant-producer strain FB4_88 identified it as Anoxybacillus flavithermus, with a high genetic and functional diversity potential for biotechnological applications. These initial culturomic and genomic data suggest that thermophilic bacteria from this Antarctic volcano have potential applications in the petroleum industry, for bioremediation in extreme environments and for microbial enhanced oil recovery (MEOR) in reservoirs. In addition, recovery of small-subunit rRNA from metagenomes of Deception Island showed that Firmicutes is not among the dominant phyla, indicating that these low-abundance microorganisms may be important for hydrocarbon degradation and biosurfactant production in the Deception Island volcanic sediments.

Transcriptomic analysis of chloride tolerance in Leptospirillum ferriphilum DSM 14647 adapted to NaCl

Chloride ions are toxic for most acidophilic microorganisms. In this study, the chloride tolerance mechanisms in the acidophilic iron-oxidizing bacterium Leptospirillum ferriphilum DSM 14647 adapted to 180 mM NaCl were investigated by a transcriptomic approach. Results showed that 99 genes were differentially expressed in the adapted versus the non-adapted cultures, of which 69 and 30 were significantly up-regulated or down-regulated, respectively. Genes that were up-regulated include carbonic anhydrase, cytochrome c oxidase (ccoN) and sulfide:quinone reductase (sqr), likely involved in intracellular pH regulation. Towards the same end, the cation/proton antiporter CzcA (czcA) was down-regulated. Adapted cells showed a higher oxygen consumption rate (2.2 x 10⁻⁹ ppm O₂ s^-1cell^-1) than non-adapted cells (1.2 x 10⁻⁹ ppm O₂ s^-1cell^-1). Genes coding for the antioxidants flavohemoprotein and cytochrome c peroxidase were also up-regulated. Measurements of the intracellular reactive oxygen species (ROS) level revealed that adapted cells had a lower level than non-adapted cells, suggesting that detoxification of ROS could be an important strategy to withstand NaCl. In addition, data analysis revealed the up-regulation of genes for Fe-S cluster biosynthesis (iscR), metal reduction (merA) and activation of a cellular response mediated by diffusible signal factors (DSFs) and the second messenger c-di-GMP. Several genes related to the synthesis of lipopolysaccharide and peptidoglycan were consistently down-regulated. Unexpectedly, the genes ectB, ectC and ectD involved in the biosynthesis of the compatible solutes (hydroxy)ectoine were also down-regulated. In line with these findings, although hydroxyectoine reached 20 nmol mg^-1 of wet biomass in non-adapted cells, it was not detected in L. ferriphilum adapted to NaCl, suggesting that this canonical osmotic stress response was dispensable for salt adaptation. Differentially expressed transcripts and experimental validations suggest that adaptation to chloride in acidophilic microorganisms involves a multifactorial response that is different from the response in other bacteria studied.

Draft Genomes of Six Philippine Erwinia mallotivora Isolates: Comparative Genomics and Genome-Wide Analysis of Candidate Secreted Proteins

Erwinia mallotivora is one of the most important bacterial pathogens of papaya and causes bacterial crown rot disease in the Philippines. In this paper, we present the draft genome sequences of six Philippine E. mallotivora isolates to provide insights into the genes involved in host-pathogen interactions and compare their genomes to other Erwinia species. The genomes were sequenced using Illumina Miseq platform. The draft whole-genome assemblies of the E. mallotivora isolates are composed of 36-64 contigs with N50 value ranging from 285 to 332 kbp and cover 96.2-100% of the estimated genome size. Structural genome annotation of these assemblies has predicted 4489-4749 protein-coding genes. Comparative genomic analysis using orthologous gene sets led to the identification of conserved genes within the genus and species-specific gene orthologous groups, which collectively provide a baseline for functional genomic studies to determine genes affecting virulence and host specificity. Secreted proteins of E. mallotivora were also predicted and characterized to unravel putative genes involved in plant-pathogen interactions. This study provides the first draft whole-genome sequences of Philippine isolates of E. mallotivora, thus expanding the genomic knowledge for this species in comparison with other members of the genus Erwinia.

Genome-Wide Survey and Development of the First Microsatellite Markers Database (AnCorDB) in Anemone coronaria L

Anemone coronaria L. (2n = 2x = 16) is a perennial, allogamous, highly heterozygous plant marketed as a cut flower or in gardens. Due to its large genome size, limited efforts have been made in order to develop species-specific molecular markers. We obtained the first draft genome of the species by Illumina sequencing an androgenetic haploid plant of the commercial line “MISTRAL^® Magenta”. The genome assembly was obtained by applying the MEGAHIT pipeline and consisted of 2 × 10⁶ scaffolds. The SciRoKo SSR (Simple Sequence Repeats)-search module identified 401.822 perfect and 188.987 imperfect microsatellites motifs. Following, we developed a user-friendly “Anemone coronaria Microsatellite DataBase” (AnCorDB), which incorporates the Primer3 script, making it possible to design couples of primers for downstream application of the identified SSR markers. Eight genotypes belonging to eight cultivars were used to validate 62 SSRs and a subset of markers was applied for fingerprinting each cultivar, as well as to assess their intra-cultivar variability. The newly developed microsatellite markers will find application in Breeding Rights disputes, developing genetic maps, marker assisted breeding (MAS) strategies, as well as phylogenetic studies.

Investigating the Functional Role of Hypothetical Proteins From an Antarctic Bacterium Pseudomonas sp. Lz4W: Emphasis on Identifying Proteins Involved in Cold Adaptation

Exploring the molecular mechanisms behind bacterial adaptation to extreme temperatures has potential biotechnological applications. In the present study, Pseudomonas sp. Lz4W, a Gram-negative psychrophilic bacterium adapted to survive in Antarctica, was selected to decipher the molecular mechanism underlying the cold adaptation. Proteome analysis of the isolates grown at 4°C was performed to identify the proteins and pathways that are responsible for the adaptation. However, many proteins from the expressed proteome were found to be hypothetical proteins (HPs), whose function is unknown. Investigating the functional roles of these proteins may provide additional information in the biological understanding of the bacterial cold adaptation. Thus, our study aimed to assign functions to these HPs and understand their role at the molecular level. We used a structured insilico workflow combining different bioinformatics tools and databases for functional annotation. Pseudomonas sp. Lz4W genome (CP017432, version 1) contains 4493 genes and 4412 coding sequences (CDS), of which 743 CDS were annotated as HPs. Of these, from the proteome analysis, 61 HPs were found to be expressed consistently at the protein level. The amino acid sequences of these 61 HPs were submitted to our workflow and we could successfully assign a function to 18 HPs. Most of these proteins were predicted to be involved in biological mechanisms of cold adaptations such as peptidoglycan metabolism, cell wall organization, ATP hydrolysis, outer membrane fluidity, catalysis, and others. This study provided a better understanding of the functional significance of HPs in cold adaptation of Pseudomonas sp. Lz4W. Our approach emphasizes the importance of addressing the “hypothetical protein problem” for a thorough understanding of mechanisms at the cellular level, as well as, provided the assessment of integrating proteomics methods with various annotation and curation approaches to characterize hypothetical or uncharacterized protein data. The MS proteomics data generated from this study has been deposited to the ProteomeXchange through PRIDE with the dataset identifier–PXD029741.

The Space-Exposed Kombucha Microbial Community Member Komagataeibacter oboediens Showed Only Minor Changes in Its Genome After Reactivation on Earth

Komagataeibacter is the dominant taxon and cellulose-producing bacteria in the Kombucha Microbial Community (KMC). This is the first study to isolate the K. oboediens genome from a reactivated space-exposed KMC sample and comprehensively characterize it. The space-exposed genome was compared with the Earth-based reference genome to understand the genome stability of K. oboediens under extraterrestrial conditions during a long time. Our results suggest that the genomes of K. oboediens IMBG180 (ground sample) and K. oboediens IMBG185 (space-exposed) are remarkably similar in topology, genomic islands, transposases, prion-like proteins, and number of plasmids and CRISPR-Cas cassettes. Nonetheless, there was a difference in the length of plasmids and the location of cas genes. A small difference was observed in the number of protein coding genes. Despite these differences, they do not affect any genetic metabolic profile of the cellulose synthesis, nitrogen-fixation, hopanoid lipids biosynthesis, and stress-related pathways. Minor changes are only observed in central carbohydrate and energy metabolism pathways gene numbers or sequence completeness. Altogether, these findings suggest that K. oboediens maintains its genome stability and functionality in KMC exposed to the space environment most probably due to the protective role of the KMC biofilm. Furthermore, due to its unaffected metabolic pathways, this bacterial species may also retain some promising potential for space applications.

MEDUSA: A Pipeline for Sensitive Taxonomic Classification and Flexible Functional Annotation of Metagenomic Shotgun Sequences

Metagenomic studies unravel details about the taxonomic composition and the functions performed by microbial communities. As a complete metagenomic analysis requires different tools for different purposes, the selection and setup of these tools remain challenging. Furthermore, the chosen toolset will affect the accuracy, the formatting, and the functional identifiers reported in the results, impacting the results interpretation and the biological answer obtained. Thus, we surveyed state-of-the-art tools available in the literature, created simulated datasets, and performed benchmarks to design a sensitive and flexible metagenomic analysis pipeline. Here we present MEDUSA, an efficient pipeline to conduct comprehensive metagenomic analyses. It performs preprocessing, assembly, alignment, taxonomic classification, and functional annotation on shotgun data, supporting user-built dictionaries to transfer annotations to any functional identifier. MEDUSA includes several tools, as fastp, Bowtie2, DIAMOND, Kaiju, MEGAHIT, and a novel tool implemented in Python to transfer annotations to BLAST/DIAMOND alignment results. These tools are installed via Conda, and the workflow is managed by Snakemake, easing the setup and execution. Compared with MEGAN 6 Community Edition, MEDUSA correctly identifies more species, especially the less abundant, and is more suited for functional analysis using Gene Ontology identifiers.

Hybrid Assembly Improves Genome Quality and Completeness of Trametes villosa CCMB561 and Reveals a Huge Potential for Lignocellulose Breakdown

Trametes villosa is a wood-decaying fungus with great potential to be used in the bioconversion of agro-industrial residues and to obtain high-value-added products, such as biofuels. Nonetheless, the lack of high-quality genomic data hampers studies investigating genetic mechanisms and metabolic pathways in T. villosa, hindering its application in industry. Herein, applying a hybrid assembly pipeline using short reads (Illumina HiSeq) and long reads (Oxford Nanopore MinION), we obtained a high-quality genome for the T. villosa CCMB561 and investigated its genetic potential for lignocellulose breakdown. The new genome possesses 143 contigs, N50 of 1,009,271 bp, a total length of 46,748,415 bp, 14,540 protein-coding genes, 22 secondary metabolite gene clusters, and 426 genes encoding Carbohydrate-Active enzymes. Our CAZome annotation and comparative genomic analyses of nine Trametes spp. genomes revealed T. villosa CCMB561 as the species with the highest number of genes encoding lignin-modifying enzymes and a wide array of genes encoding proteins for the breakdown of cellulose, hemicellulose, and pectin. These results bring to light the potential of this isolate to be applied in the bioconversion of lignocellulose and will support future studies on the expression, regulation, and evolution of genes, proteins, and metabolic pathways regarding the bioconversion of lignocellulosic residues.

Integration of Aspergillus niger transcriptomic profile with metabolic model identifies potential targets to optimise citric acid production from lignocellulosic hydrolysate

BACKGROUND

Citric acid is typically produced industrially by Aspergillus niger-mediated fermentation of a sucrose-based feedstock, such as molasses. The fungus Aspergillus niger has the potential to utilise lignocellulosic biomass, such as bagasse, for industrial-scale citric acid production, but realising this potential requires strain optimisation. Systems biology can accelerate strain engineering by systematic target identification, facilitated by methods for the integration of omics data into a high-quality metabolic model. In this work, we perform transcriptomic analysis to determine the temporal expression changes during fermentation of bagasse hydrolysate and develop an evolutionary algorithm to integrate the transcriptomic data with the available metabolic model to identify potential targets for strain engineering.

RESULTS

The novel integrated procedure matures our understanding of suboptimal citric acid production and reveals potential targets for strain engineering, including targets consistent with the literature such as the up-regulation of citrate export and pyruvate carboxylase as well as novel targets such as the down-regulation of inorganic diphosphatase.

CONCLUSIONS

In this study, we demonstrate the production of citric acid from lignocellulosic hydrolysate and show how transcriptomic data across multiple timepoints can be coupled with evolutionary and metabolic modelling to identify potential targets for further engineering to maximise productivity from a chosen feedstock. The in silico strategies employed in this study can be applied to other biotechnological goals, assisting efforts to harness the potential of microorganisms for bio-based production of valuable chemicals.

Functional Analysis of Hypothetical Proteins of Vibrio parahaemolyticus Reveals the Presence of Virulence Factors and Growth-Related Enzymes With Therapeutic Potential

Vibrio parahaemolyticus, an aquatic pathogen, is a major concern in the shrimp aquaculture industry. Several strains of this pathogen are responsible for causing acute hepatopancreatic necrosis disease as well as other serious illness, both of which result in severe economic losses. The genome sequence of two pathogenic strains of V. parahaemolyticus, MSR16 and MSR17, isolated from Bangladesh, have been reported to gain a better understanding of their diversity and virulence. However, the prevalence of hypothetical proteins (HPs) makes it challenging to obtain a comprehensive understanding of the pathogenesis of V. parahaemolyticus. The aim of the present study is to provide a functional annotation of the HPs to elucidate their role in pathogenesis employing several in silico tools. The exploration of protein domains and families, similarity searches against proteins with known function, gene ontology enrichment, along with protein-protein interaction analysis of the HPs led to the functional assignment with a high level of confidence for 656 proteins out of a pool of 2631 proteins. The in silico approach used in this study was important for accurately assigning function to HPs and inferring interactions with proteins with previously described functions. The HPs with function predicted were categorized into various groups such as enzymes involved in small-compound biosynthesis pathway, iron binding proteins, antibiotics resistance proteins, and other proteins. Several proteins with potential druggability were identified among them. In addition, the HPs were investigated in search of virulent factors, which led to the identification of proteins that have the potential to be exploited as vaccine candidate. The findings of the study will be effective in gaining a better understanding of the molecular mechanisms of bacterial pathogenesis. They may also provide an insight into the process of evaluating promising targets for the development of drugs and vaccines against V. parahaemolyticus.

Genome-scale protein interaction network construction and topology analysis of functional hypothetical proteins in Helicobacter pylori divulges novel therapeutic targets

The emergence and spread of multi-drug resistance among Helicobacter pylori (H. pylori) strain raise more stakes for genetic research for discovering new drugs. The quantity of uncharacterized hypothetical proteins in the genome may provide an opportunity to explore their property and promulgation could act as a platform for designing the drugs, making them an intriguing genetic target. In this context, the present study aims to identify the key hypothetical proteins (HPs) and their biological regulatory processes in H. pylori. This investigation could provide a foundation to establish the molecular connectivity among the pathways using topological analysis of the protein interaction networks (PINs). The giant network derived from the extended network has 374 nodes connected via 925 edges. A total of 43 proteins with high betweenness centrality (BC), 54 proteins with a large degree, and 23 proteins with high BC and large degrees have been identified. HP 1479, HP 0056, HP 1481, HP 1021, HP 0043, HP 1019, gmd, flgA, HP 0472, HP 1486, HP 1478, and HP 1473 are categorized as hub nodes because they have a higher number of direct connections and are potentially more important in understanding HP's molecular interactions. The pathway enrichment analysis of the network clusters revealed significant involvement of HPs in pathways such as flagellar assembly, bacterial chemotaxis and lipopolysaccharide biosynthesis. This comprehensive computational study revealed HP's functional role and its druggability characteristics, which could be useful in the development of drugs to combat H. pylori infections.

Proteins from the core genome of Corynebacterium ulcerans respond for pathogenicity and reveal promising vaccine targets for diphtheria

Corynebacterium ulcerans is an emerging pathogen able to transmit the acute infection diphtheria to humans. Although there is a well-established vaccine based on the toxin produced by Corynebacterium diphtheriae, another species of this genus known to cause the disease, there is still no vaccine formulations described for C. ulcerans; this fact contributes to the increase in cases of infection that has been observed. In this study, we want to provide information at the genomic level of this bacterium in order to suggest proteins as possible vaccine targets. We carried out an in silico prospection of vaccine candidates through reverse vaccinology for targets that exhibit antigenic potential against diphtheria. We found important virulence factors, such as adhesion-related ones, that are responsible for pathogen-host interaction after infection, but we did not find the diphtheria toxin, which is the main component of the currently available vaccine. This study provides detailed information about the exoproteome and hypothetical proteins from the core genome of C. ulcerans, suggesting vaccine targets to be further tested in vitro for the development of a new vaccine against diphtheria.

Host Colonization as a Major Evolutionary Force Favoring the Diversity and the Emergence of the Worldwide Multidrug-Resistant Escherichia coli ST131

The emergence of multidrug-resistant Escherichia coli ST131 is a major worldwide public health problem in humans. According to the “one health” approach, this study investigated animal reservoirs of ST131, their relationships with human strains, and the genetic features associated with host colonization. High-quality genomes originating from human, avian, and canine hosts were classified on the basis of their accessory gene content using pangenomic. Pangenomic clusters and subclusters were specifically and significantly associated with hosts. The functions of clustering accessory genes were mainly enriched in functions involved in DNA acquisition, interactions, and virulence (e.g., pathogenesis, response to biotic stimulus and interaction between organisms). Accordingly, networks of cooccurrent host interaction factors were significantly associated with the pangenomic clusters and the originating hosts. The avian strains exhibited a specific content in virulence factors. Rarely found in humans, they corresponded to pathovars responsible for severe human infections. An emerging subcluster significantly associated with both human and canine hosts was evidenced. This ability to significantly colonize canine hosts in addition to humans was associated with a specific content in virulence factors (VFs) and metabolic functions encoded by a new pathogenicity island in ST131 and an improved fitness that is probably involved in its emergence. Overall, VF content, unlike the determinants of antimicrobial resistance, appeared as a key actor of bacterial host adaptation. The host dimension emerges as a major driver of genetic evolution that shapes ST131 genome, enhances its diversity, and favors its dissemination.

An Integrated Database of Small RNAs and Their Interplay With Transcriptional Gene Regulatory Networks in Corynebacteria

Small RNAs (sRNAs) are one of the key players in the post-transcriptional regulation of bacterial gene expression. These molecules, together with transcription factors, form regulatory networks and greatly influence the bacterial regulatory landscape. Little is known concerning sRNAs and their influence on the regulatory machinery in the genus Corynebacterium, despite its medical, veterinary and biotechnological importance. Here, we expand corynebacterial regulatory knowledge by integrating sRNAs and their regulatory interactions into the transcriptional regulatory networks of six corynebacterial species, covering four human and animal pathogens, and integrate this data into the CoryneRegNet database. To this end, we predicted sRNAs to regulate 754 genes, including 206 transcription factors, in corynebacterial gene regulatory networks. Amongst them, the sRNA Cd-NCTC13129-sRNA-2 is predicted to directly regulate ydfH, which indirectly regulates 66 genes, including the global regulator glxR in C. diphtheriae. All of the sRNA-enriched regulatory networks of the genus Corynebacterium have been made publicly available in the newest release of CoryneRegNet(www.exbio.wzw.tum.de/coryneregnet/) to aid in providing valuable insights and to guide future experiments.

Mantis: flexible and consensus-driven genome annotation

Background

The rapid development of the (meta-)omics fields has produced an unprecedented amount of high-resolution and high-fidelity data. Through the use of these datasets we can infer the role of previously functionally unannotated proteins from single organisms and consortia. In this context, protein function annotation can be described as the identification of regions of interest (i.e., domains) in protein sequences and the assignment of biological functions. Despite the existence of numerous tools, challenges remain in terms of speed, flexibility, and reproducibility. In the big data era, it is also increasingly important to cease limiting our findings to a single reference, coalescing knowledge from different data sources, and thus overcoming some limitations in overly relying on computationally generated data from single sources.

Results

We implemented a protein annotation tool, Mantis, which uses database identifiers intersection and text mining to integrate knowledge from multiple reference data sources into a single consensus-driven output. Mantis is flexible, allowing for the customization of reference data and execution parameters, and is reproducible across different research goals and user environments. We implemented a depth-first search algorithm for domain-specific annotation, which significantly improved annotation performance compared to sequence-wide annotation. The parallelized implementation of Mantis results in short runtimes while also outputting high coverage and high-quality protein function annotations.

Conclusions

Mantis is a protein function annotation tool that produces high-quality consensus-driven protein annotations. It is easy to set up, customize, and use, scaling from single genomes to large metagenomes. Mantis is available under the MIT license at https://github.com/PedroMTQ/mantis.

Characteristics of the AT-Hook Motif Containing Nuclear Localized (AHL) Genes in Carrot Provides Insight into Their Role in Plant Growth and Storage Root Development

The AT-hook motif containing nuclear localized (AHL) gene family, controlling various developmental processes, is conserved in land plants. They comprise Plant and Prokaryote Conserved (PPC) domain and one or two AT-hook motifs. DcAHLc1 has been proposed as a candidate gene governing the formation of the carrot storage root. We identified and in-silico characterized carrot AHL proteins, performed phylogenetic analyses, investigated their expression profiles and constructed gene coexpression networks. We found 47 AHL genes in carrot and grouped them into two clades, A and B, comprising 29 and 18 genes, respectively. Within Clade-A, we distinguished three subclades, one of them grouping noncanonical AHLs differing in their structure (two PPC domains) and/or cellular localization (not nucleus). Coexpression network analysis attributed AHLs expressed in carrot roots into four of the 72 clusters, some of them showing a large number of interactions. Determination of expression profiles of AHL genes in various tissues and samples provided basis to hypothesize on their possible roles in the development of the carrot storage root. We identified a group of rapidly evolving noncanonical AHLs, possibly differing functionally from typical AHLs, as suggested by their expression profiles and their predicted cellular localization. We pointed at several AHLs likely involved in the development of the carrot storage root.

Genomic Resources for the North American Water Vole (Microtus richardsoni) and the Montane Vole (Microtus montanus)

Voles of the genus Microtus are important research organisms, yet genomic resources are lacking. Such resources would benefit future studies of immunology, phylogeography, cryptic diversity, and more. We sequenced and assembled nuclear genomes from two subspecies of water vole (Microtus richardsoni) and from the montane vole (Microtus montanus). The water vole genomes were sequenced with Illumina and 10× Chromium plus Illumina sequencing, resulting in assemblies with ∼1600,000 and ∼30,000 scaffolds, respectively. The montane vole was also assembled into ∼13,000 scaffolds using Illumina sequencing. Mitochondrial genome assemblies were also performed for both species. Structural and functional annotation for the best water vole nuclear genome resulted in ∼24,500 annotated genes, with 83% of these having functional annotations. Assembly quality statistics for our nuclear assemblies fall within the range of genomes previously published in the genus Microtus, making the water vole and montane vole genomes useful additions to currently available genomic resources.

Bioinformatic analysis and validation of microRNA-508-3p as a protective predictor by targeting NR4A3/MEK axis in pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) featured a debilitating progressive disorder. Here, we intend to determine diagnosis‐valuable biomarkers for PAH and decode the fundamental mechanisms of the biological function of these markers. Two mRNA microarray profiles (GSE70456 and GSE117261) and two microRNA microarray profiles (GSE55427 and GSE67597) were mined from the Gene Expression Omnibus platform. Then, we identified the differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs), respectively. Besides, we investigated online miRNA prediction tools to screen the target gene of DEMs. In this study, 185 DEGs and three common DEMs were screened as well as 1266 target genes of the three DEMs were identified. Next, 16 overlapping dysregulated genes from 185 DEGs and 1266 target gene were obtained. Meanwhile, we constructed the miRNA gene regulatory network and determined miRNA‐508‐3p‐NR4A3 pair for deeper exploring. Experiment methods verified the functional expression of miR‐508‐3p in PAH and its signalling cascade. We observed that ectopic miR‐508‐3p expression promotes proliferation and migration of pulmonary artery smooth muscle cell (PASMC). Bioinformatic, dual‐luciferase assay showed NR4A3 represents directly targeted gene of miR‐508‐3p. Mechanistically, we demonstrated that down‐regulation of miR‐508‐3p advances PASMC proliferation and migration via inducing NR4A3 to activate MAPK/ERK kinase signalling pathway. Altogether, our research provides a promising diagnosis of predictor and therapeutic avenues for patients in PAH.

Probiogenomics of Lactobacillus delbrueckii subsp. lactis CIDCA 133: In Silico, In Vitro, and In Vivo Approaches

Lactobacillus delbrueckii subsp. lactis CIDCA 133 (CIDCA 133) has been reported as a potential probiotic strain, presenting immunomodulatory properties. This study investigated the possible genes and molecular mechanism involved with a probiotic profile of CIDCA 133 through a genomic approach associated with in vitro and in vivo analysis. Genomic analysis corroborates the species identification carried out by the classical microbiological method. Phenotypic assays demonstrated that the CIDCA 133 strain could survive acidic, osmotic, and thermic stresses. In addition, this strain shows antibacterial activity against Salmonella Typhimurium and presents immunostimulatory properties capable of upregulating anti-inflammatory cytokines Il10 and Tgfb1 gene expression through inhibition of Nfkb1 gene expression. These reported effects can be associated with secreted, membrane/exposed to the surface and cytoplasmic proteins, and bacteriocins-encoding genes predicted in silico. Furthermore, our results showed the genes and the possible mechanisms used by CIDCA 133 to produce their beneficial host effects and highlight its use as a probiotic microorganism.

Identification of novel biomarkers involved in pulmonary arterial hypertension based on multiple-microarray analysis

Pulmonary arterial hypertension (PAH) is a life-threatening chronic cardiopulmonary disorder. However, studies providing PAH-related gene expression profiles are scarce. To identify hub genes involved in PAH, we investigate two microarray data sets from gene expression omnibus (GEO). A total of 150 differentially expressed genes (DEGs) were identified by limma package. Enriched Gene Ontology (GO) annotations and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of DEGs mostly included mitotic nuclear division, ATPase activity, and Herpes simplex virus one infection. Ten hub genes from three significant modules were ascertained by Cytoscape (CytoHubba). Gene set enrichment analysis (GSEA) plots showed that transcription elongation factor complex was the most significantly enriched gene set positively correlated with the PAH group. At the same time, solute proton symporter activity was the most significantly enriched gene set positively correlated with the control group. Correlation analysis between hub genes suggested that SMC4, TOP2A, SMC2, KIF11, KIF23, ANLN, ARHGAP11A, SMC3, SMC6 and RAD50 may involve in the pathogenesis of PAH. Then, the miRNA-target genes regulation network was performed to unveil the underlying complex association among them. Finally, RNA extracted from monocrotaline (MCT)-induced Rat-PAH model lung artery tissues were to conduct quantitative real-time PCR (qRT-PCR) to validate these hub genes. In conclusion, our study offers new evidence for the underlying molecular mechanisms of PAH as well as attractive targets for diagnosis and treatment of PAH.

Mutation load at a mimicry supergene sheds new light on the evolution of inversion polymorphisms

Chromosomal inversions are ubiquitous in genomes and often coordinate complex phenotypes, such as the covariation of behavior and morphology in many birds, fishes, insects or mammals^1-11. However, why and how inversions become associated with polymorphic traits remains obscure. Here we show that despite a strong selective advantage when they form, inversions accumulate recessive deleterious mutations that generate frequency-dependent selection and promote their maintenance at intermediate frequency. Combining genomics and in vivo fitness analyses in a model butterfly for wing-pattern polymorphism, Heliconius numata, we reveal that three ecologically advantageous inversions have built up a heavy mutational load from the sequential accumulation of deleterious mutations and transposable elements. Inversions associate with sharply reduced viability when homozygous, which prevents them from replacing ancestral chromosome arrangements. Our results suggest that other complex polymorphisms, rather than representing adaptations to competing ecological optima, could evolve because chromosomal rearrangements are intrinsically prone to carrying recessive harmful mutations.

Biodiversity of magnetotactic bacteria in the tropical marine environment of Singapore revealed by metagenomic analysis

Most studies on the diversity of magnetotactic bacteria (MTB) have been conducted on samples obtained from the Northern or the Southern hemispheres. The diversity of MTB in tropical Asia near the geo-equator, with a close-to-zero geomagnetic inclination, weak magnetic field and constantly high seawater temperature has never been explored. This study aims to decipher the diversity of MTB in the marine environment of Singapore through shotgun metagenomics. Although MTB has been acknowledged to be ubiquitous in aquatic environments, we did not observe magnetotactic behaviour in the samples. However, we detected the presence and determined the diversity of MTB through bioinformatic analyses. Metagenomic analysis suggested majority of the MTB in the seafloor sediments represents novel MTB taxa that cannot be classified at the species level. The relative abundance of MTB (~0.2-1.69%) in the samples collected from the marine environment of Singapore was found to be substantially lower than studies for other regions. In contrast to other studies, the genera Magnetovibrio and Desulfamplus, but not Magnetococcus, were the dominant MTB. Additionally, we recovered 3 MTB genomic bins that are unclassified at the species level, with Magnetovibrio blakemorei being the closest-associated genome. All the recovered genomic bins contain homologs of at least 5 of the 7 mam genes but lack homologs for mamI, a membrane protein suggested to take part in the magenetosome invagination. This study fills in the knowledge gap of MTB biodiversity in the tropical marine environment near the geo-equator. Our findings will facilitate future research efforts aiming to unravel the ecological roles of MTB in the tropical marine environments as well as to bioprospecting novel MTB that have been adapted to tropical marine environments for biotechnological applications.

TOA: A software package for automated functional annotation in non-model plant species

The increase of sequencing capacity provided by high-throughput platforms has made it possible to routinely obtain large sets of genomic and transcriptomic sequences from model and non-model organisms. Subsequent genomic analysis and gene discovery in next-generation sequencing experiments are, however, bottlenecked by functional annotation. One common way to perform functional annotation of sets of sequences obtained from next-generation sequencing experiments, is by searching for homologous sequences and accessing the related functional information deposited in genomic databases. Functional annotation is especially challenging for non-model organisms, like many plant species. In such cases, existing free and commercial general-purpose applications may not offer complete and accurate results. We present TOA (Taxonomy-oriented annotation), a Python-based user-friendly open source application designed to establish functional annotation pipelines geared towards non-model plant species that can run in Linux/Mac computers, HPCs and cloud servers. TOA performs homology searches against proteins stored in the PLAZA databases, NCBI RefSeq Plant, Nucleotide Database and Non-Redundant Protein Sequence Database, and outputs functional information from several ontology systems: Gene Ontology, InterPro, EC, KEGG, Mapman and MetaCyc. The software performance was validated by comparing the runtimes, total number of annotated sequences and accuracy of the functional information obtained for several plant benchmark data sets with TOA and other functional annotation solutions. TOA outperformed the other software in terms of number of annotated sequences and accuracy of the annotation and constitutes a good alternative to improve functional annotation in plants. TOA is especially recommended for gymnosperms or for low quality sequence data sets of non-model plants.

Host-symbiont coevolution, cryptic structure, and bleaching susceptibility, in a coral species complex (Scleractinia; Poritidae)

The 'species' is a key concept for conservation and evolutionary biology, yet the lines between population and species-level variation are often blurred, especially for corals. The 'Porites lobata species complex' consists of branching and mounding corals that form reefs across the Pacific. We used reduced representation meta-genomic sequencing to examine genetic relationships within this species complex and to identify candidate loci associated with colony morphology, cryptic genetic structure, and apparent bleaching susceptibility. We compared existing Porites data with bleached and unbleached colonies of the branching coral P. compressa collected in Kāne'ohe Bay Hawai'i during the 2015 coral bleaching event. Loci that mapped to coral, symbiont, and microbial references revealed genetic structure consistent with recent host-symbiont co-evolution. Cryptic genetic clades were resolved that previous work has associated with distance from shore, but no genetic structure was associated with bleaching. We identified many candidate loci associated with morphospecies, including candidate host and symbiont loci with fixed differences between branching and mounding corals. We also found many loci associated with cryptic genetic structure, yet relatively few loci associated with bleaching. Recent host-symbiont co-evolution and rapid diversification suggests that variation and therefore the capacity of these corals to adapt may be underappreciated.

Plasma Proteome Responses in Salmonid Fish Following Immunization

Vaccination plays a critical role in the protection of humans and other animals from infectious diseases. However, the same vaccine often confers different protection levels among individuals due to variation in genetics and/or immunological histories. While this represents a well-recognized issue in humans, it has received little attention in fish. Here we address this knowledge gap in a proteomic study of rainbow trout (Oncorhynchus mykiss, Walbaum), using non-lethal repeated blood sampling to establish the plasma protein response of individual fish following immunization. Six trout were immunized with adjuvanted hen egg-white lysozyme (HEL) and peripheral blood sampled at ten time points from day 0 to day 84 post-injection. We confirm that an antigen-specific antibody response to HEL was raised, showing differences in timing and magnitude among individuals. Using label-free liquid chromatography-mass spectrometry, we quantified the abundance of 278 plasma proteins across the timecourse. As part of the analysis, we show that this approach can distinguish many (but not all) duplicated plasma proteins encoded by paralogous genes retained from the salmonid-specific whole genome duplication event. Global variation in the plasma proteome was predominantly explained by individual differences among fish. However, sampling day explained a major component of variation in abundance for a statistically defined subset of 41 proteins, representing 15% of those detected. These proteins clustered into five groups showing distinct temporal responses to HEL immunization at the population level, and include classical immune (e.g. complement system members) and acute phase molecules (e.g. apolipoproteins, haptoglobins), several enzymes and other proteins supporting the immune response, in addition to evolutionarily conserved molecules that are as yet uncharacterized. Overall, this study improves our understanding of the fish plasma proteome, provides valuable marker proteins for different phases of the immune response, and has implications for vaccine development and the design of immune challenge experiments.

Review on the Computational Genome Annotation of Sequences Obtained by Next-Generation Sequencing

Next-Generation Sequencing (NGS) has made it easier to obtain genome-wide sequence data and it has shifted the research focus into genome annotation. The challenging tasks involved in annotation rely on the currently available tools and techniques to decode the information contained in nucleotide sequences. This information will improve our understanding of general aspects of life and evolution and improve our ability to diagnose genetic disorders. Here, we present a summary of both structural and functional annotations, as well as the associated comparative annotation tools and pipelines. We highlight visualization tools that immensely aid the annotation process and the contributions of the scientific community to the annotation. Further, we discuss quality-control practices and the need for re-annotation, and highlight the future of annotation.

Keep Calm and Survive: Adaptation Strategies to Energy Crisis in Fruit Trees under Root Hypoxia

Plants are permanently facing challenges imposed by the environment which, in the context of the current scenario of global climate change, implies a constant process of adaptation to survive and even, in the case of crops, at least maintain yield. O₂ deficiency at the rhizosphere level, i.e., root hypoxia, is one of the factors with the greatest impact at whole-plant level. At cellular level, this O₂ deficiency provokes a disturbance in the energy metabolism which has notable consequences on the yield of plant crops. In this sense, although several physiological studies describe processes involved in plant adaptation to root hypoxia in woody fruit trees, with emphasis on the negative impacts on photosynthetic rate, there are very few studies that include -omics strategies for specifically understanding these processes in the roots of such species. Through a de novo assembly approach, a comparative transcriptome study of waterlogged Prunus spp. genotypes contrasting in their tolerance to root hypoxia was revisited in order to gain a deeper insight into the reconfiguration of pivotal pathways involved in energy metabolism. This re-analysis describes the classically altered pathways seen in the roots of woody fruit trees under hypoxia, but also routes that link them to pathways involved with nitrogen assimilation and the maintenance of cytoplasmic pH and glycolytic flow. In addition, the effects of root hypoxia on the transcription of genes related to the mitochondrial oxidative phosphorylation system, responsible for providing adenosine triphosphate (ATP) to the cell, are discussed in terms of their roles in the energy balance, reactive oxygen species (ROS) metabolism and aerenchyma formation. This review compiles key findings that help to explain the trait of tolerance to root hypoxia in woody fruit species, giving special attention to their strategies for managing the energy crisis. Finally, research challenges addressing less-explored topics in recovery and stress memory in woody fruit trees are pointed out.

In silico functional prediction of hypothetical proteins from the core genome of Corynebacterium pseudotuberculosis biovar ovis

Corynebacterium pseudotuberculosis is a pathogen of veterinary relevance diseases, being divided into two biovars: equi and ovis; causing ulcerative lymphangitis and caseous lymphadenitis, respectively. The isolation and sequencing of C. pseudotuberculosis biovar ovis strains in the Northern and Northeastern regions of Brazil exhibited the emergence of this pathogen, which causes economic losses to small ruminant producers, and condemnation of carcasses and skins of animals. Through the pan-genomic approach, it is possible to determine and analyze genes that are shared by all strains of a species—the core genome. However, many of these genes do not have any predicted function, being characterized as hypothetical proteins (HP). In this study, we considered 32 C. pseudotuberculosis biovar ovis genomes for the pan-genomic analysis, where were identified 172 HP present in a core genome composed by 1255 genes. We are able to functionally annotate 80 sequences previously characterized as HP through the identification of structural features as conserved domains and families. Furthermore, we analyzed the physicochemical properties, subcellular localization and molecular function. Additionally, through RNA-seq data, we investigated the differential gene expression of the annotated HP. Genes inserted in pathogenicity islands had their virulence potential evaluated. Also, we have analyzed the existence of functional associations for their products based on protein–protein interaction networks, and perform the structural prediction of three targets. Due to the integration of different strategies, this study can underlie deeper in vitro researches in the characterization of these HP and the search for new solutions for combat this pathogen.

Genomic Variation and Arsenic Tolerance Emerged as Niche Specific Adaptations by Different Exiguobacterium Strains Isolated From the Extreme Salar de Huasco Environment in Chilean - Altiplano

Polyextremophilic bacteria can thrive in environments with multiple stressors such as the Salar de Huasco (SH). Microbial communities in SH are exposed to low atmospheric pressure, high UV radiation, wide temperature ranges, salinity gradient and the presence of toxic compounds such as arsenic (As). In this work we focus on arsenic stress as one of the main adverse factors in SH and bacteria that belong to the Exiguobacterium genus due to their plasticity and ubiquity. Therefore, our aim was to shed light on the effect of niche conditions pressure (particularly arsenic), on the adaptation and divergence (at genotypic and phenotypic levels) of Exiguobacterium strains from five different SH sites. Also, to capture greater diversity in this genus, we use as outgroup five As(III) sensitive strains isolated from Easter Island (Chile) and The Great Salt Lake (United States). For this, samples were obtained from five different SH sites under an arsenic gradient (9 to 321 mg/kg: sediment) and isolated and sequenced the genomes of 14 Exiguobacterium strains, which had different arsenic tolerance levels. Then, we used comparative genomic analysis to assess the genomic divergence of these strains and their association with phenotypic differences such as arsenic tolerance levels and the ability to resist poly-stress. Phylogenetic analysis showed that SH strains share a common ancestor. Consequently, populations were separated and structured in different SH microenvironments, giving rise to multiple coexisting lineages. Hence, this genotypic variability is also evidenced by the COG (Clusters of Orthologous Groups) composition and the size of their accessory genomes. Interestingly, these observations correlate with physiological traits such as growth patterns, gene expression, and enzyme activity related to arsenic response and/or tolerance. Therefore, Exiguobacterium strains from SH are adapted to physiologically overcome the contrasting environmental conditions, like the arsenic present in their habitat.

Co-Expression Networks for Causal Gene Identification Based on RNA-Seq Data of Corynebacterium pseudotuberculosis

Corynebacterium pseudotuberculosis is a Gram-positive bacterium that causes caseous lymphadenitis, a disease that predominantly affects sheep, goat, cattle, buffalo, and horses, but has also been recognized in other animals. This bacterium generates a severe economic impact on countries producing meat. Gene expression studies using RNA-Seq are one of the most commonly used techniques to perform transcriptional experiments. Computational analysis of such data through reverse-engineering algorithms leads to a better understanding of the genome-wide complexity of gene interactomes, enabling the identification of genes having the most significant functions inferred by the activated stress response pathways. In this study, we identified the influential or causal genes from four RNA-Seq datasets from different stress conditions (high iron, low iron, acid, osmosis, and PH) in C. pseudotuberculosis, using a consensus-based network inference algorithm called miRsigand next identified the causal genes in the network using the miRinfluence tool, which is based on the influence diffusion model. We found that over 50% of the genes identified as influential had some essential cellular functions in the genomes. In the strains analyzed, most of the causal genes had crucial roles or participated in processes associated with the response to extracellular stresses, pathogenicity, membrane components, and essential genes. This research brings new insight into the understanding of virulence and infection by C. pseudotuberculosis.

Novel Target Exploration from Hypothetical Proteins of Klebsiella pneumoniae MGH 78578 Reveals a Protein Involved in Host-Pathogen Interaction

The opportunistic pathogen Klebsiella pneumoniae is a causative agent of several hospital-acquired infections. It has become resistant to a wide range of currently available antibiotics, leading to high mortality rates among patients; this has further led to a demand for novel therapeutic intervention to treat such infections. Using a series of in silico analyses, the present study aims to explore novel drug/vaccine candidates from the hypothetical proteins of K. pneumoniae. A total of 540 proteins were found to be hypothetical in this organism. Analysis of these 540 hypothetical proteins revealed 30 pathogen-specific proteins essential for pathogen survival. A motifs/domain family analysis, similarity search against known proteins, gene ontology, and protein–protein interaction analysis of the shortlisted 30 proteins led to functional assignment for 17 proteins. They were mainly cataloged as enzymes, lipoproteins, stress-induced proteins, transporters, and other proteins (viz., two-component proteins, skeletal proteins and toxins). Among the annotated proteins, 16 proteins, located in the cytoplasm, periplasm, and inner membrane, were considered as potential drug targets, and one extracellular protein was considered as a vaccine candidate. A druggability analysis indicated that the identified 17 drug/vaccine candidates were “novel”. Furthermore, a host–pathogen interaction analysis of these identified target candidates revealed a betaine/carnitine/choline transporters (BCCT) family protein showing interactions with five host proteins. Structure prediction and validation were carried out for this protein, which could aid in structure-based inhibitor design.

Toward a chimeric vaccine against multiple isolates of Mycobacteroides - An integrative approach

AIM

Nontuberculous mycobacterial (NTM) infection such as endophthalmitis, dacryocystitis, and canaliculitis are pervasive across the globe and are currently managed by antibiotics. However, the recent cases of Mycobacteroides developing drug resistance reported along with the improper practice of medicine intrigued us to explore its genomic and proteomic canvas at a global scale and develop a chimeric vaccine against Mycobacteroides.

MAIN METHODS

We carried out a vivid genomic study on five recently sequenced strains of Mycobacteroides and explored their Pan-core genome/proteome in three different phases. The promiscuous antigenic proteins were identified via a subtractive proteomics approach that qualified for virulence causation, resistance and essentiality factors for this notorious bacterium. An integrated pipeline was developed for the identification of B-Cell, MHC (Major histocompatibility complex) class I and II epitopes.

KEY FINDINGS

Phase I identified the shreds of evidence of reductive evolution and propensity of the Pan-genome of Mycobacteroides getting closed soon. Phase II and Phase III produced 8 vaccine constructs. Our final vaccine construct, V6 qualified for all tests such as absence for allergenicity, presence of antigenicity, etc. V6 contains β-defensin as an adjuvant, linkers, Lysosomal-associated membrane protein 1 (LAMP1) signal peptide, and PADRE (Pan HLA-DR epitopes) amino acid sequence. Besides, V6 also interacts with a maximum number of MHC molecules and the TLR4/MD2 (Toll-like receptor 4/Myeloid differentiation factor 2) complex confirmed by docking and molecular dynamics simulation studies.

SIGNIFICANCE

The knowledge harnessed from the current study can help improve the current treatment regimens or in an event of an outbreak and propel further related studies.

Gene expression variation in the brains of harvester ant foragers is associated with collective behavior

Natural selection on collective behavior acts on variation among colonies in behavior that is associated with reproductive success. In the red harvester ant (Pogonomyrmex barbatus), variation among colonies in the collective regulation of foraging in response to humidity is associated with colony reproductive success. We used RNA-seq to examine gene expression in the brains of foragers in a natural setting. We find that colonies differ in the expression of neurophysiologically-relevant genes in forager brains, and a fraction of these gene expression differences are associated with two colony traits: sensitivity of foraging activity to humidity, and forager brain dopamine to serotonin ratio. Loci that were correlated with colony behavioral differences were enriched in neurotransmitter receptor signaling & metabolic functions, tended to be more central to coexpression networks, and are evolving under higher protein-coding sequence constraint. Natural selection may shape colony foraging behavior through variation in gene expression.

Grafting alters tomato transcriptome and enhances tolerance to an airborne virus infection

Grafting of commercial tomato varieties and hybrids on the tomato ecotype Manduria resulted in high levels of tolerance to the infection of Sw5 resistance-breaking strains of tomato spotted wilt virus and of severe cucumber mosaic virus strains supporting hypervirulent satellite RNAs that co-determine stunting and necrotic phenotypes in tomato. To decipher the basis of such tolerance, here we used a RNAseq analysis to study the transcriptome profiles of the Manduria ecotype and of the susceptible variety UC82, and of their graft combinations, exposed or not to infection of the potato virus Y recombinant strain PVY^C-to. The analysis identified graft- and virus-responsive mRNAs differentially expressed in UC82 and Manduria, which led to an overall suitable level of tolerance to viral infection confirmed by the appearance of a recovery phenotype in Manduria and in all graft combinations. The transcriptome analysis suggested that graft wounding and viral infection had diverging effects on tomato transcriptome and that the Manduria ecotype was less responsive than the UC82 to both graft wounding and potyviral infection. We propose that the differential response to the two types of stress could account for the tolerance to viral infection observed in the Manduria ecotype as well as in the susceptible tomato variety UC82 self-grafted or grafted on the Manduria ecotype.