Comprehensive genomic analysis of the DUF4228 gene family in land plants and expression profiling of ATDUF4228 under abiotic stresses

AtPADRE13 Negatively Regulates Salt Stress Tolerance in Arabidopsis thaliana

The PADRE (Pathogen and abiotic stress response, cadmium tolerance, disordered region-containing) family of genes, which contains the structural DUF4228 domain of unknown function (DUF), has been reported to be associated with plant responses to abiotic stresses. However, the specific functions of this family in the salt stress response remain unknown. AtPADRE13 is induced by salt stress and ABA (abscisic acid). After the overexpression of AtPADRE13 in Arabidopsis, seeds were found to be insensitive to ABA treatment. After salt treatment, the overexpression lines presented a significantly lower survival rate, increased MDA (Malondialdehyde) content, and reduced antioxidant enzyme activities compared with the wild-type, and were more sensitive to salt stress. Transcriptome data analysis further revealed that AtPADRE13 overexpression resulted in different degrees of down-regulation for a series of positive regulators related to ABA catabolism, transport, and their mediated plant responses to salt stress. In addition, the expression of genes related to ROS (reactive oxygen species) scavenging was down-regulated. In conclusion, AtPADRE13 plays a negative regulatory role in the response to salt stress in Arabidopsis.

Dietary resveratrol improves the flesh quality of Siberian sturgeon (Acipenser baerii) by enhancing myofiber growth, nutrient accumulation and antioxidant capacity

BACKGROUND

In aquaculture, sturgeons are generally maintained in the confined spaces, which not only hinders sturgeon movement, but also threatens their flesh quality that seriously concerned by aquaculture industry. As a typical antioxidant, resveratrol can improve the flesh quality of livestock and poultry. However, the mechanism of resveratrol's effect on the muscle of Siberian sturgeon is still unclear.

RESULTS

In this study, the dietary resveratrol increased the myofiber diameter, the content of the amino acids, antioxidant capacity markers (CAT, LDH and SOD) levels and the expression levels of mTORC1 and MYH9 in muscle of Siberian sturgeon. Further transcriptome analysis displayed that ROS production-related pathways ("Oxidative phosphorylation" and "Chemical carcinogenes-reactive oxygen species") were enriched in KEGG analysis, and the expression levels of genes related to the production of ROS (COX4, COX6A, ATPeF1A, etc.) in mitochondria were significantly down-regulated, while the expression levels of genes related to scavenging ROS (SOD1) were up-regulated.

CONCLUSIONS

In summary, this study reveals that resveratrol may promote the flesh quality of Siberian sturgeon probably by enhancing myofiber growth, nutritional value and the antioxidant capacity of muscle, which has certain reference significance for the development of a new type of feed for Siberian sturgeon.

Domain of unknown function (DUF) proteins in plants: function and perspective

Domains of unknown function (DUFs), which are deposited in the protein family database (Pfam), are protein domains with conserved amino acid sequences and uncharacterized functions. Proteins with the same DUF were classified as DUF families. Although DUF families are generally not essential for the survival of plants, they play roles in plant development and adaptation. Characterizing the functions of DUFs is important for deciphering biological puzzles. DUFs were generally studied through forward and reverse genetics. Some novelty approaches, especially the determination of crystal structures and interaction partners of the DUFs, should attract more attention. This review described the identification of DUF genes by genome-wide and transcriptome-wide analyses, summarized the function of DUF-containing proteins, and addressed the prospects for future studies in DUFs in plants.

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

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

Genomic wide association study and selective sweep analysis identify genes associated with improved yield under drought in Turkish winter wheat germplasm

