Genetic determinants of mate recognition in Brachionus manjavacas (Rotifera).

BACKGROUND: Mate choice is of central importance to most animals, influencing population structure, speciation, and ultimately the survival of a species. Mating behavior of male brachionid rotifers is triggered by the product of a chemosensory gene, a glycoprotein on the body surface of females called the mate recognition pheromone. The mate recognition pheromone has been biochemically characterized, but little was known about the gene(s). We describe the isolation and characterization of the mate recognition pheromone gene through protein purification, N-terminal amino acid sequence determination, identification of the mate recognition pheromone gene from a cDNA library, sequencing, and RNAi knockdown to confirm the functional role of the mate recognition pheromone gene in rotifer mating. RESULTS: A 29 kD protein capable of eliciting rotifer male circling was isolated by high-performance liquid chromatography. Two transcript types containing the N-terminal sequence were identified in a cDNA library; further characterization by screening a genomic library and by polymerase chain reaction revealed two genes belonging to each type. Each gene begins with a signal peptide region followed by nearly perfect repeats of an 87 to 92 codon motif with no codons between repeats and the final motif prematurely terminated by the stop codon. The two Type A genes contain four and seven repeats and the two Type B genes contain three and five repeats, respectively. Only the Type B gene with three repeats encodes a peptide with a molecular weight of 29 kD. Each repeat of the Type B gene products contains three asparagines as potential sites for N-glycosylation; there are no asparagines in the Type A genes. RNAi with Type A double-stranded RNA did not result in less circling than in the phosphate-buffered saline control, but transfection with Type B double-stranded RNA significantly reduced male circling by 17%. The very low divergence between repeat units, even at synonymous positions, suggests that the repeats are kept nearly identical through a process of concerted evolution. Information-rich molecules like surface glycoproteins are well adapted for chemical communication and aquatic animals may have evolved signaling systems based on these compounds, whereas insects use cuticular hydrocarbons. CONCLUSION: Owing to its critical role in mating, the mate recognition pheromone gene will be a useful molecular marker for exploring the mechanisms and rates of selection and the evolution of reproductive isolation and speciation using rotifers as a model system. The phylogenetic variation in the mate recognition pheromone gene can now be studied in conjunction with the large amount of ecological and population genetic data being gathered for the Brachionus plicatilis species complex to understand better the evolutionary drivers of cryptic speciation.

Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch.

Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA's secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent (approximately 1 micros total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex.

Enzymatic ligation strategies for the preparation of purine riboswitches with site-specific chemical modifications.

One of the most versatile riboswitch classes refers to purine nucleoside metabolism. In the cell, purine riboswitches of the respective mRNAs either act at the transcriptional or translational level and off- or on-regulate genes upon binding to their dedicated ligands. Biophysical studies on ligand-induced folding of these RNA domains in vitro contribute to understanding their regulation mechanisms in vivo. For such studies, in particular, for approaches using fluorescence spectroscopy, the preparation of large RNAs with site-specific chemical modifications is required. Here, we describe a strategy for the preparation of riboswitch aptamers and aptamers adjoined to their expression platforms by chemical synthesis and enzymatic ligation. The modular design enables fast access to a large number of purine riboswitch derivatives with the modification of interest at any strand position. We exemplarily provide a detailed protocol for the preparation of adenosine deaminase (add) A-riboswitch variants with 2-aminopurine (AP) modifications at the 40-nmol scale.

Application of fluorescent measurements for characterization of riboswitch-ligand interactions.

