HMMER web server: 2015 update

Communication leads to bacterial heterogeneous adaptation to changing conditions in partial nitrification reactor

Recently, it is reported that bacterial communication coordinates the whole consortia to jointly resist the adverse environments. Here, we found the bacterial communication inevitably distinguished bacterial adaptation among different species in partial nitrification reactor under decreasing temperatures. We operated a partial nitrification reactor under temperature gradient from 30 °C to 5 °C and found the promotion of bacterial communication on adaptation of ammonia-oxidizing bacteria (AOB) was greater than that of nitrite-oxidizing bacteria (NOB). Signal pathways with single-component sensing protein in AOB can regulate more genes involved in bacterial adaptation than that with two-component sensing protein in NOB. The negative effects of bacterial communication, which were seriously ignored, have been highlighted, and Clp regulator downstream diffusible signal factor (DSF) based signal pathways worked as transcription activators and inhibitors of adaptation genes in AOB and NOB respectively. Bacterial communication can induce differential adaptation through influencing bacterial interactions. AOB inclined to cooperate with DSF synthesis bacteria as temperature declined, however, cooperation between NOB and DSF synthesis bacteria inclined to get weakening. According to the regulatory effects of signal pathways, bacterial survival strategies for self-protection were revealed. This study hints a potential way to govern niche differentiation in the microbiota by bacterial communication, contributing to forming an efficient artificial ecosystem.

The phage shock protein (PSP) envelope stress response: discovery of novel partners and evolutionary history

Bacterial phage shock protein (PSP) systems stabilize the bacterial cell membrane and protect against envelope stress. These systems have been associated with virulence, but despite their critical roles, PSP components are not well characterized outside proteobacteria. Using comparative genomics and protein sequence-structure-function analyses, we systematically identified and analyzed PSP homologs, phyletic patterns, domain architectures, and gene neighborhoods. This approach underscored the evolutionary significance of the system, revealing that its core protein PspA (Snf7 in ESCRT outside bacteria) was present in the last universal common ancestor and that this ancestral functionality has since diversified into multiple novel, distinct PSP systems across life. Several novel partners of the PSP system were identified: (i) the Toastrack domain, likely facilitating assembly of sub-membrane stress-sensing and signaling complexes, (ii) the newly defined HTH-associated α-helical signaling domain-PadR-like transcriptional regulator pair system, and (iii) multiple independent associations with ATPase, CesT/Tir-like chaperone, and Band-7 domains in proteins thought to mediate sub-membrane dynamics. Our work also uncovered links between the PSP components and other domains, such as novel variants of SHOCT-like domains, suggesting roles in assembling membrane-associated complexes of proteins with disparate biochemical functions. Results are available at our interactive web app, https://jravilab.org/psp.IMPORTANCEPhage shock proteins (PSP) are virulence-associated, cell membrane stress-protective systems. They have mostly been characterized in Proteobacteria and Firmicutes. We now show that a minimal PSP system was present in the last universal common ancestor that evolved and diversified into newly identified functional contexts. Recognizing the conservation and evolution of PSP systems across bacterial phyla contributes to our understanding of stress response mechanisms in prokaryotes. Moreover, the newly discovered PSP modularity will likely prompt new studies of lineage-specific cell envelope structures, lifestyles, and adaptation mechanisms. Finally, our results validate the use of domain architecture and genetic context for discovery in comparative genomics.

Selective haematological cancer eradication with preserved haematopoiesis

Haematopoietic stem cell (HSC) transplantation (HSCT) is the only curative treatment for a broad range of haematological malignancies, but the standard of care relies on untargeted chemotherapies and limited possibilities to treat malignant cells after HSCT without affecting the transplanted healthy cells1. Antigen-specific cell-depleting therapies hold the promise of much more targeted elimination of diseased cells, as witnessed in the past decade by the revolution of clinical practice for B cell malignancies2. However, target selection is complex and limited to antigens expressed on subsets of haematopoietic cells, resulting in a fragmented therapy landscape with high development costs2-5. Here we demonstrate that an antibody-drug conjugate (ADC) targeting the pan-haematopoietic marker CD45 enables the antigen-specific depletion of the entire haematopoietic system, including HSCs. Pairing this ADC with the transplantation of human HSCs engineered to be shielded from the CD45-targeting ADC enables the selective eradication of leukaemic cells with preserved haematopoiesis. The combination of CD45-targeting ADCs and engineered HSCs creates an almost universal strategy to replace a diseased haematopoietic system, irrespective of disease aetiology or originating cell type. We propose that this approach could have broad implications beyond haematological malignancies.

Genome-wide identification and expression analysis of Ubiquitin-specific protease gene family in maize (Zea mays L.)

BACKGROUND

Ubiquitin-specific proteases (UBPs) are a large family of deubiquitinating enzymes (DUBs). They are widespread in plants and are critical for plant growth, development, and response to external stresses. However, there are few studies on the functional characteristics of the UBP gene family in the important staple crop, maize (Zea mays L.).

RESULTS

In this study, we performed a bioinformatic analysis of the entire maize genome and identified 45 UBP genes. Phylogenetic analysis indicated that 45 ZmUBP genes can be divided into 15 subfamilies. Analysis of evolutionary patterns and divergence levels indicated that ZmUBP genes were present before the isolation of dicotyledons, were highly conserved and subjected to purifying selection during evolution. Most ZmUBP genes exhibited different expression levels in different tissues and developmental stages. Based on transcriptome data and promoter element analysis, we selected eight ZmUBP genes whose promoters contained a large number of plant hormones and stress response elements and were up-regulated under different abiotic stresses for RT-qPCR analysis, results showed that these genes responded to abiotic stresses and phytohormones to varying degrees, indicating that they play important roles in plant growth and stress response.

CONCLUSIONS

In this study, the structure, location and evolutionary relationship of maize UBP gene family members were analyzed for the first time, and the ZmUBP genes that may be involved in stress response and plant growth were identified by combining promoter element analysis, transcriptome data and RT-qPCR analysis. This study informs research on the involvement of maize deubiquitination in stress response.

Genomic evidence for the widespread presence of GH45 cellulases among soil invertebrates

Lignocellulose is a major component of vascular plant biomass. Its decomposition is crucial for the terrestrial carbon cycle. Microorganisms are considered primary decomposers, but evidence increases that some invertebrates may also decompose lignocellulose. We investigated the taxonomic distribution and evolutionary origins of GH45 hydrolases, important enzymes for the decomposition of cellulose and hemicellulose, in a collection of soil invertebrate genomes. We found that these genes are common in springtails and oribatid mites. Phylogenetic analysis revealed that cellulase genes were acquired early in the evolutionary history of these groups. Domain architectures and predicted 3D enzyme structures indicate that these cellulases are functional. Patterns of presence and absence of these genes across different lineages prompt further investigation into their evolutionary and ecological benefits. The ubiquity of cellulase genes suggests that soil invertebrates may play a role in lignocellulose decomposition, independently or in synergy with microorganisms. Understanding the ecological and evolutionary implications might be crucial for understanding soil food webs and the carbon cycle.

Molecular characterization of the complete genome of a novel ormycovirus infecting the ectomycorrhizal fungus Hortiboletus rubellus

Advancements in high-throughput sequencing and the development of new bioinformatics tools for large-scale data analysis play a crucial role in uncovering virus diversity and enhancing our understanding of virus evolution. The discovery of the ormycovirus clades, a group of RNA viruses that are phylogenetically distinct from all known Riboviria members and are found in fungi, highlights the value of these tools for the discovery of novel viruses. The aim of this study was to examine viral populations in fungal hosts to gain insights into the diversity, evolution, and classification of these viruses. Here, we report the molecular characterization of a newly discovered ormycovirus, which we have named "Hortiboletus rubellus ormycovirus 1" (HrOMV1), that was found in the ectomycorrhizal fungus Hortiboletus rubellus. The bipartite genome of HrOMV1, whose nucleotide sequence was determined by HTS and RLM-RACE, consists of two RNA segments (RNA1 and RNA2) that exhibit similarity to those of previously studied ormycoviruses in their organization and the proteins they encode. The presence of upstream, in-frame AUG triplets in the 5' termini of both RNA segments suggests that HrOMV1, like certain other ormycoviruses, employs a non-canonical translation initiation strategy. Phylogenetic analysis showed that HrOMV1 is positioned within the gammaormycovirus clade. Its putative RNA-dependent RNA polymerase (RdRp) exhibits sequence similarity to those of other gammaormycovirus members, the most similarity to that of Termitomyces ormycovirus 1, with 33.05% sequence identity. This protein was found to contain conserved motifs that are crucial for RNA replication, including the distinctive GDQ catalytic triad observed in gammaormycovirus RdRps. The results of this study underscore the significance of investigating the ecological role of mycoviruses in mycorrhizal fungi. This is the first report of an ormycovirus infecting a member of the ectomycorrhizal genus Hortiboletus.

StCoExpNet: a global co-expression network analysis facilitates identifying genes underlying agronomic traits in potatoes

KEY MESSAGE

