ExPASy: The proteomics server for in-depth protein knowledge and analysis

Molecular and Functional Characterization of Pheromone Binding Protein 2 from Cyrtotrachelus buqueti (Coleoptera: Curculionidae)

Pheromone-binding proteins (PBPs) play important roles in binding and transporting sex pheromones. However, the PBP genes identified in coleopteran insects and their information sensing mechanism are largely unknown. Cyrtotrachelus buqueti (Coleoptera: Curculionidae) is a major insect pest of bamboo plantations. In this study, a novel PBP gene, CbuqPBP2, from C. buqueti was functionally characterized. CbuqPBP2 was more abundantly expressed in the antennae of both sexes than other body parts, and its expression level was significantly male-biased. Fluorescence competitive binding assays showed that CbuqPBP2 exhibited the strongest binding affinity to dibutyl phthalate (Ki = 6.32 μM), followed by styrene (Ki = 11.37 μM), among twelve C. buqueti volatiles. CbuqPBP2, on the other hand, showed high binding affinity to linalool (Ki = 10.55), the main volatile of host plant Neosinocalamus affinis. Furthermore, molecular docking also demonstrated the strong binding ability of CbuqPBP2 to dibutyl phthalate, styrene, and linalool, with binding energy values of -5.7, -6.6, and -6.0 kcal/mol, respectively, and hydrophobic interactions were the prevailing forces. The knockdown of CbuqPBP2 expression via RNA interference significantly reduced the electroantennography (EAG) responses of male adults to dibutyl phthalate and styrene. In conclusion, these results will be conducive to understanding the olfactory mechanisms of C. buqueti and promoting the development of novel strategies for controlling this insect pest.

Taxonomic Distribution and Molecular Evolution of Mytilectins

R-type lectins are a widespread group of sugar-binding proteins found in nearly all domains of life, characterized by the presence of a carbohydrate-binding domain that adopts a β-trefoil fold. Mytilectins represent a recently described subgroup of β-trefoil lectins, which have been functionally characterized in a few mussel species (Mollusca, Bivalvia) and display attractive properties, which may fuel the development of artificial lectins with different biotechnological applications. The detection of different paralogous genes in mussels, together with the description of orthologous sequences in brachiopods, supports the formal description of mytilectins as a gene family. However, to date, an investigation of the taxonomic distribution of these lectins and their molecular diversification and evolution was still lacking. Here, we provide a comprehensive overview of the evolutionary history of mytilectins, revealing an ancient monophyletic evolutionary origin and a very broad but highly discontinuous taxonomic distribution, ranging from heteroscleromorphan sponges to ophiuroid and crinoid echinoderms. Moreover, the overwhelming majority of mytilectins display a chimera-like architecture, which combines the β-trefoil carbohydrate recognition domain with a C-terminal pore-forming domain, suggesting that the simpler structure of most functionally characterized mytilectins derives from a secondary domain loss.

COVID-19 Genomic Surveillance in Bangui (Central African Republic) Reveals a Landscape of Circulating Variants Linked to Validated Antiviral Targets of SARS-CoV-2 Proteome

Since its outbreak, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) spread rapidly, causing the Coronavirus Disease 19 (COVID-19) pandemic. Even with the vaccines' administration, the virus continued to circulate due to inequal access to prevention and therapeutic measures in African countries. Information about COVID-19 in Africa has been limited and contradictory, and thus regional studies are important. On this premise, we conducted a genomic surveillance study about COVID-19 lineages circulating in Bangui, Central African Republic (CAR). We collected 2687 nasopharyngeal samples at four checkpoints in Bangui from 2 to 22 July 2021. Fifty-three samples tested positive for SARS-CoV-2, and viral genomes were sequenced to look for the presence of different viral strains. We performed phylogenetic analysis and described the lineage landscape of SARS-CoV-2 circulating in the CAR along 15 months of pandemics and in Africa during the study period, finding the Delta variant as the predominant Variant of Concern (VoC). The deduced aminoacidic sequences of structural and non-structural genes were determined and compared to reference and reported isolates from Africa. Despite the limited number of positive samples obtained, this study provides valuable information about COVID-19 evolution at the regional level and allows for a better understanding of SARS-CoV-2 circulation in the CAR.

Codon optimization of a gene encoding DNA polymerase from Pyrococcus furiosus and its expression in Escherichia coli

BACKGROUND

DNA polymerase is an essential component in PCR assay for DNA synthesis. Improving DNA polymerase with characteristics indispensable for a powerful assay is crucial because it can be used in wide-range applications. Derived from Pyrococcus furiosus, Pfu DNA polymerase (Pfu pol) is one of the excellent polymerases due to its high fidelity. Therefore, we aimed to develop Pfu pol from a synthetic gene with codon optimization to increase its protein yield in Escherichia coli.

RESULTS

Recombinant Pfu pol was successfully expressed and purified with a two-step purification process using nickel affinity chromatography, followed by anion exchange chromatography. Subsequently, the purified Pfu pol was confirmed by Western blot analysis, resulting in a molecular weight of approximately 90 kDa. In the final purification process, we successfully obtained a large amount of purified enzyme (26.8 mg/L). Furthermore, the purified Pfu pol showed its functionality and efficiency when tested for DNA amplification using the standard PCR.

CONCLUSIONS

Overall, a high-level expression of recombinant Pfu pol was achieved by employing our approach in the present study. In the future, our findings will be useful for studies on synthesizing recombinant DNA polymerase in E. coli expression system.

The role of WRKY transcription factors in exogenous potassium (K+) response to NaCl stress in Tamarix ramosissima

Introduction: Soil salinization poses a significant challenge to plant growth and vitality. Plants like Tamarix ramosissima Ledeb (T. ramosissima), which are halophytes, are often integrated into planting schemes tailored for saline environments. Yet, the role of WRKY transcription factors in T. ramosissima, especially under sodium chloride (NaCl) stress mitigated by exogenous K+ application, is not well-understood. This research endeavors to bridge this knowledge gap. Methods: Using Pfam protein domain prediction and physicochemical property analysis, we delved into the WRKY genes in T. ramosissima roots that are implicated in counteracting NaCl stress when aided by exogenous K+ applications. By observing shifts in the expression levels of WRKY genes annotated to the KEGG pathway under NaCl stress at 0, 48, and 168 h, we aimed to identify potential key WRKY genes. Results: We found that the expression of 56 WRKY genes in T. ramosissima roots responded to exogenous K+ application during NaCl stress at the indicated time points. Particularly, the expression levels of these genes were primarily upregulated within 168 h. From these, 10 WRKY genes were found to be relevant in the KEGG pathways. Moreover, six genes, namely Unigene0024962, Unigene0024963, Unigene0010090, Unigene0007135, Unigene0070215, and Unigene0077293, were annotated to the Plant-pathogen interaction pathway or the MAPK signaling pathway in plants. These genes exhibited dynamic expression regulation at 48 h with the application of exogenous K+ under NaCl stress. Discussion: Our research highlights that WRKY transcription factors can modulate the activation or inhibition of related genes during NaCl stress with the application of exogenous K+. This regulation enhances the plant's adaptability to saline environments and mitigates the damage induced by NaCl. These findings provide valuable gene resources for future salt-tolerant Tamarix breeding and expand our understanding of the molecular mechanisms of WRKY transcription factors in alleviating NaCl toxicity.

A Genome-Wide Identification and Expression Analysis of the Casparian Strip Membrane Domain Protein-like Gene Family in Pogostemon cablin in Response to p-HBA-Induced Continuous Cropping Obstacles

Casparian strip membrane domain protein-like (CASPL) genes are key genes for the formation and regulation of the Casparian strip and play an important role in plant abiotic stress. However, little research has focused on the members, characteristics, and biological functions of the patchouli PatCASPL gene family. In this study, 156 PatCASPL genes were identified at the whole-genome level. Subcellular localization predicted that 75.6% of PatCASPL proteins reside on the cell membrane. A phylogenetic analysis categorized PatCASPL genes into five subclusters alongside Arabidopsis CASPL genes. In a cis-acting element analysis, a total of 16 different cis-elements were identified, among which the photo-responsive element was the most common in the CASPL gene family. A transcriptome analysis showed that p-hydroxybenzoic acid, an allelopathic autotoxic substance, affected the expression pattern of PatCASPLs, including a total of 27 upregulated genes and 30 down-regulated genes, suggesting that these PatCASPLs may play an important role in the regulation of patchouli continuous cropping obstacles by affecting the formation and integrity of Casparian strip bands. These results provided a theoretical basis for exploring and verifying the function of the patchouli PatCASPL gene family and its role in continuous cropping obstacles.

Comprehensive characterization and expression analysis of enzymatic antioxidant gene families in passion fruit (Passiflora edulis)

