RNA Recognition Motif Domain Research Articles

Pseudouridine-modified tRNA fragments repress aberrant protein synthesis and predict leukaemic progression in myelodysplastic syndrome

Transfer RNA-derived fragments (tRFs) are emerging small noncoding RNAs that, although commonly altered in cancer, have poorly defined roles in tumorigenesis1. Here we show that pseudouridylation (Ψ) of a stem cell-enriched tRF subtype2, mini tRFs containing a 5′ terminal oligoguanine (mTOG), selectively inhibits aberrant protein synthesis programmes, thereby promoting engraftment and differentiation of haematopoietic stem and progenitor cells (HSPCs) in patients with myelodysplastic syndrome (MDS). Building on evidence that mTOG-Ψ targets polyadenylate-binding protein cytoplasmic 1 (PABPC1), we employed isotope exchange proteomics to reveal critical interactions between mTOG and functional RNA-recognition motif (RRM) domains of PABPC1. Mechanistically, this hinders the recruitment of translational co-activator PABPC1-interacting protein 1 (PAIP1)3 and strongly represses the translation of transcripts sharing pyrimidine-enriched sequences (PES) at the 5′ untranslated region (UTR), including 5′ terminal oligopyrimidine tracts (TOP) that encode protein machinery components and are frequently altered in cancer4. Significantly, mTOG dysregulation leads to aberrantly increased translation of 5′ PES messenger RNA (mRNA) in malignant MDS-HSPCs and is clinically associated with leukaemic transformation and reduced patient survival. These findings define a critical role for tRFs and Ψ in difficult-to-treat subsets of MDS characterized by high risk of progression to acute myeloid leukaemia (AML).

Molecular characterization and expression profiling oftransformer 2 and fruitless-like homologs in the black tiger shrimp, Penaeus monodon.

Transformer 2 (tra 2) and fruitless (fru) genes have been proven to play a key role in sex determination pathways in many Arthropods, including insects and crustaceans. In this study, a paralog of P. monodon tra 2 (Pmtra 2), P. monodon ovarian associated transformer 2 (PmOvtra 2) and 2 isoforms of P. monodon fruitless-like gene (Pmfru-1 and Pmfru-2) were identified and characterized. The full cDNA sequence of PmOvtra 2 consisted of 1,774 bp with the longest open reading frame (ORF) of 744 bp encoding for 247 amino acids. The PmOvtra 2 exhibited a predicted RNA-recognition motif (RRM) domain and two arginine-serine (RS) regions, suggesting its function in RNA splicing. The full cDNA sequence of Pmfru-1 consisted of 1,306 bp with 1,182 bp ORF encoding for 393 amino acids, whereas the full cDNA sequence of Pmfru-2 consisted of 1,858 bp with 1,437 bp ORF encoding 478 amino acids. The deduced amino acid sequences of Pmfru-1 and Pmfru-2 exhibited highly conserved domains of Fru proteins, including Broad-complex, Tramtrack and Bric-a-brac (BTB), and zinc finger (ZF) domains. In addition, Pmfru-1 and Pmfru-2 were suggestively originated from the same single genomic locus by genomic sequence analysis. Specifically, Pmfru pre-mRNA was alternatively spliced for Pmfru-1 and Pmfru-2 to include mutually exclusive exon 7 and exon 6, respectively. Temporal and spatial expression of PmOvtra 2, Pmfru-1, and Pmfru-2 were also investigated by qPCR. The results showed that all were expressed in early developmental stages with undifferentiated gonads starting from nauplius until postlarvae. The expression of PmOvtra 2 started at nauplius stage and gradually increased from mysis to postlarvae (PL) 1. However, the expression of Pmfru-1 was low at the nauplii stage and slightly increased from protozoea to PL5, whereas the expression of Pmfru-2 maintained a low level from nauplius to mysis and then gradually increased at the PL stages. Expressions of PmOvtra 2, Pmfru-1, and Pmfru-2 were detected in various tissues including nervous tissue, gill, heart, hepatopancreas, gut, and gonads. Interestingly, the sexually dimorphic expression of PmOvtra 2, Pmfru-1, and Pmfru-2 was demonstrated in fully developed gonads in which the ovary showed significantly higher expressions than the testis. The great difference in the expression pattern of PmOvtra 2, Pmfru-1, and Pmfru-2 in the ovary and testis suggested their roles in the female sex determination in P. monodon.