Riboswitches are recently discovered messenger RNA motifs involved in gene regulation. They modulate gene expression at various levels, such as transcription, translation, splicing, and mRNA degradation. Because riboswitches exhibit relatively complex structures, they are able to form highly complex ligand-binding sites, which enable the specific recognition of target metabolites in a complex cellular environment. Practically in all studied cases, riboswitches use ligand-induced conformational changes to control gene expression. To monitor the structural reorganization of riboswitches, we use the local fluorescent reporter 2-aminopurine (2AP), which is a structural analog of adenine. The 2AP fluorescence is strongly quenched when the fluorophore is involved in stacking interactions with surrounding bases, and can, therefore, be used to monitor local structural rearrangements. Here, we show specific examples in which 2AP fluorescence can be used to monitor structural changes in the aptamer domain of the S-adenosyl methionine (SAM) riboswitch and where it can be used as a ligand for the guanine riboswitch.

Nonlocal helix formation is key to understanding S-adenosylmethionine-1 riboswitch function.

Riboswitches are noncoding RNAs that regulate gene expression in response to changing concentrations of specific metabolites. Switching activity is affected by the interplay between the aptamer domain and expression platform of the riboswitch. The aptamer domain binds the metabolite, locking the riboswitch in a ligand-bound conformation. In absence of the metabolite, the expression platform forms an alternative secondary structure by sequestering the 3' end of a nonlocal helix called P1. We use all-atom structure-based simulations to characterize the folding, unfolding, and metabolite binding of the aptamer domain of the S-adenosylmethionine-1 (SAM-1) riboswitch. Our results suggest that folding of the nonlocal helix (P1) is rate-limiting in aptamer domain formation. Interestingly, SAM assists folding of the P1 helix by reducing the associated free energy barrier. Because the 3' end of the P1 helix is sequestered by an alternative helix in the absence of metabolites, this observed ligand-control of P1 formation provides a mechanistic explanation of expression platform regulation.

Transcriptome-wide analysis of uncapped mRNAs in Arabidopsis reveals regulation of mRNA degradation.

The composition of the transcriptome is determined by a balance between mRNA synthesis and degradation. An important route for mRNA degradation produces uncapped mRNAs, and this decay process can be initiated by decapping enzymes, endonucleases, and small RNAs. Although uncapped mRNAs are an important intermediate for mRNA decay, their identity and abundance have never been studied on a large scale until recently. Here, we present an experimental method for transcriptome-wide profiling of uncapped mRNAs that can be used in any eukaryotic system. We applied the method to study the prevalence of uncapped transcripts during the early stages of Arabidopsis thaliana flower development. Uncapped transcripts were identified for the majority of expressed genes, although at different levels. By comparing uncapped RNA levels with steady state overall transcript levels, our study provides evidence for widespread mRNA degradation control in numerous biological processes involving genes of varied molecular functions, implying that uncapped mRNA levels are dynamically regulated. Sequence analyses identified structural features of transcripts and cis-elements that were associated with different levels of uncapping. These transcriptome-wide profiles of uncapped mRNAs will aid in illuminating new regulatory mechanisms of eukaryotic transcriptional networks.

[Analysis of correlation of local GC level in human protein coding genes].

GC level is an important feature of genomic composition, which significantly improve our understanding of structure, function and evolution of genes. In this paper, the nonredundant DNA sequence of 7,992 human protein coding genes were retrieved from public database and the local GC level of different sequence regions and correlation between GC levels were analyzed.. The results showed that the GC levels of different sequence regions were strikingly nonuniform. 5' untranslated regions were of richest GC, with average GC content being 62.5%. 3'-untranslated regions were of poorest GC, with average GC content being 43.97%. GC contents of 3' flanking sequences profoundly matched the GC levels of DNA large fragments where the genes were located. Although the GC contents of open reading frames (ORFs) were higher than that of intron, 3' non-translated region and 3' flanking sequences, high correlation existed among the GC contents of the four regions. Average GC content of the third codon position (GC3) was 58.9%, higher than that of the fist and second position, and showed high correlation to GC contents of ORFs, with correlation coefficients being 0.91, besides of its significant association with GC contents of intron, 3'-untranslated region and 3' flanking sequences. Moreover, the linear regression of GC3 against GC contents of 3' flanking sequences yielded a slope of 1.25. Thus, GC3 was a sensitive indicator for GC change of local genome. As for 5' flanking sequences, 5' untranslated regions, fist and second codon position, however, their GC level exhibited weaker correlation with that of other regions. These results suggest that the third codon positions, introns, 3'-untranslated regions and 3' flanking sequences may evolve similarly while first and second codon positions, 5' flanking sequences and 5' untranslated region were expected to bear more selective stress for holding their functions.

