The Pfam protein families database

Genomic and epigenomic insights into the mechanism of cold response in upland cotton (Gossypium hirsutum)

Functional genome research, including gene transcriptional and posttranslational modifications of histones, can benefit greatly from a high-quality genome assembly. Histone modification plays a significant role in modulating the responses to abiotic stress in plants. However, there are limited reports on the involvement of dynamic changes in histone modification in cold stress response in upland cotton. In this study, the genome of an elite accession, YM11, with considerable cold stress tolerance was de novo assembled, which yielded a genome of 2343.06 Mb with a contig N50 of 88.96 Mb, and a total of 73,821 protein-coding gene models were annotated. Comparisons among YM11 and five Gossypium allopolyploid cotton assemblies highlighted a large amount of structural variations and presence/absence variations. We analyzed transcriptome and metabolome changes in YM11 seedlings subjected to cold stress. Using the CUT&Tag method, genome-wide H3K4me3 and H3K9ac modification patterns and effect of histone changes on gene expression were profiled during cold stress. Significant and consistently changing histone modifications and the gene expressions were screened, of which transcription factors (TFs) were highlighted. Our results suggest a positive correlation between the changes in H3K4me3, H3K9ac modifications and cold stress-responsive gene activation. This genome assembly and comprehensive analysis of genome-wide histone modifications and gene expression provide insights into the genomic variation and epigenetic responses to cold stress in upland cotton.

Putative novel hydrogen- and iron-oxidizing sheath-producing Zetaproteobacteria thrive at the Fåvne deep-sea hydrothermal vent field

IMPORTANCE

Knowledge on microbial iron oxidation is important for understanding the cycling of iron, carbon, nitrogen, nutrients, and metals. The current study yields important insights into the niche sharing, diversification, and Fe(III) oxyhydroxide morphology of Ghiorsea, an iron- and hydrogen-oxidizing Zetaproteobacteria representative belonging to Zetaproteobacteria operational taxonomic unit 9. The study proposes that Ghiorsea exhibits a more extensive morphology of Fe(III) oxyhydroxide than previously observed. Overall, the results increase our knowledge on potential drivers of Zetaproteobacteria diversity in iron microbial mats and can eventually be used to develop strategies for the cultivation of sheath-forming Zetaproteobacteria.

Taxonomic distribution of metabolic functions in bacteria associated with Trichodesmium consortia

IMPORTANCE

Colonies of the cyanobacteria Trichodesmium act as a biological hotspot for the usage and recycling of key resources such as C, N, P, and Fe within an otherwise oligotrophic environment. While Trichodesmium colonies are known to interact and support a unique community of algae and particle-associated microbes, our understanding of the taxa that populate these colonies and the gene functions they encode is still limited. Characterizing the taxa and adaptive strategies that influence consortium physiology and its concomitant biogeochemistry is critical in a future ocean predicted to have increasingly resource-depleted regions.

Probiotic Potential and Safety Assessment of Lactiplantibacillus plantarum cqf-43 and Whole-Genome Sequence Analysis

This study reports the whole-genome sequence of Lactiplantibacillus plantarum cqf-43 isolated from healthy sow feces. Based on genomic analysis, we performed a comprehensive safety assessment of strain cqf-43, using both in vitro and in vivo experiments, and explored its probiotic potential. The total genome length measures 3,169,201 bp, boasting a GC content of 44.59%. Through phylogenetic analyses, leveraging both 16S rRNA gene and whole-genome sequences, we confidently categorize strain cqf-43 as a member of Lactiplantibacillus. Genome annotation using Prokka unveiled a total of 3141 genes, encompassing 2990 protein-coding sequences, 71 tRNAs, 16 rRNAs, and 1 tmRNA. Functional annotations derived from COG and KEGG databases highlighted a significant abundance of genes related to metabolism, with a notable emphasis on carbohydrate utilization. The genome also revealed the presence of prophage regions and CRISPR-Cas regions while lacking virulence and toxin genes. Screening for antibiotic resistance genes via the CARD database yielded no detectable transferable resistance genes, effectively eliminating the potential for harmful gene transfer. It is worth highlighting that the virulence factors identified via the VFDB database primarily contribute to bolstering pathogen resilience in hostile environments. This characteristic is particularly advantageous for probiotics. Furthermore, the genome is devoid of menacing genes such as hemolysin, gelatinase, and biogenic amine-producing genes. Our investigation also unveiled the presence of three unannotated secondary metabolite biosynthetic gene clusters, as detected by the online tool antiSMASH, suggesting a great deal of unknown potential for this strain. Rigorous in vitro experiments confirmed tolerance of strain cqf-43 in the intestinal environment, its antimicrobial efficacy, sensitivity to antibiotics, absence of hemolysis and gelatinase activity, and its inability to produce biogenic amines. In addition, a 28-day oral toxicity test showed that the strain cqf-43 did not pose a health hazard in mice, further establishing it as a safe strain.

A reference transcriptome for walnut anthracnose pathogen, Ophiognomonia leptostyla, guides the discovery of candidate virulence genes

Despite the economic losses due to the walnut anthracnose, Ophiognomonia leptostyla is an orphan fungus with respect to genomic resources. In the present study, the transcriptome of O. leptostyla was assembled for the first time. RNA sequencing was conducted for the fungal mycelia grown in a liquid media, and the inoculated leaf samples of walnut with the fungal conidia sampled at 48, 96 and 144 h post inoculation (hpi). The completeness, correctness, and contiguity of the de novo transcriptome assemblies generated with Trinity, Oases, SOAPdenovo-Trans and Bridger were compared to identify a single superior reference assembly. In most of the assessment criteria including N50, Transrate score, number of ORFs with known description in gene bank, the percentage of reads mapped back to the transcript (RMBT), BUSCO score, Swiss-Prot coverage bin and RESM-EVAL score, the Bridger assembly was the superior and thus used as a reference for profiling the O. leptostyla transcriptome in liquid media vs. during walnut infection. The k-means clustering of transcripts resulted in four distinct transcription patterns across the three sampling time points. Most of the detected CAZy transcripts had elevated transcription at 96 hpi that is hypothetically concurrent with the start of intracellular growth. The in-silico analysis revealed 103 candidate effectors of which six were members of Necrosis and Ethylene Inducing Like Protein (NLP) gene family belonging to three distinct k-means clusters. This study provided a complex and temporal pattern of the CAZys and candidate effectors transcription during six days post O. leptostyla inoculation on walnut leaves, introducing a list of candidate virulence genes for validation in future studies.

Comparative Transcriptomic Assessment of Chemosensory Genes in Adult and Larval Olfactory Organs of Cnaphalocrocis medinalis

The rice leaf folder, Cnaphalocrocis medinalis (Lepidoptera: Pyralidae), is a notorious pest of rice in Asia. The larvae and adults of C. medinalis utilize specialized chemosensory systems to adapt to different environmental odors and physiological behaviors. However, the differences in chemosensory genes between the olfactory organs of these two different developmental stages remain unclear. Here, we conducted a transcriptome analysis of larvae heads, male antennae, and female antennae in C. medinalis and identified 131 putative chemosensory genes, including 32 OBPs (8 novel OBPs), 23 CSPs (2 novel CSPs), 55 ORs (17 novel ORs), 19 IRs (5 novel IRs) and 2 SNMPs. Comparisons between larvae and adults of C. medinalis by transcriptome and RT-qPCR analysis revealed that the number and expression of chemosensory genes in larval heads were less than that of adult antennae. Only 17 chemosensory genes (7 OBPs and 10 CSPs) were specifically or preferentially expressed in the larval heads, while a total of 101 chemosensory genes (21 OBPs, 9 CSPs, 51 ORs, 18 IRs, and 2 SNMPs) were specifically or preferentially expressed in adult antennae. Our study found differences in chemosensory gene expression between larvae and adults, suggesting their specialized functions at different developmental stages of C. medinalis. These results provide a theoretical basis for screening chemosensory genes as potential molecular targets and developing novel management strategies to control C. medinalis.

Chromosome-scale genomes of commercially important mahoganies, Swietenia macrophylla and Khaya senegalensis

Mahogany species (family Meliaceae) are highly valued for their aesthetic and durable wood. Despite their economic and ecological importance, genomic resources for mahogany species are limited, hindering genetic improvement and conservation efforts. Here we perform chromosome-scale genome assemblies of two commercially important mahogany species: Swietenia macrophylla and Khaya senegalensis. By combining 10X sequencing and Hi-C data, we assemble high-quality genomes of 274.49 Mb (S. macrophylla) and 406.50 Mb (K. senegalensis), with scaffold N50 lengths of 8.51 Mb and 7.85 Mb, respectively. A total of 99.38% and 98.05% of the assembled sequences are anchored to 28 pseudo-chromosomes in S. macrophylla and K. senegalensis, respectively. We predict 34,129 and 31,908 protein-coding genes in S. macrophylla and K. senegalensis, respectively, of which 97.44% and 98.49% are functionally annotated. The chromosome-scale genome assemblies of these mahogany species could serve as a vital genetic resource, especially in understanding the properties of non-model woody plants. These high-quality genomes could support the development of molecular markers for breeding programs, conservation efforts, and the sustainable management of these valuable forest resources.