An interactome analysis reveals that Arabidopsis thaliana GRDP2 interacts with proteins involved in post-transcriptional processes.

The Arabidopsis thaliana glycine-rich domain protein 2 (AtGRDP2) gene encodes a protein of unknown function that is involved in plant growth and salt stress tolerance. The AtGRDP2 protein (787 aa, At4g37900) is constituted by three domains: a DUF1399 located at the N-terminus, a potential RNA Recognition Motif (RRM) in the central region, and a short glycine-rich domain at the C-terminus. Herein, we analyzed the subcellular localization of AtGRDP2 protein as a GFP translational fusion and found it was localized in the cytosol and the nucleus of tobacco leaf cells. Truncated versions of AtGRDP2 showed that the DUF1399 or the RRM domains were sufficient for nuclear localization. In addition, we performed a yeast two-hybrid split-ubiquitin assay (Y2H) to identify potential interactors for AtGRDP2 protein. The Y2H assay identified proteins associated with RNA binding functions such as PABN3 (At5g65260), EF-1α (At1g07920), and CL15 (At3g25920). Heterodimeric associations in planta between AtGRDP2 and its interactors were carried out by Bimolecular Fluorescence Complementation (BiFC) assays. The data revealed heterodimeric interactions between AtGRDP2 and PABN3 in the nucleus and AtGRDP2 with EF-1α in the cytosol, while AtGRDP2-CL15 associations occurred only in the chloroplasts. Finally, functional characterization of the protein-protein interaction regions revealed that both DUF1399 and RRM domains were key for heterodimerization with its interactors. The AtGRDP2 interaction with these proteins in different compartments suggests that this glycine-rich domain protein is involved in post-transcriptional processes.

1H, 13C and 15N resonance assignments of stress granule key component G3BP1 RRM domain.

G3BP1 exists as a helicase and one of the core components in stress granules, which are associated with a variety of neurodegenerative diseases. Its RNA recognition motif (RRM) domain performs the paramount function of binding mRNA. Here we report the resonance assignment of human G3BP1 RRM domain to understand its structure–function relationship.

Characterization of Sequence-Specific Binding of LARP6 to the 5′ Stem-Loop of Type I Collagen mRNAs and Implications for Rational Design of Antifibrotic Drugs

Excessive synthesis of type I collagen is a hallmark of fibrotic diseases. Binding of La-related protein 6 (LARP6) to the 5′ stem-loop (5′SL) of collagen mRNAs regulates their translation leading to an unnaturally elevated rate of collagen biosynthesis in fibrosis. Previous work suggested that LARP6 needs two domains to form stable complex with 5′SL RNA, the La domain and the juxtaposed RNA recognition motif (RRM), jointly called the La-module. Here we describe that La domain of LARP6 is necessary and sufficient for recognition of 5′SL in RNA sequence specific manner. A three-amino-acid motif located in the flexible loop connecting the second α-helix to the β-sheet of the La domain, called the RNK-motif, is critical for binding. Mutation of any of these three amino acids abolishes the binding of the La domain to 5′SL. The major site of crosslinking of LARP6 to 5′SL RNA was mapped to this motif, as well. The RNK-motif is not found in other LARPs, which cannot bind 5′SL. Presence of RRM increases the stability of complex between La domain and 5′SL RNA and RRM domain does not make extensive contacts with 5′SL RNA. We propose a model in which the initial recognition of 5′SL by LARP6 is mediated by the RNK epitope and further stabilized by the RRM domain. This discovery suggests that the interaction between LARP6 and collagen mRNAs can be blocked by small molecules that target the RNK epitope and will help rational design of the LARP6 binding inhibitors as specific antifibrotic drugs.

1H, 13C and 15N resonance assignments and solution structures of the two RRM domains of Matrin-3.

