Protein diversification through post-translational modifications, alternative splicing, and gene duplication

Selenoprotein K at the intersection of cellular pathways

Selenoprotein K (selenok) is linked to the integrated stress response, which helps cells combat stressors and regain normal function. The selenoprotein contains numerous protein interaction hubs and post-translational modification sites and is involved in protein palmitoylation, vesicle trafficking, and the resolution of ER stress. Anchored to the endoplasmic reticulum (ER) membrane, selenok interacts with protein partners to influence their stability, localization, and trafficking, impacting various cellular functions such as calcium homeostasis, cellular migration, phagocytosis, gene expression, and immune response. Consequently, selenok expression level is linked to cancer and neurodegenerative diseases. Because it contains the reactive amino acid selenocysteine, selenok is likely to function as an enzyme. However, highly unusual for enzymes, the protein segment containing the selenocysteine lacks a stable secondary or tertiary structure, yet it includes multiple interaction sites for protein partners and post-translational modifications. Currently, the reason(s) for the presence of the rare selenocysteine in selenok is not known. Furthermore, of selenok's numerous interaction sites, only some have been sufficiently characterized, leaving many of selenok's potential protein partners to be discovered. In this review, we explore selenok's role in various cellular pathways and its impact on human health, thereby highlighting the links between its diverse cellular functions.

The Role of E3 Ubiquitin Ligase Gene FBK in Ubiquitination Modification of Protein and Its Potential Function in Plant Growth, Development, Secondary Metabolism, and Stress Response

As a crucial post-translational modification (PTM), protein ubiquitination mediates the breakdown of particular proteins, which plays a pivotal role in a large number of biological processes including plant growth, development, and stress response. The ubiquitin-proteasome system (UPS) consists of ubiquitin (Ub), ubiquitinase, deubiquitinating enzyme (DUB), and 26S proteasome mediates more than 80% of protein degradation for protein turnover in plants. For the ubiquitinases, including ubiquitin-activating enzyme (E1), ubiquitin-conjugating enzyme (E2), and ubiquitin ligase (E3), the FBK (F-box Kelch repeat protein) is an essential component of multi-subunit E3 ligase SCF (Skp1-Cullin 1-F-box) involved in the specific recognition of target proteins in the UPS. Many FBK genes have been identified in different plant species, which regulates plant growth and development through affecting endogenous phytohormones as well as plant tolerance to various biotic and abiotic stresses associated with changes in secondary metabolites such as phenylpropanoid, phenolic acid, flavonoid, lignin, wax, etc. The review summarizes the significance of the ubiquitination modification of protein, the role of UPS in protein degradation, and the possible function of FBK genes involved in plant growth, development, secondary metabolism, and stress response, which provides a systematic and comprehensive understanding of the mechanism of ubiquitination and potential function of FBKs in plant species.

Exploring gene expression, alternative splicing events and RNA-binding proteins changes in PBMC from patients with hyperuricemia

AIM

The objective of this study was to examine the transcriptomic profile changes in hyperuricemia (HUA) and to investigate the pathogenic mechanisms and biomarkers of HUA from a transcriptomic perspective.

METHODS

In this study, three patients with HUA were randomly selected and matched with three healthy controls. Six participants provided peripheral blood mononuclear cells (PBMCs) for analysis. RNA sequencing (RNA-seq) was used to identify differentially expressed genes (DEGs) and alternative splicing events (ASEs). Gene Ontology (GO) biological processes and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to identify the functions and pathways of the DEGs and ASEs. Additionally, a co-expression network was constructed to analyze the regulation of DEGs and ASEs by RNA-binding protein (RBP) genes. In addition, important DEGs and ASEs were validated using quantitative real-time PCR (qPCR).

RESULTS

There were 633 DEGs identified, 348 up-regulated DEGs and 285 down-regulated DEGs, including RGS18, CAVIN2, GZMH, GNLY and MT-TV, which were mainly enriched in inflammatory and immune-related biological processes. A total of 1542 ASEs were significantly differentially expressed in HUA, of which LTB4R and ENTPD4 were closely associated with the development of HUA. In addition, 15 RBP genes were detected to be differentially expressed in HUA. Three RBP genes (IFIT1, IFFIT2, and IFIT3) were highly associated with immunoinflammation and affected HUA by modulating downstream immune responses, inflammatory response-associated DEGs, and ASEs. The selected five DEGs and two ASEs were verified by qPCR, which was consistent with the results of RNA sequencing.

CONCLUSIONS

In summary, the findings indicate that HUA is associated with significant changes in inflammatory and immune response-related genes (RGS18, CAVIN2, GZMH, GNLY, MT-TV, LTB4R, ENTPD4, IFIT1, IFFIT2, and IFIT3). These findings suggest potential biomarkers and therapeutic targets.

