SMART, a simple modular architecture research tool: identification of signaling domains

Etiology, clinical characteristics, genetic profile, and outcomes of children with refractory rickets at a referral center in India: a cohort study

BACKGROUND

Limited research exists regarding the genetic profile, clinical characteristics, and outcomes of refractory rickets in children from India.

METHODS

Patients with refractory rickets aged ≤ 18 years were enrolled. Data regarding clinical features, etiology, genotype-phenotype correlation, and estimated glomerular filtration rate (eGFR) were recorded.

RESULTS

Seventy-two patients with refractory rickets (non-nutritional, with normal kidney function at presentation) from 65 families attending the pediatric nephrology clinic from 2005-2024 were included. Median (IQR) age at first presentation was 2 (1, 4) years. Clinical features included failure-to-thrive (49 [68.1%]), polyuria (37 [51.4%]), nephrocalcinosis (33 [45.8%]), fractures (10 [13.9%]), and hypokalemic paralysis (4 [5.6%]). Major etiologies included distal renal tubular acidosis (dRTA) [34(47.2%)], hereditary hypophosphatemic rickets (11 [15.3%]), cystinosis (9 [12.5%]), Lowe syndrome (3 [4.2%]), vitamin D-dependent rickets (4 [5.5%]), and Fanconi-Bickel syndrome (3 [4.2%]). Next-generation sequencing identified 61 variants among 71 children tested (85.9%), of which 56 variants (among 55 children) were pathogenic (P)/likely-pathogenic (LP) (77.5% diagnostic-yield). P/LP variants included SLC4A1 (n = 14), CTNS (n = 9), PHEX (n = 8), WDR72 (n = 5), OCRL (n = 2), SLC2A2 (n = 3), ATP6V0A4 (n = 4), VDR (n = 3), CLDN16 (n = 2), ATP6V1B1 (n = 1), SLC12A1 (n = 1), CLCN5 (n = 1), SLC34A3 (n = 1), ATP7B (n = 1), and KCNJ1 (n = 1). Fifteen novel P/LP variants and five novel variants-of-uncertain-significance (VUS) were identified. c.2573C > A in exon 19 among SLC4A1-dRTA (n = 14) was a recurrent mutation. Five patients with cystinosis, two patients with SLC4A1-dRTA, two with WDR72-dRTA, and two with Bartter syndrome showed progression to CKD stage 2 or greater during follow-up.

CONCLUSIONS

dRTA, X-linked hypophosphatemic rickets, and cystinosis were common causes of refractory rickets. The c.2573C > A variant in exon 19 was a recurrent mutation in SLC4A1-dRTA.

Cuticular proteins (crusticuls) affect 3D chitin bundle nanostructure

The crustacean exoskeleton features a micrometric, three-dimensional chitin scaffold. The intricate organization of this structure makes it an ideal model for investigating scaffold proteins at the nanoscale. Periodic exoskeleton replacement during a rapid and punctual molt cycle involves proteins that govern exoskeleton formation. Relying on binary expression pattern analysis of a molt-related transcriptomic library generated from the cuticle-forming epithelium of the crayfish Cherax quadricarinatus, a family of crustacean cuticle structural proteins termed 'crusticuls' was discovered and shown to present an exoskeleton formation-related expression pattern. All nine crusticuls include a chitin-binding domain bordered by two acidic residue-rich regions, putative functional domains related to exoskeletal formation and biomineralization. Crusticuls knock-down via RNAi resulted in over 95% reduced relative expression in treated versus control crayfish, with phenotypic effects ranging from prolonged molt cycles to lethality. Crusticuls were largely absent from newly formed cuticles following knockdown, resulting in exoskeletal deformities in the three-dimensional organization of chitinous bundles at the micro- and nanometric scales. These structural alterations were phenotypically translated into changes in cuticular hardness and elasticity. The identification of crusticuls as being key for proper nanometric three-dimensional organization of cuticular chitinous scaffolds opens new avenues for synthetic scaffold bio-mimetic applications.

Multi-omics insights into growth and fruiting body development in the entomopathogenic fungus Cordycepsblackwelliae

Cordycepsblackwelliae is an entomopathogenic fungus with significant potential for research and development due to its ease of cultivation. However, the lack of omics-based studies has limited our understanding of the molecular mechanisms governing its growth and fruiting body development. This study employed a multi-omics approach, integrating genomic, transcriptomic and metabolomic analyses. Utilising both Illumina and Nanopore sequencing technologies, we assembled a 31.06 Mb nuclear genome comprising 11 scaffolds, with telomere presence at one or both ends in eight scaffolds and annotated 8,138 identified genes (8,136 from genome prediction and two from local BLAST searches). Transcriptomic analysis identified 2,078 differentially expressed genes across three developmental stages: liquid culture mycelia, wheat culture mycelia and fruiting bodies. Amongst these, 745 genes were up-regulated in fruiting bodies, primarily associated with biosynthetic and catabolic pathways. Metabolomic analysis identified 1,161 metabolites, with 1,014 showing significant variations across developmental stages. Integrated transcriptomic and metabolomic analyses uncovered 17 genes positively correlated with 34 metabolites, which are likely crucial regulators of fruiting body development. These findings provide new insights into the molecular networks underlying C.blackwelliae growth and fruiting body formation.

Genomic and Methylomic Signatures Associated With the Maintenance of Genome Stability and Adaptive Evolution in Two Closely Allied Wolf Spiders

Pardosa spiders, belonging to the wolf spider family Lycosidae, play a vital role in maintaining the health of forest and agricultural ecosystems due to their function in pest control. This study presents chromosome-level genome assemblies for two allied Pardosa species, P. laura and P. agraria. Both species' genomes show a notable expansion of helitron transposable elements, which contributes to their large genome sizes. Methylome analysis indicates that P. laura has higher overall DNA methylation levels compared to P. agraria. DNA methylation may not only aids in transposable element-driven genome expansion but also positively affects the three-dimensional organisation of P. laura after transposon amplification, thereby potentially enhancing genome stability. Genes associated with hyper-differentially methylated regions in P. laura (compared to P. agraria) are enriched in functions related to mRNA processing and energy production. Furthermore, combined transcriptome and methylome profiling has uncovered a complex regulatory interplay between DNA methylation and gene expression, emphasising the important role of gene body methylation in the regulation of gene expression. Comparative genomic analysis shows a significant expansion of cuticle protein and detoxification-related gene families in P. laura, which may improve its adaptability to various habitats. This study provides essential genomic and methylomic insights, offering a deeper understanding of the relationship between transposable elements and genome stability, and illuminating the adaptive evolution and species differentiation among allied spiders.

Cadherin-16 regulates acoustic sensory gating in zebrafish through endocrine signaling

Sensory thresholds enable animals to regulate their behavioral responses to environmental threats. Despite the importance of sensory thresholds for animal behavior and human health, we do not yet have a full appreciation of the underlying molecular-genetic and circuit mechanisms. The larval zebrafish acoustic startle response provides a powerful system to identify molecular mechanisms underlying establishment of sensory thresholds and plasticity of thresholds through mechanisms like habituation. Using this system, we identify Cadherin-16 as a previously undescribed regulator of sensory gating. We demonstrate that Cadherin-16 regulates sensory thresholds via an endocrine organ, the corpuscle of Stannius (CS), which is essential in zebrafish for regulating Ca2+ homeostasis. We further show that Cadherin-16 regulates whole-body calcium and ultimately behavior through the hormone Stanniocalcin 1l (Stc1l), and the IGF-regulatory metalloprotease, Papp-aa. Finally, we demonstrate the importance of the CS through ablation experiments that reveal its role in promoting normal acoustic sensory gating. Together, our results uncover a previously undescribed brain non-autonomous pathway for the regulation of behavior and underscore Ca2+ homeostasis as a critical process underlying sensory gating in vivo.

Time-series RNA sequencing reveals that immulectin-1 (IML-1) is a potential target for biological control of fall armyworm, Spodoptera frugiperda

Insects mount robust immune responses to reduce the efficacy of biopesticides. Identification of critical immune targets for biocontrol strategies is thus essential. However, immune targets in the fall armyworm, Spodoptera frugiperda (Lepidoptera: Noctuidae), a major pest of corn, remain largely unexplored. Here, we investigated the immune response of S. frugiperda to the entomopathogen Beauveria bassiana (Hypocreales: Cordycipitaceae). Time-series RNA sequencing revealed 863, 717, and 165 differentially expressed genes (DEGs) at 24, 48, and 60 h post-infection, respectively. GO and KEGG enrichment analyses revealed that DEGs were significantly enriched in immune-related terms and pathways. Among the DEGs, the gene encoding immulectin SfIML-1, a C-type lectin with dual carbohydrate recognition domains, was significantly upregulated during the early and intermediate stages of infection. Genome-wide analysis identified nine IML genes, with SfIML-1 as the only one highly responsive to B. bassiana. Recombinant SfIML-1 binds pathogen-associated molecular patterns and pathogens, enhancing cellular responses. RNAi of SfIML-1 decreased survival after fungal infection, while injection of recombinant SfIML-1 increased survival. Suppression of the NF-κB and JAK-STAT pathways diminished SfIML-1 expression in both fat body and Sf9 cells. These results identify SfIML-1 as a key pattern recognition receptor and a promising target for developing novel biocontrol strategies.

