Evidence that BmTXK beta-BmKCT cDNA from Chinese scorpion Buthus martensii Karsch is an artifact generated in the reverse transcription process

Artifacts and biases of the reverse transcription reaction in RNA sequencing

RNA sequencing has spurred a significant number of research areas in recent years. Most protocols rely on synthesizing a more stable complementary DNA (cDNA) copy of the RNA molecule during the reverse transcription reaction. The resulting cDNA pool is often wrongfully assumed to be quantitatively and molecularly similar to the original RNA input. Sadly, biases and artifacts confound the resulting cDNA mixture. These issues are often overlooked or ignored in the literature by those that rely on the reverse transcription process. In this review, we confront the reader with intra- and intersample biases and artifacts caused by the reverse transcription reaction during RNA sequencing experiments. To fight the reader's despair, we also provide solutions to most issues and inform on good RNA sequencing practices. We hope the reader can use this review to their advantage, thereby contributing to scientifically sound RNA studies.

A Survey of Virus Recombination Uncovers Canonical Features of Artificial Chimeras Generated During Deep Sequencing Library Preparation

Chimeric reads can be generated by in vitro recombination during the preparation of high-throughput sequencing libraries. Our attempt to detect biological recombination between the genomes of dengue virus (DENV; +ssRNA genome) and its mosquito host using the Illumina Nextera sequencing library preparation kit revealed that most, if not all, detected host-virus chimeras were artificial. Indeed, these chimeras were not more frequent than with control RNA from another species (a pillbug), which was never in contact with DENV RNA prior to the library preparation. The proportion of chimera types merely reflected those of the three species among sequencing reads. Chimeras were frequently characterized by the presence of 1-20 bp microhomology between recombining fragments. Within-species chimeras mostly involved fragments in opposite orientations and located less than 100 bp from each other in the parental genome. We found similar features in published datasets using two other viruses: Ebola virus (EBOV; -ssRNA genome) and a herpesvirus (dsDNA genome), both produced with the Illumina Nextera protocol. These canonical features suggest that artificial chimeras are generated by intra-molecular template switching of the DNA polymerase during the PCR step of the Nextera protocol. Finally, a published Illumina dataset using the Flock House virus (FHV; +ssRNA genome) generated with a protocol preventing artificial recombination revealed the presence of 1-10 bp microhomology motifs in FHV-FHV chimeras, but very few recombining fragments were in opposite orientations. Our analysis uncovered sequence features characterizing recombination breakpoints in short-read sequencing datasets, which can be helpful to evaluate the presence and extent of artificial recombination.

A multiplex RNA quantification method to determine the absolute amounts of mRNA without reverse transcription

We have developed a highly sensitive microarray-based method that determines the absolute amounts of mRNA in a total RNA sample in a multiplex manner without reverse transcription. This direct mRNA measurement promotes high-throughput testing and reduces bias in transcriptome analyses. Furthermore, quantification of the absolute amount of mRNA allows transcriptome analysis without common controls or additional, complicated normalization. The method, called Photo-DEAN, was validated using chemically synthesized RNAs of known quantities and mouse liver total RNA samples. We found that the absolute amounts of mRNA were successfully measured without the cDNA synthesis step, with a sensitivity of 15 zmol achieved in 7 h.

Gene activity in primary T cells infected with HIV89.6: intron retention and induction of genomic repeats

Background

HIV infection has been reported to alter cellular gene activity, but published studies have commonly assayed transformed cell lines and lab-adapted HIV strains, yielding inconsistent results. Here we carried out a deep RNA-Seq analysis of primary human T cells infected with the low passage HIV isolate HIV_89.6.

Results

