SMART: recent updates, new developments and status in 2020

Positive selection and functional diversification of transcription factor Cmr1 homologs in Alternaria

Transcription factor Cmr1 (Colletotrichum melanin regulation 1) and its homologs in several plant fungal pathogens are the regulators of the 1,8-dihydroxynaphthalene (DHN)-melanin biosynthesis pathway and have evolved functional diversification in morphology and pathogenicity. The fungal genus Alternaria comprises the group of "black fungi" that are rich in DHN-melanin in the primary cell wall and septa of the conidia. Some Alternaria species cause many economically important plant diseases worldwide. However, the evolution and function of Cmr1 homologs in Alternaria remain poorly understood. Here, we identified a total of forty-two Cmr1 homologs from forty-two Alternaria spp. and all contained one additional diverse fungal specific transcription factor motif. Phylogenetic analysis indicated the division of these homologs into five major clades and three branches. Dated phylogeny showed the A and D clades diverged latest and earliest, respectively. Molecular evolutionary analyses revealed that three amino acid sites of Cmr1 homologs in Alternaria were the targets of positive selection. Asmr1, the homolog of Cmr1 in the potato early blight pathogen, Alternaria solani was amplified and displayed the sequence conservation at the amino acid level in different A. solani isolates. Asmr1 was further confirmed to have the transcriptional activation activity and was upregulated during the early stage of potato infection. Deletion of asmr1 led to the decreased melanin content and pathogenicity, deformed conidial morphology, and responses to cell wall and fungicide stresses in A. solani. These results suggest positive selection and functional divergence have played a role in the evolution of Cmr1 homologs in Alternaria. KEY POINTS: • Cmr1 homologs were under positive selection in Alternaria species • Asmr1 is a functional transcription factor, involved in spore development, melanin biosynthesis, pathogenicity, and responses to cell wall and fungicide stresses in A. solani • Cmr1 might be used as a potential taxonomic marker of the genus Alternaria.

Genome-wide identification and analysis of DEAD-box RNA helicases in Gossypium hirsutum

DEAD-box RNA helicases, a prominent subfamily within the RNA helicase superfamily 2 (SF2), play crucial roles in the growth, development, and abiotic stress responses of plants. This study identifies 146 DEAD-box RNA helicase genes (GhDEADs) and categorizes them into four Clades (Clade A-D) through phylogenetic analysis. Promoter analysis reveals cis-acting elements linked to plant responses to light, methyl jasmonate (MeJA), abscisic acid (ABA), low temperature, and drought. RNA-seq data demonstrate that Clade C GhDEADs exhibit elevated and ubiquitous expression across different tissues, validating their connection to leaf development through real-time quantitative polymerase chain reaction (RT-qPCR) analysis. Notably, over half of GhDEADs display up-regulation in the leaves of virus-induced gene silencing (VIGS) plants of GhVIR-A/D (members of m6A methyltransferase complex, which regulate leaf morphogenesis). In conclusion, this study offers a comprehensive insight into GhDEADs, emphasizing their potential involvement in leaf development.

Cloning and expression analysis of Janus activated kinase family genes from spotted seabass (Lateolabrax maculatus)

Janus kinases (JAKs) are important components of the JAK-STAT signaling pathway and play vital roles in innate immunity, autoimmune diseases, and inflammation. However, information about JAKs remains largely unknown in the spotted seabass, a fish species of Perciformes with great commercial value in the aquaculture industry. The aims of this study are to obtain the complete cDNA sequences of JAKs (JAK1, JAK2A, JAK2B, JAK3 and TYK2) from spotted seabass and to investigate their roles upon stimulation with lipopolysaccharides (LPS) and Edwardsiella tarda, using RT-PCR, PCR and qRT-PCR methods. All five JAK genes from the spotted seabass, each encode more than 1100 amino acids residues. JAK1 and JAK3 consist of 24 exons and 23 introns, whereas JAK2A, JAK2B and TYK2 consist of 23 exons and 22 introns. Furthermore, these five spotted seabass JAKs share high sequence identities with those of other fish species in protein domain analysis, synteny analysis, and phylogenetic analysis. Moreover, these five JAK genes were ubiquitously expressed in all tissues examined from healthy fish, and inducible expressions of JAKs were observed in the intestine, gill, head kidney, and spleen following LPS treatment or E. tarda infection. These findings indicate that all these JAK genes are involved in the antibacterial immunity of the spotted seabass and provide a basis for further understanding the mechanism of JAKs antibacterial response in the spotted sea bass.

Exploring the genetic diversity pattern of PvEBP/DBP2: A promising candidate for an effective Plasmodium vivax vaccine

Malaria remains a public health challenge. Since many control strategies have proven ineffective in eradicating this disease, new strategies are required, among which the design of a multivalent vaccine stands out. However, the effectiveness of this strategy has been hindered, among other reasons, by the genetic diversity observed in parasite antigens. In Plasmodium vivax, the Erythrocyte Binding Protein (PvEBP, also known as DBP2) is an alternate ligand to Duffy Binding Protein (DBP); given its structural resemblance to DBP, EBP/DBP2 is proposed as a promising antigen for inclusion in vaccine design. However, the extent of genetic diversity within the locus encoding this protein has not been comprehensively assessed. Thus, this study aimed to characterize the genetic diversity of the locus encoding the P. vivax EBP/DBP2 protein and to determine the evolutionary mechanisms modulating this diversity. Several intrapopulation genetic variation parameters were estimated using 36 gene sequences of PvEBP/DBP2 from Colombian P. vivax clinical isolates and 186 sequences available in databases. The study then evaluated the worldwide genetic structure and the evolutionary forces that may influence the observed patterns of genetic variation. It was found that the PvEBP/DBP2 gene exhibits one of the lowest levels of genetic diversity compared to other vaccine-candidate antigens. Four major haplotypes were shared worldwide. Analysis of the protein's 3D structure and epitope prediction identified five regions with potential antigenic properties. The results suggest that the PvEBP/DBP2 protein possesses ideal characteristics to be considered when designing a multivalent effective antimalarial vaccine against P. vivax.

Comprehensive identification and analysis of DUF640 genes associated with rice growth

Protein domains with conserved amino acid sequences and uncharacterized functions are called domains of unknown function (DUF). The DUF640 gene family plays a crucial role in plant growth, particularly in light regulation, floral organ development, and fruit development. However, there exists a lack of systematic understanding of the evolutionary relationships and functional differentiation of DUF640 within the Oryza genus. In this study, 61 DUF640 genes were identified in the Oryza genus. The expression of DUF640s is induced by multiple hormonal stressors including abscisic acid (ABA), cytokinin (CK), ethylene (ETH), and indole-3-acetic acid (IAA). Specifically, OiDUF640-10 expression significantly increased after ETH treatment. Transgenic experiments showed that overexpressing OiDUF640-10 lines were sensitive to ETH, and seedling length was obstructed. Evolutionary analysis revealed differentiation of the OiDUF640-10 gene in O. sativa ssp. indica and japonica varieties, likely driven by natural selection during the domestication of cultivated rice. These results indicate that OiDUF640-10 plays a vital role in the regulation of rice seedling length.

Discovery and characterization of the α-amylases cDNAs from Enchytraeus albidus shed light on the evolution of "Enchytraeus-Eisenia type" Amy homologs in Annelida

Although enchytraeids have gained popularity in scientific research, fundamental questions regarding their feeding ecology and biology remain largely unexplored. This study investigates α-amylases, major digestive enzymes responsible for hydrolyzing starch and similar polysaccharides into sugars, in Enchytraeus albidus. Genetic data related to α-amylases is currently lacking for the family Enchytraeidae but also for the entire Annelida. To detect and identify coding sequences of the expressed α-amylase genes in COI-monohaplotype culture (PL-A strain) of E. albidus, we used classical "gene fishing" and transcriptomic approaches. We also compared coding sequence variants of α-amylase retrieved from transcriptomic data related to freeze-tolerant strains. Our results reveal that E. albidus possesses two distinct α-amylase genes (Amy I and Amy II) that are homologs to earthworm Eisenia fetida Ef-Amy genes. Different strains of E. albidus possess distinctive alleles of α-amylases with unique SNP patterns specific to a particular strain. Unlike Amy II, Amy I seems to be a highly polymorphic and multicopy gene. The domain architecture of the putative Amy proteins was found the same as for classical animal α-amylases with ABC-domains. A characteristic feature of Amy II is the lack of GHGA motif in the flexible loop region, similarly to many insect amylases. We identified "Enchytraeus-Eisenia type" α-amylase homologs in other clitellates and polychaetes, indicating the ancestral origin of Amy I/II proteins in Annelida. This study provides the first insight into the endogenous non-proteolytic digestive enzyme genes in potworms, discusses the evolution of Amy α-amylases in Annelida, and explores phylogenetic implications.

TaTPS11 enhances wheat cold resistance by regulating source-sink factor