The GlnE protein of Azorhizobium caulinodans ORS571 plays a crucial role in the nodulation process of the legume host Sesbania rostrata

The GlnE enzyme, functioning as an adenylyltransferase/adenylyl-removing enzyme, plays a crucial role in reversible adenylylation of glutamine synthetase (GS), which in turn regulates bacterial nitrogen assimilation. Genomic analysis of Azorhizobium caulinodans ORS571 revealed an open reading frame encoding a GlnE protein, whose function in the free-living and symbiotic states remains to be elucidated. A glnE deletion mutant retained high GS activity even under nitrogen-rich conditions. However, a reduction in growth was observed for the mutant strain at lower NH4+ concentrations than for the wild-type strain. Furthermore, the ΔglnE mutant strain showed reduced motility on ammonium-containing media. Inactivation of GlnE led to an increase in root adhesion, biofilm formation, and nodulation on Sesbania rostrata. Nevertheless, the nodules induced by the glnE mutant strain were ineffective. In addition, A. caulinodans GlnE played a significant role in enhancing resistance against environmental stresses, such as heat, heavy metals, and cumene hydroperoxide. This study demonstrates that GlnE plays multiple regulatory roles in A. caulinodans beyond nitrogen metabolism and is essential for establishing symbiotic relationships with host plants.

MibR and LibR are involved in the transcriptional regulation of the ipdC gene in Azospirillum brasilense Sp7

Azospirillum brasilense is a PGPR that produces the phytohormone IAA, a signaling molecule involved in bacteria-plant interaction processes. IAA biosynthesis in Azospirillum is mainly carried out via the IPyA pathway in which the enzyme phenylpyruvate decarboxylase encoded by the ipdC gene is the main. The promoter region of ipdC gene contains cis elements that are highly conserved among different Azospirillum strains. In this work, we identified two proteins that interact with the promoter region of the ipdC gene, named MibR and LibR that belong to the MarR and LuxR transcriptional regulators family, respectively. Both proteins have an HTH domain, and in the case of LibR, it has a REC domain, with aspartic acid residue conserved in positions 7, 8 and 54, this last as a possible phosphorylation target. To explore their participation in the regulation of the ipdC gene, mutants of libR, mibR, and libR-mibR double mutant were generated. The results showed a decrease in IAA biosynthesis that was related to the observed decrease in ipdC gene expression mostly in the doble mutant compared with the wild type. In this work we suggest that ipdC transcription is regulated by LibR and MibR, providing new findings insight into the mechanism employed by A. brasilense to control IAA biosynthesis.

β-Catenin localization in the ctenophore Mnemiopsis leidyi suggests an ancestral role in cell adhesion and nuclear function

BACKGROUND

The emergence of multicellularity in animals marks a pivotal evolutionary event, which was likely enabled by molecular innovations in the way cells adhere and communicate with one another. β-Catenin is significant to this transition due to its dual role as both a structural component in the cadherin-catenin complex and as a transcriptional coactivator involved in the Wnt/β-catenin signaling pathway. However, our knowledge of how this protein functions in ctenophores, one of the earliest diverging metazoans, is limited.

RESULTS

To study β-catenin function in the ctenophore Mnemiopsis leidyi, we generated affinity-purified polyclonal antibodies targeting Mlβ-catenin. We then used this tool to observe β-catenin protein localization in developing Mnemiopsis embryos. In this article, we provide evidence of consistent β-catenin protein enrichment at cell-cell interfaces in Mnemiopsis embryos. Additionally, we found β-catenin enrichment in some nuclei, particularly restricted to the oral pole around the time of gastrulation. The Mlβ-catenin affinity-purified antibodies now provide us with a powerful reagent to study the ancestral functions of β-catenin in cell adhesion and transcriptional regulation.

CONCLUSIONS

The localization pattern of embryonic Mlβ-catenin suggests that this protein had an ancestral role in cell adhesion and may have a nuclear function as well.

Primary regulatory T cell activator FOXP3 is present across Amphibia

The overall structure of the immune system is highly conserved across jawed vertebrates, but characterization and description of the immune system is heavily biased toward mammals. One arm of the vertebrate immune system, the adaptive immune system, mounts pathogen-specific responses that tend to be robust and effective at clearing pathogens. This system requires selection against self-recognition and modulation of the immune response. One of the mechanisms of immune modulation is the presence of regulatory T cells that suppress other effector immune cells. Regulatory T cells and their primary activator forkhead box protein P3 (FOXP3) have been well characterized in mammalian models but unexplored in most other vertebrate taxa. Amphibians are a good focal group for the characterization of FOXP3 due to their phylogenetic position on the vertebrate tree of life, and their susceptibility to emerging pathogens. In this study, we mined available transcriptomic and genomic data to confirm the presence of FOXP3 across the amphibian tree of life. We find that FOXP3 is present in all major clades of amphibians. We also test whether selection on FOXP3 shows signatures of intensification among the three main clades of amphibians, which may reflect shifts in the stringency of natural selection on this gene. Our findings provide insights into the evolutionary history of the vertebrate immune system and confirm the conservation of vertebrate immune genes within amphibians.

Evolution and expression of TLR5a and TLR5b in lamprey (Lampetra japonica)

The lamprey serves as a key model organism for studying the origin and evolution of species, embryonic development, and the immune system. The immune system primarily relies on pattern recognition receptors, including Toll-like receptors (TLRs), with Toll-like receptor 5 (TLR5) having a particularly complex evolutionary history. Currently, although TLR5 is being identified in an expanding array of taxonomic groups, a comprehensive study on its evolutionary aspects is yet to be conducted. In this study, we identified Lj-TLR5a and Lj-TLR5b in Lampetra japonica and examined their distribution across various tissues in this species. Furthermore, we conducted preliminary investigations into their immune functions and discovered that, as primitive genes, they are highly sensitive to various pathogens. Upon recognizing flagellar proteins, both Lj-TLR5a and Lj-TLR5b work together; however, these TLRs may function independently in response to other stimuli. Subsequently, we performed comprehensive structural and evolutionary analyses of the TLR5 family, incorporating TLR5 data from various species at different evolutionary stages. Our findings revealed that TLR5a and TLR5b in lamprey are in a relatively primitive evolutionary state. Meanwhile, TLR5L differentiated during the early stages of evolution and exhibits a pseudogenic trend throughout this evolutionary process; notably, this TLR is currently preserved only in certain amphibian and reptile species. In cartilaginous fish, only one type of TLR5 is usually retained, whereas bony fish typically possess both TLR5 and TLR5S. TLR5S has a relatively simple structure, likely arising from repetitive whole-genome events in bony fishes. Among bony fishes, flesh-finned fish were found to retain only one TLR5, which eventually evolved into the TLR5 found in quadrupeds. In summary, this study provides significant insights into the origin and evolution of the TLR5 family by analyzing the evolutionary status and immune functions of Lj-TLR5a and Lj-TLR5b in Japanese lampreys.

Ethylene signals through an ethylene receptor to modulate biofilm formation and root colonization in a beneficial plant-associated bacterium

Ethylene is a plant hormone involved in many aspects of plant growth and development as well as responses to stress. The role of ethylene in plant-microbe interactions has been explored from the perspective of plants. However, only a small number of studies have examined the role of ethylene in microbes. We demonstrated that Azospirillum brasilense contains a functional ethylene receptor that we call Azospirillum Ethylene Response1 (AzoEtr1) after the nomenclature used in plants. AzoEtr1 directly binds ethylene with high affinity. Treating cells with ethylene or disrupting the receptor reduces biofilm formation and colonization of plant root surfaces. Additionally, RNA sequencing and untargeted metabolomics showed that ethylene causes wide-spread metabolic changes that affect carbon and nitrogen metabolism. One result is the accumulation of poly-hydroxybutyrate. Our data suggests a model in which ethylene from host plants alters the density of colonization by A. brasilense and re-wires its metabolism, suggesting that the bacterium implements an adaptation program upon sensing ethylene. These data provide potential new targets to regulate beneficial plant-microbe interactions.

Protein phosphatase PP2C19 controls hypocotyl phototropism through the phosphorylation modification of NONPHOTOTROPIC HYPOCOTYL3 in Arabidopsis