Ensemble Learning with Supervised Methods Based on Large-Scale Protein Language Models for Protein Mutation Effects Prediction

Machine learning has been increasingly utilized in the field of protein engineering, and research directed at predicting the effects of protein mutations has attracted increasing attention. Among them, so far, the best results have been achieved by related methods based on protein language models, which are trained on a large number of unlabeled protein sequences to capture the generally hidden evolutionary rules in protein sequences, and are therefore able to predict their fitness from protein sequences. Although numerous similar models and methods have been successfully employed in practical protein engineering processes, the majority of the studies have been limited to how to construct more complex language models to capture richer protein sequence feature information and utilize this feature information for unsupervised protein fitness prediction. There remains considerable untapped potential in these developed models, such as whether the prediction performance can be further improved by integrating different models to further improve the accuracy of prediction. Furthermore, how to utilize large-scale models for prediction methods of mutational effects on quantifiable properties of proteins due to the nonlinear relationship between protein fitness and the quantification of specific functionalities has yet to be explored thoroughly. In this study, we propose an ensemble learning approach for predicting mutational effects of proteins integrating protein sequence features extracted from multiple large protein language models, as well as evolutionarily coupled features extracted in homologous sequences, while comparing the differences between linear regression and deep learning models in mapping these features to quantifiable functional changes. We tested our approach on a dataset of 17 protein deep mutation scans and indicated that the integrated approach together with linear regression enables the models to have higher prediction accuracy and generalization. Moreover, we further illustrated the reliability of the integrated approach by exploring the differences in the predictive performance of the models across species and protein sequence lengths, as well as by visualizing clustering of ensemble and non-ensemble features.

Genome sequencing of Syzygium cumini (jamun) reveals adaptive evolution in secondary metabolism pathways associated with its medicinal properties

Syzygium cumini, also known as jambolan or jamun, is an evergreen tree widely known for its medicinal properties, fruits, and ornamental value. To understand the genomic and evolutionary basis of its medicinal properties, we sequenced S. cumini genome for the first time from the world's largest tree genus Syzygium using Oxford Nanopore and 10x Genomics sequencing technologies. We also sequenced and assembled the transcriptome of S. cumini in this study. The tetraploid and highly heterozygous draft genome of S. cumini had a total size of 709.9 Mbp with 61,195 coding genes. The phylogenetic position of S. cumini was established using a comprehensive genome-wide analysis including species from 18 Eudicot plant orders. The existence of neopolyploidy in S. cumini was evident from the higher number of coding genes and expanded gene families resulting from gene duplication events compared to the other two sequenced species from this genus. Comparative evolutionary analyses showed the adaptive evolution of genes involved in the phenylpropanoid-flavonoid (PF) biosynthesis pathway and other secondary metabolites biosynthesis such as terpenoid and alkaloid in S. cumini, along with genes involved in stress tolerance mechanisms, which was also supported by leaf transcriptome data generated in this study. The adaptive evolution of secondary metabolism pathways is associated with the wide range of pharmacological properties, specifically the anti-diabetic property, of this species conferred by the bioactive compounds that act as nutraceutical agents in modern medicine.

ATP synthase evolution on a cross-braced dated tree of life

The timing of early cellular evolution, from the divergence of Archaea and Bacteria to the origin of eukaryotes, is poorly constrained. The ATP synthase complex is thought to have originated prior to the Last Universal Common Ancestor (LUCA) and analyses of ATP synthase genes, together with ribosomes, have played a key role in inferring and rooting the tree of life. We reconstruct the evolutionary history of ATP synthases using an expanded taxon sampling set and develop a phylogenetic cross-bracing approach, constraining equivalent speciation nodes to be contemporaneous, based on the phylogenetic imprint of endosymbioses and ancient gene duplications. This approach results in a highly resolved, dated species tree and establishes an absolute timeline for ATP synthase evolution. Our analyses show that the divergence of ATP synthase into F- and A/V-type lineages was a very early event in cellular evolution dating back to more than 4 Ga, potentially predating the diversification of Archaea and Bacteria. Our cross-braced, dated tree of life also provides insight into more recent evolutionary transitions including eukaryogenesis, showing that the eukaryotic nuclear and mitochondrial lineages diverged from their closest archaeal (2.67-2.19 Ga) and bacterial (2.58-2.12 Ga) relatives at approximately the same time, with a slightly longer nuclear stem-lineage.

Chromosome-level genome assembly of the deep-sea snail Phymorhynchus buccinoides provides insights into the adaptation to the cold seep habitat

BACKGROUND

The deep-sea snail Phymorhynchus buccinoides belongs to the genus Phymorhynchus (Neogastropoda: Raphitomidae), and it is a dominant specie in the cold seep habitat. As the environment of the cold seep is characterized by darkness, hypoxia and high concentrations of toxic substances such as hydrogen sulfide (H2S), exploration of the diverse fauna living around cold seeps will help to uncover the adaptive mechanisms to this unique habitat. In the present study, a chromosome-level genome of P. buccinoides was constructed and a series of genomic and transcriptomic analyses were conducted to explore its molecular adaptation mechanisms to the cold seep environments.

RESULTS

The assembled genome size of the P. buccinoides was approximately 2.1 Gb, which is larger than most of the reported snail genomes, possibly due to the high proportion of repetitive elements. About 92.0% of the assembled base pairs of contigs were anchored to 34 pseudo-chromosomes with a scaffold N50 size of 60.0 Mb. Compared with relative specie in the shallow water, the glutamate regulative and related genes were expanded in P. buccinoides, which contributes to the acclimation to hypoxia and coldness. Besides, the relatively high mRNA expression levels of the olfactory/chemosensory genes in osphradium indicate that P. buccinoides might have evolved a highly developed and sensitive olfactory organ for its orientation and predation. Moreover, the genome and transcriptome analyses demonstrate that P. buccinoides has evolved a sulfite-tolerance mechanism by performing H2S detoxification. Many genes involved in H2S detoxification were highly expressed in ctenidium and hepatopancreas, suggesting that these tissues might be critical for H2S detoxification and sulfite tolerance.

CONCLUSIONS

In summary, our report of this chromosome-level deep-sea snail genome provides a comprehensive genomic basis for the understanding of the adaptation strategy of P. buccinoides to the extreme environment at the deep-sea cold seeps.

Exploring the genetic landscape of nitrogen uptake in durum wheat: genome-wide characterization and expression profiling of NPF and NRT2 gene families

Nitrate uptake by plants primarily relies on two gene families: Nitrate transporter 1/peptide transporter (NPF) and Nitrate transporter 2 (NRT2). Here, we extensively characterized the NPF and NRT2 families in the durum wheat genome, revealing 211 NPF and 20 NRT2 genes. The two families share many Cis Regulatory Elements (CREs) and Transcription Factor binding sites, highlighting a partially overlapping regulatory system and suggesting a coordinated response for nitrate transport and utilization. Analyzing RNA-seq data from 9 tissues and 20 cultivars, we explored expression profiles and co-expression relationships of both gene families. We observed a strong correlation between nucleotide variation and gene expression within the NRT2 gene family, implicating a shared selection mechanism operating on both coding and regulatory regions. Furthermore, NPF genes showed highly tissue-specific expression profiles, while NRT2s were mainly divided in two co-expression modules, one expressed in roots (NAR2/NRT3 dependent) and the other induced in anthers and/ovaries during maturation. Our evidences confirmed that the majority of these genes were retained after small-scale duplication events, suggesting a neo- or sub-functionalization of many NPFs and NRT2s. Altogether, these findings indicate that the expansion of these gene families in durum wheat could provide valuable genetic variability useful to identify NUE-related and candidate genes for future breeding programs in the context of low-impact and sustainable agriculture.

Combination of linkage and association mapping with genomic prediction to infer QTL regions associated with gray leaf spot and northern corn leaf blight resistance in tropical maize