Matrin-3 is a multifunctional protein that binds to both DNA and RNA. Its DNA-binding activity is linked to the formation of the nuclear matrix and transcriptional regulation, while its RNA-binding activity is linked to mRNA metabolism including splicing, transport, stabilization, and degradation. Correspondingly, Matrin-3 has two zinc finger domains for DNA binding and two consecutive RNA recognition motif (RRM) domains for RNA binding. Matrin-3 has been reported to cause amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) when its disordered region contains pathogenic mutations. Simultaneously, it has been shown that the RNA-binding activity of Matrin-3 mediated by its RRM domains, affects the formation of insoluble cytoplasmic granules, which are related to the pathogenic mechanism of ALS/FTD. Thus, the effect of the RRM domains on the phase separation of condensed protein/RNA mixtures has to be clarified for a comprehensive understanding of ALS/FTD. Here, we report the 1H, 15N, and 13C resonance assignments of the two RNA binding domains and their solution structures. The resonance assignments and the solution structures obtained in this work will contribute to the elucidation of the molecular basis of Matrin-3 in the pathogenic mechanism of ALS and/or FTD.

Structural and biochemical insights into CRISPR RNA processing by the Cas5c ribonuclease SMU1763 from Streptococcus mutans

The cariogenic pathogen Streptococcus mutans contains two CRISPR systems (type I-C and type II-A) with the Cas5c protein (SmuCas5c) involved in processing of long CRISPR RNA transcripts (pre-crRNA) containing repeats and spacers to mature crRNA guides. In this study, we determined the crystal structure of SmuCas5c at a resolution of 1.72 Å, which revealed the presence of an N-terminal modified RNA recognition motif and a C-terminal twisted β-sheet domain with four bound sulphate molecules. Analysis of surface charge and residue conservation of the SmuCas5c structure suggested the location of an RNA-binding site in a shallow groove formed by the RNA recognition motif domain with several conserved positively charged residues (Arg39, Lys52, Arg109, Arg127, and Arg134). Purified SmuCas5c exhibited metal-independent ribonuclease activity against single-stranded pre-CRISPR RNAs containing a stem–loop structure with a seven-nucleotide stem and a pentaloop. We found SmuCas5c cleaves substrate RNA within the repeat sequence at a single cleavage site located at the 3′-base of the stem but shows significant tolerance to substrate sequence variations downstream of the cleavage site. Structure-based mutational analysis revealed that the conserved residues Tyr50, Lys120, and His121 comprise the SmuCas5c catalytic residues. In addition, site-directed mutagenesis of positively charged residues Lys52, Arg109, and Arg134 located near the catalytic triad had strong negative effects on the RNase activity of this protein, suggesting that these residues are involved in RNA binding. Taken together, our results reveal functional diversity of Cas5c ribonucleases and provide further insight into the molecular mechanisms of substrate selectivity and activity of these enzymes.

RNA binding motif 47 (RBM47): emerging roles in vertebrate development, RNA editing and cancer.

RNA-binding proteins (RBPs) are critical players in the post-transcriptional regulation of gene expression and are associated with each event in RNA metabolism. The term 'RNA-binding motif' (RBM) is assigned to novel RBPs with one or more RNA recognition motif (RRM) domains that are mainly involved in the nuclear processing of RNAs. RBM47 is a novel RBP conserved in vertebrates with three RRM domains whose contributions to various aspects of cellular functions are as yet emerging. Loss of RBM47 function affects head morphogenesis in zebrafish embryos and leads to perinatal lethality in mouse embryos, thereby assigning it to be an essential gene in early development of vertebrates. Its function as an essential cofactor for APOBEC1 in C to U RNA editing of several targets through substitution for A1CF in the A1CF-APOBEC1 editosome, established a new paradigm in the field. Recent advances in the understanding of its involvement in cancer progression assigned RBM47 to be a tumor suppressor that acts by inhibiting EMT and Wnt/[Formula: see text]-catenin signaling through post-transcriptional regulation. RBM47 is also required to maintain immune homeostasis, which adds another facet to its regulatory role in cellular functions. Here, we review the emerging roles of RBM47 in various biological contexts and discuss the current gaps in our knowledge alongside future perspectives for the field.

Insights into the aggregation mechanism of RNA recognition motif domains in TDP-43: a theoretical exploration.

