The Pfam protein families database

Identification of type VI secretion system effector-immunity pairs using structural bioinformatics

The type VI secretion system (T6SS) is an important mediator of microbe-microbe and microbe-host interactions. Gram-negative bacteria use the T6SS to inject T6SS effectors (T6Es), which are usually proteins with toxic activity, into neighboring cells. Antibacterial effectors have cognate immunity proteins that neutralize self-intoxication. Here, we applied novel structural bioinformatic tools to perform systematic discovery and functional annotation of T6Es and their cognate immunity proteins from a dataset of 17,920 T6SS-encoding bacterial genomes. Using structural clustering, we identified 517 putative T6E families, outperforming sequence-based clustering. We developed a logistic regression model to reliably quantify protein-protein interaction of new T6E-immunity pairs, yielding candidate immunity proteins for 231 out of the 517 T6E families. We used sensitive structure-based annotation which yielded functional annotations for 51% of the T6E families, again outperforming sequence-based annotation. Next, we validated four novel T6E-immunity pairs using basic experiments in E. coli. In particular, we showed that the Pfam domain DUF3289 is a homolog of Colicin M and that DUF943 acts as its cognate immunity protein. Furthermore, we discovered a novel T6E that is a structural homolog of SleB, a lytic transglycosylase, and identified a specific glutamate that acts as its putative catalytic residue. Overall, this study applies novel structural bioinformatic tools to T6E-immunity pair discovery, and provides an extensive database of annotated T6E-immunity pairs.

Unlocking the power of AI models: exploring protein folding prediction through comparative analysis

Protein structure determination has made progress with the aid of deep learning models, enabling the prediction of protein folding from protein sequences. However, obtaining accurate predictions becomes essential in certain cases where the protein structure remains undescribed. This is particularly challenging when dealing with rare, diverse structures and complex sample preparation. Different metrics assess prediction reliability and offer insights into result strength, providing a comprehensive understanding of protein structure by combining different models. In a previous study, two proteins named ARM58 and ARM56 were investigated. These proteins contain four domains of unknown function and are present in Leishmania spp. ARM refers to an antimony resistance marker. The study's main objective is to assess the accuracy of the model's predictions, thereby providing insights into the complexities and supporting metrics underlying these findings. The analysis also extends to the comparison of predictions obtained from other species and organisms. Notably, one of these proteins shares an ortholog with Trypanosoma cruzi and Trypanosoma brucei, leading further significance to our analysis. This attempt underscored the importance of evaluating the diverse outputs from deep learning models, facilitating comparisons across different organisms and proteins. This becomes particularly pertinent in cases where no previous structural information is available.

On knowing a gene: A distributional hypothesis of gene function

As words can have multiple meanings that depend on sentence context, genes can have various functions that depend on the surrounding biological system. This pleiotropic nature of gene function is limited by ontologies, which annotate gene functions without considering biological contexts. We contend that the gene function problem in genetics may be informed by recent technological leaps in natural language processing, in which representations of word semantics can be automatically learned from diverse language contexts. In contrast to efforts to model semantics as "is-a" relationships in the 1990s, modern distributional semantics represents words as vectors in a learned semantic space and fuels current advances in transformer-based models such as large language models and generative pre-trained transformers. A similar shift in thinking of gene functions as distributions over cellular contexts may enable a similar breakthrough in data-driven learning from large biological datasets to inform gene function.

Polysaccharide utilization loci encoded DUF1735 likely functions as membrane-bound spacer for carbohydrate active enzymes

Proteins featuring the Domain of Unknown Function 1735 are frequently found in Polysaccharide Utilization Loci, yet their role remains unknown. The domain and vicinity analyzer programs we developed mine the Kyoto Encyclopedia of Genes and Genomes and UniProt to enhance the functional prediction of DUF1735. Our datasets confirmed the exclusive presence of DUF1735 in Bacteroidota genomes, with Bacteroidetes thetaiotaomicron harboring 46 copies. Notably, 97.8% of DUF1735 are encoded in PULs, and 89% are N-termini of multimodular proteins featuring C-termini like Laminin_G_3, F5/8-typeC, and GH18 domains. Predominantly possessing a predicted lipoprotein signal peptide and sharing an immunoglobulin-like β-sandwich fold with the BACON domain and the N-termini of SusE/F, DUF1735 likely functions as N-terminal, membrane-bound spacer for diverse C-termini involved in PUL-mediated carbohydrate utilization.

Chromosome-level genome assembly of the diploid oat species Avena longiglumis

Diploid wild oat Avena longiglumis has nutritional and adaptive traits which are valuable for common oat (A. sativa) breeding. The combination of Illumina, Nanopore and Hi-C data allowed us to assemble a high-quality chromosome-level genome of A. longiglumis (ALO), evidenced by contig N50 of 12.68 Mb with 99% BUSCO completeness for the assembly size of 3,960.97 Mb. A total of 40,845 protein-coding genes were annotated. The assembled genome was composed of 87.04% repetitive DNA sequences. Dotplots of the genome assembly (PI657387) with two published ALO genomes were compared to indicate the conservation of gene order and equal expansion of all syntenic blocks among three genome assemblies. Two recent whole-genome duplication events were characterized in genomes of diploid Avena species. These findings provide new knowledge for the genomic features of A. longiglumis, give information about the species diversity, and will accelerate the functional genomics and breeding studies in oat and related cereal crops.

Genome-wide analysis of the mulberry (Morus abla L.) GH9 gene family and the functional characterization of MaGH9B6 during the development of the abscission zone

Plant glycoside hydrolase family 9 genes (GH9s) are widely distributed in plants and involved in a variety of cellular and physiological processes. In the current study, nine GH9 genes were identified in the mulberry and were divided into two subfamilies based on the phylogenetic analysis. Conserved motifs and gene structure analysis suggested that the evolution of the two subfamilies is relatively conserved and the glycoside hydrolase domain almost occupy the entire coding region of the GH9s gene. Only segmental duplication has played a role in the expansion of gene family. Collinearity analysis showed that mulberry GH9s had the closest relationship with poplar GH9s. MaGH9B1, MaGH9B6, MaGH9B5, and MaGH9B3 were detected to have transcript accumulation in the stalk of easy-to drop mature fruit drop, suggesting that these could play a role in mulberry fruit drop. Multiple cis-acting elements related to plant hormones and abiotic stress responses were found in the mulberry GH9 promoter regions and showed different activities under exogenous abscisic acid (ABA) and 2,4- dichlorophenoxyacetic acid (2,4-D) stresses. We found that the lignin content in the fruit stalk decreased with the formation of the abscission zone (AZ), which could indirectly reflect the formation process of the AZ. These results provide a theoretical basis for further research on the role of GH9s in mulberry abscission.

Transcriptome analysis of polyploid giant cancer cells and their progeny reveals a functional role for p21 in polyploidization and depolyploidization

Polyploid giant cancer cells (PGCC) are frequently detected in tumors and are increasingly recognized for their roles in chromosomal instability and associated genome evolution that leads to cancer recurrence. We previously reported that therapy stress promotes polyploidy, and that acid ceramidase plays a role in depolyploidization. In this study, we used an RNA-seq approach to gain a better understanding of the underlying transcriptomic changes that occur as cancer cells progress through polyploidization and depolyploidization. Our results revealed gene signatures that are associated with disease-free and/or overall survival in several cancers and identified the cell cycle inhibitor CDKN1A/p21 as the major hub in PGCC and early progeny. Increased expression of p21 in PGCC was limited to the cytoplasm. We previously demonstrated that the sphingolipid enzyme acid ceramidase is dispensable for polyploidization upon therapy stress but plays a crucial role in depolyploidization. The current study demonstrates that treatment of cells with ceramide is not sufficient for p53-independent induction of p21 and that knockdown of acid ceramidase, which hydrolyzes ceramide, does not interfere with upregulation of p21. In contrast, blocking the expression of p21 with UC2288 prevented the induction of acid ceramidase and inhibited both the formation of PGCC from parental cells as well as the generation of progeny from PGCC. Taken together, our data suggest that p21 functions upstream of acid ceramidase and plays an important role in polyploidization and depolyploidization.

Effect of a probiotic and an antibiotic on the mobilome of the porcine microbiota