Among the diseases threatening maize production in Africa are gray leaf spot (GLS) caused by Cercospora zeina and northern corn leaf blight (NCLB) caused by Exserohilum turcicum. The two pathogens, which have high genetic diversity, reduce the photosynthesizing ability of susceptible genotypes and, hence, reduce the grain yield. To identify population-based quantitative trait loci (QTLs) for GLS and NCLB resistance, a biparental population of 230 lines derived from the tropical maize parents CML511 and CML546 and an association mapping panel of 239 tropical and sub-tropical inbred lines were phenotyped across multi-environments in western Kenya. Based on 1,264 high-quality polymorphic single-nucleotide polymorphisms (SNPs) in the biparental population, we identified 10 and 18 QTLs, which explained 64.2% and 64.9% of the total phenotypic variance for GLS and NCLB resistance, respectively. A major QTL for GLS, qGLS1_186 accounted for 15.2% of the phenotypic variance, while qNCLB3_50 explained the most phenotypic variance at 8.8% for NCLB resistance. Association mapping with 230,743 markers revealed 11 and 16 SNPs significantly associated with GLS and NCLB resistance, respectively. Several of the SNPs detected in the association panel were co-localized with QTLs identified in the biparental population, suggesting some consistent genomic regions across genetic backgrounds. These would be more relevant to use in field breeding to improve resistance to both diseases. Genomic prediction models trained on the biparental population data yielded average prediction accuracies of 0.66-0.75 for the disease traits when validated in the same population. Applying these prediction models to the association panel produced accuracies of 0.49 and 0.75 for GLS and NCLB, respectively. This research conducted in maize fields relevant to farmers in western Kenya has combined linkage and association mapping to identify new QTLs and confirm previous QTLs for GLS and NCLB resistance. Overall, our findings imply that genetic gain can be improved in maize breeding for resistance to multiple diseases including GLS and NCLB by using genomic selection.

Genomic analysis of the emergent aquatic plant Sparganium stoloniferum provides insights into its clonality, local adaptation and demographic history

Clonal propagation and extensive dispersal of seeds and asexual propagules are two important features of aquatic plants that help them adapt to aquatic environments. Accurate measurements of clonality and effective clonal dispersal are essential for understanding the evolution of aquatic plants. Here, we first assembled a high-quality chromosome-level genome of a widespread emergent aquatic plant Sparganium stoloniferum to provide a reference for its population genomic study. We then performed high-depth resequencing of 173 individuals from 20 populations covering different basins across its range in China. Population genomic analyses revealed three genetic lineages reflecting the northeast (NE), southwest (SW) and northwest (NW) of its geographical distribution. The NE lineage diverged in the middle Pleistocene while the SW and NW lineages diverged until about 2400 years ago. Clonal relationship analyses identified nine populations as monoclonal population. Dispersal of vegetative propagules was identified between five populations covering three basins in the NE lineage, and dispersal distance was up to 1041 km, indicating high dispersibility in emergent aquatic plant species. We also identified lineage-specific positively selected genes that are likely to be involved in adaptations to saline wetlands and high-altitude environments. Our findings accurately measure the clonality, determine the dispersal range and frequency of vegetative propagules, and detect genetic signatures of local adaptation in a widespread emergent aquatic plant species, providing new perspectives on the evolution of aquatic plants.

Comparative genomics profiling revealed multi-stress responsive roles of the CC-NBS-LRR genes in three mango cultivars

The nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family is the largest group of disease resistance (R) genes in plants and is active in response to viruses, bacteria, and fungi usually involved in effector-triggered immunity (ETI). Pangenome-wide studies allow researchers to analyze the genetic diversity of multiple species or their members simultaneously, providing a comprehensive understanding of the evolutionary relationships and diversity present among them. The draft pan-genome of three Mangifera indica cultivars (Alphonso, Hong Xiang Ya, and Tommy atkins) was constructed and Presence/absence variants (PAVs) were filtered through the ppsPCP pipeline. As a result, 2823 genes and 5907 PAVs from H. Xiang Ya, and 1266 genes and 2098 PAVs from T. atkins were added to the reference genome. For the identification of CC-NBS-LRR (CNL) genes in these mango cultivars, this draft pan-genome study has successfully identified 47, 27, and 36 members in Alphonso, H. Xiang Ya, and T. atkins respectively. The phylogenetic analysis divided MiCNL proteins into four distinct subgroups. All MiCNL genes are unevenly distributed on chromosomes. Both tandem and segmental duplication events played a significant role in the expansion of the CNL gene family. These genes contain cis-elements related to light, stress, hormone, and development. The analysis of protein-protein interactions (PPI) revealed that MiCNL proteins interacted with other defense-responsive proteins. Gene Ontology (GO) analysis indicated that MiCNL genes play a role in defense mechanisms within the organism. The expression level of the identified genes in fruit peel was observed under disease and cold stress which showed that Mi_A_CNL13 and 14 were up-regulated while Mi_A_CNL15, 25, 30, 31, and 40 were down-regulated in disease stress. On the other hand, Mi_A_CNL2, 14, 41, and 45 were up-regulated and Mi_A_CNL47 is down-regulated in cold stress. Subsequently, the Random Forest (RF) classifier was used to assess the multi-stress response of MiCNLs. It was found that Mi_A_CNL14 is a gene that responds to multiple stress conditions. The CNLs have similar protein structures which show that they are involved in the same function. The above findings provide a foundation for a deeper understanding of the functional characteristics of the mango CNL gene family.

Genome-wide identification of nucleotide-binding domain leucine-rich repeat (NLR) genes and their association with green peach aphid (Myzus persicae) resistance in peach

Resistance genes (R genes) are a class of genes that are immune to a wide range of diseases and pests. In planta, NLR genes are essential components of the innate immune system. Currently, genes belonging to NLR family have been found in a number of plant species, but little is known in peach. Here, 286 NLR genes were identified on peach genome by using their homologous genes in Arabidopsis thaliana as queries. These 286 NLR genes contained at least one NBS domain and LRR domain. Phylogenetic and N-terminal domain analysis showed that these NLRs could be separated into four subfamilies (I-IV) and their promoters contained many cis-elements in response to defense and phytohormones. In addition, transcriptome analysis showed that 22 NLR genes were up-regulated after infected by Green Peach Aphid (GPA), and showed different expression patterns. This study clarified the NLR gene family and their potential functions in aphid resistance process. The candidate NLR genes might be useful in illustrating the mechanism of aphid resistance in peach.

Unraveling the evolutionary dynamics of the TPS gene family in land plants

Terpenes and terpenoids are key natural compounds for plant defense, development, and composition of plant oil. The synthesis and accumulation of a myriad of volatile terpenoid compounds in these plants may dramatically alter the quality and flavor of the oils, which provide great commercial utilization value for oil-producing plants. Terpene synthases (TPSs) are important enzymes responsible for terpenic diversity. Investigating the differentiation of the TPS gene family could provide valuable theoretical support for the genetic improvement of oil-producing plants. While the origin and function of TPS genes have been extensively studied, the exact origin of the initial gene fusion event - it occurred in plants or microbes - remains uncertain. Furthermore, a comprehensive exploration of the TPS gene differentiation is still pending. Here, phylogenetic analysis revealed that the fusion of the TPS gene likely occurred in the ancestor of land plants, following the acquisition of individual C- and N- terminal domains. Potential mutual transfer of TPS genes was observed among microbes and plants. Gene synteny analysis disclosed a differential divergence pattern between TPS-c and TPS-e/f subfamilies involved in primary metabolism and those (TPS-a/b/d/g/h subfamilies) crucial for secondary metabolites. Biosynthetic gene clusters (BGCs) analysis suggested a correlation between lineage divergence and potential natural selection in structuring terpene diversities. This study provides fresh perspectives on the origin and evolution of the TPS gene family.

Comprehensive Analysis of GH3 Gene Family in Potato and Functional Characterization of StGH3.3 under Drought Stress

As an important hormone response gene, Gretchen Hagen 3 (GH3) maintains hormonal homeostasis by conjugating excess auxin with amino acids during plant stress-related signaling pathways. GH3 genes have been characterized in many plant species, but they are rarely reported in potato. Here, 19 StGH3 genes were isolated and characterized. Phylogenetic analysis indicated that StGH3s were divided into two categories (group I and group III). Analyses of gene structure and motif composition showed that the members of a specific StGH3 subfamily are relatively conserved. Collinearity analysis of StGH3 genes in potato and other plants laid a foundation for further exploring the evolutionary characteristics of the StGH3 genes. Promoter analysis showed that most StGH3 promoters contained hormone and abiotic stress response elements. Multiple transcriptome studies indicated that some StGH3 genes were responsive to ABA, water deficits, and salt treatments. Moreover, qRT-PCR analysis indicated that StGH3 genes could be induced by phytohormones (ABA, SA, and MeJA) and abiotic stresses (water deficit, high salt, and low temperature), although with different patterns. Furthermore, transgenic tobacco with transient overexpression of the StGH3.3 gene showed positive regulation in response to water deficits by increasing proline accumulation and reducing the leaf water loss rate. These results suggested that StGH3 genes may be involved in the response to abiotic stress through hormonal signal pathways. Overall, this study provides useful insights into the evolution and function of StGH3s and lays a foundation for further study on the molecular mechanisms of StGH3s in the regulation of potato drought resistance.

Insights into bacterial community metatranscriptome and metabolome in river water influenced by palm oil mill effluent final discharge

AIMS