A panel comprising of 84 Turkish winter wheat landraces (LR) and 73 modern varieties (MV) was analyzed with genome wide association study (GWAS) to identify genes/genomic regions associated with increased yield under favorable and drought conditions. In addition, selective sweep analysis was conducted to detect signatures of selection in the winter wheat genome driving the differentiation between LR and MV, to gather an understanding of genomic regions linked to adaptation and yield improvement. The panel was genotyped with 25 K wheat SNP array and phenotyped for agronomic traits for two growing seasons (2018 and 2019) in Konya, Turkey. Year 2018 was treated as drought environment due to very low precipitation prior to heading whereas year 2019 was considered as a favorable season. GWAS conducted with SNPs and haplotype blocks using mixed linear model identified 18 genomic regions in the vicinities of known genes i.e., TaERF3-3A, TaERF3-3B, DEP1-5A, FRIZZY PANICLE-2D, TaSnRK23-1A, TaAGL6-A, TaARF12-2A, TaARF12-2B, WAPO1, TaSPL16-7D, TaTGW6-A1, KAT-2B, TaOGT1, TaSPL21-6B, TaSBEIb, trs1/WFZP-A, TaCwi-A1-2A and TaPIN1-7A associated with grain yield (GY) and yield related traits. Haplotype-based GWAS identified five haplotype blocks (H1A-42, H2A-71, H4A-48, H7B-123 and H7B-124), with the favorable haplotypes showing a yield increase of > 700 kg/ha in the drought season. SNP-based GWAS, detected only one larger effect genomic region on chromosome 7B, in common with haplotype-based GWAS. On an average, the percentage variation (PV) explained by haplotypes was 8.0% higher than PV explained by SNPs for all the investigated traits. Selective sweep analysis detected 39 signatures of selection between LR and MV of which 15 were within proximity of known functional genes controlling flowering (PRR-A1, PPR-D1, TaHd1-6B), GY and GY components (TaSus2-2B, TaGS2-B1, AG1-1A/WAG1-1A, DUO-A1, DUO-B1, AG2-3A/WAG2-3A, TaLAX1, TaSnRK210-4A, FBP, TaLAX1, TaPIL1 and AP3-1-7A/WPA3-7A) and 10 regions underlying various transcription factors and regulatory genes. The study outcomes contribute to utilization of LR in breeding winter wheat.

Genome-Wide Identification, Characterization and Expression Analysis of the TaDUF724 Gene Family in Wheat (Triticum aestivum)

Unknown functional domain (DUF) proteins constitute a large number of functionally uncharacterized protein families in eukaryotes. DUF724s play crucial roles in plants. However, the insight understanding of wheat TaDUF724s is currently lacking. To explore the possible function of TaDUF724s in wheat growth and development and stress response, the family members were systematically identified and characterized. In total, 14 TaDUF724s were detected from a wheat reference genome; they are unevenly distributed across the 11 chromosomes, and, according to chromosome location, they were named TaDUF724-1 to TaDUF724-14. Evolution analysis revealed that TaDUF724s were under negative selection, and fragment replication was the main reason for family expansion. All TaDUF724s are unstable proteins; most TaDUF724s are acidic and hydrophilic. They were predicted to be located in the nucleus and chloroplast. The promoter regions of TaDUF724s were enriched with the cis-elements functionally associated with growth and development, as well as being hormone-responsive. Expression profiling showed that TaDUF724-9 was highly expressed in seedings, roots, leaves, stems, spikes and grains, and strongly expressed throughout the whole growth period. The 12 TaDUF724 were post-transcription regulated by 12 wheat MicroRNA (miRNA) through cleavage and translation. RT-qPCR showed that six TaDUF724s were regulated by biological and abiotic stresses. Conclusively, TaDUF724s were systematically analyzed using bioinformatics methods, which laid a theoretical foundation for clarifying the function of TaDUF724s in wheat.

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

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

Unraveling the Diverse Roles of Neglected Genes Containing Domains of Unknown Function (DUFs): Progress and Perspective

Domain of unknown function (DUF) is a general term for many uncharacterized domains with two distinct features: relatively conservative amino acid sequence and unknown function of the domain. In the Pfam 35.0 database, 4795 (24%) gene families belong to the DUF type, yet, their functions remain to be explored. This review summarizes the characteristics of the DUF protein families and their functions in regulating plant growth and development, generating responses to biotic and abiotic stress, and other regulatory roles in plant life. Though very limited information is available about these proteins yet, by taking advantage of emerging omics and bioinformatic tools, functional studies of DUF proteins could be utilized in future molecular studies.

