The Clostridium-infecting filamentous phage CAK1 genome analysis allows to define a new potential clade of Tubulavirales

What we know about Tubulavirales, i.e. filamentous phages, essentially comes from Gram-negative-infecting Inoviridae. However, metagenomics recently suggests filamentous phages are much more widespread and diverse. Here, we report the complete sequence and functional annotation of CAK1, a 6.6 kb filamentous phage that was shown to chronically infect Clostridium beijerinckii 30 years ago and only represents the second filamentous phage cultivated on a Gram-positive bacterium. CAK1 has a typical filamentous phage modular genome with no homologs in databases and we were interested to compare it with a pig gut filamentous phage metagenomics dataset that we previously assembled and for which many filamentous phages were predicted to infect Clostridium species by bioinformatics means. CAK1 is distantly related to nine of these sequences, two of which have been predicted as Clostridium-associated. In itself, this small cluster of CAK1-connected sequences sheds light on the diversity of filamentous phages that putatively infect Clostridium species, and probably many other Gram-positive genera.

Transcriptome profiling disclosed the effect of single and combined drought and heat stress on reprogramming of genes expression in barley flag leaf

Despite numerous studies aimed at unraveling the genetic background of barley's response to abiotic stress, the modulation of the transcriptome induced by combinatorial drought and increased temperature remains largely unrecognized. Very limited studies were done, especially on the flag leaf, which plays an important role in grain filling in cereals. In the present study, transcriptome profiles, along with chlorophyll fluorescence parameters and yield components, were compared between barley genotypes with different flag leaf sizes under single and combined drought and heat stress. High-throughput mRNA sequencing revealed 2,457 differentially expressed genes, which were functionally interpreted using Gene Ontology term enrichment analysis. The transcriptomic signature under double stress was more similar to effects caused by drought than by elevated temperature; it was also manifested at phenotypic and chlorophyll fluorescence levels. Both common and stress-specific changes in transcript abundance were identified. Genes regulated commonly across stress treatments, determining universal stress responses, were associated, among others, with responses to drought, heat, and oxidative stress. In addition, changes specific to the size of the flag leaf blade were found. Our study allowed us to identify sets of genes assigned to various processes underlying the response to drought and heat, including photosynthesis, the abscisic acid pathway, and lipid transport. Genes encoding LEA proteins, including dehydrins and heat shock proteins, were especially induced by stress treatments. Some association between genetic composition and flag leaf size was confirmed. However, there was no general coincidence between SNP polymorphism of genotypes and differential expression of genes induced by stress factors. This research provided novel insight into the molecular mechanisms of barley flag leaf that determine drought and heat response, as well as their co-occurrence.

Comparison between the Gametophyte and the Sporophyte Transcriptomes of the Endangered Fern Vandenboschia speciosa

Genomic resources are essential to understanding the evolution and functional biology of organisms. Nevertheless, generating genomic resources from endangered species may be challenging due to the scarcity of available specimens and sampling difficulties. In this study, we compare the transcriptomes of the sporophyte and the gametophyte of the endangered fern Vandenboschia speciosa. After Illumina sequencing and de novo transcriptome assembly of the gametophyte, annotation proved the existence of cross-species contamination in the gametophyte sample. Thus, we developed an in silico decontamination step for the gametophyte sequences. Once the quality check of the decontaminated reads passed, we produced a de novo assembly with the decontaminated gametophyte reads (with 43,139 contigs) and another combining the sporophyte and in silico decontaminated gametophyte reads (with 42,918 contigs). A comparison of the enriched GO terms from the top 1000 most expressed transcripts from both tissues showed that the gametophyte GO term set was enriched in sequences involved in development, response to stress, and plastid organization, while the sporophyte GO term set had a larger representation of more general metabolic functions. This study complements the available genomic resources on the life cycle of the endangered fern Vandenboschia speciosa.

Liver-Specific Polygenic Risk Score Is Associated with Alzheimer's Disease Diagnosis

BACKGROUND

Our understanding of the pathophysiology underlying Alzheimer's disease (AD) has benefited from genomic analyses, including those that leverage polygenic risk score (PRS) models of disease. The use of functional annotation has been able to improve the power of genomic models.

OBJECTIVE

We sought to leverage genomic functional annotations to build tissue-specific AD PRS models and study their relationship with AD and its biomarkers.

METHODS

We built 13 tissue-specific AD PRS and studied the scores' relationships with AD diagnosis, cerebrospinal fluid (CSF) amyloid, CSF tau, and other CSF biomarkers in two longitudinal cohort studies of AD.

RESULTS

The AD PRS model that was most predictive of AD diagnosis (even without APOE) was the liver AD PRS: n = 1,115; odds ratio = 2.15 (1.67-2.78), p = 3.62×10-9. The liver AD PRS was also statistically significantly associated with cerebrospinal fluid biomarker evidence of amyloid-β (Aβ42:Aβ40 ratio, p = 3.53×10-6) and the phosphorylated tau:amyloid-β ratio (p = 1.45×10-5).

CONCLUSION

These findings provide further evidence of the role of the liver-functional genome in AD and the benefits of incorporating functional annotation into genomic research.

Molecular characteristics and transcriptional regulatory of spermatogenesis-related gene RFX2 in adult Banna mini-pig inbred line (BMI)

RFX2 plays critical roles in mammalian spermatogenesis and cilium maturation. Here, the testes of 12-month-old adult boars of Banna mini-pig inbred line (BMI) were subjected to whole-transcriptome sequencing. The results indicated that the average expression (raw count) of RFX2 gene in BMI testes was 16138.25, and the average expression value of the corresponding transcript ENSSSCT00000043271.2 was 123.1898. The CDS of RFX2 obtained from BMI testes was 2,817 bp (GenBank accession number: OL362242). Gene structure analysis showed that RFX2 was located on chromosome 2 of the pig genome with 19 exons. Protein structure analysis indicated that RFX2 contains 728 amino acids with two conserved domains. Phylogenetic analysis revealed that RFX2 was highly conserved with evolutionary homologies among mammalian species. Other analyses, including PPI networks, KEGG, and GO, indicated that BMI RFX2 had interactions with 43 proteins involving various functions, such as in cell cycle, spermatid development, spermatid differentiation, cilium assembly, and cilium organization, etc. Correlation analysis between these proteins and the transcriptome data implied that RFX2 was significantly associated with FOXJ1, DNAH9, TMEM138, E2F7, and ATR, and particularly showed the highest correlation with ATR, demonstrating the importance of RFX2 and ART in spermatogenesis. Functional annotation implied that RFX2 was involved in 17 GO terms, including three cellular components (CC), six molecular functions (MF), and eight biological processes (BP). The analysis of miRNA-gene targeting indicated that BMI RFX2 was mainly regulated by two miRNAs, among which four lncRNAs and five lncRNAs competitively bound ssc-miR-365-5p and ssc-miR-744 with RFX2, respectively. Further, the dual-luciferase report assay indicated that the ssc-miR-365-5p and ssc-miR-744 significantly reduced luciferase activity of RFX2 3'UTR in the 293T cells, suggesting that these two miRNAs regulated the expression of RFX2. Our results revealed the important role of RFX2 in BMI spermatogenesis, making it an intriguing candidate for follow-up studies.

State-of-the-art Application of Artificial Intelligence to Transporter-centered Functional and Pharmaceutical Research

Protein transporters not only have essential functions in regulating the transport of endogenous substrates and remote communication between organs and organisms, but they also play a vital role in drug absorption, distribution, and excretion and are recognized as major determinants of drug safety and efficacy. Understanding transporter function is important for drug development and clarifying disease mechanisms. However, the experimental-based functional research on transporters has been challenged and hinged by the expensive cost of time and resources. With the increasing volume of relevant omics datasets and the rapid evolution of artificial intelligence (AI) techniques, next-generation AI is becoming increasingly prevalent in the functional and pharmaceutical research of transporters. Thus, a comprehensive discussion on the state-of-the-art application of AI in three cutting-edge directions was provided in this review, which included (a) transporter classification and function annotation, (b) structure discovery of membrane transporters, and (c) drug-transporter interaction prediction. This study provides a panoramic view of AI algorithms and tools applied to the field of transporters. It is expected to guide a better understanding and utilization of AI techniques for in-depth studies of transporter-centered functional and pharmaceutical research.