Passion fruit, a valuable tropical fruit, faces climate-related growth challenges. Antioxidant enzymes are vital for both stress protection and growth regulation in plants. We first provided systemic analysis of enzymatic antioxidant gene families in passion fruit, identifying 90 members including 11 PeSODs, 45 PeAPXs, 8 PeCATs, 7 PeGPXs, 6 PeMDHARs, 8 PeDHARs, and 5 PeGRs. Gene members in each gene family with same subcellular localization showed closer phylogenetic relationship. Many antioxidant genes exhibited tissue- or developmental stage-specific expression patterns during floral and fruit development, with some widely expressed. Their co-expressed genes were linked to photosynthesis and energy metabolism, suggesting roles in protecting highly proliferating tissues from oxidative damage. Potential genes for enhancing temperature stress resistance were identified. The involvement of diverse regulatory factors including miRNAs, transcription factors, and CREs might contribute to the complex roles of antioxidant genes. This study informs future research on antioxidant genes and passion fruit breeding.

Retracing the Rapid Evolution of an Herbicide-Degrading Enzyme by Protein Engineering

The mechanisms underlying the rapid evolution of novel enzymatic activities from promiscuous side activities are poorly understood. Recently emerged enzymes catalyzing the catabolic degradation of xenobiotic substances that have been spread out into the environment during the last decades provide an exquisite opportunity to study these mechanisms. A prominent example is the herbicide atrazine (2-chloro-4-ethylamino-6-isopropylamino-1,3,5-triazine), which is degraded through a number of enzymatic reactions constituting the Atz pathway. Here, we analyzed the evolution of the hydroxyatrazine ethylaminohydrolase AtzB, a Zn(II)-dependent metalloenzyme that adopts the amidohydrolase fold and catalyzes the second step of the Atz pathway. We searched for promiscuous side activities of AtzB, which might point to the identity of its progenitor. These investigations revealed that AtzB has low promiscuous guanine deaminase activity. Furthermore, we found that the two closest AtzB homologues, which have not been functionally annotated up to now, are guanine deaminases with modest promiscuous hydroxyatrazine hydrolase activity. Based on sequence comparisons with the closest AtzB homologues, the guanine deaminase activity of AtzB could be increased by three orders of magnitude through the introduction of only four active site mutations. Interestingly, introducing the inverse four mutations into the AtzB homologues significantly enhanced their hydroxyatrazine hydrolase activity, and in one case is even equivalent to that of wild-type AtzB. Molecular dynamics simulations elucidated the structural and molecular basis for the mutation-induced activity changes. The example of AtzB highlights how novel enzymes with high catalytic proficiency can evolve from low promiscuous side activities by only few mutational events within a short period of time.

A Genome-Wide Analysis of Serine Protease Inhibitors in Cydia pomonella Provides Insights into Their Evolution and Expression Pattern

Serine protease inhibitors (serpins) appear to be ubiquitous in almost all living organisms, with a conserved structure and varying functions. Serpins can modulate immune responses by negatively regulating serine protease activities strictly and precisely. The codling moth, Cydia pomonella (L.), a major invasive pest in China, can cause serious economic losses. However, knowledge of serpin genes in this insect remain largely unknown. In this study, we performed a systematic analysis of the serpin genes in C. pomonella, obtaining 26 serpins from the C. pomonella genome. Subsequently, their sequence features, evolutionary relationship, and expression pattern were characterized. Comparative analysis revealed the evolution of a number of serpin genes in Lepidoptera. Importantly, the evolutionary relationship and putative roles of serpin genes in C. pomonella were revealed. Additionally, selective pressure analysis found amino acid sites with strong evidence of positive selection. Interestingly, the serpin1 gene possessed at least six splicing isoforms with distinct reactive-center loops, and these isoforms were experimentally validated. Furthermore, we observed a subclade expansion of serpins, and these genes showed high expression in multiple tissues, suggesting their important roles in C. pomonella. Overall, this study will enrich our knowledge of the immunity of C. pomonella and help to elucidate the role of serpins in the immune response.

Genome wide identification, structural characterization and phylogenetic analysis of High-Affinity potassium (HAK) ion transporters in common bean (Phaseolus vulgaris L.)

BACKGROUND

High-Affinity Potassium ions represent one of the most important and large group of potassium transporters. Although HAK genes have been studied in a variety of plant species, yet, remain unexplored in common bean.

RESULTS

In the current study, 20 HAK genes were identified in common bean genome. Super-family "K_trans" domain was found in all PvHAK genes. Signals for localization of PvHAK proteins were detected in cell membrane. Fifty three HAKs genes, across diverse plant species, were divided into 5 groups based on sequential homology. Twelve pairs of orthologs genes were found in various plant species. PvHAKs genes were distributed unequally on 7 chromosomes with maximum number (7) mapped on chromosome 2 while only 1 PvHAK found on each chromosome 1, 4, and 6. Tandem gene duplication was witnessed in 2 paralog pairs while 1 pair exhibited segmental gene duplication. Five groups were made in PvHAK gene family based on Phylogeny. Maximum PvHAKs (10) were detected in Group-V while group-II composed of only 1 PvHAK gene. Variation was witnessed in number and size of motifs, and structure of PvHAKs associated with different groups. Light and hormone responsive elements contributed 57 and 24% share, respectively, to cis regulatory elements. qRT-PCR based results revealed significant increase in expression of all 4 PvHAK genes under low-potassium stress.

CONCLUSION

The current study provides valuable information for further functional characterization and uncovering the molecular mechanism associated with Potassium transportation in plants.

Molecular docking, simulations of animal peptides against the envelope protein of Dengue virus

Peptides are biologically active, small molecules with high specificity in its mode of action that can be effective at nanomolar concentrations. Peptide-based antiviral medicines have already been licensed and used to treat human immunodeficiency virus (HIV), influenza virus and hepatitis C virus. So far, no peptide drug has been approved for antiviral treatment against Dengue virus, and many are under clinical trials. Therefore, developing a reasonable peptide against the Dengue virus Envelope protein structure will be a successful strategy for treating Dengue. Hence, we investigated protein-protein docking interactions between 215 peptides retrieved from the AVP database against the envelope protein using Cluspro and HADDOCK followed by the evaluation of their allegenicity, toxicity and physicochemical characteristics investigation. Further validation of the protein-peptide complexes was performed with Molecular dynamics simulations and Molecular Mechanics Poisson-Boltzmann surface area (MMPBSA) analysis on the hit inhibitors. This study revealed that Indolicidin (-75.026 ± 1.54 KJ/mol) and Human Neutrophil peptide-1 (-71.6551 ± 2.112 KJ/mol) shows higher negative ΔG binding implicating their relative stabilization in the protein-peptide interactions during 100 ns of dynamic simulations. Also, both the peptides exhibited desirable physicochemical properties and were nonallergenic. Hence, we further aim to test these peptides by in vitro and in vivo studies to confirm their efficacy against Dengue virus.Communicated by Ramaswamy H. Sarma.

Genome-Wide Identification and Characterization of CLAVATA3/EMBRYO SURROUNDING REGION (CLE) Gene Family in Foxtail Millet (Setaria italica L.)

The CLAVATA3/EMBRYO-SURROUNDING REGION (CLE) genes encode signaling peptides that play important roles in various developmental and physiological processes. However, the systematic identification and characterization of CLE genes in foxtail millet (Setaria italica L.) remain limited. In this study, we identified and characterized 41 SiCLE genes in the foxtail millet genome. These genes were distributed across nine chromosomes and classified into four groups, with five pairs resulting from gene duplication events. SiCLE genes within the same phylogenetic group shared similar gene structure and motif patterns, while 34 genes were found to be single-exon genes. All SiCLE peptides harbored the conserved C-terminal CLE domain, with highly conserved positions in the CLE core sequences shared among foxtail millet, Arabidopsis, rice, and maize. The SiCLE genes contained various cis-elements, including five plant hormone-responsive elements. Notably, 34 SiCLE genes possessed more than three types of phytohormone-responsive elements on their promoters. Comparative analysis revealed higher collinearity between CLE genes in maize and foxtail millet, which may be because they are both C4 plants. Tissue-specific expression patterns were observed, with genes within the same group exhibiting similar and specific expression profiles. SiCLE32 and SiCLE41, classified in Group D, displayed relatively high expression levels in all tissues except panicles. Most SiCLE genes exhibited low expression levels in young panicles, while SiCLE6, SiCLE24, SiCLE25, and SiCLE34 showed higher expression in young panicles, with SiCLE24 down-regulated during later panicle development. Greater numbers of SiCLE genes exhibited higher expression in roots, with SiCLE7, SiCLE22, and SiCLE36 showing the highest levels and SiCLE36 significantly down-regulated after abscisic acid (ABA) treatment. Following treatments with ABA, 6-benzylaminopurine (6-BA), and gibberellic acid 3 (GA3), most SiCLE genes displayed down-regulation followed by subsequent recovery, while jasmonic acid (JA) and indole-3-acetic acid (IAA) treatments led to upregulation at 30 min in leaves. Moreover, identical hormone treatments elicited different expression patterns of the same genes in leaves and stems. This comprehensive study enhances our understanding of the SiCLE gene family and provides a foundation for further investigations into the functions and evolution of SiCLE genes in foxtail millet.

Molecular modeling, docking and dynamics studies of fenugreek (Trigonella foenum-graecum) α-amylase

