RNA Recognition Motif Domain Research Articles

Transcriptome studies of the floodwater mosquito, Aedes vexans (Diptera: Culicidae) with potential as secondary vectors using Illumina HiSeq 4,000 sequencing

AbstractAedes vexans is the inland floodwater mosquito widely distributed in every continent excepting Antarctica and South America. They are opportunistic feeders preferring blood meal of larger animals including cattle, horses, deer, and humans. Further, the mosquito species is a compatible vector of several diseases, including West Nile virus and dog heartworm. In this study, we performed transcriptome characterization of Ae. vexans using Illumina HiSeq 4,000 sequencing and assembly of sequenced reads using Trinity. A total of 55,813,852 raw read and 54,630,771 clean reads (97.88% of raw reads) were obtained after Illumina paired‐end sequencing and pre‐processing steps. After Trinity de novo assembly, TransDecoder and TGICL clustering, a total of 37,111 unigenes were obtained. Out of the total unigenes count, 28,733, 17,893, 14,626, and 17,055 showed homologous matches against the PANM, UniGene, SwissProt, and KOG databases. A total of 9,483 unigenes were assigned to Gene Ontology (GO) terms, and 3,741 unigenes were mapped to 483 KEGG pathways. The zinc finger (C2H2‐type), reverse transcriptase, integrase (catalytic core), protein kinase, and RNA recognition motif domain among others showed as the top InterProScan domains. The obtained datasets serves as a basis for future studies towards understanding ecology, metabolism, and parasitism potential of Ae. vexans.

Translational regulation of Chk1 expression by eIF3a via interaction with the RNA-binding protein HuR.

eIF3a is a putative subunit of the eukaryotic translation initiation factor 3 complex. Accumulating evidence suggests that eIF3a may have a translational regulatory function by suppressing translation of a subset of mRNAs while accelerating that of other mRNAs. Albeit the suppression of mRNA translation may derive from eIF3a binding to the 5'-UTRs of target mRNAs, how eIF3a may accelerate mRNA translation remains unknown. In this study, we show that eIF3a up-regulates translation of Chk1 but not Chk2 mRNA by interacting with HuR, which binds directly to the 3'-UTR of Chk1 mRNA. The interaction between eIF3a and HuR occurs at the 10-amino-acid repeat domain of eIF3a and the RNA recognition motif domain of HuR. This interaction may effectively circularize Chk1 mRNA to form an end-to-end complex that has recently been suggested to accelerate mRNA translation. Together with previous findings, we conclude that eIF3a may regulate mRNA translation by directly binding to the 5'-UTR to suppress or by interacting with RNA-binding proteins at 3'-UTRs to accelerate mRNA translation.

Structural basis of UCUU RNA motif recognition by splicing factor RBM20.

The vertebrate splicing factor RBM20 (RNA binding motif protein 20) regulates protein isoforms important for heart development and function, with mutations in the gene linked to cardiomyopathy. Previous studies have identified the four nucleotide RNA motif UCUU as a common element in pre-mRNA targeted by RBM20. Here, we have determined the structure of the RNA Recognition Motif (RRM) domain from mouse RBM20 bound to RNA containing a UCUU sequence. The atomic details show that the RRM domain spans a larger region than initially proposed in order to interact with the complete UCUU motif, with a well-folded C-terminal helix encoded by exon 8 critical for high affinity binding. This helix only forms upon binding RNA with the final uracil, and removing the helix reduces affinity as well as specificity. We therefore find that RBM20 uses a coupled folding-binding mechanism by the C-terminal helix to specifically recognize the UCUU RNA motif.

An autoinhibitory intramolecular interaction proof-reads RNA recognition by the essential splicing factor U2AF2