Identification and functional annotation of hypothetical proteins of uropathogenic Escherichia coli strain CFT073 towards designing antimicrobial drug targets

Urinary tract infections are a serious health concern worldwide, especially in developing countries. Escherichia coli strain CFT073 is a highly virulent pathogenic bacterial strain. CFT073 proteome contains 4897 proteins, out of which 992 have been classified as hypothetical proteins. Identification and characterization of hypothetical proteins can aid in the selection of targets for drug design. In this study, we studied the hypothetical proteins from the UPEC strain CFT073 using various computational tools. By NCBI-CDD, 376 protein sequences showed conserved domains. Based on the functional motifs in their primary sequences, we classified these 376 hypothetical proteins into 7 functional categories. Further KEGG database was used to find the roles of these hypothetical proteins in several pathways. Protein interaction network analysis of hypothetical proteins identified 53 proteins as highly interacting metabolic proteins. Virulence factor analysis of the proteins identified 8 proteins as virulent. We conducted a non-homology search for the identified proteins of UPEC in the available human proteome. We observed that 35 proteins are non-homologous to humans and hence could be selected for drug designing targets. Qualitative characterization of the selected 35 non-homologous hypothetical proteins including essentiality analysis and evaluation of druggability by similarity search against drug bank database was performed. Out of these 35 proteins, three-dimensional structures of six proteins (NP_752562.1, NP_756345.1, NP_754893.1, NP_756600.2, NP_755264.1 and NP_752994.1) could be successfully modelled. These new annotations can help to better understand disease mechanisms at the molecular level, as well as provide new targets for drug development against the UPEC strain CFT073.Communicated by Ramaswamy H. Sarma.

Antisense-oligonucleotide-mediated perturbation of long non-coding RNA reveals functional features in stem cells and across cell types

Within the scope of the FANTOM6 consortium, we perform a large-scale knockdown of 200 long non-coding RNAs (lncRNAs) in human induced pluripotent stem cells (iPSCs) and systematically characterize their roles in self-renewal and pluripotency. We find 36 lncRNAs (18%) exhibiting cell growth inhibition. From the knockdown of 123 lncRNAs with transcriptome profiling, 36 lncRNAs (29.3%) show molecular phenotypes. Integrating the molecular phenotypes with chromatin-interaction assays further reveals cis- and trans-interacting partners as potential primary targets. Additionally, cell-type enrichment analysis identifies lncRNAs associated with pluripotency, while the knockdown of LINC02595, CATG00000090305.1, and RP11-148B6.2 modulates colony formation of iPSCs. We compare our results with previously published fibroblasts phenotyping data and find that 2.9% of the lncRNAs exhibit a consistent cell growth phenotype, whereas we observe 58.3% agreement in molecular phenotypes. This highlights that molecular phenotyping is more comprehensive in revealing affected pathways.

Functional and Taxonomic Traits of the Gut Microbiota in Type 1 Diabetes Children at the Onset: A Metaproteomic Study

Type 1 diabetes (T1D) is a chronic autoimmune metabolic disorder with onset in pediatric/adolescent age, characterized by insufficient insulin production, due to a progressive destruction of pancreatic β-cells. Evidence on the correlation between the human gut microbiota (GM) composition and T1D insurgence has been recently reported. In particular, 16S rRNA-based metagenomics has been intensively employed in the last decade in a number of investigations focused on GM representation in relation to a pre-disease state or to a response to clinical treatments. On the other hand, few works have been published using alternative functional omics, which is more suitable to provide a different interpretation of such a relationship. In this work, we pursued a comprehensive metaproteomic investigation on T1D children compared with a group of siblings (SIBL) and a reference control group (CTRL) composed of aged matched healthy subjects, with the aim of finding features in the T1D patients' GM to be related with the onset of the disease. Modulated metaproteins were found either by comparing T1D with CTRL and SIBL or by stratifying T1D by insulin need (IN), as a proxy of β-cells damage, showing some functional and taxonomic traits of the GM, possibly related to the disease onset at different stages of severity.

Genome-level analyses resolve an ancient lineage of symbiotic ascomycetes

Ascomycota account for about two-thirds of named fungal species.1 Over 98% of known Ascomycota belong to the Pezizomycotina, including many economically important species as well as diverse pathogens, decomposers, and mutualistic symbionts.2 Our understanding of Pezizomycotina evolution has until now been based on sampling traditionally well-defined taxonomic classes.3,4,5 However, considerable diversity exists in undersampled and uncultured, putatively early-diverging lineages, and the effect of these on evolutionary models has seldom been tested. We obtained genomes from 30 putative early-diverging lineages not included in recent phylogenomic analyses and analyzed these together with 451 genomes covering all available ascomycete genera. We show that 22 of these lineages, collectively representing over 600 species, trace back to a single origin that diverged from the common ancestor of Eurotiomycetes and Lecanoromycetes over 300 million years BP. The new clade, which we recognize as a more broadly defined Lichinomycetes, includes lichen and insect symbionts, endophytes, and putative mycorrhizae and encompasses a range of morphologies so disparate that they have recently been placed in six different taxonomic classes. To test for shared hidden features within this group, we analyzed genome content and compared gene repertoires to related groups in Ascomycota. Regardless of their lifestyle, Lichinomycetes have smaller genomes than most filamentous Ascomycota, with reduced arsenals of carbohydrate-degrading enzymes and secondary metabolite gene clusters. Our expanded genome sample resolves the relationships of numerous "orphan" ascomycetes and establishes the independent evolutionary origins of multiple mutualistic lifestyles within a single, morphologically hyperdiverse clade of fungi.

Functional Characterization, Mechanism, and Mode of Action of Putative Streptomycin Adenylyltransferase from Serratia marcescens

Nosocomial infections are serious threats to the entire world in healthcare settings. The major causative agents of nosocomial infections are bacterial pathogens, among which Enterobacteriaceae family member Serratia marcescens plays a crucial role. It is a gram-negative opportunistic pathogen, predominantly affecting patients in intensive-care units. The presence of intrinsic genes in S. marcescens led to the development of resistance to antibiotics for survival. Complete scanning of the proteome, including hypothetical and partially annotated proteins, paves the way for a better understanding of potential drug targets. The targeted protein expressed in E. coli BL21 (DE3) pLysS cells has shown complete resistance to aminoglycoside antibiotic streptomycin (>256 MCG). The recombinant protein was purified using affinity and size-exclusion chromatography and characterized using SDS-PAGE, western blotting, and MALDI-TOF analysis. Free phosphate bound to malachite green was detected at 620 nm, evident of the conversion of adenosine triphosphate to adenosine monophosphate during the adenylation process. Similarly, in the chromatographic assay, adenylated streptomycin absorbed at 260 nm in AKTA (FPLC), confirming the enzyme-catalyzed adenylation of streptomycin. Further, the adenylated product of streptomycin was confirmed through HPLC and mass spectrometry analysis. In conclusion, our characterization studies identified the partially annotated hypothetical protein as streptomycin adenylyltransferase.

[Single-cell transcriptome analysis of multigrade glioma heterogeneity and immune microenvironment revealed potential prognostic biomarkers]

