Sequencing and analysis for the full-length genome RNA of foot-and-mouth disease virus China/99

The complete nucleotide sequence of genomic RNA of foot and mouth disease virus (FMDV) strain China/99 from infected bovine tongue epithelium is presented. The nucleotide sequence extending from the 5' end of the genomic RNA to the 5' end of poly (A) tail contains 8173 nucleotides (nt). Its open reading frame, which encodes a single polypeptide of 2332 amino acids, encompasses 6999 nt starting from the initiation codon AUG and terminating at the UAA codon 93 bases upstream from the 5' end of poly (A) tract. The 5' untranslated region (UTR) is composed of 1081 nt. The consensus of the 1d gene of FMDV strain China/99 compared with that of UKG/6/2001, UKG/12/2001, China/99HN4 and China/3/Tibet is over 97%. The result showed the stains belong to the members of the Pan-Asia family. There is a remarkable differentiation in the function-unknown (FUR), p2 and p3 regions between FMDV isolates from infected cattle and swine, especially in 3a gene. No deletion was found in genes /, 1a, 1b, 2a, 2c, 3b, and 3d. Thes     

口蹄疫病毒O/NY00株基因组L片段的克隆及其基因特征

采用RT—PCR方法对口蹄疫病毒O／NY00株基因组L片段进行了分子克隆和序列测定。结果表明：获得的O／NY00株基因组L片段长7805nt,其中包括715nt的5’非编码区和6999nt的多聚蛋白编码区。将O／NY00株与其它参考毒株进行序列比较,结果显示：O／NY00株与O／SKR／2000、O／1JKG／3／2001遗传关系最近,而与其它毒株遗传关系较远,属于0型口蹄疫ME-SA拓扑型Pan／ksia株的成员。与Cathay拓扑型的代表毒株O／TAW／97比较,两者在主要抗原位点1和位点3呈现出明显的以拓扑型为特征的氨基酸差别,提示两拓扑型毒株之间有较高的抗原差异。     

Detection and characterization of defective interfering RNAs associated with the cocksfoot mottle sobemovirus

A new RNA of about 900 nt was found in the virions of cocksfoot mottle virus (CfMV) and in infected plants by RNA hybridization and RT-PCR. Structural features suggested that this RNA is a defective interfering RNA (diRNA). The CfMV diRNA was shown to consist of a 35-nt 5′-terminal genomic region, which formed a hairpin, and a 3′-terminal genomic region, which included the coat protein (CP) gene lacking the first 120 nt.In vitro translation of the diRNA started at the third Met codon to produce truncated CP. The CfMV diRNA was assumed totrans-activate synthesis of the CP subgenomic RNA (sgRNA).     

Expression Enhancement of a Rice Polyubiquitin Gene Promoter

An 808 bp promoter from a rice polyubiquitin gene, rubi3, has been isolated. The rubi3 gene contained an open reading frame of 1140 bp encoding a pentameric polyubiquitin arranged as five tandem, head-to-tail repeats of 76 aa. The 1140 bp 5′ UTR intron of the gene enhanced its promoter activity in transient expression assays by 20-fold. Translational fusion of the GUS reporter gene to the coding sequence of the ubiquitin monomer enhanced GUS enzyme activity in transient expression assays by 4.3-fold over the construct containing the original rubi3 promoter (including the 5′ UTR intron) construct. The enhancing effect residing in the ubiquitin monomer coding sequence has been narrowed down to the first 9 nt coding for the first three amino acid residues of the ubiquitin protein. Mutagenesis at the third nucleotide of this 9 nt sequence still maintains the enhancing effect, but leads to translation of the native GUS protein rather than a fusion protein. The resultant 5′ regulatory sequence, consisting of the rubi3 promoter, 5′ UTR exon and intron, and the mutated first 9 nt coding sequence, has an activity nearly 90-fold greater than the rubi3 promoter only (without the 5′ UTR intron), and 2.2-fold greater than the maize Ubi1 gene promoter (including its 5′ UTR intron). The newly created expression vector is expected to enhance transgene expression in monocot plants. Considering the high conservation of the polyubiquitin gene structure in higher plants, the observed enhancement in gene expression may apply to 5′ regulatory sequences of other plant polyubiquitin genes.     

Characterization of proteins associating with 5’ terminus of PGHS-1 mRNA