This study aimed to investigate the effect of palm oil mill effluent (POME) final discharge on the active bacterial composition, gene expression, and metabolite profiles in the receiving rivers to establish a foundation for identifying potential biomarkers for monitoring POME pollution in rivers.

METHODS AND RESULTS

The POME final discharge, upstream (unpolluted by POME), and downstream (effluent receiving point) parts of the rivers from two sites were physicochemically characterized. The taxonomic and gene profiles were then evaluated using de novo metatranscriptomics, while the metabolites were detected using qualitative metabolomics. A similar bacterial community structure in the POME final discharge samples from both sites was recorded, but their composition varied. Redundancy analysis showed that several families, particularly Comamonadaceae and Burkholderiaceae [Pr(>F) = 0.028], were positively correlated with biochemical oxygen demand (BOD5) and chemical oxygen demand (COD). The results also showed significant enrichment of genes regulating various metabolisms in the POME-receiving rivers, with methane, carbon fixation pathway, and amino acids among the predominant metabolisms identified (FDR < 0.05, PostFC > 4, and PPDE > 0.95). This was further validated through qualitative metabolomics, whereby amino acids were detected as the predominant metabolites.

CONCLUSIONS

The results suggest that genes regulating amino acid metabolism have significant potential for developing effective biomonitoring and bioremediation strategies in river water influenced by POME final discharge, fostering a sustainable palm oil industry.

Revisiting evolutionary trajectories and the organization of the Pleolipoviridae family

Archaeal pleomorphic viruses belonging to the Pleolipoviridae family represent an enigmatic group as they exhibit unique genomic features and are thought to have evolved through recombination with different archaeal plasmids. However, most of our understanding of the diversity and evolutionary trajectories of this clade comes from a handful of isolated representatives. Here we present 164 new genomes of pleolipoviruses obtained from metagenomic data of Australian hypersaline lakes and publicly available metagenomic data. We perform a comprehensive analysis on the diversity and evolutionary relationships of the newly discovered viruses and previously described pleolipoviruses. We propose to classify the viruses into five genera within the Pleolipoviridae family, with one new genus represented only by virus genomes retrieved in this study. Our data support the current hypothesis that pleolipoviruses reshaped their genomes through recombining with multiple different groups of plasmids, which is reflected in the diversity of their predicted replication strategies. We show that the proposed genus Epsilonpleolipovirus has evolutionary ties to pRN1-like plasmids from Sulfolobus, suggesting that this group could be infecting other archaeal phyla. Interestingly, we observed that the genome size of pleolipoviruses is correlated to the presence or absence of an integrase. Analyses of the host range revealed that all but one virus exhibit an extremely narrow range, and we show that the predicted tertiary structure of the spike protein is strongly associated with the host family, suggesting a specific adaptation to the host S-layer glycoprotein organization.

Genome-wide identification, phylogenetic relationships and expression patterns of the NOD-like receptor (NLR) gene family in flounder (Paralichthys olivaceus)

NOD-like receptors (NLRs) are one of the pattern recognition receptors which have been widely known for identifying pathogens and regulating innate immunity in mammals, but the functions of the NLR gene family in teleost fish remain poorly understood. In this study, we conducted a comprehensive identification and analysis of the flounder (Paralichthys olivaceus) NLR gene family, including bioinformatics information, evolutionary relationships, gene structures, conserved motifs, domain composition, expression patterns and protein-protein interaction (PPI). We identified 22 NLRs in flounder (flNLRs) which were clustered into three subfamilies according to their domain organizations and phylogenetic features, i.e., NLR-A (6 members) resembling mammalian NODs, NLR-B (1 member) resembling mammalian NLRPs, and NLR-C (15 members) unique to teleost fish. All flNLRs shared a conserved NACHT domain including an N-terminal nucleotide-binding domain, a middle helical domain 1, and a winged helix domain. Gene structure analysis displayed that flNLRs were significantly different, with exon numbers from 1 to 52. Conserved domain analysis showed that the N-terminus of flNLRs possessed different characteristics of the domains including CARD domain, PYRIN domain, RING domain, and fish-specific FISNA domain, and the C-terminus of seven NLR-C members contained an extra B30.2 domain, named NLRC-B30.2 group. Notably, flNLRs were expressed in all nine tested tissues, showing higher expressions in the systemic and mucosal immune tissues (e.g., kidney, spleen, hindgut, gills, skin, liver) in healthy flounder, and significant responses to intraperitoneal injection and immersion immunization of inactivated Vibrio anguillarum in mucosal tissues, especially the NLR-C members. In addition, PPI analysis demonstrated that some flNLRs of NLR-A and NLR-C shared the same interacting proteins such as RIPK2, TRAF6, MAVS, CASP, ASC, and ATG5, suggesting they might play crucial roles in host defense, antiviral innate immunity, inflammation, apoptosis and autophagy. This study for the first time characterized the NLR gene family of flounder at the genome-wide level, and the results provided a better understanding of the evolution of the NLR gene family and their immune functions in innate immunity in fish.

Substrate Specificity of Biofilms Proximate to Historic Shipwrecks

The number of built structures on the seabed, such as shipwrecks, energy platforms, and pipelines, is increasing in coastal and offshore regions. These structures, typically composed of steel or wood, are substrates for microbial attachment and biofilm formation. The success of biofilm growth depends on substrate characteristics and local environmental conditions, though it is unclear which feature is dominant in shaping biofilm microbiomes. The goal of this study was to understand the substrate- and site-specific impacts of built structures on short-term biofilm composition and functional potential. Seafloor experiments were conducted wherein steel and wood surfaces were deployed for four months at distances extending up to 115 m away from three historic (>50 years old) shipwrecks in the Gulf of Mexico. DNA from biofilms on the steel and wood was extracted, and metagenomes were sequenced on an Illumina NextSeq. A bioinformatics analysis revealed that the taxonomic composition was significantly different between substrates and sites, with substrate being the primary determining factor. Regardless of site, the steel biofilms had a higher abundance of genes related to biofilm formation, and sulfur, iron, and nitrogen cycling, while the wood biofilms showed a higher abundance of manganese cycling and methanol oxidation genes. This study demonstrates how substrate composition shapes biofilm microbiomes and suggests that marine biofilms may contribute to nutrient cycling at depth. Analyzing the marine biofilm microbiome provides insight into the ecological impact of anthropogenic structures on the seabed.

Mechanism of herpesvirus UL24 protein regulating viral immune escape and virulence

Herpesviruses have evolved a series of abilities involved in the process of host infection that are conducive to virus survival and adaptation to the host, such as immune escape, latent infection, and induction of programmed cell death for sustainable infection. The herpesvirus gene UL24 encodes a highly conserved core protein that plays an important role in effective viral infection. The UL24 protein can inhibit the innate immune response of the host by acting on multiple immune signaling pathways during virus infection, and it also plays a key role in the proliferation and pathogenicity of the virus in the later stage of infection. This article reviews the mechanism by which the UL24 protein mediates herpesvirus immune escape and its effects on viral proliferation and virulence by influencing syncytial formation, DNA damage and the cell cycle. Reviewing these studies will enhance our understanding of the pathogenesis of herpesvirus infection and provide evidence for new strategies to combat against viral infection.

Lysine lactylation regulates metabolic pathways and biofilm formation in Streptococcus mutans

In eukaryotes, lactate produced during glycolysis is involved in regulating multiple metabolic processes through lysine lactylation (Kla). To explore the potential link between metabolism and Kla in prokaryotes, we investigated the distribution of Kla in the cariogenic bacterium Streptococcus mutans during planktonic growth in low-sugar conditions and in biofilm-promoting, high-sugar conditions. We identified 1869 Kla sites in 469 proteins under these two conditions, with the biofilm growth state showing a greater number of lactylated sites and proteins. Although high sugar increased Kla globally, it reduced lactylation of RNA polymerase subunit α (RpoA) at Lys173. Lactylation at this residue inhibited the synthesis of extracellular polysaccharides, a major constituent of the cariogenic biofilm. The Gcn5-related N-acetyltransferase (GNAT) superfamily enzyme GNAT13 exhibited lysine lactyltransferase activity in cells and lactylated Lys173 in RpoA in vitro. Either GNAT13 overexpression or lactylation of Lys173 in RpoA inhibited biofilm formation. These results provide an overview of the distribution and potential functions of Kla and improve our understanding of the role of lactate in the metabolic regulation of prokaryotes.

Feeding during hibernation shifts gene expression toward active season levels in brown bears (Ursus arctos)