The presence of sugar in plant tissue can lead to an increase in the osmotic pressure within cells, a decrease in the freezing point of plants, and protection against ice crystal damage to the tissue. Trehalose is closely related to sucrose, which comprises the largest proportion of sugar and has become a hot topic of research in recent years. Our previous studies have confirmed that a key trehalose synthesis gene, TaTPS11, from the cold-resistant winter wheat DM1, could enhance the cold resistance of plants by increasing sugar content. However, the underlying mechanism behind this phenomenon remains unclear. In this study, we cloned TaTPS11-6D, edited TaTPS11-6D using CRISPR/Cas9 technology and transformed 'Fielder' to obtain T2 generation plants. We screened out OE3-3 and OE8-7 lines with significantly higher cold resistance than that of 'Fielder' and Cri 4-3 edited lines with significantly lower cold resistance than that of 'Fielder'. Low temperature storage limiting factors were measured for OE3-3, OE8-7 and Cri 4-3 treated at different temperatures.The results showed that TaTPS11-6D significantly increased the content of sugar in plants and the transfer of sugar from source to storage organs under cold conditions. The TaTPS11-6D significantly increased the levels of salicylic, jasmonic, and abscisic acids while also significantly decreasing the level of gibberellic acid. Our research improves the model of low temperature storage capacity limiting factor.

Genome-wide identification and expression pattern analysis of WRKY transcription factors in response to biotic and abiotic stresses in tea plants (Camellia sinensis)

Plants would encounter various biotic and abiotic stresses during the growth and development. WRKY transcription factors (TFs) as plant-specific TFs, play an important role in responding to various adverse circumstances. Despite some advances were achieved in functional studies of WRKY TFs in tea plants, systematic analysis of the involvement of CsWRKY TFs when facing cold, salt, drought stresses and pathogen and insect attack was lacked. In present study, a total of 78 CsWRKY TFs were identified following the genomic and transcript databases. The expression patterns of CsWRKYs in various organs of tea plants and the expression profiles in response to biotic and abiotic stresses were investigated by examining representative RNA-seq data. Moreover, the effects of hormone treatments (SA and MeJA) on the transcription levels of WRKY TFs were also investigated. The phylogenetic tree of CsWRKY TFs from different species indicated the functional diversity of WRKY TFs was not closely related to their protein classification. Concurrently, CsWRKY70-2 TF was identified as a positive regulator in response to drought stress. This study provided solid and valuable information, helping us better understand the functional diversity of CsWRKY TFs, and laid the foundation for further research on the function of key WRKY genes in tea plants.

KIF1A novel frameshift variant p.(Ser887Profs*64) exhibits clinical heterogeneity in a Pakistani family with hereditary sensory and autonomic neuropathy type IIC

Background: Hereditary sensory and autonomic neuropathies (HSANs) are rare heterogeneous group of neurological disorders caused by peripheral nerve deterioration. The HSANs sub-clinical classes have clinical and genetic overlap which often lead to misdiagnosis. In the present study a Pakistani family with five affected members suffering from severe neuropathy were genetically analyzed to identify the disease causative element in the family.

Methods: Genome wide high-density single nucleotide polymorphism (SNP) microarray analysis was carried out followed by whole exome sequencing of the affected proband and another affected sibling. Shared homozygous regions in all severely affected members were identified through homozygosity mapping approach.

Results: The largest homozygous region of 14.1 Mb shared by the five severely affected members of the family was identified on chromosome 2. Subsequent exome sequencing identified a novel single nucleotide deletion c.2658del; p.(Ser887Profs*64) in KIF1A. Segregation analysis revealed that this mutation was homozygous in all five affected individuals of the family with severe clinical manifestation, while members of the family that were heterozygous carriers shared abnormal skin features (scaly skin) only with the homozygous affected members.

Conclusions: A novel frameshift mutation p.(Ser887Profs*64) in KIF1A is the potential cause of severe HSANIIC in a Pakistani family along with incomplete penetrance in mutation carriers. We demonstrate that using a combination of different techniques not only strengthens the gene finding approach but also helps in proper sub-clinical characterization along with identification of mutated alleles exhibiting incomplete penetrance leading to intrafamilial clinical variability in HSAN group of inherited diseases.

The novel regulator HdrR controls the transcription of the heterodisulfide reductase operon hdrBCA in Methanosarcina barkeri

Methanogenic archaea play a key role in the global carbon cycle because these microorganisms remineralize organic compounds in various anaerobic environments. The microorganism Methanosarcina barkeri is a metabolically versatile methanogen, which can utilize acetate, methanol, and H2/CO2 to synthesize methane. However, the regulatory mechanisms underlying methanogenesis for different substrates remain unknown. In this study, RNA-seq analysis was used to investigate M. barkeri growth and gene transcription under different substrate regimes. According to the results, M. barkeri showed the best growth under methanol, followed by H2/CO2 and acetate, and these findings corresponded well with the observed variations in genes transcription abundance for different substrates. In addition, we identified a novel regulator, MSBRM_RS03855 (designated as HdrR), which specifically activates the transcription of the heterodisulfide reductase hdrBCA operon in M. barkeri. HdrR was able to bind to the hdrBCA operon promoter to regulate transcription. Furthermore, the structural model analyses revealed a helix-turn-helix domain, which is likely involved in DNA binding. Taken together, HdrR serves as a model to reveal how certain regulatory factors control the expression of key enzymes in the methanogenic pathway.IMPORTANCEThe microorganism Methanosarcina barkeri has a pivotal role in the global carbon cycle and contributes to global temperature homeostasis. The consequences of biological methanogenesis are far-reaching, including impacts on atmospheric methane and CO2 concentrations, agriculture, energy production, waste treatment, and human health. As such, reducing methane emissions is crucial to meeting set climate goals. The methanogenic activity of certain microorganisms can be drastically reduced by inhibiting the transcription of the hdrBCA operon, which encodes heterodisulfide reductases. Here, we provide novel insight into the mechanisms regulating hdrBCA operon transcription in the model methanogen M. barkeri. The results clarified that HdrR serves as a regulator of heterodisulfide reductase hdrBCA operon transcription during methanogenesis, which expands our understanding of the unique regulatory mechanisms that govern methanogenesis. The findings presented in this study can further our understanding of how genetic regulation can effectively reduce the methane emissions caused by methanogens.

Bioinformatics approach for structure modeling, vaccine design, and molecular docking of Brucella candidate proteins BvrR, OMP25, and OMP31

Brucellosis is a zoonotic disease with significant economic and healthcare costs. Despite the eradication efforts, the disease persists. Vaccines prevent disease in animals while antibiotics cure humans with limitations. This study aims to design vaccines and drugs for brucellosis in animals and humans, using protein modeling, epitope prediction, and molecular docking of the target proteins (BvrR, OMP25, and OMP31). Tertiary structure models of three target proteins were constructed and assessed using RMSD, TM-score, C-score, Z-score, and ERRAT. The best models selected from AlphaFold and I-TASSER due to their superior performance according to CASP 12 - CASP 15 were chosen for further analysis. The motif analysis of best models using MotifFinder revealed two, five, and five protein binding motifs, however, the Motif Scan identified seven, six, and eight Post-Translational Modification sites (PTMs) in the BvrR, OMP25, and OMP31 proteins, respectively. Dominant B cell epitopes were predicted at (44-63, 85-93, 126-137, 193-205, and 208-237), (26-46, 52-71, 98-114, 142-155, and 183-200), and (29-45, 58-82, 119-142, 177-198, and 222-251) for the three target proteins. Additionally, cytotoxic T lymphocyte epitopes were detected at (173-181, 189-197, and 202-210), (61-69, 91-99, 159-167, and 181-189), and (3-11, 24-32, 167-175, and 216-224), while T helper lymphocyte epitopes were displayed at (39-53, 57-65, 150-158, 163-171), (79-87, 95-108, 115-123, 128-142, and 189-197), and (39-47, 109-123, 216-224, and 245-253), for the respective target protein. Furthermore, structure-based virtual screening of the ZINC and DrugBank databases using the docking MOE program was followed by ADMET analysis. The best five compounds of the ZINC database revealed docking scores ranged from (- 16.8744 to - 15.1922), (- 16.0424 to - 14.1645), and (- 14.7566 to - 13.3222) for the BvrR, OMP25, and OMP31, respectively. These compounds had good ADMET parameters and no cytotoxicity, while DrugBank compounds didn't meet Lipinski's rule criteria. Therefore, the five selected compounds from the ZINC20 databases may fulfill the pharmacokinetics and could be considered lead molecules for potentially inhibiting Brucella's proteins.

Identification of the mutual gliding locus as a factor for gut colonization in non-native bee hosts using the ARTP mutagenesis

BACKGROUND

The gut microbiota and their hosts profoundly affect each other's physiology and evolution. Identifying host-selected traits is crucial to understanding the processes that govern the evolving interactions between animals and symbiotic microbes. Current experimental approaches mainly focus on the model bacteria, like hypermutating Escherichia coli or the evolutionary changes of wild stains by host transmissions. A method called atmospheric and room temperature plasma (ARTP) may overcome the bottleneck of low spontaneous mutation rates while maintaining mild conditions for the gut bacteria.

RESULTS

We established an experimental symbiotic system with gnotobiotic bee models to unravel the molecular mechanisms promoting host colonization. By in vivo serial passage, we tracked the genetic changes of ARTP-treated Snodgrassella strains from Bombus terrestris in the non-native honeybee host. We observed that passaged isolates showing genetic changes in the mutual gliding locus have a competitive advantage in the non-native host. Specifically, alleles in the orphan mglB, the GTPase activating protein, promoted colonization potentially by altering the type IV pili-dependent motility of the cells. Finally, competition assays confirmed that the mutations out-competed the ancestral strain in the non-native honeybee gut but not in the native host.

CONCLUSIONS

Using the ARTP mutagenesis to generate a mutation library of gut symbionts, we explored the potential genetic mechanisms for improved gut colonization in non-native hosts. Our findings demonstrate the implication of the cell mutual-gliding motility in host association and provide an experimental system for future study on host-microbe interactions. Video Abstract.