Glioma, the most common intrinsic tumor of the central nervous system, is characterized by its high incidence and poor prognosis. The aim of this study was to identify differentially expressed genes (DEGs) between glioblastoma multiforme (GBM) and low-grade glioma (LGG) to explore prognostic factors of different grades of gliomas. Single-cell transcriptome sequencing data of gliomas were collected from the NCBI Gene Expression Omnibus (GEO), which included a total of 29 097 cell samples from three datasets. For the analysis of human gliomas of different grades, 21 071 cells were obtained by filtering, and 70 genes were screened from differentially expressed genes by gene ontology (GO) analysis, Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis, from which the gene DLL3 was focused by reviewing the literature. The TCGA-based gene expression profiling interactive analysis (GEPIA) database was used to explore the survival curves of genes in LGG and GBM, and the gene expression profiling interactive analysis and tumor immune estimation resource (TIMER) database was used to study the expression of key genes in gliomas of different grades, predicting biomarkers that were closely related to immunotherapy. The cBioPortal database was used to explore the relationship between DLL3 expression and 25 immune checkpoints. Gene set enrichment analysis (GSEA) further identified pathways associated with central genes. Finally, the efficacy of biomarkers in prognosis and prediction was validated in the Chinese glioma genome atlas (CGGA). These results demonstrated that prognostic genes are associated with tumor proliferation and progression. Analysis of biological information and survival suggested that these genes might serve as a promising prognostic biomarker and as new targets for selecting therapeutic strategies.

Transcriptome and Proteome Co-Profiling Offers an Understanding of Pre-Harvest Sprouting (PHS) Molecular Mechanisms in Wheat (Triticum aestivum)

While wheat (Triticum aestivum L.) is a widely grown and enjoyed crop, the diverse and complex global situation and climate are exacerbating the instability of its supply. In particular, pre-harvest sprouting (PHS) is one of the major abiotic stresses that frequently occurs due to irregular climate conditions, causing serious damage to wheat and its quality. In this study, transcriptomic analysis with RNA-seq and proteomic analysis with LC-MS/MS were performed in PHS-treated spikes from two wheat cultivars presenting PHS sensitivity and tolerance, respectively. A total of 13,154 differentially expressed genes (DEGs) and 706 differentially expressed proteins (DEPs) were identified in four comparison groups between the susceptible/tolerant cultivars. Gene function and correlation analysis were performed to determine the co-profiled genes and proteins affected by PHS treatment. In the functional annotation of each comparative group, similar functions were confirmed in each cultivar under PHS treatment; however, in Keumgang PHS+7 (K7) vs. Woori PHS+7 (W7), functional annotations presented clear differences in the "spliceosome" and "proteasome" pathways. In addition, our results indicate that alternative splicing and ubiquitin-proteasome support the regulation of germination and seed dormancy. This study provides an advanced understanding of the functions involved in transcription and translation related to PHS mechanisms, thus enabling specific proposals for the further analysis of germination and seed dormancy mechanisms and pathways in wheat.

The Fluorescent Enzyme Cascade Detects Low Abundance Protein Modifications Suitable for the Assembly of Functionally Annotated Modificatome Databases

Pathophysiological functions of proteins critically depend on both their chemical composition, including post-translational modifications, and their three-dimensional structure, commonly referred to as structure-activity relationship. Current analytical methods, like capillary electrophoresis or mass spectrometry, suffer from limitations, such as the detection of unexpected modifications at low abundance and their insensitivity to conformational changes. Building on previous enzyme-based analytical methods, we here introduce a fluorescence-based enzyme cascade (fEC), which can detect diverse chemical and conformational variations in protein samples and assemble them into digital databases. Together with complementary analytical methods an automated fEC analysis established unique modification-function relationships, which can be expanded to a proteome-wide scale, i. e. a functionally annotated modificatome. The fEC offers diverse applications, including hypersensitive biomarker detection in complex samples.

Forty Years without Family: Three Novel Bacteriophages with High Similarity to SPP1 Reveal Decades of Evolutionary Stasis since the Isolation of Their Famous Relative

SPP1, an extensively studied bacteriophage of the Gram-positive Bacillus subtilis, is a model system for the study of phage-host interactions. Despite progress in the isolation and characterization of Bacillus phages, no previously fully sequenced phages have shared more than passing genetic similarity to SPP1. Here, we describe three virulent phages very similar to SPP1; SPP1 has greater than 80% nucleotide sequence identity and shares more that 85% of its protein coding genes with these phages. This is remarkable, given more than 40 years between the isolation of SPP1 and these phages. All three phages have somewhat larger genomes and more genes than SPP1. We identified a new putative gene in SPP1 based on a conserved sequence found in all phages. Gene conservation connotes purifying selection and is observed in structural genes and genes involved in DNA metabolism, but also in genes of unknown function, suggesting an important role in phage survival independent of the environment. Patterns of divergence point to genes or gene domains likely involved in adaptation to diverse hosts or different environments. Ultimately, comparative genomics of related phages provides insight into the long-term selective pressures that affect phage-bacteria interactions and alter phage genome content.

Comparative Genomics of Six Lytic Bacillus subtilis Phages from the Southwest United States

Background

Despite their importance to microbial dynamics involving Bacillus subtilis, we have a limited understanding of the diversity of phages that can lyse this model organism.

Materials and Methods

Phages were isolated from soil samples collected from various sites in the southwest U.S. deserts on a wild B. subtilis strain. Their genomes were assembled, characterized, and bioinformatically compared.

Results

Six Siphoviruses with high nucleotide and amino acid similarity to each other (>80%) but very limited similarity to phages currently in GenBank were isolated. These phages have double-stranded DNA genomes (55,312 to 56,127 bp) with 86-91 putative protein coding genes, and a low GC content. Comparative genomics reveal differences in loci encoding proteins that are putatively involved in bacterial adsorption with evidence for genomic mosaicism and a possible role for small genes.

Conclusions

A comparative approach provides insights into phage evolution, including the role of indels in protein folding.

Short-Term Test for Toxicogenomic Analysis of Ecotoxic Modes of Action in Lemna minor

In the environmental risk assessment of substances, toxicity to aquatic plants is evaluated using, among other methods, the 7 dayLemna sp. growth inhibition test following the OECD TG 221. So far, the test is not applicable for short-term screening of toxicity, nor does it allow evaluation of toxic modes of action (MoA). The latter is also complicated by the lack of knowledge of gene functions in the test species. Using ecotoxicogenomics, we developed a time-shortened 3 day assay inLemna minor which allows discrimination of ecotoxic MoA. By examining the changes in gene expression induced by low effect concentrations of the pharmaceutical atorvastatin and the herbicide bentazon at the transcriptome and proteome levels, we were able to identify candidate biomarkers for the respective MoA. We developed a homology-based functional annotation pipeline for the reference genome ofL. minor, which allowed overrepresentation analysis of the gene ontologies affected by both test compounds. Genes affected by atorvastatin mainly influenced lipid synthesis and metabolism, whereas the bentazon-responsive genes were mainly involved in light response. Our approach is therefore less time-consuming but sensitive and allows assessment of MoA in L. minor. Using this shortened assay, investigation of expression changes of the identified candidate biomarkers may allow the development of MoA-specific screening approaches in the future.

Functional Annotation of Hypothetical Proteins From the Enterobacter cloacae B13 Strain and Its Association With Pathogenicity

Enterobacter cloacae B13 strain is a rod-shaped gram-negative bacterium that belongs to the Enterobacteriaceae family. It can cause respiratory and urinary tract infections, and is responsible for several outbreaks in hospitals. E. cloacae has become an important pathogen and an emerging global threat because of its opportunistic and multidrug resistant ability. However, little knowledge is present about a large portion of its proteins and functions. Therefore, functional annotation of the hypothetical proteins (HPs) can provide an improved understanding of this organism and its virulence activity. The workflow in the study included several bioinformatic tools which were utilized to characterize functions, family and domains, subcellular localization, physiochemical properties, and protein-protein interactions. The E. cloacae B13 strain has overall 604 HPs, among which 78 were functionally annotated with high confidence. Several proteins were identified as enzymes, regulatory, binding, and transmembrane proteins with essential functions. Furthermore, 23 HPs were predicted to be virulent factors. These virulent proteins are linked to pathogenesis with their contribution to biofilm formation, quorum sensing, 2-component signal transduction or secretion. Better knowledge about the HPs’ characteristics and functions will provide a greater overview of the proteome. Moreover, it will help against E. cloacae in neonatal intensive care unit (NICU) outbreaks and nosocomial infections.

Global Identification of White Lupin lncRNAs Reveals Their Role in Cluster Roots under Phosphorus Deficiency