Plants exhibit shoot growth in the direction of the light source to facilitate photosynthesis, known as positive phototropism. In Arabidopsis hypocotyl phototropism, it is thought that a gradient of the signal intensity of the blue light (BL) photoreceptor phototropin1 (phot1) between the light-irradiated and shaded sides leads to the differential growth of hypocotyls. The intensity of phot1 signal is regulated not only by the protein kinase activity of phot1 but also by the phosphorylation status of the NONPHOTOTROPIC HYPOCOTYL3 (NPH3) protein, which has a dark form and a BL form of the phosphorylation modification. Previous studies have shown that phot1 drives the forward reaction from the dark form to the BL form of NPH3. However, the molecular mechanism underlying the reverse reaction remains unknown. Here, we show that protein phosphatase PP2C19 controls the reverse reaction that converts the BL form of NPH3 to the dark form of NPH3. The PP2C19 protein possesses the protein phosphatase type 2C (PP2C) domain, two cyclic nucleoside monophosphate (cNMP)-binding domains, and the protein kinase domain. Similar to phot1 and NPH3, PP2C19 localizes to the plasma membrane, and its PP2C domain is necessary and sufficient for PP2C19 function in hypocotyl phototropism. The pp2c19 mutants show abnormalities in second positive hypocotyl phototropism with a delay in the reverse reaction of NPH3 phosphorylation modification. The present study suggests that continuous BL irradiation induces an equilibrium state of the reversible reaction of NPH3 phosphorylation, which acts as a phot1 signaling gradient with phot1 kinase activity to induce the second positive phototropism.

Computational analysis of congenital heart disease associated SNPs: unveiling their impact on the gene regulatory system

Congenital heart disease (CHD) is a prevalent condition characterized by defective heart development, causing premature death and stillbirths among infants. Genome-wide association studies (GWASs) have provided insights into the role of genetic variants in CHD pathogenesis through the identification of a comprehensive set of single-nucleotide polymorphisms (SNPs). Notably, 90-95% of these variants reside in the noncoding genome, complicating the understanding of their underlying mechanisms. Here, we developed a systematic computational pipeline for the identification and analysis of CHD-associated SNPs spanning both coding and noncoding regions of the genome. Initially, we curated a thorough dataset of SNPs from GWAS-catalog and ClinVar database and filtered them based on CHD-related traits. Subsequently, these CHD-SNPs were annotated and categorized into noncoding and coding regions based on their location. To study the functional implications of noncoding CHD-SNPs, we cross-validated them with enhancer-specific histone modification marks from developing human heart across 9 Carnegie stages and identified potential cardiac enhancers. This approach led to the identification of 2,056 CHD-associated putative enhancers (CHD-enhancers), 38.9% of them overlapping with known enhancers catalogued in human enhancer disease database. We identified heart-related transcription factor binding sites within these CHD-enhancers, offering insights into the impact of SNPs on TF binding. Conservation analysis further revealed that many of these CHD-enhancers were highly conserved across vertebrates, suggesting their evolutionary significance. Utilizing heart-specific expression quantitative trait loci data, we further identified a subset of 63 CHD-SNPs with regulatory potential distributed across various cardiac tissues. Concurrently, coding CHD-SNPs were represented as a protein interaction network and its subsequent binding energy analysis focused on a pair of proteins within this network, pinpointed a deleterious coding CHD-SNP, rs770030288, located in C2 domain of MYBPC3 protein. Overall, our findings demonstrate that SNPs have the potential to disrupt gene regulatory systems, either by affecting enhancer sequences or modulating protein-protein interactions, which can lead to abnormal developmental processes contributing to CHD pathogenesis.

Whole genome sequencing, assembly and annotation of the Southern Ground Hornbill - Bucorvus leadbeateri

The Southern Ground Hornbill (SGH - Bucorvus leadbeateri) is one of the largest hornbill species worldwide, known for its complex social structures and breeding behaviours. This bird has been of great interest due to its declining population and disappearance from historic ranges in southern Africa. Despite being the focus of numerous conservation efforts, with research forming an integral part of these initiatives, there is still a substantial lack of knowledge regarding the molecular biology aspects of this bird species. In this study, whole genome sequencing of the SGH was achieved using Illumina short-read (NovaSeq 6000) and Pacific Biosciences long-read technologies. A hybrid de novo genome assembly followed by reference-based refinement produced a 1.16 Gb high-quality draft genome assembly of the SGH comprised of 1,672 contigs (N50 value of 40.45 Mb). The availability of this genome will aid in improving our current understanding of this bird at a genomic level and provide a pivotal foundation for comparative genomic analysis with other hornbill species.

Computational Analysis of CC2D1A Missense Mutations: Insight into Protein Structure and Interaction Dynamics

CC2D1A is implicated in a range of conditions, including autism spectrum disorder, intellectual disability, seizures, autosomal recessive nonsyndromic intellectual disability, heterotaxy, and ciliary dysfunction. In order to understand the molecular mechanisms underlying these conditions, we focused on the structural and dynamic activity consequences of mutations within this gene. In this study, whole exome sequencing identified the c.1552G > A (GLU518LYS) missense mutation in the CC2D1A in an 18-year-old male, linking it to intellectual disability and autism. In addition to the GLU518LYS mutation, we conducted a comprehensive analysis of other predefined missense mutations (i.e., PRO192LEU, GLN506ARG, PRO532LEU, GLY781VAL, and GLY781GLU) found within the CC2D1A. Utilizing all-atom molecular dynamics (MD) simulations and neighborhood interaction analyses, we delve into the impact of these mutations on protein structure and function at an atomic level, aiming to shed light on their contribution to the pathogenesis of related diseases. The results suggest that GLU518LYS, GLY781VAL, and GLY781GLU mutations did not significantly alter overall global protein structure compared to the wild type, while PRO192LEU, GLN506ARG, and PRO532LEU exhibited slightly higher protein root-mean-square deviation (RMSD) values, which may indicate potential impacts on whole protein stability. Moreover, neighborhood interaction analysis indicated that ASP85 emerges as a unique interaction partner specifically associated with the GLU518LYS mutation, whereas LYS75, which interacts with the ASP85 in the mutated form, is absent in the wild type. This alteration signifies a crucial reconfiguration in the local interaction network at the site of the mutation.

Comprehensive Analysis of the NHX Gene Family and Its Regulation Under Salt and Drought Stress in Quinoa (Chenopodium quinoa Willd.)

Background/Objectives: Abiotic stresses such as salinity and drought significantly constrain crop cultivation and affect productivity. Quinoa (Chenopodium quinoa Willd.), a facultative halophyte, exhibits remarkable tolerance to drought and salinity stresses, making it a valued model for understanding stress adaptation mechanisms. The objective of this study was to identify and characterize Sodium/Hydrogen antiporter (NHX) genes from the quinoa genome and study their role in stress tolerance. Methods: We identified and characterized 10 NHX genes from the quinoa genome, which belong to the monovalent cation/proton antiporter 1 (CPA1) superfamily. Comprehensive analysis, including phylogenetic relationships, motif patterns, and structural characteristics, was performed to classify these genes into three subfamilies. Physicochemical properties such as isoelectric point (pI), GRAVY, and transmembrane domains were examined. Promoter analysis was conducted to identify cis-elements linked to abiotic stress responses, phytohormone signalling, and light regulation. qPCR analysis was used to assess the differential expression patterns of CqNHX genes under salt and drought stress. Results: The analysis revealed that the NHX genes were divided into three subfamilies localized to vacuolar, plasma, and endosomal membranes. These genes exhibited structural and functional diversity. Promoter analysis indicated the presence of cis-elements associated with abiotic stress responses, phytohormone signalling, and light regulation, suggesting diverse regulatory roles. qPCR analysis revealed differential expression patterns of CqNHX genes under salt and drought stress, with vacuolar NHXs showing higher induction in leaf tissues under salinity. This underscores their critical role in sodium sequestration and ion homeostasis. Evolutionary analysis indicated a high degree of conservation within subfamilies, alongside evidence of purifying selection. Conclusions: The findings enhance our understanding of the molecular basis of stress tolerance in quinoa and provide valuable targets for genetic engineering to improve crop resilience to environmental challenges.

Genome-wide identification of SAP family genes and characterization of TaSAP6-A1 to improve Cd tolerance in Triticum aestivum L

Stress Associated Proteins (SAPs) contain A20/AN1 zinc finger domains and, have been proposed to function in various physiological processes such as cold, salinity, drought, heavy metals, damage, and flooding resistance in plants. Here, a total of 131 SAP genes were identified, including T. aestivum (60), T. urartu (10), Ae. Tauschii (16), T. dicoccoides (13), O. sativa (18), and A. thaliana (14). A phylogenetic analysis revealed that the SAPs are clustered into two subfamilies. The TaSAP genes in the collinear region comprised 34 pairs of duplicated genes formed through segmental duplication events. Overexpressing TaSAP6-A1 in wheat enhanced Cd tolerance, whereas knock-down of this gene increased Cd sensitivity. Yeast two-hybrid (Y2H) and bimolecular fluorescent complementation assays (BiFC) demonstrated interaction between TaSAP6-A1 and phenylalanine ammonia-lyase (TaPAL), the first enzyme in the phenylpropanoid pathway. This study provides a valuable reference for further investigations into the functional and molecular mechanisms of the SAP gene family.

Genetic variation mining of the Chinese mitten crab (Eriocheir sinensis) based on transcriptome data from public databases