GpEF1A: a novel lysine methyltransferase gene from Gypsophila perfoliata L. involved in boron homeostasis

Rapid accumulation of boron (B) leads to toxicity in plant tissues, and the narrow gap between deficiency and toxicity makes it difficult to adjust essential B levels in soil for plant productivity. Therefore, understanding different aspects of B tolerance is necessary to provide new and valid solutions to B toxicity. Gypsophila perfoliata stands out as a remarkable example of a B-tolerant plant, with a natural propensity to thrive in environments such as B mines and soils enriched with high levels of B. In this study, a yeast functional screening experiment was conducted using cDNA libraries from G. perfoliata leaf and root cells for B tolerance. Ten colonies from the leaf library grew in 80 mm boric acid, while none emerged from the root library. Analysis of isolated cDNAs showed identical sequences and a unique motif related to B tolerance. The gene GpEF1A was identified in the tolerant yeast colonies, with predicted structural features suggesting its role, and RT-qPCR indicating increased expression under B stress. A regulatory role for EF1A lysine methylation was proposed in mammalian cells and fungi because of its dynamic and inducible nature under environmental constraints. This could also be relevant for plant cells, as the high similarity of the GpEF1A gene in some salt-tolerant plants might indicate the upregulation of EF1A as a conserved way to cope with abiotic stress conditions. This report represents the first instance of involvement of GpEF1A in B tolerance, and further detailed studies are necessary to understand other components of this tolerance mechanism.

Discovery of fungal onoceroid triterpenoids through domainless enzyme-targeted global genome mining

Genomics-guided methodologies have revolutionized the discovery of natural products. However, a major challenge in the field of genome mining is determining how to selectively extract biosynthetic gene clusters (BGCs) for untapped natural products from numerous available genome sequences. In this study, we developed a fungal genome mining tool that extracts BGCs encoding enzymes that lack a detectable protein domain (i.e., domainless enzymes) and are not recognized as biosynthetic proteins by existing bioinformatic tools. We searched for BGCs encoding a homologue of Pyr4-family terpene cyclases, which are representative examples of apparently domainless enzymes, in approximately 2000 fungal genomes and discovered several BGCs with unique features. The subsequent characterization of selected BGCs led to the discovery of fungal onoceroid triterpenoids and unprecedented onoceroid synthases. Furthermore, in addition to the onoceroids, a previously unreported sesquiterpene hydroquinone, of which the biosynthesis involves a Pyr4-family terpene cyclase, was obtained. Our genome mining tool has broad applicability in fungal genome mining and can serve as a beneficial platform for accessing diverse, unexploited natural products.

Genome-wide identification and examination of the wheat glycosyltransferase family 43 regulation during Fusarium graminearum infection

In Arabidopsis and rice, the glycosyltransferase (GT) 43 family is involved in xylan synthesis. However, there have been limited reports on the study of the TaGT43 family in wheat. In this study, 28 TaGT43 family members were identified in wheat (Triticum aestivum L.) and clustered into three major groups based on the similarity of amino acid sequences. The results of the TaGT43 family's conserved motif and gene structure analyses agree with this result. Collinearity analysis revealed segmental duplications mainly promoted TaGT43 family expansion. cis-Acting element analysis revealed that the TaGT43 genes were involved in the light response, phytohormone response, abiotic/biotic stress response, and growth and development. The TaGT43 family showed a tissue-specific expression pattern, primarily expressed in roots and stems. Besides, the transcriptional and expression levels of multiple TaGT43 genes were upregulated during the infection of F. graminearum. According to metabolomics studies, F. graminearum infection affected the phenylpropanoid biosynthesis pathway in wheat, a critical factor in cell wall construction. Furthermore, GO enrichment analysis indicated that the TaGT43 genes play a significant role in cell wall organization. Based on these results, it may be concluded that the TaGT43 family mediates cell wall organization in response to F. graminearum infection.

Understanding Genetic Risks: Computational Exploration of Human β-Synuclein nsSNPs and their Potential Impact on Structural Alteration

Synucleins are pivotal in neurodegenerative conditions. Beta-synuclein (β-synuclein) is part of the synuclein protein family alongside alpha-synuclein (α-synuclein) and gamma-synuclein (γ-synuclein). These proteins, found mainly in brain tissue and cancers, are soluble and unstructured. β-synuclein shares significant similarity with α-synuclein, especially in their N-terminus, with a 90% match. However, their aggregation tendencies differ significantly. While α-synuclein aggregation is believed to be counteracted by β-synuclein, which occurs in conditions like Parkinson's disease, β-synuclein may counteract α-synuclein's toxic effects on the nervous system, offering potential treatment for neurodegenerative diseases. Under normal circumstances, β-synuclein may guard against disease by interacting with α-synuclein. Yet, in pathological environments with heightened levels or toxic substances, it might contribute to disease. Our research aims to explore potential harmful mutations in the β-synuclein using computational tools to predict their destabilizing impact on protein structure. Consensus analysis revealed rs1207608813 (A63P), rs1340051870 (S72F), and rs1581178262 (G36C) as deleterious. These findings highlight the intricate relationship between nsSNPs and protein function, shedding light on their potential implications in disease pathways. Understanding the structural consequences of nsSNPs is crucial for elucidating their role in pathogenesis and developing targeted therapeutic interventions. Our results offer a robust computational framework for identifying neurodegenerative disorder-related mutations from SNP datasets, potentially reducing the costs associated with experimental characterization.

Molecular characterization revealed the role of thaumatin-like proteins of Rhizoctonia solani AG4-JY in inducing maize disease resistance

The gene family of thaumatin-like proteins (TLPs) plays a crucial role in the adaptation of organisms to environmental stresses. In recent years, fungal secreted proteins (SP) with inducing disease resistance activity in plants have emerged as important elicitors in the control of fungal diseases. Identifying SPs with inducing disease resistance activity and studying their mechanisms are crucial for controlling sheath blight. In the present study, 10 proteins containing the thaumatin-like domain were identified in strain AG4-JY of Rhizoctonia solani and eight of the 10 proteins had signal peptides. Analysis of the TLP genes of the 10 different anastomosis groups (AGs) showed that the evolutionary relationship of the TLP gene was consistent with that between different AGs of R. solani. Furthermore, it was found that RsTLP3, RsTLP9 and RsTLP10 were regarded as secreted proteins for their signaling peptides exhibited secretory activity. Prokaryotic expression and enzyme activity analysis revealed that the three secreted proteins possess glycoside hydrolase activity, suggesting they belong to the TLP family. Additionally, spraying the crude enzyme solution of the three TLP proteins could enhance maize resistance to sheath blight. Further analysis showed that genes associated with the salicylic acid and ethylene pathways were up-regulated following RsTLP3 application. The results indicated that RsTLP3 had a good application prospect in biological control.

UBC Gene Family Analysis in Salvia castanea and Roles of ScUBC2/5 Genes under Abiotic Stress

Salvia castanea Diels, a relative of the medicinal plant Salvia miltiorrhiza Bunge, belongs to the genus Salvia and family Lamiaceae. Ubiquitin-conjugating enzyme E2 (UBC) is an important ubiquitin-binding enzyme in protein ubiquitination. This study aimed to analyze the regulatory role of UBC genes, particularly ScUBC2/5, on the growth and adaptation of S. castanea to extreme environments including cold or drought stress. We identified nine UBC genes in S. castanea and found that these genes were extremely stable and more highly expressed in the roots than other tissues. This suggested that UBC genes might play a role in promoting root adaptation to cold and dry environments. Further analysis of UBC gene expression in hairy roots under cold (4 °C) and UV stress also confirmed their importance under stress. The contents of tanshinone and salvianolic acid in hairy roots with the overexpression of ScUBC2/5 were increased compared to non-transgenic wild type, and the cold and UV resistance of hairy roots was increased compared with that of wild type. Together, these findings highlighted the role of ScUBC2/5 in enhancing secondary metabolite accumulation and regulation in response to cold and ultraviolet stress in S. castanea, providing a new perspective for genetic improvement in its phytochemistry.

Molecular cloning and characterization of a salt overly sensitive3 (SOS3) gene from the halophyte Pongamia

A high concentration of sodium (Na+) is the primary stressor for plants in high salinity environments. The Salt Overly Sensitive (SOS) pathway is one of the best-studied signal transduction pathways, which confers plants the ability to export too much Na+ out of the cells or translocate the cytoplasmic Na+ into the vacuole. In this study, the Salt Overly Sensitive3 (MpSOS3) gene from Pongamia (Millettia pinnata Syn. Pongamia pinnata), a semi-mangrove, was isolated and characterized. The MpSOS3 protein has canonical EF-hand motifs conserved in other calcium-binding proteins and an N-myristoylation signature sequence. The MpSOS3 gene was significantly induced by salt stress, especially in Pongamia roots. Expression of the wild-type MpSOS3 but not the mutated nonmyristoylated MpSOS3-G2A could rescue the salt-hypersensitive phenotype of the Arabidopsis sos3-1 mutant, which suggested the N-myristoylation signature sequence of MpSOS3 was required for MpSOS3 function in plant salt tolerance. Heterologous expression of MpSOS3 in Arabidopsis accumulated less H2O2, superoxide anion radical (O2-), and malondialdehyde (MDA) than wild-type plants, which enhanced the salt tolerance of transgenic Arabidopsis plants. Under salt stress, MpSOS3 transgenic plants accumulated a lower content of Na+ and a higher content of K+ than wild-type plants, which maintained a better K+/Na+ ratio in transgenic plants. Moreover, no development and growth discrepancies were observed in the MpSOS3 heterologous overexpression plants compared to wild-type plants. Our results demonstrated that the MpSOS3 pathway confers a conservative salt-tolerant role and provided a foundation for further study of the SOS pathway in Pongamia.