The recognition of cis-regulatory RNA motifs in human transcripts by RNA binding proteins (RBPs) is essential for gene regulation. The molecular features that determine RBP specificity are often poorly understood. Here, we combined NMR structural biology with high-throughput iCLIP approaches to identify a regulatory mechanism for U2AF2 RNA recognition. We found that the intrinsically disordered linker region connecting the two RNA recognition motif (RRM) domains of U2AF2 mediates autoinhibitory intramolecular interactions to reduce nonproductive binding to weak Py-tract RNAs. This proofreading favors binding of U2AF2 at stronger Py-tracts, as required to define 3' splice sites at early stages of spliceosome assembly. Mutations that impair the linker autoinhibition enhance the affinity for weak Py-tracts result in promiscuous binding of U2AF2 along mRNAs and impact on splicing fidelity. Our findings highlight an important role of intrinsically disordered linkers to modulate RNA interactions of multidomain RBPs.

Roles of the RGG Domain and RNA Recognition Motif of Nucleolin in G-Quadruplex Stabilization.

G-quadruplexes have important biologic functions that are regulated by G-quadruplex-binding proteins. In particular, G-quadruplex structures are folded or unfolded by their binding proteins and affect transcription and other biologic functions. Here, we investigated the effect of the RNA recognition motif (RRM) and arginine–glycine–glycine repeat (RGG) domain of nucleolin on G-quadruplex formation. Our findings indicate that Phe in the RGG domain of nucleolin is responsible for G-quadruplex binding and folding. Moreover, the RRM of nucleolin potentially binds to a guanine-rich single strand and folds the G-quadruplex with a 5′-terminal and 3′-terminal single strand containing guanine. Our findings contribute to our understanding of how the RRM and RGG domains contribute to G-quadruplex folding and unfolding.

The C-Terminal RRM/ACT Domain Is Crucial for Fine-Tuning the Activation of 'Long' RelA-SpoT Homolog Enzymes by Ribosomal Complexes.

The (p)ppGpp-mediated stringent response is a bacterial stress response implicated in virulence and antibiotic tolerance. Both synthesis and degradation of the (p)ppGpp alarmone nucleotide are mediated by RelA-SpoT Homolog (RSH) enzymes which can be broadly divided in two classes: single-domain ‘short’ and multi-domain ‘long’ RSH. The regulatory ACT (Aspartokinase, Chorismate mutase and TyrA)/RRM (RNA Recognition Motif) domain is a near-universal C-terminal domain of long RSHs. Deletion of RRM in both monofunctional (synthesis-only) RelA as well as bifunctional (i.e., capable of both degrading and synthesizing the alarmone) Rel renders the long RSH cytotoxic due to overproduction of (p)ppGpp. To probe the molecular mechanism underlying this effect we characterized Escherichia coli RelA and Bacillus subtilis Rel RSHs lacking RRM. We demonstrate that, first, the cytotoxicity caused by the removal of RRM is counteracted by secondary mutations that disrupt the interaction of the RSH with the starved ribosomal complex – the ultimate inducer of (p)ppGpp production by RelA and Rel – and, second, that the hydrolytic activity of Rel is not abrogated in the truncated mutant. Therefore, we conclude that the overproduction of (p)ppGpp by RSHs lacking the RRM domain is not explained by a lack of auto-inhibition in the absence of RRM or/and a defect in (p)ppGpp hydrolysis. Instead, we argue that it is driven by misregulation of the RSH activation by the ribosome.

Exploring the proteome associated with the mRNA encoding the D1 reaction center protein of Photosystem II in plant chloroplasts