At present, public databases house an extensive repository of transcriptome data, with the volume continuing to grow at an accelerated pace. Utilizing these data effectively is a shared interest within the scientific community. In this study, we introduced a novel strategy that harnesses SNPs and InDels identified from transcriptome data, combined with sample metadata from databases, to effectively screen for molecular markers correlated with traits. We utilized 228 transcriptome datasets of Eriocheir sinensis from the NCBI database and employed the Genome Analysis Toolkit software to identify 96 388 SNPs and 20 645 InDels. Employing the genome-wide association study analysis, in conjunction with the gender information from databases, we identified 3456 sex-biased SNPs and 639 sex-biased InDels. The KOG and KEGG annotations of the sex-biased SNPs and InDels revealed that these genes were primarily involved in the metabolic processes of E. sinensis. Combined with SnpEff annotation and PCR experimental validation, a highly sex-biased SNP located in the Kelch domain containing 4 (Klhdc4) gene, CHR67-6415071, was found to alter the splicing sites of Klhdc4, generating two splice variants, Klhdc4_a and Klhdc4_b. Additionally, Klhdc4 exhibited robust expression across the ovaries, testes, and accessory glands. The sex-biased SNPs and InDels identified in this study are conducive to the development of unisexual cultivation methods for E. sinensis, and the alternative splicing event caused by the sex-biased SNP in Klhdc4 may serve as a potential mechanism for sex regulation in E. sinensis. The analysis strategy employed in this study represents a new direction for the rational exploitation and utilization of transcriptome data in public databases.

NOVEL INVENTION OF SPORE INDUCTION IN A SISTER SPECIES TO GROUP 4 DICTYOSTELIA

Background

Dictyostelia are soil amoebas that aggregate to form fruiting bodies with spores and stalk cells in response to starvation. Where known, species across the dictyostelid phylogeny use secreted cAMP, detected by cAMP receptors (cARs) to induce the differentiation of spores and to organize fruiting body construction. However, recent deletion of the single cAR of Polyspondylium violaceum (Pvio) left both its fruiting bodies and spores intact.

Methods

To investigate whether Pvio sporulation can occur in the absence of secreted cAMP and to explore alternative inducers in a bioassay , three prespore genes were identified and gene fusions of their promoters with the LacZ reporter gene were transformed into Pvio cells. After assessing the spatial expression pattern of the genes and the stage at which prespore gene expression initiated, the effect of cAMP and other Dictyostelium discoideum ( Ddis) signal molecules were tested on prespore gene expression in vitro.

Results

Pvio genes g4562 (psp1), g2696 (psp2) and g2380 (psp3) were identified as homologs of Ddis spore coat genes. They were first expressed around 4 h of starvation in aggregation centres and later in the posterior 4/5 th of emerging sorogens and the spore head of early fruiting bodies. Cells from dissociated 4 h aggregates and shaken in suspension for 6 h increased prespore- LacZ reporter activity 4-fold for psp1 and 6-fold for psp2, but this increase was at least 5-fold higher when cells were plated on solid substratum for 6 h to develop normally. cAMP had no effect on prespore gene induction and neither had the Pvio chemoattractant glorin nor the Ddis chemoattractants and differentiation inducers folate, c-di-GMP, DIF-1, GABA, cGMP and 8Br-cAMP.

Conclusions

The Pvio lineage uniquely evolved a novel genetic network for synthesis, detection and processing of the signal that triggers its main survival strategy.

The extracellular segment of CroS is not required for sensing but fine-tunes the magnitude of CroS signaling to regulate cephalosporin resistance in Enterococcus faecalis

Enterococci are Gram-positive bacteria that colonize the gastrointestinal tract. Clinically relevant enterococci are intrinsically resistant to antibiotics in the cephalosporin family, and prior therapy with cephalosporins is a major risk factor for the acquisition of an enterococcal infection. One important determinant of intrinsic cephalosporin resistance in enterococci is the two-component signal transduction system CroS/R. The CroS sensor kinase senses cephalosporin-induced cell wall stress to become activated and phosphorylates its cognate response regulator CroR, thereby enhancing CroR-dependent gene expression to drive cephalosporin resistance. CroS possesses a short (~30 amino acids) extracellular segment between its two transmembrane domains near the N-terminus, but whether this extracellular segment is important for sensing cephalosporin stress, or possesses any other function, has remained unknown. Here, we explored the role of the CroS extracellular segment through mutagenesis and functional studies. We found that mutations in the CroS extracellular segment biased CroS to adopt a more active state during ceftriaxone stress, which led to an increase in CroR-dependent gene expression and hyper-resistance to ceftriaxone. Importantly, these mutants still responded to ceftriaxone-mediated stress by enhancing CroS activity, indicating that the extracellular segment of CroS does not directly bind a regulatory ligand. Overall, our results suggest that although the extracellular segment of CroS does not directly bind a regulatory ligand, it can modulate the magnitude of CroS signaling for phosphorylation of CroR to regulate cephalosporin resistance through the resulting changes in CroR-dependent gene expression.

IMPORTANCE

Clinically relevant enterococci are intrinsically resistant to antibiotics in the cephalosporin family. The CroS sensor kinase senses cephalosporin-induced cell wall stress to trigger signaling that drives cephalosporin resistance, but the mechanism by which CroS senses stress is unknown. We report the first functional characterization of the CroS extracellular segment, revealing that mutations in the extracellular segment did not prevent CroS from responding to cell wall stress but instead biased CroS to adopt a more active state during cephalosporin stress that led to an increase in CroR-dependent gene expression and hyper-resistance to ceftriaxone. Overall, our results suggest that the extracellular segment of CroS does not directly bind to a regulatory ligand but that it can modulate the magnitude of CroS signaling.

Discovery and characterization of a novel Lepidoptera-specific antimicrobial peptide from the fall armyworm, Spodoptera frugiperda (Lepidoptera: Noctuidae)

Antimicrobial peptides (AMPs) are critical components of innate immunity in diverse organisms, including plants, vertebrates, and insects. This study identified and characterized a novel Lepidoptera-specific AMP, named lepidoptin, from the invasive pest Spodoptera frugiperda (Lepidoptera: Noctuidae). Lepidoptin is a 116-amino acid protein containing a signal peptide and a novel β-sandwich domain that is distinct from previously reported AMPs. Temporal and spatial expression analyses revealed a significant upregulation of the lepidoptin gene in vivo and in cultured SF9 cells in response to pathogens. Molecular docking analysis identified a specific binding cavity. Enzyme-linked immunosorbent assay and binding assays confirmed that lepidoptin can bind to pathogen-associated molecular patterns, bacteria, and fungi. Recombinant lepidoptin exhibited potent antibacterial activity by inducing bacterial agglutination, inhibiting bacterial growth, increasing bacterial membrane permeability, and preventing biofilm formation. Lepidoptin also showed antifungal activity against the entomopathogenic fungus Beauveria bassiana by inhibiting spore germination, increasing fungal cell permeability, and increasing reactive oxygen species. Injection of recombinant lepidoptin into S. frugiperda larvae increased survival after B. bassiana infection, whereas knockdown of lepidoptin by RNA interference decreased larval survival. In addition, lepidoptin showed antimicrobial activity against the plant pathogen Fusarium graminearum by inhibiting spore germination and alleviating disease symptoms in wheat seedlings and cherry tomatoes. This study demonstrates the remarkable dual functionality of lepidoptin in enhancing S. frugiperda immunity and controlling plant pathogens, making it a promising candidate for biocontrol strategies in both pest management and plant disease prevention.

Genome-wide identification and analysis of ERF transcription factors related to abiotic stress responses in Nelumbo nucifera

BACKGROUND

Ethylene-responsive factor (ERF) transcription factors belong to the APETALA2/ERF (AP2/ERF) superfamily, and play crucial roles in plant development process and stress responses. However, the function of ERF proteins (especially for their role in response to abiotic stresses) remains scarce in Nelumbo nucifera, which is an important aquatic plant with high ornamental, economic, and ecological values.

RESULTS

A total of 107 ERF genes were identified from the N. nucifera genome, and phylogenetic analysis classified these genes into 11 groups. The NnERF genes in the same group exhibited similar gene structure and conserved motifs, and they were unevenly distributed across the 8 chromosomes, with three pairs of tandem duplications and 21 pairs of segmental duplications. Synteny analysis revealed 44 and 39 of NnERF genes were orthologous to those in Arabidopsis thaliana and Oryza sativa, respectively. Tissue-specific expression patterns analysis of NnERF showed that 26 NnERF genes were expressed in all tested tissues, in which five genes exhibited high expression levels. Furthermore, 16 NnERF genes were selected for exploring their responses to different abiotic stresses, including cold, salt, drought, and Cd stresses. qRT-PCR analysis revealed that all these 16 investigated genes were regulated by at least one stress treatment, and 12 genes responded to all the stress treatments with different expression patterns or levels, suggesting their potential roles in diverse abiotic stress tolerance of N. nucifera. Additionally, two representative stress-related NnERFs (Nn3g19628 and Nn1g06033) were confirmed to be nuclear-localized proteins and displayed transcriptional activation.

CONCLUSIONS

In this study, we conducted a genome-wide identification and analysis of NnERF gene family related to abiotic stress responses in N. nucifera, which provides valuable information for further functional validation of these genes in stress responses, and forms a foundation for stress tolerance breeding in N. nucifera and other aquatic ornamental plants.