Phosphorus (P) deficiency heterogeneously affected plant nutritional status and physiological performance, ultimately leading to a severe yield reduction. A few putative long non-coding RNAs (lncRNAs) responding to P-starvation in the model crops Arabidopsis thaliana and Oryza sativa have been characterized. White lupin (Lupinus albus) is of prime importance, and is a legume with increasing agronomic value as a protein crop as it exhibits extreme tolerance to nutrient deficiency, particularly P deficiency. Despite its adapted nature to P deficiency, nothing is known about low P-induced lncRNAs in white lupin roots. To address this issue, we identified 39,840 mRNA and 2028 lncRNAs in the eight developmental stages of white lupin root (S0–S7 and lateral root, LR) grown under P deficiency. From these 2028 lncRNAs, 1564 were intergenic and 464 natural antisense intergenic transcript (NAT) lncRNAs. We further predicted six potential targets of miRNAs with twelve lncRNAs, which may regulate P-deficiency-related processes. Moreover, the weighted gene co-expression network analysis (WGCNA) revealed seven modules that were correlated with the expression pattern of lncRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed 606 GO terms and 27 different pathways including signal transduction, energy synthesis, detoxification, and Pi transport. In addition, we screened 13 putative lncRNAs that showed a distinct expression pattern in each root, indicating their role in the P deficiency regulatory network. Therefore, white lupin may be a reference legume to characterize P-deficiency-responsive novel lncRNAs, which would highlight the role of lncRNAs in the regulation of plant responses to P deficiency.

De novo genome assembly of Bradysia cellarum (Diptera: Sciaridae), a notorious pest in traditional special vegetables in China

Bradysia cellarum (Diptera: Sciaridae) is a destructive vegetable insect pest infesting more than 30 species of host plants from seven families in Asia and Europe. B. cellarum causes grave problems in Chinese chive, which originated in China and is cultivated widely in East Asia. The B. cellarum infestation results in economic losses and subsequent severe food safety problems in farm productions, insecticide resistance and environmental pollution. The genomic and molecular information of B. cellarum to delineate the biological features, insecticide resistance, evolution remains poorly understood. Herein, we decode the whole genome of B. cellarum to delineate the underlying molecular mechanisms causing insecticide resistance. We constructed a highly reliable genome for B. cellarum using PacBio, Illumina and 10X Genomics sequencing platforms. The genome size of B. cellarum was 375.91 Mb with a contig N50 of 1.57 Mb. A total of 16,231 genes were identified, among which 93.8% were functionally annotated, and 42.06% were repeat sequences. According to phylogenetic analysis, B. cellarum diverged from the common ancestor of Drosophila melanogaster and Musca domestica ~139.3-191.0 million years ago. Moreover, some important genes responsible for significant insecticide resistance, such as cytochrome P450s, ABC transporters and those involved in glutathione metabolism, were expanded in B. cellarum. We assembled a high-quality B. cellarum genome to provide valuable insights into their life history strategies, insecticide resistance and biological behaviours. It also lays the foundation for exploring gene structure and functional evolution, as well as comparative genomics of B. cellarum and other model insect species.

InterCellDB: A User-Defined Database for Inferring Intercellular Networks

Recent advances in single cell RNA sequencing (scRNA-seq) empower insights into cell-cell crosstalk within specific tissues. However, customizable data analysis tools that decipher intercellular communication from gene expression in association with biological functions are lacking. The authors have developed InterCellDB, a platform that allows a user-defined analysis of intercellular communication using scRNA-seq datasets in combination with protein annotation information, including cellular localization and functional classification, and protein interaction properties. The application of InterCellDB in tumor microenvironment research is exemplified using two independent scRNA-seq datasets from human and mouse and it is demonstrated that InterCellDB-inferred cell-cell interactions and ligand-receptor pairs are experimentally valid.

RNA m6A Modification Changes in Postmortem Nucleus Accumbens of Subjects with Alcohol Use Disorder: A Pilot Study

BACKGROUND

The nucleus accumbens (NAc) is a key brain structure mediating the rewarding effect of alcohol and drug abuse. Chronic alcohol consumption may alter RNA methylome (or epitranscriptome) in the NAc, leading to altered gene expression and thus behavioral neuroadaptation to alcohol.

METHODS

This pilot study profiled the epitranscriptomes of mRNAs, long noncoding RNAs (lncRNAs), and microRNAs (miRNAs) in postmortem NAc of three male Caucasian subjects with alcohol use disorder (AUD) and three matched male Caucasian control subjects using Arraystar's m6A-mRNA&lncRNA Epitranscriptomic Microarray assay. Differentially methylated (DM) RNAs and the function of DM RNAs were analyzed by biostatistics and bioinformatics programs.

RESULTS

26 mRNAs were hypermethylated and three mRNAs were hypomethylated in the NAc of AUD subjects (≥2-fold changes and p ≤ 0.05). Most of these 29 DM mRNAs are involved in immune-related pathways (e.g., IL-17 signaling). Moreover, four lncRNAs were hypermethylated and one lncRNA was hypomethylated in the NAc of AUD subjects (≥2-fold changes and p ≤ 0.05). Additionally, three miRNAs were hypermethylated in the NAc of AUD subjects (≥2-fold changes and p ≤ 0.05).

CONCLUSIONS

This study revealed RNA methylomic changes in the NAc of AUD subjects, suggesting that chronic alcohol consumption may lead to AUD through epitranscriptomic RNA modifications. Our findings need to be replicated in a larger sample.

Extracting Information from Gene Coexpression Networks of Rhizobium leguminosarum

Nitrogen uptake in legumes is facilitated by bacteria such as Rhizobium leguminosarum. For this bacterium, gene expression data are available, but functional gene annotation is less well developed than for other model organisms. More annotations could lead to a better understanding of the pathways for growth, plant colonization, and nitrogen fixation in R. leguminosarum. In this study, we present a pipeline that combines novel scores from gene coexpression network analysis in a principled way to identify the genes that are associated with certain growth conditions or highly coexpressed with a predefined set of genes of interest. This association may lead to putative functional annotation or to a prioritized list of genes for further study.

MiRLog and dbmiR: prioritization and functional annotation tools to study human microRNA sequence variants

The recent identification of noncoding variants with pathogenic effects suggests that these variations could underlie a significant number of undiagnosed cases. Several computational methods have been developed to predict the functional impact of noncoding variants, but they exhibit only partial concordance and are not integrated with functional annotation resources, making the interpretation of these variants still challenging. MicroRNAs (miRNAs) are small noncoding RNA molecules that act as fine regulators of gene expression and play crucial functions in several biological processes, such as cell proliferation and differentiation. An increasing number of studies demonstrate a significant impact of miRNA single nucleotide variants (SNVs) both in Mendelian diseases and complex traits. To predict the functional effect of miRNA SNVs, we implemented a new meta‐predictor, MiRLog, and we integrated it into a comprehensive database, dbmiR, which includes a precompiled list of all possible miRNA allelic SNVs, providing their biological annotations at nucleotide and miRNA levels. MiRLog and dbmiR were used to explore the genetic variability of miRNAs in 15,708 human genomes included in the gnomAD project, finding several ultra‐rare SNVs with a potentially deleterious effect on miRNA biogenesis and function representing putative contributors to human phenotypes.

Effect of Prenatal Opioid Exposure on the Human Placental Methylome

Prenatal exposure to addictive drugs can lead to placental epigenetic modifications, but a methylome-wide evaluation of placental DNA methylation changes after prenatal opioid exposure has not yet been performed. Placental tissue samples were collected at delivery from 19 opioid-exposed and 20 unexposed control full-term pregnancies. Placental DNA methylomes were profiled using the Illumina Infinium HumanMethylationEPIC BeadChip. Differentially methylated CpG sites associated with opioid exposure were identified with a linear model using the ‘limma’ R package. To identify differentially methylated regions (DMRs) spanning multiple CpG sites, the ‘DMRcate’ R package was used. The functions of genes mapped by differentially methylated CpG sites and DMRs were further annotated using Enrichr. Differentially methylated CpGs (n = 684, unadjusted p < 0.005 and |∆β| ≥ 0.05) were mapped to 258 genes (including PLD1, MGAM, and ALCS2). Differentially methylated regions (n = 199) were located in 174 genes (including KCNMA1). Enrichment analysis of the top differentially methylated CpG sites and regions indicated disrupted epigenetic regulation of genes involved in synaptic structure, chemical synaptic transmission, and nervous system development. Our findings imply that placental epigenetic changes due to prenatal opioid exposure could result in placental dysfunction, leading to abnormal fetal brain development and the symptoms of opioid withdrawal in neonates.