Hibernation in bears involves a suite of metabolical and physiological changes, including the onset of insulin resistance, that are driven in part by sweeping changes in gene expression in multiple tissues. Feeding bears glucose during hibernation partially restores active season physiological phenotypes, including partial resensitization to insulin, but the molecular mechanisms underlying this transition remain poorly understood. Here, we analyze tissue-level gene expression in adipose, liver, and muscle to identify genes that respond to midhibernation glucose feeding and thus potentially drive postfeeding metabolical and physiological shifts. We show that midhibernation feeding stimulates differential expression in all analyzed tissues of hibernating bears and that a subset of these genes responds specifically by shifting expression toward levels typical of the active season. Inferences of upstream regulatory molecules potentially driving these postfeeding responses implicate peroxisome proliferator-activated receptor gamma (PPARG) and other known regulators of insulin sensitivity, providing new insight into high-level regulatory mechanisms involved in shifting metabolic phenotypes between hibernation and active states.

Genome of Phyllanthus emblica: the medicinal plant Amla with super antioxidant properties

Phyllanthus emblica or Indian gooseberry, commonly known as amla, is an important medicinal horticultural plant used in traditional and modern medicines. It bears stone fruits with immense antioxidant properties due to being one of the richest natural sources of vitamin C and numerous flavonoids. This study presents the first genome sequencing of this species performed using 10x Genomics and Oxford Nanopore Technology. The draft genome assembly was 519 Mbp in size and consisted of 4,384 contigs, N50 of 597 Kbp, 98.4% BUSCO score, and 37,858 coding sequences. This study also reports the genome-wide phylogeny of this species with 26 other plant species that resolved the phylogenetic position of P. emblica. The presence of three ascorbate biosynthesis pathways including L-galactose, galacturonate, and myo-inositol pathways was confirmed in this genome. A comprehensive comparative evolutionary genomic analysis including gene family expansion/contraction and identification of multiple signatures of adaptive evolution provided evolutionary insights into ascorbate and flavonoid biosynthesis pathways and stone fruit formation through lignin biosynthesis. The availability of this genome will be beneficial for its horticultural, medicinal, dietary, and cosmetic applications and will also help in comparative genomics analysis studies.

Epigenetic mechanisms of lncRNA in response to thermal stress during embryogenesis of allotetraploid Cyprinus carpio

Embryogenesis and epigenetic mechanisms of lncRNA may play an important role in the formation of temperature tolerance in allotetraploid Cyprinus carpio. To investigate the response of lncRNA to thermal stress during embryogenesis of C. carpio, transcriptome sequencing was performed on 81 embryo or larva samples from different early development stages and temperatures. We identified 45,097 lncRNAs and analyzed transcriptome variation during embryogenesis. Stage-specific and temperature-specific DE lncRNAs and DEGs were screened. GO and KEGG analysis identified numerous pathways involved in thermal stress. Temperature-specific regulation of cis-/trans-/antisense lncRNAs was analyzed. Interaction network analysis identified 6 hub lncRNAs and many hub genes, such as cdk1 and hsf1. Decreased expression of many essential genes regulated by lncRNAs may lead to the death of embryos at 33 °C. Our findings provide new insights into the regulation of lncRNA in thermal stress response during embryogenesis and contribute to the understanding of environmental adaptation of allotetraploid species.

Characterization of a novel polyextremotolerant fungus, Exophiala viscosa, with insights into its melanin regulation and ecological niche

Black yeasts are polyextremotolerant fungi that contain high amounts of melanin in their cell wall and maintain a primar yeast form. These fungi grow in xeric, nutrient depletes environments which implies that they require highly flexible metabolisms and have been suggested to contain the ability to form lichen-like mutualisms with nearby algae and bacteria. However, the exact ecological niche and interactions between these fungi and their surrounding community are not well understood. We have isolated 2 novel black yeasts from the genus Exophiala that were recovered from dryland biological soil crusts. Despite notable differences in colony and cellular morphology, both fungi appear to be members of the same species, which has been named Exophiala viscosa (i.e. E. viscosa JF 03-3 Goopy and E. viscosa JF 03-4F Slimy). A combination of whole genome sequencing, phenotypic experiments, and melanin regulation experiments have been performed on these isolates to fully characterize these fungi and help decipher their fundamental niche within the biological soil crust consortium. Our results reveal that E. viscosa is capable of utilizing a wide variety of carbon and nitrogen sources potentially derived from symbiotic microbes, can withstand many forms of abiotic stresses, and excretes melanin which can potentially provide ultraviolet resistance to the biological soil crust community. Besides the identification of a novel species within the genus Exophiala, our study also provides new insight into the regulation of melanin production in polyextremotolerant fungi.

Characterizing interactions in E-cadherin assemblages

Cadherin intermolecular interactions are critical for cell-cell adhesion and play essential roles in tissue formation and the maintenance of tissue structures. In this study, we focus on E-cadherin, a classical cadherin that connects epithelial cells, to understand how they interact in cis and trans conformations when attached to the same cell or opposing cells. We employ coevolutionary sequence analysis and molecular dynamics simulations to confirm previously known interaction sites as well as to identify new interaction sites. The sequence coevolutionary results yield a surprising result indicating that there are no strongly favored intermolecular interaction sites, which is unusual and suggests that many interaction sites may be possible, with none being strongly preferred over others. By using molecular dynamics, we test the persistence of these interactions and how they facilitate adhesion. We build several types of cadherin assemblages, with different numbers and combinations of cis and trans interfaces to understand how these conformations act to facilitate adhesion. Our results suggest that, in addition to the established interaction sites on the EC1 and EC2 domains, an additional plausible cis interface at the EC3-EC5 domain exists. Furthermore, we identify specific mutations at cis/trans binding sites that impair adhesion within E-cadherin assemblages.

Chromosome-scale genomes of commercial timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis)

Wood is the most important natural and endlessly renewable source of energy. Despite the ecological and economic importance of wood, many aspects of its formation have not yet been investigated. We performed chromosome-scale genome assemblies of three timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis) which exhibit different wood properties such as wood density, hardness, growth rate, and fiber cell wall thickness. The combination of 10X, stLFR, Hi-Fi sequencing and HiC data led us to assemble high-quality genomes evident by scaffold N50 length of 55.97 Mb (O. pyramidale), 22.37 Mb (M. ferrea) and 14.55 Mb (T. grandis) with >97% BUSCO completeness of the assemblies. A total of 35774, 24027, and 44813 protein-coding genes were identified in M. ferrea, T. grandis and O. pyramidale, respectively. The data generated in this study is anticipated to serve as a valuable genetic resource and will promote comparative genomic analyses, and it is of practical importance in gaining a further understanding of the wood properties in non-model woody species.

High-quality genome assemblies provide clues on the evolutionary advantage of blue peafowl over green peafowl

An intriguing example of differential adaptability is the case of two Asian peafowl species, Pavo cristatus (blue peafowl) and Pavo muticus (green peafowl), where the former has a "Least Concern" conservation status and the latter is an "Endangered" species. To understand the genetic basis of this differential adaptability of the two peafowl species, a comparative analysis of these species is much needed to gain the genomic and evolutionary insights. Thus, we constructed a high-quality genome assembly of blue peafowl with an N50 value of 84.81 Mb (pseudochromosome-level assembly), and a high-confidence coding gene set to perform the genomic and evolutionary analyses of blue and green peafowls with 49 other avian species. The analyses revealed adaptive evolution of genes related to neuronal development, immunity, and skeletal muscle development in these peafowl species. Major genes related to axon guidance such as NEO1 and UNC5, semaphorin (SEMA), and ephrin receptor showed adaptive evolution in peafowl species. However, blue peafowl showed the presence of 42% more coding genes compared to the green peafowl along with a higher number of species-specific gene clusters, segmental duplicated genes and expanded gene families, and comparatively higher evolution in neuronal and developmental pathways. Blue peafowl also showed longer branch length compared to green peafowl in the species phylogenetic tree. These genomic insights obtained from the high-quality genome assembly of P. cristatus constructed in this study provide new clues on the superior adaptability of the blue peafowl over green peafowl despite having a recent species divergence time.

Linked candidate genes of different functions for white mold resistance in common bean (Phaseolus vulgaris L) are identified by multiple QTL mapping approaches

White mold (WM) is a major disease in common bean (Phaseolus vulgaris L.), and its complex quantitative genetic control limits the development of WM resistant cultivars. WM2.2, one of the nine meta-QTL with a major effect on WM tolerance, explains up to 35% of the phenotypic variation and was previously mapped to a large genomic interval on Pv02. Our objective was to narrow the interval of this QTL using combined approach of classic QTL mapping and QTL-based bulk segregant analysis (BSA), and confirming those results with Khufu de novo QTL-seq. The phenotypic and genotypic data from two RIL populations, 'Raven'/I9365-31 (R31) and 'AN-37'/PS02-029C-20 (Z0726-9), were used to select resistant and susceptible lines to generate subpopulations for bulk DNA sequencing. The QTL physical interval was determined by considering overlapping interval of the identified QTL or peak region in both populations by three independent QTL mapping analyses. Our findings revealed that meta-QTL WM2.2 consists of three regions, WM2.2a (4.27-5.76 Mb; euchromatic), WM 2.2b (12.19 to 17.61 Mb; heterochromatic), and WM2.2c (23.01-25.74 Mb; heterochromatic) found in both populations. Gene models encoding for gibberellin 2-oxidase 8, pentatricopeptide repeat, and heat-shock proteins are the likely candidate genes associated with WM2.2a resistance. A TIR-NBS-LRR class of disease resistance protein (Phvul.002G09200) and LRR domain containing family proteins are potential candidate genes associated with WM2.2b resistance. Nine gene models encoding disease resistance protein [pathogenesis-related thaumatin superfamily protein and disease resistance-responsive (dirigent-like protein) family protein etc] found within the WM2.2c QTL interval are putative candidate genes. WM2.2a region is most likely associated with avoidance mechanisms while WM2.2b and WM2.2c regions trigger physiological resistance based on putative candidate genes.