Establishment of the auxin inducible degron system for Babesia duncani: a conditional knockdown tool to study precise protein regulation in Babesia spp

BACKGROUND

Babesia duncani is a pathogen within the phylum Apicomplexa that causes human babesiosis. It poses a significant threat to public health, as it can be transmitted not only through tick bites but also via blood transfusion. Consequently, an understanding of the gene functions of this pathogen is necessary for the development of drugs and vaccines. However, the absence of conditional gene knockdown tools has hindered the research on this pathogen. The auxin-inducible degron (AID) system is a rapid, reversible conditional knockdown system widely used in gene function studies. Thus, there is an urgent need to establish the AID system in B. duncani to study essential gene functions.

METHODS

The endogenous genes of the Skp1-Cullin-F-box (SCF) complex in B. duncani were identified and confirmed through multiple sequence alignment and conserved domain analysis. The expression of the F-box protein TIR1 from Oryza sativa (OsTIR1) was achieved by constructing a transgenic parasite strain using a homologous recombination strategy. Polymerase chain reaction (PCR), western blot, and indirect immunofluorescence assay (IFA) were used to confirm the correct monoclonal parasite strain. The degradation of enhanced green fluorescent protein (eGFP) tagged with an AID degron was detected through western blot and live-cell fluorescence microscopy after treatment of indole-3-acetic acid (IAA).

RESULTS

In this study, Skp1, Cul1, and Rbx1 of the SCF complex in B. duncani were identified through sequence alignment and domain analysis. A pure BdTIR1 strain with expression of the OsTIR1 gene was constructed through homologous recombination and confirmed. This strain showed no significant differences from the wild type (WT) in terms of growth rate and proportions of different parasite forms. The eGFP tagged with an AID degron was successfully induced for degradation using 500 μM IAA. Grayscale analysis of western blot indicated a 61.3% reduction in eGFP expression levels, while fluorescence intensity analysis showed a 77.5% decrease in fluorescence intensity. Increasing the IAA concentration to 2 mM accelerated eGFP degradation and enhanced the extent of degradation.

CONCLUSIONS

This study demonstrated the functionality of the AID system in regulating protein levels by inducing rapid degradation of eGFP using IAA, providing an important research tool for studying essential gene functions related to invasion, egress, and virulence of B. duncani. Moreover, it also offers a construction strategy for apicomplexan parasites that have not developed an AID system.

Exploring the naturally acquired response to Pvs47 gametocyte antigen

Malaria represents a challenging global public health task, with Plasmodium vivax being the predominant parasite in Brazil and the most widely distributed species throughout the world. Developing a vaccine against P. vivax malaria demands innovative strategies, and targeting gametocyte antigens shows promise for blocking transmission prevention. Among these antigens, Pvs47, expressed in gametocytes, has shown remarkable efficacy in transmission blocking. However, remains underexplored in vaccine formulations. This study employed in silico methods to comprehensively characterize the physicochemical properties, structural attributes, epitope presence, and conservation profile of Pvs47. Additionally, we assessed its antigenicity in individuals exposed to malaria in endemic Brazilian regions. Recombinant protein expression occurred in a eukaryotic system, and antigenicity was evaluated using immunoenzymatic assays. The responses of naturally acquired IgM, total IgG, and IgG subclasses were analyzed in three groups of samples from Amazon region. Notably, all samples exhibited anti-Pvs47 IgM and IgG antibodies, with IgG3 predominating. Asymptomatic patients demonstrated stronger IgG responses and more diverse subclass responses. Anti-Pvs47 IgM and IgG responses in symptomatic individuals decrease over time. Furthermore, we observed a negative correlation between anti-Pvs47 IgM response and gametocytemia in samples of symptomatic patients, indicating a gametocyte-specific response. Additionally, negative correlation was observed among anti-Pvs47 antibody response and hematocrit levels. Furthermore, comparative analysis with widely characterized blood antigens, PvAMA1 and PvMSP119, revealed that Pvs47 was equally or more recognized than both proteins. In addition, there is positive correlation between P. vivax blood asexual and sexual stage immune responses. In summary, our study unveils a significant prevalence of anti-Pvs47 antibodies in diverse Amazonian samples and the importance of IgM response for gametocytes depuration. These findings regarding the in silico characterization and antigenicity of Pvs47 provide crucial insights for potential integration into P. vivax vaccine formulations.

AnnoDUF: A Web-Based Tool for Annotating Functions of Proteins Having Domains of Unknown Function

The rapid expansion of biological sequence databases due to high-throughput genomic and proteomic sequencing methods has left a considerable number of identified protein sequences with unclear or incomplete functional annotations. Domains of unknown function (DUFs) are protein domains that lack functional annotations but are present in numerous proteins. To address the challenge of finding functional annotations for DUFs, we have developed a computational method that efficiently identifies and annotates these enigmatic protein domains by utilizing the position-specific iterative basic local alignment search tool (PSI-BLAST) and data mining techniques. Our pipeline identifies putative potential functionalities of DUFs, thereby decreasing the gap between known sequences and functions. The tool can also take user input sequences to annotate. We executed our pipeline on 5111 unique DUF sequences obtained from Pfam, resulting in putative annotations for 2007 of these. These annotations were subsequently incorporated into a comprehensive database and interfaced with a web-based server named "AnnoDUF". AnnoDUF is freely accessible to both academic and industrial users, via the World Wide Web at the link http://bts.ibab.ac.in/annoduf.php. All scripts used in this study are uploaded to the GitHub repository, and these can be accessed from https://github.com/BioToolSuite/AnnoDUF.

Identification of significant SNPs and candidate loci for blast disease resistance via GWAS and population structure analysis in ARC panel of Oryza sativa

Pyricularia (syn. Magnaporthe) oryzae is responsible for the blast disease in rice resulting in a greater extent of yield loss. However, some of the cultivars of rice have the ability to survive this devastating infection due to the presence of R (resistance) genes. Therefore, genome wide association study (GWAS) was undertaken using a panel of 400 rice landraces (ARC panel) and a set of filtered 38,723 single nucleotide polymorphisms (SNPs). The highest SNPs were mapped to chromosome 1 with a number of 4332 SNPs and lowest (2252) in chromosome 12. The ARC panel was evaluated phenotypically which revealed that 6% of the selected cultivars has resistance to rice blast disease with SES score of 1. The majority of the resistant cultivars belong to the group Asra of the panel. The population structure analysis was executed wherein three genetic subpopulations were identified namely RC1, RC2, RC3 and an admixture population constituting 48 accessions. Further, GWAS detected 15 significant association signal with P value in the range of 1.03E-05 to 1.03E-04, effect ranged from - 1.18 to 1.06, phenotypic variance explained was from 0 to 7.14%, R2 of 0.047 to 0.058, and minor allele frequency of 0.107 to 0.444. Eleven (Os01g39980, Os01g56130, Os01g67100, Os01g67110, Os03g41030, Os04g33310, Os07g42104, Os09g06464, Os09g08920, Os09g38800, Os12g37680) out of these 15 significant associations were identified as the candidate loci for the blast resistance in rice that will serve as an important genetic resistance source to be introgressed into an elite rice line in future breeding programs for deciphering blast resistance in rice. The GWAS study presented in this article helped to uncover significant gene regions which encode proteins to resist blast infection in rice plant. This is the first report on the GWAS analysis for blast resistance in unique landraces of rice from Northeast India employing single nucleotide polymorphism.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01518-6.

Strenuous expression of porcine epidemic diarrhea virus ORF3 protein suggests host resistance

Porcine epidemic diarrhea virus is attenuated upon adaptation to cell culture. Exclusively genomic mutations have been traced to the ORF3 gene of the laboratory strains. Previous attempts to express the protein were unsuccessful. We sought to express the ORF3 protein in both mammalian and bacteria cells as a prerequisite for investigation of the protein's role. For prokaryotic expression, two vector systems, pET28-a(+) and pGEX-4T-1 were constructed and expressed in Escherichia coli cells. For eukaryotic analyses, ORF3/pEGFP-C1 vector constructs were expressed in human embryonic, green monkey kidney and mouse fibrous cells. Intriguingly, there was minimal expression of the ORF3 gene. Following a documented hint that truncated ORF3 revealed higher expression, ORF3 gene was truncated. The simple modular architecture research tool analysis predicted two transmembrane domains between amino acid (aa) 41-63 and aa 76-98. Consequently, we generated two fragments; ORF-N (aa 1-98) inclusive of transmembrane domains and ORF3-C (aa 99-224). These truncated sequences were constructed as the whole gene here referred to as ORF3 wild type (wt). Coomassie blue stained gels revealed bands of ORF3-C expressed as a fusion protein of 17.5 and 39 kDa in pET28-a(+) and pGEX-4T-1 vectors, respectively. In contrast, ORF3-N was not. Additionally, ORF3-N induction decreased total cellular proteins suggesting inhibition of protein synthesis or metabolism. Solubility tests carried out at 30 °C, 25 °C and 18 °C showed that ORF3 formed inclusion bodies. Similar findings were observed in mammalian cells. Noteworthy, morphological distortions appeared in mammalian cells expressing ORF3 protein or its truncated mutants suggesting significance in host viability.

