Identification of high-confidence human poly(A) RNA isoform scaffolds using nanopore sequencing

Nanopore sequencing devices read individual RNA strands directly. This facilitates identification of exon linkages and nucleotide modifications; however, using conventional direct RNA nanopore sequencing, the 5' and 3' ends of poly(A) RNA cannot be identified unambiguously. This is due in part to RNA degradation in vivo and in vitro that can obscure transcription start and end sites. In this study, we aimed to identify individual full-length human RNA isoforms among ∼4 million nanopore poly(A)-selected RNA reads. First, to identify RNA strands bearing 5' m⁷G caps, we exchanged the biological cap for a modified cap attached to a 45-nt oligomer. This oligomer adaptation method improved 5' end sequencing and ensured correct identification of the 5' m⁷G capped ends. Second, among these 5'-capped nanopore reads, we screened for features consistent with a 3' polyadenylation site. Combining these two steps, we identified 294,107 individual high-confidence full-length RNA scaffolds from human GM12878 cells, most of which (257,721) aligned to protein-coding genes. Of these, 4876 scaffolds indicated unannotated isoforms that were often internal to longer, previously identified RNA isoforms. Orthogonal data for m⁷G caps and open chromatin, such as CAGE and DNase-HS seq, confirmed the validity of these high-confidence RNA scaffolds.

An allelic series at the KARRIKIN INSENSITIVE 2 locus of Arabidopsis thaliana decouples ligand hydrolysis and receptor degradation from downstream signalling

Karrikins are butenolide compounds present in post-fire environments that can stimulate seed germination in many species, including Arabidopsis thaliana. Plants also produce endogenous butenolide compounds that serve as hormones, namely strigolactones (SLs). The receptor for karrikins (KARRIKIN INSENSITIVE 2; KAI2) and the receptor for SLs (DWARF14; D14) are homologous proteins that share many similarities. The mode of action of D14 as a dual enzyme receptor protein is well established, but the nature of KAI2-dependent signalling and its function as a receptor are not fully understood. To expand our knowledge of how KAI2 operates, we screened ethyl methanesulphonate (EMS)-mutagenized populations of A. thaliana for mutants with kai2-like phenotypes and isolated 13 new kai2 alleles. Among these alleles, kai2-10 encoded a D184N protein variant that was stable in planta. Differential scanning fluorimetry assays indicated that the KAI2 D184N protein could interact normally with bioactive ligands. We developed a KAI2-active version of the fluorescent strigolactone analogue Yoshimulactone Green to show that KAI2 D184N exhibits normal rates of ligand hydrolysis. KAI2 D184N degraded in response to treatment with exogenous ligands, suggesting that receptor degradation is a consequence of ligand binding and hydrolysis, but is insufficient for signalling activity. Remarkably, KAI2 D184N degradation was hypersensitive to karrikins, but showed a normal response to strigolactone analogues, implying that these butenolides may interact differently with KAI2. These results demonstrate that the enzymatic and signalling functions of KAI2 can be decoupled, and provide important insights into the mechanistic events that underpin butenolide signalling in plants.

Polyadenylation Site-Based Analysis of Transcript Expression by 3'READS

Deep sequencing of the 3' end region of poly(A)+ RNA identifies the cleavage and polyadenylation site (PAS) and measures transcript abundance. However, mispriming at internal A-rich regions by the oligo-dT oligo in reverse transcription can lead to falsely identified PASs. This problem can be resolved by direct ligation of an adapter to the 3' end of RNA. However, ligation-based methods are often inefficient. Here, we describe 3'READS+, an accurate and sensitive method for deep sequencing of the 3' end of poly(A)+ RNA. Through partial digestion by RNase H of the poly(A) tail bound to a locked nucleic acid (LNA)/DNA hybrid oligo, this method sequences an optimal number of terminal A's, which balances sequencing quality and accurate identification of PAS in A-rich regions. With efficient ligation steps, 3'READS+ is amenable to small amounts of input RNA. 3'READS+ can also be readily used as a cost-effective method for gene expression analysis.

Genome-wide dynamics of alternative polyadenylation in rice

Alternative polyadenylation (APA), in which a transcript uses one of the poly(A) sites to define its 3'-end, is a common regulatory mechanism in eukaryotic gene expression. However, the potential of APA in determining crop agronomic traits remains elusive. This study systematically tallied poly(A) sites of 14 different rice tissues and developmental stages using the poly(A) tag sequencing (PAT-seq) approach. The results indicate significant involvement of APA in developmental and quantitative trait loci (QTL) gene expression. About 48% of all expressed genes use APA to generate transcriptomic and proteomic diversity. Some genes switch APA sites, allowing differentially expressed genes to use alternate 3' UTRs. Interestingly, APA in mature pollen is distinct where differential expression levels of a set of poly(A) factors and different distributions of APA sites are found, indicating a unique mRNA 3'-end formation regulation during gametophyte development. Equally interesting, statistical analyses showed that QTL tends to use APA for regulation of gene expression of many agronomic traits, suggesting a potential important role of APA in rice production. These results provide thus far the most comprehensive and high-resolution resource for advanced analysis of APA in crops and shed light on how APA is associated with trait formation in eukaryotes.

A novel TP53 variant (rs78378222 A > C) in the polyadenylation signal is associated with increased cancer susceptibility: evidence from a meta-analysis