We constructed a gene expression atlas and co-expression network for potatoes and identified several novel genes associated with various agronomic traits. This resource will accelerate potato genetics and genomics research. Potato (Solanum tuberosum L.) is the world's most crucial non-cereal food crop and ranks third in food production after wheat and rice. Despite the availability of several potato transcriptome datasets at public databases like NCBI SRA, an effort has yet to be put into developing a global transcriptome atlas and a co-expression network for potatoes. The objectives of our study were to construct a global expression atlas for potatoes using publicly available transcriptome datasets, identify housekeeping and tissue-specific genes, construct a global co-expression network and identify co-expression clusters, investigate the transcriptional complexity of genes involved in various essential biological processes related to agronomic traits, and provide a web server (StCoExpNet) to easily access the newly constructed expression atlas and co-expression network to investigate the expression and co-expression of genes of interest. In this study, we used data from 2299 publicly available potato transcriptome samples obtained from 15 different tissues to construct a global transcriptome atlas. We found that roughly 87% of the annotated genes exhibited detectable expression in at least one sample. Among these, we identified 281 genes with consistent and stable expression levels, indicating their role as housekeeping genes. Conversely, 308 genes exhibited marked tissue-specific expression patterns. We exemplarily linked some co-expression clusters to important agronomic traits of potatoes, such as self-incompatibility, anthocyanin biosynthesis, tuberization, and defense responses against multiple pathogens. The dataset compiled here constitutes a new resource (StCoExpNet), which can be accessed at https://stcoexpnet.julius-kuehn.de . This transcriptome atlas and the co-expression network will accelerate potato genetics and genomics research.

Shaping human gut community assembly and butyrate production by controlling the arginine dihydrolase pathway

The arginine dihydrolase pathway (arc operon) present in a subset of diverse human gut species enables arginine catabolism. This specialized metabolic pathway can alter environmental pH and nitrogen availability, which in turn could shape gut microbiota inter-species interactions. By exploiting synthetic control of gene expression, we investigated the role of the arc operon in probiotic Escherichia coli Nissle 1917 on human gut community assembly and health-relevant metabolite profiles in vitro and in the murine gut. By stabilizing environmental pH, the arc operon reduced variability in community composition across different initial pH perturbations. The abundance of butyrate producing bacteria were altered in response to arc operon activity and butyrate production was enhanced in a physiologically relevant pH range. While the presence of the arc operon altered community dynamics, it did not impact production of short chain fatty acids. Dynamic computational modeling of pH-mediated interactions reveals the quantitative contribution of this mechanism to community assembly. In sum, our framework to quantify the contribution of molecular pathways and mechanism modalities on microbial community dynamics and functions could be applied more broadly.

Curation, nomenclature, and topological classification of receptor-like kinases from 528 plant species for novel domain discovery and functional inference

