The InterPro protein families database: the classification resource after 15 years

Water wisteria genome reveals environmental adaptation and heterophylly regulation in amphibious plants

Heterophylly is a phenomenon whereby an individual plant dramatically changes leaf shape in response to the surroundings. Hygrophila difformis (Acanthaceae; water wisteria), has recently emerged as a model plant to study heterophylly because of its striking leaf shape variation in response to various environmental factors. When submerged, H. difformis often develops complex leaves, but on land it develops simple leaves. Leaf complexity is also influenced by other factors, such as light density, humidity, and temperature. Here, we sequenced and assembled the H. difformis chromosome-level genome (scaffold N50: 60.43 Mb, genome size: 871.92 Mb), which revealed 36 099 predicted protein-coding genes distributed over 15 pseudochromosomes. H. difformis diverged from its relatives during the Oligocene climate-change period and expanded gene families related to its amphibious habit. Genes related to environmental stimuli, leaf development, and other pathways were differentially expressed in submerged and terrestrial conditions, possibly modulating morphological and physiological acclimation to changing environments. We also found that auxin plays a role in H. difformis heterophylly. Finally, we discovered candidate genes that respond to different environmental conditions and elucidated the role of LATE MERISTEM IDENTITY 1 (LMI1) in heterophylly. We established H. difformis as a model for studying interconnections between environmental adaptation and morphogenesis.

Coupling CO2 Reduction and Acetyl-CoA Formation: The Role of a CO Capturing Tunnel in Enzymatic Catalysis

The bifunctional CO-dehydrogenase/acetyl-CoA synthase (CODH/ACS) complex couples the reduction of CO2 to the condensation of CO with a methyl moiety and CoA to acetyl-CoA. Catalysis occurs at two sites connected by a tunnel transporting the CO. In this study, we investigated how the bifunctional complex and its tunnel support catalysis using the CODH/ACS from Carboxydothermus hydrogenoformans as a model. Although CODH/ACS adapted to form a stable bifunctional complex with a secluded substrate tunnel, catalysis and CO transport is even more efficient when two monofunctional enzymes are coupled. Efficient CO channeling appears to be ensured by hydrophobic binding sites for CO, which act in a bucket-brigade fashion rather than as a simple tube. Tunnel remodeling showed that opening the tunnel increased activity but impaired directed transport of CO. Constricting the tunnel impaired activity and CO transport, suggesting that the tunnel evolved to sequester CO rather than to maximize turnover.

Comparative transcriptomics identifies genes underlying growth performance of the Pacific black-lipped pearl oyster Pinctada margaritifera

BACKGROUND

In bivalves, the rate at which organisms grow is a major functional trait underlying many aspects of their commercial production. Growth is a highly polygenic trait, which is typically regulated by many genes with small to moderate effects. Due to its complexity, growth variability in such shellfish remains poorly understood. In this study, we aimed to investigate differential gene expression among spat of the pearl oyster Pinctada margaritifera with distinct growth phenotypes.

RESULTS

We selected two groups of P. margaritifera spat belonging to the same F2 cohort based on their growth performance at 5.5 months old. Transcriptome profile analysis identified a total of 394 differentially expressed genes between these Fast-growing (F) and Slow-growing (S) phenotypes. According to functional enrichment analysis, S oysters overexpressed genes associated with stress-pathways and regulation of innate immune responses. In contrast, F oysters up-regulated genes associated with cytoskeleton activity, cell proliferation, and apoptosis. Analysis of genome polymorphism identified 16 single nucleotide polymorphisms (SNPs) significantly associated with the growth phenotypes. SNP effect categorization revealed one SNP identified for high effect and annotated for a stop codon gained mutation. Interestingly, this SNP is located within a gene annotated for scavenger receptor class F member 1 (SRF1), which is known to modulate apoptosis. Our analyses also revealed that all F oysters showed up-regulation for this gene and were homozygous for the stop-codon mutation. Conversely, S oysters had a heterozygous genotype and a reduced expression of this gene.

CONCLUSIONS

Altogether, our findings suggest that differences in growth among the same oyster cohort may be explained by contrasted metabolic allocation between regulatory pathways for growth and the immune system. This study provides a valuable contribution towards our understanding of the molecular components associated with growth performance in the pearl oyster P. margaritifera and bivalves in general.

Major changes in domain arrangements are associated with the evolution of termites

Domains as functional protein units and their rearrangements along the phylogeny can shed light on the functional changes of proteomes associated with the evolution of complex traits like eusociality. This complex trait is associated with sterile soldiers and workers, and long-lived, highly fecund reproductives. Unlike in Hymenoptera (ants, bees, and wasps), the evolution of eusociality within Blattodea, where termites evolved from within cockroaches, was accompanied by a reduction in proteome size, raising the question of whether functional novelty was achieved with existing rather than novel proteins. To address this, we investigated the role of domain rearrangements during the evolution of termite eusociality. Analysing domain rearrangements in the proteomes of three solitary cockroaches and five eusocial termites, we inferred more than 5,000 rearrangements over the phylogeny of Blattodea. The 90 novel domain arrangements that emerged at the origin of termites were enriched for several functions related to longevity, such as protein homeostasis, DNA repair, mitochondrial activity, and nutrient sensing. Many domain rearrangements were related to changes in developmental pathways, important for the emergence of novel castes. Along with the elaboration of social complexity, including permanently sterile workers and larger, foraging colonies, we found 110 further domain arrangements with functions related to protein glycosylation and ion transport. We found an enrichment of caste-biased expression and splicing within rearranged genes, highlighting their importance for the evolution of castes. Furthermore, we found increased levels of DNA methylation among rearranged compared to non-rearranged genes suggesting fundamental differences in their regulation. Our findings indicate the importance of domain rearrangements in the generation of functional novelty necessary for termite eusociality to evolve.

Chromosome-level genome assembly of Helwingia omeiensis: the first genome in the family Helwingiaceae

Helwingia, a shrub of the monotypic cosmopolitan family Helwingiaceae, is distinguished by its inflorescence, in which flowers are borne on the midrib of the leaf-a trait not commonly observed in related plant families. Previous studies have investigated the development of this unusual structure using comparative anatomical methods. However, the scarcity of genomic data has hindered our understanding of the origins and evolutionary history of this uncommon trait at the molecular level. Here, we report the first high-quality genome of the family Helwingiaceae. Assembled using HiFi sequencing and Hi-C technologies, the genome of H. omeiensis is anchored to 19 chromosomes, with a total length of 2.75 Gb and a contig N50 length of 6.78 Mb. The BUSCO completeness score of the assembled genome was 98.2%. 53,951 genes were identified, of which 99.7% were annotated in at least one protein database. The high-quality reference genome of H. omeiensis provides an essential genetic resource and sheds light on the phylogeny and evolution of specific traits in the family Helwingiaceae.

Chromosome-level genome assembly of Cnidium monnieri, a highly demanded traditional Chinese medicine

Cnidium monnieri, a medicinal herb of the Cnidium genus and the Apiaceae family, is among the most important traditional Chinese medicines and is widely distributed in China. However, to date, no C. monnieri-related genomic information has been described. In this study, we assembled the C. monnieri genome of approximately 1210.23 Mb with a contig N50 of 83.14 Mb. Using PacBio HiFi and Hi-C sequencing data, we successfully anchored 93.86% of the assembled sequences to 10 pseudochromosomes (2n = 20). We predicted a total of 37,460 protein-coding genes, with 97.02% of them being functionally annotated in Non-Redundant, Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and other databases. In addition, we identified 2,778 tRNAs, 4,180 rRNAs, 258 miRNAs, and 1,700 snRNAs in the genome. This is the first reported C. monnieri genome. Hopefully, the availability of this chromosome-level reference genome provides a significant basis for upcoming natural product-related biosynthetic pathway assessment in C. monnieri.

SAFPred: synteny-aware gene function prediction for bacteria using protein embeddings

MOTIVATION