Polymorphisms in TP53 are involved in the progression of different types of cancer. A rare novel TP53 variant (rs78378222 A > C allele) was found via whole-genome sequencing in 2011. This variant was shown to significantly increase the risk of glioma, colorectal adenoma and prostate cancer. Functional analysis further revealed that this variant hindered TP53 expression and its downstream effect on apoptosis. Several studies have investigated the relationship between rs78378222 and cancer susceptibility. However, the results were not consistent. We conducted the first meta-analysis to give a more credible assessment on the association about this variant and cancer risk. Our meta-analysis included 34 studies consisting of 36599 cases and 91272 controls. These studies were mostly on the basis of high-grade data from Genome-wide association studies (GWASs). The results indicated that TP53 rs78378222 was significantly associated with an increased risk of overall cancer (AC vs. AA: OR = 1.511, 95% CI = 1.285-1.777). Furthermore, stratified analyses indicated that rs78378222 increased the risk of nervous system cancer, skin cancer and other cancer. To summarize, this meta-analysis suggested that rs78378222 C allele is a potent risk factor for overall cancer.

Mutated polyadenylation signals for controlling expression levels of multiple genes in mammalian cells

A set of mutated SV40 early polyadenylation signals (SV40pA) with varying strengths is generated by mutating the AATAAA sequence in the wild-type SV40pA. They are shown to control the expression level of a gene over a 10-fold range using luciferase reporter genes in transient transfection assays. The relative strength of these SV40pA variants remains similar under three commonly used mammalian promoters and in five mammalian cell lines. Application of SV40pA variants for controlling expression level of multiple genes is demonstrated in a study of monoclonal antibody (mAb) synthesis in mammalian cells. By using SV40pA variants of different strengths, the expression of light chain (LC) and heavy chain (HC) genes encoded in a single vector is independently altered which results in different ratios of LC to HC expression spanning a range from 0.24 to 16.42. The changes in gene expression are determined by measuring mRNA levels and intracellular LC and HC polypeptides. It is found that a substantial decrease of HC expression, which increases the LC/HC mRNA ratio, only slightly reduces mAb production. However, reducing the LC expression by a similar magnitude, which decreases the LC/HC mRNA ratio results in a sharp decline of mAb production to trace amounts. This set of SV40pA variants offers a new tool for accurate control of the relative expression levels of multiple genes. It will have wide-ranging applications in fields related to the study of biosynthesis of multi-subunit proteins, proteomic research on protein interactions, and multi-gene metabolic engineering.

Multiple RNAs from the mouse carboxypeptidase M locus: functional RNAs or transcription noise?

BACKGROUND

A major effort of the scientific community has been to obtain complete pictures of the genomes of many organisms. This has been accomplished mainly by annotation of structural and functional elements in the genome sequence, a process that has been centred in the gene concept and, as a consequence, biased toward protein coding sequences. Recently, the explosion of transcriptome data generated and the discovery of many functional non-protein coding RNAs have painted a more detailed and complex scenario for the genome. Here we analyzed the mouse carboxypeptidase M locus in this broader perspective in order to define the mouse CPM gene structure and evaluate the existence of other transcripts from the same genomic region.

RESULTS

Bioinformatic analysis of nucleotide sequences that map to the mouse CPM locus suggests that, in addition to the mouse CPM mRNA, it expresses at least 33 different transcripts, many of which seem to be non-coding RNAs. We randomly chose to evaluate experimentally four of these extra transcripts. They are expressed in a tissue specific manner, indicating that they are not artefacts or transcriptional noise. Furthermore, one of these four extra transcripts shows expression patterns that differed considerably from the other ones and from the mouse CPM gene, suggesting that there may be more than one transcriptional unit in this locus. In addition, we have confirmed the mouse CPM gene RefSeq sequence by rapid amplification of cDNA ends (RACE) and directional cloning.

CONCLUSION

This study supports the recent view that the majority of the genome is transcribed and that many of the resulting transcripts seem to be non-coding RNAs from introns of genes or from independent transcriptional units. Although some of the information on the transcriptome of many organisms may actually be artefacts or transcriptional noise, we argue that it can be experimentally evaluated and used to find and define biological functional elements on the genome. Furthermore, the transcription of other functional RNAs besides the protein coding RNA from a specific genomic locus imposes extra care when designing and interpreting experiments involving genetic manipulations or expression detection and quantification.

Specific trans-acting proteins interact with auxiliary RNA polyadenylation elements in the COX-2 3'-UTR

Two cyclooxygenase (COX) enzymes, COX-1 and COX-2, are present in human cells. While COX-1 is constitutively expressed, COX-2 is inducible and up-regulated in response to many signals. Since increased transcriptional activity accounts for only part of COX-2 up-regulation, we chose to explore other RNA processing mechanisms in the regulation of this gene. Previously, we showed that COX-2 is regulated by alternative polyadenylation, and that the COX-2 proximal polyadenylation signal contains auxiliary upstream sequence elements (USEs) that are very important in efficient polyadenylation. To explore trans-acting protein factors interacting with these cis-acting RNA elements, we performed pull-down assays with HeLa nuclear extract and biotinylated RNA oligonucleotides representing COX-2 USEs. We identified PSF, p54(nrb), PTB, and U1A as proteins specifically bound to the COX-2 USEs. We further explored their participation in polyadenylation using MS2 phage coat protein-MS2 RNA binding site tethering assays, and found that tethering any of these four proteins to the COX-2 USE mutant RNA can compensate for these cis-acting elements. Finally, we suggest that these proteins (p54(nrb), PTB, PSF, and U1A) may interact as a complex since immunoprecipitations of the transfected MS2 fusion proteins coprecipitate the other proteins.