α-Amylase catalyses the hydrolysis of glucosidic bonds in polysaccharides such as starch, glycogen and their degradation products. In the present study, the three-dimensional structure of fenugreek (Trigonella foenum-graecum) α-amylase was determined using a homology modeling-based technique. The best predicted model was deposited in PMDB server with PMDB ID PM0084364. The phylogenetic tree was created using the UPGMA method with 8 homologous protein sequences, Trigonella foenum-graecum was utilized as the target protein. Alignment of the phylogenetic tree identified two primary functional groupings (A and B). α-Amylase from the target genome Trigonella foenum-graecum (Acc. No: GHNA01022531.1) was clustered with Medicago truncatula (Acc. No: XP003589186.1), Cicer arietinum (Acc. No: XP004499059.1), Cajanus cajan (Acc. No: XP020231823.1), Vigna angularis (Acc. No: NP001316768.1) and Vigna mungo (Acc. No: P17859.1), in group A cluster, while Hordeum vulgare (Acc. No: Q40015) and Oryza sativa (PDB ID: 3WN6) were in cluster B. The molecular dynamics simulations were performed to understand the molecular basis and mode of action of Trigonella foenum-graecum α-amylase. Additionally, a geometry-based molecular docking technique was used to evaluate potential binding interactions between the modeled structure of α-amylase and maltose. The results show that Trp228, Glu226, Arg199, His308, Tyr165, Asp309, Phe202 and Asp201 from Trigonella foenum-graecum α-amylase enzyme is involved in the binding to the substrate maltose. Our study provides a 3D model of Trigonella foenum-graecum α-amylase and aids in understanding the atomic level molecular underpinnings of the mechanism of α-amylase interaction with substrate maltose. Ca2+ are essential for the stability of domain B since they are connected to it. Ca2+ site ligands are Asp139, Glu130, Thr133, Asp135 and Gly131 residues. HIGHLIGHTSIn silico analysis, gene prediction of α-amylase was carried from Trigonella foenum-graecum.Analysis of the structure of α-amylase was carried out using homology modelling.Calcium binding sites and their interactions with α-amylase were visualised using BIOVIA DISCOVERY STUDIO 2019.The molecular interaction between Trigonella foenum-graecum α-amylase and maltose was studied in silico using a molecular docking-based method.To give the required simulation parameters, RMSD, RMSF, and Total Energy were calculated using BIOVIA DISCOVERY STUDIO 2019.[Figure: see text]Communicated by Ramaswamy H. Sarma.

Genome sequencing of cucumber mosaic virus (CMV) isolates infecting chilli and its interaction with host ferredoxin protein of different host for causing mosaic symptoms

Chilli (Capsicum annuum L.) is an important vegetable crop grown in the Indian sub-continent and is prone to viral infections under field conditions. During the field survey, leaf samples from chilli plants showing typical symptoms of disease caused by cucumber mosaic virus (CMV) such as mild mosaic, mottling and leaf distortion were collected. DAC-ELISA analysis confirmed the presence of CMV in 71 out of 100 samples, indicating its widespread prevalence in the region. Five CMV isolates, named Gu1, Gu2, BA, Ho, and Sal were mechanically inoculated onto cucumber and Nicotiana glutinosa plants to study their virulence. Inoculated plants expressed the characteristic symptoms of CMV such as chlorotic spots followed by mild mosaic and leaf distortion. Complete genomes of the five CMV isolates were amplified, cloned, and sequenced, revealing RNA1, RNA2, and RNA3 sequences with 3358, 3045, and 2220 nucleotides, respectively. Phylogenetic analysis classified the isolates as belonging to the CMV-IB subgroup, distinguishing them from subgroup IA and II CMV isolates. Recombination analysis showed intra and interspecific recombination in all the three RNA segments of these isolates. In silico protein-protein docking approach was used to decipher the mechanism behind the production of mosaic symptoms during the CMV-host interaction in 13 host plants. Analysis revealed that the production of mosaic symptoms could be due to the interaction between the coat protein (CP) of CMV and chloroplast ferredoxin proteins. Further, in silico prediction was validated in 13 host plants of CMV by mechanical sap inoculation. Twelve host plants produced systemic symptoms viz., chlorotic spot, chlorotic ringspot, chlorotic local lesion, mosaic and mild mosaic and one host plant, Solanum lycopersicum produced mosaic followed by shoestring symptoms.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03777-8.

Wild lime psyllid Leuronota fagarae Burckhardt (Hemiptera: Psylloidea) picorna-like virus full genome annotation and classification

Picorna-like viruses of the order Picornavirales are a poorly defined group of positive-sense, single-stranded RNA viruses that include numerous pathogens known to infect plants, animals, and insects. A new picorna-like viral species was isolated from the wild lime psyllid (WLP), Leuronota fagarae, in the state of Florida, USA, and labelled: Leuronota fagarae picorna-like virus isolate FL (LfPLV-FL). The virus was found to have homology to a picorna-like virus identified in the Asian Citrus Psyllid (ACP), Diaphorina citri, collected in the state of Florida. Computational analysis of RNA extracts from WLP adult heads identified a 10,006-nucleotide sequence encoding a 2,942 amino acid polyprotein with similar functional domain structure to polyproteins of both Dicistroviridae and Iflaviridae. Sequence comparisons of nucleic acid and amino acid translations of the conserved RNA-dependent RNA polymerase, along with the entire N-terminal nonstructural coding region, provided insight into an evolutionary relationship of LfPLV-FL to insect-infecting iflaviruses. Viruses belonging to the family Iflaviridae encode a polyprotein of around 3000 amino acids in length that is processed post-translationally to produce components necessary for replication. The classification of a novel picorna-like virus in L. fagarae, with evolutionary characteristics similar to picorna-like viruses infecting Bactericera cockerelli and D. citri, provides an opportunity to examine virus host specificity, as well as identify critical components of the virus' genome required for successful transmission, infection, and replication. This bioinformatic classification allows for further insight into a novel virus species, and aids in the research of a closely related virus of the invasive psyllid, D. citri, a major pest of Floridian citriculture. The potential use of viral pathogens as expression vectors to manage the spread D. citri is an area that requires additional research; however, it may bring forth an effective control strategy to reduce the transmission of Candidatus Liberibacter asiaticus (CLas), the causative agent of Huanglongbing (HLB).

A new member of the Nudiviridae from the Florida stone crab (Menippe mercenaria)

Menippe mercenaria, the Florida stone crab, supports an unconventional fishery across the southern USA and Caribbean that involves claw-removal and the return of de-clawed animals to the sea. We provide pathological, ultrastructural, and genomic detail for a novel hepatopancreatic, nucleus-specific virus - Menippe mercenaria nudivirus (MmNV) - isolated from M. mercenaria, captured during fisheries-independent monitoring. The virus has a genome of 99,336 bp and encodes 84 predicted protein coding genes and shows greatest similarity to Aratus pisonii nudivirus (ApNV) (<60% protein similarity and 31 shared genes of greatest similarity), collected from the Florida Keys, USA. MmNV is a member of the Gammanudivirus genus (Naldaviricetes: Lefavirales: Nudiviridae). Comparisons of virus genome size, preferred host environment, and gene number revealed no clear associations between the viral traits and phylogenetic position. Evolution of the virus alongside the diversification of host taxa, with the potential for host-switching, remain more likely evolutionary pathways.

Structural Evolution of Gene Promoters Driven by Primate-Specific KRAB Zinc Finger Proteins

Krüppel-associated box (KRAB) zinc finger proteins (KZNFs) recognize and repress transposable elements (TEs); TEs are DNA elements that are capable of replicating themselves throughout our genomes with potentially harmful consequences. However, genes from this family of transcription factors have a much wider potential for genomic regulation. KZNFs have become integrated into gene-regulatory networks through the control of TEs that function as enhancers and gene promoters; some KZNFs also bind directly to gene promoters, suggesting an additional, more direct layer of KZNF co-option into gene-regulatory networks. Binding site analysis of ZNF519, ZNF441, and ZNF468 suggests the structural evolution of KZNFs to recognize TEs can result in coincidental binding to gene promoters independent of TE sequences. We show a higher rate of sequence turnover in gene promoter KZNF binding sites than neighboring regions, implying a selective pressure is being applied by the binding of a KZNF. Through CRISPR/Cas9 mediated genetic deletion of ZNF519, ZNF441, and ZNF468, we provide further evidence for genome-wide co-option of the KZNF-mediated gene-regulatory functions; KZNF knockout leads to changes in expression of KZNF-bound genes in neuronal lineages. Finally, we show that the opposite can be established upon KZNF overexpression, further strengthening the support for the role of KZNFs as bona-fide gene regulators. With no eminent role for ZNF519 in controlling its TE target, our study may provide a snapshot into the early stages of the completed co-option of a KZNF, showing the lasting, multilayered impact that retrovirus invasions and host response mechanisms can have upon the evolution of our genomes.

In Silico Analysis: Genome-Wide Identification, Characterization and Evolutionary Adaptations of Bone Morphogenetic Protein (BMP) Gene Family in Homo sapiens