Identifying Calmodulin and Calmodulin-like Protein Members in Canavalia rosea and Exploring Their Potential Roles in Abiotic Stress Tolerance

Calmodulins (CaMs) and calmodulin-like proteins (CMLs) belong to families of calcium-sensors that act as calcium ion (Ca2+) signal-decoding proteins and regulate downstream target proteins. As a tropical halophyte, Canavalia rosea shows great resistance to multiple abiotic stresses, including high salinity/alkalinity, extreme drought, heat, and intense sunlight. However, investigations of calcium ion signal transduction involved in the stress responses of C. rosea are limited. The CaM and CML gene families have been identified and characterized in many other plant species. Nevertheless, there is limited available information about these genes in C. rosea. In this study, a bioinformatic analysis, including the gene structures, conserved protein domains, phylogenetic relationships, chromosome distribution, and gene synteny, was comprehensively performed to identify and characterize CrCaMs and CrCMLs. A spatio-temporal expression assay in different organs and environmental conditions was then conducted using the RNA sequencing technique. Additionally, several CrCaM and CrCML members were then cloned and functionally characterized using the yeast heterogeneous expression system, and some of them were found to change the tolerance of yeast to heat, salt, alkalinity, and high osmotic stresses. The results of this study provide a foundation for understanding the possible roles of the CrCaM and CrCML genes, especially for halophyte C. rosea's natural ecological adaptability for its native habitats. This study also provides a theoretical basis for further study of the physiological and biochemical functions of plant CaMs and CMLs that are involved in tolerance to multiple abiotic stresses.

Genome-Wide Identification and Expression Analysis of the Cys2His2 Zinc Finger Protein Gene Family in Flammulina filiformis

Zinc finger proteins (ZFPs) are essential transcription factors in eukaryotes, particularly the extensively studied C2H2 family, which is known for its involvement in various biological processes. This research provides a thorough examination and analysis of the C2H2-ZFP gene family in Flammulina filiformis. Using bioinformatics tools, 58 FfC2H2-ZFP genes spread across 11 chromosomes were identified and scrutinized in detail for their gene structures, protein characteristics, and phylogenetic relationships. The study of phylogenetics and synteny sheds light on the evolutionary relationships among C2H2-ZFPs in F. filiformis and other fungi, revealing a complex evolutionary past. The identification of conserved cis-regulatory elements in the gene promoter regions suggests intricate functionalities, particularly in the developmental and stress response pathways. By utilizing RNA-seq and qRT-PCR techniques, the expression patterns of these genes were explored across different developmental stages and tissues of F. filiformis, unveiling distinct expression profiles. Notably, significant expression variations were observed in the stipe elongation region and pilei of various sizes, indicating potential roles in fruiting body morphogenesis. This study enhances our knowledge of the C2H2-ZFP gene family in F. filiformis and lays the groundwork for future investigations into their regulatory mechanisms and applications in fungal biology and biotechnology.

The role of protein post-translational modifications in prostate cancer

Involving addition of chemical groups or protein units to specific residues of the target protein, post-translational modifications (PTMs) alter the charge, hydrophobicity, and conformation of a protein, which in turn influences protein function, protein-protein interaction, and protein aggregation. These alterations, which include phosphorylation, glycosylation, ubiquitination, methylation, acetylation, lipidation, and lactylation, are significant biological events in the development of cancer, and play vital roles in numerous biological processes. The processes behind essential functions, the screening of clinical illness signs, and the identification of therapeutic targets all depend heavily on further research into the PTMs. This review outlines the influence of several PTM types on prostate cancer (PCa) diagnosis, therapy, and prognosis in an effort to shed fresh light on the molecular causes and progression of the disease.

Genetic and selective constraints on the optimization of gene product diversity

RNA and protein expressed from the same gene can have diverse isoforms due to various post-transcriptional and post-translational modifications. For the vast majority of alternative isoforms, It is unknown whether they are adaptive or simply biological noise. As we cannot experimentally probe the function of each isoform, we can ask whether the distribution of isoforms across genes and across species is consistent with expectations from different evolutionary processes. However, there is currently no theoretical framework that can generate such predictions. To address this, we developed a mathematical model where isoform abundances are determined collectively by cis-acting loci, trans-acting factors, gene expression levels, and isoform decay rates to predict isoform abundance distributions across species and genes in the face of mutation, genetic drift, and selection. We found that factors beyond selection, such as effective population size and the number of cis-acting loci, significantly influence evolutionary outcomes. Notably, suboptimal phenotypes are more likely to evolve when the population is small and/or when the number of cis-loci is large. We also explored scenarios where modification processes have both beneficial and detrimental effects, revealing a non-monotonic relationship between effective population size and optimization, demonstrating how opposing selection pressures on cis- and trans-acting loci can constrain the optimization of gene product diversity. As a demonstration of the power of our theory, we compared the expected distribution of A-to-I RNA editing levels in coleoids and found this to be largely consistent with non-adaptive explanations.