Introduction: To consider the growing health issues caused by antibiotic resistance from a "one health" perspective, the contribution of meat production needs to be addressed. While antibiotic resistance is naturally present in microbial communities, the treatment of farm animals with antibiotics causes an increase in antibiotic resistance genes (ARG) in the gut microbiome. Pigs are among the most prevalent animals in agriculture; therefore, reducing the prevalence of antibiotic-resistant bacteria in the pig gut microbiome could reduce the spread of antibiotic resistance. Probiotics are often studied as a way to modulate the microbiome and are, therefore, an interesting way to potentially decrease antibiotic resistance. Methods: To assess the efficacy of a probiotic to reduce the prevalence of ARGs in the pig microbiome, six pigs received either treatment with antibiotics (tylvalosin), probiotics (Pediococcus acidilactici MA18/5M; Biopower® PA), or a combination of both. Their faeces and ileal digesta were collected and DNA was extracted for whole genome shotgun sequencing. The reads were compared with taxonomy and ARG databases to identify the taxa and resistance genes in the samples. Results: The results showed that the ARG profiles in the faeces of the antibiotic and combination treatments were similar, and both were different from the profiles of the probiotic treatment (p < 0.05). The effects of the treatments were different in the digesta and faeces. Many macrolide resistance genes were detected in a higher proportion in the microbiome of the pigs treated with antibiotics or the combination of probiotics and antibiotics. Resistance-carrying conjugative plasmids and horizontal transfer genes were also amplified in faeces samples for the antibiotic and combined treatments. There was no effect of treatment on the short chain fatty acid content in the digesta or the faeces. Conclusion: There is no positive effect of adding probiotics to an antibiotic treatment when these treatments are administered simultaneously.

Chromosome-level genome assembly of Platycarya strobilacea

Platycarya strobilacea belongs to the walnut family (Juglandaceae), is commonly known as species endemic to East Asia, and is an ecologically important, wind pollinated, woody deciduous tree. To facilitate this ancient tree for the ecological value and conservation of this ancient tree, we report a new high-quality genome assembly of P. strobilacea. The genome size was 677.30 Mb, with a scaffold N50 size of 45,791,698 bp, and 98.43% of the assembly was anchored to 15 chromosomes. We annotated 32,246 protein-coding genes in the genome, of which 96.30% were functionally annotated in six databases. This new high-quality assembly of P. strobilacea provide valuable resource for the phylogenetic and evolutionary analysis of the walnut family and angiosperm.

Engineering serine hydroxymethyltransferases for efficient synthesis of L-serine in Escherichia coli

L-serine is a high-value amino acid widely used in the food, medicine, and cosmetic industries. However, the low yield of L-serine has limited its industrial production. In this study, a cellular factory for efficient synthesis of L-serine was obtained by engineering the serine hydroxymethyltransferases (SHMT). Firstly, after screening the SHMT from Alcanivorax dieselolei by genome mining, a mutant AdSHMTE266M with high thermal stability was identified through rational design. Subsequently, an iterative saturating mutant library was constructed by using coevolutionary analysis, and a mutant AdSHMTE160L/E193Q with enzyme activity 1.35 times higher than AdSHMT was identified. Additionally, the target protein AdSHMTE160L/E193Q/E266M was efficiently overexpressed by improving its mRNA stability. Finally, combining the substrate addition strategy and system optimization, the optimized strain BL21/pET28a-AdSHMTE160L/E193Q/E266M-5'UTR-REP3S16 produced 106.06 g/L L-serine, which is the highest production to date. This study provides new ideas and insights for the engineering design of SHMT and the industrial production of L-serine.

Verticillium dahliae VdPBP1 Transcription Factor Is Required for Hyphal Growth, Virulence, and Microsclerotia Formation

Verticillium dahliae, a fungal pathogen that affects more than 200 plant species, including tomatoes, requires specific proteins for its early steps in plant infection. One such crucial protein, VdPBP1, exhibits high expression in the presence of tomato roots. Its 313-amino acid C-terminal section restores adhesion in nonadhesive Saccharomyces cerevisiae strains. To uncover its role, we employed a combination of bioinformatics, genetics, and morphological analyses. Our findings underscore the importance of VdPBP1 in fungal growth and pathogenesis. Bioinformatic analysis revealed that the VdPBP1 gene consists of four exons and three introns, encoding a 952-codon reading frame. The protein features a 9aaTAD domain, LsmAD, and PAB1 DNA-binding sites, as well as potential nuclear localization and transmembrane helix signals. Notably, the deletion of a 1.1 kb fragment at the gene's third end impedes microsclerotia formation and reduces pathogenicity. Mutants exhibit reduced growth and slower aerial mycelial development compared to the wild type. The VdPBP1 deletion strain does not induce disease symptoms in tomato plants. Furthermore, VdPBP1 deletion correlates with downregulated microsclerotia formation-related genes, and promoter analysis reveals regulatory elements, including sites for Rfx1, Mig1, and Ste12 proteins. Understanding the regulation and target genes of VdPBP1 holds promise for managing Verticillium wilt disease and related fungal pathogens.

Genome-Wide Identification of Trehalose-6-phosphate Synthase (TPS) Gene Family Reveals the Potential Role in Carbohydrate Metabolism in Peach

Trehalose-6-phosphate synthase (TPS) is essential for plant growth and development, linking trehalose-6-phosphate (T6P) to carbon metabolism. However, little is known about the TPS gene family in peaches and their potential roles in regulating carbohydrates in peach fruit. In this study, nine TPS genes were identified in the peach genome and named according to the homologous genes in Arabidopsis. Phylogenetic analysis showed that three subfamilies were identified, including TPSI, TPSII-1, and TPSII-2, which were also consistent with gene structure analysis. Considerable cis-elements were enriched in the promoters, including plant hormone-related elements. Tissue-specific analysis showed that these TPS genes were mainly expressed in leaves, stems, and fruit, showing different expression patterns for each gene. In addition, during fruit development, the content of trehalose-6-phosphate (T6P) was positively correlated with the expression of PpTPS7a and negatively with sucrose non-fermenting-1-related kinase 1 (SnRK1) activity. Transient overexpression and silencing of PpTPS7a in peach fruit validated its function in regulating T6P content and SnRK1 activity.

Polyploid Genome Assembly Provides Insights into Morphological Development and Ascorbic Acid Accumulation of Sauropus androgynus

Sauropus androgynus (S. androgynus) (2n = 4x = 52) is one of the most popular functional leafy vegetables in South and Southeast Asia. With its rich nutritional and pharmaceutical values, it has traditionally had widespread use for dietary and herbal purposes. Here, the genome of S. androgynus was sequenced and assembled, revealing a genome size of 1.55 Gb with 26 pseudo-chromosomes. Phylogenetic analysis traced back the divergence of Sauropus from Phyllanthus to approximately 29.67 million years ago (Mya). Genome analysis revealed that S. androgynus polyploidized around 20.51 Mya and shared a γ event about 132.95 Mya. Gene function analysis suggested that the expansion of pathways related to phloem development, lignin biosynthesis, and photosynthesis tended to result in the morphological differences among species within the Phyllanthaceae family, characterized by varying ploidy levels. The high accumulation of ascorbic acid in S. androgynus was attributed to the high expression of genes associated with the L-galactose pathway and recycling pathway. Moreover, the expanded gene families of S. androgynus exhibited multiple biochemical pathways associated with its comprehensive pharmacological activity, geographic adaptation and distinctive pleasurable flavor. Altogether, our findings represent a crucial genomic asset for S. androgynus, casting light on the intricate ploidy within the Phyllanthaceae family.

Unraveling the dataset transcriptomic response of Hydrangea macrophylla stem to mechanical stimulation: De novo assembly and functional annotation

A crucial attribute of potted ornamental plants is compactness characterized by well branched plants with rather short stems bearing numerous flowers. To gain plant compactness, producers use plant growth regulators (PGRs), in particular growth retardants during culture. However, due to their negative environmental impacts, growth retardants are progressively withdrawn from the market. As a response, eco-friendly alternative methods to chemicals need to be developed. One method consists in mimicking mechanical stimulation (MS) imposed by wind on plants which causes reduction in stem elongation, an increase in stem diameter and an increase in branching, all contributing to plant compactness. So far, few plant species were studied under MS and little is known on molecular response mechanisms to MS. This first transcriptomic data after MS in Hydrangea macrophylla will contribute unravelling how plants respond to mechanical stimuli. RNAseq data were obtained from total mRNA of stems collected 15 min before MS and 1, 3, 24 and 72 h after MS treatment. RNA from non-MS treated plants were used as control. MS treatment consisted in 12 consecutive bendings (i.e. 6 forth and 6 back) applied at 9 a.m. during 1 h and for a single day. From RNAseq data a de novo assembly of the transcriptome was produced and 78,398 transcripts functionally annotated. These transcriptomic data also contribute to a better knowledge of how outdoor crop respond to the increasing frequency of strong harmful winds under climate change.