The transactive response DNA-binding protein 43 (TDP-43) is associated with several diseases such as amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD) due to pathogenic aggregations. In this work, we examined the dimer, tetramer and hexamer models built from the RRM domains of TDP-43 using molecular dynamics simulations in combination with the protein–protein docking. Our results showed that the formations of the dimer models are mainly achieved by the interactions of the RRM1 domains. The parallel β-sheet layers between the RRM1 domains provide most of the binding sites in these oligomer models, and thus play an important role in the aggregation process. The approaching of the parallel β-sheet layers from small oligomer models gradually expand to large ones through the allosteric communication between the α1/α2 helices of the RRM1 domains, which maintains the binding affinities and interactions in the larger oligomer models. Using the repeatable-superimposing method based on the tetramer models, we proposed a new aggregation mechanism of RRM domains in TDP-43, which could well characterize the formation of the large aggregation models with the repeated, helical and rope-like structures. These new insights help to understand the amyloid-like aggregation phenomena of TDP-43 protein in ALS and FTLD diseases.

Further evidence for de novo variants in SYNCRIP as the cause of a neurodevelopmental disorder.

SYNCRIP encodes for the Synaptotagmin-binding cytoplasmic RNA-interacting protein, involved in RNA-binding and regulation of multiple cellular pathways. It has been proposed as a candidate gene for neurodevelopmental disorders (NDDs) with autism spectrum disorder (ASD), intellectual disability (ID), and epilepsy. We ascertained genetic, clinical, and neuroradiological data of three additional individuals with novel de novo SYNCRIP variants. All individuals had ID. Autistic features were observed in two. One individual showed myoclonic-atonic epilepsy. Neuroradiological features comprised periventricular nodular heterotopia and widening of subarachnoid spaces. Two frameshift variants in the more severely affected individuals, likely result in haploinsufficiency. The third missense variant lies in the conserved RNA recognition motif (RRM) 2 domain likely affecting RNA-binding. Our findings support the importance of RRM domains for SYNCRIP functionality and suggest genotype-phenotype correlations. Our study provides further evidence for a SYNCRIP-associated NDD characterized by ID and ASD sporadically accompanied by malformations of cortical development and myoclonic-atonic epilepsy.

A p53/CPEB2 negative feedback loop regulates renal cancer cell proliferation and migration

The tumor suppressor p53 transactivates the expression of multiple genes to exert its multifaceted functions and ultimately maintains genome stability. Thus, cancer cells develop various mechanisms to diminish p53 expression and bypass the cell cycle checkpoint. In this study, we identified the gene encoding RNA-binding protein cytoplasmic polyadenylation element-binding protein 2 (CPEB2) as a p53 target. In turn, CPEB2 decreases p53 messenger RNA stability and translation to fine-tune p53 level. Specifically, we showed that CPEB2 binds the cytoplasmic polyadenylation elements in the p53 3′-untranslated region, and the RNA recognition motif and zinc finger (ZF) domains of CPEB2 are required for this binding. Furthermore, we found that CPEB2 was upregulated in renal cancer tissues and promotes the renal cancer cell proliferation and migration. The oncogenic effect of CPEB2 is partially dependent on negative feedback regulation of p53. Overall, we identify a novel regulatory feedback loop between p53 and CPEB2 and demonstrate that CPEB2 promotes tumor progression by inactivating p53, suggesting that CPEB2 is a potential therapeutic target in human renal cancer.

Human DND1-RRM2 forms a non-canonical domain swapped dimer.

RNA recognition motif (RRM) being the most abundant RNA binding domain in eukaryotes, is a major player in cellular regulation. Several variations in the canonical βαββαβ topology have been observed. We have determined the 2.3Å crystal structure of the human DND1-RRM2 domain. The structure revealed an interesting non-canonical RRM fold, which is maintained by the formation of a 3D domain swapped dimer between β1 and β4 strands across protomers. We have delineated the structural basis of the stable domain swapped dimer formation using the residue level dynamics of protein explored by NMR spectroscopy and MD simulations. Our structural and dynamics studies substantiate major determinants and molecular basis for domain swapped dimerization observed in the RRM domain.

Molecular Cloning, Expression and Characterization of Peroxisome Proliferators-Activated Receptors Gamma in the Sea Urchin (Strongylocentrotus intermedius)