bHLH transcription factor family identification, phylogeny, and its response to abiotic stress in Chenopodium quinoa

The second-largest transcription factor superfamily in plants is that of the basic helix-loop-helix (bHLH) family, which plays an important complex physiological role in plant growth, tissue development, and environmental adaptation. Systematic research on the Chenopodium quinoa bHLH family will enable a better understanding of this species. Herein, authors used a variety of bioinformatics methods and quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) to explore the evolution and function of the 218 CqbHLH genes identified. A total of 218 CqbHLH transcription factor genes were identified in the whole genome, located on 18 chromosomes. A phylogenetic tree was constructed using the CqbHLH and AtbHLH proteins to determine their homology, and the members were divided into 20 subgroups and one unclustered gene. Authors also analyzed 218 CqbHLH genes, conservative motifs, chromosome diffusion, and gene replication. The author constructed one Neighbor-Joining (NJ) tree and a collinearity analysis map of the bHLH family in C. quinoa and six other plant species to study the evolutionary relationship and homology among multiple species. In addition, the expression levels of 20 CqbHLH members from different subgroups in various tissues, different fruit developmental stages, and six abiotic stresses were analyzed. Authors identified 218 CqbHLH genes and studied their biological functions, providing a basis for better understanding and further studying the bHLH family in quinoa.

Identification and Evolutionary Analysis of Cotton (Gossypium hirsutum) WOX Family Genes and Their Potential Function in Somatic Embryogenesis

WUSCHEL-related homeobox (WOX) proteins participate profoundly in plant development and stress responses. As the difficulty of somatic embryogenesis severely constrains cotton genetic modification, in this study, we identified and comprehensively analyzed WOX genes in cotton. As a result, 40 WOX genes were identified in the upland cotton genome. All these cotton WOX genes were classified into three clades, ancient, intermediate, and modern clades, based on the phylogenetic analysis of previous studies. The majority (24) of the cotton WOX genes belonged to the modern clade, in which all gene members contain the vital functional domain WUS-box, which is necessary for plant stem cell regulation and maintenance. Collinearity analysis indicated that the WOX gene family in cotton expanded to some degree compared to Arabidopsis, especially in the modern clade. Genome duplication and segmental duplication may greatly contribute to expansion. Hormone-response- and abiotic-stress-response-related cis-acting regulatory elements were widely distributed in the promoter regions of cotton WOX genes, suggesting that the corresponding functions of stress responses and the participation of development processes were involved in hormone responses. By RNA sequencing, we profiled the expression patterns of cotton WOX genes in somatic embryogenesis. Only about half of cotton WOX genes were actively expressed during somatic embryogenesis; different cotton WOX genes may function in different development stages. The most representative, GhWOX4 and GhWOX13, may function in almost all stages of somatic embryogenesis; GhWOX2 and GhWOX9 function in the late stages of embryo patterning and embryo development during cotton somatic embryogenesis. Co-expression analysis showed that the cotton WOXs co-expressed with genes involved in extensive genetic information processing, including DNA replication, DNA repair, homologous recombination, RNA transport, protein processing, and several signaling and metabolism pathways, in which plant hormones signal transduction, MAPK signaling pathways, phosphatidylinositol signaling systems, and ABC transporters, as well as the metabolism of fatty acid; valine, leucine, and isoleucine biosynthesis; and cutin, suberine, and wax biosynthesis, were most significantly enriched. Taken together, the present study provides useful information and new insights into the functions of cotton WOX genes during somatic embryogenesis. The specific regulatory roles of some WOX genes in somatic embryogenesis are worthy of further functional research.

Metapangenomic investigation provides insight into niche differentiation of methanogenic populations from the subsurface serpentinizing environment, Samail Ophiolite, Oman

Serpentinization reactions produce highly reduced waters that have hyperalkaline pH and that can have high concentrations of H₂ and CH₄. Putatively autotrophic methanogenic archaea have been identified in the subsurface waters of the Samail Ophiolite, Sultanate of Oman, though the strategies to overcome hyperalkaline pH and dissolved inorganic carbon limitation remain to be fully understood. Here, we recovered metagenome assembled genomes (MAGs) and applied a metapangenomic approach to three different Methanobacterium populations to assess habitat-specific functional gene distribution. A Type I population was identified in the fluids with neutral pH, while a Type II and "Mixed" population were identified in the most hyperalkaline fluids (pH 11.63). The core genome of all Methanobacterium populations highlighted potential DNA scavenging techniques to overcome phosphate or nitrogen limitation induced by environmental conditions. With particular emphasis on the Mixed and Type II population found in the most hyperalkaline fluids, the accessory genomes unique to each population reflected adaptation mechanisms suggesting lifestyles that minimize niche overlap. In addition to previously reported metabolic capability to utilize formate as an electron donor and generate intracellular CO₂, the Type II population possessed genes relevant to defense against antimicrobials and assimilating potential osmoprotectants to provide cellular stability. The accessory genome of the Mixed population was enriched in genes for multiple glycosyltransferases suggesting reduced energetic costs by adhering to mineral surfaces or to other microorganisms, and fostering a non-motile lifestyle. These results highlight the niche differentiation of distinct Methanobacterium populations to circumvent the challenges of serpentinization impacted fluids through coexistence strategies, supporting our ability to understand controls on methanogenic lifestyles and adaptations within the serpentinizing subsurface fluids of the Samail Ophiolite.

Production of recombinant lytic polysaccharide monooxygenases and evaluation effect of its addition into Aspergillus fumigatus var. niveus cocktail for sugarcane bagasse saccharification

Lignocellulosic biomass is a promising alternative for producing biofuels, despite its recalcitrant nature. There are microorganisms in nature capable of efficiently degrade biomass, such as the filamentous fungi. Among them, Aspergillus fumigatus var. niveus (AFUMN) has a wide variety of carbohydrate-active enzymes (CAZymes), especially hydrolases, but a low number of oxidative enzymes in its genome. To confirm the enzymatic profile of this fungus, this study analyzed the secretome of AFUMN cultured in sugarcane bagasse as the sole carbon source. As expected, the secretome showed a predominance of hydrolytic enzymes compared to oxidative activity. However, it is known that hydrolytic enzymes act in synergy with oxidative proteins to efficiently degrade cellulose polymer, such as the Lytic Polysaccharide Monooxygenases (LPMOs). Thus, three LPMOs from the fungus Thermothelomyces thermophilus (TtLPMO9D, TtLPMO9H, and TtLPMO9O) were selected, heterologous expressed in Aspergillus nidulans, purified, and used to supplement the AFUMN secretome to evaluate their effect on the saccharification of sugarcane bagasse. The saccharification assay was carried out using different concentrations of AFUMN secretome supplemented with recombinant T. thermophilus LPMOs, as well as ascorbic acid as reducing agent for oxidative enzymes. Through a statistic design created by Design-Expert software, we were able to analyze a possible cooperative effect between these components. The results indicated that, in general, the addition of TtLPMO9D and ascorbic acid did not favor the conversion process in this study, while TtLPMO9O had a highly significant cooperative effect in bagasse saccharification compared to the control using only AFUMN secretome.

Identification of BBX gene family and its function in the regulation of microtuber formation in yam

BBX proteins play important roles in all of the major light-regulated developmental processes. However, no systematic analysis of BBX gene family regarding the regulation of photoperiodic microtuber formation has been previously performed in yam. In this study, a systematic analysis on the BBX gene family was conducted in three yam species, with the results, indicating that this gene plays a role in regulating photoperiodic microtuber formation. These analyses included identification the BBX gene family in three yam species, their evolutionary relationships, conserved domains, motifs, gene structure, cis-acting elements, and expressional patterns. Based on these analyses, DoBBX2/DoCOL5 and DoBBX8/DoCOL8 showing the most opposite pattern of expression during microtuber formation were selected as candidate genes for further investigation. Gene expression analysis showed DoBBX2/DoCOL5 and DoBBX8/DoCOL8 were highest expressed in leaves and exhibited photoperiod responsive expression patterns. Besides, the overexpression of DoBBX2/DoCOL5 and DoBBX8/DoCOL8 in potato accelerated tuber formation under short-day (SD) conditions, whereas only the overexpression of DoBBX8/DoCOL8 enhanced the accelerating effect of dark conditions on tuber induction. Tuber number was increased in DoBBX8/DoCOL8 overexpressing plants under dark, as well as in DoBBX2/DoCOL5 overexpressing plants under SD. Overall, the data generated in this study may form the basis of future functional characterizations of BBX genes in yam, especially regarding their regulation of microtuber formation via the photoperiodic response pathway.