We systematically analyzed BMP gene family in H. sapiens to elucidate genetic structure, phylogenetic relationships, adaptive evolution and tissue-specific expression pattern. Total of 13 BMPs genes were identified in the H. sapiens genome. Bone morphogenetic proteins (BMPs) are composed of a variable number of exons ranging from 2 to 21. They exhibit a molecular weight ranging from 31,081.81 to 82,899.61 Da. These proteins possess hydrophilic characteristics, display thermostability, and exhibit a pH range from acidic to basic. We identified four segmental and two tandem duplication events in BMP gene family of H. sapiens. All of the vertebrate species that were studied show the presence of BMPs 1, 2, 3, 4, 5, 6, 7, 8A, and 15, however only Homo sapiens demonstrated the presence of BMP9 and BMP11. The pathway and process enrichment analysis of BMPs genes showed that these were considerably enriched in positive regulation of pathway-restricted SMAD protein phosphorylation (92%) and cartilage development (77%) biological processes. These genes exhibited positive selection signals that were shown to be conserved across vertebrate lineages. The results showed that BMP2/3/5/6/8a/15 proteins underwent adaptive selection at many amino acid locations and increased positive selection was detected in TGF-β propeptide and TGF-β super family domains which were involved in dorso-ventral patterning, limb bud development. More over the expression pattern of BMP genes revealed that BMP1 and BMP5; BMP4 and BMP6 exhibited substantially identical expression patterns in all tissues while BMP10, BMP15, and BMP3 showed tissue-specific expression.

Genome-wide identification and expression analysis of Na+/H+antiporter (NHX) genes in tomato under salt stress

Plant Na +/H + antiporter (NHX) genes enhance salt tolerance by preventing excessive Na+ accumulation in the cytosol through partitioning of Na+ ions into vacuoles or extracellular transport across the plasma membrane. However, there is limited detailed information regarding the salt stress responsive SlNHXs in the most recent tomato genome. We investigated the role of this gene family's expression patterns in the open flower tissues under salt shock in Solanum lycopersicum using a genome-wide approach. A total of seven putative SlNHX genes located on chromosomes 1, 4, 6, and 10 were identified, but no ortholog of the NHX5 gene was identified in the tomato genome. Phylogenetic analysis revealed that these genes are divided into three different groups. SlNHX proteins with 10-12 transmembrane domains were hypothetically localized in vacuoles or cell membranes. Promoter analysis revealed that SlNHX6 and SlNHX8 are involved with the stress-related MeJA hormone in response to salt stress signaling. The structural motif analysis of SlNHX1, -2, -3, -4, and -6 proteins showed that they have highly conserved amiloride binding sites. The protein-protein network revealed that SlNHX7 and SlNHX8 interact physically with Salt Overly Sensitive (SOS) pathway proteins. Transcriptome analysis demonstrated that the SlNHX2 and SlNHX6 genes were substantially expressed in the open flower tissues. Moreover, quantitative PCR analysis indicated that all SlNHX genes, particularly SlNHX6 and SlNHX8, are significantly upregulated by salt shock in the open flower tissues. Our results provide an updated framework for future genetic research and development of breeding strategies against salt stress in the tomato.

Dysregulation of MicroRNAs in Adult Osteogenesis Imperfecta: The miROI Study

As epigenetic regulators of gene expression, circulating micro-RiboNucleic Acids (miRNAs) have been described in several bone diseases as potential prognostic markers. The aim of our study was to identify circulating miRNAs potentially associated with the severity of osteogenesis imperfecta (OI) in three steps. We have screened by RNA sequencing for the miRNAs that were differentially expressed in sera of a small group of OI patients versus controls and then conducted a validation phase by RT-qPCR analysis of sera of a larger patient population. In the first phase of miROI, we found 79 miRNAs that were significantly differentially expressed. We therefore selected 19 of them as the most relevant. In the second phase, we were able to validate the significant overexpression of 8 miRNAs in the larger OI group. Finally, we looked for a relationship between the level of variation of the validated miRNAs and the clinical characteristics of OI. We found a significant difference in the expression of two microRNAs in those patients with dentinogenesis imperfecta. After reviewing the literature, we found 6 of the 8 miRNAs already known to have a direct action on bone homeostasis. Furthermore, the use of a miRNA-gene interaction prediction model revealed a 100% probability of interaction between 2 of the 8 confirmed miRNAs and COL1A1 and/or COL1A2. This is the first study to establish the miRNA signature in OI, showing a significant modification of miRNA expression potentially involved in the regulation of genes involved in the physiopathology of OI. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

IgE response to Aed al 13 and Aed al 14 recombinant allergens from Aedes albopictus saliva in humans

Background

Mosquito bite is normally associated with mild allergic responses, but severe localized or systemic reactions are also possible. Reliable tools for the diagnosis of mosquito allergy are still unavailable. Here, we investigated the IgE response to 3 potential salivary allergens identified in the saliva of the tiger mosquito Aedes albopictus.

Methods

Serum from 55 adult individuals (28 controls and 27 allergic people), were analysed using an in-house Enzyme Linked ImmunoSorbent Assay (ELISA) against the Salivary Gland Extract (SGE) and the recombinant proteins albD7l2 (Aed al 2), albAntigen5-3 (Aed al 13) and albLIPS-2 (Aed al 14).

Results

Fifteen of the 27 (56%) individuals having hypersensitive reactions to mosquito bites had IgE serum levels recognizing SGE. Negative sera did not show detectable levels of IgE targeting the SGE from the most common sympatric mosquito Culex pipiens. Among the positive individuals, 2 subjects displayed IgE targeting Aed al 2 (13%), while IgE recognizing Aed al 13 and Aed al 14 were detected in ten (67%) and seven (47%) individuals, respectively. Two sera from non-hypersensitive subjects had detectable levels of IgE targeting Aed al 13, suggesting possible cross-reaction with the homologue salivary proteins of multiple mosquito species or, more generally, of hematophagous insects.

Conclusions

Our results indicate that Aed al 13 and Aed al 14 hold the potential to be developed as tools for the diagnosis of allergy to Ae. albopictus bites. Such tools would facilitate epidemiological studies on tiger mosquito allergy in humans and might foster the development of further protein-based assays to investigate cross-species allergies.

Deficiency of ligase IV leads to reduced NHEJ, accumulation of DNA damage, and can sensitize cells to cancer therapeutics

Ligase IV is a key enzyme involved during DNA double-strand breaks (DSBs) repair through nonhomologous end joining (NHEJ). However, in contrast to Ligase IV deficient mouse cells, which are embryonic lethal, Ligase IV deficient human cells, including pre-B cells, are viable. Using CRISPR-Cas9 mediated genome editing, we have generated six different LIG4 mutants in cervical cancer and normal kidney epithelial cell lines. While the LIG4 mutant cells showed a significant reduction in NHEJ, joining mediated through microhomology-mediated end joining (MMEJ) and homologous recombination (HR) were significantly high. The reduced NHEJ joining activity was restored by adding purified Ligase IV/XRCC4. Accumulation of DSBs and reduced cell viability were observed in LIG4 mutant cells. LIG4 mutant cells exhibited enhanced sensitivity towards DSB-inducing agents such as ionizing radiation (IR) and etoposide. More importantly, the LIG4 mutant of cervical cancer cells showed increased sensitivity towards FDA approved drugs such as Carboplatin, Cisplatin, Paclitaxel, Doxorubicin, and Bleomycin used for cervical cancer treatment. These drugs, in combination with IR showed enhanced cancer cell death in the background of LIG4 gene mutation. Thus, our study reveals that mutation in LIG4 results in compromised NHEJ, leading to sensitization of cervical cancer cells towards currently used cancer therapeutics.

RAS gene mutations and histomorphometric measurements in oral squamous cell carcinoma

Members of the RAS gene family frequently are mutated in cancers including oral squamous cell carcinoma (OSCC). We investigated the correlation of histological characteristics of OSCC with RAS gene mutations. We graded tumors and extracted genomic DNA from OSCC. The first two exons of KRAS, HRAS and NRAS genes were subjected to PCR amplification and DNA sequencing followed by bioinformatic analysis to explore the structural and functional impact of the mutations on encoding of proteins. Cellular and nuclear diameters in histological sections were varied in all grades of cancer. Using sequence analysis, we identified nonsynonymous mutations in both HRAS (G12S, G15C, D54H, Q61H, Q61L, E62D, E63D, Q70E, Q70V) and NRAS (Q22P, K88R). Stop codon mutations, however, were observed in KRAS. Spatial orientation of substituted amino acids was observed despite conservation of overall structure of variant proteins. Our findings suggest that KRAS may be mutated more frequently in OSCC compared to HRAS and NRAS. Also, the histological features of nuclear and cellular diameter differed significantly between the KRAS mutated and unmutated cases.

Reactive architecture profiling with a methyl acyl phosphate electrophile