FvMbp1-Swi6 complex regulates vegetative growth, stress tolerance, and virulence in Fusarium verticillioides

The mycotoxigenic fungus Fusarium verticillioides is a common pathogen of grain and medicine that contaminates the host with fumonisin B1 (FB1) mycotoxin, poses serious threats to human and animal health. Therefore, it is crucial to unravel the regulatory mechanisms of growth, and pathogenicity of F. verticillioides. Mbp1 is a component of the MluI cell cycle box binding factor complex and acts as an APSES-type transcription factor that regulates cell cycle progression. However, no information is available regarding its role in F. verticillioides. In this study, we demonstrate that FvMbp1 interacts with FvSwi6 that acts as the cell cycle transcription factor, to form the heteromeric transcription factor complexes in F. verticillioides. Our results show that ΔFvMbp1 and ΔFvSwi6 both cause a severe reduction of vegetative growth, conidiation, and increase tolerance to diverse environmental stresses. Moreover, ΔFvMbp1 and ΔFvSwi6 dramatically decrease the virulence of the pathogen on the stalk and ear of maize. Transcriptome profiling show that FvMbp1-Swi6 complex co-regulates the expression of genes associated with multiple stress responses. These results indicate the functional importance of the FvMbp1-Swi6 complex in the filamentous fungi F. verticillioides and reveal a potential target for the effective prevention and control of Fusarium diseases.

Exome sequencing implicates ancestry-related Mendelian variation at SYNE1 in childhood-onset essential hypertension

Childhood-onset essential hypertension (COEH) is an uncommon form of hypertension that manifests in childhood or adolescence and, in the United States, disproportionately affects children of African ancestry. The etiology of COEH is unknown, but its childhood onset, low prevalence, high heritability, and skewed ancestral demography suggest the potential to identify rare genetic variation segregating in a Mendelian manner among affected individuals and thereby implicate genes important to disease pathogenesis. However, no COEH genes have been reported to date. Here, we identify recessive segregation of rare and putatively damaging missense variation in the spectrin domain of spectrin repeat containing nuclear envelope protein 1 (SYNE1), a cardiovascular candidate gene, in 3 of 16 families with early-onset COEH without an antecedent family history. By leveraging exome sequence data from an additional 48 COEH families, 1,700 in-house trios, and publicly available data sets, we demonstrate that compound heterozygous SYNE1 variation in these COEH individuals occurred more often than expected by chance and that this class of biallelic rare variation was significantly enriched among individuals of African genetic ancestry. Using in vitro shRNA knockdown of SYNE1, we show that reduced SYNE1 expression resulted in a substantial decrease in the elasticity of smooth muscle vascular cells that could be rescued by pharmacological inhibition of the downstream RhoA/Rho-associated protein kinase pathway. These results provide insights into the molecular genetics and underlying pathophysiology of COEH and suggest a role for precision therapeutics in the future.

Identification and characterization of the CRK gene family in the wheat genome and analysis of their expression profile in response to high temperature-induced male sterility

Cysteine-rich receptor-like kinases (CRKs) play many important roles during plant development, including defense responses under both biotic and abiotic stress, reactive oxygen species (ROS) homeostasis, callose deposition and programmed cell death (PCD). However, there are few studies on the involvement of the CRK family in male sterility due to heat stress in wheat (Triticum aestivum L.). In this study, a genome-wide characterization of the CRK family was performed to investigate the structural and functional attributes of the wheat CRKs in anther sterility caused by heat stress. A total of 95 CRK genes were unevenly distributed on 18 chromosomes, with the most genes distributed on chromosome 2B. Paralogous homologous genes with Ka/Ks ratios less than 1 may have undergone strong purifying selection during evolution and are more functionally conserved. The collinearity analysis results of CRK genes showed that wheat and Arabidopsis (A. thaliana), foxtail millet, Brachypodium distachyon (B. distachyon), and rice have three, 12, 15, and 11 pairs of orthologous genes, respectively. In addition, the results of the network interactions of genes and miRNAs showed that five miRNAs were in the hub of the interactions map, namely tae-miR9657b-5p, tae-miR9780, tae-miR9676-5p, tae-miR164, and tae-miR531. Furthermore, qRT-PCR validation of the six TaCRK genes showed that they play key roles in the development of the mononuclear stage anthers, as all six genes were expressed at highly significant levels in heat-stressed male sterile mononuclear stage anthers compared to normal anthers. We hypothesized that the TaCRK gene is significant in the process of high-temperature-induced sterility in wheat based on the combination of anther phenotypes, paraffin sections, and qRT-PCR data. These results improve our understanding of their relationship.

Genome-Wide Analysis of Transcription Factor R2R3-MYB Gene Family and Gene Expression Profiles during Anthocyanin Synthesis in Common Walnut (Juglans regia L.)

The R2R3-MYB gene family, encoding plant transcriptional regulators, participates in many metabolic pathways of plant physiology and development, including flavonoid metabolism and anthocyanin synthesis. This study proceeded as follows: the JrR2R3-MYB gene family was analyzed genome-wide, and the family members were identified and characterized using the high-quality walnut reference genome "Chandler 2.0". All 204 JrR2R3-MYBs were established and categorized into 30 subgroups via phylogenetic analysis. JrR2R3-MYBs were unevenly distributed over 16 chromosomes. Most JrR2R3-MYBs had similar structures and conservative motifs. The cis-acting elements exhibit multiple functions of JrR2R3-MYBs such as light response, metabolite response, and stress response. We found that the expansion of JrR2R3-MYBs was mainly caused by WGD or segmental duplication events. Ka/Ks analysis indicated that these genes were in a state of negative purifying selection. Transcriptome results suggested that JrR2R3-MYBs were widely entangled in the process of walnut organ development and differentially expressed in different colored varieties of walnuts. Subsequently, we identified 17 differentially expressed JrR2R3-MYBs, 9 of which may regulate anthocyanin biosynthesis based on the results of a phylogenetic analysis. These genes were present in greater expression levels in 'Zijing' leaves than in 'Lvling' leaves, as revealed by the results of qRT-PCR experiments. These results contributed to the elucidation of the functions of JrR2R3-MYBs in walnut coloration. Collectively, this work provides a foundation for exploring the functional characteristics of the JrR2R3-MYBs in walnuts and improving the nutritional value and appearance quality of walnuts.

Assessment of the hypothetical protein BB0616 in the murine infection of Borrelia burgdorferi

bb0616 of Borrelia burgdorferi, the Lyme disease pathogen, encodes a hypothetical protein of unknown function. In this study, we showed that BB0616 was not surface-exposed or associated with the membrane through localization analyses using proteinase K digestion and cell partitioning assays. The expression of bb0616 was influenced by a reduced pH but not by growth phases, elevated temperatures, or carbon sources during in vitro cultivation. A transcriptional start site for bb0616 was identified by using 5' rapid amplification of cDNA ends, which led to the identification of a functional promoter in the 5' regulatory region upstream of bb0616. By analyzing a bb0616-deficient mutant and its isogenic complemented counterparts, we found that the infectivity potential of the mutant was significantly attenuated. The inactivation of bb0616 displayed no effect on borrelial growth in the medium or resistance to oxidative stress, but the mutant was significantly more susceptible to osmotic stress. In addition, the production of global virulence regulators such as BosR and RpoS as well as virulence-associated outer surface lipoproteins OspC and DbpA was reduced in the mutant. These phenotypes were fully restored when gene mutation was complemented with a wild-type copy of bb0616. Based on these findings, we concluded that the hypothetical protein BB0616 is required for the optimal infectivity of B. burgdorferi, potentially by impacting B. burgdorferi virulence gene expression as well as survival of the spirochete under stressful conditions.

Identification and analysis of MAPK cascade gene families of Camellia oleifera and their roles in response to cold stress

BACKGROUND

Low-temperature severely limits the growth and development of Camellia oleifera (C. oleifera). The mitogen-activated protein kinase (MAPK) cascade plays a key role in the response to cold stress.

METHODS AND RESULTS

Our study aims to identify MAPK cascade genes in C. oleifera and reveal their roles in response to cold stress. In our study, we systematically identified and analyzed the MAPK cascade gene families of C. oleifera, including their physical and chemical properties, conserved motifs, and multiple sequence alignments. In addition, we characterized the interacting networks of MAPKK kinase (MAPKKK)-MAPK kinase (MAPKK)-MAPK in C. oleifera. The molecular mechanism of cold stress resistance of MAPK cascade genes in wild C. oleifera was analyzed by differential gene expression and real-time quantitative reverse transcription-PCR (qRT-PCR).

CONCLUSION

In this study, 21 MAPKs, 4 MAPKKs and 55 MAPKKKs genes were identified in the leaf transcriptome of C. oleifera. According to the phylogenetic results, MAPKs were divided into 4 groups (A, B, C and D), MAPKKs were divided into 3 groups (A, B and D), and MAPKKKs were divided into 2 groups (MEKK and Raf). Motif analysis showed that the motifs in each subfamily were conserved, and most of the motifs in the same subfamily were basically the same. The protein interaction network based on Arabidopsis thaliana (A. thaliana) homologs revealed that MAPK, MAPKK, and MAPKKK genes were widely involved in C. oleifera growth and development and in responses to biotic and abiotic stresses. Gene expression analysis revealed that the CoMAPKKK5/CoMAPKKK43/CoMAPKKK49-CoMAPKK4-CoMAPK8 module may play a key role in the cold stress resistance of wild C. oleifera at a high-elevation site in Lu Mountain (LSG). This study can facilitate the mining and utilization of genetic resources of C. oleifera with low-temperature tolerance.