The poly(A)-dependent transcriptional pause is mediated by CPSF acting on the body of the polymerase

Eukaryotic poly(A) signals direct mRNA 3'-end processing and also pausing and termination of transcription. We show that pausing and termination require the processing factor CPSF, which binds the AAUAAA hexamer of the mammalian poly(A) signal. Pausing does not require the RNA polymerase II C-terminal domain (CTD) or the cleavage stimulation factor, CstF, that binds the CTD. Pull-down experiments show that CPSF binds, principally through its 30-kDa subunit, to the body of the polymerase. CPSF can also bind CstF, but this seems to be mutually exclusive with polymerase binding. We suggest that CPSF, while binding the body of the polymerase, scans for hexamers in the extruding RNA. Any encounter with a hexamer triggers pausing. If the hexamer is part of a functional poly(A) signal, CstF is recruited and binds CPSF, causing it to release the polymerase body and move (with CstF) to the CTD.

Comparative analysis of the full genome sequence of European bat lyssavirus type 1 and type 2 with other lyssaviruses and evidence for a conserved transcription termination and polyadenylation motif in the G-L 3' non-translated region

We report the first full-length genomic sequences for European bat lyssavirus type-1 (EBLV-1) and type-2 (EBLV-2). The EBLV-1 genomic sequence was derived from a virus isolated from a serotine bat in Hamburg, Germany, in 1968 and the EBLV-2 sequence was derived from a virus isolate from a human case of rabies that occurred in Scotland in 2002. A long-distance PCR strategy was used to amplify the open reading frames (ORFs), followed by standard and modified RACE (rapid amplification of cDNA ends) techniques to amplify the 3' and 5' ends. The lengths of each complete viral genome for EBLV-1 and EBLV-2 were 11 966 and 11 930 base pairs, respectively, and follow the standard rhabdovirus genome organization of five viral proteins. Comparison with other lyssavirus sequences demonstrates variation in degrees of homology, with the genomic termini showing a high degree of complementarity. The nucleoprotein was the most conserved, both intra- and intergenotypically, followed by the polymerase (L), matrix and glyco- proteins, with the phosphoprotein being the most variable. In addition, we have shown that the two EBLVs utilize a conserved transcription termination and polyadenylation (TTP) motif, approximately 50 nt upstream of the L gene start codon. All available lyssavirus sequences to date, with the exception of Pasteur virus (PV) and PV-derived isolates, use the second TTP site. This observation may explain differences in pathogenicity between lyssavirus strains, dependent on the length of the untranslated region, which might affect transcriptional activity and RNA stability.

PolyA_DB 2: mRNA polyadenylation sites in vertebrate genes

Polyadenylation of nascent transcripts is one of the key mRNA processing events in eukaryotic cells. A large number of human and mouse genes have alternative polyadenylation sites, or poly(A) sites, leading to mRNA variants with different protein products and/or 3′-untranslated regions (3′-UTRs). PolyA_DB 2 contains poly(A) sites identified for genes in several vertebrate species, including human, mouse, rat, chicken and zebrafish, using alignments between cDNA/ESTs and genome sequences. Several new features have been added to the database since its last release, including syntenic genome regions for human poly(A) sites in seven other vertebrates and cis-element information adjacent to poly(A) sites. Trace sequences are used to provide additional evidence for poly(A/T) tails in cDNA/ESTs. The updated database is intended to broaden poly(A) site coverage in vertebrate genomes, and provide means to assess the authenticity of poly(A) sites identified by bioinformatics. The URL for this database is .

Systematic variation in mRNA 3'-processing signals during mouse spermatogenesis

Gene expression and processing during mouse male germ cell maturation (spermatogenesis) is highly specialized. Previous reports have suggested that there is a high incidence of alternative 3′-processing in male germ cell mRNAs, including reduced usage of the canonical polyadenylation signal, AAUAAA. We used EST libraries generated from mouse testicular cells to identify 3′-processing sites used at various stages of spermatogenesis (spermatogonia, spermatocytes and round spermatids) and testicular somatic Sertoli cells. We assessed differences in 3′-processing characteristics in the testicular samples, compared to control sets of widely used 3′-processing sites. Using a new method for comparison of degenerate regulatory elements between sequence samples, we identified significant changes in the use of putative 3′-processing regulatory sequence elements in all spermatogenic cell types. In addition, we observed a trend towards truncated 3′-untranslated regions (3′-UTRs), with the most significant differences apparent in round spermatids. In contrast, Sertoli cells displayed a much smaller trend towards 3′-UTR truncation and no significant difference in 3′-processing regulatory sequences. Finally, we identified a number of genes encoding mRNAs that were specifically subject to alternative 3′-processing during meiosis and postmeiotic development. Our results highlight developmental differences in polyadenylation site choice and in the elements that likely control them during spermatogenesis.

An exploration of 3'-end processing signals and their tissue distribution in Oryza sativa