Activity-based protein profiling has facilitated the study of the activity of enzymes in proteomes, inhibitor development, and identification of enzymes that share mechanistic and active-site architectural features. Since methyl acyl phosphate monoesters act as electrostatically selective anionic electrophiles for the covalent modification of nucleophiles that reside adjacent to cationic sites in proteins, we synthesized methyl hex-5-ynoyl phosphate (MHP) to broadly target such protein architectures. After treating the soluble proteome of Paucimonas lemoignei with MHP, biotinylating the resulting acylated proteins using click chemistry, enriching the protein adducts using streptavidin, and analyzing the proteins by LC-MS/MS, a set of 240 enzymes and 132 non-enzyme proteins were identified for a wide spectrum of biological processes and from all 7 enzyme classes. Among those enzymes identified, β-hydroxybutyrate dehydrogenase (PlHBDH) and CTP synthase (E. coli orthologue, EcCTPS) were purified as recombinant enzymes and their rates of inactivation and sites of modification by MHP and methyl acetyl phosphate (MAP) were characterized. MHP reacted more slowly with these proteins than MAP but exhibited greater specificity, despite its lack of multiple binding determinants. Generally, MAP modified more surface residues than MHP. MHP specifically modified Ser 146, Lys 156, and Lys 163 at the active site of PlHBDH. MHP and MAP modified numerous residues of EcCTPS with CTP furnishing the greatest level of protection against MHP- and MAP-dependent modification and inactivation, respectively, followed by ATP and glutamine. Overall, MHP served as an effective probe to identify proteins that are potentially amenable to inhibition by methyl acyl phosphates.

HSP60 mimetic peptides from Mycobacterium leprae as new antigens for immunodiagnosis of Leprosy

The early diagnosis of leprosy serves as an important tool to reduce the incidence of this disease in the world. Phage display (PD) technology can be used for mapping new antigens to the development of immunodiagnostic platforms. Our objective was to identify peptides that mimic Mycobacterium leprae proteins as serological markers using phage display technology. The phages were obtained in the biopanning using negative and positive serum from household contacts and leprosy patients, respectively. Then, the peptides were synthesized and validated in silico and in vitro for detection of IgG from patients and contacts. To characterize the native protein of M. leprae, scFv antibodies were selected against the synthetic peptides by PD. The scFv binding protein was obtained by immunocapture and confirmed using mass spectrometry. We selected two phase-fused peptides, MPML12 and MPML14, which mimic the HSP60 protein from M. leprae. The peptides MPML12 and MPML14 obtained 100% and 92.85% positivity in lepromatous patients. MPML12 and MPM14 detect IgG, especially in the multibacillary forms. The MPML12 and MPML14 peptides had positivity of 11.1% and 16.6% in household contacts, respectively. There was no cross-reaction in patient's samples with visceral leishmaniasis, tuberculosis and other mycobacteriosis for both peptides. Given these results and the easy obtainment of mimetic antigens, our peptides are promising markers for application in the diagnosis of leprosy, especially in endemic and hyperendemic regions.

Comprehensive Genome-Wide Identification of the RNA-Binding Glycine-Rich Gene Family and Expression Profiling under Abiotic Stress in Brassica oleracea

The RNA-binding glycine-rich proteins (RBGs) of the glycine-rich protein family play vital roles in regulating gene expression both at the transcriptional and post-transcriptional levels. However, the members and functions in response to abiotic stresses of the RBG gene family remain unclear in Brassica oleracea. In this study, a total of 19 BoiRBG genes were identified through genome-wide analysis in broccoli. The characteristics of BoiRBG sequences and their evolution were examined. An analysis of synteny indicated that the expansion of the BoiRBG gene family was primarily driven by whole-genome duplication and tandem duplication events. The BoiRBG expression patterns revealed that these genes are involved in reaction to diverse abiotic stress conditions (i.e., simulated drought, salinity, heat, cold, and abscisic acid) and different organs. In the present research, the up-regulation of BoiRBGA13 expression was observed when subjected to both NaCl-induced and cold stress conditions in broccoli. Moreover, the overexpression of BoiRBGA13 resulted in a noteworthy reduction in taproot lengths under NaCl stress, as well as the inhibition of seed germination under cold stress in broccoli, indicating that RBGs play different roles under various stresses. This study provides insights into the evolution and functions of BoiRBG genes in Brassica oleracea and other Brassicaceae family plants.

In Silico Evaluation, Phylogenetic Analysis, and Structural Modeling of the Class II Hydrophobin Family from Different Fungal Phytopathogens

The class II hydrophobin group (HFBII) is an extracellular group of proteins that contain the HFBII domain and eight conserved cysteine residues. These proteins are exclusively secreted by fungi and have multiple functions with a probable role as effectors. In the present study, a total of 45 amino acid sequences of hydrophobin class II proteins from different phytopathogenic fungi were retrieved from the NCBI database. We used the integration of well-designed bioinformatic tools to characterize and predict their physicochemical parameters, novel motifs, 3D structures, multiple sequence alignment (MSA), evolution, and functions as effector proteins through molecular docking. The results revealed new features for these protein members. The ProtParam tool detected the hydrophobicity properties of all proteins except for one hydrophilic protein (KAI3335996.1). Out of 45 proteins, six of them were detected as GPI-anchored proteins by the PredGPI server. Different 3D structure templates with high pTM scores were designed by Multifold v1, AlphaFold2, and trRosetta. Most of the studied proteins were anticipated as apoplastic effectors and matched with the ghyd5 gene of Fusarium graminearum as virulence factors. A protein-protein interaction (PPI) analysis unraveled the molecular function of this group as GTP-binding proteins, while a molecular docking analysis detected a chitin-binding effector role. From the MSA analysis, it was observed that the HFBII sequences shared conserved 2 Pro (P) and 2 Gly (G) amino acids besides the known eight conserved cysteine residues. The evolutionary analysis and phylogenetic tree provided evidence of episodic diversifying selection at the branch level using the aBSREL tool. A detailed in silico analysis of this family and the present findings will provide a better understanding of the HFBII characters and evolutionary relationships, which could be very useful in future studies.

Apple CRISPR-Cas9-A Recipe for Successful Targeting of AGAMOUS-like Genes in Domestic Apple

Fruit trees and other fruiting hardwood perennials are economically valuable, and there is interest in developing improved varieties. Both conventional breeding and biotechnology approaches are being utilized towards the goal of developing advanced cultivars. Increased knowledge of the effectiveness and efficiency of biotechnology approaches can help guide use of the CRISPR gene-editing technology. Here, we examined CRISPR-Cas9-directed genome editing in the valuable commodity fruit tree Malus x domestica (domestic apple). We transformed two cultivars with dual CRISPR-Cas9 constructs designed to target two AGAMOUS-like genes simultaneously. The main goal was to determine the effectiveness of this approach for achieving target gene changes. We obtained 6 Cas9 control and 38 independent CRISPR-Cas9 events. Of the 38 CRISPR-Cas9 events, 34 (89%) had gene edits and 14 (37%) showed changes to all alleles of both target genes. The most common change was large deletions, which were present in 59% of all changed alleles, followed by small deletions (21%), small insertions (12%), and a combination of small insertions and deletions (8%). Overall, a high rate of successful gene alterations was found. Many of these changes are predicted to cause frameshifts and alterations to the predicted peptides. Future work will include monitoring the floral development and floral form.

Genome-Wide Identification and Analysis of MYB Transcription Factors in Pyropia yezoensis

MYB transcription factors are one of the largest transcription factor families in plants, and they regulate numerous biological processes. Red algae are an important taxonomic group and have important roles in economics and research. However, no comprehensive analysis of the MYB gene family in any red algae, including Pyropia yezoensis, has been conducted. To identify the MYB gene members of Py. yezoensis, and to investigate their family structural features and expression profile characteristics, a study was conducted. In this study, 3 R2R3-MYBs and 13 MYB-related members were identified in Py. yezoensis. Phylogenetic analysis indicated that most red algae MYB genes could be clustered with green plants or Glaucophyta MYB genes, inferring their ancient origins. Synteny analysis indicated that 13 and 5 PyMYB genes were orthologous to Pyropia haitanensis and Porphyra umbilicalis, respectively. Most Bangiaceae MYB genes contain several Gly-rich motifs, which may be the result of an adaptation to carbon limitations and maintenance of important regulatory functions. An expression profile analysis showed that PyMYB genes exhibited diverse expression profiles. However, the expression patterns of different members appeared to be diverse, and PyMYB5 was upregulated in response to dehydration, low temperature, and Pythium porphyrae infection. This is the first comprehensive study of the MYB gene family in Py. Yezoensis and it provides vital insights into the functional divergence of MYB genes.

Triosephosphate Isomerase: The Crippling Effect of the P168A/I172A Substitution at the Heart of an Enzyme Active Site