Genomic insights into cytokinin oxidase/dehydrogenase (CKX) gene family, identification, phylogeny and synteny analysis for its possible role in regulating seed number in Pigeonpea (Cajanus cajan (L.) Millsp.)

Cytokinin oxidase/dehydrogenase (CKX) regulates cytokinin levels in plants which are vital for plant growth and development. However, there is a paucity of evidence regarding their role in controlling embryo/seed development in pigeonpea. This comprehensive study provides information on the identification and characterization of CKX genes in pigeonpea. A genome-wide analysis identified 18 CKX genes, each with distinct structure, expression patterns, and possible diverse functions. Domain analysis revealed the presence of the sequences including FAD and CK-Binding domain, and subcellular localization analysis showed that almost 50 % of them reside within the nucleus. They were observed to be located unevenly on chromosome numbers 2, 4, 6, 7, and 11 with a majority of them present on the scaffolds. The 8 homologous pairs and various orthologous gene pairs provided further insights into their evolution pattern. Further, SNP/Indels variation in CKX genes and haplotype groups among contrasting genotypes for SNPP (seed number per pod) were analyzed. Spatiotemporal expression analysis revealed the significant expression pattern of CcCKX15, CcCKX17, and CcCKX2 in genotypes carrying low SNPP reiterating their possible role as negative regulators. These genes can be potential targets to undertake seed and biomass improvement in pigeonpea.

Bovine Transcription Factor POU Class 2 Homeobox 1 (POU2F1/Oct1) Protein Promotes BoHV-1 Replication in MDBK Cells

Bovine herpesvirus type 1 (BoHV-1) causes severe diseases in bovine species and great economic burden to the cattle industry worldwide. Due to its complex life cycle, many host factors that affect BoHV-1 replication remain to be explored. To understand the possible roles that the Oct1 cellular protein could play in this process, we first created Oct1-deficient MDBK cells using CRISPR/Cas9-mediated genome editing. Upon infection, the absence of Oct1 in MDBK cells significantly impacted BoHV-1 replication, a phenotype rescued by over-expressing the wild-type Oct1 protein in the deficient cells. We further found that the expression of all three classes of temporal genes, including essential and non-essential viral genes, were significantly reduced in Oct1 knockout MDBK cells, following both high and low multiplicity of infection. In summary, our findings confirm that the bovine Oct1 protein acts as a pro-viral factor for BoHV-1 replication by promoting its viral gene transcription in MDBK cells.

Unusual photodynamic characteristics of the light-oxygen-voltage domain of phototropin linked to terrestrial adaptation of Klebsormidium nitens

Phototropin (Phot), a blue light-sensing LOV domain protein, mediates blue light responses and is evolutionarily conserved across the green lineage. Klebsormidium nitens, a green terrestrial alga, presents a valuable opportunity to study adaptive responses from aquatic to land habitat transitions. We determined the crystal structure of Klebsormidium nitens Phot LOV1 domain (KnLOV1) in the dark and engineered different mutations (R60K, Q122N, and D33N) to modulate the lifetime of the photorecovery cycle. We observed unusual, slow recovery kinetics in the wild-type KnLOV1 domain (τ = 41 ± 3 min) compared to different mutants (R60K: τ = 2.0 ± 0.1 min, Q122N: τ = 1.7 ± 0.1 min, D33N: τ = 9.6 ± 0.1 min). Crystal structures of wild-type KnLOV1 and mutants revealed subtle but critical changes near the protein chromophore that is responsible for modulating protein dark recovery time. Our findings shed light on the unique structural and biochemical characteristics of the newly studied KnLOV1 and its evolutionary importance for phototropin-mediated physiology.

Cadherin 16 promotes sensory gating via the endocrine corpuscles of Stannius

Sensory thresholds enable animals to regulate their behavioral responses to environmental threats. Despite the importance of sensory thresholds for animal behavior and human health, we do not yet have a full appreciation of the underlying molecular-genetic and circuit mechanisms. The larval zebrafish acoustic startle response provides a powerful system to identify molecular mechanisms underlying establishment of sensory thresholds and plasticity of thresholds through mechanisms like habituation. Using this system, we identify Cadherin 16 as a previously undescribed regulator of sensory gating. We demonstrate that Cadherin 16 regulates sensory thresholds via an endocrine organ, the corpuscle of Stannius (CS), which is essential in zebrafish for regulating Ca2+ homeostasis. We further show that Cadherin 16 regulates whole-body calcium and ultimately behavior through the hormone Stanniocalcin 1L, and the IGF-regulatory metalloprotease, Papp-aa. Finally, we demonstrate the importance of the CS through ablation experiments that reveal its role in promoting normal acoustic sensory gating. Together, our results uncover a previously undescribed brain non-autonomous pathway for the regulation of behavior and establish Ca2+ homeostasis as a critical process underlying sensory gating in vivo.

Biomolecular Condensation of Trypsin Prevents Autolysis and Promotes Ca2+-Mediated Activation of Esterase Activity

The presence of Ca2+ ions is known to facilitate the activity of trypsin-like serine proteases via structural stabilization against thermal denaturation and autolysis. Herein, we report a new and hidden regulatory role of Ca2+ in the catalytic pathways of trypsin and α-chymotrypsin under physiological conditions. We discovered that macromolecular crowding promotes spontaneous homotypic condensation of trypsin via liquid-liquid phase separation to yield membraneless condensates over a broad range of concentrations, pH, and temperature, which are stabilized by multivalent hydrophobic interactions. Interestingly, we found that Ca2+ binding in the calcium binding loop reversibly regulates the condensation of trypsin and α-chymotrypsin. Spontaneous condensation effectively prevents autolysis of trypsin and preserves its native-like esterase activity for a prolonged period of time. It has also been found that phase-separated trypsin responds to Ca2+-dependent activation of its esterase activity even after 14 days of storage while free trypsin failed to do so. The present study highlights an important physiological aspect by which cells can spatiotemporally regulate the biocatalytic efficacy of trypsin-like serine proteases via Ca2+-signaling.

Functional complementation of two splicing variants of Gustavus in Neocaridina denticulata sinensis during ovarian maturation

Gustavus, a positive regulator in arthropod reproduction, features a conserved SPRY and a C-terminal SOCS box domain and belongs to the SPSB protein family. The SPSB family, encompassing SPSB1 to SPSB4, plays pivotal roles in higher animals, including immune response, apoptosis, growth, and stress responses. In Neocaridina denticulata sinensis, alternative splicing yielded two NdGustavus isoforms, NdGusX1 and NdGusX2, with distinct expression patterns-high in ovaries and muscles, respectively, and across all ovarian germ cells. These isoforms showed similar expression dynamics during embryogenesis and significant upregulation post-copper ion exposure (P < 0.05). The in situ hybridization result elucidated that NdGusX1 and NdGusX2 were expressed across the germ cell spectrum in the ovary, with NdGusX1 showing enhanced expression in oogonia and primary oocytes. In addition, RNA interference revealed functional complementation in ovaries and potential functional differentiation in muscles. Knockdown of NdGusX1 and NdGusX2 potentially disrupted endogenous vitellogenin synthesis, regulating vitellogenesis and reducing mature oocyte volume, affecting follicular cavity occupation. This study provides a theoretical framework for understanding the biological functions of the SPSB family in crustacean ovarian maturation.

The growth-blocking peptide is a dual regulator of development and immunity with biocontrol potential in Spodoptera frugiperda (Lepidoptera: Noctuidae)

Insect growth-blocking peptides (GBPs) are a family of cytokines found in several insect orders and are known for their roles in regulating development, paralysis, cell proliferation, and immune responses. Despite their diverse functions, the potential of GBPs as biocontrol targets against the pest Spodoptera frugiperda (Lepidoptera: Noctuidae) has not been fully explored. In this study, S. frugiperda GBP (SfGBP) was identified and functionally characterized. SfGBP is synthesized as a 146 amino acid proprotein with a 24 amino acid C-terminal active peptide (Glu123-Gly146). Predominant expression of SfGBP occurs in fourth to sixth instar larvae and in the larval fat body, with significant upregulation in response to pathogens and pathogen-associated molecular patterns. Injection of the synthetic active peptide into larvae induced growth retardation, delayed pupation, and increased survival against Beauveria bassiana infection. Conversely, RNA interference-mediated knockdown of SfGBP resulted in accelerated growth, earlier pupation, and decreased survival against B. bassiana infection. Further analysis revealed that SfGBP promoted SF9 cell proliferation and spreading, enhanced bacteriostatic activity of larval hemolymph, and directly inhibited germination of B. bassiana conidia. In addition, SfGBP enhanced humoral responses, such as upregulation of immunity-related genes and generation of reactive oxygen species, and cellular responses, such as nodulation, phagocytosis, and encapsulation. These results highlight the dual regulatory role of SfGBP in development and immune responses and establish it as a promising biocontrol target for the management of S. frugiperda.

Stress pathway outputs are encoded by pH-dependent clustering of kinase components