Today, we know the function of only a small fraction of the protein sequences predicted from genomic data. This problem is even more salient for bacteria, which represent some of the most phylogenetically and metabolically diverse taxa on Earth. This low rate of bacterial gene annotation is compounded by the fact that most function prediction algorithms have focused on eukaryotes, and conventional annotation approaches rely on the presence of similar sequences in existing databases. However, often there are no such sequences for novel bacterial proteins. Thus, we need improved gene function prediction methods tailored for bacteria. Recently, transformer-based language models-adopted from the natural language processing field-have been used to obtain new representations of proteins, to replace amino acid sequences. These representations, referred to as protein embeddings, have shown promise for improving annotation of eukaryotes, but there have been only limited applications on bacterial genomes.

RESULTS

To predict gene functions in bacteria, we developed SAFPred, a novel synteny-aware gene function prediction tool based on protein embeddings from state-of-the-art protein language models. SAFpred also leverages the unique operon structure of bacteria through conserved synteny. SAFPred outperformed both conventional sequence-based annotation methods and state-of-the-art methods on multiple bacterial species, including for distant homolog detection, where the sequence similarity to the proteins in the training set was as low as 40%. Using SAFPred to identify gene functions across diverse enterococci, of which some species are major clinical threats, we identified 11 previously unrecognized putative novel toxins, with potential significance to human and animal health.

AVAILABILITY AND IMPLEMENTATION

https://github.com/AbeelLab/safpred.

Mutual annotation-based prediction of protein domain functions with Domain2GO

Identifying unknown functional properties of proteins is essential for understanding their roles in both health and disease states. The domain composition of a protein can reveal critical information in this context, as domains are structural and functional units that dictate how the protein should act at the molecular level. The expensive and time-consuming nature of wet-lab experimental approaches prompted researchers to develop computational strategies for predicting the functions of proteins. In this study, we proposed a new method called Domain2GO that infers associations between protein domains and function-defining gene ontology (GO) terms, thus redefining the problem as domain function prediction. Domain2GO uses documented protein-level GO annotations together with proteins' domain annotations. Co-annotation patterns of domains and GO terms in the same proteins are examined using statistical resampling to obtain reliable associations. As a use-case study, we evaluated the biological relevance of examples selected from the Domain2GO-generated domain-GO term mappings via literature review. Then, we applied Domain2GO to predict unknown protein functions by propagating domain-associated GO terms to proteins annotated with these domains. For function prediction performance evaluation and comparison against other methods, we employed Critical Assessment of Function Annotation 3 (CAFA3) challenge datasets. The results demonstrated the high potential of Domain2GO, particularly for predicting molecular function and biological process terms, along with advantages such as producing interpretable results and having an exceptionally low computational cost. The approach presented here can be extended to other ontologies and biological entities to investigate unknown relationships in complex and large-scale biological data. The source code, datasets, results, and user instructions for Domain2GO are available at https://github.com/HUBioDataLab/Domain2GO. Additionally, we offer a user-friendly online tool at https://huggingface.co/spaces/HUBioDataLab/Domain2GO, which simplifies the prediction of functions of previously unannotated proteins solely using amino acid sequences.

A comprehensive review and comparison of existing computational methods for protein function prediction

Protein function prediction is critical for understanding the cellular physiological and biochemical processes, and it opens up new possibilities for advancements in fields such as disease research and drug discovery. During the past decades, with the exponential growth of protein sequence data, many computational methods for predicting protein function have been proposed. Therefore, a systematic review and comparison of these methods are necessary. In this study, we divide these methods into four different categories, including sequence-based methods, 3D structure-based methods, PPI network-based methods and hybrid information-based methods. Furthermore, their advantages and disadvantages are discussed, and then their performance is comprehensively evaluated and compared. Finally, we discuss the challenges and opportunities present in this field.

Protein function prediction through multi-view multi-label latent tensor reconstruction

BACKGROUND

In last two decades, the use of high-throughput sequencing technologies has accelerated the pace of discovery of proteins. However, due to the time and resource limitations of rigorous experimental functional characterization, the functions of a vast majority of them remain unknown. As a result, computational methods offering accurate, fast and large-scale assignment of functions to new and previously unannotated proteins are sought after. Leveraging the underlying associations between the multiplicity of features that describe proteins could reveal functional insights into the diverse roles of proteins and improve performance on the automatic function prediction task.

RESULTS

We present GO-LTR, a multi-view multi-label prediction model that relies on a high-order tensor approximation of model weights combined with non-linear activation functions. The model is capable of learning high-order relationships between multiple input views representing the proteins and predicting high-dimensional multi-label output consisting of protein functional categories. We demonstrate the competitiveness of our method on various performance measures. Experiments show that GO-LTR learns polynomial combinations between different protein features, resulting in improved performance. Additional investigations establish GO-LTR's practical potential in assigning functions to proteins under diverse challenging scenarios: very low sequence similarity to previously observed sequences, rarely observed and highly specific terms in the gene ontology.

IMPLEMENTATION

The code and data used for training GO-LTR is available at https://github.com/aalto-ics-kepaco/GO-LTR-prediction .

Examining the impacts of salt specificity on freshwater microbial community and functional potential following salinization

The degradation of freshwater systems by salt pollution is a threat to global freshwater resources. Salinization is commonly identified by increased specific conductance (conductivity), a proxy for salt concentrations. However, conductivity fails to account for the diversity of salts entering freshwaters and the potential implications this has on microbial communities and functions. We tested 4 types of salt pollution-MgCl2, MgSO4, NaCl, and Na2SO4-on bacterial taxonomic and functional α-, β-diversity of communities originating from streams in two distinct localities (Nebraska [NE] and Ohio [OH], USA). Community responses depended on the site of origin, with NE and OH exhibiting more pronounced decreases in community diversity in response to Na2SO4 and MgCl2 than other salt amendments. A closer examination of taxonomic and functional diversity metrics suggests that core features of communities are more resistant to induced salt stress and that marginal features at both a population and functional level are more likely to exhibit significant structural shifts based on salt specificity. The lack of uniformity in community response highlights the need to consider the compositional complexities of salinization to accurately identify the ecological consequences of instances of salt pollution.

Chromosome-level genome assembly of the yellow-cheek carp Elopichthys bambusa

Yellow-cheek carp (Elopichthys bambusa) is a typical large and ferocious carnivorous fish endemic to East Asia, with high growth rate, nutritional value and economic value. In this study, a chromosome-level genome of yellow-cheek carp was generated by combining PacBio reads, Illumina reads and Hi-C data. The genome size is 827.63 Mb with a scaffold N50 size of 33.65 Mb, and 99.51% (823.61 Mb) of the assembled sequences were anchored to 24 pseudo-chromosomes. The genome is predicted to contain 24,153 protein-coding genes, with 95.54% having functional annotations. Repeat elements account for approximately 55.17% of the genomic landscape. The completeness of yellow-cheek carp genome assembly is highlighted by a BUSCO score of 98.4%. This genome will help us understand the genetic diversity of yellow-cheek carp and facilitate its conservation planning.

SNP and Structural Study of the Notch Superfamily Provides Insights and Novel Pharmacological Targets against the CADASIL Syndrome and Neurodegenerative Diseases

The evolutionary conserved Notch signaling pathway functions as a mediator of direct cell-cell communication between neighboring cells during development. Notch plays a crucial role in various fundamental biological processes in a wide range of tissues. Accordingly, the aberrant signaling of this pathway underlies multiple genetic pathologies such as developmental syndromes, congenital disorders, neurodegenerative diseases, and cancer. Over the last two decades, significant data have shown that the Notch signaling pathway displays a significant function in the mature brains of vertebrates and invertebrates beyond neuronal development and specification during embryonic development. Neuronal connection, synaptic plasticity, learning, and memory appear to be regulated by this pathway. Specific mutations in human Notch family proteins have been linked to several neurodegenerative diseases including Alzheimer's disease, CADASIL, and ischemic injury. Neurodegenerative diseases are incurable disorders of the central nervous system that cause the progressive degeneration and/or death of brain nerve cells, affecting both mental function and movement (ataxia). There is currently a lot of study being conducted to better understand the molecular mechanisms by which Notch plays an essential role in the mature brain. In this study, an in silico analysis of polymorphisms and mutations in human Notch family members that lead to neurodegenerative diseases was performed in order to investigate the correlations among Notch family proteins and neurodegenerative diseases. Particular emphasis was placed on the study of mutations in the Notch3 protein and the structure analysis of the mutant Notch3 protein that leads to the manifestation of the CADASIL syndrome in order to spot possible conserved mutations and interpret the effect of these mutations in the Notch3 protein structure. Conserved mutations of cysteine residues may be candidate pharmacological targets for the potential therapy of CADASIL syndrome.