Overexpression of a DUF740 family gene (LOC_Os04g59420) imparts enhanced climate resilience through multiple stress tolerance in rice

Functional characterization of stress-responsive genes through the analysis of transgenic plants is a standard approach to comprehend their role in climate resilience and subsequently exploit them for sustainable crop improvement. In this study, we investigated the function of LOC_Os04g59420, a gene of DUF740 family (OsSRDP-Oryza sativa Stress Responsive DUF740 Protein) from rice, which showed upregulation in response to abiotic stress in the available global expression data, but is yet to be functionally characterized. Transgenic plants of the rice OsSRDP gene, driven by a stress-inducible promoter AtRd29A, were developed in the background of cv. Pusa Sugandh 2 (PS2) and their transgene integration and copy number were confirmed by molecular analysis. The three independent homozygous transgenic plants (AtRd29A::OsSRDP rice transformants) showed better resilience to drought, salinity, and cold stresses, but not heat stress, as compared to the non-transformed PS2, which corresponded with their respective relative transcript abundance for OsSRDP. Transgenic plants maintained higher RWC, photosynthetic pigments, and proline accumulation under drought and salinity stresses. Furthermore, they exhibited less accumulation of reactive oxygen species (ROS) than PS2 under drought stress, as seen from the transcript abundance studies of the ROS genes. Under cold stress, OsSRDP transgenic lines illustrated minimal cell membrane injury compared to PS2. Additionally, the transgenic plants showed resistance to a virulent strain of rice blast fungus, Magnaporthe oryzae (M. oryzae). The promoter analysis of the gene in N22 and PS2 revealed the presence of multiple abiotic and biotic stress-specific motif elements supporting our observation on multiple stress tolerance. Based on bioinformatics studies, we identified four potential candidate interaction partners for LOC_Os04g59420, of which two genes (LOC_Os05g09640 and LOC_Os06g50370) showed co-expression under biotic and drought stress along with OsSRDP. Altogether, our findings established that stress-inducible expression of OsSRDP can significantly enhance tolerance to multiple abiotic stresses and a biotic stress.

Genome-Wide Identification of DUF668 Gene Family and Expression Analysis under Drought and Salt Stresses in Sweet Potato [Ipomoea batatas (L.) Lam]

The domain of unknown function 668 (DUF668) is a gene family that plays a vital role in responses to adversity coercion stresses in plant. However, the function of the DUF668 gene family is not fully understood in sweet potato. In this study, bioinformatics methods were used to analyze the number, physicochemical properties, evolution, structure, and promoter cis-acting elements of the IbDUF668 family genes, and RNA-seq and qRT-PCR were performed to detect gene expression and their regulation under hormonal and abiotic stress. A total of 14 IbDUF668 proteins were identified in sweet potato, distributed on nine chromosomes. By phylogenetic analysis, IbDUF668 proteins can be divided into two subfamilies. Transcriptome expression profiling revealed that many genes from DUF668 in sweet potato showed specificity and differential expression under cold, heat, drought, salt and hormones (ABA, GA3 and IAA). Four genes (IbDUF668-6, 7, 11 and 13) of sweet potato were significantly upregulated by qRT-PCR under ABA, drought and NaCl stress. Results suggest that the DUF668 gene family is involved in drought and salt tolerance in sweet potato, and it will further provide the basic information of DUF668 gene mechanisms in plants.

A Comprehensive Analysis of the DUF4228 Gene Family in Gossypium Reveals the Role of GhDUF4228-67 in Salt Tolerance