Identification of Phytophthora cinnamomi CRN effectors and their roles in manipulating cell death during Persea americana infection

The oomycete Phytophthora cinnamomi is a devastating plant pathogen with a notably broad host range. It is the causal agent of Phytophthora root rot (PRR), arguably the most economically important yield-limiting disease in Persea americana (avocado). Despite this, our understanding of the mechanisms P. cinnamomi employs to infect and successfully colonize avocado remains limited, particularly regarding the pathogen's ability to maintain its biotrophic and necrotrophic lifestyles during infection. The pathogen utilises a large repertoire of effector proteins which function in facilitating and establishing disease in susceptible host plants. Crinkling and necrosis effectors (CRN/Crinklers) are suspected to manipulate cell death to aid in maintenance of the pathogens biotrophic and necrotrophic lifestyles during different stages of infection. The current study identified 25 P. cinnamomi CRN effectors from the GKB4 genome using an HMM profile and assigned putative function to them as either cell death inducers or suppressors. Function was assigned to 10 PcinCRNs by analysing their RNA-seq expression profiles, relatedness to other functionally characterised Phytophthora CRNs and tertiary protein predictions. The full-length coding sequences for these PcinCRNs were confirmed by Sanger sequencing, six of which were found to have two divergent alleles. The presence of alleles indicates that the proteins encoded may perform contradicting functions in cell death manipulation, or function in different host plant species. Overall, this study provides a foundation for future research on P. cinnamomi infection and cell death manipulation mechanisms.

Asymmetric evolution of ISG15 homologs and the immune adaptation to LBUSV infection in spotted seabass (Lateolabrax maculatus)

The battle between host and viral is ubiquitous across all ecosystems. Despite this, research is scarce on the antiviral characteristics of fish, particularly in those that primarily rely on innate immune responses. This study, comprehensively explored the genetic and antiviral features of ISG15 in spotted seabass, focusing on its response to largemouth bass ulcerative syndrome virus (LBUSV). Through whole-genome BLAST and PCR cloning, two ISG15 homologs, namely LmISG15a and LmISG15b, were identified in spotted seabass, both encoding highly conserved proteins. However, a distinctive contrast emerged in their expression patterns, with LmISG15a exhibiting high expression in immune organs while LmISG15b remained largely silent across various organs. Regulatory elements analysis indicated an asymmetric evolution of the two ISG15s, with the minimal expression of LmISG15b may attribute to the loss of a necessary ISRE and an additional instability "ATTTA" motif. Association analysis demonstrated a significant correlation between LmISG15a expression and LBUSV infection. Subsequent antiviral activity detection revealed that LmISG15a interacted with LBUSV, inhibiting its replication by activating ISGylation and downstream pro-inflammatory mediators. In summary, this study unveils a distinct evolutionary strategy of fish antiviral gene ISG15 and delineates its kinetic characteristics in response to LBUSV infection.

VARista: a free web platform for streamlined whole-genome variant analysis across T2T, hg38, and hg19

With the increasing importance of genomic data in understanding genetic diseases, there is an essential need for efficient and user-friendly tools that simplify variant analysis. Although multiple tools exist, many present barriers such as steep learning curves, limited reference genome compatibility, or costs. We developed VARista, a free web-based tool, to address these challenges and provide a streamlined solution for researchers, particularly those focusing on rare monogenic diseases. VARista offers a user-centric interface that eliminates much of the technical complexity typically associated with variant analysis. The tool directly supports VCF files generated using reference genomes hg19, hg38, and the emerging T2T, with seamless remapping capabilities between them. Features such as gene summaries and links, tissue and cell-specific gene expression data for both adults and fetuses, as well as automated PCR design and integration with tools such as SpliceAI and AlphaMissense, enable users to focus on the biology and the case itself. As we demonstrate, VARista proved effective in narrowing down potential disease-causing variants, prioritizing them effectively, and providing meaningful biological context, facilitating rapid decision-making. VARista stands out as a freely available and comprehensive tool that consolidates various aspects of variant analysis into a single platform that embraces the forefront of genomic advancements. Its design inherently supports a shift in focus from technicalities to critical thinking, thereby promoting better-informed decisions in genetic disease research. Given its unique capabilities and user-centric design, VARista has the potential to become an essential asset for the genomic research community. https://VARista.link.

Characterization and Hydrolysis Mechanism Analysis of a Cold-Adapted Trypsin-Like Protease from Antarctic Krill

Cold-adapted proteases are capable of efficient protein hydrolysis at reduced temperatures, which offer significant potential applications in the area of low temperature food processing. In this paper, we attempted to characterize cold-adapted proteases from Antarctic krill. Antarctic krill possesses an extremely active autolytic enzyme system in their bodies, and the production of peptides and free amino acids accompanies the rapid breakdown of muscle proteins following the death. The crucial role of trypsin in this process is recognized. A cold-adapted trypsin named OUC-Pp-20 from Antarctic krill genome was cloned and expressed in Pichia pastoris. Recombinant trypsin is a monomeric protein of 26.8 ± 1.0 kDa with optimum reaction temperature at 25 °C. In addition, the catalytic specificity of OUC-Pp-20 was assessed by identifying its hydrolysis sites through LC-MS/MS. OUC-Pp-20 appeared to prefer Gln and Asn at the P1 position, which is an amino acid with an amide group in its side chain. Hydrolysis reactions on milk and shrimp meat revealed that it can effectively degrade allergenic components in milk and arginine kinase in shrimp meat. These findings update the current knowledge of cold-adapted trypsin and demonstrate the potential application of OUC-Pp-20 in low temperature food processing.

A chromosome-scale genome of Rhus chinensis Mill. provides new insights into plant-insect interaction and gallotannins biosynthesis

Rhus chinensis Mill., an economically valuable Anacardiaceae species, is parasitized by the galling aphid Schlechtendalia chinensis, resulting in the formation of the Chinese gallnut (CG). Here, we report a chromosomal-level genome assembly of R. chinensis, with a total size of 389.40 Mb and scaffold N50 of 23.02 Mb. Comparative genomic and transcriptome analysis revealed that the enhanced structure of CG and nutritional metabolism contribute to improving the adaptability of R. chinensis to S. chinensis by supporting CG and galling aphid growth. CG was observed to be abundant in hydrolysable tannins (HT), particularly gallotannin and its isomers. Tandem repeat clusters of dehydroquinate dehydratase/shikimate dehydrogenase (DQD/SDH) and serine carboxypeptidase-like (SCPL) and their homologs involved in HT production were determined as specific to HT-rich species. The functional differentiation of DQD/SDH tandem duplicate genes and the significant contraction in the phenylalanine ammonia-lyase (PAL) gene family contributed to the accumulation of gallic acid and HT while minimizing the production of shikimic acid, flavonoids, and condensed tannins in CG. Furthermore, we identified one UDP glucosyltransferase (UGT84A), three carboxylesterase (CXE), and six SCPL genes from conserved tandem repeat clusters that are involved in gallotannin biosynthesis and hydrolysis in CG. We then constructed a regulatory network of these genes based on co-expression and transcription factor motif analysis. Our findings provide a genomic resource for the exploration of the underlying mechanisms of plant-galling insect interaction and highlight the importance of the functional divergence of tandem duplicate genes in the accumulation of secondary metabolites.

Two maize homologs of mammalian proton-coupled folate transporter, ZmMFS_1-62 and ZmMFS_1-73, are essential to salt and drought tolerance

Folates are essential to the maintenance of normal life activities in almost all organisms. Proton-coupled folate transporter (PCFT), belonging to the major facilitator superfamily, is one of the three major folate transporter types widely studied in mammals. However, information about plant PCFTs is limited. Here, a genome-wide identification of maize PCFTs was performed, and two PCFTs, ZmMFS_1-62 and ZmMFS_1-73, were functionally investigated. Both proteins contained the typical 12 transmembrane helixes with N- and C-termini located in the cytoplasm, and were localized in the plasma membrane. Molecular docking analysis indicated their binding activity with folates via hydrogen bonding. Interference with ZmMFS_1-62 and ZmMFS_1-73 in maize seedlings through virus-induced gene silencing disrupted folate homeostasis, mainly in the roots, and reduced tolerance to drought and salt stresses. Moreover, a molecular chaperone protein, ZmHSP20, was found to interact with ZmMFS_1-62 and ZmMFS_1-73, and interference with ZmHSP20 in maize seedlings also led to folate disruption and increased sensitivity to drought and salt stresses. Overall, this is the first report of functional identification of maize PCFTs, which play essential roles in salt and drought stress tolerance, thereby linking folate metabolism with abiotic stress responses in maize.

Recurrent missense variant identified in two unrelated families with MPZL2-related hearing loss, expanding the variant spectrum associated with DFNB111