In-Silico Functional Annotation of Plasmodium falciparum Hypothetical Proteins to Identify Novel Drug Targets

Plasmodium falciparum is one of the plasmodium species responsible for the majority of life-threatening malaria cases. The current antimalarial therapies are becoming less effective due to growing drug resistance, leading to the urgent requirement for alternative and more effective antimalarial drugs or vaccines. To facilitate the novel drug discovery or vaccine development efforts, recent advances in sequencing technologies provide valuable information about the whole genome of the parasite, yet a lot more needs to be deciphered due to its incomplete proteome annotation. Surprisingly, out of the 5,389 proteins currently annotated in the Plasmodium falciparum 3D7 strain, 1,626 proteins (∼30% data) are annotated as hypothetical proteins. In parasite genomic studies, the challenge to annotate hypothetical proteins is often ignored, which may obscure the crucial information related to the pathogenicity of the parasite. In this study, we attempt to characterize hypothetical proteins of the parasite to identify novel drug targets using a computational pipeline. The study reveals that out of the overall pool of the hypothetical proteins, 266 proteins have conserved functional signatures. Furthermore, the pathway analysis of these proteins revealed that 23 proteins have an essential role in various biochemical, signalling and metabolic pathways. Additionally, all the proteins (266) were subjected to computational structure analysis. We could successfully model 11 proteins. We validated and checked the structural stability of the models by performing molecular dynamics simulation. Interestingly, eight proteins show stable conformations, and seven proteins are specific for Plasmodium falciparum, based on homology analysis. Lastly, mapping the seven shortlisted hypothetical proteins on the Plasmodium falciparum protein-protein interaction network revealed 3,299 nodes and 2,750,692 edges. Our study revealed interesting functional details of seven hypothetical proteins of the parasite, which help learn more about the less-studied molecules and their interactions, providing valuable clues to unravel the role of these proteins via future experimental validation.

Comprehensive analysis of miRNA-mRNA interactions in ovaries of aged mice

Advanced maternal age and ovarian aging are deleterious to the quantity and quality of oocytes and epigenetic modifications, which can affect the health of offspring. However, relatively little is known about the regulation of microRNA-mediated transcription during ovarian aging. We therefore aimed to identify age-related mRNA and microRNA changes and their interactions in the ovaries of aged mice. We performed QuantSeq 3'mRNA and small RNA sequencing to compare their expression patterns in post-ovulation ovaries from young (12-week-old) and old (44-week-old) mice. Functional annotation and integrative analyses were performed to identify the potential functions of differentially expressed genes and identify binding sites for critical microRNAs. We found 343 differentially expressed genes and 9 microRNAs in our comparison of the two mouse groups, with fold changes >2.0 (P < 0.01). Furthermore, we identified possible direct interactions between 24 differentially expressed mRNAs and 8 microRNAs. The differentially expressed genes are involved in fat digestion and absorption, the PI3K-Akt signaling pathway, serotonergic synapse, and ovarian steroidogenesis, which are important for folliculogenesis and oocyte growth. During ovarian aging, changes in gene expression induce alterations in folliculogenesis, oocyte growth, and steroidogenesis, resulting in decreased oocyte quality and reproductive outcomes.

Machine Learning and Deep Learning Applications in Metagenomic Taxonomy and Functional Annotation

Shotgun sequencing of environmental DNA (i.e., metagenomics) has revolutionized the field of environmental microbiology, allowing the characterization of all microorganisms in a sequencing experiment. To identify the microbes in terms of taxonomy and biological activity, the sequenced reads must necessarily be aligned on known microbial genomes/genes. However, current alignment methods are limited in terms of speed and can produce a significant number of false positives when detecting bacterial species or false negatives in specific cases (virus, plasmids, and gene detection). Moreover, recent advances in metagenomics have enabled the reconstruction of new genomes using de novo binning strategies, but these genomes, not yet fully characterized, are not used in classic approaches, whereas machine and deep learning methods can use them as models. In this article, we attempted to review the different methods and their efficiency to improve the annotation of metagenomic sequences. Deep learning models have reached the performance of the widely used k-mer alignment-based tools, with better accuracy in certain cases; however, they still must demonstrate their robustness across the variety of environmental samples and across the rapid expansion of accessible genomes in databases.

Bacterial Communities of Forest Soils along Different Elevations: Diversity, Structure, and Functional Composition with Potential Impacts on CO2 Emission

Soil bacteria are important components of forest ecosystems, there compostion structure and functions are sensitive to environmental conditions along elevation gradients. Using 16S rRNA gene amplicon sequencing followed by FAPROTAX function prediction, we examined the diversity, composition, and functional potentials of soil bacterial communities at three sites at elevations of 1400 m, 1600 m, and 2200 m in a temperate forest. We showed that microbial taxonomic composition did not change with elevation (p = 0.311), though soil bacterial α-diversities did. Proteobacteria, Acidobacteria, Actinobacteria, and Verrucomicrobia were abundant phyla in almost all soil samples, while Nitrospirae, closely associated with soil nitrogen cycling, was the fourth most abundant phylum in soils at 2200 m. Chemoheterotrophy and aerobic chemoheterotrophy were the two most abundant functions performed in soils at 1400 m and 1600 m, while nitrification (25.59% on average) and aerobic nitrite oxidation (19.38% on average) were higher in soils at 2200 m. Soil CO2 effluxes decreased (p < 0.050) with increasing elevation, while they were positively correlated (r = 0.55, p = 0.035) with the abundances of bacterial functional groups associated with carbon degradation. Moreover, bacterial functional composition, rather than taxonomic composition, was significantly associated with soil CO2 effluxes, suggesting a decoupling of taxonomy and function, with the latter being a better predictor of ecosystem functions. Annual temperature, annual precipitation, and pH shaped (p < 0.050) both bacterial taxonomic and functional communities. By establishing linkages between bacterial taxonomic communities, abundances of bacterial functional groups, and soil CO2 fluxes, we provide novel insights into how soil bacterial communities could serve as potential proxies of ecosystem functions.

Hypothetical protein predicted to be tumor suppressor: a protein functional analysis

Litorilituus sediminis is a Gram-negative, aerobic, novel bacterium under the family of Colwelliaceae, has a stunning hypothetical protein containing domain called von Hippel-Lindau that has significant tumor suppressor activity. Therefore, this study was designed to elucidate the structure and function of the biologically important hypothetical protein EMK97_00595 (QBG34344.1) using several bioinformatics tools. The functional annotation exposed that the hypothetical protein is an extracellular secretory soluble signal peptide and contains the von Hippel-Lindau (VHL; VHL beta) domain that has a significant role in tumor suppression. This domain is conserved throughout evolution, as its homologs are available in various types of the organism like mammals, insects, and nematode. The gene product of VHL has a critical regulatory activity in the ubiquitous oxygen-sensing pathway. This domain has a significant role in inhibiting cell proliferation, angiogenesis progression, kidney cancer, breast cancer, and colon cancer. At last, the current study depicts that the annotated hypothetical protein is linked with tumor suppressor activity which might be of great interest to future research in the higher organism.

Genomic sequence and transcriptome analysis of the medicinal fungus Keithomyces neogunnii