Soil salinization conditions seriously restrict cotton yield and quality. Related studies have shown that the DUF4228 proteins are pivotal in plant resistance to abiotic stress. However, there has been no systematic identification and analysis of the DUF4228 gene family in cotton and their role in abiotic stress. In this study, a total of 308 DUF4228 genes were identified in four Gossypium species, which were divided into five subfamilies. Gene structure and protein motifs analysis showed that the GhDUF4228 proteins were conserved in each subfamily. In addition, whole genome duplication (WGD) events and allopolyploidization might play an essential role in the expansion of the DUF4228 genes. Besides, many stress-responsive (MYB, MYC) and hormone-responsive (ABA, MeJA) related cis-elements were detected in the promoters of the DUF4228 genes. The qRT-PCR results showed that GhDUF4228 genes might be involved in the response to abiotic stress. VIGS assays and the measurement of relative water content (RWC), Proline content, POD activity, and malondialdehyde (MDA) content indicated that GhDUF4228-67 might be a positive regulator of cotton response to salt stress. The results in this study systematically characterized the DUF4228s in Gossypium species and will provide helpful information to further research the role of DUF4228s in salt tolerance.

Genome-Wide Identification and Expression Profiling of DUF221 Gene Family Provides New Insights Into Abiotic Stress Responses in Potato

The domain of the unknown function 221 proteins regulate several processes in plants, including development, growth, hormone transduction mechanism, and abiotic stress response. Therefore, a comprehensive analysis of the potato genome was conducted to identify the deafness-dystonia peptide (DDP) proteins' role in potatoes. In the present study, we performed a genome-wide analysis of the potato domain of the unknown function 221 (DUF221) genes, including phylogenetic inferences, chromosomal locations, gene duplications, gene structures, and expression analysis. In our results, we identified 10 DDP genes in the potato genome. The phylogenetic analysis results indicated that StDDPs genes were distributed in all four clades, and clade IV was the largest clade. The gene duplication under selection pressure analysis indicated various positive and purifying selections in StDDP genes. The putative stu-miRNAs from different families targeting StDDPs were also predicted in the present study. Promoter regions of StDDP genes contain different cis-acting components involved in multiple stress responses, such as phytohormones and abiotic stress-responsive factors. The analysis of the tissue-specific expression profiling indicated the StDDPs gene expression in stem, root, and leaf tissues. We subsequently observed that StDDP4, StDDP5, and StDDP8 showed higher expressions in roots, stems, and leaves. StDDP5 exhibited high expression against heat stress response, and StDDP7 showed high transcript abundance against salt stress in potatoes. Under abscisic acid (ABA) and indole acetic acid (IAA) treatments, seven StDDP genes' expressions indicated that ABA and IAA performed important roles in immunity response. The expression profiling and real-time qPCR of stems, roots, and leaves revealed StDDPs' significant role in growth and development. These expression results of DDPs are primary functional analysis and present basic information for other economically important crops.

Genome-Wide Identification and Transcriptional Analysis of Arabidopsis DUF506 Gene Family

The Domain of unknown function 506 (DUF506) family, which belongs to the PD-(D/E)XK nuclease superfamily, has not been functionally characterized. In this study, 266 DUF506 domain-containing genes were identified from algae, mosses, and land plants showing their wide occurrence in photosynthetic organisms. Bioinformatics analysis identified 211 high-confidence DUF506 genes across 17 representative land plant species. Phylogenetic modeling classified three groups of plant DUF506 genes that suggested functional preservation among the groups based on conserved gene structure and motifs. Gene duplication and Ka/Ks evolutionary rates revealed that DUF506 genes are under purifying positive selection pressure. Subcellular protein localization analysis revealed that DUF506 proteins were present in different organelles. Transcript analyses showed that 13 of the Arabidopsis DUF506 genes are ubiquitously expressed in various tissues and respond to different abiotic stresses and ABA treatment. Protein-protein interaction network analysis using the STRING-DB, AtPIN (Arabidopsis thaliana Protein Interaction Network), and AI-1 (Arabidopsis Interactome-1) tools indicated that AtDUF506s potentially interact with iron-deficiency response proteins, salt-inducible transcription factors, or calcium sensors (calmodulins), implying that DUF506 genes have distinct biological functions including responses to environmental stimuli, nutrient-deficiencies, and participate in Ca(2+) signaling. Current results provide insightful information regarding the molecular features of the DUF506 family in plants, to support further functional characterizations.

Identification of RNA helicases in Medicago truncatula and their expression patterns under abiotic stress