Going through phages: a computational approach to revealing the role of prophage in Staphylococcus aureus

Prophages have important roles in virulence, antibiotic resistance, and genome evolution in Staphylococcus aureus . Rapid growth in the number of sequenced S. aureus genomes allows for an investigation of prophage sequences at an unprecedented scale. We developed a novel computational pipeline for phage discovery and annotation. We combined PhiSpy, a phage discovery tool, with VGAS and PROKKA, genome annotation tools to detect and analyse prophage sequences in nearly 10 011  S . aureus genomes, discovering thousands of putative prophage sequences with genes encoding virulence factors and antibiotic resistance. To our knowledge, this is the first large-scale application of PhiSpy on a large-scale set of genomes (10 011  S . aureus ). Determining the presence of virulence and resistance encoding genes in prophage has implications for the potential transfer of these genes/functions to other bacteria via transduction and thus can provide insight into the evolution and spread of these genes/functions between bacterial strains. While the phage we have identified may be known, these phages were not necessarily known or characterized in S. aureus and the clustering and comparison we did for phage based on their gene content is novel. Moreover, the reporting of these genes with the S. aureus genomes is novel.

Characterization and Coexpression Analysis of the TIFY Family Genes in Euryale ferox Related to Leaf Development

TIFYs are plant-specific transcription factors that contain the TIFY structural domain and play an important role in plant leaf growth and development. However, the role played by TIFY in E. ferox (Euryale ferox Salisb.) leaf development has not been investigated. In this study, 23 TIFY genes were identified in E. ferox. Phylogenetic analyses of the TIFY genes showed clustering into three groups (JAZ, ZIM, and PPD). The TIFY domain was shown to be conserved. JAZ was mainly expanded via wholegenome triplication (WGT) in E. ferox. Based on analyses of the TIFY genes in nine species, we found that JAZ has a closer relationship with PPD, in addition to appearing the most recently and expanding most rapidly, leading to the rapid expansion of TIFYs in Nymphaeaceae. In addition, their different evolution types were discovered. Different gene expressions showed the distinct and corresponsive expression patterns of the EfTIFYs in different stages of tissue and leaf development. Finally, The qPCR analysis revealed that the expression of EfTIFY7.2 and EfTIFY10.1 showed an upward trend and high expression throughout leaf development. Further co-expression analysis indicated that EfTIFY7.2 might be more important for the development of E. ferox leaves. This information will be valuable when exploring the molecular mechanisms of EfTIFYs in plants.

Genome-Wide Identification of B-Box Gene Family and Expression Analysis Suggest Its Roles in Responses to Cercospora Leaf Spot in Sugar Beet (Beta Vulgaris L.)

The B-box (BBX) protein, which is a zinc-finger protein containing one or two B-box domains, plays a crucial role in the growth and development of plants. Plant B-box genes are generally involved in morphogenesis, the growth of floral organs, and various life activities in response to stress. In this study, the sugar beet B-box genes (hereafter referred to as BvBBXs) were identified by searching the homologous sequences of the Arabidopsis thaliana B-box gene family. The gene structure, protein physicochemical properties, and phylogenetic analysis of these genes were systematically analyzed. In this study, 17 B-box gene family members were identified from the sugar beet genome. A B-box domain can be found in all sugar beet BBX proteins. BvBBXs encode 135 to 517 amino acids with a theoretical isoelectric point of 4.12 to 6.70. Chromosome localization studies revealed that BvBBXs were dispersed across nine sugar beet chromosomes except chromosomes 5 and 7. The sugar beet BBX gene family was divided into five subfamilies using phylogenetic analysis. The gene architectures of subfamily members on the same evolutionary tree branch are quite similar. Light, hormonal, and stress-related cis-acting elements can be found in the promoter region of BvBBXs. The BvBBX gene family was differently expressed in sugar beet following Cercospora leaf spot infection, according to RT-qPCR data. It is shown that the BvBBX gene family may influence how the plant reacts to a pathogen infection.

Worker Reproduction and Caste Polymorphism Impact Genome Evolution and Social Genes Across the Ants

Eusocial insects are characterized by several traits, including reproductive division of labor and caste polymorphisms, which likely modulate genome evolution. Concomitantly, evolution may act on specific genes and pathways underlying these novel, sociality-associated phenotypes. Reproductive division of labor should increase the magnitude of genetic drift and reduce the efficacy of selection by reducing effective population size. Caste polymorphism has been associated with relaxed selection and may facilitate directional selection on caste-specific genes. Here, we use comparative analyses of 22 ant genomes to test how reproductive division of labor and worker polymorphism influence positive selection and selection intensity across the genome. Our results demonstrate that worker reproductive capacity is associated with a reduction in the degree of relaxed selection but is not associated with any significant change to positive selection. We find decreases in positive selection in species with polymorphic workers, but no increase in the degree of relaxed selection. Finally, we explore evolutionary patterns in specific candidate genes associated with our focal traits in eusocial insects. Two oocyte patterning genes previously implicated in worker sterility evolve under intensified selection in species with reproductive workers. Behavioral caste genes generally experience relaxed selection associated with worker polymorphism, whereas vestigial and spalt, both associated with soldier development in Pheidole ants, experience intensified selection in worker polymorphic species. These findings expand our understanding of the genetic mechanisms underlying elaborations of sociality. The impacts of reproductive division of labor and caste polymorphisms on specific genes illuminate those genes' roles in generating complex eusocial phenotypes.

Identification of a covert evolutionary pathway between two protein folds

Although homologous protein sequences are expected to adopt similar structures, some amino acid substitutions can interconvert α-helices and β-sheets. Such fold switching may have occurred over evolutionary history, but supporting evidence has been limited by the: (1) abundance and diversity of sequenced genes, (2) quantity of experimentally determined protein structures, and (3) assumptions underlying the statistical methods used to infer homology. Here, we overcome these barriers by applying multiple statistical methods to a family of ~600,000 bacterial response regulator proteins. We find that their homologous DNA-binding subunits assume divergent structures: helix-turn-helix versus α-helix + β-sheet (winged helix). Phylogenetic analyses, ancestral sequence reconstruction, and AlphaFold2 models indicate that amino acid substitutions facilitated a switch from helix-turn-helix into winged helix. This structural transformation likely expanded DNA-binding specificity. Our approach uncovers an evolutionary pathway between two protein folds and provides a methodology to identify secondary structure switching in other protein families.

Genome-wide identification and expression profiling of two-component system (TCS) genes in Brassica oleracea in response to shade stress

The Two-component system (TCS) consists of Histidine kinases (HKs), Phosphotransfers (HPs), and response regulator (RR) proteins. It has an important role in signal transduction to respond to a wide variety of abiotic stresses and hence in plant development. Brassica oleracea (cabbage) is a leafy vegetable, which is used for food and medicinal purposes. Although this system was identified in several plants, it had not been identified in Brassica oleracea yet. This genome-wide study identified 80 BoTCS genes consisting of 21 HKs, 8 HPs, 39 RRs, and 12 PRRs. This classification was done based on conserved domains and motif structure. Phylogenetic relationships of BoTCS genes with Arabidopsis thaliana, Oryza sativa, Glycine max, and Cicer arietinum showed conservation in TCS genes. Gene structure analysis revealed that each subfamily had conserved introns and exons. Both tandem and segmental duplication led to the expansion of this gene family. Almost all of the HPs and RRs were expanded through segmental duplication. Chromosomal analysis showed that BoTCS genes were dispersed across all nine chromosomes. The promoter regions of these genes were found to contain a variety of cis-regulatory elements. The 3D structure prediction of proteins also confirmed the conservation of structure within subfamilies. MicroRNAs (miRNAs) involved in the regulation of BoTCSs were also predicted and their regulatory roles were also evaluated. Moreover, BoTCSs were docked with abscisic acid to evaluate their binding. RNA-seq-based expression analysis and validation by qRT-PCR showed significant variation of expression for BoPHYs, BoERS1.1, BoERS2.1, BoERS2.2, BoRR10.2, and BoRR7.1 suggesting their importance in stress response. These genes showing unique expression can be further used in manipulating the plant's genome to make the plant more resistant the environmental stresses which will ultimately help in the increase of plant's yield. More specifically, these genes have altered expression in shade stress which clearly indicates their importance in biological functions. These findings are important for future functional characterization of TCS genes in generating stress-responsive cultivars.