The 3' untranslated regions deeply affect many properties of eukaryotic mRNA. In plants, the polyadenine control signals contained in these regions seem to be more variable than of mammals. Three cDNA libraries derived from the leaf, endosperm and stem tissues of rice were sequenced from the 3'-end. Of the 9911 transcripts analyzed, 5723 unique transcripts were identified from the leaf sequences, 2934 from the endosperm and 1254 from the stem. The information entropy and two statistical methods were used to compile a list of rice poly(A) control signals. Based on their distribution, these signals can be roughly grouped into far-upstream element (FUE), near-upstream element (NUE), T-rich region (TRE) and downstream element (DE). The distribution of rice conserved regions is similar to the previous model from Arabidopsis and yeast, with a few differences in word constructions. Interestingly, we also found the word distributions were diverse in the cleavage site of downstream sequences of different rice tissues. The signal bias in downstream sequences may lead mRNA to be differently cleaved in different rice tissues.

Expression of TPK1 and TPK2 genes encoding PKA catalytic subunits during growth and morphogenesis in Candida albicans

The transcript levels of Candida albicans TPK1 and TPK2 genes, encoding PKA catalytic subunits, as well as phosphotransferase activity, were measured in the parental strain CAI4 and in homozygous tpk1Delta and tpk2Delta mutants during vegetative growth and during yeast-to-mycelial transition in N-acetylglucosamine liquid inducing medium at 37 degrees C. We observed two TPK2 transcripts, a major one of 1.8 kb and a minor one of 1.4 kb, and established by 3'-RACE that they originate from the recognition of the three polyadenylation signals present in the 3' untranslated region of the gene. During vegetative growth of CAI4 strain, the expression profiles of TPK1 and TPK2 varied similarly, reaching maximal expression at the late logarithmic phase. TPK1 mRNA levels were lower than those of TPK2 at all stages measured. In the corresponding homozygous tpk mutants, mRNA levels and the expression patterns of TPK1 and TPK2 were similar to those of CAI4, suggesting that the loss of one catalytic isoform is not compensated by overexpression of the other. Changes in PKA specific activity roughly correlated with fluctuations of mRNA expression levels. During yeast-to-mycelial transition, a sharp increase in TPK1 mRNA levels and in PKA-specific activity correlated with the onset of germ-tube formation in strain tpk2Delta. We also showed that tpk1Delta strain exhibited a delayed morphogenetic shift in comparison with CAI4 and tpk2Delta strains in several liquid inducing media, reinforcing the idea that Tpk1p is important for faster germ-tube appearance.

Characterization of an abundant COL9A1 transcript in the cochlea with a novel 3' UTR: Expression studies and detection of miRNA target sequence

EST N66408 represents one of several large unique clusters expressed in the Morton human fetal cochlear cDNA library. N66408 is 575 bp in size and initial BLAST analysis of this sequence showed no homology to any known genes or expressed sequence tags (ESTs) from other organs or tissues. Sequence of the original cochlear clone from which N66408 was derived revealed that the corresponding cDNA was about 700 bp in size, including 125 bp at its 5' end with homology to the 3' end of COL9A1 in addition to 575 bp of novel sequence. RT-PCR analysis using primers specific to COL9A1 isoforms 1 and 2 detected expression of both isoforms in human fetal cochlea. Tissue in situ hybridization using the novel 3' UTR sequence as probe showed abundant expression in spiral limbus and spiral ligament, and a moderate level of expression in the organ of Corti. dbEST analysis of ESTs specific to the 3' UTR of COL9A1 showed 19 ESTs derived from various tissues; three polyadenylation sites were identified and the majority of these ESTs were derived from overlapping polyadenylation signals at the second site (position 749-758). Comparison of the 3' UTR of human COL9A1 with its orthologs as well as with dbEST uncovered a highly conserved region around the overlapping polyadenylation signals at position 749-758 in mammals. A search of the microRNA database revealed a highly conserved target sequence for miR-9 immediately preceding the overlapping polyadenylation signals in the novel 3' UTR of COL9A1, suggesting its role in posttranscriptional regulation of COL9A1.

Characteristics of transcription-regulatory elements for gene expression from plasmid vectors in human trophoblast cell lines

Nonviral gene delivery systems are useful for basic research in trophoblasts. In these systems, gene expression is regulated by a cassette of regulatory elements within the plasmid, and the transcriptional activity differs among cell lines. In the present study, we used BeWo and JAR human trophoblast cell lines to systematically compare the transcriptional activities of several expression cassettes and those of a control plasmid made up of a simian virus 40 (SV40) promoter, a polyadenylation (PA) signal, and an enhancer. We also found that insertion of intron elements enhanced transcriptional activities in the following order: intron A>hybrid beta-globin-immunoglobin intron>no intron. Of several PA signals tested including those from SV40, bovine growth hormone, and the minimal rabbit beta-globin, the latter had the highest transcriptional activities (3.9- and 26-fold over control plasmid in BeWo and JAR cells, respectively). Addition of a second enhancer increased the transcriptional activity in these cells. We also found that gene expression level can be controlled by selecting the expression cassette. These results should be useful for further transgene experiments in BeWo and JAR cells.

The hepatitis C virus 3'-untranslated region or a poly(A) tract promote efficient translation subsequent to the initiation phase