Selective recognition of a DNA G-quadruplex by an engineered antibody.

Particular guanine rich nucleic acid sequences can fold into stable secondary structures called G-quadruplexes. These structures have been identified in various regions of the genome that include the telomeres, gene promoters and UTR regions, raising the possibility that they may be associated with biological function(s). Computational analysis has predicted that intramolecular G-quadruplex forming sequences are prevalent in the human genome, thus raising the desire to differentially recognize genomic G-quadruplexes. We have employed antibody phage display and competitive selection techniques to generate a single-chain antibody that shows >1000-fold discrimination between G-quadruplex and duplex DNA, and furthermore >100-fold discrimination between two related intramolecular parallel DNA G-quadruplexes. The amino acid sequence composition at the antigen binding site shows conservation within the light and heavy chains of the selected scFvs, suggesting sequence requirements for G-quadruplex recognition. Circular dichroism (CD) spectroscopic data showed that the scFv binds to the prefolded G-quadruplex and does not induce G-quadruplex structure formation. This study demonstrates the strongest discrimination that we are aware of between two intramolecular genomic G-quadruplexes.

Molecular cloning and characterization of the chicken cationic amino acid transporter-2 gene.

The system y(+) cationic amino acid transporters (CATs) play an important role in the regulation of lysine and arginine utilization. We have characterized the genomic organization and tissue transcription of the chicken CAT-2 (cCAT-2) isoforms from pectoralis muscle. The primary cCAT-2 transcript is alternatively spliced within the same position of the mRNA to produce cCAT-2A and cCAT-2B isoforms that share high nucleotide homology with their mammalian counterparts. We also identified a novel third CAT-2 isoform, cCAT-2C. This isoform contains a premature termination codon and studies in the chicken LMH cell line show that cCAT-2C mRNA is degraded by the nonsense-mediated mRNA decay pathway. All three cCAT-2 isoforms are transcribed in liver, pectoralis, gastrocnemius and heart. Chicken cCAT-2A was the predominant cCAT-2 isoform in liver, pectoralis and gastrocnemius. Analysis of the 5'-untranslated region of cCAT-2 identified multiple cCAT-2 promoters. Chicken CAT-2 promoter usage was tissue dependent and was not responsible for cCAT-2 isoform production. Analysis of genomic sequence upstream cCAT-2 promoter regions revealed binding sites for transcription factors involved in amino acid sensing, hormone signaling and immune function. These results provide insight into the role of the cCAT-2 gene and its regulation of lysine and arginine utilization in aves.

Comparative analysis of full genomic sequences among different genotypes of dengue virus type 3.

BACKGROUND: Although the previous study demonstrated the envelope protein of dengue viruses is under purifying selection pressure, little is known about the genetic differences of full-length viral genomes of DENV-3. In our study, complete genomic sequencing of DENV-3 strains collected from different geographical locations and isolation years were determined and the sequence diversity as well as selection pressure sites in the DENV genome other than within the E gene were also analyzed. RESULTS: Using maximum likelihood and Bayesian approaches, our phylogenetic analysis revealed that the Taiwan's indigenous DENV-3 isolated from 1994 and 1998 dengue/DHF epidemics and one 1999 sporadic case were of the three different genotypes - I, II, and III, each associated with DENV-3 circulating in Indonesia, Thailand and Sri Lanka, respectively. Sequence diversity and selection pressure of different genomic regions among DENV-3 different genotypes was further examined to understand the global DENV-3 evolution. The highest nucleotide sequence diversity among the fully sequenced DENV-3 strains was found in the nonstructural protein 2A (mean +/- SD: 5.84 +/- 0.54) and envelope protein gene regions (mean +/- SD: 5.04 +/- 0.32). Further analysis found that positive selection pressure of DENV-3 may occur in the non-structural protein 1 gene region and the positive selection site was detected at position 178 of the NS1 gene. CONCLUSION: Our study confirmed that the envelope protein is under purifying selection pressure although it presented higher sequence diversity. The detection of positive selection pressure in the non-structural protein along genotype II indicated that DENV-3 originated from Southeast Asia needs to monitor the emergence of DENV strains with epidemic potential for better epidemic prevention and vaccine development.