The P168 and I172 side chains sit at the heart of the active site of triosephosphate isomerase (TIM) and play important roles in the catalysis of the isomerization reaction. The phosphodianion of substrate glyceraldehyde 3-phosphate (GAP) drives a conformational change at the TIM that creates a steric interaction with the P168 side chain that is relieved by the movement of P168 that carries the basic E167 side chain into a clamp that consists of the hydrophobic I172 and L232 side chains. The P168A/I172A substitution at TIM from Trypanosoma brucei brucei (TbbTIM) causes a large 120,000-fold decrease in kcat for isomerization of GAP that eliminates most of the difference in the reactivity of TIM compared to the small amine base quinuclidinone for deprotonation of catalyst-bound GAP. The I172A substitution causes a > 2-unit decrease in the pKa of the E167 carboxylic acid in a complex to the intermediate analog PGA, but the P168A substitution at the I172A variant has no further effect on this pKa. The P168A/I172A substitutions cause a 5-fold decrease in Km for the isomerization of GAP from a 0.9 kcal/mol stabilization of the substrate Michaelis complexes. The results show that the P168 and I172 side chains play a dual role in destabilizing the ground-state Michaelis complex to GAP and in promoting stabilization of the transition state for substrate isomerization. This is consistent with an important role for these side chains in an induced fit reaction mechanism [Richard, J. P. (2022) Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution. Biochemistry 61, 1533-1542].

Identification and characterization of the cupin_1 domain-containing proteins in ma bamboo (Dendrocalamus latiflorus) and their potential role in rhizome sprouting

Cupin_1 domain-containing protein (CDP) family, which is a member of the cupin superfamily with the most diverse functions in plants, has been found to be involved in hormone pathways that are closely related to rhizome sprouting (RS), a vital form of asexual reproduction in plants. Ma bamboo is a typical clumping bamboo, which mainly reproduces by RS. In this study, we identified and characterized 53 Dendrocalamus latiflorus CDP genes and divided them into seven subfamilies. Comparing the genetic structures among subfamilies showed a relatively conserved gene structure within each subfamily, and the number of cupin_1 domains affected the conservation among D. latiflorus CDP genes. Gene collinearity results showed that segmental duplication and tandem duplication both contributed to the expansion of D. latiflorus CDP genes, and lineage-specific gene duplication was an important factor influencing the evolution of CDP genes. Expression patterns showed that CDP genes generally had higher expression levels in germinating underground buds, indicating that they might play important roles in promoting shoot sprouting. Transcription factor binding site prediction and co-expression network analysis indicated that D. latiflorus CDPs were regulated by a large number of transcription factors, and collectively participated in rhizome buds and shoot development. This study significantly provided new insights into the evolutionary patterns and molecular functions of CDP genes, and laid a foundation for further studying the regulatory mechanisms of plant rhizome sprouting.

The Cytotoxic Mycobacteriophage Protein Phaedrus gp82 Interacts with and Modulates the Activity of the Host ATPase, MoxR

Approximately 70% of bacteriophage-encoded proteins are of unknown function. Elucidating these protein functions represents opportunities to discover new phage-host interactions and mechanisms by which the phages modulate host activities. Here, we describe a pipeline for prioritizing phage-encoded proteins for structural analysis and characterize the gp82 protein encoded by mycobacteriophage Phaedrus. Structural and solution studies of gp82 show it is a trimeric protein containing two domains. Co-precipitation studies with the host Mycobacterium smegmatis identified the ATPase MoxR as an interacting partner protein. Phaedrus gp82-MoxR interaction requires the presence of a loop sequence within gp82 that is highly exposed and disordered in the crystallographic structure. We show that Phaedrus gp82 overexpression in M. smegmatis retards the growth of M. smegmatis on solid medium, resulting in a small colony phenotype. Overexpression of gp82 containing a mutant disordered loop or the overexpression of MoxR both rescue this phenotype. Lastly, we show that recombinant gp82 reduces levels of MoxR-mediated ATPase activity in vitro that is required for its chaperone function, and that the disordered loop plays an important role in this phenotype. We conclude that Phaedrus gp82 binds to and reduces mycobacterial MoxR activity, leading to reduced function of host proteins that require MoxR chaperone activity for their normal activity.

Immunoinformatics aided designing of a next generation poly-epitope vaccine against uropathogenic Escherichia coli to combat urinary tract infections

Urinary tract infections (UTIs) are the second most prevalent bacterial infections and uropathogenic Escherichia coli (UPEC) stands among the primary causative agents of UTIs. The usage of antibiotics is the routine therapy being used in various countries to treat UTIs but becoming ineffective because of increasing antibiotic resistance among UPEC strains. Thus, there must be the development of some alternative treatment strategies such as vaccine development against UPEC. In the following study, pan-genomics along with reverse vaccinology approaches is used under the framework of bioinformatics for the identification of core putative vaccine candidates, employing 307 UPEC genomes (complete and draft), available publicly. A total of nine T-cell epitopes (derived from B-cells) of both MHC classes (I and II), were prioritized among three potential protein candidates. These epitopes were then docked together by using linkers (GPGPG and AAY) and an adjuvant (Cholera Toxin B) to form a poly-valent vaccine construct. The chimeric vaccine construct was undergone by molecular modelling, further refinement and energy minimization. We predicted positive results of the vaccine construct in immune simulations with significantly high levels of immune cells. The protein-protein docking analysis of vaccine construct with toll-like receptors predicted efficient binding, which was further validated by molecular dynamics simulation of vaccine construct with TLR-2 and TLR-4 at 120 ns, resulting in stable complexes' conformation throughout the simulation run. Overall, the vaccine construct demonstrated positive antigenic response. In future, this chimeric vaccine construct or the identified epitopes could be experimentally validated for the development of UPEC vaccines against UTIs.Communicated by Ramaswamy H. Sarma.

Zucchini Yellow Mosaic Virus (ZYMV) as a Serious Biotic Stress to Cucurbits: Prevalence, Diversity, and Its Implications for Crop Sustainability

Zucchini yellow mosaic virus (ZYMV) is a severe threat to cucurbit crops worldwide, including Pakistan. This study was pursued to evaluate the prevalence, geographic distribution, and molecular diversity of ZYMV isolates infecting cucurbits in Pakistan's Pothwar region. Almost all the plant viruses act as a biotic stress on the host plants, which results in a yield loss. These viruses cause losses in single-infection or in mixed-infection cucurbit crops, and we have found a number of mixed-infected samples belonging to the Curubitaceae family. Serological detection of the tested potyviruses in the collected cucurbit samples revealed that ZYMV was the most prevalent virus, with a disease incidence (DI) at 35.2%, followed by Papaya ringspot virus (PRSV) with an incidence of 2.2%, and Watermelon mosaic virus (WMV) having an incidence as little as 0.5% in 2016. In the year 2017, a relatively higher disease incidence of 39.7%, 2.4%, and 0.3% for ZYMV, WMV, and PRSV, respectively, was recorded. ZYMV was the most prevalent virus with the highest incidence in Attock, Rawalpindi, and Islamabad, while PRSV was observed to be the highest in Islamabad and Jhelum. WMV infection was observed only in Rawalpindi and Chakwal. Newly detected Pakistani ZYMV isolates shared 95.8-97.0% nucleotide identities among themselves and 77.1-97.8% with other isolates retrieved from GenBank. Phylogenetic relationships obtained using different ZYMV isolates retrieved from GenBank and validated by in silico restriction analysis revealed that four Pakistani isolates clustered with other ZYMV isolates in group IIb with Chinese, Italian, Polish, and French isolates, while another isolate (MK848239) formed a separate minor clade within IIb. The isolate MK8482490, reported to infect bitter gourd in Pakistan, shared a minor clade with a Chinese isolate (KX884570). Recombination analysis revealed that the recently found ZYMV isolate (MK848239) is most likely a recombinant of Pakistani (MK848237) and Italian (MK956829) isolates, with a recombinant breakpoint between 266 and 814 nucleotide positions. Local isolate comparison and recombination detection may aid in the development of a breeding program that identifies resistant sources against recombinant isolates because the ZYMV is prevalent in a few cucurbit species grown in the surveyed areas and causes heavy losses and economic damage to the agricultural community.

Genome-wide identification and comprehensive analysis of WRKY transcription factor family in safflower during drought stress

The WRKY family is an important family of transcription factors in plant development and stress response. Currently, there are few reports on the WRKY gene family in safflower (Carthamus tinctorius L.). In this study, a total of 82 CtWRKY genes were identified from the safflower genome and could be classified into 3 major groups and 5 subgroups based on their structural and phylogenetic characteristics. The results of gene structure, conserved domain and motif analyses indicated that CtWRKYs within the same subfamily maintained a consistent exon/intron organization and composition. Chromosomal localization and gene duplication analysis results showed that CtWRKYs were randomly localized on 12 chromosomes and that fragment duplication and purification selection may have played an important role in the evolution of the WRKY gene family in safflower. Promoter cis-acting element analysis revealed that the CtWRKYs contain many abiotic stress response elements and hormone response elements. Transcriptome data and qRT-PCR analyses revealed that the expression of CtWRKYs showed tissue specificity and a strong response to drought stress. Notably, the expression level of the CtWRKY55 gene rapidly increased more than eightfold under drought treatment and rehydration, indicating that it may be a key gene in response to drought stress. These results provide useful insights for investigating the regulatory function of the CtWRKY gene in safflower growth and development, as well as identifying key genes for future molecular breeding programmes.