Induction of Prostaglandin Endoperoxide H Synthase-1 (PGHS-1) gene has been previously documented in a few studies during events such as development and cellular differentiation. However, molecular mechanisms governing the regulation of PGHS-1 gene expression and contributing to changes in protein levels are poorly understood. Using the MEG-01 cell model of PGHS-1 gene induction, our laboratory has previously demonstrated that the 5’UTR and the first two exons of PGHS-1 mRNA had a significant impact on decreasing the translational efficiency of a reporter gene and suggested that the presence of a secondary structure is required for conservation of this activity. This 5’end of PGHS-1 mRNA sequence has also been shown to associate with nucleolin protein. In the current study, we set to investigate the protein composition of the mRNP (messenger ribonucleoprotein) associating with the 5’end of PGHS-1 mRNA and to identify its protein members. RNA/protein binding assays coupled with LC-MS analysis identified serpin B1 and NF45 (nuclear factor 45) proteins as potential members of PGHS-1 mRNP complex. Immunoprecipitation experiments using MEG-01 protein extracts validated mass spectrometry data and confirmed binding of nucleolin, serpin B1, NF45 and NF90. The RNA fraction was extracted from immunoprecipitated mRNP complexes and association of RNA binding proteins, serpin B1, NF45 and NF90, to PGHS-1 mRNA target sequence was confirmed by RT-PCR. Together these data suggest that serpin B1, NF45 and NF90 associate with PGHS-1 mRNA and can potentially participate in the formation a single or a number of PGHS-1 ribonucleoprotein complexes, through nucleolin that possibly serves as a docking base for other protein complex members.     

Evolution of Prokaryotic Genes by Shift of Stop Codons

De novo origin of coding sequence remains an obscure issue in molecular evolution. One of the possible paths for addition (subtraction) of DNA segments to (from) a gene is stop codon shift. Single nucleotide substitutions can destroy the existing stop codon, leading to uninterrupted translation up to the next stop codon in the gene’s reading frame, or create a premature stop codon via a nonsense mutation. Furthermore, short indels-caused frameshifts near gene’s end may lead to premature stop codons or to translation past the existing stop codon. Here, we describe the evolution of the length of coding sequence of prokaryotic genes by change of positions of stop codons. We observed cases of addition of regions of 3′UTR to genes due to mutations at the existing stop codon, and cases of subtraction of C-terminal coding segments due to nonsense mutations upstream of the stop codon. Many of the observed stop codon shifts cannot be attributed to sequencing errors or rare deleterious variants segregating within bacterial populations. The additions of regions of 3′UTR tend to occur in those genes in which they are facilitated by nearby downstream in-frame triplets which may serve as new stop codons. Conversely, subtractions of coding sequence often give rise to in-frame stop codons located nearby. The amino acid composition of the added region is significantly biased, compared to the overall amino acid composition of the genes. Our results show that in prokaryotes, shift of stop codon is an underappreciated contributor to functional evolution of gene length.     

The factors dictating the codon usage variation among the genes in the genome of Burkholderia pseudomallei

Burkholderia pseudomallei is a recognized biothreat agent and the causative agent of melioidosis. Codon usage biases of all protein-coding genes (length greater than or equal to 300 bp) from the complete genome of B. pseudomallei K96243 have been analyzed. As B. pseudomallei is a GC-rich organism (68.5%), overall codon usage data analysis indicates that indeed codons ending in G and/or C are predominant in this organism. But multivariate statistical analysis indicates that there is a single major trend in the codon usage variation among the genes in this organism, which has a strong positively correlation with the expressivities of the genes. The majority of the lowly expressed genes are scattered towards the negative end of the major axis whereas the highly expressed genes are clustered towards the positive end. At the same time, from the results that there were two significant correlations between axis 1 coordinates and the GC, GC3s content at silent sites of each sequence, and clearly significant negatively correlations between the ‘Effective Number of Codons’ values and GC, GC3s content, we inferred that codon usage bias was affected by gene nucleotide composition also. In addition, some other factors such as the lengths of the genes as well as the hydrophobicity of genes also influence the codon usage variation among the genes in this organism in a minor way. At the same time, notably, 21 codons have been defined as ‘optimal codons’ of the B. pseudomallei. In summary, our work have provided a basic understanding of the mechanisms for codon usage bias and some more useful information for improving the expression of target genes in vivo and in vitro. Sheng Zhao and Qin Zhang contributed equally to this work.     

Differential mRNA stability to endogenous ribonucleases of the coding region and 3′ untranslated regions of wheat (Triticum aestivum L.) manganese superoxide dismutase genes

The sequences of the 3′ untranslated region (UTR) of the manganese superoxide dismutase (MnSOD) genes in wheat (Triticum aestivum) were found to be quite variable with different predicted thermostabilities. The degradation rates of the 3′ UTR variants and the coding region were measured following exposure to endogenous nucleases. The degradation rates of the 3′ UTR variants for 15 min were not significantly different, meaning the degradation rates of the 3′ UTR variants were not directly related to the thermostabilities. However, the degradation rate of the coding region was significantly faster than those of the 3’ UTR variants. Further investigation revealed the coding region seemed to have specific sites for degradation, indicating a possibility of increasing MnSOD expression by the degradation site alteration.     