Seventeen percent of cellular genes showed altered activity 48 h after infection. In a meta-analysis including four other studies, our data differed from studies of HIV infection in cell lines but showed more parallels with infections of primary cells. We found a global trend toward retention of introns after infection, suggestive of a novel cellular response to infection. HIV_89.6 infection was also associated with activation of several human endogenous retroviruses (HERVs) and retrotransposons, of interest as possible novel antigens that could serve as vaccine targets. The most highly activated group of HERVs was a subset of the ERV-9. Analysis showed that activation was associated with a particular variant of ERV-9 long terminal repeats that contains an indel near the U3-R border. These data also allowed quantification of >70 splice forms of the HIV_89.6 RNA and specified the main types of chimeric HIV_89.6-host RNAs. Comparison to over 100,000 integration site sequences from the same infected cell populations allowed quantification of authentic versus artifactual chimeric reads, showing that 5′ read-in, splicing out of HIV_89.6 from the D4 donor and 3′ read-through were the most common HIV_89.6-host cell chimeric RNA forms.

Conclusions

Analysis of RNA abundance after infection of primary T cells with the low passage HIV_89.6 isolate disclosed multiple novel features of HIV-host interactions, notably intron retention and induction of transcription of retrotransposons and endogenous retroviruses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-015-0205-1) contains supplementary material, which is available to authorized users.

The role of canonical and noncanonical pre-mRNA splicing in plant stress responses

Plants are sessile organisms capable of adapting to various environmental constraints, such as high or low temperatures, drought, soil salinity, or pathogen attack. To survive the unfavorable conditions, plants actively employ pre-mRNA splicing as a mechanism to regulate expression of stress-responsive genes and reprogram intracellular regulatory networks. There is a growing evidence that various stresses strongly affect the frequency and diversity of alternative splicing events in the stress-responsive genes and lead to an increased accumulation of mRNAs containing premature stop codons, which in turn have an impact on plant stress response. A number of studies revealed that some mRNAs involved in plant stress response are spliced counter to the traditional conception of alternative splicing. Such noncanonical mRNA splicing events include trans-splicing, intraexonic deletions, or variations affecting multiple exons and often require short direct repeats to occur. The noncanonical alternative splicing, along with common splicing events, targets the spliced transcripts to degradation through nonsense-mediated mRNA decay or leads to translation of truncated proteins. Investigation of the diversity, biological consequences, and mechanisms of the canonical and noncanonical alternative splicing events will help one to identify those transcripts which are promising for using in genetic engineering and selection of stress-tolerant plants.

Application of metatranscriptomics to soil environments

The activities of soil microbial communities are of critical importance to terrestrial ecosystem functioning. The mechanisms that determine the interactions between soil microorganisms, their environment and neighbouring organisms, however, are poorly understood. Due to advances in sequencing technologies, an increasing number of metagenomics studies are being conducted on samples from diverse environments including soils. This information has not only increased our awareness of the functional potential of soil microbial communities, but also constitutes powerful reference material for soil metatranscriptomics studies. Metatranscriptomics provides a snapshot of transcriptional profiles that correspond to discrete populations within a microbial community at the time of sampling. This information can indicate the potential activities of complex microbial communities and the mechanisms that regulate them. Here we summarise the technical challenges for metatranscriptomics applied to soil environments and discuss approaches for gaining biologically meaningful insight into these datasets.

Characterization of BmKbpp, a multifunctional peptide from the Chinese scorpion Mesobuthus martensii Karsch: gaining insight into a new mechanism for the functional diversification of scorpion venom peptides

BmKbpp is a novel cationic and α-helical peptide from the Chinese scorpion Mesobuthus martensii Karsch, of which function or biological activity has not been characterized so far. Here we showed that BmKbpp possesses strong antimicrobial activity against both Gram-positive and Gram-negative bacteria with a MIC range from 2.3 μM to 68.2 μM for the majority of tested bacteria. BmKbpp also inhibits the growth of tested fungi with an IC50 range from 0.2 μM to 3.1 μM. Because BmKbpp potently inhibits the growth of some antibiotics-resistant pathogens, and shows very weak hemolytic activity, it has considerable potentials for therapeutic applications. Moreover, we found that BmKbpp markedly inhibits the superoxide production in granulocytes or HL-60 cells at the concentrations of submicromolar level; this suggests that BmKbpp can act as a signaling molecule involving innate immune regulation at low concentrations. The C-terminal region of BmKbpp (BmKbpp-C) shows 72% similarity to the peptide K-12, a bradykinin-potentiating peptide. We found that both BmKbpp and BmKbpp-C possess bradykinin-potentiating activity, and the activity of BmKbpp-C is stronger than that of BmKbpp. PCR amplification for the genomic gene of BmBpp showed that it is not a continuous sequence in the genome; it suggests that BmKbpp could come from a recombination event in transcript level. Taken together, our data suggest that multi-functionalization of a single peptide, which is probably mediated by trans-splicing, could be a new mechanism for the functional diversification of scorpion venom peptides.