MPZL2-related hearing loss is a rare form of autosomal recessive hearing loss characterized by progressive, mild sloping to severe sensorineural hearing loss. Thirty-five previously reported patients had biallelic truncating variants in MPZL2, with the exception of one patient with a missense variant of uncertain significance and a truncating variant. Here, we describe the clinical characteristics and genotypes of five patients from four families with confirmed MPZL2-related hearing loss. A rare missense likely pathogenic variant [NM_005797.4(MPZL2):c.280C>T,p.(Arg94Trp)] located in exon 3 was confirmed to be in trans with a recurrent pathogenic truncating variant that segregated with hearing loss in three of the patients from two unrelated families. This is the first recurrent likely pathogenic missense variant identified in MPZL2. Apparently milder or later-onset hearing loss associated with rare missense variants in MPZL2 indicates that some missense variants in this gene may cause a milder phenotype than that resulting from homozygous or compound heterozygous truncating variants. This study, along with the identification of truncating loss of function and missense MPZL2 variants in several diverse populations, suggests that MPZL2-related hearing loss may be more common than previously appreciated and demonstrates the need for MPZL2 inclusion in hearing loss testing panels.

Genome-wide identification and characterization of the sweet orange (Citrus sinensis) basic helix-loop-helix (bHLH) family reveals a role for CsbHLH085 as a regulator of citrus bacterial canker resistance

Citrus bacterial canker (CBC) is a harmful bacterial disease caused by Xanthomonas citri subsp. citri (Xcc), negatively impacting citrus production worldwide. The basic helix-loop-helix (bHLH) transcription factor family plays crucial roles in plant development and stress responses. This study aimed to identify and annotate bHLH proteins encoded in the Citrus sinensis genome and explore their involvement and functional importance in regulating CBC resistance. A total of 135 putative CsbHLHs TFs were identified and categorized into 16 subfamilies. Their chromosomal locations, collinearity, and phylogenetic relationships were comprehensively analyzed. Upon Xcc strain YN1 infection, certain CsbHLHs were differentially regulated in CBC-resistant and CBC-sensitive citrus varieties. Among these, CsbHLH085 was selected for further functional characterization. CsbHLH085 was upregulated in the CBC-resistant citrus variety, was localized in the nucleus, and had a transcriptional activation activity. CsbHLH085 overexpression in Citrus significantly enhanced CBC resistance, accompanied by increased levels of salicylic acid (SA), jasmonic acid (JA), reactive oxygen species (ROS), and decreased levels of abscisic acid (ABA) and antioxidant enzymes. Conversely, CsbHLH085 virus-induced gene silencing resulted in opposite phenotypic and biochemical responses. CsbHLH085 silencing also affected the expression of phytohormone biosynthesis and signaling genes involved in SA, JA, and ABA signaling. These findings highlight the crucial role of CsbHLH085 in regulating CBC resistance, suggesting its potential as a target for biotechnological-assisted breeding citrus varieties with improved resistance against phytopathogens.

Two aminopeptidase I homologs convergently contribute to pathobiology of fungal entomopathogen Beauveria bassiana via divergent physiology-dependent autophagy pathways for vacuolar targeting

INTRODUCTION

In yeast, the cytoplasm-to-vacuole targeting (Cvt) pathway acts as a biosynthetic autophagy-related process, in which vacuolar targeting of hydrolase is mediated by the machineries involved in the selective autophagy. However, the mechanistic insights into vacuolar targeting of hydrolases through the selective autophagy pathway still remain enigmatic in filamentous fungi.

OBJECTIVES

Our study aims to investigate the mechanisms involved in vacuolar targeting of hydrolases in filamentous fungi.

METHODS

The filamentous entomopathogenic fungus Beauveria bassiana was used as a representative of filamentous fungi. We identified the homologs of yeast aminopeptidase I (Ape1) in B. bassiana by bioinformatic analyses and characterized their physiological roles by gene function analyses. Pathways for vacuolar targeting of hydrolases were investigated via molecular trafficking analyses.

RESULTS

B. bassiana has two homologs of yeast aminopeptidase I (Ape1) which are designated as BbApe1A and BbApe1B. The two homologs of yeast Ape1 contribute to starvation tolerance, development, and virulence in B. bassiana. Significantly, BbNbr1 acts as a selective autophagy receptor to mediate the vacuolar targeting of the two Ape1 proteins, in which BbApe1B interacts with BbNbr1 also directly interacting with BbAtg8, and BbApe1A has an additional requirement of the scaffold protein BbAtg11 that interacts with BbNbr1 and BbAtg8. Protein processing occurs at both terminuses of BbApe1A and only at carboxyl terminus of BbApe1B, which is also dependent on the autophagy-related proteins. Together, the functions and translocation processes of the two Ape1 proteins are associated with autophagy in fungal lifecycle.

CONCLUSION

This study reveals the functions and translocation processes for vacuolar hydrolases in the insect-pathogenic fungi and improves our understandings of the Nbr1-mediated vacuolar targeting pathway in the filamentous fungi.

Genome‑wide identification and expression analysis of the UBC gene family in wheat (Triticum aestivum L.)

BACKGROUND

Ubiquitination is an important regulatory step of selective protein degradation in the plant UPS (ubiquitin-proteasome system), which is involved in various biological processes in eukaryotes. Ubiquitin-conjugating enzymes play an intermediate role in the process of protein ubiquitination reactions and thus play an essential role in regulating plant growth and response to adverse environmental conditions. However, a genome-wide analysis of the UBC gene family in wheat (Triticum aestivum L.) has not yet been performed.

RESULTS

In this study, the number, physiochemical properties, gene structure, collinearity, and phylogenetic relationships of TaUBC family members in wheat were analyzed using bioinformatics methods. The expression pattern of TaUBC genes in different tissues/organs and developmental periods, as well as the transcript levels under abiotic stress treatment, were analyzed using RNA-Seq data and qRT-PCR. Meanwhile, favorable haplotypes of TaUBC25 were investigated based on wheat resequencing data of 681 wheat cultivars from the Wheat Union Database. The analyses identified a total of 93 TaUBC family members containing a UBC domain in wheat genome. These genes were unevenly distributed across 21 chromosomes, and numerous duplication events were observed between gene members. Based on phylogenetic analysis, the TaUBC family was divided into 13 E2 groups and a separate UEV group. We investigated the expression of TaUBC family genes under different tissue/organ and stress conditions by quantitative real-time PCR (qRT-PCR) analysis. The results showed that some TaUBC genes were specifically expressed in certain tissues/organs and that most TaUBC genes responded to NaCl, PEG6000, and ABA treatment with different levels of expression. In addition, we performed association analysis for the two haplotypes based on key agronomic traits such as thousand-kernel weight (TKW), kernel length (KL), kernel weight (KW), and kernel thickness (KT), examining 122 wheat accessions at three environmental sites. The results showed that TaUBC25-Hap II had significantly higher TKW, KL, KW, and KT than TaUBC25-Hap I. The distribution analysis of haplotypes showed that TaUBC25-Hap II was preferred in the natural population of wheat.

CONCLUSION

Our results identified 93 members of the TaUBC family in wheat, and several genes involved in grain development and abiotic stress response. Based on the SNPs detected in the TaUBC sequence, two haplotypes, TaUBC25-Hap I and TaUBC25-Hap II, were identified among wheat cultivars, and their potential value for wheat breeding was validated by association analysis. The above results provide a theoretical basis for elucidating the evolutionary relationships of the TaUBC gene family and lay the foundation for studying the functions of family members in the future.

Trophic Position of the White Worm (Enchytraeus albidus) in the Context of Digestive Enzyme Genes Revealed by Transcriptomics Analysis

To assess the impact of Enchytraeidae (potworms) on the functioning of the decomposer system, knowledge of the feeding preferences of enchytraeid species is required. Different food preferences can be explained by variations in enzymatic activities among different enchytraeid species, as there are no significant differences in the morphology or anatomy of their alimentary tracts. However, it is crucial to distinguish between the contribution of microbial enzymes and the animal's digestive capacity. Here, we computationally analyzed the endogenous digestive enzyme genes in Enchytraeus albidus. The analysis was based on RNA-Seq of COI-monohaplotype culture (PL-A strain) specimens, utilizing transcriptome profiling to determine the trophic position of the species. We also corroborated the results obtained using transcriptomics data from genetically heterogeneous freeze-tolerant strains. Our results revealed that E. albidus expresses a wide range of glycosidases, including GH9 cellulases and a specific digestive SH3b-domain-containing i-type lysozyme, previously described in the earthworm Eisenia andrei. Therefore, E. albidus combines traits of both primary decomposers (primary saprophytophages) and secondary decomposers (sapro-microphytophages/microbivores) and can be defined as an intermediate decomposer. Based on assemblies of publicly available RNA-Seq reads, we found close homologs for these cellulases and i-type lysozymes in various clitellate taxa, including Crassiclitellata and Enchytraeidae.

The repertoire and structure of adhesion GPCR transcript variants assembled from publicly available deep-sequenced human samples

Alternative splicing and multiple transcription start and termination sites can produce a diverse repertoire of mRNA transcript variants from a given gene. While the full picture of the human transcriptome is still incomplete, publicly available RNA datasets have enabled the assembly of transcripts. Using publicly available deep sequencing data from 927 human samples across 48 tissues, we quantified known and new transcript variants, provide an interactive, browser-based application Splice-O-Mat and demonstrate its relevance using adhesion G protein-coupled receptors (aGPCRs) as an example. On average, 24 different transcript variants were detected for each of the 33 human aGPCR genes, and several dominant transcript variants were not yet annotated. Variable transcription starts and complex exon-intron structures encode a flexible protein domain architecture of the N- and C termini and the seven-transmembrane helix domain (7TMD). Notably, we discovered the first GPCR (ADGRG7/GPR128) with eight transmembrane helices. Both the N- and C terminus of this aGPCR were intracellularly oriented, anchoring the N terminus in the plasma membrane. Moreover, the assessment of tissue-specific transcript variants, also for other gene classes, in our application may change the evaluation of disease-causing mutations, as their position in different transcript variants may explain tissue-specific phenotypes.