IRES-driven translation is stimulated separately by the FMDV 3'-NCR and poly(A) sequences

The 3′ end region of foot-and-mouth disease virus (FMDV) consists of two distinct elements, a 90 nt untranslated region (3′-NCR) and a poly(A) tract. Removal of either the poly(A) tract or both the 3′-NCR and the poly(A) tract abrogated infectivity in susceptible cells in the context of a full-length cDNA clone. We have addressed the question of whether the impairment of RNA infectivity is related to defects at the translation level using a double approach. First, compared to the full-length viral RNA, removal of the 3′ sequences reduced the efficiency of translation in vitro. Secondly, a stimulatory effect of the 3′ end sequences on IRES-dependent translation was found in vivo using bicistronic constructs. RNAs carrying the FMDV 3′ end sequences linked to the second cistron showed a significant stimulation of IRES-dependent translation, whereas cap-dependent translation was not affected. Remarkably, IRES-dependent stimulation exerted by the poly(A) tract or the 3′-NCR seems to be the result of two separate events, as the 3′-NCR alone enhanced IRES activity on its own. Under conditions of FMDV Lb protease-induced translation shut-off, the stimulation of IRES activity reached values above 6-fold in living cells. A northern blot analysis indicated that IRES stimulation was not the consequence of a change in the stability of the bicistronic RNA produced in transfected cells. Analysis of the RNA-binding proteins interacting with a mixture of 3′ end and IRES probes showed an additive pattern. Altogether, our results strongly suggest that individual signals in the viral 3′ end ensure stimulation of FMDV translation.     

Sequences within Both the 5′ UTR and Gag Are Required for Optimal In Vivo Packaging and Propagation of Mouse Mammary Tumor Virus (MMTV) Genomic RNA

Background

This study mapped regions of genomic RNA (gRNA) important for packaging and propagation of mouse mammary tumor virus (MMTV). MMTV is a type B betaretrovirus which preassembles intracellularly, a phenomenon distinct from retroviruses that assemble the progeny virion at cell surface just before budding such as the type C human and feline immunodeficiency viruses (HIV and FIV). Studies of FIV and Mason-Pfizer monkey virus (MPMV), a type D betaretrovirus with similar intracellular virion assembly processes as MMTV, have shown that the 5′ untranslated region (5′ UTR) and 5′ end of gag constitute important packaging determinants for gRNA.

Methodology

Three series of MMTV transfer vectors containing incremental amounts of gag or 5′ UTR sequences, or incremental amounts of 5′ UTR in the presence of 400 nucleotides (nt) of gag were constructed to delineate the extent of 5′ sequences that may be involved in MMTV gRNA packaging. Real time PCR measured the packaging efficiency of these vector RNAs into MMTV particles generated by co-transfection of MMTV Gag/Pol, vesicular stomatitis virus envelope glycoprotein (VSV-G Env), and individual transfer vectors into human 293T cells. Transfer vector RNA propagation was monitored by measuring transduction of target HeLaT4 cells following infection with viral particles containing a hygromycin resistance gene expression cassette on the packaged RNA.

Principal Findings

MMTV requires the entire 5′ UTR and a minimum of ∼120 nucleotide (nt) at the 5′ end of gag for not only efficient gRNA packaging but also propagation of MMTV-based transfer vector RNAs. Vector RNAs without the entire 5′ UTR were defective for both efficient packaging and propagation into target cells.

Conclusions/Significance

These results reveal that the 5′ end of MMTV genome is critical for both gRNA packaging and propagation, unlike the recently delineated FIV and MPMV packaging determinants that have been shown to be of bipartite nature.     

The LOX-1 3′UTR188CT polymorphism and coronary artery disease in Turkish patients

In coronary artery disease (CAD), a potentially reversible factor leading to cardiac death is left ventricular hypertrophy (LVH). The 3′untranslated region (3′UTR) 188CT polymorphism of lectin-like oxidized low-density lipoproteins receptor-1 (LOX-1) gene has been associated with an increased risk for CAD. We aim to investigate, in a Turkish population, whether 3′UTR188CT variation could affect the development of LVH in CAD patients. In a population-based case–control study, we compared 83 cases with CAD and 99 healthy controls for this polymorphism. The LOX-1 3′UTR188CT genotypes were determined by PCR–RFLP technique. LOX-1 3′UTR188 TT genotype was associated with significantly increased systolic blood pressure (P = 0.047) and risk of LVH (P = 0.014, OR: 3.541) when compared with the C allele carriers. In addition, the TT genotype was positively associated with decreased levels of HDL-cholesterol in the control subjects (P = 0.031) and increased levels of VLDL-C in the patient group (P = 0.009). The LOX-1 3′UTR188CT gene polymorphism may predispose to the development of LVH in CAD patients, dependent on blood pressure.     