Learning sequence, structure, and function representations of proteins with language models

The sequence-structure-function relationships that ultimately generate the diversity of extant observed proteins is complex, as proteins bridge the gap between multiple informational and physical scales involved in nearly all cellular processes. One limitation of existing protein annotation databases such as UniProt is that less than 1% of proteins have experimentally verified functions, and computational methods are needed to fill in the missing information. Here, we demonstrate that a multi-aspect framework based on protein language models can learn sequence-structure-function representations of amino acid sequences, and can provide the foundation for sensitive sequence-structure-function aware protein sequence search and annotation. Based on this model, we introduce a multi-aspect information retrieval system for proteins, Protein-Vec, covering sequence, structure, and function aspects, that enables computational protein annotation and function prediction at tree-of-life scales.

VDR promotes testosterone synthesis in mouse Leydig cells via regulation of cholesterol side chain cleavage cytochrome P450 (Cyp11a1) expression

BACKGROUND

The vitamin D receptor (VDR) mediates the pleiotropic biological actions that include osteoporosis, immune responses and androgen synthesis.VDR is widely expressed in testis cells such as Leydig cells, Sertoli cells, and sperm. The levels of steroids are critical for sexual development. In the early stage of steroidogenesis, cholesterol is converted to pregnenolone (precursor of most steroid hormones) by cholesterol side-chain lyase (CYP11A1), which eventually synthesizes the male hormone testosterone.

OBJECTIVE

This study aims to reveal how VDR regulates CYP11A1 expression and affects testosterone synthesis in murine Leydig cells.

METHODS

The levels of VDR, CYP11A1 were determined by quantitative real-time polymerase chain reaction (RT-qPCR) or western blot. Targeted relationship between VDR and Cyp11a1 was evaluated by dual-luciferase reporter assay. The levels of testosterone concentrations in cell culture media serum by enzyme-linked immunosorbent assay (ELISA).

RESULTS

Phylogenetic and motif analysis showed that the Cyp11a1 family had sequence loss, which may have special biological functions during evolution. The results of promoter prediction showed that vitamin D response element (VDRE) existed in the upstream promoter region of murine Cyp11a1. Dual-luciferase assay confirmed that VDR could bind candidate VDREs in upstream region of Cyp11a1, and enhance gene expression. Tissue distribution and localizatio analysis showed that Cyp11a1 was mainly expressed in testis, and dominantly existed in murine Leydig cells. Furthermore, over-expression VDR and CYP11A1 significantly increased testosterone synthesis in mice Leydig cells.

CONCLUSIONS

Active vitamin D3 (VD3) and Vdr interference treatment showed that VD3/VDR had a positive regulatory effect on Cyp11a1 expression and testosterone secretion. VDR promotes testosterone synthesis in male mice by up-regulating Cyp11a1 expression, which played an important role for male reproduction.

Systematic comparison of genome information processing and boundary recognition tools used for genomic island detection

Genomic islands are fragments of foreign DNA that are found in bacterial and archaeal genomes, and are typically associated with symbiosis or pathogenesis. While numerous genomic island detection methods have been proposed, there has been limited evaluation of the efficiency of the genome information processing and boundary recognition tools. In this study, we conducted a review of the statistical methods involved in genomic signatures, host signature extraction, informative signature selection, divergence measures, and boundary detection steps in genomic island prediction. We compared the performances of these methods on simulated experiments using alien fragments obtained from both artificial and real genomes. Our results indicate that among the nine genomic signatures evaluated, genomic signature frequency and full probability performed the best. However, their performance declined when normalized to their expectations and variances, such as Z-score and composition vector. Based on our experiments of the E. coli genome, we found that the confidence intervals of the window variances achieved the best performance in the signature extraction of the host, with the best confidence interval being 1.5-2 times the standard error. Ordered kurtosis was most effective in selecting informative signatures from a single genome, without requiring prior knowledge from other datasets. Among the three divergence measures evaluated, the two-sample t-test was the most successful, and a non-overlapping window with a small eye window (size 2) was best suited for identifying compositionally distinct regions. Finally, the maximum of the Markovian Jensen-Shannon divergence score, in terms of GC-content bias, was found to make boundary detection faster while maintaining a similar error rate.

Comparative analysis of KNOX genes and their expression patterns under various treatments in Dendrobium huoshanense

Introduction

KNOX plays a pivotal role in governing plant growth, development, and responses to diverse abiotic and biotic stresses. However, information on the relationship between the KNOX gene family and expression levels under different treatments in Dendrobium is still limited.

Methods

To address this problem, we first used bioinformatics methods and revealed the presence of 19 KNOX genes distributed among 13 chromosomes in the Dendrobium huoshanense genome. Through an analysis of phylogenetic relationships, these genes were classified into three distinct clades: class I, class II, and class M. Our investigation included promoter analysis, revealing various cis-acting elements associated with hormones, growth and development, and abiotic stress responses. Additionally, qRT-PCR experiments were conducted to assess the expression patterns of DhKNOX genes under different treatments, including ABA, MeJA, SA, and drought.

Results

The results demonstrated differential expression of DhKNOX genes in response to these treatments, thereby highlighting their potential roles in stress adaptation.

Discussion

Overall, our results contribute important insights for further investigations into the functional characterization of the Dendrobium KNOX gene family, shedding light on their roles in plant development and stress responses.

Pangenomic Analysis of Nucleo-Cytoplasmic Large DNA Viruses. I: The Phylogenetic Distribution of Conserved Oxygen-Dependent Enzymes Reveals a Capture-Gene Process

The Nucleo-Cytoplasmic Large DNA Viruses (NCLDVs) infect a wide range of eukaryotic species, including amoeba, algae, fish, amphibia, arthropods, birds, and mammals. This group of viruses has linear or circular double-stranded DNA genomes whose size spans approximately one order of magnitude, from 100 to 2500 kbp. The ultimate origin of this peculiar group of viruses remains an open issue. Some have argued that NCLDVs' origin may lie in a bacteriophage ancestor that increased its genome size by subsequent recruitment of eukaryotic and bacterial genes. Others have suggested that NCLDVs families originated from cells that underwent an irreversible process of genome reduction. However, the hypothesis that a number of NCLDVs sequences have been recruited from the host genomes has been largely ignored. In the present work, we have performed pangenomic analyses of each of the seven known NCLDVs families. We show that these families' core- and shell genes have cellular homologs, supporting possible escaping-gene events as part of its evolution. Furthermore, the detection of sequences that belong to two protein families (small chain ribonucleotide reductase and Erv1/Air) and to one superfamily [2OG-Fe(II) oxygenases] that are for distribution in all NCLDVs core and shell clusters encoding for oxygen-dependent enzymes suggests that the highly conserved core these viruses originated after the Proterozoic Great Oxidation Event that transformed the terrestrial atmosphere 2.4-2.3 Ga ago.

Characterization of a GH10 extremely thermophilic xylanase from the metagenome of hot spring for prebiotic production

A xylanase gene (named xyngmqa) was identified from the metagenomic data of the Gumingquan hot spring (92.5 °C, pH 9.2) in Tengchong City, Yunnan Province, southwest China. It showed the highest amino acid sequence identity (82.70%) to endo-1,4-beta-xylanase from Thermotoga caldifontis. A constitutive expression plasmid (denominated pSHY211) and double-layer plate (DLP) method were constructed for cloning, expression, and identification of the XynGMQA gene. The XynGMQA gene was synthesized and successfully expressed in Escherichia coli DH5α. XynGMQA exhibited optimal activity at 90 °C and pH 4.6, being thermostable by maintaining 100% of its activity after 2 h incubated at 80 °C. Interestingly, its enzyme activity was enhanced by high temperatures (70 and 80 °C) and low pH (3.0-6.0). About 150% enzyme activity was detected after incubation at 70 °C for 20 to 60 min or 80 °C for 10 to 40 min, and more than 140% enzyme activity after incubation at pH 3.0 to 6.0 for 12 h. Hydrolytic products of beechwood xylan with XynGMQA were xylooligosaccharides, including xylobiose (X2), xylotriose (X3), and xylotetraose (X4). These properties suggest that XynGMQA as an extremely thermophilic xylanase, may be exploited for biofuel and prebiotic production from lignocellulosic biomass.

Identification of Dioscorea opposite Thunb. CDPK gene family reveals that DoCDPK20 is related to heat resistance