Receptor-like kinases (RLKs) are the most numerous signal transduction components in plants and play important roles in determining how different plants adapt to their ecological environments. Research on RLKs has focused mainly on a small number of typical RLK members in a few model plants. There is an urgent need to study the composition, distribution, and evolution of RLKs at the holistic level to increase our understanding of how RLKs assist in the ecological adaptations of different plant species. In this study, we collected the genome assemblies of 528 plant species and constructed an RLK dataset. Using this dataset, we identified and characterized 524 948 RLK family members. Each member underwent systematic topological classification and was assigned a gene ID based on a unified nomenclature system. Furthermore, we identified two novel extracellular domains in some RLKs, designated Xiao and Xiang. Evolutionary analysis of the RLK family revealed that the RLCK-XVII and RLCK-XII-2 classes were present exclusively in dicots, suggesting that diversification of RLKs between monocots and dicots may have led to differences in downstream cytoplasmic responses. We also used an interaction proteome to help empower data mining for inference of new RLK functions from a global perspective, with the ultimate goal of understanding how RLKs shape the adaptation of different plants to the environments/ecosystems. The assembled RLK dataset, together with annotations and analytical tools, forms an integrated foundation of multiomics data that is publicly accessible via the metaRLK web portal (http://metaRLK.biocloud.top).

Molecular and phylogenetic evidence of parallel expansion of anion channels in plants

Aluminum-activated malate transporters (ALMTs) and slow anion channels (SLACs) are important in various physiological processes in plants, including stomatal regulation, nutrient uptake, and in response to abiotic stress such as aluminum toxicity. To understand their evolutionary history and functional divergence, we conducted phylogenetic and expression analyses of ALMTs and SLACs in green plants. Our findings from phylogenetic studies indicate that ALMTs and SLACs may have originated from green algae and red algae, respectively. The ALMTs of early land plants and charophytes formed a monophyletic clade consisting of three subgroups. A single duplication event of ALMTs was identified in vascular plants and subsequent duplications into six clades occurred in angiosperms, including an identified clade, 1-1. The ALMTs experienced gene number losses in clades 1-1 and 2-1 and expansions in clades 1-2 and 2-2b. Interestingly, the expansion of clade 1-2 was also associated with higher expression levels compared to genes in clades that experienced apparent loss. SLACs first diversified in bryophytes, followed by duplication in vascular plants, giving rise to three distinct clades (I, II, and III), and clade II potentially associated with stomatal control in seed plants. SLACs show losses in clades II and III without substantial expansion in clade I. Additionally, ALMT clade 2-2 and SLAC clade III contain genes specifically expressed in reproductive organs and roots in angiosperms, lycophytes, and mosses, indicating neofunctionalization. In summary, our study demonstrates the evolutionary complexity of ALMTs and SLACs, highlighting their crucial role in the adaptation and diversification of vascular plants.

First characterization of the nuclear receptor superfamily in the Mediterranean mussel Mytilus galloprovincialis: developmental expression dynamics and potential susceptibility to environmental chemicals

Endocrine-disrupting chemicals (EDCs) represent a global threat to human health and the environment. In vertebrates, lipophilic EDCs primarily act by mimicking endogenous hormones, thus interfering with the transcriptional activity of nuclear receptors (NRs). The demonstration of the direct translation of these mechanisms into perturbation of NR-mediated physiological functions in invertebrates, however, has rarely proven successful, as the modes of action of EDCs in vertebrates and invertebrates seem to be distinct. In the present work, we investigated the members of the NR superfamily in a bivalve mollusk, the Mediterranean mussel Mytilus galloprovincialis. In addition to annotating the M. galloprovincialis NR complement, we assessed the potential developmental functions and susceptibility to EDC challenge during early development by gene expression analyses. Our results indicate that a majority of mussel NRs are dynamically expressed during early development, including receptors characterized by a potential susceptibility to EDCs. This study thus indicates that NRs are major regulators of early mussel development and that NR-mediated endocrine disruption in the mussel could be occurring at a larger scale and at earlier stages of the life cycle than previously anticipated. Altogether, these findings will have significant repercussions for our understanding of the stability of natural mussel populations. This article is part of the theme issue 'Endocrine responses to environmental variation: conceptual approaches and recent developments'.

Fungal microtubule organizing centers are evolutionarily unstable structures

For most Eukaryotic species the requirements of cilia formation dictate the structure of microtubule organizing centers (MTOCs). In this study we find that loss of cilia corresponds to loss of evolutionary stability for fungal MTOCs. We used iterative search algorithms to identify proteins homologous to those found in Saccharomyces cerevisiae, and Schizosaccharomyces pombe MTOCs, and calculated site-specific rates of change for those proteins that were broadly phylogenetically distributed. Our results indicate that both the protein composition of MTOCs as well as the sequence of MTOC proteins are poorly conserved throughout the fungal kingdom. To begin to reconcile this rapid evolutionary change with the rigid structure and essential function of the S. cerevisiae MTOC we further analyzed how structural interfaces among proteins influence the rates of change for specific residues within a protein. We find that a more stable protein may stabilize portions of an interacting partner where the two proteins are in contact. In summary, while the protein composition and sequences of the MTOC may be rapidly changing the proteins within the structure have a stabilizing effect on one another. Further exploration of fungal MTOCs will expand our understanding of how changes in the functional needs of a cell have affected physical structures, proteomes, and protein sequences throughout fungal evolution.

Comparative genomics of N-acetyl-5-methoxytryptamine members in four Prunus species with insights into bud dormancy and abiotic stress responses in Prunus avium

KEY MESSAGE

This study provides novel insights into the evolution, diversification, and functions of melatonin biosynthesis genes in Prunus species, highlighting their potential role in regulating bud dormancy and abiotic stresses. The biosynthesis of melatonin (MEL) in plants is primarily governed by enzymatic reactions involving key enzymes such as serotonin N-acetyltransferase (SNAT), tryptamine 5-hydroxylase (T5H), N-acetylserotonin methyltransferase (ASMT) and tryptophan decarboxylase (TDC). In this study, we analyzed Melatonin genes in four Prunus species such as Prunus avium (Pavi), Prunus pusilliflora (Ppus), Prunus serulata (Pser), and Prunus persica (Pper) based on comparative genomics approach. Among the four Prunus species, a total of 29 TDCs, 998 T5Hs, 16 SNATs, and 115 ASMTs within the genome of four Prunus genomes. A thorough investigation of melatonin-related genes was carried out using systematic biological methods and comparative genomics. Through phylogenetic analysis, orthologous clusters, Go enrichment, syntenic relationship, and gene duplication analysis, we discovered both similarities and variations in Melatonin genes among these Prunus species. Additionally, our study revealed the existence of unique subgroup members in the Melatonin genes of these species, which were distinct from those found in Arabidopsis genes. Furthermore, the transcriptomic expression analysis revealed the potential significance of melatonin genes in bud dormancy regulation and abiotic stresses. Our extensive results offer valuable perspectives on the evolutionary patterns, intricate expansion, and functions of PavMEL genes. Given their promising attributes, PavTDCs, PavT5H, PavNAT, and three PavASMT genes warrant in-depth exploration as prime candidates for manipulating dormancy in sweet cherry. This was done to lay the foundation for future explorations into the structural and functional aspects of these factors in Prunus species. This study offers significant insights into the functions of ASMT, SNAT, T5H, and TDC genes and sheds light on their roles in Prunus avium. Moreover, it established a robust foundation for further exploration functional characterization of melatonin genes in fruit species.

A chromosome level genome assembly of Pseudoroegneria Libanotica reveals a key Kcs gene involves in the cuticular wax elongation for drought resistance

BACKGROUND

The genus Pseudoroegneria (Nevski) Löve (Triticeae, Poaceae), whose genome symbol was designed as "St", accounts for more than 60% of perennial Triticeae species. The diploid species Psudoroegneria libanotica (2n = 14) contains the most ancient St genome, exhibited strong drought resistance, and was morphologically covered by cuticular wax on the aerial part. Therefore, the St-genome sequencing data could provide fundamental information for studies of genome evolution and reveal its mechanisms of cuticular wax and drought resistance.

RESULTS

In this study, we reported the chromosome-level genome assembly for the St genome of Pse. libanotica, with a total size of 2.99 Gb. 46,369 protein-coding genes annotated and 71.62% was repeat sequences. Comparative analyses revealed that the genus Pseudoroegneria diverged during the middle and late Miocene. During this period, unique genes, gene family expansion, and contraction in Pse. libanotica were enriched in biotic and abiotic stresses, such as fatty acid biosynthesis which may greatly contribute to its drought adaption. Furthermore, we investigated genes associated with the cuticular wax formation and water deficit and found a new Kcs gene evm.TU.CTG175.54. It plays a critical role in the very long chain fatty acid (VLCFA) elongation from C18 to C26 in Pse. libanotica. The function needs more evidence to be verified.

CONCLUSIONS

We sequenced and assembled the St genome in Triticeae and discovered a new KCS gene that plays a role in wax extension to cope with drought. Our study lays a foundation for the genome diversification of Triticeae species and deciphers cuticular wax formation genes involved in drought resistance.

The crystal structure of mycothiol disulfide reductase (Mtr) provides mechanistic insight into the specific low-molecular-weight thiol reductase activity of Actinobacteria

Low-molecular-weight (LMW) thiols are involved in many processes in all organisms, playing a protective role against reactive species, heavy metals, toxins and antibiotics. Actinobacteria, such as Mycobacterium tuberculosis, use the LMW thiol mycothiol (MSH) to buffer the intracellular redox environment. The NADPH-dependent FAD-containing oxidoreductase mycothiol disulfide reductase (Mtr) is known to reduce oxidized mycothiol disulfide (MSSM) to MSH, which is crucial to maintain the cellular redox balance. In this work, the first crystal structures of Mtr are presented, expanding the structural knowledge and understanding of LMW thiol reductases. The structural analyses and docking calculations provide insight into the nature of Mtrs, with regard to the binding and reduction of the MSSM substrate, in the context of related oxidoreductases. The putative binding site for MSSM suggests a similar binding to that described for the homologous glutathione reductase and its respective substrate glutathione disulfide, but with distinct structural differences shaped to fit the bulkier MSSM substrate, assigning Mtrs as uniquely functioning reductases. As MSH has been acknowledged as an attractive antitubercular target, the structural findings presented in this work may contribute towards future antituberculosis drug development.

Comprehensive analysis of the aldehyde dehydrogenase gene family in Phaseolus vulgaris L. and their response to saline-alkali stress

Background

Aldehyde dehydrogenase (ALDH) scavenges toxic aldehyde molecules by catalyzing the oxidation of aldehydes to carboxylic acids. Although ALDH gene family members in various plants have been extensively studied and were found to regulate plant response to abiotic stress, reports on ALDH genes in the common bean (Phaseolus vulgaris L.) are limited. In this study, we aimed to investigate the effects of neutral (NS) and basic alkaline (AS) stresses on growth, physiological and biochemical indices, and ALDH activity, ALDH gene expression of common bean. In addition, We used bioinformatics techniques to analyze the physical and chemical properties, phylogenetic relationships, gene replication, collinearity, cis-acting elements, gene structure, motifs, and protein structural characteristics of PvALDH family members.

Results

We found that both NS and AS stresses weakened the photosynthetic performance of the leaves, induced oxidative stress, inhibited common bean growth, and enhanced the antioxidative system to scavenge reactive oxygen species. Furthermore, we our findings revealed that ALDH in the common bean actively responds to NS or AS stress by inducing the expression of PvALDH genes. In addition, using the established classification criteria and phylogenetic analysis, 27 PvALDHs were identified in the common bean genome, belonging to 10 ALDH families. The primary expansion mode of PvALDH genes was segmental duplication. Cis-acting elemental analysis showed that PvALDHs were associated with abiotic stress and phytohormonal responses. Gene expression analysis revealed that the PvALDH gene expression was tissue-specific. For instance, PvALDH3F1 and PvALDH3H1 were highly expressed in flower buds and flowers, respectively, whereas PvALDH3H2 and PvALDH2B4 were highly expressed in green mature pods and young pods, respectively. PvALDH22A1 and PvALDH11A2 were highly expressed in leaves and young trifoliates, respectively; PvALDH18B2 and PvALDH18B3 were highly expressed in stems and nodules, respectively; and PvALDH2C2 and PvALDH2C3 were highly expressed in the roots. PvALDHs expression in the roots responded positively to NS-AS stress, and PvALDH2C3, PvALDH5F1, and PvALDH10A1 were significantly (P < 0.05) upregulated in the roots.

Conclusion

These results indicate that AS stress causes higher levels of oxidative damage than NS stress, resulting in weaker photosynthetic performance and more significant inhibition of common bean growth. The influence of PvALDHs potentially modulates abiotic stress response, particularly in the context of saline-alkali stress. These findings establish a basis for future research into the potential roles of ALDHs in the common bean.

The Mediterranean mussel Mytilus galloprovincialis: a novel model for developmental studies in mollusks

A model organism in developmental biology is defined by its experimental amenability and by resources created for the model system by the scientific community. For the most powerful invertebrate models, the combination of both has already yielded a thorough understanding of developmental processes. However, the number of developmental model systems is still limited, and their phylogenetic distribution heavily biased. Members of one of the largest animal lineages, the Spiralia, for example, have long been neglected. In order to remedy this shortcoming, we have produced a detailed developmental transcriptome for the bivalve mollusk Mytilus galloprovincialis, and have expanded the list of experimental protocols available for this species. Our high-quality transcriptome allowed us to identify transcriptomic signatures of developmental progression and to perform a first comparison with another bivalve mollusk: the Pacific oyster Crassostrea gigas. To allow co-labelling studies, we optimized and combined protocols for immunohistochemistry and hybridization chain reaction to create high-resolution co-expression maps of developmental genes. The resources and protocols described here represent an enormous boost for the establishment of Mytilus galloprovincialis as an alternative model system in developmental biology.

Chromosome-scale genome assembly of Lepus oiostolus (Lepus, Leporidae)

Lepus oiostolus (L. oiostolus) is a species endemic to the Qinghai-Tibet Plateau. However, the absence of a reference genome limits genetic studies. Here, we reported a high-quality L. oiostolus genome assembly, with scaffolds anchored to 24 chromosomes and a total assembled length of 2.80 Gb (contig N50 = 64.25 Mb). Genomic annotation uncovered 22,295 protein-coding genes and identified 49.84% of the sequences as transposable elements. Long interspersed nuclear elements (LINEs) constitute a high proportion of the genome. Our study is at the first time to report the chromosome-scale genome for the species of the L. oiostolus. It provides a valuable genomic resource for future research on the evolution of the Leporidae.

Molecular Insights into the Synergistic Effects of Putrescine and Ammonium on Dinoflagellates

Ammonium and polyamines are essential nitrogen metabolites in all living organisms. Crosstalk between ammonium and polyamines through their metabolic pathways has been demonstrated in plants and animals, while no research has been directed to explore this relationship in algae or to investigate the underlying molecular mechanisms. Previous research demonstrated that high concentrations of ammonium and putrescine were among the active substances in bacteria-derived algicide targeting dinoflagellates, suggesting that the biochemical inter-connection and/or interaction of these nitrogen compounds play an essential role in controlling these ecologically important algal species. In this research, putrescine, ammonium, or a combination of putrescine and ammonium was added to cultures of three dinoflagellate species to explore their effects. The results demonstrated the dose-dependent and species-specific synergistic effects of putrescine and ammonium on these species. To further explore the molecular mechanisms behind the synergistic effects, transcriptome analysis was conducted on dinoflagellate Karlodinium veneficum treated with putrescine or ammonium vs. a combination of putrescine and ammonium. The results suggested that the synergistic effects of putrescine and ammonium disrupted polyamine homeostasis and reduced ammonium tolerance, which may have contributed to the cell death of K. veneficum. There was also transcriptomic evidence of damage to chloroplasts and impaired photosynthesis of K. veneficum. This research illustrates the molecular mechanisms underlying the synergistic effects of the major nitrogen metabolites, ammonium and putrescine, in dinoflagellates and provides direction for future studies on polyamine biology in algal species.

Genome-wide molecular evolution analysis of the GRF and GIF gene families in Plantae (Archaeplastida)

BACKGROUND

Plant growth-regulating factors (GRFs) and GRF-interacting factors (GIFs) interact with each other and collectively have important regulatory roles in plant growth, development, and stress responses. Therefore, it is of great significance to explore the systematic evolution of GRF and GIF gene families. However, our knowledge and understanding of the role of GRF and GIF genes during plant evolution has been fragmentary.

RESULTS

In this study, a large number of genomic and transcriptomic datasets of algae, mosses, ferns, gymnosperms and angiosperms were used to systematically analyze the evolution of GRF and GIF genes during the evolution of plants. The results showed that GRF gene first appeared in the charophyte Klebsormidium nitens, whereas the GIF genes originated relatively early, and these two gene families were mainly expanded by segmental duplication events after plant terrestrialization. During the process of evolution, the protein sequences and functions of GRF and GIF family genes are relatively conservative. As cooperative partner, GRF and GIF genes contain the similar types of cis-acting elements in their promoter regions, which enables them to have similar transcriptional response patterns, and both show higher levels of expression in reproductive organs and tissues and organs with strong capacity for cell division. Based on protein-protein interaction analysis and verification, we found that the GRF-GIF protein partnership began to be established in pteridophytes and is highly conserved across different terrestrial plants.

CONCLUSIONS

These results provide a foundation for further exploration of the molecular evolution and biological functions of GRF and GIF genes.

Structure of a tripartite protein complex that targets toxins to the type VII secretion system

Type VII secretion systems are membrane-embedded nanomachines used by Gram-positive bacteria to export effector proteins from the cytoplasm to the extracellular environment. Many of these effectors are polymorphic toxins comprised of an N-terminal Leu-x-Gly (LXG) domain of unknown function and a C-terminal toxin domain that inhibits the growth of bacterial competitors. In recent work, it was shown that LXG effectors require two cognate Lap proteins for T7SS-dependent export. Here, we present the 2.6 Å structure of the LXG domain of the TelA toxin from the opportunistic pathogen Streptococcus intermedius in complex with both of its cognate Lap targeting factors. The structure reveals an elongated α-helical bundle within which each Lap protein makes extensive hydrophobic contacts with either end of the LXG domain. Remarkably, despite low overall sequence identity, we identify striking structural similarity between our LXG complex and PE-PPE heterodimers exported by the distantly related ESX type VII secretion systems of Mycobacteria implying a conserved mechanism of effector export among diverse Gram-positive bacteria. Overall, our findings demonstrate that LXG domains, in conjunction with their cognate Lap targeting factors, represent a tripartite secretion signal for a widespread family of T7SS toxins.

Altered microbiota, antimicrobial resistance genes, and functional enzyme profiles in the rumen of yak calves fed with milk replacer

IMPORTANCE

Yaks, as ruminants inhabiting high-altitude environments, possess a distinct rumen microbiome and are resistant to extreme living conditions. This study investigated the microbiota, resistome, and functional gene profiles in the rumen of yaks fed milk or milk replacer (MR), providing insights into the regulation of the rumen microbiome and the intervention of antimicrobial resistance in yaks through dietary methods. The abundance of Prevotella members increased significantly in response to MR. Tetracycline resistance was the most predominant. The rumen of yaks contained multiple antimicrobial resistance genes (ARGs) originating from different bacteria, which could be driven by MR, and these ARGs displayed intricate and complex interactions. MR also induced changes in functional genes. The enzymes associated with fiber degradation and butyrate metabolism were activated and showed close correlations with Prevotella members and butyrate concentration. This study allows us to deeply understand the ruminal microbiome and ARGs of yaks and their relationship with rumen bacteria in response to different milk sources.

Multiple Origins of Bioluminescence in Beetles and Evolution of Luciferase Function

Bioluminescence in beetles has long fascinated biologists, with diverse applications in biotechnology. To date, however, our understanding of its evolutionary origin and functional variation mechanisms remains poor. To address these questions, we obtained high-quality reference genomes of luminous and nonluminous beetles in 6 Elateroidea families. We then reconstructed a robust phylogenetic relationship for all luminous families and related nonluminous families. Comparative genomic analyses and biochemical functional experiments suggested that gene evolution within Elateroidea played a crucial role in the origin of bioluminescence, with multiple parallel origins observed in the luminous beetle families. While most luciferase-like proteins exhibited a conserved nonluminous amino acid pattern (TLA346 to 348) in the luciferin-binding sites, luciferases in the different luminous beetle families showed divergent luminous patterns at these sites (TSA/CCA/CSA/LVA). Comparisons of the structural and enzymatic properties of ancestral, extant, and site-directed mutant luciferases further reinforced the important role of these sites in the trade-off between acyl-CoA synthetase and luciferase activities. Furthermore, the evolution of bioluminescent color demonstrated a tendency toward hypsochromic shifts and variations among the luminous families. Taken together, our results revealed multiple parallel origins of bioluminescence and functional divergence within the beetle bioluminescent system.

Genome-wide analysis of the SMXL gene family in common bean and identification of karrikin-responsive PvSMXL2 as a negative regulator of PEG-induced drought stress

Common bean (Phaseolus vulgaris L.) is a major legume crop worldwide, but its growth and development frequently face challenges due to abiotic stresses, particularly drought. Proper supplement of copper could mitigate the adverse effects of drought, but excessive accumulation of this metal in plants can be harmful. The suppressor of MAX2 1-like (SMXL) gene family, which plays important roles in various plant processes, including stress responses, remains poorly understood in common bean. In this study, we identified nine orthologues of SMXL genes in common bean, which are located on six chromosomes and classified into four subgroups. Basic molecular properties, including theoretical isoelectric point (PI), molecular weight (MW), grand average of hydropathicity (GVIO), gene structure, and conserved motifs were characterized, and numerous cis-elements in promoters were predicted. The expression patterns of PvSMXL genes were found to be distinct under 10% polyethylene glycol (PEG)-induced drought stress and 200 μM Cu treatments. Most PvSMXLs showed reduced expression in response to Cu treatment, whereas nearly half PvSMXLs exhibited inducible expression under drought stress. PvSMXL2, which exhibited a rapid response to karrikin 1 (KAR1), an active form of the plant growth regulators newly found in the smoke of burning plant material, was down-regulated by both PEG-induced drought and Cu stresses. Transient silencing of PvSMXL2 resulted in enhanced drought stress tolerance without conferring Cu tolerance. These findings provide valuable insights into the functions of SMXL genes in common bean under abiotic stress conditions.

Genome-wide identification and comprehensive analysis heat shock transcription factor (Hsf) members in asparagus (Asparagus officinalis) at the seeding stage under abiotic stresses

Heat shock transcription factors (Hsf) are pivotal as essential transcription factors. They function as direct transcriptional activators of genes regulated by thermal stress and are closely associated with various abiotic stresses. Asparagus (Asparagus officinalis) is a vegetable of considerable economic and nutritional significance, abundant in essential vitamins, minerals, and dietary fiber. Nevertheless, asparagus is sensitive to environmental stresses, and specific abiotic stresses harm its yield and quality. In this context, Hsf members have been discerned through the reference genome, and a comprehensive analysis encompassing physical and chemical attributes, evolutionary aspects, motifs, gene structure, cis-acting elements, collinearity, and expression patterns under abiotic stresses has been conducted. The findings identified 18 members, categorized into five distinct subgroups. Members within each subgroup exhibited analogous motifs, gene structures, and cis-acting elements. Collinearity analysis unveiled a noteworthy pattern, revealing that Hsf members within asparagus shared one, two, and three pairs with counterparts in Arabidopsis, Oryza sativa, and Glycine max, respectively.Furthermore, members displayed tissue-specific expression during the seedling stage, with roots emerging as viable target tissue. Notably, the expression levels of certain members underwent modification under the influence of abiotic stresses. This study establishes a foundational framework for understanding Hsf members and offers valuable insights into the potential application of molecular breeding in the context of asparagus cultivation.

Duplication and neofunctionalization of a horizontally-transferred xyloglucanase as a facet of the red queen co-evolutionary dynamic

Oomycetes are heterotrophic protists that share phenotypic similarities with fungi, including the ability to cause plant diseases, but branch in a separate and distant region of the eukaryotic tree of life. It has been suggested that multiple horizontal gene transfers (HGTs) from fungi-to-oomycetes contributed to the evolution of plant-pathogenic traits. These HGTs are predicted to include secreted proteins that degrade plant cell walls. This is a key trait in the pathology of many oomycetes, as the plant cell wall represents a primary barrier to pathogen invasion and a rich source of carbohydrates. Many of the HGT gene families identified have undergone multiple rounds of duplication. Using a combination of phylogenomic analysis and functional assays, we investigate the diversification of a horizontally-transferred xyloglucanase gene family in the model oomycete species Phytophthora sojae. Our analyses detect 11 genes retained in P. sojae among a complex pattern of gene duplications and losses. Using a phenotype assay, based on heterologous expression in yeast, we show that eight of these paralogs have xyloglucanase function, including variants with distinct protein characteristics, such as a long-disordered C-terminal extension that can increase xyloglucanase activity. The functional xyloglucanase variants analysed subtend an ancestral node close to the fungi-oomycetes gene transfer, suggesting the horizontally-transferred gene was a bona fide xyloglucanase. Expression of xyloglucanase paralogs in Nicotiana benthamiana triggers distinct patterns of reactive oxygen species (ROS) generation, demonstrating that enzyme variants differentially stimulate pattern-triggered immunity in plants. Mass spectrometry of detectable enzymatic products demonstrates that some paralogs catalyze production of variant breakdown profiles, suggesting that secretion of multiple xyloglucanase variants increases efficiency of xyloglucan breakdown, as well as potentially diversifying the range of Damage-Associated Molecular Patterns (DAMPs) released during pathogen attack. We suggest that such patterns of protein neofunctionalization, and variant host responses, represent an aspect of the Red Queen host-pathogen co-evolutionary dynamic.

Significance Statement

The oomycetes are a diverse group of eukaryotic microbes that include some of the most devastating pathogens of plants. Oomycetes perceive, invade, and colonize plants in similar ways to fungi, in part because they acquired the genes to attack and feed on plants from fungi. These genes are predicted to be useful to oomycete plant pathogens because they have undergone multiple rounds of gene duplication. One key enzyme for attacking plant cell wall structures is called xyloglucanase. Xyloglucanase in the oomycetes has undergone multiple rounds of gene duplication, leading to variants including an enzyme with a C-terminal extension that increases activity. Some xyloglucanase variants trigger unique patterns of reactive oxygen species (ROS) in planta, and generate different profiles of cell wall breakdown products - such outcomes could act to mystify and increase the workload of the plant immune system, allowing successful pathogens to proliferate.

Comparative transcriptome analysis of two Daphnia galeata genotypes displaying contrasting phenotypic variation induced by fish kairomones in the same environment of the Han River, Korea

BACKGROUND

Phenotypic plasticity is a crucial adaptive mechanism that enables organisms to modify their traits in response to changes in their environment. Predator-induced defenses are an example of phenotypic plasticity observed across a wide range of organisms, from single-celled organisms to vertebrates. In addition to morphology and behavior, these responses also affect life-history traits. The crustacean Daphnia galeata is a suitable model organism for studying predator-induced defenses, as it exhibits life-history traits changes under predation risk. To get a better overview of their phenotypic plasticity under predation stress, we conducted RNA sequencing on the transcriptomes of two Korean Daphnia galeata genotypes, KE1, and KB11, collected in the same environment.

RESULTS

When exposed to fish kairomones, the two genotypes exhibited phenotypic variations related to reproduction and growth, with opposite patterns in growth-related phenotypic variation. From both genotypes, a total of 135,611 unigenes were analyzed, of which 194 differentially expressed transcripts (DETs) were shared among the two genotypes under predation stress, which showed consistent, or inconsistent expression patterns in both genotypes. Prominent DETs were related to digestion and reproduction and consistently up-regulated in both genotypes, thus associated with changes in life-history traits. Among the inconsistent DETs, transcripts encode vinculin (VINC) and protein obstructor-E (OBST-E), which are associated with growth; these may explain the differences in life-history traits between the two genotypes. In addition, genotype-specific DETs could explain the variation in growth-related life-history traits between genotypes, and could be associated with the increased body length of genotype KE1.

CONCLUSIONS

The current study allows for a better understanding of the adaptation mechanisms related to reproduction and growth of two Korean D. galeata genotypes induced by predation stress. However, further research is necessary to better understand the specific mechanisms by which the uncovered DETs are related with the observed phenotypic variation in each genotype. In the future, we aim to unravel the precise adaptive mechanisms underlying predator-induced responses.

Screening and Functional Analyses of Novel Cecropins from Insect Transcriptome

Antibiotic resistance is a significant and growing threat to global public health. However, antimicrobial peptides (AMPs) have shown promise as they exhibit a broad spectrum of antibacterial activities with low potential for resistance development. Insects, which inhabit a wide range of environments and are incredibly diverse, remain largely unexplored as a source of novel AMPs. To address this, we conducted a screening of the representative transcriptomes from the 1000 Insect Transcriptome Evolution (1KITE) dataset, focusing on the homologous reference genes of Cecropins, the first identified AMPs in insects known for its high efficiency. Our analysis identified 108 Cecropin genes from 105 insect transcriptomes, covering all major hexapod lineages. We validated the gene sequences and synthesized mature peptides for three identified Cecropin genes. Through minimal inhibition concentration and agar diffusion assays, we confirmed that these peptides exhibited antimicrobial activities against Gram-negative bacteria. Similar to the known Cecropin, the three Cecropins adopted an alpha-helical conformation in membrane-like environments, efficiently disrupting bacterial membranes through permeabilization. Importantly, none of the three Cecropins demonstrated cytotoxicity in erythrocyte hemolysis tests, suggesting their safety in real-world applications. Overall, this newly developed methodology provides a high-throughput bioinformatic pipeline for the discovery of AMP, taking advantage of the expanding genomic resources available for diverse organisms.

Predicting and Interpreting Protein Developability Via Transfer of Convolutional Sequence Representation

Engineered proteins have emerged as novel diagnostics, therapeutics, and catalysts. Often, poor protein developability─quantified by expression, solubility, and stability─hinders utility. The ability to predict protein developability from amino acid sequence would reduce the experimental burden when selecting candidates. Recent advances in screening technologies enabled a high-throughput (HT) developability dataset for 105 of 1020 possible variants of protein ligand scaffold Gp2. In this work, we evaluate the ability of neural networks to learn a developability representation from a HT dataset and transfer this knowledge to predict recombinant expression beyond observed sequences. The model convolves learned amino acid properties to predict expression levels 44% closer to the experimental variance compared to a non-embedded control. Analysis of learned amino acid embeddings highlights the uniqueness of cysteine, the importance of hydrophobicity and charge, and the unimportance of aromaticity, when aiming to improve the developability of small proteins. We identify clusters of similar sequences with increased recombinant expression through nonlinear dimensionality reduction and we explore the inferred expression landscape via nested sampling. The analysis enables the first direct visualization of the fitness landscape and highlights the existence of evolutionary bottlenecks in sequence space giving rise to competing subpopulations of sequences with different developability. The work advances applied protein engineering efforts by predicting and interpreting protein scaffold expression from a limited dataset. Furthermore, our statistical mechanical treatment of the problem advances foundational efforts to characterize the structure of the protein fitness landscape and the amino acid characteristics that influence protein developability.

Analysis of Delta(9) fatty acid desaturase gene family and their role in oleic acid accumulation in Carya cathayensis kernel

Carya cathayensis, commonly referred to as Chinese hickory, produces nuts that contain high-quality edible oils, particularly oleic acid (18:1). It is known that stearoyl-ACP desaturase (SAD) is the first key step converting stearic acid (C18:0, SA) to oleic acid (C18:1, OA) in the aminolevulinic acid (ALA) biosynthetic pathway and play an important role in OA accumulation. Thus far, there is little information about SAD gene family in C. cathayensis and the role of individual members in OA accumulation. This study searched the Chinese Hickory Genome Database and identified five members of SAD genes, designated as CcSADs, at the whole genome level through the comparison with the homologous genes from Arabidopsis. RNA-Seq analysis showed that CcSSI2-1, CcSSI2-2, and CcSAD6 were highly expressed in kernels. The expression pattern of CcSADs was significantly correlated with fatty acid accumulation during the kernel development. In addition, five full-length cDNAs encoding SADs were isolated from the developing kernel of C. cathayensis. CcSADs-green fluorescent protein (GFP) fusion construct was infiltrated into tobacco epidermal cells, and results indicated their chloroplast localization. The catalytic function of these CcSADs was further analyzed by heterologous expression in Saccharomyces cerevisiae, Nicotiana benthamiana, and walnut. Functional analysis demonstrated that all CcSADs had fatty acid desaturase activity to catalyze oleic acid biosynthesis. Some members of CcSADs also have strong substrate specificity for 16:0-ACP to synthesize palmitoleic acid (C16:1, PA). Our study documented SAD gene family in C. cathayensis and the role of CcSSI2-1, CcSSI2-2, and CcSAD6 in OA accumulation, which could be important for future improvement of OA content in this species via genetic manipulation.

A systematic approach to classify and characterize genomic islands driven by conjugative mobility using protein signatures

Genomic islands (GIs) play a crucial role in the spread of antibiotic resistance, virulence factors and antiviral defense systems in a broad range of bacterial species. However, the characterization and classification of GIs are challenging due to their relatively small size and considerable genetic diversity. Predicting their intercellular mobility is of utmost importance in the context of the emerging crisis of multidrug resistance. Here, we propose a large-scale classification method to categorize GIs according to their mobility profile and, subsequently, analyze their gene cargo. We based our classification decision scheme on a collection of mobility protein motif definitions available in publicly accessible databases. Our results show that the size distribution of GI classes correlates with their respective structure and complexity. Self-transmissible GIs are usually the largest, except in Bacillota and Actinomycetota, accumulate antibiotic and phage resistance genes, and favour the use of a tyrosine recombinase to insert into a host's replicon. Non-mobilizable GIs tend to use a DDE transposase instead. Finally, although tRNA genes are more frequently targeted as insertion sites by GIs encoding a tyrosine recombinase, most GIs insert in a protein-encoding gene. This study is a stepping stone toward a better characterization of mobile GIs in bacterial genomes and their mechanism of mobility.

Genome-Wide Analysis and Expression Profiling of DUF668 Genes in Glycine max under Salt Stress

The DUF668 gene performs a critical role in mitigating the impact of abiotic stress factors. In this study, we identified 30 DUF668 genes in a soybean genome, distributed across fifteen chromosomes. The phylogenetic analysis classified the DUF668 genes into three groups (group I, group II, and group III). Interestingly, gene structure analysis illustrated that several GmDUF668 genes were without introns. Furthermore, the subcellular localization results suggested that GmDUF668 proteins were present in the nucleus, mitochondria, cytoplasm, and plasma membrane. GmDUF668 promoters were analyzed in silico to gain insight into the presence of regulatory sequences for TFs binding. The expression profiling illustrated that GmDUF668 genes showed expression in leaves, roots, nodules, and flowers. To investigate their response to salt stress, we utilized the RNA sequencing data of GmDUF668 genes. The results unveiled that GmDUF668-8, GmDUF668-20, and GmDUF668-30 genes were upregulated against salt stress treatment. We further validated these findings using qRT-PCR analysis. These findings provide a scientific basis to explore the functions of GmDUF668 genes against different stress conditions.

Antibiotic potentiation and inhibition of cross-resistance in pathogens associated with cystic fibrosis

Critical Gram-negative pathogens, like Pseudomonas, Stenotrophomonas and Burkholderia, have become resistant to most antibiotics. Complex resistance profiles together with synergistic interactions between these organisms increase the likelihood of treatment failure in distinct infection settings, for example in the lungs of cystic fibrosis patients. Here, we discover that cell envelope protein homeostasis pathways underpin both antibiotic resistance and cross-protection in CF-associated bacteria. We find that inhibition of oxidative protein folding inactivates multiple species-specific resistance proteins. Using this strategy, we sensitize multi-drug resistant Pseudomonas aeruginosa to β-lactam antibiotics and demonstrate promise of new treatment avenues for the recalcitrant pathogen Stenotrophomonas maltophilia. The same approach also inhibits cross-protection between resistant S. maltophilia and susceptible P. aeruginosa, allowing eradication of both commonly co-occurring CF-associated organisms. Our results provide the basis for the development of next-generation strategies that target antibiotic resistance, while also impairing specific interbacterial interactions that enhance the severity of polymicrobial infections.

The genetic structure and demographic history revealed by whole-genome resequencing provide insights into conservation of critically endangered Artocarpus nanchuanensis

Introduction

Whole-genome resequencing technology covers almost all nucleotide variations in the genome, which makes it possible to carry out conservation genomics research on endangered species at the whole-genome level.

Methods

In this study, based on the whole-genome resequencing data of 101 critically endangered Artocarpus nanchuanensis individuals, we evaluated the genetic diversity and population structure, inferred the demographic history and genetic load, predicted the potential distributions in the past, present and future, and classified conservation units to propose targeted suggestions for the conservation of this critically endangered species.

Results

Whole-genome resequencing for A. nanchuanensis generated approximately 2 Tb of data. Based on abundant mutation sites (25,312,571 single nucleotide polymorphisms sites), we revealed that the average genetic diversity (nucleotide diversity, π) of different populations of A. nanchuanensis was relatively low compared with other trees that have been studied. And we also revealed that the NHZ and QJT populations harboured unique genetic backgrounds and were significantly separated from the other five populations. In addition, positive genetic selective signals, significantly enriched in biological processes related to terpene synthesis, were identified in the NHZ population. The analysis of demographic history of A. nanchuanensis revealed the existence of three genetic bottleneck events. Moreover, abundant genetic loads (48.56% protein-coding genes) were identified in Artocarpus nanchuanensis, especially in genes related to early development and immune function of plants. The predication analysis of suitable habitat areas indicated that the past suitable habitat areas shifted from the north to the south due to global temperature decline. However, in the future, the actual distribution area of A. nanchuanensis will still maintain high suitability.

Discussion

Based on total analyses, we divided the populations of A. nanchuanensis into four conservation units and proposed a number of practical management suggestions for each conservation unit. Overall, our study provides meaningful guidance for the protection of A. nanchuanensis and important insight into conservation genomics research.

A Candidate Bacterial-Type Amino Acid Decarboxylase Is Essential for Male Gamete Exflagellation and Mosquito Transmission of the Malaria Parasite

A frequent side effect of chemotherapy against malaria parasite blood infections is a dramatic induction of the sexual blood stages, thereby enhancing the risk of future malaria transmissions. The polyamine biosynthesis pathway has been suggested as a candidate target for transmission-blocking anti-malarial drug development. Herein, we describe the role of a bacterial-type amino acid decarboxylase (AAD) in the life cycle of the malaria model parasite Plasmodium yoelii. Hallmarks of AAD include a conserved catalytic lysine residue and high-level homology to arginine/lysine/ornithine decarboxylases of pathogenic bacteria. By targeted gene deletion, we show that AAD plays an essential role in the exflagellation of microgametes, resulting in complete absence of sporozoites in the mosquito vector. These data highlight the central role of the biosysthesis of polyamines in the final steps of male gamete sexual development of the malaria parasite and, hence, onward transmission to mosquitoes.

F-Type Pyocins Are Diverse Noncontractile Phage Tail-Like Weapons for Killing Pseudomonas aeruginosa

Most Pseudomonas aeruginosa strains produce bacteriocins derived from contractile or noncontractile phage tails known as R- and F-type pyocins, respectively. These bacteriocins possess strain-specific bactericidal activity against P. aeruginosa and likely increase evolutionary fitness through intraspecies competition. R-type pyocins have been studied extensively and show promise as alternatives to antibiotics. Although they have similar therapeutic potential, experimental studies on F-type pyocins are limited. Here, we provide a bioinformatic and experimental investigation of F-type pyocins. We introduce a systematic naming scheme for genes found in R- and F-type pyocin operons and identify 15 genes invariably found in strains producing F-type pyocins. Five proteins encoded at the 3' end of the F-type pyocin cluster are divergent in sequence and likely determine bactericidal specificity. We use sequence similarities among these proteins to define eleven distinct F-type pyocin groups, five of which had not been previously described. The five genes encoding the variable proteins associate in two modules that have clearly reassorted independently during the evolution of these operons. These proteins are considerably more diverse than the specificity-determining tail fibers of R-type pyocins, suggesting that F-type pyocins may have emerged earlier. Experimental studies on six F-type pyocin groups show that each displays a distinct spectrum of bactericidal activity. This activity is strongly influenced by the lipopolysaccharide O-antigen type, but other factors also play a role. F-type pyocins appear to kill as efficiently as R-type pyocins. These studies set the stage for the development of F-type pyocins as antibacterial therapeutics. IMPORTANCE Pseudomonas aeruginosa is an opportunistic pathogen that causes antibiotic-resistant infections with high mortality rates, particularly in immunocompromised individuals and cystic fibrosis patients. Due to the increasing frequency of multidrug-resistant P. aeruginosa infections, there is great need for the development of alternative therapeutics. In this study, we investigate one such potential therapeutic: F-type pyocins, which are bacteriocins naturally produced by P. aeruginosa that resemble noncontractile phage tails. We show that they are potent killers of P. aeruginosa and identify their probable bactericidal specificity determinants, which opens up the possibility of engineering them to precisely target strains of pathogenic bacteria. The resemblance of F-type pyocins to well-characterized phage tails will greatly facilitate their development into effective antibacterials.

Allele-aware chromosome-level genome assembly of the autohexaploid Diospyros kaki Thunb

Artificially improving persimmon (Diospyros kaki Thunb.), one of the most important fruit trees, remains challenging owing to the lack of reference genomes. In this study, we generated an allele-aware chromosome-level genome assembly for the autohexaploid persimmon 'Xiaoguotianshi' (Chinese-PCNA type) using PacBio CCS and Hi-C technology. The final assembly contained 4.52 Gb, with a contig N50 value of 5.28 Mb and scaffold N50 value of 44.01 Mb, of which 4.06 Gb (89.87%) of the assembly were anchored onto 90 chromosome-level pseudomolecules comprising 15 homologous groups with 6 allelic chromosomes in each. A total of 153,288 protein-coding genes were predicted, of which 98.60% were functionally annotated. Repetitive sequences accounted for 64.02% of the genome; and 110,480 rRNAs, 12,297 tRNAs, 1,483 miRNAs, and 3,510 snRNA genes were also identified. This genome assembly fills the knowledge gap in the autohexaploid persimmon genome, which is conducive in the study on the regulatory mechanisms underlying the major economically advantageous traits of persimmons and promoting breeding programs.

Identification of nsLTP family in Chinese white pear (Pyrus bretschneideri) reveals its potential roles in russet skin formation

MAIN CONCLUSION

Identification of PbLTP genes in pear and functional characterization of PbLTP4 in the transport of suberin monomers of russet skin formation. Non-specific lipid-transfer protein (nsLTP) is an abundant and diverse alkaline small molecule protein in the plant kingdom with complex and diverse biophysiological functions, such as transfer of phospholipids, reproductive development, pathogen defence and abiotic stress response. Up to now, only a tiny fraction of nsLTPs have been functionally identified, and the distribution of nsLTPs in pear (Pyrus bretschneideri) (PbLTPs) has not been fully characterized. In this study, the genome-wide analysis of the nsLTP gene family in the pear genome identified 67 PbLTP proteins, which could be divided into six types (1, 2, C, D, E, and G). Similar intron/exon structural patterns were observed in the same type, strongly supporting their close evolutionary relationship. In addition, PbLTP4 was highly expressed in russet pear skin compared with green skin, which was located in the plasma membrane. Coexpression network analysis showed that PbLTP4 closely related to suberin biosynthetic genes. The biological function of PbLTP4 in promoting suberification has been demonstrated by overexpression in Arabidopsis. Identification of suberin monomers showed that PbLTP4 promotes suberification by regulating 9,12-octadecadienoic acid and hexadecanoic acid transport. These results provide helpful insights into the characteristics of PbLTP genes and their biological function in the transport of suberin monomers of russet skin formation.

The stability of transcription factor PfSPL1 participates in the response to phytoplasma stress in Paulownia fortunei

The current understanding of the pathogenesis of phytoplasma is still very limited and challenging. Here, ceRNA regulatory network and degradome sequencing identified a PfmiR156f-PfSPL regulatory module in Paulownia fortunei infected by phytoplasma, and RLM-5'RACE and dual luciferase analyses verified the relationship. The PfmiR156 cleavage site was located at 1104 nt and 1177 nt of PfSPL1 and PfSPL10, respectively. MG132 and epoxomicin, two 26S proteasome inhibitors, significantly increased the accumulation of PfSPL1. PfSPL1 was also the attack target of phytoplasma effectors (Pawb 3/9/16/37/51) after the phytoplasma invaded Paulownia. Moreover, molecular docking implied that the effectors may interact with the conserved SBP domain of the target protein PfSPL1. Basically, these results indicated that the stability of PfSPL1 was regulated by PfmiR156 cleavage activity and/or the 26S proteasome pathway at the post-translation level. The PfSPL1, which is a transcription factor, was also the one of the targets of multiple effectors attacking Paulownia. This study provides a good scope to understand the paulownia phytoplasma infecting mechanism.

Polyphenol oxidases regulate pollen development through modulating flavonoids homeostasis in tobacco

Pollen development is critical in plant reproduction. Polyphenol oxidases (PPOs) genes encode defense-related enzymes, but the role of PPOs in pollen development remains largely unexplored. Here, we characterized NtPPO genes, and then investigated their function in pollen via creating NtPPO9/10 double knockout mutant (cas-1), overexpression 35S::NtPPO10 (cosp) line and RNAi lines against all NtPPOs in Nicotiana tabacum. NtPPOs were abundantly expressed in the anther and pollen (especially NtPPO9/10). The pollen germination, polarity ratio and fruit weights were significantly reduced in the NtPPO-RNAi and cosp lines, while they were normal in cas-1 likely due to compensation by other NtPPO isoforms. Comparisons of metabolites and transcripts between the pollen of WT and NtPPO-RNAi, or cosp showed that decreased enzymatic activity of NtPPOs led to hyper-accumulation of flavonoids. This accumulation might reduce the content of ROS. Ca²⁺ and actin levels also decreased in pollen of the transgenic lines.Thus, the NtPPOs regulate pollen germination through the flavonoid homeostasis and ROS signal pathway. This finding provides novel insights into the native physiological functions of PPOs in pollen during reproduction.

The first chromosome-level Fallopia multiflora genome assembly provides insights into stilbene biosynthesis

Fallopia multiflora (Thunb.) Harald, a vine belonging to the Polygonaceae family, is used in traditional medicine. The stilbenes contained in it have significant pharmacological activities in anti-oxidation and anti-aging. This study describes the assembly of the F. multiflora genome and presents its chromosome-level genome sequence containing 1.46 gigabases of data (with a contig N50 of 1.97 megabases), 1.44 gigabases of which was assigned to 11 pseudochromosomes. Comparative genomics confirmed that F. multiflora shared a whole-genome duplication event with Tartary buckwheat and then underwent different transposon evolution after separation. Combining genomics, transcriptomics, and metabolomics data to map a network of associated genes and metabolites, we identified two FmRS genes responsible for the catalysis of one molecule of p-coumaroyl-CoA and three molecules of malonyl-CoA to resveratrol in F. multiflora. These findings not only serve as the basis for revealing the stilbene biosynthetic pathway but will also contribute to the development of tools for increasing the production of bioactive stilbenes through molecular breeding in plants or metabolic engineering in microbes. Moreover, the reference genome of F. multiflora is a useful addition to the genomes of the Polygonaceae family.

Identification of F-box gene family in Brassica oleracea and expression analysis in response to low-temperature stress

Unfavorable climatic conditions, such as low temperatures, often hinder the growth and production of crops worldwide. The F-box protein-encoding gene family performs an essential role in plant stress resistance. However, a comprehensive analysis of the F-box gene family in cabbage (Brassica oleracea var capitata L.) has not been reported yet. In this study, genome-wide characterization of F-box proteins in cabbage yielded 303 BoFBX genes and 224 BoFBX genes unevenly distributed on 9 chromosomes of cabbage. Phylogenetic analysis of 303 BoFBX genes was classified into nine distinct subfamily groups (GI-GIX). Analysis of the gene structure of BoFBX genes indicated that most genes within the same clade are highly conserved. In addition, tissue-specific expression analysis revealed that six F-box genes in cabbage showed the highest expression in rosette leaves, followed by roots and stems and the lowest expression was observed in the BoFBX156 gene. In contrast, the expression of the other five genes, BoFBX100, BoFBX117, BoFBX136, BoFBX137 and BoFBX213 was observed to be upregulated in response to low-temperature stress. Moreover, we found that the expression level of the BoFBX gene in the cold-tolerant cultivar "ZG" was higher than that in cold-sensitive "YC" with the extension of stress duration, while expression levels of each gene in "ZG" were higher than "YC" at 24 h. Knowledge of the various functions provided by BoFBXs genes and their expression patterns provides a firm theoretical foundation for explaining the functions of BoFBXs, thereby contributing to the molecular breeding process of cabbage.

Mutating novel interaction sites in NRP1 reduces SARS-CoV-2 spike protein internalization

The global pandemic of coronavirus disease 2019 caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus has become a severe global health problem because of its rapid spread. Both Ace2 and NRP1 provide initial viral binding sites for SARS-CoV-2. Here, we show that cysteine residues located in the vestigial plasminogen-apple-nematode (PAN) domain of NRP1 are necessary for SARS-CoV-2 spike protein internalization. Mutating novel cysteine residues in the PAN altered NRP1 stability and downstream activation of extracellular signal-regulated kinase (ERK) signaling pathway and impaired its interaction with the spike protein. This resulted in a significant reduction in spike protein abundance in Vero-E6 cells for the original, alpha, and delta SARS-CoV-2 variants even in the presence of the Ace2. Moreover, mutating these cysteine residues in NRP1 significantly lowered its association with Plexin-A1. As the spike protein is a critical component for targeted therapy, our biochemical study may represent a distinct mechanism to develop a path for future therapeutic discovery.

Transcriptomic analysis reveals key genes and pathways corresponding to Cd and Pb in the hyperaccumulator Arabis paniculata

Soil and water are increasingly at risk of contamination from the toxic heavy metals lead (Pb) and cadmium (Cd). Arabis paniculata (Brassicaceae) is a hyperaccumulator of heavy metals (HMs) found widely distributed in areas impacts by mining activities. However, the mechanism by which A. paniculata tolerates HMs is still uncharacterized. For this experiment, we employed RNA sequencing (RNA-seq) in order to find Cd (0.25 mM)- and Pb (2.50 mM)-coresponsive genes A. paniculata. In total, 4490 and 1804 differentially expressed genes (DEGs) were identified in root tissue, and 955 and 2209 DEGs were identified in shoot tissue, after Cd and Pb exposure, respectively. Interestingly in root tissue, gene expression corresponded similarly to both Cd and Pd exposure, of which 27.48% were co-upregulated and 41.00% were co-downregulated. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses showed that the co-regulated genes were predominantly involved in transcription factors (TFs), cell wall biosynthesis, metal transport, plant hormone signal transduction, and antioxidant enzyme activity. Many critical Pb/Cd-induced DEGs involved in phytohormone biosynthesis and signal transduction, HM transport, and transcription factors were also identified. Especially the gene ABCC9 was co-downregulated in root tissues but co-upregulated in shoot tissues. The co-downregulation of ABCC9 in the roots prevented Cd and Pb from entering the vacuole rather than the cytoplasm for transporting HMs to shoots. While in shoots, the ABCC9 co-upregulated results in vacuolar Cd and Pb accumulation, which may explain why A. paniculata is a hyperaccumulator. These results will help to reveal the molecular and physiological processes underlying tolerance to HM exposure in the hyperaccumulator A. paniculata, and aid in future efforts to utilize this plant in phytoremediation.

The involvement of Ein3-binding F-box protein PbrEBF3 in regulating ethylene signaling during Cuiguan pear fruit ripening

Ein3-binding F-box (EBF) proteins have been determined to modulate ethylene response processes by regulating EIN3/EIL protein degradation in Arabidopsis and tomato. However, the function of pear PbrEBFs in ethylene-dependent responses during fruit ripening remains unclear. In this study, PbrEBF1, PbrEBF2, and PbrEBF3 display contrasting expression patterns in response to ethylene and 1-MCP treatment. PbrEBF3 displayed potential fruit ripening-associated function in a transient expression experiment. Yeast two-hybrid (Y2H) and Firefly luciferase complementation imaging (LCI) assays indicated that PbrEBF3 interacts with PbrEIL1, PbrEIL2, and PbrEIL3 proteins. In turn, the transcription of PbrEBF3 is directly regulated by PbrEILs via a feedback loop. PbrEILs trigger a transcriptional cascade of PbrERF24 and finally affect ethylene synthesis. Overall, PbrEBF3 plays a central role in pear fruit ripening through mediation of the ethylene signaling pathway.

Revealing evolution of tropane alkaloid biosynthesis by analyzing two genomes in the Solanaceae family

Tropane alkaloids (TAs) are widely distributed in the Solanaceae, while some important medicinal tropane alkaloids (mTAs), such as hyoscyamine and scopolamine, are restricted to certain species/tribes in this family. Little is known about the genomic basis and evolution of TAs biosynthesis and specialization in the Solanaceae. Here, we present chromosome-level genomes of two representative mTAs-producing species: Atropa belladonna and Datura stramonium. Our results reveal that the two species employ a conserved biosynthetic pathway to produce mTAs despite being distantly related within the nightshade family. A conserved gene cluster combined with gene duplication underlies the wide distribution of TAs in this family. We also provide evidence that branching genes leading to mTAs likely have evolved in early ancestral Solanaceae species but have been lost in most of the lineages, with A. belladonna and D. stramonium being exceptions. Furthermore, we identify a cytochrome P450 that modifies hyoscyamine into norhyoscyamine. Our results provide a genomic basis for evolutionary insights into the biosynthesis of TAs in the Solanaceae and will be useful for biotechnological production of mTAs via synthetic biology approaches.

Pangenomic analysis identifies structural variation associated with heat tolerance in pearl millet

Pearl millet is an important cereal crop worldwide and shows superior heat tolerance. Here, we developed a graph-based pan-genome by assembling ten chromosomal genomes with one existing assembly adapted to different climates worldwide and captured 424,085 genomic structural variations (SVs). Comparative genomics and transcriptomics analyses revealed the expansion of the RWP-RK transcription factor family and the involvement of endoplasmic reticulum (ER)-related genes in heat tolerance. The overexpression of one RWP-RK gene led to enhanced plant heat tolerance and transactivated ER-related genes quickly, supporting the important roles of RWP-RK transcription factors and ER system in heat tolerance. Furthermore, we found that some SVs affected the gene expression associated with heat tolerance and SVs surrounding ER-related genes shaped adaptation to heat tolerance during domestication in the population. Our study provides a comprehensive genomic resource revealing insights into heat tolerance and laying a foundation for generating more robust crops under the changing climate.

CrAssphage May Be Viable Markers of Contamination in Pristine and Contaminated River Water

Viruses are the most biologically abundant entities and may be ideal indicators of fecal pollutants in water. Anthropogenic activities have triggered drastic ecosystem changes in rivers, leading to substantial shifts in chemical and biological attributes. Here, we evaluate the viability of using the presence of crAssphage as indicators of fecal contamination in South African rivers. Shotgun analysis revealed diverse crAssphage viruses in these rivers, which are impacted by chemical and biological pollution. Overall, the diversity and relative abundances of these viruses was higher in contaminated sites compared to pristine locations. In contrast to fecal coliform counts, crAssphage sequences were detected in pristine rivers, supporting the assertion that the afore mentioned marker may be a more accurate indicator of fecal contamination. Our data demonstrate the presence of diverse putative hosts which includes members of the phyla Bacteroidota, Pseudomonadota, Verrucomicrobiota, and Bacillota. Phylogenetic analysis revealed novel subfamilies, suggesting that rivers potentially harbor distinct and uncharacterized clades of crAssphage. These data provide the first insights regarding the diversity, distribution, and functional roles of crAssphage in rivers. Taken together, the results support the potential application of crAssphage as viable markers for water quality monitoring. IMPORTANCE Rivers support substantial populations and provide important ecosystem services. Despite the application of fecal coliform tests and other markers, we lack rapid and reproducible approaches for determining fecal contamination in rivers. Waterborne viral outbreaks have been reported even after fecal indicator bacteria (FIB) were suggested to be absent or below regulated levels of coliforms. This indicates a need to develop and apply improved indicators of pollutants in aquatic ecosystems. Here, we evaluate the viability of crAssphage as indicators of fecal contamination in two South African rivers. We assess the abundance, distribution, and diversity of these viruses in sites that had been predicted pristine or contaminated by FIB analysis. We show that crAssphage are ideal and sensitive markers for fecal contamination and describe novel clades of crAss-like phages. Known crAss-like subfamilies were unrepresented in our data, suggesting that the diversity of these viruses may reflect geographic locality and dependence.

Coevolution-based prediction of key allosteric residues for protein function regulation

Allostery is fundamental to many biological processes. Due to the distant regulation nature, how allosteric mutations, modifications, and effector binding impact protein function is difficult to forecast. In protein engineering, remote mutations cannot be rationally designed without large-scale experimental screening. Allosteric drugs have raised much attention due to their high specificity and possibility of overcoming existing drug-resistant mutations. However, optimization of allosteric compounds remains challenging. Here, we developed a novel computational method KeyAlloSite to predict allosteric site and to identify key allosteric residues (allo-residues) based on the evolutionary coupling model. We found that protein allosteric sites are strongly coupled to orthosteric site compared to non-functional sites. We further inferred key allo-residues by pairwise comparing the difference of evolutionary coupling scores of each residue in the allosteric pocket with the functional site. Our predicted key allo-residues are in accordance with previous experimental studies for typical allosteric proteins like BCR-ABL1, Tar, and PDZ3, as well as key cancer mutations. We also showed that KeyAlloSite can be used to predict key allosteric residues distant from the catalytic site that are important for enzyme catalysis. Our study demonstrates that weak coevolutionary couplings contain important information of protein allosteric regulation function. KeyAlloSite can be applied in studying the evolution of protein allosteric regulation, designing and optimizing allosteric drugs, and performing functional protein design and enzyme engineering.

Staying below the Radar: Unraveling a New Family of Ubiquitous "Cryptic" Non-Tailed Temperate Vibriophages and Implications for Their Bacterial Hosts

Bacteriophages are the most abundant biological entities in the oceans and play key roles in bacterial activity, diversity and evolution. While extensive research has been conducted on the role of tailed viruses (Class: Caudoviricetes), very little is known about the distribution and functions of the non-tailed viruses (Class: Tectiliviricetes). The recent discovery of the lytic Autolykiviridae family demonstrated the potential importance of this structural lineage, emphasizing the need for further exploration of the role of this group of marine viruses. Here, we report the novel family of temperate phages under the class of Tectiliviricetes, which we propose to name "Asemoviridae" with phage NO16 as a main representative. These phages are widely distributed across geographical regions and isolation sources and found inside the genomes of at least 30 species of Vibrio, in addition to the original V. anguillarum isolation host. Genomic analysis identified dif-like sites, suggesting that NO16 prophages recombine with the bacterial genome based on the XerCD site-specific recombination mechanism. The interactions between the NO16 phage and its V. anguillarum host were linked to cell density and phage-host ratio. High cell density and low phage predation levels were shown to favor the temperate over the lytic lifestyle for NO16 viruses, and their spontaneous induction rate was highly variable between different V. anguillarum lysogenic strains. NO16 prophages coexist with the V. anguillarum host in a mutualistic interaction by rendering fitness properties to the host, such as increased virulence and biofilm formation through lysogenic conversion, likely contributing to their global distribution.