The filamentous fungus Keithomyces neogunnii can infect the larvae of Lepidoptera (Hepialus sp.) and form an insect–fungi complex, which is utilized as an important traditional Chinese medicine. As a valuable medicinal fungus, K. neogunnii produces diverse bioactive substances (e.g., polysaccharide, vitamins, cordycepic acid, and adenosine) under cultivation conditions. Herein, we report the first high-quality genome of the K. neogunnii single-spore isolate Cg7.2a using single-molecule real-time sequencing technology in combination with Illumina sequencing. The assembled genome was 32.6 Mb in size, containing 8,641 predicted genes and having a GC content of 52.16%. RNA sequencing analysis revealed the maximum number of differentially expressed genes in the fungus during the stroma formation stage compared with those during the mycelium stage. These data are valuable to enhance our understanding of the biology, development, evolution, and physiological metabolism of K. neogunnii.

Identification of Iron-Sulfur (Fe-S) Cluster and Zinc (Zn) Binding Sites Within Proteomes Predicted by DeepMind's AlphaFold2 Program Dramatically Expands the Metalloproteome

DeepMind's AlphaFold2 software has ushered in a revolution in high quality, 3D protein structure prediction. In very recent work by the DeepMind team, structure predictions have been made for entire proteomes of twenty-one organisms, with >360,000 structures made available for download. Here we show that thousands of novel binding sites for iron-sulfur (Fe-S) clusters and zinc (Zn) ions can be identified within these predicted structures by exhaustive enumeration of all potential ligand-binding orientations. We demonstrate that AlphaFold2 routinely makes highly specific predictions of ligand binding sites: for example, binding sites that are comprised exclusively of four cysteine sidechains fall into three clusters, representing binding sites for 4Fe-4S clusters, 2Fe-2S clusters, or individual Zn ions. We show further: (a) that the majority of known Fe-S cluster and Zn binding sites documented in UniProt are recovered by the AlphaFold2 structures, (b) that there are occasional disputes between AlphaFold2 and UniProt with AlphaFold2 predicting highly plausible alternative binding sites, (c) that the Fe-S cluster binding sites that we identify in E. coli agree well with previous bioinformatics predictions, (d) that cysteines predicted here to be part of ligand binding sites show little overlap with those shown via chemoproteomics techniques to be highly reactive, and (e) that AlphaFold2 occasionally appears to build erroneous disulfide bonds between cysteines that should instead coordinate a ligand. These results suggest that AlphaFold2 could be an important tool for the functional annotation of proteomes, and the methodology presented here is likely to be useful for predicting other ligand-binding sites.

In-silico characterization and structure-based functional annotation of a hypothetical protein from Campylobacter jejuni involved in propionate catabolism

Campylobacter jejuni is one of the most prevalent organisms associated with foodborne illness across the globe causing campylobacteriosis and gastritis. Many proteins of C. jejuni are still unidentified. The purpose of this study was to determine the structure and function of a non-annotated hypothetical protein (HP) from C. jejuni. A number of properties like physiochemical characteristics, 3D structure, and functional annotation of the HP (accession No. CAG2129885.1) were predicted using various bioinformatics tools followed by further validation and quality assessment. Moreover, the protein-protein interactions and active site were obtained from the STRING and CASTp server, respectively. The hypothesized protein possesses various characteristics including an acidic pH, thermal stability, water solubility, and cytoplasmic distribution. While alpha-helix and random coil structures are the most prominent structural components of this protein, most of it is formed of helices and coils. Along with expected quality, the 3D model has been found to be novel. This study has identified the potential role of the HP in 2-methylcitric acid cycle and propionate catabolism. Furthermore, protein-protein interactions revealed several significant functional partners. The in-silico characterization of this protein will assist to understand its molecular mechanism of action better. The methodology of this study would also serve as the basis for additional research into proteomic and genomic data for functional potential identification.

A Glance into MTHFR Deficiency at a Molecular Level

MTHFR deficiency still deserves an investigation to associate the phenotype to protein structure variations. To this aim, considering the MTHFR wild type protein structure, with a catalytic and a regulatory domain and taking advantage of state-of-the-art computational tools, we explore the properties of 72 missense variations known to be disease associated. By computing the thermodynamic ΔΔG change according to a consensus method that we recently introduced, we find that 61% of the disease-related variations destabilize the protein, are present both in the catalytic and regulatory domain and correspond to known biochemical deficiencies. The propensity of solvent accessible residues to be involved in protein-protein interaction sites indicates that most of the interacting residues are located in the regulatory domain, and that only three of them, located at the interface of the functional protein homodimer, are both disease-related and destabilizing. Finally, we compute the protein architecture with Hidden Markov Models, one from Pfam for the catalytic domain and the second computed in house for the regulatory domain. We show that patterns of disease-associated, physicochemical variation types, both in the catalytic and regulatory domains, are unique for the MTHFR deficiency when mapped into the protein architecture.

Whole Genome Sequencing and Annotation of Naematelia aurantialba (Basidiomycota, Edible-Medicinal Fungi)

Naematelia aurantialba is a rare edible fungus with both nutritional and medicinal values and especially rich in bioactive polysaccharides. However, due to the lack of genomic information, researches on the mining of active compounds, artificial breeding and cultivation, genetics, and molecular biology are limited. To facilitate the medicinal and food applications of N. aurantialba, we sequenced and analyzed the whole genome of N. aurantialba for the first time. The 21-Mb genome contained 15 contigs, and a total of 5860 protein-coding genes were predicted. The genome sequence shows that 296 genes are related to polysaccharide synthesis, including 15 genes related to nucleoside-activated sugar synthesis and 11 genes related to glucan synthesis. The genome also contains genes and gene clusters for the synthesis of other active substances, including terpenoids, unsaturated fatty acids, and bioactive proteins. In addition, it was also found that N. aurantialba was more closely related to Naematelia encephala than to Tremella fuciformis. In short, this study provides a reference for molecular cognition of N. aurantialba and related researches.

AMPing Up the Search: A Structural and Functional Repository of Antimicrobial Peptides for Biofilm Studies, and a Case Study of Its Application to Corynebacterium striatum, an Emerging Pathogen

Antimicrobial peptides (AMPs) have been recognized for their ability to target processes important for biofilm formation. Given the vast array of AMPs, identifying potential anti-biofilm candidates remains a significant challenge, and prompts the need for preliminary in silico investigations prior to extensive in vitro and in vivo studies. We have developed Biofilm-AMP (B-AMP), a curated 3D structural and functional repository of AMPs relevant to biofilm studies. In its current version, B-AMP contains predicted 3D structural models of 5544 AMPs (from the DRAMP database) developed using a suite of molecular modeling tools. The repository supports a user-friendly search, using source, name, DRAMP ID, and PepID (unique to B-AMP). Further, AMPs are annotated to existing biofilm literature, consisting of a vast library of over 10,000 articles, enhancing the functional capabilities of B-AMP. To provide an example of the usability of B-AMP, we use the sortase C biofilm target of the emerging pathogen Corynebacterium striatum as a case study. For this, 100 structural AMP models from B-AMP were subject to in silico protein-peptide molecular docking against the catalytic site residues of the C. striatum sortase C protein. Based on docking scores and interacting residues, we suggest a preference scale using which candidate AMPs could be taken up for further in silico, in vitro and in vivo testing. The 3D protein-peptide interaction models and preference scale are available in B-AMP. B-AMP is a comprehensive structural and functional repository of AMPs, and will serve as a starting point for future studies exploring AMPs for biofilm studies. B-AMP is freely available to the community at https://b-amp.karishmakaushiklab.com and will be regularly updated with AMP structures, interaction models with potential biofilm targets, and annotations to biofilm literature.

BGvar: A comprehensive resource for blood group immunogenetics

BACKGROUND

Blood groups form the basis of effective and safe blood transfusion. There are about 43 well-recognised human blood group systems presently known. Blood groups are molecularly determined by the presence of specific antigens on the red blood cells and are genetically determined and inherited following Mendelian principles. The lack of a comprehensive, relevant, manually compiled and genome-ready dataset of red cell antigens limited the widespread application of genomic technologies to characterise and interpret the blood group complement of an individual from genomic datasets.

MATERIALS AND METHODS