Peroxisome proliferators-activated receptor gamma (PPARγ) plays important regulatory roles in adipocyte differentiation. In this study, we cloned the full-length sequence of the PPARγ gene and analyzed its expression profile in different developmental stages and tissues of Strongylocentrotus intermedius. The full-length cDNA of PPARγ contains 2286 base pairs (bp) with a putative open reading frame of 1755 bp, and the gene encodes encoding a polypeptide of 584 amino acid residues. The predicted molecular mass of the protein is 67.27 kDa, and its theoretical isoelectric point (pI) is 10.07. The protein contains conserved motifs, including an RRM (RNA recognition motif) domain. PPARγ expression with the highest level was observed in the gonad, and the lowest was observed in the tube feet of S. intermedius. Time-course expression measurements at different developmental stages showed that the highest expression level of PPARγ is in the eggs and its weakest expression level is in the 32-cells stage. Knock-down of PPARγ by specific siRNA revealed that UCP2 expression is significantly decreased in the gonads and intestines 48 h post-transfection, indicating that the UCP2 is a downstream target gene of PPARγ. This finding suggests that PPARγ and UCP2 have positive regulatory effects in regulating adipocyte differentiation. Changes in fatty acid levels in the gonads before and after PPARγ interference were assessed, and decreased C18:2(trans, n™6) and C20:3(n−6) levels were observed 48 h after siRNA transfection. The results showed the function of PPARγ in fatty acid anabolism, The data are helpful to improve the current understanding of the fatty acid synthesis pathways and regulatory mechanisms in S. intermedius. They also provide an experimental basis for improving fatty acid synthesis in sea urchins, which is important for cultivating sea urchins with high nutritional value.

Cyp33 binds AU-rich RNA motifs via an extended interface that competitively disrupts the gene repressive Cyp33-MLL1 interaction in vitro.

Cyp33 is an essential human cyclophilin prolyl isomerase that plays myriad roles in splicing and chromatin remodeling. In addition to a canonical cyclophilin (Cyp) domain, Cyp33 contains an RNA-recognition motif (RRM) domain, and RNA-binding triggers proline isomerase activity. One prominent role for Cyp33 is through a direct interaction with the mixed lineage leukemia protein 1 (MLL1, also known as KMT2A) complex, which is a histone methyltransferase that serves as a global regulator of human transcription. MLL activity is regulated by Cyp33, which isomerizes a key proline in the linker between the PHD3 and Bromo domains of MLL1, acting as a switch between gene activation and repression. The direct interaction between MLL1 and Cyp33 is critical, as deletion of the MLL1-PHD3 domain responsible for this interaction results in oncogenesis. The Cyp33 RRM is central to these activities, as it binds both the PHD3 domain and RNA. To better understand how RNA binding drives the action of Cyp33, we performed RNA-SELEX against full-length Cyp33 accompanied by deep sequencing. We have identified an enriched Cyp33 binding motif (AAUAAUAA) broadly represented in the cellular RNA pool as well as tightly binding RNA aptamers with affinities comparable and competitive with the Cyp33 MLL1-PHD3 interaction. RNA binding extends beyond the canonical RRM domain, but not to the Cyp domain, suggesting an indirect mechanism of interaction. NMR chemical shift mapping confirms an overlapping, but not identical, interface on Cyp33 for RNA and PHD3 binding. This finding suggests RNA can disrupt the gene repressive Cyp33-MLL1 complex providing another layer of regulation for chromatin remodeling by MLL1.

A combined NMR and EPR investigation on the effect of the disordered RGG regions in the structure and the activity of the RRM domain of FUS

Structural disorder represents a key feature in the mechanism of action of RNA-binding proteins (RBPs). Recent insights revealed that intrinsically disordered regions (IDRs) linking globular domains modulate their capability to interact with various sequences of RNA, but also regulate aggregation processes, stress-granules formation, and binding to other proteins. The FET protein family, which includes FUS (Fused in Sarcoma), EWG (Ewing Sarcoma) and TAF15 (TATA binding association factor 15) proteins, is a group of RBPs containing three different long IDRs characterized by the presence of RGG motifs. In this study, we present the characterization of a fragment of FUS comprising two RGG regions flanking the RNA Recognition Motif (RRM) alone and in the presence of a stem-loop RNA. From a combination of EPR and NMR spectroscopies, we established that the two RGG regions transiently interact with the RRM itself. These interactions may play a role in the recognition of stem-loop RNA, without a disorder-to-order transition but retaining high dynamics.

Resonance assignments of La protein RRM domain from Trypanosoma brucei.