Signal processing by intracellular kinases controls near all biological processes but how signal pathway functions evolve with changed cellular context is poorly understood. Functional specificity of c-Jun N-terminal Kinases (JNK) are partly encoded by signal strength. Here we reveal that intracellular pH (pHi) is a significant component of the JNK network and defines signal response to specific stimuli. We show pHi regulates JNK activity in response to cell stress, with the relationship between pHi and JNK activity dependent on specific stimuli and upstream kinases activated. Using the optogenetic clustering tag CRY2, we show that an increase in pHi promotes the light-induced phase transition of ASK1 to augment JNK activation. While increased pHi similarly promoted CRY2-tagged JNK2 to form light-induced condensates, this attenuated JNK activity. Mathematical modelling of feedback signalling incorporating pHi and differential contributions by ASK1 and JNK2 condensates was sufficient to delineate signal responses to specific stimuli. Taking pHi and ASK1/JNK2 signal contributions into consideration may delineate oncogenic versus tumour suppressive JNK functions and cancer cell drug responses.

ALA Promotes Sucrose Accumulation in Early Peach Fruit by Regulating SPS Activity

5-Aminolevulinic acid (ALA), as a novel plant growth regulator, is a critical precursor for the biosynthesis of porphyrin compounds in all organisms. Many studies have reported that exogenous ALA treatment could improve fruit sweetness. However, the mechanism by which ALA promotes the increase in sugar content in fruit remains unclear. In this study, we found that ALA significantly promoted sucrose accumulation and SPS (sucrose phosphate synthase) activity in peach fruit. At 14, 28, 42, 50 and 60 days after ALA treatment, sucrose content of fruit was increased by 23%, 43%, 37%, 40% and 16%, respectively, compared with control treatment, and SPS enzyme activity was increased by 21%, 28%, 47%, 37% and 29%, respectively. Correlation analysis showed that the sucrose content of peach fruit under ALA treatment was significantly positively correlated with SPS activity. Subsequently, bioinformatics was used to identify SPS gene family members in peach fruit, and it was found that there were four members of the PpSPS gene family, distributed on chromosomes 1, 7 and 8, named PpSPS1, PpSPS2, PpSPS3 and PpSPS4, respectively. The results of qRT-PCR showed that PpSPS2 and PpSPS3 were highly expressed in response to ALA during fruit development, and the expression of PpSPS2 was positively correlated with SPS activity and sucrose accumulation in peach fruit. The results of tobacco subcellular localization showed that PpSPS2 was mainly distributed in the cytoplasm and nucleus, while PpSPS3 was mainly distributed in the nucleus. The results of this study will lay the foundation for further study on the functions of PpSPS and the regulation of sugar metabolism during the development and ripening of peach fruit by ALA.

CANDy: Automated analysis of domain architectures in carbohydrate-active enzymes

Carbohydrate-active enzymes (CAZymes) can be found in all domains of life and play a crucial role in metabolic and physiological processes. CAZymes often possess a modular structure, comprising not only catalytic domains but also associated domains such as carbohydrate-binding modules (CBMs) and linker domains. By exploring the modular diversity of CAZy families, catalysts with novel properties can be discovered and further insight in their biological functions and evolutionary relationships can be obtained. Here we present the carbohydrate-active enzyme domain analysis tool (CANDy), an assembly of several novel scripts, tools and databases that allows users to analyze the domain architecture of all protein sequences in a given CAZy family. CANDy's usability is shown on glycoside hydrolase family 48, a small yet underexplored family containing multi-domain enzymes. Our analysis reveals the existence of 35 distinct domain assemblies, including eight known architectures, with the remaining assemblies awaiting characterization. Moreover, we substantiate the occurrence of horizontal gene transfer from prokaryotes to insect orthologs and provide evidence for the subsequent removal of auxiliary domains, likely through a gene fission event. CANDy is available at https://github.com/PyEED/CANDy.

A rho-type GTPase activating protein affects the growth and development of Cordyceps cicadae

Cordyceps cicadae is recognized for its medicinal properties, attributed to bioactive constituents like polysaccharides and adenosine, which have been shown to improve kidney and liver functions and possess anti-tumor properties. Rho GTPase activating proteins (Rho GAPs) serve as inhibitory regulators of Rho GTPases in eukaryotic cells by accelerating the GTP hydrolysis of Rho GTPases, leading to their inactivation. In this study, we explored the function of the CcRga8 gene in C. cicadae, which encodes a Rho-type GTPase activating protein. Our study found that the knockout of CcRga8 resulted in a decrease in polysaccharide levels and an increase in adenosine concentration. Furthermore, the mutants exhibited altered spore yield and morphology, fruiting body development, decreased infectivity, reduced resistance to hyperosmotic stress, oxidative conditions, and cell wall inhibitors. These findings suggest that CcRga8 plays a crucial role in the development, stress response, and bioactive compound production of C. cicadae.

Genome-scale identification, expression and evolution analysis of B-box members in Dendrobium huoshanense

B-box (BBX) proteins have been recognized as vital determinants in plant development, morphogenesis, and adaptive responses to a myriad of environmental stresses. These zinc-finger proteins play a pivotal role in various biological processes. Their influence spans photomorphogenesis, the regulation of flowering, and imparting resilience to a wide array of challenges, encompassing both biotic and abiotic factors. Chromosome localization, gene structure and conserved motifs, phylogenetic analysis, collinearity analysis, expression profiling, fluorescence quantitative analysis, and tobacco transient transformation methods were used for functional localization and expression pattern analysis of the DhBBX gene. A total of 23 DhBBX members were identified from Dendrobium huoshanense. Subsequent phylogenetic evaluations effectively segregated these genes into five discrete evolutionary subsets. The predictions of subcellular localizations revealed that all these proteins were localized in the nucleus. The genetic composition and patterns showed that the majority of these genes consisted of several exons, with a few variations that could be attributed to transposon insertion. A comprehensive analysis using qRT-PCR was conducted to unravel the expression patterns of these genes in D. huoshanense, with a specific concentration on their responses to various hormone treatments and cold stress. Subcellular localization reveals that DhBBX21 and DhBBX9 are located in the nucleus. Our results provide a deep comprehension of the complex regulatory mechanisms of BBXs in response to various environmental and hormonal stimuli. These discoveries encourage further detailed and focused investigations into the operational dynamics of the BBX gene family in a wider range of plant species.

Genome-wide identification and mining elite allele variation of the Monoacylglycerol lipase (MAGL) gene family in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Monoacylglycerol lipase (MAGL) genes belong to the alpha/beta hydrolase superfamily, catalyze the terminal step of triglyceride (TAG) hydrolysis, converting monoacylglycerol (MAG) into free fatty acids and glycerol.

RESULTS

In this study, 30 MAGL genes in upland cotton have been identified, which have been classified into eight subgroups. The duplication of GhMAGL genes in upland cotton was predominantly influenced by segmental duplication events, as revealed through synteny analysis. Furthermore, all GhMAGL genes were found to contain light-responsive elements. Through comprehensive association and haplotype analyses using resequencing data from 355 cotton accessions, GhMAGL3 and GhMAGL6 were detected as key genes related to lipid hydrolysis processes, suggesting a negative regulatory effect.

CONCLUSIONS

In summary, MAGL has never been studied in upland cotton previously. This study provides the genetic mechanism foundation for the discover of new genes involved in lipid metabolism to improve cottonseed oil content, which will provide a strategic avenue for marker-assisted breeding aimed at incorporating desirable traits into cultivated cotton varieties.

Identifying components of the Shewanella phage LambdaSo lysis system

Phage-induced lysis of Gram-negative bacterial hosts usually requires a set of phage lysis proteins, a holin, an endopeptidase, and a spanin system, to disrupt each of the three cell envelope layers. Genome annotations and previous studies identified a gene region in the Shewanella oneidensis prophage LambdaSo, which comprises potential holin- and endolysin-encoding genes but lacks an obvious spanin system. By a combination of candidate approaches, mutant screening, characterization, and microscopy, we found that LambdaSo uses a pinholin/signal-anchor-release (SAR) endolysin system to induce proton leakage and degradation of the cell wall. Between the corresponding genes, we found that two extensively nested open-reading frames encode a two-component spanin module Rz/Rz1. Unexpectedly, we identified another factor strictly required for LambdaSo-induced cell lysis, the phage protein Lcc6. Lcc6 is a transmembrane protein of 65 amino acid residues with hitherto unknown function, which acts at the level of holin in the cytoplasmic membrane to allow endolysin release. Thus, LambdaSo-mediated cell lysis requires at least four protein factors (pinholin, SAR endolysin, spanin, and Lcc6). The findings further extend the known repertoire of phage proteins involved in host lysis and phage egress.

IMPORTANCE

Lysis of bacteria can have multiple consequences, such as the release of host DNA to foster robust biofilm. Phage-induced lysis of Gram-negative cells requires the disruption of three layers, the outer and inner membranes and the cell wall. In most cases, the lysis systems of phages infecting Gram-negative cells comprise holins to disrupt or depolarize the membrane, thereby releasing or activating endolysins, which then degrade the cell wall. This, in turn, allows the spanins to become active and fuse outer and inner membranes, completing cell envelope disruption and allowing phage egress. Here, we show that the presence of these three components may not be sufficient to allow cell lysis, implicating that also in known phages, further factors may be required.