Enhancement of eukaryotic messenger RNA (mRNA) translation initiation by the 3′ poly(A) tail is mediated through interaction of poly(A)-binding protein with eukaryotic initiation factor (eIF) 4G, bridging the 5′ terminal cap structure. In contrast to cellular mRNA, translation of the uncapped, non-polyadenylated hepatitis C virus (HCV) genome occurs independently of eIF4G and a role for 3′-untranslated sequences in modifying HCV gene expression is controversial. Utilizing cell-based and in vitro translation assays, we show that the HCV 3′-untranslated region (UTR) or a 3′ poly(A) tract of sufficient length interchangeably stimulate translation dependent upon the HCV internal ribosomal entry site (IRES). However, in contrast to cap-dependent translation, the rate of initiation at the HCV IRES was unaffected by 3′-untranslated sequences. Analysis of post-initiation events revealed that the 3′ poly(A) tract and HCV 3′-UTR improve translation efficiency by enabling termination and possibly ribosome recycling for successive rounds of translation.

GRSDB: a database of quadruplex forming G-rich sequences in alternatively processed mammalian pre-mRNA sequences

Guanine-rich nucleic acids are known to form highly stable G-quadruplex structures, also known as G-quartets. Recently, there has been a tremendous amount of interest in studying G-quadruplexes owing to the realization of their biological importance. G-rich sequences (GRSs) capable of forming G-quadruplexes are found in the vicinity of polyadenylation regions and are involved in regulating 3′ end processing of mammalian pre-mRNAs. G-rich motifs are also known to play an important role in alternative, tissue-specific splicing by interacting with hnRNP H protein subfamily. Whether quadruplex structure directly plays a role in regulating RNA processing events requires further investigation. To date there has not been a comprehensive effort to study G-quadruplexes near RNA processing sites. We have applied a computational approach to map putative Quadruplex forming GRSs within the transcribed regions of a large number of alternatively processed human and mouse gene sequences that were obtained as fully annotated entries from GenBank and RefSeq. We have used the computed data to build the GRSDB database that provides a unique avenue for studying G-quadruplexes in the context of RNA processing sites. GRSDB website offers visual comparison of G-quadruplex distribution patterns among all the alternative RNA products of a gene with the help of dynamic graphics. At present, GRSDB contains data from 1310 human and mouse genes, of which 1188 are alternatively processed. It has a total of 379 223 predicted G-quadruplexes, of which 54 252 are near RNA processing sites. GRSDB is a good resource for researchers interested in investigating the functional relevance of G-quadruplexes, especially in the context of alternative RNA processing. It can be accessed at .

Identification of post-transcriptionally regulated Xenopus tropicalis maternal mRNAs by microarray

Cytoplasmic control of the adenylation state of mRNAs is a critical post-transcriptional process involved in the regulation of mRNAs stability and translational efficiency. The early development of Xenopus laevis has been a major model for the study of such regulations. We describe here a microarray analysis to identify mRNAs that are regulated by changes in their adenylation state during oogenesis and early development of the diploid frog Xenopus tropicalis. The microarray data were validated using qRT–PCR and direct analysis of the adenylation state of endogenous maternal mRNAs during the period studied. We identified more than 500 mRNAs regulated at the post-transcriptional level among the 3000 mRNAs potentially detected by the microarray. The mRNAs were classified into nine different adenylation behavior categories. The various adenylation profiles observed during oocyte maturation and early development and the analyses of 3′-untranslated region sequences suggest that previously uncharacterized sequence elements control the adenylation behavior of the newly identified mRNAs. These data should prove useful in identifying mRNAs with important functions during oocyte maturation and early development.

Molecular Characterization of Major Histocompatibility Complex Class I (B-F) mRNA Variants from Chickens Differing in Resistance to Marek's Disease

In this study, the relative distributions of two alternatively polyadenylated chicken major histocompatibility complex (MHC) mRNA isoforms of approximately 1.5 and 1.9 kb were analysed in spleen cells from chickens homozygous for the MHC haplotypes B21 and B19v1 as well as in heterozygous B19v1/B21 birds. Both isoforms are likely to encode classical MHC class I (B-F) alpha chains. The B19v1 and B21 MHC haplotypes confer different levels of protection against Marek's disease (MD), which is caused by infection with MD virus (MDV). In spleen cells, MD-resistant B21 birds were shown to have the highest percentage of the 1.5 kb variant relative to the total MHC class I expression, MD-susceptible B19v1 birds the lowest and B19v1/B21 birds an intermediate percentage. Infection of 4-week-old chickens with the GA strain of MDV was shown to cause a significant increase in the relative amount of 1.5 kb transcripts in B21 birds 32 days postinfection (dpi). Alternatively polyadenylated mRNA isoforms may encode identical proteins, but differences in the 3' untranslated region (UTR) can influence polyadenylation, mRNA stability, intracellular localization and translation efficiency. It was shown that the increased 1.5 kb percentage in B21 birds 32 days postinfection may be a result of a change in the choice of poly(A) site rather than a locus-specific upregulated transcription of the BF1 gene that preferentially expresses the 1.5 kb variant. Furthermore, the 3' end of the 1.5 kb mRNA variants deriving from B19v1 and B21 chickens was characterized by Rapid Amplification of cDNA Ends (RACE) and sequencing. No potentially functional elements were identified in the 3' UTR of the RACE products corresponding to this short isoform. However, variation in polyadenylation site was observed between the BF1 and BF2 mRNA transcripts and alternative splicing-out of the sequence (exon 7) encoding the second segment of the cytoplasmic part of the mature BF2*19 molecules. This alternative exon 7 splice variant was also detected in other MD-susceptible haplotypes, but not in the MD-resistant B21 and B21-like haplotypes, suggesting a potential role of exon 7 in MHC-related MD resistance.