Comparative full genome analysis revealed E1: A226V shift in 2007 Indian Chikungunya virus isolates.

The resurgence of Chikungunya virus (CHIKV) in the form of unprecedented explosive epidemic after a gap of 32 years in India is a point of major public health concern. In 2007 again there was outbreak in Kerala, India, affecting more than 25,000 cases with many reported mortalities. To understand the molecular basis of this high virulence and its implication in large-scale epidemic, a detailed systematic serological, virological and molecular investigation was undertaken with the epidemic samples of Kerala-2007. The comparative analysis of full genome sequence of Chikungunya virus isolate of 2007 with 2006 revealed three unique substitutions in structural and non-structural genes of 2007 isolate [two in E1 region (V14A and A226V) and one in Nsp1 (M184T)]. Our finding further substantiates the association of A226V shift in E1 gene with evolutionary success possibly due to adaptation in the mosquito vector with progression of epidemic, as observed in Reunion Island. This A226V shift which was absent in all 2006 Indian isolates, is found to be present in the four 2007 isolates, analysed in this study. These unique molecular features of the 2007 isolates with the progression of the epidemic from 2005 to 2007 demonstrate their high evolutionary and epidemic potential and thereby suggesting possible implication in higher magnitude and virulence of this outbreak.

Bioinformatics approaches for studying untranslated regions of mRNAs.

Interactions between RNA-binding proteins and cis-acting elements in the 5'- and 3'-untranslated regions (UTRs) of transcripts are responsible for regulating essential biological activities, such as mRNA localization, mRNA turnover, and translation efficiency. This chapter introduces some of the publicly available free bioinformatics resources, including software tools and databases, which can be used for predicting, mapping, and characterizing regulatory motifs found in the eukaryotic mRNA-untranslated regions.

Complex riboswitches.

Using simple biochemical tricks, metabolite-binding riboswitches take on gene control functions that have long been thought to be the work of protein factors. Although modern riboswitches might be the last holdouts of primitive genetic elements, some are capable of sensory and regulatory feats that are competitive with their protein counterparts.

Genomic sequence of mandarin fish rhabdovirus with an unusual small non-transcriptional ORF.

The complete genome of mandarin fish Siniperca chuatsi rhabdovirus (SCRV) was cloned and sequenced. It comprises 11,545 nucleotides and contains five genes encoding the nucleoprotein N, the phosphoprotein P, the matrix protein M, the glycoprotein G, and the RNA-dependent RNA polymerase protein L. At the 3' and 5' termini of SCRV genome, leader and trailer sequences show inverse complementarity. The N, P, M and G proteins share the highest sequence identities (ranging from 14.8 to 41.5%) with the respective proteins of rhabdovirus 903/87, the L protein has the highest identity with those of vesiculoviruses, especially with Chandipura virus (44.7%). Phylogenetic analysis of L proteins showed that SCRV clustered with spring vireamia of carp virus (SVCV) and was most closely related to viruses in the genus Vesiculovirus. In addition, an overlapping open reading frame (ORF) predicted to encode a protein similar to vesicular stomatitis virus C protein is present within the P gene of SCRV. Furthermore, an unoverlapping small ORF downstream of M ORF within M gene is predicted (tentatively called orf4). Therefore, the genomic organization of SCRV can be proposed as 3' leader-N-P/C-M-(orf4)-G-L-trailer 5'. Orf4 transcription or translation products could not be detected by northern or Western blot, respectively, though one similar mRNA band to M mRNA was found. This is the first report on one small unoverlapping ORF in M gene of a fish rhabdovirus.