A novel O-methyltransferase Cp4MP-OMT catalyses the final step in the biosynthesis of the volatile 1,4-dimethoxybenzene in pumpkin (Cucurbita pepo) flowers

BACKGROUND

Floral scents play a crucial role in attracting insect pollinators. Among the compounds attractive to pollinators is 1,4-dimethoxybenzene (1,4-DMB). It is a significant contributor to the scent profile of plants from various genera, including economically important Cucurbita species. Despite its importance, the biosynthetic pathway for the formation of 1,4-DMB was not elucidated so far.

RESULTS

In this study we showed the catalysis of 1,4-DMB in the presence of 4-methoxyphenol (4-MP) by protein extract from Styrian oil pumpkin (Cucurbita pepo) flowers. Based on this finding, we identified a novel O-methyltransferase gene, Cp4MP-OMT, whose expression is highly upregulated in the volatile-producing tissue of pumpkin flowers when compared to vegetative tissues. OMT activity was verified by purified recombinant Cp4MP-OMT, illustrating its ability to catalyse the methylation of 4-MP to 1,4-DMB in the presence of cofactor SAM (S-(5'-adenosyl)-L-methionine).

CONCLUSIONS

Cp4MP-OMT is a novel O-methyltransferase from C. pepo, responsible for the final step in the biosynthesis of the floral scent compound 1,4-DMB. Considering the significance of 1,4-DMB in attracting insects for pollination and in the further course fruit formation, enhanced understanding of its biosynthetic pathways holds great promise for both ecological insights and advancements in plant breeding initiatives.

In silico mutation of aromatic with aliphatic amino acid residues in Clostridium perfringens epsilon toxin (ETX) reduces its binding efficiency to Caprine Myelin and lymphocyte (MAL) protein receptors

Enterotoxaemia (ET) is a severe disease that affects domestic ruminants, including sheep and goats, and is caused by Clostridium perfringens type B and D strains. The disease is characterized by the production of Epsilon toxin (ETX), which has a significant impact on the farming industry due to its high lethality. The binding of ETX to the host cell receptor is crucial, but still poorly understood. Therefore, the structural features of goat Myelin and lymphocytic (MAL) protein were investigated and defined in this study. We induced the mutations in aromatic amino acid residues of ETX and substituted them with aliphatic residues at domains I and II. Subsequently, protein-protein interactions (PPI) were performed between ETX (wild)-MAL and ETX (mutated)-MAL protein predicting the domain sites of ETX structure. Further, molecular dynamics (MD) simulation studies were performed for both complexes to investigate the dynamic behavior of the proteins. The binding efficiency between 'ETX (wild)-MAL protein' and 'ETX (mutated)-MAL protein complex' interactions were compared and showed that the former had stronger interactions and binding efficiency due to the higher stability of the complex. The MD analysis showed destabilization and higher fluctuations in the PPI of the mutated heterodimeric ETX-MAL complex which is otherwise essential for its functional conformation. Such kind of interactions with mutated functional domains of ligands provided much-needed clarity in understanding the pre-pore complex formation of epsilon toxin with the MAL protein receptor of goats. The findings from this study would provide an impetus for designing a novel vaccine for Enterotoxaemia in goats.Communicated by Ramaswamy H. Sarma.

Designing a T-cell epitope-based vaccine using in silico approaches against the Sal k 1 allergen of Salsola kali plant

Allergens originated from Salsola kali (Russian thistle) pollen grains are one of the most important sources of aeroallergens causing pollinosis in desert and semi-desert regions. T-cell epitope-based vaccines (TEV) are more effective among different therapeutic approaches developed to alleviate allergic diseases. The physicochemical properties, and B as well as T cell epitopes of Sal k 1 (a major allergen of S. kali) were predicted using immunoinformatic tools. A TEV was constructed using the linkers EAAAK, GPGPG and the most suitable CD4+ T cell epitopes. RS04 adjuvant was added as a TLR4 agonist to the amino (N) and carboxyl (C) terminus of the TEV protein. The secondary and tertiary structures, solubility, allergenicity, toxicity, stability, physicochemical properties, docking with immune receptors, BLASTp against the human and microbiota proteomes, and in silico cloning of the designed TEV were assessed using immunoinformatic analyses. Two CD4+ T cell epitopes of Sal k1 that had high affinity with different alleles of MHC-II were selected and used in the TEV. The molecular docking of the TEV with HLADRB1, and TLR4 showed TEV strong interactions and stable binding pose to these receptors. Moreover, the codon optimized TEV sequence was cloned between NcoI and XhoI restriction sites of pET-28a(+) expression plasmid. The designed TEV can be used as a promising candidate in allergen-specific immunotherapy against S. kali. Nonetheless, effectiveness of this vaccine should be validated through immunological bioassays.

A systematic screen for co-option of transposable elements across the fungal kingdom

How novel protein functions are acquired is a central question in molecular biology. Key paths to novelty include gene duplications, recombination or horizontal acquisition. Transposable elements (TEs) are increasingly recognized as a major source of novel domain-encoding sequences. However, the impact of TE coding sequences on the evolution of the proteome remains understudied. Here, we analyzed 1237 genomes spanning the phylogenetic breadth of the fungal kingdom. We scanned proteomes for evidence of co-occurrence of TE-derived domains along with other conventional protein functional domains. We detected more than 13,000 predicted proteins containing potentially TE-derived domain, of which 825 were identified in more than five genomes, indicating that many host-TE fusions may have persisted over long evolutionary time scales. We used the phylogenetic context to identify the origin and retention of individual TE-derived domains. The most common TE-derived domains are helicases derived from Academ, Kolobok or Helitron. We found putative TE co-options at a higher rate in genomes of the Saccharomycotina, providing an unexpected source of protein novelty in these generally TE depleted genomes. We investigated in detail a candidate host-TE fusion with a heterochromatic transcriptional silencing function that may play a role in TE and gene regulation in ascomycetes. The affected gene underwent multiple full or partial losses within the phylum. Overall, our work establishes a kingdom-wide view of putative host-TE fusions and facilitates systematic investigations of candidate fusion proteins.

Modeling methyl-sensitive transcription factor motifs with an expanded epigenetic alphabet

BACKGROUND

Transcription factors bind DNA in specific sequence contexts. In addition to distinguishing one nucleobase from another, some transcription factors can distinguish between unmodified and modified bases. Current models of transcription factor binding tend not to take DNA modifications into account, while the recent few that do often have limitations. This makes a comprehensive and accurate profiling of transcription factor affinities difficult.

RESULTS

Here, we develop methods to identify transcription factor binding sites in modified DNA. Our models expand the standard A/C/G/T DNA alphabet to include cytosine modifications. We develop Cytomod to create modified genomic sequences and we also enhance the MEME Suite, adding the capacity to handle custom alphabets. We adapt the well-established position weight matrix (PWM) model of transcription factor binding affinity to this expanded DNA alphabet. Using these methods, we identify modification-sensitive transcription factor binding motifs. We confirm established binding preferences, such as the preference of ZFP57 and C/EBPβ for methylated motifs and the preference of c-Myc for unmethylated E-box motifs.

CONCLUSIONS

Using known binding preferences to tune model parameters, we discover novel modified motifs for a wide array of transcription factors. Finally, we validate our binding preference predictions for OCT4 using cleavage under targets and release using nuclease (CUT&RUN) experiments across conventional, methylation-, and hydroxymethylation-enriched sequences. Our approach readily extends to other DNA modifications. As more genome-wide single-base resolution modification data becomes available, we expect that our method will yield insights into altered transcription factor binding affinities across many different modifications.

The Digital World of Cytogenetic and Cytogenomic Web Resources

The dynamic growth of technological capabilities at the cellular and molecular level has led to a rapid increase in the amount of data on the genes and genomes of organisms. In order to store, access, compare, validate, classify, and understand the massive data generated by different researchers, and to promote effective communication among research communities, various genome and cytogenetic online databases have been established. These data platforms/resources are essential not only for computational analyses and theoretical syntheses but also for helping researchers select future research topics and prioritize molecular targets. Furthermore, they are valuable for identifying shared recurrent genomic patterns related to human diseases and for avoiding unnecessary duplications among different researchers. The website interface, menu, graphics, animations, text layout, and data from databases are displayed by a front end on the screen of a monitor or smartphone. A database front-end refers to the user interface or application that enables accessing tabular, structured, or raw data stored in the database. The Internet makes it possible to reach a greater number of users around the world and gives them quick access to information stored in databases. The number of ways of presenting this data by front-ends increases as well. This requires unifying the ways of operating and presenting information by front-ends and ensuring contextual switching between front-ends of different databases. This chapter aims to present selected cytogenetic and cytogenomic Internet resources in terms of obtaining the needed information and to indicate how to increase the efficiency of access to stored information. Through a brief introduction of these databases and by providing examples of their usage in cytogenetic analyses, we aim to bridge the gap between cytogenetics and molecular genomics by encouraging their utilization.

Bioinformatic Workflows for Metaproteomics

The strong influence of microbiomes on areas such as ecology and human health has become widely recognized in the past years. Accordingly, various techniques for the investigation of the composition and function of microbial community samples have been developed. Metaproteomics, the comprehensive analysis of the proteins from microbial communities, allows for the investigation of not only the taxonomy but also the functional and quantitative composition of microbiome samples. Due to the complexity of the investigated communities, methods developed for single organism proteomics cannot be readily applied to metaproteomic samples. For this purpose, methods specifically tailored to metaproteomics are required. In this work, a detailed overview of current bioinformatic solutions and protocols in metaproteomics is given. After an introduction to the proteomic database search, the metaproteomic post-processing steps are explained in detail. Ten specific bioinformatic software solutions are focused on, covering various steps including database-driven identification and quantification as well as taxonomic and functional assignment.

Chromosomal scale assembly reveals localized structural variants in avian caecal coccidian parasite Eimeria tenella

Eimeria tenella is a major cause of caecal coccidiosis in commercial poultry chickens worldwide. Here, we report chromosomal scale assembly of Eimeria tenella strain APU2, a strain isolated from commercial broiler chickens in the U.S. We obtained 100× sequencing Oxford Nanopore Technology (ONT) and more than 800× Coverage of Illumina Next-Seq. We created the assembly using the hybrid approach implemented in MaSuRCA, achieving a contiguous 51.34 Mb chromosomal-scale scaffolding enabling identification of structural variations. The AUGUSTUS pipeline predicted 8060 genes, and BUSCO deemed the genomes 99% complete; 6278 (78%) genes were annotated with Pfam domains, and 1395 genes were assigned GO-terms. Comparing E. tenella strains (APU2, US isolate and Houghton, UK isolate) derived Houghton strain of E. tenella revealed 62,905 high stringency differences, of which 45,322 are single nucleotide polymorphisms (SNPs) (0.088%). The rate of transitions/transversions among the SNPs are 1.63 ts/tv. The strains possess conserved gene order but have profound sequence heterogeneity in a several chromosomal segments (chr 2, 11 and 15). Genic and intergenic variation in defined gene families was evaluated between the two strains to possibly identify sequences under selection. The average genic nucleotide diversity of 2.8 with average 2 kb gene length (0.145%) at genic level. We examined population structure using available E. tenella sequences in NCBI, revealing that the two E. tenella isolates from the U.S. (E. tenella APU2 and Wisconsin, "ERR296879") share a common maternal inheritance with the E. tenella Houghton. Our chromosomal level assembly promotes insight into Eimeria biology and evolution, hastening drug discovery and vaccine development.

Genome-wide identification and expression analysis of the JMJ-C gene family in melon (Cucumis melo L.) reveals their potential role in fruit development

BACKGROUND

Proteins with the jumonji (JMJ)-C domain belong to the histone demethylase family and contribute to reverse histone methylation. Although JMJ-C family genes have an essential role in regulating plant growth and development, the characterization of the JMJ-C family genes in melon has not been uncovered.

RESULTS

In this study, a total of 17 JMJ-C proteins were identified in melon (Cucumis melo L.). CmJMJs were categorized into five subfamilies based on the specific conserved domain: KDM4/JHDM3, KDM5/JARID1, JMJD6, KDM3/JHDM2, and JMJ-C domain-only. The chromosome localization analyses showed that 17 CmJMJs were distributed on nine chromosomes. Cis-acting element analyses of the 17 CmJMJ genes showed numerous hormone, light, and stress response elements distributed in the promoter region. Covariance analysis revealed one pair of replicated fragments (CmJMJ3a and CmJMJ3b) in 17 CmJMJ genes. We investigated the expression profile of 17 CmJMJ genes in different lateral organs and four developmental stages of fruit by RNA-seq transcriptome analysis and RT-qPCR. The results revealed that most CmJMJ genes were prominently expressed in female flowers, ovaries, and developing fruits, suggesting their active role in melon fruit development. Subcellular localization showed that the fruit-related CmJMJ5a protein is specifically localized in the cell nucleus.

CONCLUSIONS

This study provides a comprehensive understanding of the gene structure, classification, and evolution of JMJ-C in melon and supports the clarification of the JMJ-C functions in further research.

Comprehensive Evolution and Expression anaLysis of PHOSPHATE 1 Gene Family in Allotetraploid Brassica napus and Its Diploid Ancestors

The members of PHOSPHATE 1 (PHO1) family play important roles in plant phosphate (Pi) transport and adaptation to Pi deficiency. The functions of PHO1 family proteins have been reported in several plant species, with the exception of Brassica species. Here, we identified 23, 23, and 44 putative PHO1 family genes in Brassica rapa, Brassica oleracea, and Brassica napus by whole genome analysis, respectively. The phylogenetic analysis divided PHO1 family proteins into eight groups, which represented the orthologous relationships among PHO1 members. The gene structure and the conserved motif analysis indicated that the most PHO1 family genes had similar gene structures and the PHO1 proteins shared mutual conserved motifs. The chromosome distribution analysis showed that the majority of BnPHO1 family genes distributed analogously at chromosomes with BrPHO1 and BoPHO1 family genes. The data showed that PHO1 family genes were highly conserved during evolution from diploid to tetraploid. Furthermore, the expression analysis showed that PHO1 family genes had different expression patterns in plant tissues, suggesting the diversity of gene functions in Brassica species. Meanwhile, the expression analysis also revealed that some PHO1 family genes were significantly responsive to Pi deficiency, suggesting that PHO1 family genes play critical roles in Pi uptake and homeostasis under low Pi stress. Altogether, the characteristics of PHO1 family genes provide a reliable groundwork for further dissecting their functions in Brassica species.

Exploring complexity of class-A Beta-lactamase family using physiochemical-based multiplex networks

The Beta-lactamase protein family is vital in countering Beta-lactam antibiotics, a widely used antimicrobial. To enhance our understanding of this family, we adopted a novel approach employing a multiplex network representation of its multiple sequence alignment. Each network layer, derived from the physiochemical properties of amino acids, unveils distinct insights into the intricate interactions among nodes, thereby enabling the identification of key motifs. Nodes with identical property signs tend to aggregate, providing evidence of the presence of consequential functional and evolutionary constraints shaping the Beta-lactamase family. We further investigate the distribution of evolutionary links across various layers. We observe that polarity manifests the highest number of unique links at lower thresholds, followed by hydrophobicity and polarizability, wherein hydrophobicity exerts dominance at higher thresholds. Further, the combinations of polarizability and volume, exhibit multiple simultaneous connections at all thresholds. The combination of hydrophobicity, polarizability, and volume uncovers shared links exclusive to these layers, implying substantial evolutionary impacts that may have functional or structural implications. By assessing the multi-degree of nodes, we unveil the hierarchical influence of properties at each position, identifying crucial properties responsible for the protein's functionality and providing valuable insights into potential targets for modulating enzymatic activity.

Evolutionary Insights into the Relationship of Frogs, Salamanders, and Caecilians and Their Adaptive Traits, with an Emphasis on Salamander Regeneration and Longevity

The extant amphibians have developed uncanny abilities to adapt to their environment. I compared the genes of amphibians to those of other vertebrates to investigate the genetic changes underlying their unique traits, especially salamanders' regeneration and longevity. Using the well-supported Batrachia tree, I found that salamander genomes have undergone accelerated adaptive evolution, especially for development-related genes. The group-based comparison showed that several genes are under positive selection, rapid evolution, and unexpected parallel evolution with traits shared by distantly related species, such as the tail-regenerative lizard and the longer-lived naked mole rat. The genes, such as EEF1E1, PAFAH1B1, and OGFR, may be involved in salamander regeneration, as they are involved in the apoptotic process, blastema formation, and cell proliferation, respectively. The genes PCNA and SIRT1 may be involved in extending lifespan, as they are involved in DNA repair and histone modification, respectively. Some genes, such as PCNA and OGFR, have dual roles in regeneration and aging, which suggests that these two processes are interconnected. My experiment validated the time course differential expression pattern of SERPINI1 and OGFR, two genes that have evolved in parallel in salamanders and lizards during the regeneration process of salamander limbs. In addition, I found several candidate genes responsible for frogs' frequent vocalization and caecilians' degenerative vision. This study provides much-needed insights into the processes of regeneration and aging, and the discovery of the critical genes paves the way for further functional analysis, which could open up new avenues for exploiting the genetic potential of humans and improving human well-being.

Domain-PFP allows protein function prediction using function-aware domain embedding representations

Domains are functional and structural units of proteins that govern various biological functions performed by the proteins. Therefore, the characterization of domains in a protein can serve as a proper functional representation of proteins. Here, we employ a self-supervised protocol to derive functionally consistent representations for domains by learning domain-Gene Ontology (GO) co-occurrences and associations. The domain embeddings we constructed turned out to be effective in performing actual function prediction tasks. Extensive evaluations showed that protein representations using the domain embeddings are superior to those of large-scale protein language models in GO prediction tasks. Moreover, the new function prediction method built on the domain embeddings, named Domain-PFP, substantially outperformed the state-of-the-art function predictors. Additionally, Domain-PFP demonstrated competitive performance in the CAFA3 evaluation, achieving overall the best performance among the top teams that participated in the assessment.

Identification of potential novel therapeutic drug target against Elizabethkingia anophelis by integrative pan and subtractive genomic analysis: An in silico approach

Elizabethkingia anophelis is a human pathogen responsible for severe nosocomial infections in neonates and immunocompromised patients. The significantly higher mortality rate from E. anophelis infections and the lack of available regimens highlight the critical need to explore novel drug targets. The current study investigated effective novel drug targets by employing a comprehensive in silico subtractive genomic approach integrated with pangenomic analysis of E. anophelis strains. A total of 2809 core genomic proteins were found by pangenomic analysis of non-paralogous proteins. Subsequently, 156 pathogen-specific, 442 choke point, 202 virulence factor, 53 antibiotic resistant and 119 host-pathogen interacting proteins were identified in E. anophelis. By subtractive genomic approach, at first 791 proteins were found to be indispensable for the survival of E. anophelis. 558 and 315 proteins were detected as non-homologous to human and gut microflora respectively. Following that 245 cytoplasmic, 245 novel, and 23 broad-spectrum targets were selected and finally four proteins were considered as potential therapeutic targets of E. anophelis based on highest degree score in PPI network. Among those, three proteins were subjected to molecular docking and subsequent MD simulation as one protein did not contain a plausible binding pocket with sufficient surface area and volume. All the complexes were found to be stable and compact in 100 ns molecular dynamics simulation studies as measured by RMSD, RMSF, and Rg. These three short-listed targets identified in this study may lead to the development of novel antimicrobials capable of curing infections and pave the way to prevent and control the disease progression caused by the deadly agent E. anophelis.

Long-term biodegradation of crude oil in high-arctic backshore sediments: The Baffin Island Oil Spill (BIOS) after nearly four decades

With an on-going disproportional warming of the Arctic Ocean and the reduction of the sea ice cover, the risk of an accidental oil spill from ships or future oil exploration is increasing. It is hence important to know how crude oil weathers in this environment and what factors affect oil biodegradation in the Arctic. However, this topic is currently poorly studied. In the 1980s, the Baffin Island Oil Spill (BIOS) project carried out a series of simulated oil spills in the backshore zone of beaches located on Baffin Island in the Canadian High Arctic. In this study two BIOS sites were re-visited, offering the unique opportunity to study the long-term weathering of crude oil under Arctic conditions. Here we show that residual oil remains present at these sites even after almost four decades since the original oiling. Oil at both BIOS sites appears to have attenuated very slowly with estimated loss rates of 1.8-2.7% per year. The presence of residual oil continues to significantly affect sediment microbial communities at the sites as manifested by a significantly decreased diversity, differences in the abundance of microorganisms and an enrichment of putative oil-degrading bacteria in oiled sediments. Reconstructed genomes of putative oil degraders suggest that only a subset is specifically adapted for growth under psychrothermic conditions, further reducing the time for biodegradation during the already short Arctic summers. Altogether, this study shows that crude oil spilled in the Arctic can persist and significantly affect the Arctic ecosystem for a long time, in the order of several decades.

Genomes of cultivated and wild Capsicum species provide insights into pepper domestication and population differentiation

Pepper (Capsicum spp.) is one of the earliest cultivated crops and includes five domesticated species, C. annuum var. annuum, C. chinense, C. frutescens, C. baccatum var. pendulum and C. pubescens. Here, we report a pepper graph pan-genome and a genome variation map of 500 accessions from the five domesticated Capsicum species and close wild relatives. We identify highly differentiated genomic regions among the domesticated peppers that underlie their natural variations in flowering time, characteristic flavors, and unique resistances to biotic and abiotic stresses. Domestication sweeps detected in C. annuum var. annuum and C. baccatum var. pendulum are mostly different, and the common domestication traits, including fruit size, shape and pungency, are achieved mainly through the selection of distinct genomic regions between these two cultivated species. Introgressions from C. baccatum into C. chinense and C. frutescens are detected, including those providing genetic sources for various biotic and abiotic stress tolerances.

Domain-PFP: Protein Function Prediction Using Function-Aware Domain Embedding Representations

Domains are functional and structural units of proteins that govern various biological functions performed by the proteins. Therefore, the characterization of domains in a protein can serve as a proper functional representation of proteins. Here, we employ a self-supervised protocol to derive functionally consistent representations for domains by learning domain-Gene Ontology (GO) co-occurrences and associations. The domain embeddings we constructed turned out to be effective in performing actual function prediction tasks. Extensive evaluations showed that protein representations using the domain embeddings are superior to those of large-scale protein language models in GO prediction tasks. Moreover, the new function prediction method built on the domain embeddings, named Domain-PFP, significantly outperformed the state-of-the-art function predictors. Additionally, Domain-PFP demonstrated competitive performance in the CAFA3 evaluation, achieving overall the best performance among the top teams that participated in the assessment.

Potential Role of Lysine Acetylation and Autophagy in Brown Film Formation and Postripening of Lentinula edodes Mycelium

Lentinula edodes is one of the most widely cultivated edible mushrooms in the world. When cultivated in sawdust, the surface mycelium of L. edodes needs a long postripening stage wherein it forms a brown film (BF) by secreting and accumulating pigments. BF formation is critical for the high quality and yield of fruiting bodies. Protein lysine acetylation (KAC) is an important post-translational modification that regulates growth and development. Previous studies have shown that deacetylase levels are significantly increased during BF formation in the postripening stage of L. edodes. The aim of this study was to assess the role of protein acetylation during BF formation. To this end, we compared the acetylome of L. edodes mycelia before and after BF formation using anti-acetyl antibody-based label-free quantitative proteomics. We identified 5,613 acetylation sites in 1,991 proteins, and quantitative information was available for 4,848 of these sites in 1,815 proteins. Comparative acetylome analysis showed that the modification of 699 sites increased and that of 562 sites decreased during BF formation. Bioinformatics analysis of the differentially acetylated proteins showed significant enrichment in the tricarboxylic acid (TCA) cycle and proteasome pathways. Furthermore, functional assays showed that BF formation is associated with significant changes in the activities of proteasome, citrate synthase, and isocitrate dehydrogenase. Consistent with this hypothesis, the lysine deacetylase inhibitor trichostatin (TSA) delayed autophagy and BF formation in L. edodes. Taken together, KAC and autophagy play important roles in the mycelial BF formation and postripening stage of L. edodes. IMPORTANCE Mycelial BF formation and postripening of L. edodes affects the quality and quantity of its edible fruiting bodies. In this study, we explored the role of protein KAC in this biological process, with the aim of optimizing the cultivation and yield of L. edodes.

The Proteome and Lipidome of Extracellular Vesicles from Haemonchus contortus to Underpin Explorations of Host-Parasite Cross-Talk

Many parasitic worms have a major adverse impact on human and animal populations worldwide due to the chronicity of their infections. There is a growing body of evidence indicating that extracellular vesicles (EVs) are intimately involved in modulating (suppressing) inflammatory/immune host responses and parasitism. As one of the most pathogenic nematodes of livestock animals, Haemonchus contortus is an ideal model system for EV exploration. Here, employing a multi-step enrichment process (in vitro culture, followed by ultracentrifugation, size exclusion and filtration), we enriched EVs from H. contortus and undertook the first comprehensive (qualitative and quantitative) multi-omic investigation of EV proteins and lipids using advanced liquid chromatography-mass spectrometry and informatics methods. We identified and quantified 561 proteins and 446 lipids in EVs and compared these molecules with those of adult worms. We identified unique molecules in EVs, such as proteins linked to lipid transportation and lipid species (i.e., sphingolipids) associated with signalling, indicating the involvement of these molecules in parasite-host cross-talk. This work provides a solid starting point to explore the functional roles of EV-specific proteins and lipids in modulating parasite-host cross-talk, and the prospect of finding ways of disrupting or interrupting this relationship to suppress or eliminate parasite infection.

Biotechnological Tools to Elucidate the Mechanism of Plant and Nematode Interactions

Plant-parasitic nematodes (PPNs) pose a threat to global food security in both the developed and developing worlds. PPNs cause crop losses worth a total of more than USD 150 billion worldwide. The sedentary root-knot nematodes (RKNs) also cause severe damage to various agricultural crops and establish compatible relationships with a broad range of host plants. This review aims to provide a broad overview of the strategies used to identify the morpho-physiological and molecular events that occur during RKN parasitism. It describes the most current developments in the transcriptomic, proteomic, and metabolomic strategies of nematodes, which are important for understanding compatible interactions of plants and nematodes, and several strategies for enhancing plant resistance against RKNs. We will highlight recent rapid advances in molecular strategies, such as gene-silencing technologies, RNA interference (RNAi), and small interfering RNA (siRNA) effector proteins, that are leading to considerable progress in understanding the mechanism of plant-nematode interactions. We also take into account genetic engineering strategies, such as targeted genome editing techniques, the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR associated protein 9 (Cas9) (CRISPR/Cas-9) system, and quantitative trait loci (QTL), to enhance the resistance of plants against nematodes.

Host transcriptomic plasticity and photosymbiotic fidelity underpin Pocillopora acclimatization across thermal regimes in the Pacific Ocean

Heat waves are causing declines in coral reefs globally. Coral thermal responses depend on multiple, interacting drivers, such as past thermal exposure, endosymbiont community composition, and host genotype. This makes the understanding of their relative roles in adaptive and/or plastic responses crucial for anticipating impacts of future warming. Here, we extracted DNA and RNA from 102 Pocillopora colonies collected from 32 sites on 11 islands across the Pacific Ocean to characterize host-photosymbiont fidelity and to investigate patterns of gene expression across a historical thermal gradient. We report high host-photosymbiont fidelity and show that coral and microalgal gene expression respond to different drivers. Differences in photosymbiotic association had only weak impacts on host gene expression, which was more strongly correlated with the historical thermal environment, whereas, photosymbiont gene expression was largely determined by microalgal lineage. Overall, our results reveal a three-tiered strategy of thermal acclimatization in Pocillopora underpinned by host-photosymbiont specificity, host transcriptomic plasticity, and differential photosymbiotic association under extreme warming.

A Deep Learning Framework for Predicting Protein Functions With Co-Occurrence of GO Terms

The understanding of protein functions is critical to many biological problems such as the development of new drugs and new crops. To reduce the huge gap between the increase of protein sequences and annotations of protein functions, many methods have been proposed to deal with this problem. These methods use Gene Ontology (GO) to classify the functions of proteins and consider one GO term as a class label. However, they ignore the co-occurrence of GO terms that is helpful for protein function prediction. We propose a new deep learning model, named DeepPFP-CO, which uses Graph Convolutional Network (GCN) to explore and capture the co-occurrence of GO terms to improve the protein function prediction performance. In this way, we can further deduce the protein functions by fusing the predicted propensity of the center function and its co-occurrence functions. We use Fmax and AUPR to evaluate the performance of DeepPFP-CO and compare DeepPFP-CO with state-of-the-art methods such as DeepGOPlus and DeepGOA. The computational results show that DeepPFP-CO outperforms DeepGOPlus and other methods. Moreover, we further analyze our model at the protein level. The results have demonstrated that DeepPFP-CO improves the performance of protein function prediction. DeepPFP-CO is available at https://csuligroup.com/DeepPFP/.

Chromosome-level Genome Assembly and Sex-specific Differential Transcriptome of the White-backed Planthopper, Sogatella furcifera

Background

The white-backed planthopper (WBPH), Sogatella furcifera, causes great damage to many crops (mainly rice) by direct feeding or transmitting plant viruses. The previous genome assembly was generated by second-generation sequencing technologies, with a contig N50 of only 51.5 kb, and contained a lot of heterozygous sequences.

Methods

We utilized third-generation sequencing technologies and Hi-C data to generate a high-quality chromosome-level assembly. We also provide a large amount of transcriptome data for full-length transcriptome analysis and gender differential expression analysis.

Results

The final assembly comprised 56.38 Mb, with a contig N50 of 2.20 Mb and a scaffold N50 of 45.25 Mb. Fourteen autosomes and one X chromosome were identified. More than 99.5% of the assembled bases located on the 15 chromosomes. 95.9% of the complete BUSCO Hemiptera genes were detected in the final assembly and 16,880 genes were annotated. 722 genes were relatively highly expressed in males, while 60 in the females.

Conclusion

The integrated genome, definite sex chromosomes, comprehensive transcriptome profiles, high efficiency of RNA interference and short life cycle substantially made WBPH an efficient research object for functional genomics.

Odorant-Binding and Chemosensory Proteins in Anthonomus eugenii (Coleoptera: Curculionidae) and Their Tissue Expression

The pepper weevil Anthonomus eugenii is one of the most damaging pests to the pepper crop. To offer alternative management strategies to insecticides, several studies have identified the semiochemicals that are involved in the pepper weevil's aggregation and mating behavior; however, there is no information on its perireceptor molecular mechanism, to date. In this study, bioinformatics tools were used to functionally annotate and characterize the A. eugenii head transcriptome and their probable coding proteins. We identified twenty-two transcripts belonging to families related to chemosensory processes, seventeen corresponding to odorant-binding proteins (OBP), and six to chemosensory proteins (CSP). All results matched with closely related Coleoptera: Curculionidae homologous proteins. Likewise, twelve OBP and three CSP transcripts were experimentally characterized by RT-PCR in different female and male tissues. The results by sex and tissue display the different expression patterns of the AeugOBPs and AeugCSPs; some are present in both sexes and all tissues, while others show expressions with higher specificity, which suggests diverse physiological functions in addition to chemo-detection. This study provides information to support the understanding of odor perception in the pepper weevil.

The genome of a Far Eastern isolate of Diaporthe caulivora, a soybean fungal pathogen

Diaporthe caulivora is an economically important fungal pathogen and a causal agent of soybean stem canker and seed decay. Here, the genome of a Russian Far Eastern isolate of D. caulivora was sequenced, assembled, and announced. Assembly quality was enough for advanced annotation, including prediction of potential disease-related genes encoding virulence factors and molecular determinants contributing to pathogen-host selection, interactions, and adaptation. Comparative analysis of 15 Diaporthe species was conducted regarding general genome properties, collinearity, and proteomes, and included detailed investigation of interspersed repeats. A notable feature of this analysis is a high recombinant variability of Diaporthe genomes, determined by the number and distribution of interspersed repeats, which also proved to be responsible for the diversity of GC content and genome size. This variability is assumed the main determinant of the divergence of Diaporthe genomes. A Bayesian multi-gene phylogeny was inferred for the 15 Diaporthe species on the basis of twenty thousand polymorphic sites of > 100 orthologous genes using independently adjusted evolutionary models. This allowed for the most accurate determination of evolutionary relationships and species boundaries for effective reporting about these plant pathogens. The evidence, obtained by different genome analysis techniques, implies the host-independent evolution of Diaporthe species. KEY POINTS: • The genome of a Far Eastern isolate of D. caulivora was announced. • A high degree of recombinant variability determines genomic divergence in Diaporthe genus. • The multi-gene phylogeny implies host-independent evolution of Diaporthe species.

Genome Assembly and Comparative Analysis of the Egg Parasitoid Wasp Trichogramma dendrolimi Shed Light on the Composition and Evolution of Olfactory Receptors and Venoms

Trichogramma dendrolimi is one of the most successfully industrialized Trichogramma species used to control agricultural and forestry pests in China. However, the molecular mechanisms underlying its host recognition and parasitism remain largely unknown, partially due to the limited genome information of this parasitoid wasp. Here, we present a high-quality de novo assembly of T. dendrolimi through a combination of Illumina and PacBio sequencing technologies. The final assembly had a length of 215.2 Mb and contains 316 scaffolds with a scaffold N50 size of 1.41 Mb. Repetitive sequences with a length of 63.4 Mb and 12,785 protein-coding genes were identified. Significantly expanded gene families were identified to be involved in the development and regulatory processes, while remarkably contracted gene families were involved in the transport processes in T. dendrolimi. The olfactory and venom-associated genes were identified in T. dendrolimi and 24 other hymenopteran species, using uniform methods combining BLAST and HMM profiling. The identified venom genes of T. dendrolimi were enriched in antioxidant activity, tricarboxylic acid cycle, response to oxidative stress and cell redox homeostasis. Our study provides an important resource for comparative genomics and functional studies to interpret the molecular mechanisms underlying host recognition and parasitism of Trichogramma species.

Systematic analysis reveals novel insight into the molecular determinants of function, diversity and evolution of sweet taste receptors T1R2/T1R3 in primates

Sweet taste is a primary sensation for the preference and adaption of primates to diet, which is crucial for their survival and fitness. It is clear now that the sweet perception is mediated by a G protein-coupled receptor (GPCR)-sweet taste receptor T1R2/T1R3, and many behavioral or physiological experiments have described the diverse sweet taste sensitivities in primates. However, the structure-function relationship of T1R2s/T1R3s in primates, especially the molecular basis for their species-dependent sweet taste, has not been well understood until now. In this study, we performed a comprehensive sequence, structural and functional analysis of sweet taste receptors in primates to elucidate the molecular determinants mediating their species-dependent sweet taste recognition. Our results reveal distinct taxonomic distribution and significant characteristics (interaction, coevolution and epistasis) of specific key function-related residues, which could partly account for the previously reported behavioral results of taste perception in primates. Moreover, the prosimians Lemuriformes species, which were reported to have no sensitivity to aspartame, could be proposed to be aspartame tasters based on the present analysis. Collectively, our study provides new insights and promotes a better understanding for the diversity, function and evolution of sweet taste receptors in primates.

DIAMOND + MEGAN Microbiome Analysis

Metagenomics is the study of microbiomes using DNA sequencing technologies. Basic computational tasks are to determine the taxonomic composition (who is out there?), the functional composition (what can they do?), and also to correlate changes of composition to changes in external parameters (how do they compare?). One approach to address these issues is to first align all sequences against a protein reference database such as NCBI-nr and to then perform taxonomic and functional binning of all sequences based on their alignments. The resulting classifications can then be interactively analyzed and compared. Here we illustrate how to pursue this approach using the DIAMOND+MEGAN pipeline, on two different publicly available datasets, one containing short-read samples and other containing long-read samples.

Altered ureido protein modification profiles in seminal plasma extracellular vesicles of non-normozoospermic men

INTRODUCTION

Extracellular vesicles (EVs) have been recognized as key players in numerous physiological functions. These vesicles alter their compositions attuned to the health and disease states of the organism. In men, significant changes in the proteomic composition(s) of seminal plasma EVs (sEVs) have already been found to be related to infertility.

METHODS

Methods: In this study, we analyze the posttranslational configuration of sEV proteomes from normozoospermic (NZ) men and non-normozoospermic (non-NZ) men diagnosed with teratozoospermia and/or asthenozoospermia by unbiased, discovery-driven proteomics and advanced bioinformatics, specifically focusing on citrulline (Cit) and homocitrulline (hCit) posttranscriptional residues, both considered product of ureido protein modifications.

RESULTS AND DISCUSSION

Significant increase in the proteome-wide cumulative presence of hCit together with downregulation of Cit in specific proteins related to decisive molecular functions have been encountered in sEVs of non-NZ subjects. These findings identify novel culprits with a higher chance of affecting fundamental aspects of sperm functional quality and define potential specific diagnostic and prognostic non-invasive markers for male infertility.

Wolbachia manipulates reproduction of spider mites by influencing herbivore salivary proteins

BACKGROUND

The endosymbiont Wolbachia is known for manipulating host reproduction. Wolbachia also can affect host fitness by mediating interactions between plant and herbivores. However, it remains unclear whether saliva proteins are involved in this process.

RESULTS

We found that Wolbachia infection decreased the number of deposited eggs but increased the egg hatching rate in the spider mite Tetranychus urticae Koch (Acari: Tetranychidae), a cosmopolitan pest that infects >1000 species of plants. Transcriptomic and proteomic analyses revealed that Wolbachia-infected mites upregulated the gene expression levels of many T. urticae salivary proteins including a cluster of Tetranychidae-specific, functionally uncharacterized SHOT1s (secreted host-responsive proteins of Tetranychidae). The SHOT1 genes were expressed more in the feeding stages (nymphs and adults) of mites than in eggs and highly enriched in the proterosomas. RNA interference experiments showed that knockdown of SHOT1s significantly decreased Wolbachia density, increased the number of deposited eggs and decreased the egg hatching rate.

CONCLUSION

Together, these results indicate that SHOT1s are positively correlated with Wolbachia density and account for Wolbachia-mediated phenotypes. Our results provide new evidence that herbivore salivary proteins are related to Wolbachia-mediated manipulations of host performance on plants. © 2022 Society of Chemical Industry.

Systems biology and in silico-based analysis of PCOS revealed the risk of metabolic disorders

Background

Polycystic ovarian syndrome (PCOS) is a common condition of hyperandrogenism, chronic ovulation, and polycystic ovaries in females during the reproduction and maturation of the ovum. Although PCOS has been associated with metabolic disorders, including type 2 diabetes (T2D), obesity (OBE), and cardiovascular disease (CVD), Causal connection and molecular features are still unknown.

Purpose

Therefore, we investigated the shared common differentially expressed genes (DEGs), pathways, and networks of associated proteins in PCOS and metabolic diseases with therapeutic intervention.

Methods

We have used a bioinformatics pipeline to analyze transcriptome data for the polycystic ovarian syndrome (PCOS), type 2 diabetes (T2D), obesity (OBE), and cardiovascular diseases (CVD) in female patients. Then we employed gene-disease association network, gene ontology (GO) and signaling pathway analysis, selection of hub genes from protein-protein interaction (PPI) network, molecular docking, and gold benchmarking approach to screen potential hub proteins.

Result

We discovered 2225 DEGs in PCOS patients relative to healthy controls and 34, 91, and 205 significant DEGs with T2D, Obesity, and CVD, respectively. Gene Ontology analysis revealed several significant shared and metabolic pathways from signaling pathway analysis. Furthermore, we identified ten potential hub proteins from PPI analysis that may serve as a therapeutic intervention in the future. Finally, we targeted one significant hub protein, IGF2R (PDB ID: 2V5O), out of ten hub proteins based on the Maximal clique centrality (MCC) algorithm and literature review for molecular docking study. Enzastaurin (-12.5), Kaempferol (-9.1), Quercetin (-9.0), and Coumestrol (-8.9) kcal/mol showed higher binding affinity in the molecular docking approach than 19 drug compounds. We have also found that the selected four compounds displayed favorable ADMET properties compared to the native ligand.

Conclusion

Our in-silico research findings identified a shared molecular etiology between PCOS and metabolic diseases that may suggest new therapeutic targets and warrants future experimental validation of the key targets.

Drug repositioning based on heterogeneous networks and variational graph autoencoders

Predicting new therapeutic effects (drug repositioning) of existing drugs plays an important role in drug development. However, traditional wet experimental prediction methods are usually time-consuming and costly. The emergence of more and more artificial intelligence-based drug repositioning methods in the past 2 years has facilitated drug development. In this study we propose a drug repositioning method, VGAEDR, based on a heterogeneous network of multiple drug attributes and a variational graph autoencoder. First, a drug-disease heterogeneous network is established based on three drug attributes, disease semantic information, and known drug-disease associations. Second, low-dimensional feature representations for heterogeneous networks are learned through a variational graph autoencoder module and a multi-layer convolutional module. Finally, the feature representation is fed to a fully connected layer and a Softmax layer to predict new drug-disease associations. Comparative experiments with other baseline methods on three datasets demonstrate the excellent performance of VGAEDR. In the case study, we predicted the top 10 possible anti-COVID-19 drugs on the existing drug and disease data, and six of them were verified by other literatures.

Multiple relaxases contribute to the horizontal transfer of the virulence plasmids from the tumorigenic bacterium Pseudomonas syringae pv. savastanoi NCPPB 3335

Pseudomonas syringae pv. savastanoi NCPPB 3335 is the causal agent of olive knot disease and contains three virulence plasmids: pPsv48A (pA), 80 kb; pPsv48B (pB), 45 kb, and pPsv48C (pC), 42 kb. Here we show that pB contains a complete MPF_T (previously type IVA secretion system) and a functional origin of conjugational transfer adjacent to a relaxase of the MOB_P family; pC also contains a functional oriT-MOB_P array, whereas pA contains an incomplete MPF_I (previously type IVB secretion system), but not a recognizable oriT. Plasmid transfer occurred on solid and in liquid media, and on leaf surfaces of a non-host plant (Phaseolus vulgaris) with high (pB) or moderate frequency (pC); pA was transferred only occasionally after cointegration with pB. We found three plasmid-borne and three chromosomal relaxase genes, although the chromosomal relaxases did not contribute to plasmid dissemination. The MOB_P relaxase genes of pB and pC were functionally interchangeable, although with differing efficiencies. We also identified a functional MOB_Q mobilization region in pC, which could only mobilize this plasmid. Plasmid pB could be efficiently transferred to strains of six phylogroups of P. syringae sensu lato, whereas pC could only be mobilized to two strains of phylogroup 3 (genomospecies 2). In two of the recipient strains, pB was stably maintained after 21 subcultures in liquid medium. The carriage of several relaxases by the native plasmids of P. syringae impacts their transfer frequency and, by providing functional diversity and redundancy, adds robustness to the conjugation system.

Molecular and functional characterization of chicken interleukin 1 receptor 2 (chIL-1R2)

Interleukin-1 receptor type 2 (IL1R2) is a decoy receptor for exogenous IL-1. However, its functional role in chicken immunity is poorly understood. Herein, chicken IL-1R2 (chIL-1R2) was identified and functionally characterized in vivo and in vitro. The chIL-1R2 coding sequence includes 1,236 nucleotides encoding 412 amino acids, is highly conserved, and has a close relationship with its mammalian counterpart. Its extracellular region has three Ig-like domains but no TIR domain for intracellular signaling. Using ELISA, the recombinant chIL-1R2 protein was demonstrated to specifically bind to the chicken IL-1β. ChIL-1R2 mRNA expression was shown to be higher in the spleen, lung, kidney, small intestine, and liver. The expression of chIL-1R2 and chIL-1R1 was significantly upregulated in DF-1 cells treated with poly (I:C), but significantly downregulated in the presence of NF-κB, JNK, and MEK inhibitors, indicating that the NF-κB, JNK, and MEK signaling pathways are required for the transcriptional regulation of chIL-1R1 and chIL-1R2 expression. It is worth noting that while the p30 MAPK pathway was required for chIL-1R1 expression, it was not required for chIL-1R2 expression. Furthermore, chIL-1R2 expression increased as early as day 1, and then significantly decreased until day 3, while chIL-1R1 was dramatically upregulated in four organs of chickens infected with the highly pathogenic avian influenza virus (HPAIV). These findings indicate that chIL-1R1 and chIL-1R2 may play a crucial in innate and adaptive immune responses toward HPAIV infection. In summary the present study showed that chIL-1R2 binds to chIL-1β antibody. ChIL-1R2 expression can be induced by a viral infection, and may be regulated through NF-κB/JNK/MEK-mediated signaling pathways.

Proteomic analysis of Sarcoptes scabiei reveals that proteins differentially expressed between eggs and female adult stages are involved predominantly in genetic information processing, metabolism and/or host-parasite interactions

Presently, there is a dearth of proteomic data for parasitic mites and their relationship with the host animals. Here, using a high throughput LC-MS/MS-based approach, we undertook the first comprehensive, large-scale proteomic investigation of egg and adult female stages of the scabies mite, Sarcoptes scabiei-one of the most important parasitic mites of humans and other animals worldwide. In total, 1,761 S. scabiei proteins were identified and quantified with high confidence. Bioinformatic analyses revealed differentially expressed proteins to be involved predominantly in biological pathways or processes including genetic information processing, energy (oxidative phosphorylation), nucleotide, amino acid, carbohydrate and/or lipid metabolism, and some adaptive processes. Selected, constitutively and highly expressed proteins, such as peptidases, scabies mite inactivated protease paralogues (SMIPPs) and muscle proteins (myosin and troponin), are proposed to be involved in key biological processes within S. scabiei, host-parasite interactions and/or the pathogenesis of scabies. These proteomic data will enable future molecular, biochemical and physiological investigations of early developmental stages of S. scabiei and the discovery of novel interventions, targeting the egg stage, given its non-susceptibility to acaricides currently approved for the treatment of scabies in humans.

mTOR Contributes to the Proteome Diversity through Transcriptome-Wide Alternative Splicing

The mammalian target of rapamycin (mTOR) pathway is crucial in energy metabolism and cell proliferation. Previously, we reported transcriptome-wide 3′-untranslated region (UTR) shortening by alternative polyadenylation upon mTOR activation and its impact on the proteome. Here, we further interrogated the mTOR-activated transcriptome and found that hyperactivation of mTOR promotes transcriptome-wide exon skipping/exclusion, producing short isoform transcripts from genes. This widespread exon skipping confers multifarious regulations in the mTOR-controlled functional proteomics: AS in coding regions widely affects the protein length and functional domains. They also alter the half-life of proteins and affect the regulatory post-translational modifications. Among the RNA processing factors differentially regulated by mTOR signaling, we found that SRSF3 mechanistically facilitates exon skipping in the mTOR-activated transcriptome. This study reveals a role of mTOR in AS regulation and demonstrates that widespread AS is a multifaceted modulator of the mTOR-regulated functional proteome.

Isolation and characterization of a main porin from the outer membrane of Salinibacter ruber

Salinibacter ruber is an extremophilic bacterium able to grow in high-salts environments, such as saltern crystallizer ponds. This halophilic bacterium is red-pigmented due to the production of several carotenoids and their derivatives. Two of these pigment molecules, salinixanthin and retinal, are reported to be essential cofactors of the xanthorhodopsin, a light-driven proton pump unique to this bacterium. Here, we isolate and characterize an outer membrane porin-like protein that retains salinixanthin. The characterization by mass spectrometry identified an unknown protein whose structure, predicted by AlphaFold, consists of a 8 strands beta-barrel transmembrane organization typical of porins. The protein is found to be part of a functional network clearly involved in the outer membrane trafficking. Cryo-EM micrographs showed the shape and dimensions of a particle comparable with the ones of the predicted structure. Functional implications, with respect to the high representativity of this protein in the outer membrane fraction, are discussed considering its possible role in primary functions such as the nutrients uptake and the homeostatic balance. Finally, also a possible involvement in balancing the charge perturbation associated with the xanthorhodopsin and ATP synthase activities is considered.