Large-Scale Identification and Characterization Analysis of VQ Family Genes in Plants, Especially Gymnosperms

VQ motif-containing (VQ) proteins are a class of transcription regulatory cofactors widely present in plants, playing crucial roles in growth and development, stress response, and defense. Although there have been some reports on the member identification and functional research of VQ genes in some plants, there is still a lack of large-scale identification and clear graphical presentation of their basic characterization information to help us to better understand this family. Especially in gymnosperms, the VQ family genes and their evolutionary relationships have not yet been reported. In this study, we systematically identified 2469 VQ genes from 56 plant species, including bryophytes, gymnosperms, and angiosperms, and analyzed their molecular and evolutionary features. We found that amino acids are only highly conserved in the VQ domain, while other positions are relatively variable; most VQ genes encode relatively small proteins and do not have introns. The GC content in Poaceae plants is the highest (up to 70%); these VQ proteins can be divided into nine subgroups. In particular, we analyzed the molecular characteristics, chromosome distribution, duplication events, and expression levels of VQ genes in three gymnosperms: Ginkgo biloba, Taxus chinensis, and Pinus tabuliformis. In gymnosperms, VQ genes are classified into 11 groups, with highly similar motifs in each group; most VQ proteins have less than 300 amino acids and are predicted to be located in nucleus. Tandem duplication is an important driving force for the expansion of the VQ gene family, and the evolutionary processes of most VQ genes and duplication events are relatively independent; some candidate VQ genes are preliminarily screened, and they are likely to be involved in plant growth and stress and defense responses. These results provide detailed information and powerful references for further understanding and utilizing the VQ family genes in various plants.

Himalayan upliftment and Shiwalik succession act as a cradle for divergence in Bengal monitor lizard Varanus bengalensis (Reptilia: Varanidae) in India

The Himalayan foothills and associated environment are well-known for driving the rapid diversification of many species and the formation of biodiversity hotspots. The effects of environmental change since the Miocene have accelerated species diversification, and hence are useful for studying population genetic structure, and evolutionary relationships via genetic approaches. To date, the effects of climatic fluctuations on the biogeography of large-bodied lizards have not been assessed comprehensively. Herein, we examine the diversification of Varanus bengalensis, focusing on its genetic structure to provide insights into how landscape structure and climatic fluctuations have shaped species differentiation. We confirm the existence of two distinct lineages within V. bengalensis distributed across the Himalayan foothills and the remainder of mainland India. Divergence analyses revealed the split between the Himalayan foothills and the remainder of the mainland lineages of V. bengalensis in the mid-Pliocene ~3.06 Ma, potentially as a consequence of the Siwalik broadening and climatic fluctuations across the Himalayan foothills. The results suggest recognition of a new lineage of V. bengalensis from the Himalayan foothills as a distinctive evolutionarily significant unit.

Truncated glycosyltransferase coding regions in novel ABO alleles give rise to weak A or B blood group expression and discrepant typing results

BACKGROUND

Correct ABO blood-group matching between donor and patient is crucial for safe transfusions. We investigated the underlying reason causing inconclusive ABO serology in samples referred to our laboratory.

STUDY DESIGN AND METHODS

Flow cytometric analysis, ABO genotyping, and sequencing were used to characterize ABO-discrepant blood samples (n = 13). ABO gene variants were inserted in a GFP-containing bicistronic vector to assess A/B expression following overexpression in HeLa cells.

RESULTS

Seven novel alleles with nonsense mutations predicted to truncate the encoded ABO glycosyltransferases were identified. While these variants could represent O alleles, serology showed signs of ABO glycosyltransferase activity. ABO*A1.01-related alleles displayed remarkably characteristic percentages of A-positive cells for samples with the same variant: c.42C>A (p.Cys14*; 10%), c.102C>A (p.Tyr34*; 31%-32%, n = 2), c.106dup (p.Val36Glyfs*21; 16%-17%, n = 3) or c.181_182ins (p.Leu61Argfs*21; 12%-13%, n = 2). Transfection studies confirmed significantly decreased A expression compared to wild type. The remaining variants were found on ABO*B.01 background: c.1_5dup (pGly3Trpfs*20), c.15dup (p.Arg6Alafs*51) or c.496del (p.Thr166Profs*26). Although the absence of plasma anti-B was noted overall, B antigen expression was barely detected on erythrocytes. Overexpression confirmed decreased B in two variants compared to wildtype while c.1_5dup only showed a non-significant downward trend.

CONCLUSION

Samples displaying aberrant ABO serology revealed seven principally interesting alleles. Despite the presence of truncating mutations, normally resulting in null alleles, low levels of ABO antigens were detectable where alterations affected ABO exons 1-4 but not exon 7. This is compatible with the previously proposed concept that alternative start codons in early exons can be used to initiate the translation of functional ABO glycosyltransferase.

Biochemical and X-ray analyses of the players involved in the faRel2/aTfaRel2 toxin-antitoxin operon

The aTfaRel2/faRel2 operon from Coprobacillus sp. D7 encodes a bicistronic type II toxin-antitoxin (TA) module. The FaRel2 toxin is a toxic small alarmone synthetase (toxSAS) that inhibits translation through the pyrophosphorylation of uncharged tRNAs at the 3'-CCA end. The toxin is neutralized by the antitoxin ATfaRel2 through the formation of an inactive TA complex. Here, the production, biophysical analysis and crystallization of ATfaRel2 and FaRel2 as well as of the ATfaRel2-FaRel2 complex are reported. ATfaRel2 is monomeric in solution. The antitoxin crystallized in space group P21212 with unit-cell parameters a = 53.3, b = 34.2, c = 37.6 Å, and the best crystal diffracted to a resolution of 1.24 Å. Crystals of FaRel2 in complex with APCPP, a nonhydrolysable ATP analogue, belonged to space group P21, with unit-cell parameters a = 31.5, b = 60.6, c = 177.2 Å, β = 90.6°, and diffracted to 2.6 Å resolution. The ATfaRel2-FaRel2Y128F complex forms a heterotetramer in solution composed of two toxins and two antitoxins. This complex crystallized in two space groups: F4132, with unit-cell parameters a = b = c = 227.1 Å, and P212121, with unit-cell parameters a = 51.7, b = 106.2, c = 135.1 Å. The crystals diffracted to 1.98 and 2.1 Å resolution, respectively.

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

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

Identification of the zebrafish homologues of IMPG2, a retinal proteoglycan

Photoreceptor outer segments are surrounded by a carbohydrate-rich matrix, the interphotoreceptor matrix, necessary for physiological retinal function. Few roles for molecules characterizing the interphotoreceptor matrix have been clearly defined. Recent studies have found the presence of nonsense mutations in the interphotoreceptor matrix proteoglycan 2 (IMPG2) gene in patients affected by retinal dystrophies. IMPG2 encodes for a proteoglycan synthesized by photoreceptors and secreted in the interphotoreceptor matrix. Little is known about the structure and function of this protein, we thus decided to characterize zebrafish impg2. In zebrafish there are two Impg2 proteins, Impg2a and Impg2b. We generated a phylogenetic tree based on IMPG2 protein sequence similarity among vertebrates, showing a significant similarity between humans and teleosts. The human and zebrafish proteins share conserved domains, as also shown by homology models. Expression analyses of impg2a and impg2b show a continued expression in the photoreceptor layer starting from developmental stages and continuing through adulthood. Between 1 and 6 months post-fertilization, there is a significant shift of Impg2 expression toward the outer segment region, suggesting an increase in secretion. This raises intriguing hypotheses about its possible role(s) during retinal maturation, laying the groundwork for the generation of most needed models for the study of IMPG2-related inherited retinal dystrophies.

Potential for host-symbiont communication via neurotransmitters and neuromodulators in an aneural animal, the marine sponge Amphimedon queenslandica

Interkingdom signalling within a holobiont allows host and symbionts to communicate and to regulate each other's physiological and developmental states. Here we show that a suite of signalling molecules that function as neurotransmitters and neuromodulators in most animals with nervous systems, specifically dopamine and trace amines, are produced exclusively by the bacterial symbionts of the demosponge Amphimedon queenslandica. Although sponges do not possess a nervous system, A. queenslandica expresses rhodopsin class G-protein-coupled receptors that are structurally similar to dopamine and trace amine receptors. When sponge larvae, which express these receptors, are exposed to agonists and antagonists of bilaterian dopamine and trace amine receptors, we observe marked changes in larval phototactic swimming behaviour, consistent with the sponge being competent to recognise and respond to symbiont-derived trace amine signals. These results indicate that monoamines synthesised by bacterial symbionts may be able to influence the physiology of the host sponge.

Urine-HILIC: Automated Sample Preparation for Bottom-Up Urinary Proteome Profiling in Clinical Proteomics