RNA helicase catalyzes the denaturation of DNA or the unwinding of double-stranded RNA. It is vital to RNA splicing, transport, editing, degradation and the initiation of protein translation. However, the function of RNA helicase in Medicago truncatula has rarely been reported. In this study, 170 putative RNA helicase genes were identified in the M. truncatula genome, and classified into three subfamilies based on the presence of either a DEAD-box (52 genes), DEAH-box (38 genes), or DExD/H-box (80 genes) in their coding regions. Additionally, conserved helicase_C domains and other functional domains (e.g., the HA2, DUF, and ZnF domains) were also present in these genes. Chromosomal mapping and synteny analyses showed that there were tandem and segment duplications of RNA helicase genes. Furthermore, transcriptome and real-time PCR analysis showed that the expression of 35 RNA helicase genes was affected by abiotic stress. To be specific, 17, 12 and 19 genes were regulated by salt, drought and cold stress, respectively. It is worth noting that MtDEAD8, MtDEAH3, MtDExD/H18 and MtDExD/H23 responded to all three types of stress. These results provide valuable information for understanding the RNA helicase genes in M. truncatula and their abiotic stress-related functions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01087-y.

Genome-Wide Analysis of the DUF4228 Family in Soybean and Functional Identification of GmDUF4228 -70 in Response to Drought and Salt Stresses

Domain of unknown function 4228 (DUF4228) proteins are a class of proteins widely found in plants, playing an important role in response to abiotic stresses. However, studies on the DUF4228 family in soybean (Glycine max L.) are sparse. In this study, we identified a total of 81 DUF4228 genes in soybean genome, named systematically based on their chromosome distributions. Results showed that these genes were unevenly distributed on the 20 chromosomes of soybean. The predicted soybean DUF4228 proteins were identified in three groups (Groups I-III) based on a maximum likelihood phylogenetic tree. Genetic structure analysis showed that most of the GmDUF4228 genes contained no introns. Expression profiling showed that GmDUF4228 genes were widely expressed in different organs and tissues in soybean. RNA-seq data were used to characterize the expression profiles of GmDUF4228 genes under the treatments of drought and salt stresses, with nine genes showing significant up-regulation under both drought and salt stress further functionally verified by promoter (cis-acting elements) analysis and quantitative real-time PCR (qRT-PCR). Due to its upregulation under drought and salt stresses based on both RNA-seq and qRT-PCR analyses, GmDUF4228-70 was selected for further functional analysis in transgenic plants. Under drought stress, the degree of leaf curling and wilting of the GmDUF4228-70-overexpressing (GmDUF4228-70-OE) line was lower than that of the empty vector (EV) line. GmDUF4228-70-OE lines also showed increased proline content, relative water content (RWC), and chlorophyll content, and decreased contents of malondialdehyde (MDA), H2O2, and O2-. Under salt stress, the changes in phenotypic and physiological indicators of transgenic plants were the same as those under drought stress. In addition, overexpression of the GmDUF4228-70 gene promoted the expression of marker genes under both drought and salt stresses. Taken together, the results indicated that GmDUF4228 genes play important roles in response to abiotic stresses in soybean.

Characterization, identification and expression profiling of genome-wide R-genes in melon and their putative roles in bacterial fruit blotch resistance

Background

Bacterial fruit blotch (BFB), a disease caused by Acidovorax citrulli, results in significant economic losses in melon. The causal QTLs and genes for resistance to this disease have yet to be identified. Resistance (R)-genes play vital roles in resistance to plant diseases. Since the complete genome sequence of melon is available and genome-wide identification of R-genes has been performed for this important crop, comprehensive expression profiling may lead to the identification of putative candidate genes that function in the response to BFB.

Results