Genomic analysis of G protein gamma subunits in human and mouse - the relationship between conserved gene structure and G protein betagamma dimer formation

Analysis of the genomic sequences, cDNAs and expressed sequence tags (ESTs) in human and mouse for the 12 genes of the gamma subunits of the heterotrimeric G proteins has allowed us to identify the common versus unique elements of the organization and expression of the members of this important gene family. All of the G protein gamma subunit genes are organized around two coding exons, each containing about 100 nucleotides coding for 30-40 amino acids. These two exons each correspond to a functional domain of the protein, which interestingly appears to impose constraints on both the structure of the protein and the structure of the gene. There is large variation in the intron size between these two coding exons, the number and size of 5' and 3' UTRs, and the overall size of the genes. There is general but not absolute conservation in the size and structure of these genes between humans and mice. Alternative splicing and potential differential promoter usage were detected for several Ggamma subunits, indicating possible differential regulation in expression. Only for Ggamma10, however, did we find an alternative coding transcript. This alternative transcript appears to code for a hybrid protein containing a DnaJ domain in place of its Ggamma exon 1 domain, joined to the Ggamma10 second exon domain. The predicted mRNA is expressed in humans, and the protein coded by it is readily translated in vitro. This protein does not form a functional G protein betagamma dimer, but it could generate a chaperone-like protein related to its DNA-J domain. These studies suggest that alternative splicing is not a prominent mechanism for generating G protein subunit diversity from within the human or mouse genomes. Instead, each of the known 12 gamma subunit genes generate transcripts with one prevalent protein.

In vitro correction of medium chain acyl CoA dehydrogenase deficiency with a recombinant adenoviral vector

Defects of mitochondrial beta-oxidation are a growing group of disorders with variable clinical presentations ranging from mild hypotonia to sudden infant death. Current therapy involves avoidance of fasting, dietary restrictions, and cofactor supplementation. Unfortunately, times of acute illness and noncompliance can interfere with these therapies and result in a rapid clinical decline. The development of a safe, durable, and effective gene delivery system remains an attractive alternative therapy for individuals with these disorders. To this end, a recombinant first-generation adenovirus vector (Ad/cmv-hMCAD) has been prepared that constitutively expresses the human medium chain acyl CoA dehydrogenase (MCAD) protein under the control of the CMV promoter and bovine polyadenylation signal. Characterization of human fibroblasts deficient in MCAD infected with Ad/cmv-hMCAD including Western analysis, immunohistological staining visualized with confocal microscopy, electron transfer protein (ETF) reduction assay, and palmitate loading studies was performed. Infection of MCAD deficient fibroblast with Ad/cmv-hmcad resulted in the production of a 55kDa protein that co-localized in cells with a mitochondrial marker. Extracts prepared from Ad/cmv-hMCAD infected deficient fibroblasts demonstrated correction of the block seen in the MCAD catalyzed reduction of ETF in the presence of octanoyl CoA. Finally, MCAD deficient fibroblasts infected with increasing amounts of Ad/cmv-hMCAD showed a stepwise improvement of the abnormal acylcarnitine profile exhibited by the deficient cells. Together these studies demonstrate our ability to express and monitor the expression of MCAD in treated cells and support further in vivo murine studies to assess toxicity and duration of correction with this and other MCAD recombinant vectors.

A 57-nucleotide upstream early polyadenylation element in human papillomavirus type 16 interacts with hFip1, CstF-64, hnRNP C1/C2, and polypyrimidine tract binding protein