GRSDB2 and GRS_UTRdb: databases of quadruplex forming G-rich sequences in pre-mRNAs and mRNAs.

G-quadruplex motifs in the RNA play significant roles in key cellular processes and human disease. While sequences capable of forming G-quadruplexes in the pre-mRNA are involved in regulation of polyadenylation and splicing events in mammalian transcripts, the G-quadruplex motifs in the UTRs may help regulate mRNA expression. GRSDB2 is a second-generation database containing information on the composition and distribution of putative Quadruplex-forming G-Rich Sequences (QGRS) mapped in approximately 29 000 eukaryotic pre-mRNA sequences, many of which are alternatively processed. The data stored in the GRSDB2 is based on computational analysis of NCBI Entrez Gene entries with the help of an improved version of the QGRS Mapper program. The database allows complex queries with a wide variety of parameters, including Gene Ontology terms. The data is displayed in a variety of formats with several additional computational capabilities. We have also developed a new database, GRS_UTRdb, containing information on the composition and distribution patterns of putative QGRS in the 5'- and 3'-UTRs of eukaryotic mRNA sequences. The goal of these experiments has been to build freely accessible resources for exploring the role of G-quadruplex structure in regulation of gene expression at post-transcriptional level. The databases can be accessed at the G-Quadruplex Resource Site at: http://bioinformatics.ramapo.edu/GQRS/.

Ligand-induced folding of the thiM TPP riboswitch investigated by a structure-based fluorescence spectroscopic approach.

Riboswitches are genetic control elements within non-coding regions of mRNA. They consist of a metabolite-sensitive aptamer and an adjoining expression platform. Here, we describe ligand-induced folding of a thiamine pyrophosphate (TPP) responsive riboswitch from Escherichia coli thiM mRNA, using chemically labeled variants. Referring to a recent structure determination of the TPP/aptamer complex, each variant was synthesized with a single 2-aminopurine (AP) nucleobase replacement that was selected to monitor formation of tertiary interactions of a particular region during ligand binding in real time by fluorescence experiments. We have determined the rate constants for conformational adjustment of the individual AP sensors. From the 7-fold differentiation of these constants, it can be deduced that tertiary contacts between the two parallel helical domains (P2/J3-2/P3/L3 and P4/P5/L5) that grip the ligand's ends in two separate pockets, form significantly faster than the function-critical three-way junction with stem P1 fully developed. Based on these data, we characterize the process of ligand binding by an induced fit of the RNA and propose a folding model of the TPP riboswitch aptamer. For the full-length riboswitch domain and for shorter constructs that represent transcriptional intermediates, we have additionally evaluated ligand-induced folding via AP-modified variants and provide insights into the sequential folding pathway that involves a finely balanced equilibrium of secondary structures.

Ligand recognition determinants of guanine riboswitches.

Guanine riboswitches negatively modulate transcription upon guanine binding. The aptamer domain is organized around a three-way junction which forms the ligand binding site. Using currently available 89 guanine aptamer sequences, a consensus secondary structure is deduced and reveals differences from the previously identified aptamer consensus. Three positions are found to display different nucleotide requirements. Using a 2-aminopurine binding assay, we show that variations are allowed depending on the aptamer context. However, changes at position 48 markedly decrease ligand binding in a context-independent fashion. This is consistent with previous observations with the adenine riboswitch in which position 48 was proposed to interact with position 74, which normally base pairs with the ligand. The in vivo transcriptional control of endogenous Bacillus subtilis guanine riboswitches was studied using RT-qPCR assays. The ratio of elongated/terminated transcripts is decreased in presence of a high concentration of guanine but is dependent on the riboswitch analyzed. In general, the aptamer-2AP complex affinity correlates well with the in vivo regulation efficiency of the corresponding riboswitch. These studies suggest that core variations of guanine aptamers are used to produce a spectrum of ligand binding affinities which is used in vivo by host riboswitches to perform gene regulation.