Synthesis of the D1 reaction center protein of Photosystem II is dynamically regulated in response to environmental and developmental cues. In chloroplasts, much of this regulation occurs at the post-transcriptional level, but the proteins responsible are largely unknown. To discover proteins that impact psbA expression, we identified proteins that associate with maize psbA mRNA by: (i) formaldehyde cross-linking of leaf tissue followed by antisense oligonucleotide affinity capture of psbA mRNA; and (ii) co-immunoprecipitation with HCF173, a psbA translational activator that is known to bind psbA mRNA. The S1 domain protein SRRP1 and two RNA Recognition Motif (RRM) domain proteins, CP33C and CP33B, were enriched with both approaches. Orthologous proteins were also among the enriched protein set in a previous study in Arabidopsis that employed a designer RNA-binding protein as a psbA RNA affinity tag. We show here that CP33B is bound to psbA mRNA in vivo, as was shown previously for CP33C and SRRP1. Immunoblot, pulse labeling, and ribosome profiling analyses of mutants lacking CP33B and/or CP33C detected some decreases in D1 protein levels under some conditions, but no change in psbA RNA abundance or translation. However, analogous experiments showed that SRRP1 represses psbA ribosome association in the dark, represses ycf1 ribosome association, and promotes accumulation of ndhC mRNA. As SRRP1 is known to harbor RNA chaperone activity, we postulate that SRRP1 mediates these effects by modulating RNA structures. The uncharacterized proteins that emerged from our analyses provide a resource for the discovery of proteins that impact the expression of psbA and other chloroplast genes.

Early Metastable Assembly during the Stress-Induced Formation of Worm-like Amyloid Fibrils of Nucleic Acid Binding Domains of TDP-43.

TDP-43 protein travels between the cytosol and the nucleus to perform its nucleic acid binding functions through its two tandem RNA recognition motif domains (TDP-43tRRM). When exposed to various environmental stresses, it forms abnormal aggregates in the cytosol of neurons, which are the hallmarks of amyotrophic lateral sclerosis and other TDP-43 proteinopathies. However, the nature of early structural changes upon stress sensing and the consequent steps during the course of aggregation are not well understood. In this study, we show that under low-pH conditions, mimicking starvation stress, TDP-43tRRM undergoes a conformational opening reaction linked to the protonation of buried ionizable residues and grows into a metastable oligomeric assembly (called the "low-pH form" or the "L form"). In the L form, the protein molecules have disrupted tertiary structure, solvent-exposed hydrophobic patches, and mobile side chains but the native-like secondary structure remains intact. The L form structure is held by weak interactions and has a steep dependence on ionic strength. In the presence of as little as 15 mM KCl, it fully misfolds and further oligomerizes to form a β-sheet rich "β form" in at least two distinct steps. The β form has an ordered, stable structure that resembles worm-like amyloid fibrils. The unstructured regions of the protein gain structure during L ⇌ β conversion. Our results suggest that TDP-43tRRM could function as a stress sensor and support a recent model in which stress sensing during neurodegeneration occurs by assembly of proteins into metastable assemblies that are precursors to the solid aggregates.

A Plant SMALL RNA-BINDING PROTEIN 1 Family Mediates Cell-to-Cell Trafficking of RNAi Signals

In plants, RNA interference (RNAi) plays a pivotal role in growth and development, and responses to environmental inputs, including pathogen attack. The intercellular and systemic trafficking of small interfering RNA (siRNA)/microRNA (miRNA) is a central component in this regulatory pathway. Currently, little is known with regards to the molecular agents involved in the movement of these si/miRNAs. To address this situation, we employed a biochemical approach to identify and characterize a conserved SMALL RNA-BINDING PROTEIN 1 (SRBP1) family that mediates non-cell-autonomous small RNA (sRNA) trafficking. In Arabidopsis, AtSRBP1 is a glycine-rich (GR) RNA-binding protein, also known as AtGRP7, which we show binds single-stranded siRNA. A viral vector, Zucchini yellow mosaic virus (ZYMV), was employed to functionally characterized the AtSRBP1-4 (AtGRP7/2/4/8) RNA recognition motif and GR domains. Cellular-based studies revealed the GR domain as being necessary and sufficient for SRBP1 cell-to-cell movement. Taken together, our findings provide a foundation for future research into the mechanism and function of mobile sRNA signaling agents in plants.

ATP is a cryptic binder of TDP-43 RRM domains to enhance stability and inhibit ALS/AD-associated fibrillation