The autoantigen La protein is a conserved component of eukaryotic ribonucleoprotein complexes that binds the 3' poly(U) sequences of nascent RNA polymerase III transcripts to assist folding and maturation. This specific recognition is mediated by the N-terminal domain (NTD) of La, which comprises a La motif and an RNA recognition motif (RRM). Here, we report near complete 1H, 13C and 15N backbone and sidechain assignments for the RRM domain of La protein from Trypanosoma brucei.

VdNop12, containing two tandem RNA recognition motif domains, is a crucial factor for pathogenicity and cold adaption in Verticillium dahliae.

Previous studies have reported the ability of fungi to overwinter in soil or on crop debris under different environmental conditions, but how fungi adapt to chilling is still largely unknown. In this study, we have identified and characterized the RNA binding protein (RBP) (VdNop12) by screening an Agrobacterium tumefaciens-mediated transformation-mediated insertional mutational library of Verticillium dahliae. We determined that this protein was essential to the pathogen for virulence on cotton plants. VdNop12 contains two tandem RNA recognition motif domains, and its orthologs are widely distributed in filamentous fungi. Mutants produced by disruption of VdNop12 showed defects in vegetative growth, conidiation and cell wall integrity. The mutant also showed an increase in sensitivity to low temperature, as compared to the wildtype and complementation strains. Yeast complementation assay showed that VdNop12 could functionally restore the growth phenotype of ΔScNop12 mutant of Saccharomyces cerevisiae at 15°C. We demonstrated that the VdNop12 is localized in the nucleus, and its loss resulted in the downregulated expression of several genes related to cAMP-PKA and MAPK pathways in V. dahliae. Our results demonstrated a crucial role of RBPs in the regulation of morphology, cold adaption, and pathogenic development in V. dahliae.

Identification and expression analysis of chloroplast ribonucleoproteins (cpRNPs) in Arabidopsis and rice.

Chloroplast ribonucleoproteins (cpRNPs) are implicated in splicing, editing, and stability control of chloroplast RNAs as well as in regulating development and stress tolerance. To facilitate a comprehensive understanding of their functions, we carried out a genome-wide identification, curation, and phylogenetic analysis of cpRNP genes in Oryza sativa (rice) and Arabidopsis thaliana (Arabidopsis). Ten cpRNP genes were identified in each of Arabidopsis and rice genomes based on the presence of two RRM (RNA-recognition motif) domains and an N-terminal chloroplast targeting signal peptide in the predicted proteins. These proteins are localized to chloroplasts. Gene expression analysis revealed that cpRNP genes have differential tissue expression patterns and some cpRNP genes are induced by abiotic stresses such as cold, heat, and drought. Taken together, our study provides a comprehensive annotation of the cpRNP gene family and their expression patterns in Arabidopsis and rice which will facilitate further studies on their roles in plant growth and stress responses.

Identification of binding domains and key amino acids involved in the interaction between BmLARK and G4 structure in the BmPOUM2 promoter in Bombyx mori.

It has been found that the non-B form DNA structures, like G-quadruplex (G4) and i-motif, are involved in many important biological processes. Our previous study showed that the silkworm transcription factor BmLARK binds to the G4 structure in the promoter of the transcription factor BmPOUM2 and regulates its promoter activity. However, the binding mechanism between BmLARK and BmPOUM2 G4 structure remains unclear. In this study, binding domains and key amino acid residues involved in the interaction between BmLARK and BmPOUM2 G4 were studied. The electrophoretic mobility shift assay results indicated that the two RNA-recognition motifs (RRM) of BmLARK are simultaneously required for the binding with the G4 structure. Either RRM1 or RRM2 alone could not bind with the G4 structure. The zinc-finger motif was not involved in the binding. A series of mutant proteins with specific amino acid mutations were expressed and used to identify the key amino acid residues involving the interaction. The results indicated that β sheets, especially the β1 and β3 sheets, in the RRM domains of BmLARK played critical roles in the binding with the G4 structure. Several amino acid mutations of RRM1/2 in ribonucleoprotein domain 1 (RNP1) (motif in β3 strand) and RNP2 (motif in β1 strand) caused loss of binding ability, indicating that these amino acids are the key sites for the binding. All the results suggest that RRM domains, particularly their the RNP1 and RNP2 motifs, play important roles not only in RNA recognition, but also in the G4 structure binding.

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.