On knowing a gene: A distributional hypothesis of gene function

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

Genome-wide identification, characterization and expression analysis of the DUF668 gene family in tomato

The domain of unknown function 668 (DUF668) is a gene family that may play a key role in plant growth and development as well as in responding to adversity coercion stresses. However, the DUF668 gene family has not yet been well identified and characterized in tomato. In this study, a total of nine putative SlDUF668 genes were identified in tomato, distributed on six chromosomes. Phylogenetic analyses revealed that SlDUF668 proteins were classified into two major groups. Members within the same group largely displayed analogous gene structure and conserved motif compositions. Several cis-elements were exhibited in the upstream sequences of the SlDUF668 genes, including elements implicated in plant growth and development processes, abiotic stress and hormone responses. Further, the study assessed the expression patterns of the SlDUF668 gene family in various tomato tissues, five plant hormones treatments, three abiotic stresses using qRT-PCR. The SlDUF668 genes expressed ubiquitously in various tissues, and five genes (SlDUF668-04, SlDUF668-06, SlDUF668-07, SlDUF668-08 and SlDUF668-09) showed tissue specificity. And SlDUF668 genes responded to abiotic stresses such as salt, drought and cold to varying degrees. Overall, our study provided a base for the tomato DUF668 gene family and laid a foundation for further understanding the functional characteristics of DUF668 genes in tomato plants.

Genome-wide characterization and analysis of rice DUF247 gene family

BACKGROUND

The domain of unknown function 247 (DUF247) proteins is involved in plant development and stress response. Rice is an important worldwide cereal crop, although an increasing number of DUF proteins have been identified, the understanding of DUF proteins is still very limited in rice.

RESULTS

In this study, we identified 69 genes that encode DUF247 proteins in the rice (Oryza sativa) genome by homology searches and domain prediction. All the OsDUF247 proteins were classified into four major groups (I, II, III and IV) by phylogenetic analysis. Remarkably, OsDUF247 genes clustered on the chromosomes solely show close phylogenetic relationships, suggesting that gene duplications have driven the expansion of the DUF247 gene family in the rice genome. Tissue profile analysis showed that most DUF247 genes expressed at constitutive levels in seedlings, roots, stems, and leaves, except for seven genes (LOC_Os01g21670, LOC_Os03g19700, LOC_Os05g04060, LOC_Os08g26820, LOC_Os08g26840, LOC_Os08g26850 and LOC_Os09g13410) in panicles. These seven genes were induced by various abiotic stress, including cold, drought, heat, hormone treatment, and especially salt, as demonstrated by further experimental analysis. DUF247 proteins contain transmembrane domains located on the membrane, suggesting their significant roles in rice development and adaptation to the environment.

CONCLUSIONS

These findings lay the foundation for functional characterizations of DUF247 genes to unravel their exact role in rice cultivars.

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

MOTIVATION

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

RESULTS

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

AVAILABILITY AND IMPLEMENTATION

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

Genome-wide identification, gene expression and haplotype analysis of the rhomboid-like gene family in wheat (Triticum aestivum L.)

The rhomboid-like (RBL) gene encodes serine protease, which plays an important role in the response to cell development and diverse stresses. However, genome-wide identification, expression profiles, and haplotype analysis of the RBL family genes have not been performed in wheat (Triticum aestivum L.). This study investigated the phylogeny and diversity of the RBL family genes in the wheat genome through various approaches, including gene structure analysis, evolutionary relationship analysis, promoter cis-acting element analysis, expression pattern analysis, and haplotype analysis. The 41 TaRBL genes were identified and divided into five subfamilies in the wheat genome. RBL family genes were expanded through segmented duplication and purification selection. The cis-element analysis revealed their involvement in various stress responses and plant development. The results of RNA-seq and quantitative real-time-PCR showed that TaRBL genes displayed higher expression levels in developing spike/grain and were differentially regulated under polyethylene glycol, NaCl, and abscisic acid treatments, indicating their roles in grain development and abiotic stress response. A kompetitive allele-specific PCR molecular marker was developed to confirm the single nucleotide polymorphism of TaRBL14a gene in 263 wheat accessions. We found that the elite haplotype TaRBL14a-Hap2 showed a significantly higher 1000-grain weight than TaRBL14a-Hap11 in at least three environments, and the TaRBL14a-Hap2 was positively selected in wheat breeding. The findings will provide a good insight into the evolutionary and functional characteristics of the TaRBL genes family in wheat and lay the foundation for future exploration of the regulatory mechanisms of TaRBL genes in plant growth and development, as well as their response to abiotic stresses.

Chromatin remodeling protein BPTF regulates transcriptional stability in planarian stem cells

Trimethylation of histone H3 lysine 4 (H3K4me3) correlates strongly with gene expression in many different organisms, yet the question of whether it plays a causal role in transcriptional activity remains unresolved. Although H3K4me3 does not directly affect chromatin accessibility, it can indirectly affect genome accessibility by recruiting the ATP-dependent chromatin remodeling complex NuRF (Nucleosome Remodeling Factor). The largest subunit of NuRF, BPTF/NURF301, binds H3K4me3 specifically and recruits the NuRF complex to loci marked by this modification. Studies have shown that the strength and duration of BPTF binding likely also depends on additional chromatin features at these loci, such as lysine acetylation and variant histone proteins. However, the exact details of this recruitment mechanism vary between studies and have largely been tested in vitro. Here, we use stem cells isolated directly from live planarian animals to investigate the role of BPTF in regulating chromatin accessibility in vivo. We find that BPTF operates at gene promoters and is most effective at facilitating transcription at genes marked by Set1-dependent H3K4me3 peaks, which are significantly broader than those added by the lysine methyltransferase MLL1/2. Moreover, BPTF is essential for planarian stem cell biology and its loss of function phenotype mimics that of Set1 knockdown. Together, these data suggest that BPTF and H3K4me3 are important mediators of both transcription and in vivo stem cell function.

The molecular evolution of cancer associated genes in mammals

Cancer is a disease that many multicellular organisms have faced for millions of years, and species have evolved various tumour suppression mechanisms to control oncogenesis. Although cancer occurs across the tree of life, cancer related mortality risks vary across mammalian orders, with Carnivorans particularly affected. Evolutionary theory predicts different selection pressures on genes associated with cancer progression and suppression, including oncogenes, tumour suppressor genes and immune genes. Therefore, we investigated the evolutionary history of cancer associated gene sequences across 384 mammalian taxa, to detect signatures of selection across categories of oncogenes (GRB2, FGL2 and CDC42), tumour suppressors (LITAF, Casp8 and BRCA2) and immune genes (IL2, CD274 and B2M). This approach allowed us to conduct a fine scale analysis of gene wide and site-specific signatures of selection across mammalian lineages under the lens of cancer susceptibility. Phylogenetic analyses revealed that for most species the evolution of cancer associated genes follows the species' evolution. The gene wide selection analyses revealed oncogenes being the most conserved, tumour suppressor and immune genes having similar amounts of episodic diversifying selection. Despite BRCA2's status as a key caretaker gene, episodic diversifying selection was detected across mammals. The site-specific selection analyses revealed that the two apoptosis associated domains of the Casp8 gene of bats (Chiroptera) are under opposing forces of selection (positive and negative respectively), highlighting the importance of site-specific selection analyses to understand the evolution of highly complex gene families. Our results highlighted the need to critically assess different types of selection pressure on cancer associated genes when investigating evolutionary adaptations to cancer across the tree of life. This study provides an extensive assessment of cancer associated genes in mammals with highly representative, and substantially large sample size for a comparative genomic analysis in the field and identifies various avenues for future research into the mechanisms of cancer resistance and susceptibility in mammals.

Genome-wide identification of GH28 family and insight into its contributions to pod shattering resistance in Brassica napus L

Rapeseed (Brassica napus L.), accounts for nearly 16% of vegetable oil, is the world's second produced oilseed. However, pod shattering has caused significant yield loses in rapeseed production, particularly during mechanical harvesting. The GH28 genes can promote pod shattering by changing the structure of the pod cell wall in Arabidopsis. However, the role of the GH28 gene family in rapeseed was largely unknown. Therefore, a genome-wide comprehensive analysis was conducted to classify the role of GH28 gene family on rapeseed pod shattering. A total of 37 BnaGH28 genes in the rapeseed genome were identified. These BnaGH28s can be divided into five groups (Group A-E), based on phylogenetic and synteny analysis. Protein property, gene structure, conserved motif, cis-acting element, and gene expression profile of BnaGH28 genes in the same group were similar. Specially, the expression level of genes in group A-D was gradually decreased, but increased in group E with the development of silique. Among eleven higher expressed genes in group E, two BnaGH28 genes (BnaA07T0199500ZS and BnaC06T0206500ZS) were significantly regulated by IAA or GA treatment. And the significant effects of BnaA07T0199500ZS variation on pod shattering resistance were also demonstrated in present study. These results could open a new window for insight into the role of BnaGH28 genes on pod shattering resistance in rapeseed.