A systematic analysis of the effects of splicing on the diversity of post-translational modifications in protein isoforms

Post-translational modifications (PTMs) and splicing are known to be important regulatory processes for controlling protein function and activity. However, there have been limitations in analyzing the interplay of alternative splicing and PTMs, which stems from the deep differences in genomic and proteomic databases. In this work, we bridged the protein- and genome-centric world views to map PTMs to genomic locations for subsequent projection of PTMs onto alternative isoforms. We then performed a systematic analysis of the diversification of PTMs by alternative splicing, including exploration of the modification-specific rates of inclusion across isoforms and how often the regulatory sequences directly flanking a PTM are impacted by splicing, which might indicate altered regulatory or binding interactions in the alternatively spliced isoform. We found that 6-51% of PTMs are excluded from at least one isoform, depending on the modification type. Further, approximately 2% of prospective PTM sites exhibited altered regulatory sequences surrounding the modification site, suggesting that regulatory or binding interactions might be diversified in these proteoforms. Lastly, we applied this PTM-to-isoform mapping approach to explore the impacts of disease-related splicing in prostate cancer, identifying possible new hypotheses explaining the variable consequences of ESRP1 expression in different cancers. As a part of this work, we have provided an easily implementable tool for annotating splice events identified from RNA-sequencing with PTMs and their functional consequences, called PTM-POSE.

Transcript Isoform Diversity of Y Chromosome Ampliconic Genes of Great Apes Uncovered Using Long Reads and Telomere-to-Telomere Reference Genome Assemblies

Y chromosomes of great apes harbor Ampliconic Genes (YAGs)-multi-copy gene families (BPY2, CDY, DAZ, HSFY, PRY, RBMY, TSPY, VCY, and XKRY) that encode proteins important for spermatogenesis. Previous work assembled YAG transcripts based on their targeted sequencing but not using reference genome assemblies, potentially resulting in an incomplete transcript repertoire. Here we used the recently produced gapless telomere-to-telomere (T2T) Y chromosome assemblies of great ape species (bonobo, chimpanzee, human, gorilla, Bornean orangutan, and Sumatran orangutan) and analyzed RNA data from whole-testis samples for the same species. We generated hybrid transcriptome assemblies by combining targeted long reads (Pacific Biosciences), untargeted long reads (Pacific Biosciences) and untargeted short reads (Illumina)and mapping them to the T2T reference genomes. Compared to the results from the reference-free approach, average transcript length was more than two times higher, and the total number of transcripts decreased three times, improving the quality of the assembled transcriptome. The reference-based transcriptome assemblies allowed us to differentiate transcripts originating from different Y chromosome gene copies and from their non-Y chromosome homologs. We identified two sources of transcriptome diversity-alternative splicing and gene duplication with subsequent diversification of gene copies. For each gene family, we detected transcribed pseudogenes along with protein-coding gene copies. We revealed previously unannotated gene copies of YAGs as compared to currently available NCBI annotations, as well as novel isoforms for annotated gene copies. This analysis paves the way for better understanding Y chromosome gene functions, which is important given their role in spermatogenesis.

Genome-wide analysis of the C2H2-ZFP gene family in Stevia rebaudiana reveals involvement in abiotic stress response

Stevia (Stevia rebaudiana Bertoni) is a natural sweetener plant that accumulates highly sweet steviol glycosides (SGs) especially in leaves. Stevia is native to humid areas and does not have a high tolerance to drought which is the most serious abiotic stress restricting its production worldwide. C2H2 zinc finger proteins (C2H2-ZFPs) are a group of well-known transcription factors that involves in various developmental, physiological and biochemical activities as well as in response to abiotic stresses. Here we analyzed C2H2-ZFP gene family in stevia and identified a total of 185 putative SrC2H2-ZF proteins from the genome sequence of S. rebaudiana. We further characterized the identified C2H2-ZF domains and their organization, additional domains and motifs and analyzed their physicochemical properties, localization and gene expression patterns. The cis-element analysis suggested multiple roles of SrC2H2-ZFPs in response to light, phytohormone, and abiotic stresses. In silico analysis revealed that the stevia C2H2-ZFP genes are interactively expressed in different tissues and developmental stages and some C2H2-ZFP genes are involved in response to drought stress. This study provides a background for future exploration of the functional, and regulatory aspects of the C2H2-ZFP gene family in S. rebaudiana.