A range of public datasets was used to systematically annotate the variation compendium for its functionality and allele frequencies across global populations. Details on phenotype or relevant clinical importance were collated from reported literature evidence.

RESULTS

We have compiled the Blood Group Associated Genomic Variant Resource (BGvar), a manually curated online resource comprising all known human blood group related allelic variants including a total of 1700 International Society of Blood Transfusion approved alleles and 1706 alleles predicted and curated from literature reports. This repository includes 1682 single nucleotide variations (SNVs), 310 Insertions, Deletions (InDels) and Duplications (Copy Number Variations) and about 1360 combination mutations corresponding to 43 human blood group systems and 2 transcription factors. This compendium also encompasses gene fusion and rearrangement events occurring in human blood group genes.

CONCLUSION

To the best of our knowledge, BGvar is a comprehensive and a user-friendly resource with most relevant collation of blood group alleles in humans. BGvar is accessible online at URL: http://clingen.igib.res.in/bgvar/.

eggNOG-mapper v2: Functional Annotation, Orthology Assignments, and Domain Prediction at the Metagenomic Scale

Even though automated functional annotation of genes represents a fundamental step in most genomic and metagenomic workflows, it remains challenging at large scales. Here, we describe a major upgrade to eggNOG-mapper, a tool for functional annotation based on precomputed orthology assignments, now optimized for vast (meta)genomic data sets. Improvements in version 2 include a full update of both the genomes and functional databases to those from eggNOG v5, as well as several efficiency enhancements and new features. Most notably, eggNOG-mapper v2 now allows for: 1) de novo gene prediction from raw contigs, 2) built-in pairwise orthology prediction, 3) fast protein domain discovery, and 4) automated GFF decoration. eggNOG-mapper v2 is available as a standalone tool or as an online service at http://eggnog-mapper.embl.de.

eggNOG-mapper v2: Functional Annotation, Orthology Assignments, and Domain Prediction at the Metagenomic Scale

Even though automated functional annotation of genes represents a fundamental step in most genomic and metagenomic workflows, it remains challenging at large scales. Here, we describe a major upgrade to eggNOG-mapper, a tool for functional annotation based on precomputed orthology assignments, now optimized for vast (meta)genomic data sets. Improvements in version 2 include a full update of both the genomes and functional databases to those from eggNOG v5, as well as several efficiency enhancements and new features. Most notably, eggNOG-mapper v2 now allows for: 1) de novo gene prediction from raw contigs, 2) built-in pairwise orthology prediction, 3) fast protein domain discovery, and 4) automated GFF decoration. eggNOG-mapper v2 is available as a standalone tool or as an online service at http://eggnog-mapper.embl.de.

Trait-specific Selection Signature Detection Reveals Novel Loci of Meat Quality in Large White Pigs

In past decades, meat quality traits have been shaped by human-driven selection in the process of genetic improvement programs. Exploring the potential genetic basis of artificial selection and mapping functional candidate genes for economic traits are of great significance in genetic improvement of pigs. In this study, we focus on investigating the genetic basis of five meat quality traits, including intramuscular fat content (IMF), drip loss, water binding capacity, pH at 45 min (pH45min), and ultimate pH (pH24h). Through making phenotypic gradient differential population pairs, Wright’s fixation index (F_ST) and the cross-population extended haplotype homozogysity (XPEHH) were applied to detect selection signatures for these five traits. Finally, a total of 427 and 307 trait-specific selection signatures were revealed by F_ST and XPEHH, respectively. Further bioinformatics analysis indicates that some genes, such as USF1, NDUFS2, PIGM, IGSF8, CASQ1, and ACBD6, overlapping with the trait-specific selection signatures are responsible for the phenotypes including fat metabolism and muscle development. Among them, a series of promising trait-specific selection signatures that were detected in the high IMF subpopulation are located in the region of 93544042-95179724bp on SSC4, and the genes harboring in this region are all related to lipids and muscle development. Overall, these candidate genes of meat quality traits identified in this analysis may provide some fundamental information for further exploring the genetic basis of this complex trait.

Identification of a Novel Therapeutic Target against XDR Salmonella Typhi H58 Using Genomics Driven Approach Followed Up by Natural Products Virtual Screening

Typhoid fever is caused by a pathogenic, rod-shaped, flagellated, and Gram-negative bacterium known as Salmonella Typhi. It features a polysaccharide capsule that acts as a virulence factor and deceives the host immune system by protecting phagocytosis. Typhoid fever remains a major health concern in low and middle-income countries, with an estimated death rate of ~200,000 per annum. However, the situation is exacerbated by the emergence of the extensively drug-resistant (XDR) strain designated as H58 of S. Typhi. The emergence of the XDR strain is alarming, and it poses serious threats to public health due to the failure of the current therapeutic regimen. A relatively newer computational method called subtractive genomics analyses has been widely applied to discover novel and new drug targets against pathogens, particularly drug-resistant ones. The method involves the gradual reduction of the complete proteome of the pathogen, leading to few potential and novel drug targets. Thus, in the current study, a subtractive genomics approach was applied against the Salmonella XDR strain to identify potential drug targets. The current study predicted four prioritized proteins (i.e., Colanic acid biosynthesis acetyltransferase wcaB, Shikimate dehydrogenase aroE, multidrug efflux RND transporter permease subunit MdtC, and pantothenate synthetase panC) as potential drug targets. Though few of the prioritized proteins are treated in the literature as the established drug targets against other pathogenic bacteria, these drug targets are identified here for the first time against S. Typhi (i.e., S. Typhi XDR). The current study aimed at drawing attention to new drug targets against S. Typhi that remain largely unexplored. One of the prioritized drug targets, i.e., Colanic acid biosynthesis acetyltransferase, was predicted as a unique, new drug target against S. Typhi XDR. Therefore, the Colanic acid was further explored using structure-based techniques. Additionally, ~1000 natural compounds were docked with Colanic acid biosynthesis acetyltransferase, resulting in the prediction of seven compounds as potential lead candidates against the S. Typhi XDR strain. The ADMET properties and binding energies via the docking program of these seven compounds characterized them as novel drug candidates. They may potentially be used for the development of future drugs in the treatment of Typhoid fever.

Complete Genome Sequence of the Type Strain Tepidimonas taiwanensis LMG 22826T, a Thermophilic Alkaline Protease and Polyhydroxyalkanoate Producer

Tepidimonas taiwanensis is a moderately thermophilic, Gram-negative, rod-shaped, chemoorganoheterotrophic, motile bacterium. The alkaline protease producing type strain T. taiwanensis LMG 22826^T was recently reported to also be a promising producer of polyhydroxyalkanoates (PHAs)—renewable and biodegradable polymers representing an alternative to conventional plastics. Here, we present its first complete genome sequence which is also the first complete genome sequence of the whole species. The genome consists of a single 2,915,587-bp-long circular chromosome with GC content of 68.75%. Genome annotation identified 2,764 genes in total while 2,634 open reading frames belonged to protein-coding genes. Although functional annotation of the genome and division of genes into Clusters of Orthologous Groups (COGs) revealed a relatively high number of 694 genes with unknown function or unknown COG, the majority of genes were assigned a function. Most of the genes, 406 in total, were involved in energy production and conversion, and amino acid transport and metabolism. Moreover, particular key genes involved in the metabolism of PHA were identified. Knowledge of the genome in connection with the recently reported ability to produce bioplastics from the waste stream of wine production makes T. taiwanensis LMG 22826^T, an ideal candidate for further genome engineering as a bacterium with high biotechnological potential.

Gaining insight into metabolic diseases from human genetic discoveries

Human large-scale genetic association studies have identified sequence variations at thousands of genetic risk loci that are more common in patients with diverse metabolic disease compared with healthy controls. While these genetic associations have been replicated in multiple large cohorts and sometimes can explain up to 50% of heritability, the molecular and cellular mechanisms affected by common genetic variation associated with metabolic disease remains mostly unknown. A variety of new genome-wide data types, in conjunction with novel biostatistical and computational analytical methodologies and foundational experimental technologies, are paving the way for a principled approach to systematic variant-to-function (V2F) studies for metabolic diseases, turning associated regions into causal variants, cell types and states of action, effector genes, and cellular and physiological mechanisms. Identification of new target genes and cellular programs for metabolic risk loci will improve mechanistic understanding of disease biology and identification of novel therapeutic strategies.