ATP is the universal energy currency for all cells but has cellular concentrations of 2–12 mM, much higher than required for its classic functions. RNA-recognition motif (RRM) constitutes one of the most abundant domains in eukaryotes and most heterogeneous nuclear ribonucleoproteins (hnRNP) contain RRM domains which not only mediate direct interactions with nucleic acids, but whose aggregation/fibrillation is the pathological hallmark of various human diseases. Here, by NMR and molecular docking, ATP has been decoded to bind TDP-43 two tandem RRM domains with distinctive types of interactions, thus resulting in diverse affinities. Most strikingly, the binding of ATP enhances thermodynamic stability of TDP-43 RRM domains and inhibits ALS-/AD-associated fibrillation. Together, ATP is a cryptic binder of RRM-containing proteins which generally safeguards functional phase separation from transforming into pathological aggregation/fibrillation associated with various diseases and ageing. Our study thus reveals a mechanism of ATP to control protein homeostasis by specific binding.

RBM47-regulated alternative splicing of TJP1 promotes actin stress fiber assembly during epithelial-to-mesenchymal transition.

Morphological and functional changes in cells during the epithelial-mesenchymal transition (EMT) process are known to be regulated by alternative splicing. However, only a few splicing factors involved in EMT have been reported and their underlying mechanisms remain largely unknown. Here, we showed that an isoform of tight junction protein 1 (TJP1) lacking exon 20 (TJP1-α-) is predominantly expressed in tumor tissues and in A549 cells during transforming growth factor-β (TGF-β)-induced EMT. RBM47 promoted the inclusion of exon 20 of TJP1, the alternative exon encoding the α-domain, by which RBM47 recognizes to (U)GCAUG in the downstream intronic region of exon 20. We also found that the first RNA recognition motif (RRM) domain of RBM47 is critical in the regulation of alternative splicing and its recognition to pre-mRNA of TJP1. Furthermore, we demonstrated that the TJP1-α- isoform enhances the assembly of actin stress fibers, thereby promoting cellular migration in a wound healing assay. Our results suggest the regulatory mechanism for the alternative splicing of TJP1 pre-mRNA by RBM47 during EMT, providing a basis for studies related to the modulation of EMT via alternative splicing.

Distinct RNA-binding modules in a single PUF protein cooperate to determine RNA specificity.

PUF proteins, named for Drosophila Pumilio (PUM) and Caenorhabditis elegans fem-3-binding factor (FBF), recognize specific sequences in the mRNAs they bind and control. RNA binding by classical PUF proteins is mediated by a characteristic PUM homology domain (PUM-HD). The Puf1 and Puf2 proteins possess a distinct architecture and comprise a highly conserved subfamily among fungal species. Puf1/Puf2 proteins contain two types of RNA-binding domain: a divergent PUM-HD and an RNA recognition motif (RRM). They recognize RNAs containing UAAU motifs, often in clusters. Here, we report a crystal structure of the PUM-HD of a fungal Puf1 in complex with a dual UAAU motif RNA. Each of the two UAAU tetranucleotides are bound by a Puf1 PUM-HD forming a 2:1 protein-to-RNA complex. We also determined crystal structures of the Puf1 RRM domain that identified a dimerization interface. The PUM-HD and RRM domains act in concert to determine RNA-binding specificity: the PUM-HD dictates binding to UAAU, and dimerization of the RRM domain favors binding to dual UAAU motifs rather than a single UAAU. Cooperative action of the RRM and PUM-HD identifies a new mechanism by which multiple RNA-binding modules in a single protein collaborate to create a unique RNA-binding specificity.

Small Molecule Targeting TDP-43's RNA Recognition Motifs Reduces Locomotor Defects in a Drosophila Model of Amyotrophic Lateral Sclerosis (ALS).