Lectin diversity and their positive roles in WSSV replication through regulation of calreticulin expression and inhibiting ALFs expression

In biological evolution, gene duplication (GD) generates new genes to facilitate new functions. C-type lectins (CTLs) in crayfish have been extended by GD to expand their family members. In this study, four CTL genes generated by GD were identified from Procambarus clarkii (PcLec1-4). Among these four genes, PcLec1 can also generate new isoforms with different numbers of tandem repeats through DNA slip mispairing. PcLec1-4 was widely expressed in multiple tissues. The expression levels of PcLec1-4 were upregulated in the intestine of P. clarkii upon white spot syndrome virus (WSSV) challenge at multiple time points. Further analysis indicated that GATA transcription factor regulated PcLec1-4 expression. RNA interference and recombinant PcLec1-4 protein injection experiments suggested that PcLec1-4 promoted the expression of calreticulin (PcCRT) and negatively regulated the expression of antimicrobial peptides, thereby promoting WSSV replication. This study contributes to the understanding of the function of CTLs produced by GD during WSSV invasion in crustaceans.

The Biological Action and Structural Characterization of Eryngitin 3 and 4, Ribotoxin-like Proteins from Pleurotus eryngii Fruiting Bodies

Ribotoxin-like proteins (RL-Ps) are specific ribonucleases found in mushrooms that are able to cleave a single phosphodiester bond located in the sarcin-ricin loop (SRL) of the large rRNA. The cleaved SRL interacts differently with some ribosomal proteins (P-stalk). This action blocks protein synthesis because the damaged ribosomes are unable to interact with elongation factors. Here, the amino acid sequences of eryngitin 3 and 4, RL-Ps isolated from Pleurotus eryngii fruiting bodies, were determined to (i) obtain structural information on this specific ribonuclease family from edible mushrooms and (ii) explore the structural determinants which justify their different biological and antipathogenic activities. Indeed, eryngitin 3 exhibited higher toxicity with respect to eryngitin 4 against tumoral cell lines and model fungi. Structurally, eryngitin 3 and 4 consist of 132 amino acids, most of them identical and exhibiting a single free cysteinyl residue. The amino acidic differences between the two toxins are (i) an additional phenylalanyl residue at the N-terminus of eryngitin 3, not retrieved in eryngitin 4, and (ii) an additional arginyl residue at the C-terminus of eryngitin 4, not retrieved in eryngitin 3. The 3D models of eryngitins show slight differences at the N- and C-terminal regions. In particular, the positive electrostatic surface at the C-terminal of eryngitin 4 is due to the additional arginyl residue not retrieved in eryngitin 3. This additional positive charge could interfere with the binding to the SRL (substrate) or with some ribosomal proteins (P-stalk structure) during substrate recognition.

Identification of ARF transcription factor gene family and its defense responses to bacterial infection and salicylic acid treatment in sugarcane

Auxin response factor (ARF) is a critical regulator in the auxin signaling pathway, involved in a variety of plant biological processes. Here, gene members of 24 SpapARFs and 39 SpnpARFs were identified in two genomes of Saccharum spontaneum clones AP85-441 and Np-X, respectively. Phylogenetic analysis showed that all ARF genes were clustered into four clades, which is identical to those ARF genes in maize (Zea mays) and sorghum (Sorghum bicolor). The gene structure and domain composition of this ARF family are conserved to a large degree across plant species. The SpapARF and SpnpARF genes were unevenly distributed on chromosomes 1-8 and 1-10 in the two genomes of AP85-441 and Np-X, respectively. Segmental duplication events may also contribute to this gene family expansion in S. spontaneum. The post-transcriptional regulation of ARF genes likely involves sugarcane against various stressors through a miRNA-medicated pathway. Expression levels of six representative ShARF genes were analyzed by qRT-PCR assays on two sugarcane cultivars [LCP85-384 (resistant to leaf scald) and ROC20 (susceptible to leaf scald)] triggered by Acidovorax avenae subsp. avenae (Aaa) and Xanthomonas albilineans (Xa) infections and salicylic acid (SA) treatment. ShARF04 functioned as a positive regulator under Xa and Aaa stress, whereas it was a negative regulator under SA treatment. ShARF07/17 genes played positive roles against both pathogenic bacteria and SA stresses. Additionally, ShARF22 was negatively modulated by Xa and Aaa stimuli in both cultivars, particularly LCP85-384. These findings imply that sugarcane ARFs exhibit functional redundancy and divergence against stressful conditions. This work lays the foundation for further research on ARF gene functions in sugarcane against diverse environmental stressors.