Urine provides a diverse source of information related to a patient's health status and is ideal for clinical proteomics due to its ease of collection. To date, most methods for the preparation of urine samples lack the throughput required to analyze large clinical cohorts. To this end, we developed a novel workflow, urine-HILIC (uHLC), based on an on-bead protein capture, clean-up, and digestion without the need for bottleneck processing steps such as protein precipitation or centrifugation. The workflow was applied to an acute kidney injury (AKI) pilot study. Urine from clinical samples and a pooled sample was subjected to automated sample preparation in a KingFisher™ Flex magnetic handling station using the novel approach based on MagReSyn® HILIC microspheres. For benchmarking, the pooled sample was also prepared using a published protocol based on an on-membrane (OM) protein capture and digestion workflow. Peptides were analyzed by LCMS in data-independent acquisition (DIA) mode using a Dionex Ultimate 3000 UPLC coupled to a Sciex 5600 mass spectrometer. The data were searched in Spectronaut™ 17. Both workflows showed similar peptide and protein identifications in the pooled sample. The uHLC workflow was easier to set up and complete, having less hands-on time than the OM method, with fewer manual processing steps. Lower peptide and protein coefficient of variation was observed in the uHLC technical replicates. Following statistical analysis, candidate protein markers were filtered, at ≥8.35-fold change in abundance, ≥2 unique peptides and ≤1% false discovery rate, and revealed 121 significant, differentially abundant proteins, some of which have known associations with kidney injury. The pilot data derived using this novel workflow provide information on the urinary proteome of patients with AKI. Further exploration in a larger cohort using this novel high-throughput method is warranted.

Genome-Wide Identification and Expression Analysis of the TCP Gene Family Related to Developmental and Abiotic Stress in Ginger

Ginger (Zingiber officinale Roscoe), a widely consumed edible and medicinal plant, possesses significant nutritional and economic value. Abiotic stresses such as drought and low temperatures can impact the growth and development of ginger. The plant-specific transcription factor Teosinte branched1/cycloidea/proliferating cell factor (TCP) has progressively been identified in various plants for its role in regulating plant growth and development as well as conferring resistance to abiotic stresses. However, limited information on the TCP family is available in ginger. In this study, we identified 20 TCP members in the ginger genome, which were randomly distributed across 9 chromosomes. Based on phylogenetic analysis, these ginger TCP were classified into two subfamilies: Class I (PCF) and Class II (CIN, CYC/TB). The classification of the identified ginger TCPs was supported by a multi-species phylogenetic tree and motif structure analysis, suggesting that the amplification of the ginger TCP gene family occurred prior to the differentiation of angiosperms. The promoter region of ginger TCP genes was found to contain numerous cis-acting elements associated with plant growth, development, and abiotic stress response. Among these elements, the stress response element, anaerobic induction, and MYB binding site play a dominant role in drought responsiveness. Additionally, expression pattern analysis revealed variations in the expression of ginger TCP gene among different tissues and in response to diverse abiotic stresses (drought, low temperature, heat, and salt). Our research offers a thorough examination of TCP members within the ginger plant. This analysis greatly contributes to the understanding of how TCP genes regulate tissue development and response to stress, opening up new avenues for further exploration in this field.

Ets-1 transcription factor regulates glial cell regeneration and function in planarians

Glia play multifaceted roles in nervous systems in response to injury. Depending on the species, extent of injury and glial cell type in question, glia can help or hinder the regeneration of neurons. Studying glia in the context of successful regeneration could reveal features of pro-regenerative glia that could be exploited for new human therapies. Planarian flatworms completely regenerate their nervous systems after injury - including glia - and thus provide a strong model system for exploring glia in the context of regeneration. Here, we report that planarian glia regenerate after neurons, and that neurons are required for correct glial numbers and localization during regeneration. We also identify the planarian transcription factor-encoding gene ets-1 as a key regulator of glial cell maintenance and regeneration. Using ets-1 (RNAi) to perturb glia, we show that glial loss is associated with altered neuronal gene expression, impeded animal movement and impaired nervous system architecture - particularly within the neuropil. Importantly, our work reveals the inter-relationships of glia and neurons in the context of robust neural regeneration.

Identifying novel mechanisms of biallelic TP53 loss refines poor outcome for patients with multiple myeloma

Biallelic TP53 inactivation is the most important high-risk factor associated with poor survival in multiple myeloma. Classical biallelic TP53 inactivation has been defined as simultaneous mutation and copy number loss in most studies; however, numerous studies have demonstrated that other factors could lead to the inactivation of TP53. Here, we hypothesized that novel biallelic TP53 inactivated samples existed in the multiple myeloma population. A random forest regression model that exploited an expression signature of 16 differentially expressed genes between classical biallelic TP53 and TP53 wild-type samples was subsequently established and used to identify novel biallelic TP53 samples from monoallelic TP53 groups. The model reflected high accuracy and robust performance in newly diagnosed relapsed and refractory populations. Patient survival of classical and novel biallelic TP53 samples was consistently much worse than those with mono-allelic or wild-type TP53 status. We also demonstrated that some predicted biallelic TP53 samples simultaneously had copy number loss and aberrant splicing, resulting in overexpression of high-risk transcript variants, leading to biallelic inactivation. We discovered that splice site mutation and overexpression of the splicing factor MED18 were reasons for aberrant splicing. Taken together, our study unveiled the complex transcriptome of TP53, some of which might benefit future studies targeting abnormal TP53.

Ginsentide-like Coffeetides Isolated from Coffee Waste Are Cell-Penetrating and Metal-Binding Microproteins

Coffee processing generates a huge amount of waste that contains many natural products. Here, we report the discovery of a panel of novel cell-penetrating and metal ion-binding microproteins designated coffeetide cC1a-c and cL1-6 from the husk of two popular coffee plants, Coffea canephora and Coffea liberica, respectively. Combining sequence determination and a database search, we show that the prototypic coffeetide cC1a is a 37-residue, eight-cysteine microprotein with a hevein-like cysteine motif, but without a chitin-binding domain. NMR determination of cC1a reveals a compact structure that confers its resistance to heat and proteolytic degradation. Disulfide mapping together with chemical synthesis reveals that cC1a has a ginsentide-like, and not a hevein-like, disulfide connectivity. In addition, transcriptomic analysis showed that the 98-residue micrcoproten-like coffeetide precursor contains a three-domain arrangement, like ginsentide precursors. Molecular modeling, together with experimental validation, revealed a Mg2+ and Fe3+ binding pocket at the N-terminus formed by three glutamic acids. Importantly, cC1a is amphipathic with a continuous stretch of 19 apolar amino acids, which enables its cell penetration to target intracellular proteins, despite being highly negatively charged. Our findings suggest that coffee by-products could provide a source of ginsentide-like bioactive peptides that have the potential to target intracellular proteins.

Hsp70 Gene Family in Sebastiscus marmoratus: The Genome-Wide Identification and Transcriptome Analysis under Thermal Stress

Heat shock protein 70 kDa (Hsp70) is a highly conserved heat stress protein that is important in biotic processes and responses to abiotic stress. Hsp70 genes may be important in Sebastiscus marmoratus, for it is a kind of nearshore reef fish, and habitat temperature change is more drastic during development. However, genome-wide identification and expression analysis in the Hsp70 gene family of S. marmoratus are still lacking. Here, a total of 15 Hsp70 genes in the genome of S. marmoratus are identified, and their expression patterns were investigated using transcriptomic data from thermal stress experiments. The expansion and gene duplication events of Hsp70 genes from the Hspa4, Hspa8, and Hspa12a subfamilies in S. marmoratus are revealed by phylogenetic analysis. qRT-PCR expression patterns demonstrated that seven Hsp70 genes were significantly up-regulated and none were significantly down-regulated after heat treatment. Only the hsp70 gene was significantly up-regulated after cold treatment. The selection test further showed a purifying selection on the duplicated gene pairs, suggesting that these genes underwent subfunctionalization. Our results add novel insight to aquaculture and biological research on S. marmoratus, providing important information on how Hsp70 genes are regulated in Scorpaeniformes under thermal stress.

Exploring Proteus mirabilis Methionine tRNA Synthetase Active Site: Homology Model Construction, Molecular Dynamics, Pharmacophore and Docking Validation

Currently, the treatment of Proteus mirabilis infections is considered to be complicated as the organism has become resistant to numerous antibiotic classes. Therefore, new inhibitors should be developed, targeting bacterial molecular functions. Methionine tRNA synthetase (MetRS), a member of the aminoacyl-tRNA synthetase family, is essential for protein biosynthesis offering a promising target for novel antibiotics discovery. In the context of computer-aided drug design (CADD), the current research presents the construction and analysis of a comparative homology model for P. mirabilis MetRS, enabling development of novel inhibitors with greater selectivity. Molecular Operating Environment (MOE) software was used to build a homology model for P. mirabilis MetRS using Escherichia coli MetRS as a template. The model was evaluated, and the active site of the target protein predicted from its sequence using conservation analysis. Molecular dynamic simulations were performed to evaluate the stability of the modeled protein structure. In order to evaluate the predicted active site interactions, methionine (the natural substrate of MetRS) and several inhibitors of bacterial MetRS were docked into the constructed model using MOE. After validation of the model, pharmacophore-based virtual screening for a systemically prepared dataset of compounds was performed to prove the feasibility of the proposed model, identifying possible parent compounds for further development of MetRS inhibitors against P. mirabilis.