RNA dysregulation likely contributes to disease pathogenesis of amyotrophic lateral sclerosis (ALS) and other neurodegenerative diseases. A pathological form of the transactive response (TAR) DNA binding protein (TDP-43) binds to RNA in stress granules and forms membraneless, amyloid-like TDP-43 aggregates in the cytoplasm of ALS motor neurons. In this study, we hypothesized that by targeting the RNA recognition motif (RRM) domains of TDP-43 that confer a pathogenic interaction between TDP-43 and RNA, motor neuron toxicity could be reduced. In silico docking of 50000 compounds to the RRM domains of TDP-43 identified a small molecule (rTRD01) that (i) bound to TDP-43's RRM1 and RRM2 domains, (ii) partially disrupted TDP-43's interaction with the hexanucleotide RNA repeat of the disease-linked c9orf72 gene, but not with (UG)6 canonical binding sequence of TDP-43, and (iii) improved larval turning, an assay measuring neuromuscular coordination and strength, in an ALS fly model based on the overexpression of mutant TDP-43. Our findings provide an instructive example of a chemical biology approach pivoted to discover small molecules targeting RNA-protein interactions in neurodegenerative diseases.

Unveiling the partners of the DRBD2-mRNP complex, an RBP in Trypanosoma cruzi and ortholog to the yeast SR-protein Gbp2

BackgroundRNA-binding proteins (RBPs) are well known as key factors in gene expression regulation in eukaryotes. These proteins associate with mRNAs and other proteins to form mRNP complexes that ultimately determine the fate of target transcripts in the cell. This association is usually mediated by an RNA-recognition motif (RRM). In the case of trypanosomatids, these proteins play a paramount role, as gene expression regulation is mostly posttranscriptional. Despite their relevance in the life cycle of Trypanosoma cruzi, the causative agent of Chagas’ disease, to date, few RBPs have been characterized in this parasite.ResultsWe investigated the role of DRBD2 in T. cruzi, an RBP with two RRM domains that is associated with cytoplasmic translational complexes. We show that DRBD2 is an ortholog of the Gbp2 in yeast, an SR-rich protein involved in mRNA quality control and export. We used an immunoprecipitation assay followed by shotgun proteomics and RNA-seq to assess the interaction partners of the DRBD2-mRNP complex in epimastigotes. The analysis identified mostly proteins involved in RNA metabolism and regulation, such as ALBA1, ALBA3, ALBA4, UBP1, UBP2, DRBD3, and PABP2. The RNA-seq results showed that most of the transcripts regulated by the DRBD2 complex mapped to hypothetical proteins related to multiple processes, such as to biosynthetic process, DNA metabolic process, protein modification, and response to stress.ConclusionsThe identification of regulatory proteins in the DRBD2-mRNP complex corroborates the important role of DRBD2 in gene expression regulation in T. cruzi. We consider these results an important contribution to future studies regarding gene expression regulation in T. cruzi, especially in the field of RNA-binding proteins.

A telomerase subunit homolog La protein from Trypanosomabrucei plays an essential role in ribosomal biogenesis.

The autoantigen La protein is an important component of telomerase and a predominantly nuclear phosphoprotein. As a telomerase subunit, La protein associates with the telomerase ribonucleoprotein and influences telomere length. In the reverse transcription, La protein stimulates enzymatic activity and increases repeated addition processivity of telomerase. As nuclear phosphoprotein, La protein binds the 3' poly(U)-rich elements of nascent RNA polymerase III transcripts to facilitate its correct folding and maturation. In this work, we identified a La protein homolog (TbLa) from Trypanosomabrucei (T.brucei). We revealed that TbLa interacts with ribosome-associated protein P34/P37, 40S ribosomal protein SA, and 60S ribosomal subunit L5 in T.brucei. In the interactions between TbLa protein and (P34/P37)/L5/SA, RNA recognition motif (RRM) domain of TbLa was indicated to make the major contribution to the processes. We determined the solution structure of TbLa RRM domain. NMR chemical shift perturbations revealed that the positively charged RNA-binding pocket of TbLa RRM domain is responsible for its interaction with ribosomal and ribosome-associated proteins P37/L5/SA. Furthermore, depletion of TbLa affected the maturation process of 5S rRNA and ribosomal assembly, suggesting TbLa protein might play a significant role in the ribosomal biogenesis pathway in T.brucei. Taken together, our results provide a novel insight and structural basis for better understanding the roles of TbLa and RRM domain in ribosomal biogenesis in T.brucei. DATABASE: Structural data are available in the PDB under the accession number 5ZUH.