We identified melon accessions that are resistant and susceptible to BFB through repeated bioassays and characterized all 70 R-genes in melon, including their gene structures, chromosomal locations, domain organizations, motif distributions, and syntenic relationships. Several disease resistance-related domains were identified, including NBS, TIR, LRR, CC, RLK, and DUF domains, and the genes were categorized based on the domains of their encoded proteins. In addition, we profiled the expression patterns of the genes in melon accessions with contrasting levels of BFB resistance at 12 h, 1 d, 3 d, and 6 d after inoculation with A. citrulli. Six R-genes exhibited consistent expression patterns (MELO3C023441, MELO3C016529, MELO3C022157, MELO3C022146, MELO3C025518, and MELO3C004303), with higher expression levels in the resistant vs. susceptible accession.

Conclusion

We identified six putative candidate R-genes against BFB in melon. Upon functional validation, these genes could be targeted for manipulation via breeding and biotechnological approaches to improve BFB resistance in melon in the future.

Patterns of Sequence and Expression Diversification Associate Members of the PADRE Gene Family With Response to Fungal Pathogens

Pathogen infection triggers extensive reprogramming of the plant transcriptome, including numerous genes the function of which is unknown. Due to their wide taxonomic distribution, genes encoding proteins with Domains of Unknown Function (DUFs) activated upon pathogen challenge likely play important roles in disease. In Arabidopsis thaliana, we identified thirteen genes harboring a DUF4228 domain in the top 10% most induced genes after infection by the fungal pathogen Sclerotinia sclerotiorum. Based on functional information collected through homology and contextual searches, we propose to refer to this domain as the pathogen and abiotic stress response, cadmium tolerance, disordered region-containing (PADRE) domain. Genome-wide and phylogenetic analyses indicated that PADRE is specific to plants and diversified into 10 subfamilies early in the evolution of Angiosperms. PADRE typically occurs in small single-domain proteins with a bipartite architecture. PADRE N-terminus harbors conserved sequence motifs, while its C-terminus includes an intrinsically disordered region with multiple phosphorylation sites. A pangenomic survey of PADRE genes expression upon S. sclerotiorum inoculation in Arabidopsis, castor bean, and tomato indicated consistent expression across species within phylogenetic groups. Multi-stress expression profiling and co-expression network analyses associated AtPADRE genes with the induction of anthocyanin biosynthesis and responses to chitin and to hypoxia. Our analyses reveal patterns of sequence and expression diversification consistent with the evolution of a role in disease resistance for an uncharacterized family of plant genes. These findings highlight PADRE genes as prime candidates for the functional dissection of mechanisms underlying plant disease resistance to fungi.

Analyses on clustering of the conserved residues at protein-RNA interfaces and its application in binding site identification

Background

The maintenance of protein structural stability requires the cooperativity among spatially neighboring residues. Previous studies have shown that conserved residues tend to occur clustered together within enzyme active sites and protein-protein/DNA interfaces. It is possible that conserved residues form one or more local clusters in protein tertiary structures as it can facilitate the formation of functional motifs. In this work, we systematically investigate the spatial distributions of conserved residues as well as hot spot ones within protein-RNA interfaces.

Results

The analysis of 191 polypeptide chains from 160 complexes shows the polypeptides interacting with tRNAs evolve relatively rapidly. A statistical analysis of residues in different regions shows that the interface residues are often more conserved, while the most conserved ones are those occurring at protein interiors which maintain the stability of folded polypeptide chains. Additionally, we found that 77.8% of the interfaces have the conserved residues clustered within the entire interface regions. Appling the clustering characteristics to the identification of the real interface, there are 31.1% of cases where the real interfaces are ranked in top 10% of 1000 randomly generated surface patches. In the conserved clusters, the preferred residues are the hydrophobic (Leu, Ile, Met), aromatic (Tyr, Phe, Trp) and interestingly only one positively charged Arg residues. For the hot spot residues, 51.5% of them are situated in the conserved residue clusters, and they are largely consistent with the preferred residue types in the conserved clusters.

Conclusions

The protein-RNA interface residues are often more conserved than non-interface surface ones. The conserved interface residues occur more spatially clustered relative to the entire interface residues. The high consistence of hot spot residue types and the preferred residue types in the conserved clusters has important implications for the experimental alanine scanning mutagenesis study. This work deepens the understanding of the residual organization at protein-RNA interface and is of potential applications in the identification of binding site and hot spot residues.