Absence of the JAZF1/SUZ12 chimeric transcript in the immortalized non-neoplastic endometrial stromal cell line T HESCs

Endometrial stromal sarcomas are rare malignancies, accounting for less than 10% of uterine sarcomas. The most characteristic chromosomal aberration of this tumor type is the translocation t(7;17)(p15-p21;q12-q21) leading to the fusion of two zinc finger genes, JAZF1 and SUZ12. Recently, the presence of the neoplastic JAZF1/SUZ12 fusion transcript was reported in normal cells of human endometrium. One of the positive samples for the JAZF1/SUZ12 transcript was the immortalized T HESCs cell line. This cell line was derived from the stromal cells obtained from an adult female with myomas and immortalized by transfection of a human telomerase gene. Since T HESCs has a normal karyotype and no fusion of the two genes occurs at the genomic level, the JAZF1/SUZ12 transcript was proposed to be generated by regulated trans-splicing between precursor RNAs for JAZF1 and SUZ12. However, no confirmatory reports currently exist. To determine whether the results could be reproduced, the T HESCs cell line was subjected to three different RT-PCR amplifications for the JAZF1/SUZ12 fusion transcript. RT-PCR assays did not amplify JUZF1/SUZ12 cDNA fragments in the T HESCs cell line, whereas the same assays easily generated JUZF1/SUZ12-amplified transcripts in an endometrial stromal cell sarcoma carrying the t(7;17) chromosomal aberration. Thus, the presence, if any, of a JUZF1/SUZ12 chimeric transcript in the immortalized normal T HESCs is not a constant, reproducible result.

Apparent non-canonical trans-splicing is generated by reverse transcriptase in vitro

BACKGROUND

Trans-splicing, the in vivo joining of two independently transcribed RNA molecules, is well characterized in lower eukaryotes, but was long thought absent from metazoans. However, recent bioinformatic analyses of EST sequences suggested widespread trans-splicing in mammals. These apparently spliced transcripts generally lacked canonical splice sites, leading us to question their authenticity. Particularly, the native ability of reverse transcriptase enzymes to template switch during transcription could produce apparently trans-spliced sequences.

PRINCIPAL FINDINGS

Here we report an in vitro system for the analysis of template switching in reverse transcription. Using highly purified RNA substrates, we show the reproducible occurrence of apparent trans-splicing between two RNA molecules. Other reported non-canonical splicing events such as exon shuffling and sense-antisense fusions were also readily detected. The latter caused the production of apparent antisense non-coding RNAs, which are also reported to be abundant in humans.

CONCLUSIONS

We propose that most reported examples of non-canonical splicing in metazoans arise through template switching by reverse transcriptase during cDNA preparation. We further show that the products of template switching can vary between reverse transcriptases, providing a simple diagnostic for identifying many of these experimental artifacts.

Mucroporin, the first cationic host defense peptide from the venom of Lychas mucronatus

The misuse of antibiotics has led our age to a dangerous edge, as antibiotic-resistant pathogens appear to evolve more quickly than antibiotics are invented. Thus, new agents to treat bacterial infection are badly needed. Cationic host defense peptides are on the first line of a host defense system and are thought to be good candidates for treating bacterial infection. Here, a novel cationic host defense peptide, mucroporin, was cloned and characterized from the venom of Lychas mucronatus. The MIC for Staphylococcus aureus was 25 microg/ml, including antibiotic-resistant pathogens. Based on the molecular template of mucroporin, mucroporin-M1 was designed by amino acid substitution. The MIC for S. aureus was 5 microg/ml, including the antibiotic-resistant pathogens methicillin-resistant S. aureus, methicillin-resistant coagulase-negative Staphylococcus, penicillin-resistant S. aureus, and penicillin-resistant S. epidermidis. Moreover, mucroporin-M1 also inhibited gram-negative bacteria. The modes of action of mucroporin and mucroporin-M1 were both rapid killing by disrupting the cell membrane of bacteria, and the number of surviving bacteria was reduced by about 4 to 5 orders of magnitude immediately after peptide delivery. These results showed that mucroporin could be considered a potential anti-infective drug, especially for treating antibiotic-resistant pathogens.