Origins and evolution of CRISPR-Cas systems.

CRISPR-Cas, the bacterial and archaeal adaptive immunity systems, encompass a complex machinery that integrates fragments of foreign nucleic acids, mostly from mobile genetic elements (MGE), into CRISPR arrays embedded in microbial genomes. Transcripts of the inserted segments (spacers) are employed by CRISPR-Cas systems as guide (g)RNAs for recognition and inactivation of the cognate targets. The CRISPR-Cas systems consist of distinct adaptation and effector modules whose evolutionary trajectories appear to be at least partially independent. Comparative genome analysis reveals the origin of the adaptation module from casposons, a distinct type of transposons, which employ a homologue of Cas1 protein, the integrase responsible for the spacer incorporation into CRISPR arrays, as the transposase. The origin of the effector module(s) is far less clear. The CRISPR-Cas systems are partitioned into two classes, class 1 with multisubunit effectors, and class 2 in which the effector consists of a single, large protein. The class 2 effectors originate from nucleases encoded by different MGE, whereas the origin of the class 1 effector complexes remains murky. However, the recent discovery of a signalling pathway built into the type III systems of class 1 might offer a clue, suggesting that type III effector modules could have evolved from a signal transduction system involved in stress-induced programmed cell death. The subsequent evolution of the class 1 effector complexes through serial gene duplication and displacement, primarily of genes for proteins containing RNA recognition motif domains, can be hypothetically reconstructed. In addition to the multiple contributions of MGE to the evolution of CRISPR-Cas, the reverse flow of information is notable, namely, recruitment of minimalist variants of CRISPR-Cas systems by MGE for functions that remain to be elucidated. Here, we attempt a synthesis of the diverse threads that shed light on CRISPR-Cas origins and evolution.This article is part of a discussion meeting issue ‘The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems’.

Molecular Dynamics Simulations of RNA-Recognition Motif Complexed with CAC-Containing RNA

The RNA-recognition motifs (RRMs) are the most abundant RNA-binding domains in higher vertebrates, which play diverse roles in post-transcriptional gene expression processes. In this work, the free RRM domain and its RNA-binding complex structure are studied by molecular dynamics (MD) simulations. Principal component analysis (PCA) and Molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) methods were used to explore the dynamical and recognition mechanisms of RRM domain and CAC-containing RNA. Through the motion mode comparison between free RRM domain and complex systems, it's found that the RNA-binding interface of free RRM domain is not stable, especially the C-terminal loop. The C-terminal loop in free RRM domain adopts the conformation similar to that in the complex system. The energy decomposition strategy of MM/PBSA was used to analyze the contribution of each residue or nucleotide. In complex system, 20 residues make favorable contributions for RNA binding, which are positioned on the homodimerization interface (such as Arg38, Glu39 and Lys36) or the RNA-binding interface of RRM domain (such as Phe27, Phe65, Lys100, Thr103 and Lys104). These residues are also reported by the previous experiments, which may form hydrogen bonds or intermolecular stacking interactions with RNA. Meanwhile, 11 residues in RRM domain contribute unfavorable energies for RNA binding. These residues are mainly located on the loop structure or the regions far away to the RNA binding interface. In the energy decomposition of RNA, all of 4 nucleotides make favorable energy contributions. The three nucleotides CAC contribute most of the binding energies, so these 3 nucleotides are more stable than U1 in the representative structures. This study provides some new insights into the RNA recognition mechanisms of RRM domain.