A universal protocol to generate consensus level genome sequences for foot-and-mouth disease virus and other positive-sense polyadenylated RNA viruses using the Illumina MiSeq

Background

Next-Generation Sequencing (NGS) is revolutionizing molecular epidemiology by providing new approaches to undertake whole genome sequencing (WGS) in diagnostic settings for a variety of human and veterinary pathogens. Previous sequencing protocols have been subject to biases such as those encountered during PCR amplification and cell culture, or are restricted by the need for large quantities of starting material. We describe here a simple and robust methodology for the generation of whole genome sequences on the Illumina MiSeq. This protocol is specific for foot-and-mouth disease virus (FMDV) or other polyadenylated RNA viruses and circumvents both the use of PCR and the requirement for large amounts of initial template.

Results

The protocol was successfully validated using five FMDV positive clinical samples from the 2001 epidemic in the United Kingdom, as well as a panel of representative viruses from all seven serotypes. In addition, this protocol was successfully used to recover 94% of an FMDV genome that had previously been identified as cell culture negative. Genome sequences from three other non-FMDV polyadenylated RNA viruses (EMCV, ERAV, VESV) were also obtained with minor protocol amendments. We calculated that a minimum coverage depth of 22 reads was required to produce an accurate consensus sequence for FMDV O. This was achieved in 5 FMDV/O/UKG isolates and the type O FMDV from the serotype panel with the exception of the 5′ genomic termini and area immediately flanking the poly(C) region.

Conclusions

We have developed a universal WGS method for FMDV and other polyadenylated RNA viruses. This method works successfully from a limited quantity of starting material and eliminates the requirement for genome-specific PCR amplification. This protocol has the potential to generate consensus-level sequences within a routine high-throughput diagnostic environment.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-828) contains supplementary material, which is available to authorized users.     

Strong Eukaryotic IRESs Have Weak Secondary Structure

Background

The objective of this work was to investigate the hypothesis that eukaryotic Internal Ribosome Entry Sites (IRES) lack secondary structure and to examine the generality of the hypothesis.

Methodology/Principal Findings

IRESs of the yeast and the fruit fly are located in the 5′UTR immediately upstream of the initiation codon. The minimum folding energy (MFE) of 60 nt RNA segments immediately upstream of the initiation codons was calculated as a proxy of secondary structure stability. MFE of the reverse complements of these 60 nt segments was also calculated. The relationship between MFE and empirically determined IRES activity was investigated to test the hypothesis that strong IRES activity is associated with weak secondary structure. We show that IRES activity in the yeast and the fruit fly correlates strongly with the structural stability, with highest IRES activity found in RNA segments that exhibit the weakest secondary structure.

Conclusions

We found that a subset of eukaryotic IRESs exhibits very low secondary structure in the 5′-UTR sequences immediately upstream of the initiation codon. The consistency in results between the yeast and the fruit fly suggests a possible shared mechanism of cap-independent translation initiation that relies on an unstructured RNA segment.     

Correlation Between Shine–Dalgarno Sequence Conservation and Codon Usage of Bacterial Genes

In this study, we analyzed the correlation between codon usage bias and Shine–Dalgarno (SD) sequence conservation, using complete genome sequences of nine prokaryotes. For codon usage bias, we adopted the codon adaptation index (CAI), which is based on the codon usage preference of genes encoding ribosomal proteins, elongation factors, heat shock proteins, outer membrane proteins, and RNA polymerase subunit proteins. To compute SD sequence conservation, we used SD motif sequences predicted by Tompa and systematically aligned them with 5′UTR sequences. We found that there exists a clear correlation between the CAI values and SD sequence conservation in the genomes of Escherichia coli, Bacillus subtilis, Haemophilus influenzae, Archaeoglobus fulgidus, Methanobacterium thermoautotrophicum, and Methanococcus jannaschii, and no relationship is found in M. genitalium, M. pneumoniae, and Synechocystis. That is, genes with higher CAI values tend to have more conserved SD sequences than do genes with lower CAI values in these organisms. Some organisms, such as M. thermoautotrophicum, do not clearly show the correlation. The biological significance of these results is discussed in the context of the translation initiation process and translation efficiency. Received: 22 June 2000 / Accepted: 18 October 2000