Temperature affects the growth and yield of yam (Dioscorea opposite Thunb.), and calcium-dependent protein kinases (CDPKs) play an important role in the plant stress response. However, there has been a lack of system analyses of yam's CDPK gene family. In this study, 29 CDPK transcriptome sequences with complete open reading frames (ORFs) were identified from yam RNA sequencing data. The sequences were classified into four groups (I-VI) using phylogenetic analysis. Two DoCDPK genes were randomly selected from each group and the gene patterns of yam leaves were determined using quantitative real-time PCR (qRT-PCR) under high and low temperature stress in order to show their unique functions in mediating specific responses. Among them, DoCDPK20 was significantly induced in high temperatures. The pPZP221-DoCDPK20 was transformed into tobacco leaves using an agrobacterium-mediated method. Under high temperature stress, DoCDPK20 overexpression reduced photosynthesis and improved heat tolerance in transgenic tobacco. Our research offers meaningful perspectives into CDPK genes and new avenues for the genetic engineering and molecular breeding of yam.

Protein remote homology detection and structural alignment using deep learning

Exploiting sequence-structure-function relationships in biotechnology requires improved methods for aligning proteins that have low sequence similarity to previously annotated proteins. We develop two deep learning methods to address this gap, TM-Vec and DeepBLAST. TM-Vec allows searching for structure-structure similarities in large sequence databases. It is trained to accurately predict TM-scores as a metric of structural similarity directly from sequence pairs without the need for intermediate computation or solution of structures. Once structurally similar proteins have been identified, DeepBLAST can structurally align proteins using only sequence information by identifying structurally homologous regions between proteins. It outperforms traditional sequence alignment methods and performs similarly to structure-based alignment methods. We show the merits of TM-Vec and DeepBLAST on a variety of datasets, including better identification of remotely homologous proteins compared with state-of-the-art sequence alignment and structure prediction methods.

Transcriptomic Analysis of Light-Induced Genes in Nasonia vitripennis: Possible Implications for Circadian Light Entrainment Pathways

Circadian entrainment to the environmental day-night cycle is essential for the optimal use of environmental resources. In insects, opsin-based photoreception in the compound eye and ocelli and CRYPTOCHROME1 (CRY1) in circadian clock neurons are thought to be involved in sensing photic information, but the genetic regulation of circadian light entrainment in species without light-sensitive CRY1 remains unclear. To elucidate a possible CRY1-independent light transduction cascade, we analyzed light-induced gene expression through RNA-sequencing in Nasonia vitripennis. Entrained wasps were subjected to a light pulse in the subjective night to reset the circadian clock, and light-induced changes in gene expression were characterized at four different time points in wasp heads. We used co-expression, functional annotation, and transcription factor binding motif analyses to gain insight into the molecular pathways in response to acute light stimulus and to form hypotheses about the circadian light-resetting pathway. Maximal gene induction was found after 2 h of light stimulation (1432 genes), and this included the opsin gene opblue and the core clock genes cry2 and npas2. Pathway and cluster analyses revealed light activation of glutamatergic and GABA-ergic neurotransmission, including CREB and AP-1 transcription pathway signaling. This suggests that circadian photic entrainment in Nasonia may require pathways that are similar to those in mammals. We propose a model for hymenopteran circadian light-resetting that involves opsin-based photoreception, glutamatergic neurotransmission, and gene induction of cry2 and npas2 to reset the circadian clock.

Antibiotic Resistance Mediated by Escherichia coli in Kuwait Marine Environment as Revealed through Genomic Analysis