A glycine-rich protein MoGrp1 functions as a novel splicing factor to regulate fungal virulence and growth in Magnaporthe oryzae

Glycine-rich proteins (GRPs) have diverse amino acid sequences and are involved in a variety of biological processes. The role of GRPs in plant pathogenic fungi has not been reported. In this study, we identified and functionally characterized a novel gene named MoGRP1 in Magnaporthe oryzae, which encodes a protein that has an N-terminal RNA recognition motif (RRM) and a C-terminal glycine-rich domain with four Arg-Gly-Gly (RGG) repeats. Deletion of MoGRP1 resulted in dramatic reductions in fungal virulence, mycelial growth, and conidiation. The ΔMogrp1 mutants were also defective in cell wall integrity and in their responses to different stresses. MoGrp1 was localized to the nucleus and was co-immunoprecipitated with several components of the spliceosome, including subunits of the U1 snRNP and U2 snRNP complexes. Moreover, MoGrp1 exhibited binding affinity for poly(U). Importantly, MoGrp1 was responsible for the normal splicing of genes involved in infection-related morphogenesis. Domain deletion assays showed that both the RRM domain and its two adjacent RGG repeats were essential to the full function of MoGrp1. Notably, the nine amino acids between the first and the second RGG repeats were indispensable for nuclear localization and for the biological functions of MoGrp1. Taken together, our data suggest that MoGrp1 functions as a novel splicing factor with poly(U) binding activity to regulate fungal virulence, development, and stress responses in the rice blast fungus.

Genome-Wide Identification of Cyclophilin Gene Family in Cotton and Expression Analysis of the Fibre Development in Gossypium barbadense.

Cyclophilins (CYPs) are a member of the immunophilin superfamily (in addition to FKBPs and parvulins) and play a significant role in peptidyl-prolyl cis-trans isomerase (PPIase) activity. Previous studies have shown that CYPs have important functions in plants, but no genome-wide analysis of the cotton CYP gene family has been reported, and the specific biological function of this gene is still elusive. Based on the release of the cotton genome sequence, we identified 75, 78, 40 and 38 CYP gene sequences from G. barbadense, G. hirsutum, G. arboreum, and G. raimondii, respectively; 221 CYP genes were unequally located on chromosomes. Phylogenetic analysis showed that 231 CYP genes clustered into three major groups and eight subgroups. Collinearity analysis showed that segmental duplications played a significant role in the expansion of CYP members in cotton. There were light-responsiveness, abiotic-stress and hormone-response elements upstream of most of the CYPs. In addition, the motif composition analysis revealed that 49 cyclophilin proteins had extra domains, including TPR (tetratricopeptide repeat), coiled coil, U-box, RRM (RNA recognition motif), WD40 (RNA recognition motif) and zinc finger domains, along with the cyclophilin-like domain (CLD). The expression patterns based on qRT-PCR showed that six CYP expression levels showed greater differences between Xinhai21 (long fibres, G. barbadense) and Ashmon (short fibres, G. barbadense) at 10 and 20 days postanthesis (DPA). These results signified that CYP genes are involved in the elongation stage of cotton fibre development. This study provides a valuable resource for further investigations of CYP gene functions and molecular mechanisms in cotton.

RRMdb-an evolutionary-oriented database of RNA recognition motif sequences.

RNA-recognition motif (RRM) is an RNA-interacting protein domain that plays an important role in the processes of RNA metabolism such as the splicing, editing, export, degradation, and regulation of translation. Here, we present the RNA-recognition motif database (RRMdb), which affords rapid identification and annotation of RRM domains in a given protein sequence. The RRMdb database is compiled from ~57 000 collected representative RRM domain sequences, classified into 415 families. Whenever possible, the families are associated with the available literature and structural data. Moreover, the RRM families are organized into a network of sequence similarities that allows for the assessment of the evolutionary relationships between them.