Genome-wide identification of the Hsp70 gene family in Penaeus chinensis and their response to environmental stress

The heat shock protein 70 (HSP70) gene family plays a crucial role in the response of organisms to environmental stress. However, it has not been systematically characterized in shrimp. In this study, we identified 25 PcHsp70 genes in the Penaeus chinensis genome. The encoded proteins were categorized into six subgroups based on phylogenetic relationships. Tandem duplication was the main driver of amplification in the PcHsp70 family, and the genes have experienced strong purifying selection during evolution. Transcriptome data analysis revealed that the 25 PcHsp70 members have different expression patterns in shrimp under conditions of low temperature, low salinity, and white spot syndrome virus infection. Among them, PcHsp70.11 was significantly induced under all three stress conditions, suggesting that this gene plays an important role in response to environmental stress in P. chinensis. To the best of our knowledge, this is the first study to systematically analyze the Hsp70 gene family in shrimp. The results provide important information on shrimp Hsp70s, contributing to a better understanding of the role of these genes in environmental stress and providing a basis for further functional studies.

Investigating the mechanism of rough phenotype in a naturally attenuated Brucella strain: insights from whole genome sequencing

Objective

Brucellosis, a significant zoonotic disease, not only impacts animal health but also profoundly influences the host immune responses through gut microbiome. Our research focuses on whole genome sequencing and comparative genomic analysis of these Brucella strains to understand the mechanisms of their virulence changes that may deepen our comprehension of the host immune dysregulation.

Methods

The Brucella melitensis strain CMCC55210 and its naturally attenuated variant CMCC55210a were used as models. Biochemical identification tests and in vivo experiments in mice verified the characteristics of the strain. To understand the mechanism of attenuation, we then performed de novo sequencing of these two strains.

Results

We discovered notable genomic differences between the two strains, with a key single nucleotide polymorphism (SNP) mutation in the manB gene potentially altering lipopolysaccharide (LPS) structure and influencing host immunity to the pathogen. This mutation might contribute to the attenuated strain's altered impact on the host's macrophage immune response, overing insights into the mechanisms of immune dysregulation linked to intracellular survival. Furthermore, we explore that manipulating the Type I restriction-modification system in Brucella can significantly impact its genome stability with the DNA damage response, consequently affecting the host's immune system.

Conclusion

This study not only contributes to understanding the complex relationship between pathogens, and the immune system but also opens avenues for innovative therapeutic interventions in inflammatory diseases driven by microbial and immune dysregulation.

Methyltransferase Domain-Focused Genome Mining for Fungal Polyketide Synthases

A comparison of substrate-binding site amino acid residues in the C-methyltransferase (MT) domains of fungal nonreducing polyketide synthases (NR-PKSs) suggests that these residues are correlated with the methylation modes used by the PKSs. A PKS, designated as AsbPKS, with substrate-binding site residues distinct from those of other known PKSs is focused on. The characterization of AsbPKS revealed that it yields an isocoumarin derivative, anhydrosclerotinin B (1), the biosynthesis of which involves a previously unreported methylation pattern. This study demonstrates the utility of MT domain-focused genome mining for the discovery of PKSs with new functions.

Genome-Wide Characterization and Functional Validation of the ACS Gene Family in the Chestnut Reveals Its Regulatory Role in Ovule Development

Ovule abortion significantly contributes to a reduction in chestnut yield. Therefore, an examination of the mechanisms underlying ovule abortion is crucial for increasing chestnut yield. In our previous study, we conducted a comprehensive multiomic analysis of fertile and abortive ovules and found that ACS genes in chestnuts (CmACS) play a crucial role in ovule development. Therefore, to further study the function of ACS genes, a total of seven CmACS members were identified, their gene structures, conserved structural domains, evolutionary trees, chromosomal localization, and promoter cis-acting elements were analyzed, and their subcellular localization was predicted and verified. The spatiotemporal specificity of the expression of the seven CmACS genes was confirmed via qRT-PCR analysis. Notably, CmACS7 was exclusively expressed in the floral organs, and its expression peaked during fertilization and decreased after fertilization. The ACC levels remained consistently greater in fertile ovules than in abortive ovules. The ACSase activity of CmACS7 was identified using the genetic transformation of chestnut healing tissue. Micro Solanum lycopersicum plants overexpressing CmACS7 had a significantly greater rate of seed failure than did wild-type plants. Our results suggest that ovule fertilization activates CmACS7 and increases ACC levels, whereas an overexpression of CmACS7 leads to an increase in ACC content in the ovule prior to fertilization, which can lead to abortion. In conclusion, the present study demonstrated that chestnut ovule abortion is caused by poor fertilization and not by nutritional competition. Optimization of the pollination and fertilization of female flowers is essential for increasing chestnut yield and reducing ovule abortion.

Genome-wide identification and stress response analysis of BcaCPK gene family in amphidiploid Brassica carinata

BACKGROUND

Calcium-dependent protein kinases (CPKs) are crucial for recognizing and transmitting Ca2+ signals in plant cells, playing a vital role in growth, development, and stress response. This study aimed to identify and detect the potential roles of the CPK gene family in the amphidiploid Brassica carinata (BBCC, 2n = 34) using bioinformatics methods.

RESULTS

Based on the published genomic information of B. carinata, a total of 123 CPK genes were identified, comprising 70 CPK genes on the B subgenome and 53 on the C subgenome. To further investigate the homologous evolutionary relationship between B. carinata and other plants, the phylogenetic tree was constructed using CPKs in B. carinata and Arabidopsis thaliana. The phylogenetic analysis classified 123 family members into four subfamilies, where gene members within the same subfamily exhibited similar conserved motifs. Each BcaCPK member possesses a core protein kinase domain and four EF-hand domains. Most of the BcaCPK genes contain 5 to 8 introns, and these 123 BcaCPK genes are unevenly distributed across 17 chromosomes. Among these BcaCPK genes, 120 replicated gene pairs were found, whereas only 8 genes were tandem duplication, suggesting that dispersed duplication mainly drove the family amplification. The results of the Ka/Ks analysis indicated that the CPK gene family of B. carinata was primarily underwent purification selection in evolutionary selection. The promoter region of most BcaCPK genes contained various stress-related cis-acting elements. qRT-PCR analysis of 12 selected CPK genes conducted under cadmium and salt stress at various points revealed distinct expression patterns among different family members in response to different stresses. Specifically, the expression levels of BcaCPK2.B01a, BcaCPK16.B02b, and BcaCPK26.B02 were down-regulated under both stresses, whereas the expression levels of other members were significantly up-regulated under at least one stress.

CONCLUSION

This study systematically identified the BcaCPK gene family in B. carinata, which contributes to a better understanding the CPK genes in this species. The findings also serve as a reference for analyzing stress responses, particularly in relation to cadmium and salt stress in B. carinata.

Functional domains of Acinetobacter bacteriophage tail fibers

A rapid increase in antimicrobial resistant bacterial infections around the world is causing a global health crisis. The Gram-negative bacterium Acinetobacter baumannii is categorized as a Priority 1 pathogen for research and development of new antimicrobials by the World Health Organization due to its numerous intrinsic antibiotic resistance mechanisms and ability to quickly acquire new resistance determinants. Specialized phage enzymes, called depolymerases, degrade the bacterial capsule polysaccharide layer and show therapeutic potential by sensitizing the bacterium to phages, select antibiotics, and serum killing. The functional domains responsible for the capsule degradation activity are often found in the tail fibers of select A. baumannii phages. To further explore the functional domains associated with depolymerase activity, tail-associated proteins of 71 sequenced and fully characterized phages were identified from published literature and analyzed for functional domains using InterProScan. Multisequence alignments and phylogenetic analyses were conducted on the domain groups and assessed in the context of noted halo formation or depolymerase characterization. Proteins derived from phages noted to have halo formation or a functional depolymerase, but no functional domain hits, were modeled with AlphaFold2 Multimer, and compared to other protein models using the DALI server. The domains associated with depolymerase function were pectin lyase-like (SSF51126), tailspike binding (cd20481), (Trans)glycosidases (SSF51445), and potentially SGNH hydrolases. These findings expand our knowledge on phage depolymerases, enabling researchers to better exploit these enzymes for therapeutic use in combating the antimicrobial resistance crisis.

Comprehensive virulence profiling and evolutionary analysis of specificity determinants in Staphylococcus aureus two-component systems