Molecular analysis of the 14-3-3 genes in Panax ginseng and their responses to heat stress

Background

Panax Ginseng is a perennial and semi-shady herb with tremendous medicinal value. Due to its unique botanical characteristics, ginseng is vulnerable to various abiotic factors during its growth and development, especially in high temperatures. Proteins encoded by 14-3-3 genes form a highly conserved protein family that widely exists in eukaryotes. The 14-3-3 family regulates the vital movement of cells and plays an essential role in the response of plants to abiotic stresses, including high temperatures. Currently, there is no relevant research on the 14-3-3 genes of ginseng.

Methods

The identification of the ginseng 14-3-3 gene family was mainly based on ginseng genomic data and Hidden Markov Models (HMM). We used bioinformatics-related databases and tools to analyze the gene structure, physicochemical properties, cis-acting elements, gene ontology (GO), phylogenetic tree, interacting proteins, and transcription factor regulatory networks. We analyzed the transcriptome data of different ginseng tissues to clarify the expression pattern of the 14-3-3 gene family in ginseng. The expression level and modes of 14-3-3 genes under heat stress were analyzed by quantitative real-time PCR (qRT-PCR) technology to determine the genes in the 14-3-3 gene family responding to high-temperature stress.

Results

In this study, 42 14-3-3 genes were identified from the ginseng genome and renamed PgGF14-1 to PgGF14-42. Gene structure and evolutionary relationship research divided PgGF14s into epsilon (ε) and non-epsilon (non-ε) groups, mainly located in four evolutionary branches. The gene structure and motif remained highly consistent within a subgroup. The physicochemical properties and structure of the predicted PgGF14 proteins conformed to the essential characteristics of 14-3-3 proteins. RNA-seq results indicated that the detected PgGF14s existed in different organs and tissues but differed in abundance; their expression was higher in roots, stems, leaves, and fruits but lower in seeds. The analysis of GO, cis-acting elements, interacting proteins, and regulatory networks of transcription factors indicated that PgGF14s might participate in physiological processes, such as response to stress, signal transduction, material synthesis-metabolism, and cell development. The qRT-PCR results indicated PgGF14s had multiple expression patterns under high-temperature stress with different change trends in several treatment times, and 38 of them had an apparent response to high-temperature stress. Furthermore, PgGF14-5 was significantly upregulated, and PgGF14-4 was significantly downregulated in all treatment times. This research lays a foundation for further study on the function of 14-3-3 genes and provides theoretical guidance for investigating abiotic stresses in ginseng.

Identification of Aspergillus niger Aquaporins Involved in Hydrogen Peroxide Signaling

Aspergillus niger is a robust microbial cell factory for organic acid production. However, the regulation of many industrially important pathways is still poorly understood. The regulation of the glucose oxidase (Gox) expression system, involved in the biosynthesis of gluconic acid, has recently been uncovered. The results of that study show hydrogen peroxide, a by-product of the extracellular conversion of glucose to gluconate, has a pivotal role as a signaling molecule in the induction of this system. In this study, the facilitated diffusion of hydrogen peroxide via aquaporin water channels (AQPs) was studied. AQPs are transmembrane proteins of the major intrinsic proteins (MIPs) superfamily. In addition to water and glycerol, they may also transport small solutes such as hydrogen peroxide. The genome sequence of A. niger N402 was screened for putative AQPs. Seven AQPs were found and could be classified into three main groups. One protein (AQPA) belonged to orthodox AQP, three (AQPB, AQPD, and AQPE) were grouped in aquaglyceroporins (AQGP), two (AQPC and AQPF) were in X-intrinsic proteins (XIPs), and the other (AQPG) could not be classified. Their ability to facilitate diffusion of hydrogen peroxide was identified using yeast phenotypic growth assays and by studying AQP gene knock-outs in A. niger. The X-intrinsic protein AQPF appears to play roles in facilitating hydrogen peroxide transport across the cellular membrane in both Saccharomyces cerevisiae and A. niger experiments.

Viruses of sulfur oxidizing phototrophs encode genes for pigment, carbon, and sulfur metabolisms

Viral infections modulate bacterial metabolism and ecology. Here, we investigated the hypothesis that viruses influence the ecology of purple and green sulfur bacteria in anoxic and sulfidic lakes, analogs of euxinic oceans in the geologic past. By screening metagenomes from lake sediments and water column, in addition to publicly-available genomes of cultured purple and green sulfur bacteria, we identified almost 300 high and medium-quality viral genomes. Viruses carrying the gene psbA, encoding the small subunit of photosystem II protein D1, were ubiquitous, suggesting viral interference with the light reactions of sulfur oxidizing autotrophs. Viruses predicted to infect these autotrophs also encoded auxiliary metabolic genes for reductive sulfur assimilation as cysteine, pigment production, and carbon fixation. These observations show that viruses have the genomic potential to modulate the production of metabolic markers of phototrophic sulfur bacteria that are used to identify photic zone euxinia in the geologic past.

Insights into the mechanisms of triptolide nephrotoxicity through network pharmacology-based analysis and RNA-seq

INTRODUCTION

Triptolide (TPL) is a promising plant-derived compound for clinical therapy of multiple human diseases; however, its application was limited considering its toxicity.

METHODS

To explore the underlying molecular mechanism of TPL nephrotoxicity, a network pharmacology based approach was utilized to predict candidate targets related with TPL toxicity, followed by deep RNA-seq analysis to characterize the features of three transcriptional elements include protein coding genes (PCGs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) as well as their associations with nephrotoxicity in rats with TPL treatment.

RESULTS & DISCUSSION

Although the deeper mechanisms of TPL nephrotoxcity remain further exploration, our results suggested that c-Jun is a potential target of TPL and Per1 related circadian rhythm signaling is involved in TPL induced renal toxicity.

Rapid evolution of T2/S-RNase genes in Fragaria linked to multiple transitions from self-incompatibility to self-compatibility

The T2/RNase gene family is widespread in eukaryotes, and particular members of this family play critical roles in the gametophytic self-incompatibility (GSI) system in plants. Wild diploid strawberry (Fragaria) species have diversified their sexual systems via self-incompatible and self-compatible traits, yet how these traits evolved in Fragaria remains elusive. By integrating the published and de novo assembled genomes and the newly generated RNA-seq data, members of the RNase T2 gene family were systematically identified in six Fragaria species, including three self-incompatible species (Fragaria nipponica, Fragaria nubicola, and Fragaria viridis) and three self-compatible species (Fragaria nilgerrensis, Fragaria vesca, and Fragaria iinumae). In total, 115 RNase T2 genes were identified in the six Fragaria genomes and can be classified into three classes (I-III) according to phylogenetic analysis. The identified RNase T2 genes could be divided into 22 homologous gene sets according to amino acid sequence similarity and phylogenetic and syntenic relationships. We found that extensive gene loss and pseudogenization coupled with small-scale duplications mainly accounted for variations in the RNase T2 gene numbers in Fragaria. Multiple copies of homologous genes were mainly generated from tandem and segmental duplication events. Furthermore, we newly identified five S-RNase genes in three self-incompatible Fragaria genomes, including two in F. nipponica, two in F. viridis, and one in F. nubicola, which fit for typical features of a pistil determinant, including highly pistil-specific expression, highly polymorphic proteins and alkaline isoelectric point (pI), while no S-RNase genes were found in all three self-compatible Fragaria species. Surprisingly, these T2/S-RNase genes contain at least one large intron (>10 kb). This study revealed that the rapid evolution of T2/S-RNase genes within the Fragaria genus could be associated with its sexual mode, and repeated evolution of the self-compatible traits in Fragaria was convergent via losses of S-RNase.

Distinguishing features of fold-switching proteins

Though many folded proteins assume one stable structure that performs one function, a small-but-increasing number remodel their secondary and tertiary structures and change their functions in response to cellular stimuli. These fold-switching proteins regulate biological processes and are associated with autoimmune dysfunction, severe acute respiratory syndrome coronavirus-2 infection, and more. Despite their biological importance, it is difficult to computationally predict fold switching. With the aim of advancing computational prediction and experimental characterization of fold switchers, this review discusses several features that distinguish fold-switching proteins from their single-fold and intrinsically disordered counterparts. First, the isolated structures of fold switchers are less stable and more heterogeneous than single folders but more stable and less heterogeneous than intrinsically disordered proteins (IDPs). Second, the sequences of single fold, fold switching, and intrinsically disordered proteins can evolve at distinct rates. Third, proteins from these three classes are best predicted using different computational techniques. Finally, late-breaking results suggest that single folders, fold switchers, and IDPs have distinct patterns of residue-residue coevolution. The review closes by discussing high-throughput and medium-throughput experimental approaches that might be used to identify new fold-switching proteins.