When good transcripts go bad: artifactual RT-PCR 'splicing' and genome analysis

Gene and intron prediction are essential for accurate inferences about genome evolution. Recently, two genome-wide studies searched for recent intron gains in humans, reaching very different conclusions: either of a complete absence of intron gain since early mammalian evolution, or of creation of numerous introns by genomic duplication in repetitive regions. We discuss one possible explanation: the underappreciated phenomenon of "template switching", by which reverse transcriptase may create artifactual splicing-like events in the preparation of cDNA/EST libraries, may cause complications in searches for newly gained introns in repetitive regions. We report large numbers of apparent template switching in transcript sequences from the intron-poor protists Trichomonas vaginalis and Giardia lamblia. Supplementary material for this article can be found on the BioEssays website (http://www.interscience.wiley.com/jpages/0265-9247/suppmat/index.html).

Hsp90n - An accidental product of a fortuitous chromosomal translocation rather than a regular Hsp90 family member of human proteome

Human cells express two isoforms of the Hsp90 protein, called Hsp90alpha and Hsp90beta. Although existence of the third form called Hsp90alphaDeltaN, or Hsp90N was reported in 1998, our investigation, based on the sequence analysis and attempts to reproduce previous results, demonstrate that there is no evidence that Hsp90N gene is present in human genome and no homologs of such a protein are present in other known eukaryotic genomes. We propose that Hsp90N was created as an artifact of a cDNA synthesis or that it is a chimeric protein, being a result of the chromosomal rearrangement that occurred in a single cell line, after this line was established.

The venom gland transcriptome of the Desert Massasauga rattlesnake (Sistrurus catenatus edwardsii): towards an understanding of venom composition among advanced snakes (Superfamily Colubroidea)

Background

Snake venoms are complex mixtures of pharmacologically active proteins and peptides which belong to a small number of superfamilies. Global cataloguing of the venom transcriptome facilitates the identification of new families of toxins as well as helps in understanding the evolution of venom proteomes.

Results

We have constructed a cDNA library of the venom gland of a threatened rattlesnake (a pitviper), Sistrurus catenatus edwardsii (Desert Massasauga), and sequenced 576 ESTs. Our results demonstrate a high abundance of serine proteinase and metalloproteinase transcripts, indicating that the disruption of hemostasis is a principle mechanism of action of the venom. In addition to the transcripts encoding common venom proteins, we detected two varieties of low abundance unique transcripts in the library; these encode for three-finger toxins and a novel toxin possibly generated from the fusion of two genes. We also observed polyadenylated ribosomal RNAs in the venom gland library, an interesting preliminary obsevation of this unusual phenomenon in a reptilian system.

Conclusion

The three-finger toxins are characteristic of most elapid venoms but are rare in viperid venoms. We detected several ESTs encoding this group of toxins in this study. We also observed the presence of a transcript encoding a fused protein of two well-characterized toxins (Kunitz/BPTI and Waprins), and this is the first report of this kind of fusion in a snake toxin transcriptome. We propose that these new venom proteins may have ancillary functions for envenomation. The presence of a fused toxin indicates that in addition to gene duplication and accelerated evolution, exon shuffling or transcriptional splicing may also contribute to generating the diversity of toxins and toxin isoforms observed among snake venoms. The detection of low abundance toxins, as observed in this and other studies, indicates a greater compositional similarity of venoms (though potency will differ) among advanced snakes than has been previously recognized.

Expression pattern analysis of transcribed HERV sequences is complicated by ex vivo recombination

Background

The human genome comprises numerous human endogenous retroviruses (HERVs) that formed millions of years ago in ancestral species. A number of loci of the HERV-K(HML-2) family are evolutionarily much younger. A recent study suggested an infectious HERV-K(HML-2) variant in humans and other primates. Isolating such a variant from human individuals would be a significant finding for human biology.

Results

When investigating expression patterns of specific HML-2 proviruses we encountered HERV-K(HML-2) cDNA sequences without proviral homologues in the human genome, named HERV-KX, that could very well support recently suggested infectious HML-2 variants. However, detailed sequence analysis, using the software RECCO, suggested that HERV-KX sequences were produced by recombination, possibly arising ex vivo, between transcripts from different HML-2 proviral loci.

Conclusion

As RT-PCR probably will be instrumental for isolating an infectious HERV-K(HML-2) variant, generation of "new" HERV-K(HML-2) sequences by ex vivo recombination seems inevitable. Further complicated by an unknown amount of allelic sequence variation in HERV-K(HML-2) proviruses, newly identified HERV-K(HML-2) variants should be interpreted very cautiously.

Molecular dissection of venom from Chinese scorpion Mesobuthus martensii: identification and characterization of four novel disulfide-bridged venom peptides

Scorpion venom is composed of a large repertoire of biologically active polypeptides. However, most of these peptides remain to be identified and characterized. In this paper, we report the identification and characterization of four novel disulfide-bridged venom peptides (named BmKBTx, BmKITx, BmKKx1 and BmKKx2, respectively) from the Chinese scorpion, Mesobuthus martensii (also named Buthus martensii Karsch). BmKBTx is composed of 58 amino acid residues and cross-linked by three disulfide bridges. The sequence of BmKBTx shows some similarities to that of the toxin, birtoxin, and its analogs. It is likely that BmKBTx is a beta-toxin active on Na+ channels, which is toxic to either insects or mammals. BmKITx is composed of 71 amino acid residues with four disulfide bridges. It is the longest venom peptide identified from M. martensii so far. BmKITx shows little sequence identity with scorpion alpha-toxins toxic to insects. It is likely that BmKITx is a new type of Na+ -channel specific toxin active on both insects and mammals. BmKKx1 contains 38 amino acid residues cross-linked by three disulfide bridges and shows 84% sequence identity with BmTx3, an inhibitor of A-type K+ channel and HERG currents. BmKKx1 has been classified as alpha-KTx-15.8. BmKKx2 is composed of 36 residues and stabilized by three disulfide bridges. BmKKx2 is a new member of the gamma-K+ -channel toxin subfamily (classified as gamma-KTx 2.2). The venoms of scorpions thus continue to provide novel toxins with potential novel actions on targets.

Characterization of a novel cDNA encoding a short venom peptide derived from venom gland of scorpion Buthus martensii Karsch: trans-splicing may play an important role in the diversification of scorpion venom peptides

A novel cDNA clone (named BmKT-u) which is a hybrid molecule of the 5'-terminal region of BmKT' cDNA and the 3'-terminal region of an undocumented cDNA (named BmKu), was isolated from a cDNA library made from the venom gland of scorpion Buthus martensii Karsch. BmKT-u codes for a 30 amino acid residue precursor peptide composed of a 20-residue signal sequence, and a putative 10-residue novel mature peptide. Northern blot hybridization showed BmKT-u cDNA is generated from a transcript. RT-PCR experiments excluded the possibility that BmKT-u cDNA is an artifact generated during reverse transcription. Genomic amplifications performed with three pairs of BmKT-u gene-specific primers showed the BmKT-u gene does not exist in the genome of the scorpion as a single transcriptional unit. Genomic cloning for BmKT' showed that the BmKT' gene contains an intron of 509 bp inserted into the region encoding the C-terminal region of the signal peptide. A sequence alignment comparison of the cDNA of BmKT-u with genomic BmKT' revealed that the junction site of the hybrid molecule is located at the 5'-splicing site of the intron. The data suggest that the BmKT-u transcript is a naturally occurring mature mRNA that is generated by trans-splicing. Trans-splicing may contribute to the diversity of venom peptides from venomous animals.

Reverse transcriptase template switching and false alternative transcripts

Reverse transcriptase (RT) can switch from one template to another in a homology-dependent manner. In the study of eukaryotic transcripts, this propensity of RT can produce an artificially deleted cDNA, which can be wrongly interpreted as an alternative transcript. Here, we have investigated the presence of such template-switching artifacts in cDNA databases, by scanning a collection of human splice sites (Information for the Coordinates of Exons, ICE database). We have confirmed several cases at the experimental level. Artifacts represent a significant portion of apparently spliced sequences using noncanonical splice signals but are rare in the context of the whole database. However, care should be taken in the annotation of alternative transcripts, especially when the RT used is poorly thermostable and when the putative intron is flanked by direct repeats, which are the substrate for template switching.

Intergenic transcripts in genes with phase I introns

Nonsense-mediated mRNA decay (NMD) is a eukaryotic quality-control mechanism that detects and degrades aberrant transcripts prematurely terminating translation. NMD may be elicited by intergenic transcripts that contain premature termination codons (PTCs), but chimeric mRNAs of genes that have introns of identical phase would be predicted to lack PTCs and escape NMD. We examined intron phase I-containing HLA class II genes for the presence of intergenic mRNAs and found an extraordinary diversity of correctly spliced and polyadenylated intergenic transcripts. They lacked a significant homology at the chimeric joins and had no PTCs. Their expression levels were very low and positively correlated with the expression of natural transcripts. In contrast, pair-wise mixtures of separately transcribed plasmids carrying full-length HLA-DQB1, -DQA1, -DRB1, and -DRA cDNAs produced only hybrid molecules that lacked canonical exon boundaries, had homologous chimeric joins, and occasionally contained PTCs, implicating in vitro artifacts generated by template switching of Taq polymerase and reverse transcriptase. A differential exon structure of hybrid molecules observed in vitro and in cellular RNA preparations suggests that intergenic mRNAs with canonical exon boundaries arise in vivo during exon joining and/or transcription. Since the observed intergenic mRNAs may encode mixed class II heterodimers that were previously shown to present antigens it will be interesting to determine functional properties of such molecules in future studies.

Detection of a novel sense-antisense RNA-hybrid structure by RACE experiments on endogenous troponin I antisense RNA

Conformational changes in the troponin/tropomyosin complex significantly alter the mechanical properties of cardiac muscle. Phosphorylation of cardiac troponin I, part of the troponin/tropomyosin complex, reduces calcium affinity, which leads to increased relaxation of cardiac muscle. Because cardiac troponin I plays a central role in tuning the heart to different work demands, detailed knowledge of troponin I protein regulation is required. Our group previously detected naturally occurring antisense RNA for troponin I in human and rat hearts, and here, attempt to unravel the structure of rat cardiac troponin I antisense RNA. We performed rapid amplification of cDNA ends (RACE) experiments and discovered antisense sequences identical to a copy of the sense mRNA, which led us to conclude that the antisense RNA must be transcribed from troponin I mRNA in the cytoplasm. Moreover, we isolated RNA structures comprising sense and antisense sequences in one continuous molecule. As we found no homolog structures described in the literature, we called this "hybrid RNA." Because a duplex formation was demonstrated previously we concluded that hybrid RNA is a consequence of a tight interaction between sense and antisense troponin I RNA in vivo, which we discuss in the article.

Identification and functional characterization of novel scorpion venom peptides with no disulfide bridge from Buthus martensii Karsch

The scorpion venom peptides with no disulfide bridge are rarely identified and poorly characterized so far. Here, we report the identification and characterization of four novel disulfide-bridge-free venom peptides (BmKa1, BmKa2, BmKb1 and BmKn2) from Buthus martensii Kasch. BmKa1 and BmKa2 are very acidic and hydrophilic, showing no any similarity to other proteins, whereas BmKb1 and BmKn2 both are basic, alpha-helical peptide with an amidated C-terminus, showing a little homology with other peptides. Functional tests with synthetic peptide showed that BmKn2 has strong antimicrobial activity against both Gram-positive and Gram-negative bacteria, whereas BmKb1 has weak activity in inhibiting the growth of these bacteria.