We have investigated the role of the human papillomavirus type 16 (HPV-16) early untranslated region (3' UTR) in HPV-16 gene expression. We found that deletion of the early 3' UTR reduced the utilization of the early polyadenylation signal and, as a consequence, resulted in read-through into the late region and production of late L1 and L2 mRNAs. Deletion of the U-rich 3' half of the early 3' UTR had a similar effect, demonstrating that the 57-nucleotide U-rich region acted as an enhancing upstream element on the early polyadenylation signal. In accordance with this, the newly identified hFip1 protein, which has been shown to enhance polyadenylation through U-rich upstream elements, interacted specifically with the HPV-16 upstream element. This upstream element also interacted specifically with CstF-64, hnRNP C1/C2, and polypyrimidine tract binding protein, suggesting that these factors were either enhancing or regulating polyadenylation at the HPV-16 early polyadenylation signal. Mutational inactivation of the early polyadenylation signal also resulted in increased late mRNA production. However, the effect was reduced by the activation of upstream cryptic polyadenylation signals, demonstrating the presence of additional strong RNA elements downstream of the early polyadenylation signal that direct cleavage and polyadenylation to this region of the HPV-16 genome. In addition, we identified a 3' splice site at genomic position 742 in the early region with the potential to produce E1 and E4 mRNAs on which the E1 and E4 open reading frames are preceded only by the suboptimal E6 AUG. These mRNAs would therefore be more efficiently translated into E1 and E4 than previously described HPV-16 E1 and E4 mRNAs on which E1 and E4 are preceded by both E6 and E7 AUGs.

A probabilistic model of 3' end formation in Caenorhabditis elegans

The 3' ends of mRNAs terminate with a poly(A) tail. This post-transcriptional modification is directed by sequence features present in the 3'-untranslated region (3'-UTR). We have undertaken a computational analysis of 3' end formation in Caenorhabditis elegans. By aligning cDNAs that diverge from genomic sequence at the poly(A) tract, we accurately identified a large set of true cleavage sites. When there are many transcripts aligned to a particular locus, local variation of the cleavage site over a span of a few bases is frequently observed. We find that in addition to the well-known AAUAAA motif there are several regions with distinct nucleotide compositional biases. We propose a generalized hidden Markov model that describes sequence features in C.elegans 3'-UTRs. We find that a computer program employing this model accurately predicts experimentally observed 3' ends even when there are multiple AAUAAA motifs and multiple cleavage sites. We have made available a complete set of polyadenylation site predictions for the C.elegans genome, including a subset of 6570 supported by aligned transcripts.

Secondary structure as a functional feature in the downstream region of mammalian polyadenylation signals

Secondary structure within the downstream region of mammalian polyadenylation signals has been proposed to perform important functions. The simian virus 40 late polyadenylation signal (SVLPA) forms alternate secondary structures in equilibrium. Their formation correlates with cleavage-polyadenylation efficiency (H. Hans and J. C. Alwine, Mol. Cell. Biol. 20:2926-2932, 2000; M. I. Zarudnaya, I. M. Kolomiets, A. L. Potyahaylo, and D. M. Hovorun, Nucleic Acids Res. 3:1375-1386, 2003), and oligonucleotides that disrupt the secondary structure inhibit in vitro cleavage. To define the important features of downstream secondary structure, we first minimized the SVLPA by deletion, forming a downstream region with fewer, and more stable, stem-loop structures. Specific mutagenesis showed that both stem stability and loop size are important functional features of the downstream region. Stabilization of the stem, thus minimizing alternative structures, decreased cleavage efficiency both in vitro and in vivo. This was most deleterious when the stem was stabilized at the base of the loop, constraining loop size by inhibiting breathing of the stem. The significance of loop size was supported by mutants that showed increased cleavage efficiency with increased loop size and vice versa. A loop of at least 12 nucleotides promoted cleavage; U richness in the loop also promoted cleavage and was particularly important when the stem was stabilized. A mutation designed to eliminate downstream secondary structure still formed many relatively weak alternative structures in equilibrium and retained function. The data suggest that although the downstream region is very important, its structure is quite malleable and is able to tolerate significant mutation within a wide range of primary and secondary structural features. We propose that this malleability is due to the enhanced ability of GU- and U-rich downstream elements to easily form secondary structures with surrounding sequences.

Genomic organization and transcription of the human retinol dehydrogenase 10 (RDH10) gene

A cDNA clone up-regulated in hydraulic lung edema in rabbit showed high similarity with human RDH10 mRNA, which encodes a protein involved in retinoic acid metabolism. We defined the organization of the human gene, which includes a unique transcriptional start site, a coding region with six translated exons and a 3' untranslated region containing at least two used polyadenylation sites. The two poly(A) signals are responsible for the production of the 3 and 4 kb RDH10 mRNA isoforms detected in several human tissues and cell lines.

Sequence analysis of a 5.1 kbp region of the Spodoptera frugiperda multicapsid nucleopolyhedrovirus genome that comprises a functional ecdysteroid UDP-glucosyltransferase (egt) gene

The ecdysteroid UDP-glucosyltransferase gene from the Spodoptera frugiperda multicapsid nucleopolyhedrovirus (SfMNPV) was identified using degenerate primers whose sequence were derived from conserved regions of the EGT proteins encoded by other baculoviruses. Analysis of the gene sequence revealed the presence of an open reading frame (ORF) with potential to encode a polypeptide of 525 amino acids. Promoter sequences typical of baculovirus genes were found in the 5' region of this ORF. A polyadenylation signal was identified downstream the translation stop codon. A transient expression assay showed that the product of this ORF was able to conjugate glucose from UDP-glucose with ecdysone confirming that the gene identified was indeed the SfMNPV egt gene. The SfMNPV egt gene and the sequences of other baculovirus egt genes were used to infer a phylogenetic tree. The nucleotide sequence of the entire BamHI fragment that contains the SfMNPV egt gene was determined. Search of the available sequence databases suggested that, besides the egt gene, this region contains 5 ORFs similar to the baculovirus genes gp37 (fusolin), to ptp2 and to ORFs 28, 29, and 30 of Spodoptera exigua multicapsid nucleopolyhedrovirus. Both the phylogenetic analysis of the egt genes and the gene order of the region that flanks the egt gene indicated that SfMNPV is closely related to the baculoviruses that infects S. exigua and Mamestra configurata.

ELAV inhibits 3'-end processing to promote neural splicing of ewg pre-mRNA

The embryonic lethal abnormal visual system (ELAV) is a gene-specific regulator of alternative pre-mRNA processing in neurons of Drosophila. Here we define a functional in vivo binding site for ELAV in neurons through the development of a reporter gene system in transgenic animals in combination with in vitro binding assays. ELAV binds to erect wing (ewg) RNA 3' of a polyadenylation site in the terminal intron 6. At this polyadenylation site, ELAV inhibits 3'-end processing in vitro in a dose-dependent and sequence-specific manner, and ELAV binding is necessary in vivo to promote splicing of ewg intron 6. Further, the AAUAAA poly(A) complex recognition sequence, together with ELAV, is required to regulate neural 3' splice site choice in vivo. In addition, the use of segmentally labeled RNA substrates in UV cross-linking assays suggest that ELAV does not inhibit or redirect binding of cleavage factor dCstF64 at the regulated polyadenylation site on ewg RNA. These data indicate that binding of 3'-end processing factors, together with ELAV, can regulate alternative splicing.

Comparative analysis of the effects of packaging signal, transgene orientation, promoters, polyadenylation signals, and E3 region on growth properties of first-generation adenoviruses

First-generation adenovectors have been developed for gene therapy and vaccine applications. The construction of these adenovectors has entailed the use of numerous types of expression cassettes. It has long been known that first-generation adenovectors can be rescued more easily and to higher titers with some transgenes than with others. This study has systematically shown that there can be marked differences in growth properties of recombinant adenovectors attributable to the use of promoters, the orientation of the transgene within the E1A/E1B-deleted region, and the inclusion of the E3 region. In addition, we had demonstrated the benefit of extending the packaging signal region to include elements V, VI, and VII. The effects of the complete packaging region were studied by plasmid competition studies between original and modified adenovectors. Similar competition studies between E3(+) and E3(-) adenovectors were performed and showed that the E3(+) vector had a growth advantage over its E3(-) counterpart. By making various changes, we have enhanced the growth capacity of our recombinant adenovector by more than 3-fold under serum-free and cell suspension growth conditions. Along with this enhanced growth, our adenovectors have maintained their genetic stability after 21 successive passages in cell culture. This increased robustness will be critical when adapting first-generation recombinant adenovectors to commercial production.

An uncapped RNA suggests a model for Caenorhabditis elegans polycistronic pre-mRNA processing

Polycistronic pre-mRNAs from Caenohabditis elegans operons are processed by internal cleavage and polyadenylation to create 3' ends of mature mRNAs. This is accompanied by trans-splicing with SL2 approximately 100 nucleotides downstream of the 3' end formation sites to create the 5' ends of downstream mRNAs. SL2 trans-splicing depends on a U-rich element (Ur), located approximately 70 nucleotides upstream of the trans-splice site in the intercistronic region (ICR), as well as a functional 3' end formation signal. Here we report the existence of a novel gene-length RNA, the Ur-RNA, starting just upstream of the Ur element. The expression of Ur-RNA is dependent on 3' end formation as well as on the presence of the Ur element, but does not require a trans-splice site. The Ur-RNA is not capped, and alteration of the location of the Ur element in either the 5' or 3' direction alters the location of the 5' end of the Ur-RNA. We propose that a 5' to 3' exonuclease degrades the precursor RNA following cleavage at the poly(A) site, stopping when it reaches the Ur element, presumably attributable to a bound protein. Part of the function of this protein can be performed by the MS2 coat protein. Recruitment of coat protein to the ICR in the absence of the Ur element results in accumulation of an RNA equivalent to Ur-RNA, and restores trans-splicing. Only SL1, however, is used. Therefore, coat protein is sufficient for blocking the exonuclease and thereby allowing formation of a substrate for trans-splicing, but it lacks the ability to recruit the SL2 snRNP. Our results also demonstrate that MS2 coat protein can be used as an in vivo block to an exonuclease, which should have utility in mRNA stability studies.

Characterization of the murine CD30 ligand (CD153) gene: gene structure and expression

CD153 (CD30 ligand) has been described as a 40-kDa type II transmembrane glycoprotein belonging to the TNF superfamily and is expressed primarily by activated T cells, B cells and monocytes. In this study, we have determined that the murine CD153 gene consists of four exons, with three intervening introns, spaced over approximately 26 kb of genomic sequence. Sequence analysis of the murine CD153 promoter and 5' flanking region revealed the presence of a TATA box element immediately upstream of two tsp sites, together with putative binding motifs for a variety of lymphoid-specific transcription factors. 5'RACE analysis of LPS-stimulated RAW264.7 macrophage cDNA identified at least four transcriptional start sites for murine CD153, with two sites occurring downstream of the previously predicted translation initiation codon. Additionally, 5' RACE analysis identified multiple murine CD153 polyadenylation sites. Our results indicate that primary murine CD153 transcripts may vary from 26 kb to approximately 28 kb in length.

Reexamining the polyadenylation signal: were we wrong about AAUAAA?

Polyadenylation is the process by which most eukaryotic mRNAs form their 3' ends. It was long held that polyadenylation required the sequence AAUAAA and that 90% of mRNAs had AAUAAA within 30 nucleotides of the site of poly(A) addition. More recent studies, aided by computer analysis of sequences made available in GenBank and expressed sequence tag (EST) databases, have suggested that the actual incidence of AAUAAA is much lower, perhaps as low as 50-60%. Reproductive biologists have long recognized that a large number of mRNAs in male germ cells of mammals lack AAUAAA but are otherwise normally polyadenylated. Recent research in our laboratory has uncovered a new form of an essential polyadenylation protein, tauCstF-64, that is most highly expressed in male germ cells, and to a smaller extent in the brain, and which we propose plays a significant role in AAUAAA-independent mRNA polyadenylation in germ cells.