Genetic adaptation to untranslated region-mediated enterovirus growth deficits by mutations in the nonstructural proteins 3AB and 3CD.

Both untranslated regions (UTRs) of plus-strand RNA virus genomes jointly control translation and replication of viral genomes. In the case of the Enterovirus genus of the Picornaviridae family, the 5'UTR consists of a cloverleaf-like terminus preceding the internal ribosomal entry site (IRES) and the 3' terminus is composed of a structured 3'UTR and poly(A). The IRES and poly(A) have been implicated in translation control, and all UTR structures, in addition to cis-acting genetic elements mapping to the open reading frame, have been assigned roles in RNA replication. Viral UTRs are recognized by viral and host cell RNA-binding proteins that may co-determine genome stability, translation, plus- and minus-strand RNA replication, and scaffolding of viral replication complexes within host cell substructures. In this report, we describe experiments with coxsackie B viruses with a cell type-specific propagation deficit in Sk-N-Mc neuroblastoma cells conferred by the combination of a heterologous IRES and altered 3'UTR. Serial passage of these constructs in Sk-N-Mc cells yielded genetic adaptation by mutations within the viral nonstructural proteins 3A and 3C. Our data implicate 3A and/or 3C or their precursors 3AB and/or 3CD in a functional complex with the IRES and 3'UTR that drives viral propagation. Adaptation to neuroblastoma cells suggests an involvement of cell type-specific host factors or the host cell cytoplasmic milieu in this phenomenon.

Conformational changes in the expression domain of the Escherichia coli thiM riboswitch.

The thiM riboswitch contains an aptamer domain that adaptively binds the coenzyme thiamine pyrophosphate (TPP). The binding of TPP to the aptamer domain induces structural rearrangements that are relayed to a second domain, the so-called expression domain, thereby interfering with gene expression. The recently solved crystal structures of the aptamer domains of the thiM riboswitches in complex with TPP revealed how TPP stabilizes secondary and tertiary structures in the RNA ligand complex. To understand the global modes of reorganization between the two domains upon metabolite binding the structure of the entire riboswitch in presence and absence of TPP needs to be determined. Here we report the secondary structure of the entire thiM riboswitch from Escherichia coli in its TPP-free form and its transition into the TPP-bound variant, thereby depicting domains of the riboswitch that serve as communication links between the aptamer and the expression domain. Furthermore, structural probing provides an explanation for the lack of genetic control exerted by a riboswitch variant with mutations in the expression domain that still binds TPP.

The 3a accessory protein of SARS coronavirus specifically interacts with the 5'UTR of its genomic RNA, Using a unique 75 amino acid interaction domain.

More than four years have passed since the outbreak of the severe acute respiratory syndrome (SARS) epidemic, and still very little is known about the molecular biology and pathogenesis of this deadly virus. Among the accessory proteins of the SARS coronavirus (SARS-CoV), the 3a protein has been shown to interact with the spike, envelope, and membrane glycoprotein and has recently been established to be a structural component of capsid. Recent studies suggest that the 3a protein may function as an ion channel and may promote virus release. In order to further characterize the functional properties of this protein, we initiated studies to check its RNA binding activity. Using the yeast three-hybrid system, electrophoretic mobility shift assay (EMSA), and ultraviolet (UV) cross-linking techniques, we have shown that the 3a protein is capable of binding specifically to the 5' untranslated region (5'UTR) of the SARS virus genomic RNA. Further, we have mapped the interaction domain of the 3a protein responsible for this RNA-protein interaction using a series of deletion mutants and defined it to the central 75 amino acid region. This RNA binding motif of 3a does not share homology with any other known RNA binding protein and may have an important role in viral capsid assembly and pathogenesis.