In the Staphylococcus aureus genome, a set of highly conserved two-component systems (TCSs) composed of histidine kinases (HKs) and their cognate response regulators (RRs) sense and respond to environmental stimuli, which drive the adaptation of the bacteria. This study investigates the complex interplay between TCSs in S. aureus USA300, a predominant methicillin-resistant S. aureus strain, revealing shared and unique virulence regulatory pathways and genetic variations mediating signal specificity within TCSs. Using TCS-related mutants from the Nebraska Transposon Mutant Library, we analyzed the effects of inactivated TCS HKs and RRs on the production of various virulence factors, in vitro infection abilities, and adhesion assays. We found that the TCSs' influence on virulence determinants was not associated with their phylogenetic relationship, indicating divergent functional evolution. Using the co-crystallized structure of the DesK-DesR from Bacillus subtilis and the modeled structures of the four NarL TCSs in S. aureus, we identified interacting residues, revealing specificity determinants and conservation within the same TCS, even from different strain backgrounds. The interacting residues were highly conserved within strains but varied between species due to selection pressures and the coevolution of cognate pairs. This study unveils the complex interplay and divergent functional evolution of TCSs, highlighting their potential for future experimental exploration of phosphotransfer between cognate and non-cognate recombinant HK and RRs.IMPORTANCEGiven the widespread conservation of two-component systems (TCSs) in bacteria and their pivotal role in regulating metabolic and virulence pathways, they present a compelling target for anti-microbial agents, especially in the face of rising multi-drug-resistant infections. Harnessing TCSs therapeutically necessitates a profound understanding of their evolutionary trajectory in signal transduction, as this underlies their unique or shared virulence regulatory pathways. Such insights are critical for effectively targeting TCS components, ensuring an optimized impact on bacterial virulence, and mitigating the risk of resistance emergence via the evolution of alternative pathways. Our research offers an in-depth exploration of virulence determinants controlled by TCSs in S. aureus, shedding light on the evolving specificity determinants that orchestrate interactions between their cognate pairs.

Lineage-specific gene duplication and expansion of DUF1216 gene family in Brassicaceae

Proteins containing domain of unknown function (DUF) are prevalent in eukaryotic genome. The DUF1216 proteins possess a conserved DUF1216 domain resembling to the mediator protein of Arabidopsis RNA polymerase II transcriptional subunit-like protein. The DUF1216 family are specifically existed in Brassicaceae, however, no comprehensive evolutionary analysis of DUF1216 genes have been performed. We performed a first comprehensive genome-wide analysis of DUF1216 proteins in Brassicaceae. Totally 284 DUF1216 genes were identified in 27 Brassicaceae species and classified into four subfamilies on the basis of phylogenetic analysis. The analysis of gene structure and conserved motifs revealed that DUF1216 genes within the same subfamily exhibited similar intron/exon patterns and motif composition. The majority members of DUF1216 genes contain a signal peptide in the N-terminal, and the ninth position of the signal peptide in most DUF1216 is cysteine. Synteny analysis revealed that segmental duplication is a major mechanism for expanding of DUF1216 genes in Brassica oleracea, Brassica juncea, Brassica napus, Lepidium meyneii, and Brassica carinata, while in Arabidopsis thaliana and Capsella rubella, tandem duplication plays a major role in the expansion of the DUF1216 gene family. The analysis of Ka/Ks (non-synonymous substitution rate/synonymous substitution rate) ratios for DUF1216 paralogous indicated that most of gene pairs underwent purifying selection. DUF1216 genes displayed a specifically high expression in reproductive tissues in most Brassicaceae species, while its expression in Brassica juncea was specifically high in root. Our studies offered new insights into the phylogenetic relationships, gene structures and expressional patterns of DUF1216 members in Brassicaceae, which provides a foundation for future functional analysis.

Disruption of Zfh3 abolishes mulberry-specific monophagy in silkworm larvae

Feeding behavior is critical for insect survival and fitness. Most researchers have explored the molecular basis of feeding behaviors by identifying and elucidating the function of olfactory receptors (ORs) and gustatory receptors (GRs). Other types of genes, such as transcription factors, have rarely been investigated, and little is known about their potential roles. The silkworm (Bombyx mori) is a well-studied monophagic insect which primarily feeds on mulberry leaves, but the genetic basis of its monophagy is still not understood. In this report, we focused on a transcription factor encoded by the Zfh3 gene, which is highly expressed in the silkworm central and peripheral nervous systems, including brain, antenna, and maxilla. To investigate its function, Zfh3 was abrogated using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR associated protein 9 (Cas9) mutagenesis. Since Zfh3 knockout homozygotes are not viable, we studied feeding behavior in heterozygotes, and found that disruption of Zfh3 affects both gustation and olfaction. Mutant larvae lose preference for mulberry leaves, acquire the ability to consume an expanded range of diets, and exhibit improved adaptation to the M0 artificial diet, which contains no mulberry leaves. These results provide the first demonstration that a transcription factor modulates feeding behaviors in an insect.

Intermembrane space-localized TbTim15 is an essential subunit of the single mitochondrial inner membrane protein translocase of trypanosomes

All mitochondria import >95% of their proteins from the cytosol. This process is mediated by protein translocases in the mitochondrial membranes, whose subunits are generally highly conserved. Most eukaryotes have two inner membrane protein translocases (TIMs) that are specialized to import either presequence-containing or mitochondrial carrier proteins. In contrast, the parasitic protozoan Trypanosoma brucei has a single TIM complex consisting of one conserved and five unique subunits. Here, we identify candidates for new subunits of the TIM or the presequence translocase-associated motor (PAM) using a protein-protein interaction network of previously characterized TIM and PAM subunits. This analysis reveals that the trypanosomal TIM complex contains an additional trypanosomatid-specific subunit, designated TbTim15. TbTim15 is associated with the TIM complex, lacks transmembrane domains, and localizes to the intermembrane space. TbTim15 is essential for procyclic and bloodstream forms of trypanosomes. It contains two twin CX9C motifs and mediates import of both presequence-containing and mitochondrial carrier proteins. While the precise function of TbTim15 in mitochondrial protein import is unknown, our results are consistent with the notion that it may function as an import receptor for the non-canonical trypanosomal TIM complex.

Global proteome of the saprophytic strain Leptospira biflexa and comparative analysis with pathogenic strain Leptospira interrogans uncover new pathogenesis mechanisms

Leptospira is a genus of bacteria that includes free-living saprophytic species found in water or soil, and pathogenic species, which are the etiologic agents of leptospirosis. Besides all the efforts, there are only a few proteins described as virulence factors in the pathogenic strain L. interrogans. This work aims to perform L. biflexa serovar Patoc1 strain Paris global proteome and to compare with the proteome database of pathogenic L. interrogans serovar Copenhageni strain Fiocruz L1-130. We identified a total of 2327 expressed proteins of L. biflexa by mass spectrometry. Using the Get Homologues software with the global proteome of L. biflexa and L. interrogans, we found orthologous proteins classified into conserved, low conserved, and specific proteins. Comparative bioinformatic analyses were performed to understand the biological functions of the proteins, subcellular localization, the presence of signal peptide, structural domains, and motifs using public softwares. These results lead to the selection of 182 low conserved within the saprophyte, and 176 specific proteins of L. interrogans. It is anticipated that these findings will indicate further studies to uncover virulence factors in the pathogenic strain. This work presents for the first time the global proteome of saprophytic strain L. biflexa serovar Patoc, strain Patoc1. SIGNIFICANCE: The comparative analysis established an array of specific proteins in pathogenic strain that will narrow down the identification of immune protective proteins that will help fight leptospirosis.

A bird's-eye view: exploration of the flavin-containing monooxygenase superfamily in common wheat

The Flavin Monooxygenase (FMO) gene superfamily in plants is involved in various processes most widely documented for its involvement in auxin biosynthesis, specialized metabolite biosynthesis, and plant microbial defense signaling. The roles of FMOs in defense signaling and disease resistance have recently come into focus as they may present opportunities to increase immune responses in plants including leading to systemic acquired resistance, but are not well characterized. We present a comprehensive catalogue of FMOs found in genomes across vascular plants and explore, in depth, 170 wheat TaFMO genes for sequence architecture, cis-acting regulatory elements, and changes due to Transposable Element insertions. A molecular phylogeny separates TaFMOs into three clades (A, B, and C) for which we further report gene duplication patterns, and differential rates of homoeologue expansion and retention among TaFMO subclades. We discuss Clade B TaFMOs where gene expansion is similarly seen in other cereal genomes. Transcriptome data from various studies point towards involvement of subclade B2 TaFMOs in disease responses against both biotrophic and necrotrophic pathogens, substantiated by promoter element analysis. We hypothesize that certain TaFMOs are responsive to both abiotic and biotic stresses, providing potential targets for enhancing disease resistance, plant yield and other important agronomic traits. Altogether, FMOs in wheat and other crop plants present an untapped resource to be exploited for improving the quality of crops.

Comprehensive Proteomics Analysis of the Hemolymph Composition of Sugar-Fed Aedes aegypti Female and Male Mosquitoes

In arthropods, hemolymph carries immune cells and solubilizes and transports nutrients, hormones, and other molecules that are involved in diverse physiological processes including immunity, metabolism, and reproduction. However, despite such physiological importance, little is known about its composition. We applied mass spectrometry-based label-free quantification approaches to study the proteome of hemolymph perfused from sugar-fed female and male Aedes aegypti mosquitoes. A total of 1403 proteins were identified, out of which 447 of them were predicted to be extracellular. In both sexes, almost half of these extracellular proteins were predicted to be involved in defense/immune response, and their relative abundances (based on their intensity-based absolute quantification, iBAQ) were 37.9 and 33.2%, respectively. Interestingly, among them, 102 serine proteases/serine protease-homologues were identified, with almost half of them containing CLIP regulatory domains. Moreover, proteins belonging to families classically described as chemoreceptors, such as odorant-binding proteins (OBPs) and chemosensory proteins (CSPs), were also highly abundant in the hemolymph of both sexes. Our data provide a comprehensive catalogue of A. aegypti hemolymph basal protein content, revealing numerous unexplored targets for future research on mosquito physiology and disease transmission. It also provides a reference for future studies on the effect of blood meal and infection on hemolymph composition.