Antibiotic-resistance gene elements (ARGEs) such as antibiotic-resistance genes (ARGs), integrons, and plasmids are key to the spread of antimicrobial resistance (AMR) in marine environments. Kuwait's marine area is vulnerable to sewage contaminants introduced by numerous storm outlets and indiscriminate waste disposal near recreational beaches. Therefore, it has become a significant public health issue and warrants immediate investigation. Coliforms, especially Gram-negative Escherichia coli, have been regarded as significant indicators of recent fecal pollution and carriers of ARGEs. In this study, we applied a genome-based approach to identify ARGs' prevalence in E. coli isolated from mollusks and coastal water samples collected in a previous study. In addition, we investigated the plasmids and intl1 (class 1 integron) genes coupled with the ARGs, mediating their spread within the Kuwait marine area. Whole-genome sequencing (WGS) identified genes resistant to the drug classes of beta-lactams (blaCMY-150, blaCMY-42, blaCTX-M-15, blaDHA-1, blaMIR-1, blaOKP-B-15, blaOXA-1, blaOXA-48, blaTEM-1B, blaTEM-35), trimethoprim (dfrA14, dfrA15, dfrA16, dfrA1, dfrA5, dfrA7), fluroquinolone (oqxA, oqxB, qnrB38, qnrB4, qnrS1), aminoglycoside (aadA2, ant(3'')-Ia, aph(3'')-Ib, aph(3')-Ia, aph(6)-Id), fosfomycin (fosA7, fosA_6, fosA, fosB1), sulfonamide (sul1, sul2, sul3), tetracycline (tet-A, tet-B), and macrolide (mph-A). The MFS-type drug efflux gene mdf-A is also quite common in E. coli isolates (80%). The plasmid ColRNAI was also found to be prevalent in E. coli. The integron gene intI1 and gene cassettes (GC) were reported to be in 36% and 33%, respectively, of total E. coli isolates. A positive and significant (p < 0.001) correlation was observed between phenotypic AMR-intl1 (r = 0.311) and phenotypic AMR-GC (r = 0.188). These findings are useful for the surveillance of horizontal gene transfer of AMR in the marine environments of Kuwait.

Crystal structure of domain of unknown function 507 (DUF507) reveals a new protein fold

The crystal structure of the domain of unknown function family 507 protein from Aquifex aeolicus is reported (AaDUF507, UniProt O67633, 183 residues). The structure was determined in two space groups (C2221 and P3221) at 1.9 Å resolution. The phase problem was solved by molecular replacement using an AlphaFold model as the search model. AaDUF507 is a Y-shaped α-helical protein consisting of an anti-parallel 4-helix bundle base and two helical arms that extend 30-Å from the base. The two crystal structures differ by a 25° rigid body rotation of the C-terminal arm. The tertiary structure exhibits pseudo-twofold symmetry. The structural symmetry mirrors internal sequence similarity: residues 11-57 and 102-148 are 30% identical and 53% similar with an E-value of 0.002. In one of the structures, electron density for an unknown ligand, consistent with nicotinamide or similar molecule, may indicate a functional site. Docking calculations suggest potential ligand binding hot spots in the region between the helical arms. Structure-based query of the Protein Data Bank revealed no other protein with a similar tertiary structure, leading us to propose that AaDUF507 represents a new protein fold.

SW16-7, a Novel Ackermannviridae Bacteriophage with Highly Effective Lytic Activity Targets Salmonella enterica Serovar Weltevreden

Salmonella enterica serovar Weltevreden is a foodborne pathogen commonly transmitted through fresh vegetables and seafood. In this study, a lytic phage, SW16-7, was isolated from medical sewage, demonstrating high infectivity against S. Weltevreden, S. London, S. Meleagridis, and S. Give of Group O:3. In vitro inhibition assays revealed its effective antibacterial effect for up to 12 h. Moreover, analysis using the Comprehensive Antibiotic Resistance Database (CARD) and the Virulence Factor Database (VFDB) showed that SW16-7's genome does not contain any virulence factors or antibiotic resistance genes, indicating its potential as a promising biocontrol agent against S. Weltevreden. Additionally, a TSP gene cluster was identified in SW16-7's genome, with TSP1 and TSP2 showing a high similarity to lysogenic phages ε15 and ε34, respectively, in the C-terminal region. The whole-genome phylogenetic analysis classified SW16-7 within the Ackermannviridae family and indicated a close relationship with Agtrevirus, which is consistent with the ANI results.

The Pitaya Flower Tissue's Gene Differential Expression Analysis between Self-Incompatible and Self-Compatible Varieties for the Identification of Genes Involved in Self-Incompatibility Regulation

Self-incompatible pitaya varieties have low fruit-setting rates under natural conditions, leading to higher production costs and hindering industrial prosperity. Through transcriptome sequencing, we obtained the 36,900 longest transcripts (including 9167 new transcripts) from 60 samples of flowers. Samples were collected pre- and post-pollination (at 0 h, 0.5 h, 2 h, 4 h, and 12 h) from two varieties of pitaya (self-compatible Jindu No. 1 and self-incompatible Cu Sha). Using the RNA-Seq data and comparison of reference genomes, we annotated 28,817 genes in various databases, and 1740 genes were optimized in their structure for annotation. There were significant differences in the expression of differentially expressed genes (DEGs) in the pitaya stigmas under different pollination types, especially at the late post-pollination stage, where the expression of protease genes increasedal significantly under cross-pollination. We identified DEGs involved in the ribosomal, ubiquitination-mediated, and phyto-signaling pathways that may be involved in pitaya SI regulation. Based on the available transcriptome data and bioinformatics analysis, we tentatively identified HuS-RNase2 as a candidate gynogenetic S gene in the pitaya GSI system.

The Transcriptomic Analysis of the Response of Pinus massoniana to Drought Stress and a Functional Study on the ERF1 Transcription Factor

Pinus massoniana is a major fast-growing timber tree species planted in arid areas of south China, which has a certain drought-resistant ability. However, severe drought and long-term water shortage limit its normal growth and development. Therefore, in this study, physiological indices, and the transcriptome sequencing and cloning of AP2/ERF transcription factor of P. massonsiana were determined to clarify its molecular mechanism of drought stress. The results showed that stomatal conductance (Gs) content was significantly decreased, and superoxide dismutase (SOD) activity, and malondialdehyde (MDA) and abscisic acid (ABA) content were significantly increased under drought stress. Transcriptomic analysis revealed that compared to the control, 9, 3550, and 4142 unigenes with differential expression were identified by comparing plants subjected to light, moderate or severe drought. AP2/ERF with high expression was screened out for cloning. To investigate the biological functions of ERF1, it was over-expressed in wild-type Populus davdianaand × P. bolleana via the leaf disc method. Under drought stress, compared to wild-type plants, ERF1 over-expressing poplar lines (OE) maintained a higher photosynthetic rate and growth, while the transpiration rate and stomatal conductance significantly decreased and water use efficiency was improved, indicating that drought tolerance was enhanced. This study provides an insight into the molecular mechanism of drought stress adaptation in P. massoniana.

Genome-Wide Identification and Functional Analysis of the TIFY Family Genes in Response to Abiotic Stresses and Hormone Treatments in Tartary Buckwheat (Fagopyrum tataricum)

TIFY is a plant-specific gene family with four subfamilies: ZML, TIFY, PPD, and JAZ. Recently, this family was found to have regulatory functions in hormone stimulation, environmental response, and development. However, little is known about the roles of the TIFY family in Tartary buckwheat (Fagopyrum tataricum), a significant crop for both food and medicine. In this study, 18 TIFY family genes (FtTIFYs) in Tartary buckwheat were identified. The characteristics, motif compositions, and evolutionary relationships of the TIFY proteins, as well as the gene structures, cis-acting elements, and synteny of the TIFY genes, are discussed in detail. Moreover, we found that most FtTIFYs responded to various abiotic stresses (cold, heat, salt, or drought) and hormone treatments (ABA, MeJA, or SA). Through yeast two-hybrid assays, we revealed that two FtTIFYs, FtTIFY1 and FtJAZ7, interacted with FtABI5, a homolog protein of AtABI5 involved in ABA-mediated germination and stress responses, implying crosstalk between ABA and JA signaling in Tartary buckwheat. Furthermore, the overexpression of FtJAZ10 and FtJAZ12 enhanced the heat stress tolerance of tobacco. Consequently, our study suggests that the FtTIFY family plays important roles in responses to abiotic stress and provides two candidate genes (FtJAZ10 and FtJAZ12) for the cultivation of stress-resistant crops.

Transcriptome Screening of Long Noncoding RNAs and Their Target Protein-Coding Genes Unmasks a Dynamic Portrait of Seed Coat Coloration Associated with Anthocyanins in Tibetan Hulless Barley

Many plants have the capability to accumulate anthocyanins for coloration, and anthocyanins are advantageous to human health. In the case of hulless barley (Hordeum vulgare L. var. nudum), investigation into the mechanism of anthocyanin formation is limited to the level of protein-coding genes (PCGs). Here, we conducted a comprehensive bioinformatics analysis to identify a total of 9414 long noncoding RNAs (lncRNAs) in the seed coats of purple and white hulless barley along a developmental gradient. Transcriptome-wide profiles of lncRNAs documented several properties, including GC content fluctuation, uneven length, a diverse range of exon numbers, and a wide variety of transcript classifications. We found that certain lncRNAs in hulless barley possess detectable sequence conservation with Hordeum vulgare and other monocots. Furthermore, both differentially expressed lncRNAs (DElncRNAs) and PCGs (DEPCGs) were concentrated in the later seed development stages. On the one hand, DElncRNAs could potentially cis-regulate DEPCGs associated with multiple metabolic pathways, including flavonoid and anthocyanin biosynthesis in the late milk and soft dough stages. On the other hand, there was an opportunity for trans-regulated lncRNAs in the color-forming module to affect seed coat color by upregulating PCGs in the anthocyanin pathway. In addition, the interweaving of hulless barley lncRNAs and diverse TFs may function in seed coat coloration. Notably, we depicted a dynamic portrait of the anthocyanin synthesis pathway containing hulless barley lncRNAs. Therefore, this work provides valuable gene resources and more insights into the molecular mechanisms underlying anthocyanin accumulation in hulless barley from the perspective of lncRNAs, which facilitate the development of molecular design breeding in crops.

High-quality chromosome-level scaffolds of the plant bug Pachypeltis micranthus provide insights into the availability of Mikania micrantha control

BACKGROUND

The plant bug, Pachypeltis micranthus Mu et Liu (Hemiptera: Miridae), is an effective potential biological control agent for Mikania micrantha H.B.K. (Asteraceae; one of the most notorious invasive weeds worldwide). However, limited knowledge about this species hindered its practical application and research. Accordingly, sequencing the genome of this mirid bug holds great significance in controlling M. micrantha.

RESULTS

Here, 712.72 Mb high-quality chromosome-level scaffolds of P. micranthus were generated, of which 707.51 Mb (99.27%) of assembled sequences were anchored onto 15 chromosome-level scaffolds with contig N50 of 16.84 Mb. The P. micranthus genome had the highest GC content (42.43%) and the second highest proportion of repetitive sequences (375.82 Mb, 52.73%) than the three other mirid bugs (i.e., Apolygus lucorum, Cyrtorhinus lividipennis, and Nesidiocoris tenuis). Phylogenetic analysis showed that P. micranthus clustered with other mirid bugs and diverged from the common ancestor approximately 200 million years ago. Gene family expansion and/or contraction were analyzed, and significantly expanded gene families associated with P. micranthus feeding and adaptation to M. micrantha were manually identified. Compared with the whole body, transcriptome analysis of the salivary gland revealed that most of the upregulated genes were significantly associated with metabolism pathways and peptidase activity, particularly among cysteine peptidase, serine peptidase, and polygalacturonase; this could be one of the reasons for precisely and highly efficient feeding by the oligophagous bug P. micranthus on M. micrantha.

CONCLUSION

Collectively, this work provides a crucial chromosome-level scaffolds resource to study the evolutionary adaptation between mirid bug and their host. It is also helpful in searching for novel environment-friendly biological strategies to control M. micrantha.

Genetic and Transcriptional Regulatory Mechanisms of Lipase Activity in the Plant Pathogenic Fungus Fusarium graminearum

Lipases, which catalyze the hydrolysis of long-chain triglycerides, diglycerides, and monoglycerides into free fatty acids and glycerol, participate in various biological pathways in fungi. In this study, we examined the biological functions and regulatory mechanisms of fungal lipases via two approaches. First, we performed a systemic functional characterization of 86 putative lipase-encoding genes in the plant-pathogenic fungus Fusarium graminearum. The phenotypes were assayed for vegetative growth, asexual and sexual reproduction, stress responses, pathogenicity, mycotoxin production, and lipase activity. Most mutants were normal in the assessed phenotypes, implying overlapping roles for lipases in F. graminearum. In particular, FgLip1 and Fgl1 were revealed as core extracellular lipases in F. graminearum. Second, we examined the lipase activity of previously constructed transcription factor (TF) mutants of F. graminearum and identified three TFs and one histone acetyltransferase that significantly affect lipase activity. The relative transcript levels of FgLIP1 and FGL1 were markedly reduced or enhanced in these TF mutants. Among them, Gzzc258 was identified as a key lipase regulator that is also involved in the induction of lipase activity during sexual reproduction. To our knowledge, this study is the first comprehensive functional analysis of fungal lipases and provides significant insights into the genetic and regulatory mechanisms underlying lipases in fungi. IMPORTANCE Fusarium graminearum is an economically important plant-pathogenic fungus that causes Fusarium head blight (FHB) on wheat and barley. Here, we constructed a gene knockout mutant library of 86 putative lipase-encoding genes and established a comprehensive phenotypic database of the mutants. Among them, we found that FgLip1 and Fgl1 act as core extracellular lipases in this pathogen. Moreover, several putative transcription factors (TFs) that regulate the lipase activities in F. graminearum were identified. The disruption mutants of F. graminearum-lipase regulatory TFs all showed defects in sexual reproduction, which implies a strong relationship between sexual development and lipase activity in this fungus. These findings provide valuable insights into the genetic mechanisms regulating lipase activity as well as its importance to the developmental stages of this plant-pathogenic fungus.

Wnt binding to Coatomer proteins directs secretion on exosomes independently of palmitoylation

Wnt proteins are secreted hydrophobic glycoproteins that act over long distances through poorly understood mechanisms. We discovered that Wnt7a is secreted on extracellular vesicles (EVs) following muscle injury. Structural analysis identified the motif responsible for Wnt7a secretion on EVs that we term the Exosome Binding Peptide (EBP). Addition of the EBP to an unrelated protein directed secretion on EVs. Disruption of palmitoylation, knockdown of WLS, or deletion of the N-terminal signal peptide did not affect Wnt7a secretion on purified EVs. Bio-ID analysis identified Coatomer proteins as candidates responsible for loading Wnt7a onto EVs. The crystal structure of EBP bound to the COPB2 coatomer subunit, the binding thermodynamics, and mutagenesis experiments, together demonstrate that a dilysine motif in the EBP mediates binding to COPB2. Other Wnts contain functionally analogous structural motifs. Mutation of the EBP results in a significant impairment in the ability of Wnt7a to stimulate regeneration, indicating that secretion of Wnt7a on exosomes is critical for normal regeneration in vivo . Our studies have defined the structural mechanism that mediates binding of Wnt7a to exosomes and elucidated the singularity of long-range Wnt signalling.

High and diurnally fluctuating carbon dioxide exposure produces lower mercury toxicity in a marine copepod

Coastal waters have experienced fluctuations in partial pressure of carbon dioxide (pCO₂) and mercury (Hg) pollution, yet little is known concerning how natural pCO₂ fluctuations affect Hg biotoxicity. Here, a marine copepod Tigriopus japonicus was interactively exposed to different seawater pCO₂ (ambient 400, steady elevated 1000, and fluctuating elevated 1000 ± 600 μatm) scenarios and Hg (control, 2 μg/L) treatments for 7 d. The results showed that elevated pCO₂ decreased Hg bioaccumulation, and it was even more under fluctuating elevated pCO₂ condition. We found energy depletion and oxidative stress under Hg-treated copepods, while combined exposure initiated compensatory response to alleviate Hg toxicity. Intriguingly, fluctuating acidification presented more immune defense related genes/processes in Hg-treated copepods when compared to steady acidification, probably linking with the greater decrease in Hg bioaccumulation. Collectively, understanding how fluctuating acidification interacts with Hg contaminant will become more crucial in predicting their risks to coastal biota and ecosystems.

GH20 and GH84 β-N-acetylglucosaminidases with different linkage specificities underpin mucin O-glycan breakdown capability of Bifidobacterium bifidum

Intestinal mucus layers mediate symbiosis and dysbiosis of host-microbe interactions. These interactions are influenced by the mucin O-glycan degrading ability of several gut microbes. The identities and prevalence of many glycoside hydrolyses (GHs) involved in microbial mucin O-glycan breakdown have been previously reported; however, the exact mechanisms and extent to which these GHs are dedicated to mucin O-glycan degradation pathways warrant further research. Here, using Bifidobacterium bifidum as a model mucinolytic bacterium, we revealed that two β-N-acetylglucosaminidases belonging to the GH20 (BbhI) and GH84 (BbhIV) families play important roles in mucin O-glycan degradation. Using substrate specificity analysis of natural oligosaccharides and O-glycomic analysis of porcine gastric mucin (PGM) incubated with purified enzymes or B. bifidum carrying bbhI and/or bbhIV mutations, we showed that BbhI and BbhIV are highly specific for β-(1→3)- and β-(1→6)-GlcNAc linkages of mucin core structures, respectively. Interestingly, we found that efficient hydrolysis of the β-(1→3)-linkage by BbhI of the mucin core 4 structure [GlcNAcβ1-3(GlcNAcβ1-6)GalNAcα-O-Thr] required prior removal of the β-(1→6)-GlcNAc linkage by BbhIV. Consistent with this, inactivation of bbhIV markedly decreased the ability of B. bifidum to release GlcNAc from PGM. When combined with a bbhI mutation, we observed that the growth of the strain on PGM was reduced. Finally, phylogenetic analysis suggests that GH84 members may have gained diversified functions through microbe-microbe and host-microbe horizontal gene transfer events. Taken together, these data strongly suggest GH84 family members in host glycan breakdown.

Identification and expression of the CCO family during development, ripening and stress response in banana

Plant hormone abscisic acid (ABA) plays an important role in plant growth, development and response to biotic / abiotic stressors. Thus, it is necessary to investigate the crucial genes associated with ABA synthesis. Currently, the carotenoid cleavage oxygenases (CCOs) family that function as the key step for ABA synthesis are not well understood in banana. In this study, 13 MaCCO genes and 12 MbCCO genes, divided into NCED subgroup and CCD subgroup, were identified from the banana genome, and their evolutionary relationship, protein motifs, and gene structures were also determined. Transcriptomic analysis suggested the involvement of CCO genes in banana development, ripening, and response to abiotic and biotic stressors, and homologous gene pairs showed homoeologue expression bias in the A or B subgenome. Our results identified MaNCED3A, MaCCD1, and MbNCED3B as the genes with the highest expression during fruit development and ripening. MaNCED5 / MbNCED5 and MaNCED9A might respond to abiotic stress, and MaNCED3A, 3B, 6 A, 9 A, and MbNCED9A showed transcriptional changes that could be a response to Foc4 infection. These findings may contribute to the characterization of key enzymes involved in ABA biosynthesis, as well as to identify potential targets for the genetic improvement of banana.

Streptomonospora mangrovi sp. nov., isolated from mangrove soil showing similar metabolic capabilities, but distinct secondary metabolites profiles

A novel actinomycete, designated strain S1-112 ^T, was isolated from a mangrove soil sample from Hainan, China, and characterized using a polyphasic approach. Strain S1-112 ^T showed the highest similarity of the 16S rRNA gene to Streptomonospora nanhaiensis 12A09^T (99.24%). Their close relationship was further supported by phylogenetic analyses, which placed these two strains within a stable clade. The highest values of digital DNA-DNA hybridization (dDDH, 41.4%) and average nucleotide identity (ANI, 90.55%) were detected between strain S1-112 ^T and Streptomonospora halotolerans NEAU-Jh2-17 ^T. Genotypic and phenotypic characteristics demonstrated that strain S1-112 ^T could be distinguished from its closely related relatives. We also profiled the pan-genome and metabolic features of genomic assemblies of strains belonging to the genus Streptomonospora, indicating similar functional capacities and metabolic activities. However, all of these strains showed promising potential for producing diverse types of secondary metabolites. In conclusion, strain S1-112 ^T represents a novel species of the genus Streptomonospora, for which the name Streptomonospora mangrovi sp. nov. was proposed. The type strain is S1-112 ^T (= JCM 34292 ^T).

BysR, a LysR-Type Pleiotropic Regulator, Controls Production of Occidiofungin by Activating the LuxR-Type Transcriptional Regulator AmbR1 in Burkholderia sp. Strain JP2-270

Occidiofungin is a highly effective antifungal glycopeptide produced by certain Burkholderia strains. The ocf gene cluster, responsible for occidiofungin biosynthesis, is regulated by the cluster-specific regulators encoded by an ambR homolog(s) within the same gene cluster, while the extent to which occidiofungin biosynthesis is connected with the core regulation network remains unknown. Here, we report that the LysR-type regulator BysR acts as a pleiotropic regulator and is essential for occidiofungin biosynthesis. Magnaporthe oryzae was used as an antifungal target in this study, and deletion of bysR and ocfE abolished the antagonistic activity against M. oryzae in Burkholderia sp. strain JP2-270. The ΔbysR defect can be recovered by constitutively expressing bysR or ambR1, but not ambR2. Electrophoretic mobility shift assays (EMSAs) collectively showed that BysR regulates ambR1 by directly binding to its promoter region. In addition, transcriptomic analysis revealed altered expression of 350 genes in response to bysR deletion, and the genes engaged in flagellar assembly and bacterial chemotaxis constitute the most enriched pathways. Also, 400 putative BysR-targeted loci were identified by DNA affinity purification sequencing (DAP-seq) in JP2-270. These loci include not only genes engaged in key metabolic pathways but also those involved in secondary metabolic pathways. To conclude, the occidiofungin produced by JP2-270 is the main substance inhibiting M. oryzae, and BysR controls occidiofungin production by directly targeting ambR1, an intracluster transcriptional regulatory gene that further activates the transcription of the ocf gene cluster. IMPORTANCE We report for the first time that occidiofungin production is regulated by the global transcriptional factor BysR, by directly targeting the specific regulator ambR1, which further promotes the transcription of ocf genes. BysR also acts as a pleiotropic regulator that controls various cellular processes in Burkholderia sp. strain JP2-270. This study provides insight into the regulatory mechanism of occidiofungin synthesis and enhances our understanding of the regulatory patterns of the LysR-type regulator.

Functional Protein Dynamics Directly from Sequences

The sequence correlations within a protein multiple sequence alignment are routinely being used to predict contacts within its structure, but here we point out that these data can also be used to predict a protein's dynamics directly. The elastic network protein dynamics models rely directly upon the contacts, and the normal modes of motion are obtained from the decomposition of the inverse of the contact map. To make the direct connection between sequence and dynamics, it is necessary to apply coarse-graining to the structure at the level of one point per amino acid, which has often been done, and protein coarse-grained dynamics from elastic network models has been highly successful, particularly in representing the large-scale motions of proteins that usually relate closely to their functions. The interesting implication of this is that it is not necessary to know the structure itself to obtain its dynamics and instead to use the sequence information directly to obtain the dynamics.

Pan-genome and transcriptome analyses provide insights into genomic variation and differential gene expression profiles related to disease resistance and fatty acid biosynthesis in eastern black walnut (Juglans nigra)

Walnut (Juglans) species are used as nut crops worldwide. Eastern black walnut (EBW, Juglans nigra), a diploid, horticultural important woody species is native to much of eastern North America. Although it is highly valued for its wood and nut, there are few resources for understanding EBW genetics. Here, we present a high-quality genome assembly of J. nigra based on Illumina, Pacbio, and Hi-C technologies. The genome size was 540.8 Mb, with a scaffold N50 size of 35.1 Mb, and 99.0% of the assembly was anchored to 16 chromosomes. Using this genome as a reference, the resequencing of 74 accessions revealed the effective population size of J. nigra declined during the glacial maximum. A single whole-genome duplication event was identified in the J. nigra genome. Large syntenic blocks among J. nigra, Juglans regia, and Juglans microcarpa predominated, but inversions of more than 600 kb were identified. By comparing the EBW genome with those of J. regia and J. microcarpa, we detected InDel sizes of 34.9 Mb in J. regia and 18.3 Mb in J. microcarpa, respectively. Transcriptomic analysis of differentially expressed genes identified five presumed NBS-LRR (NUCLEOTIDE BINDING SITE-LEUCINE-RICH REPEAT) genes were upregulated during the development of walnut husks and shells compared to developing embryos. We also identified candidate genes with essential roles in seed oil synthesis, including FAD (FATTY ACID DESATURASE) and OLE (OLEOSIN). Our work advances the understanding of fatty acid bioaccumulation and disease resistance in nut crops, and also provides an essential resource for conducting genomics-enabled breeding in walnut.

Whole-genome sequencing and phylogenomic analyses of a novel zearalenone-degrading Bacillus subtilis B72

Bacillus strain B72 was previously isolated as a novel zearalenone (ZEN) degradation strain from the oil field soil in Xinjiang, China. The genome of B72 was sequenced with a 400 bp paired-end using the Illumina HiSeq X Ten platform. De novo genome assembly was performed using SOAPdenovo2 assemblers. Phylogenetic analysis using 16S rRNA gene sequencing demonstrated that B72 is closely related to the novel Bacillus subtilis (B. subtilis) strain DSM 10. A phylogenetic tree based on 31 housekeeping genes, constructed with 19 strains closest at the species level, showed that B72 was closely related to B. subtilis 168, B. licheniformis PT-9, and B. tequilensis KCTC 13622. Detailed phylogenomic analysis using average nucleotide identity (ANI) and genome-to-genome distance calculator (GGDC) demonstrated that B72 might be classified as a novel B. subtilis strain. Our study demonstrated that B72 could degrade 100% of ZEN in minimal medium after 8 h of incubation, which makes it the fastest degrading strain to date. Moreover, we confirmed that ZEN degradation by B72 might involve degrading enzymes produced during the initial period of bacterial growth. Subsequently, functional genome annotation revealed that the laccase-encoding genes yfiH (gene 1743) and cotA (gene 2671) might be related to ZEN degradation in B72. The genome sequence of B. subtilis B72 reported here will provide a reference for genomic research on ZEN degradation in the field of food and feed.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03517-y.

Metatranscriptomics and metabarcoding reveal spatiotemporal shifts in fungal communities and their activities in Chinese coastal waters

Fungal communities are diverse and abundant in coastal waters, yet, their ecological roles and adaptations remain largely unknown. To address these gaps, ITS2 metabarcoding and metatranscriptomic analyses were used to capture the whole suite of fungal diversity and their metabolic potential in water column and sediments in the Yellow Sea during August and October 2019. ITS2 metabarcoding described successfully the abundance of Dikarya during August and October at the different examined habitats, but strongly underrepresented or failed to identify other fungal taxa, including zoosporic and early-diverging lineages, that were abundant in the mycobiome as uncovered by metatranscriptomes. Metatranscriptomics also revealed enriched expression of genes annotated to zoosporic fungi (e.g., chytrids) mainly in the surface water column in October. This enriched expression was correlated with the two-fold increase in chlorophyll-a intensity attributed to phytoplanktonic species which are known to be parasitized by chytrids. The concurrent high expression of genes related to calcium signalling and GTPase activity suggested that these metabolic traits facilitate the parasitic lifestyle of chytrids. Similarly, elevated expression of phagosome genes annotated to Rozellomycota, an early-diverging fungal phylum not fully detected with ITS2 metabarcoding, suggested that this taxon utilizes a suite of feeding modes, including phagotrophy in this coastal setting. Our data highlight the necessity of using combined approaches to accurately describe the community structure of coastal mycobiome. We also provide in-depth insights into the fungal ecological roles in coastal waters, and report potential metabolic mechanisms utilized by fungi to cope with environmental stresses that occur during distinct seasonal months in coastal ecosystems.

Genome-Wide Identification, Characterization and Experimental Expression Analysis of CNGC Gene Family in Gossypium

Cyclic nucleotide-gated ion channels (CNGCs) are channel proteins for calcium ions, and have been reported to play important roles in regulating survival and environmental response of various plants. However, little is known about how the CNGC family works in Gossypium. In this study, 173 CNGC genes, which were identified from two diploid and five tetraploid Gossypium species, were classified into four groups by phylogenetic analysis. The collinearity results demonstrated that CNGC genes are integrally conservative among Gossypium species, but four gene losses and three simple translocations were detected, which is beneficial to analyzing the evolution of CNGCs in Gossypium. The various cis-acting regulatory elements in the CNGCs' upstream sequences revealed their possible functions in responding to multiple stimuli such as hormonal changes and abiotic stresses. In addition, expression levels of 14 CNGC genes changed significantly after being treated with various hormones. The findings in this study will contribute to understanding the function of the CNGC family in cotton, and lay a foundation for unraveling the molecular mechanism of cotton plants' response to hormonal changes.

Insights into the Bioinformatics and Transcriptional Analysis of the Elongator Complexes (ELPs) Gene Family of Wheat: TaELPs Contribute to Wheat Abiotic Stress Tolerance and Leaf Senescence

Elongator complexes (ELPs) are the protein complexes that promote transcription through histone acetylation in eukaryotic cells and interact with elongating RNA polymerase II (RNAPII). ELPs' role in plant growth and development, signal transduction, and response to biotic and abiotic stresses have been confirmed in model plants. However, the functions of the wheat ELP genes are not well documented. The present study identified 18 members of the ELPs from the wheat genome with a homology search. Further, bioinformatics and transcription patterns in response to different stress conditions were analyzed to dissect their potential regulatory mechanisms in wheat. Gene duplication analysis showed that 18 pairs of ELP paralogous genes were derived from segmental duplication, which was divided into six clades by protein phylogenetic and cluster analysis. The orthologous analysis of wheat TaELP genes showed that TaELP genes may have evolved from orthologous genes of other plant species or closely related plants. Moreover, a variety of cis-acting regulatory elements (CAREs) related to growth and development, hormone response, and biotic and abiotic stresses were identified in the TaELPs' promoter regions. The qRT-PCR analysis showed that the transcription of TaELPs was induced under hormone, salt, and drought stress and during leaf senescence. The TaELP2 gene was silenced with BSMV-VIGS, and TaELP2 was preliminarily verified to be involved in the regulation of wheat leaf senescence. Overall, TaELP genes might be regulated by hormone signaling pathways and response to abiotic stress and leaf senescence, which could be investigated further as potential candidate genes for wheat abiotic stress tolerance and yield improvement.

Leaves of Moringa oleifera Are Potential Source of Bioactive Compound β-Carotene: Evidence from In Silico and Quantitative Gene Expression Analysis

Moringa oleifera is rich in bioactive compounds such as beta-carotene, which have high nutritional values and antimicrobial applications. Several studies have confirmed that bioactive-compound-based herbal medicines extracted from the leaves, seeds, fruits and shoots of M. oleifera are vital to cure many diseases and infections, and for the healing of wounds. The β-carotene is a naturally occurring bioactive compound encoded by zeta-carotene desaturase (ZDS) and phytoene synthase (PSY) genes. In the current study, computational analyses were performed to identify and characterize ZDS and PSY genes retrieved from Arabidopsis thaliana (as reference) and these were compared with the corresponding genes in M. oleifera, Brassica napus, Brassica rapa, Brassica oleracea and Bixa orellana. The BLAST results revealed that all the plant species considered in this study encode β-carotene genes with 80-100% similarity. The Pfam analysis on β-carotene genes of all the investigated plants confirmed that they belong to the same protein family and domain. Similarly, phylogenetic analysis revealed that β-carotene genes of M. oleifera belong to the same ancestral class. Using the ZDS and PSY genes of Arabidopsis thaliana as a reference, we conducted qRT-PCR analysis on RNA extracted from the leaves of M. oleifera, Brassica napus, Brassica rapa and Bixa orellana. It was noted that the most significant gene expression occurred in the leaves of the studied medicinal plants. We concluded that not only are the leaves of M. oleifera an effective source of bioactive compounds including beta carotene, but also the leaves of Brassica napus, Brassica rapa and Bixa orellana can be employed as antibiotics and antioxidants against bacterial or microbial infections.

Two haplotype-resolved, gap-free genome assemblies for Actinidia latifolia and Actinidia chinensis shed light on the regulatory mechanisms of vitamin C and sucrose metabolism in kiwifruit

Kiwifruit is a recently domesticated horticultural fruit crop with substantial economic and nutritional value, especially because of the high content of vitamin C in its fruit. In this study, we de novo assembled two telomere-to-telomere kiwifruit genomes from Actinidia chinensis var. 'Donghong' (DH) and Actinidia latifolia 'Kuoye' (KY), with total lengths of 608 327 852 and 640 561 626 bp for 29 chromosomes, respectively. With a burst of structural variants involving inversion, translocations, and duplications within 8.39 million years, the metabolite content of DH and KY exhibited differences in saccharides, lignans, and vitamins. A regulatory ERF098 transcription factor family has expanded in KY and Actinidia eriantha, both of which have ultra-high vitamin C content. With each assembly phased into two complete haplotypes, we identified allelic variations between two sets of haplotypes, leading to protein sequence variations in 26 494 and 27 773 gene loci and allele-specific expression of 4687 and 12 238 homozygous gene pairs. Synchronized metabolome and transcriptome changes during DH fruit development revealed the same dynamic patterns in expression levels and metabolite contents; free fatty acids and flavonols accumulated in the early stages, but sugar substances and amino acids accumulated in the late stages. The AcSWEET9b gene that exhibits allelic dominance was further identified to positively correlate with high sucrose content in fruit. Compared with wild varieties and other Actinidia species, AcSWEET9b promoters were selected in red-flesh kiwifruits that have increased fruit sucrose content, providing a possible explanation on why red-flesh kiwifruits are sweeter. Collectively, these two gap-free kiwifruit genomes provide a valuable genetic resource for investigating domestication mechanisms and genome-based breeding of kiwifruit.

Engineering processive cellulase of Clostridium thermocellum to divulge the role of the carbohydrate-binding module

The processive cellulase (CelO) is an important modular enzyme of Clostridium thermocellum. To study the effect of the carbohydrate-binding module (CBM3b) on the catalytic domain of CelO (GH5), four engineered derivatives of CelO were designed by truncation and terminal fusion of CBM3b. These are CBM at the N-terminus, native form (CelO-BC, 62 kDa); catalytic domain only (CelO-C, 42 kDa); CBM at the C-terminus (CelO-CB, 54 kDa) and CBM attached at both termini (CelO-BCB, 73 kDa). All constructs were cloned into pET22b (+) and expressed in Escherichia coli BL21 (DE3) star. The expression levels of CelO-C, CelO-CB, CelO-BC, and CelO-BCB were 35%, 35%, 30%, and 20%, respectively. The enzyme activities of CelO-C, CelO-CB, CelO-BC, and CelO-BCB against 1% regenerated amorphous cellulose (RAC) were 860, 758, 985, and 1208 units per μmole of the enzyme, respectively. The enzymes were partially purified from the lysate of E. coli cells by heat treatment followed by anion exchange FPLC purification. Against RAC, CelO-C, CelO-CB, CelO-BC, and CelO-BCB showed K_M values of 32, 33, 45, and 43 mg⋅mL^-1 and V_max values of 3571, 3846, 3571, and 4545 U⋅min^-1 , respectively. CBM3b at the N-terminus of GH5 linked through a P/T-rich linker was found to enhance the catalytic activity and thermostability of the enzyme.

Enzyme redesign and genetic code expansion

Enzyme design is an important application of computational protein design (CPD). It can benefit enormously from the additional chemistries provided by noncanonical amino acids (ncAAs). These can be incorporated into an 'expanded' genetic code, and introduced in vivo into target proteins. The key step for genetic code expansion is to engineer an aminoacyl-transfer RNA (tRNA) synthetase (aaRS) and an associated tRNA that handles the ncAA. Experimental directed evolution has been successfully used to engineer aaRSs and incorporate over 200 ncAAs into expanded codes. But directed evolution has severe limits, and is not yet applicable to noncanonical AA backbones. CPD can help address several of its limitations, and has begun to be applied to this problem. We review efforts to redesign aaRSs, studies that designed new proteins and functionalities with the help of ncAAs, and some of the method developments that have been used, such as adaptive landscape flattening Monte Carlo, which allows an enzyme to be redesigned with substrate or transition state binding as the design target.

Comprehensive Analysis of Whole-Transcriptome Profiles in Response to Acute Hypersaline Challenge in Chinese Razor Clam Sinonovacula constricta

The Chinese razor clam (Sinonovacula constricta) is an important for Chinese aquaculture marine bivalve that naturally occurs across intertidal and estuarine areas subjected to significant changes in salinity level. However, the information on the molecular mechanisms related to high salinity stress in the species remain limited. In this study, nine gill samples of S. constricta treated with 20, 30, and 40 ppt salinity for 24 h were used for whole-transcriptome RNA sequencing, and a regulatory network of competing endogenous RNAs (ceRNAs) was constructed to better understand the mechanisms responsible for adaptation of the species to high salinity. A total of 83,262 lncRNAs, 52,422 mRNAs, 2890 circRNAs, and 498 miRNAs were identified, and 4175 of them displayed differential expression pattern among the three groups examined. The KEGG analyses of differentially expressed RNAs evidenced that amino acid synthesis and membrane transport were the dominant factors involved in the adaptation of the Chinese razor clam to acute salinity increase, while lipid metabolism and signaling played only a supporting role. In addition, lncRNA/circRNA-miRNA-mRNA regulatory networks (ceRNA network) showed clearly regulatory relationships among different RNAs. Moreover, the expression of four candidate genes, including tyrosine aminotransferase (TAT), hyaluronidase 4 (HYAL4), cysteine sulfinic acid decarboxylase (CSAD), and ∆¹-pyrroline-5-carboxylate synthase (P5CS) at different challenge time were detected by qRT-PCR. The expression trend of TAT and HYAL4 was consistent with that of the ceRNA network, supporting the reliability of established network. The expression of TAT, CSAD, and P5CS were upregulated in response to increased salinity. This might be associated with increased amino acid synthesis rate, which seems to play an essential role in adaptation of the species to high salinity stress. In contrast, the expression level of HYAL4 gene decreased in response to elevated salinity level, which is associated with reduction Hyaluronan hydrolysis to help maintain water in the cell. Our findings provide a very rich reference for understanding the important role of ncRNAs in the salinity adaptation of shellfish. Moreover, the acquired information may be useful for optimization of the artificial breeding of the Chinese razor clam under aquaculture conditions.