Proteomic analysis of the colistin-resistant E. coli clinical isolate: Explorations of the resistome

BACKGROUND

Antimicrobial resistance is a worldwide problem after the emergence of colistin resistance since it was the last option left to treat carbapenemase-resistant bacterial infections. The mcr gene and its variants are one of the causes for colistin resistance. Besides mcr genes, some other intrinsic genes are also involved in colistin resistance but still need to be explored.

OBJECTIVE

The aim of this study was to investigate differential proteins expression of colistin-resistant E. coli clinical isolate and to understand their interactive partners as future drug targets.

METHODS

In this study, we have employed the whole proteome analysis through LC-MS/MS. The advance proteomics tools were used to find differentially expressed proteins in the colistin-resistant Escherichia coli clinical isolate compared to susceptible isolate. Gene ontology and STRING were used for functional annotation and protein-protein interaction networks, respectively.

RESULTS

LC-MS/MS analysis showed overexpression of 47 proteins and underexpression of 74 proteins in colistin-resistant E. coli. These proteins belong to DNA replication, transcription and translational process; defense and stress related proteins; proteins of phosphoenol pyruvate phosphotransferase system (PTS) and sugar metabolism. Functional annotation and protein-protein interaction showed translational and cellular metabolic process, sugar metabolism and metabolite interconversion.

CONCLUSION

We conclude that these protein targets and their pathways might be used to develop novel therapeutics against colistin-resistant infections. These proteins could unveil the mechanism of colistin resistance.

Decoding the Equine Genome: Lessons from ENCODE

The horse reference genome assemblies, EquCab2.0 and EquCab3.0, have enabled great advancements in the equine genomics field, from tools to novel discoveries. However, significant gaps of knowledge regarding genome function remain, hindering the study of complex traits in horses. In an effort to address these gaps and with inspiration from the Encyclopedia of DNA Elements (ENCODE) project, the equine Functional Annotation of Animal Genome (FAANG) initiative was proposed to bridge the gap between genome and gene expression, providing further insights into functional regulation within the horse genome. Three years after launching the initiative, the equine FAANG group has generated data from more than 400 experiments using over 50 tissues, targeting a variety of regulatory features of the equine genome. In this review, we examine how valuable lessons learned from the ENCODE project informed our decisions in the equine FAANG project. We report the current state of the equine FAANG project and discuss how FAANG can serve as a template for future expansion of functional annotation in the equine genome and be used as a reference for studies of complex traits in horse. A well-annotated reference functional atlas will also help advance equine genetics in the pan-genome and precision medicine era.

FA-nf: A Functional Annotation Pipeline for Proteins from Non-Model Organisms Implemented in Nextflow

Functional annotation allows adding biologically relevant information to predicted features in genomic sequences, and it is, therefore, an important procedure of any de novo genome sequencing project. It is also useful for proofreading and improving gene structural annotation. Here, we introduce FA-nf, a pipeline implemented in Nextflow, a versatile computational workflow management engine. The pipeline integrates different annotation approaches, such as NCBI BLAST+, DIAMOND, InterProScan, and KEGG. It starts from a protein sequence FASTA file and, optionally, a structural annotation file in GFF format, and produces several files, such as GO assignments, output summaries of the abovementioned programs and final annotation reports. The pipeline can be broken easily into smaller processes for the purpose of parallelization and easily deployed in a Linux computational environment, thanks to software containerization, thus helping to ensure full reproducibility.

From GWAS to Gene: Transcriptome-Wide Association Studies and Other Methods to Functionally Understand GWAS Discoveries

Since their inception, genome-wide association studies (GWAS) have identified more than a hundred thousand single nucleotide polymorphism (SNP) loci that are associated with various complex human diseases or traits. The majority of GWAS discoveries are located in non-coding regions of the human genome and have unknown functions. The valley between non-coding GWAS discoveries and downstream affected genes hinders the investigation of complex disease mechanism and the utilization of human genetics for the improvement of clinical care. Meanwhile, advances in high-throughput sequencing technologies reveal important genomic regulatory roles that non-coding regions play in the transcriptional activities of genes. In this review, we focus on data integrative bioinformatics methods that combine GWAS with functional genomics knowledge to identify genetically regulated genes. We categorize and describe two types of data integrative methods. First, we describe fine-mapping methods. Fine-mapping is an exploratory approach that calibrates likely causal variants underneath GWAS signals. Fine-mapping methods connect GWAS signals to potentially causal genes through statistical methods and/or functional annotations. Second, we discuss gene-prioritization methods. These are hypothesis generating approaches that evaluate whether genetic variants regulate genes via certain genetic regulatory mechanisms to influence complex traits, including colocalization, mendelian randomization, and the transcriptome-wide association study (TWAS). TWAS is a gene-based association approach that investigates associations between genetically regulated gene expression and complex diseases or traits. TWAS has gained popularity over the years due to its ability to reduce multiple testing burden in comparison to other variant-based analytic approaches. Multiple types of TWAS methods have been developed with varied methodological designs and biological hypotheses over the past 5 years. We dive into discussions of how TWAS methods differ in many aspects and the challenges that different TWAS methods face. Overall, TWAS is a powerful tool for identifying complex trait-associated genes. With the advent of single-cell sequencing, chromosome conformation capture, gene editing technologies, and multiplexing reporter assays, we are expecting a more comprehensive understanding of genomic regulation and genetically regulated genes underlying complex human diseases and traits in the future.

Drug Repurposing for Atopic Dermatitis by Integration of Gene Networking and Genomic Information

Atopic Dermatitis (AD) is a chronic and relapsing skin disease. The medications for treating AD are still limited, most of them are topical corticosteroid creams or antibiotics. The current study attempted to discover potential AD treatments by integrating a gene network and genomic analytic approaches. Herein, the Single Nucleotide Polymorphism (SNPs) associated with AD were extracted from the GWAS catalog. We identified 70 AD-associated loci, and then 94 AD risk genes were found by extending to proximal SNPs based on r2 > 0.8 in Asian populations using HaploReg v4.1. Next, we prioritized the AD risk genes using in silico pipelines of bioinformatic analysis based on six functional annotations to identify biological AD risk genes. Finally, we expanded them according to the molecular interactions using the STRING database to find the drug target genes. Our analysis showed 27 biological AD risk genes, and they were mapped to 76 drug target genes. According to DrugBank and Therapeutic Target Database, 25 drug target genes overlapping with 53 drugs were identified. Importantly, dupilumab, which is approved for AD, was successfully identified in this bioinformatic analysis. Furthermore, ten drugs were found to be potentially useful for AD with clinical or preclinical evidence. In particular, we identified filgotinub and fedratinib, targeting gene JAK1, as potential drugs for AD. Furthermore, four monoclonal antibody drugs (lebrikizumab, tralokinumab, tocilizumab, and canakinumab) were successfully identified as promising for AD repurposing. In sum, the results showed the feasibility of gene networking and genomic information as a potential drug discovery resource.

Disease category-specific annotation of variants using an ensemble learning framework

Understanding the impact of non-coding sequence variants on complex diseases is an essential problem. We present a novel ensemble learning framework-CASAVA, to predict genomic loci in terms of disease category-specific risk. Using disease-associated variants identified by GWAS as training data, and diverse sequencing-based genomics and epigenomics profiles as features, CASAVA provides risk prediction of 24 major categories of diseases throughout the human genome. Our studies showed that CASAVA scores at a genomic locus provide a reasonable prediction of the disease-specific and disease category-specific risk prediction for non-coding variants located within the locus. Taking MHC2TA and immune system diseases as an example, we demonstrate the potential of CASAVA in revealing variant-disease associations. A website (http://zhanglabtools.org/CASAVA) has been built to facilitate easily access to CASAVA scores.