Short hairpin RNA-mediated knockdown of protein expression in Entamoeba histolytica
© Linford et al; licensee BioMed Central Ltd. 2009
Received: 20 September 2008
Accepted: 17 February 2009
Published: 17 February 2009
Entamoeba histolytica is an intestinal protozoan parasite of humans. The genome has been sequenced, but the study of individual gene products has been hampered by the lack of the ability to generate gene knockouts. We chose to test the use of RNA interference to knock down gene expression in Entamoeba histolytica.
An episomal vector-based system, using the E. histolytica U6 promoter to drive expression of 29-basepair short hairpin RNAs, was developed to target protein-encoding genes in E. histolytica. The short hairpin RNAs successfully knocked down protein levels of all three unrelated genes tested with this system: Igl, the intermediate subunit of the galactose- and N-acetyl-D-galactosamine-inhibitable lectin; the transcription factor URE3-BP; and the membrane binding protein EhC2A. Igl levels were reduced by 72%, URE3-BP by 89%, and EhC2A by 97%.
Use of the U6 promoter to drive expression of 29-basepair short hairpin RNAs is effective at knocking down protein expression for unrelated genes in Entamoeba histolytica, providing a useful tool for the study of this parasite.
The human parasite Entamoeba histolytica (E. histolytica) is a unicellular protozoal parasite that phylogenetically is placed on one of the lowermost branches of the eukaryotic tree, closest to Dictyostelium discoideum . It is an unusual organism, having 9,938 predicted genes, with slightly less than one third (31.8%) of its predicted proteins having no homologues in GenBank . Humans are its only natural hosts, and E. histolytica is spread by ingestion of contaminated food or water via the fecal-oral route and thus tends to endemically infect people under circumstances where hygiene is poor . It has a simple life cycle, alternating between infective quadrinucleate cysts and invasive motile trophozoites . 80% of people infected with E. histolytica are colonized asymptomatically; in the remaining 20%, trophozoites invade into the intestinal epithelium, resulting in clinical disease . It is estimated that there are 50 million symptomatic cases of amebic colitis and 100,000 deaths per year worldwide due to E. histolytica .
The discovery that double-stranded RNA (dsRNA) can initiate post-transcriptional sequence-specific gene silencing of cellular genes  via translational repression or degradation of mRNA in most eukaryotic cells has become an important tool in assessing and manipulating gene function. This mechanism of RNA interference (RNAi) may have evolved as a defense against viruses and transposable elements with dsRNA intermediates [6, 7]. The small RNA intermediates in this process, short interfering RNAs (siRNAs), result from dsRNA being cleaved at 21- to 23- nucleotide intervals  by an RNase III-type protein, Dicer , and are then incorporated into the RNA-induced silencing complex (RISC), which includes Argonaute "Slicer" protein [8, 10]. The antisense strand of the siRNA is used to guide the RISC to its target mRNA, which is then cleaved by Argonaute [11, 12]. RNAi effects can be amplified by the action of RNA-dependent RNA polymerases (RdRPs). siRNAs act as primers for RdRPs, which form new dsRNAs using the target mRNA as a template, which are subsequently cleaved into siRNAs with sequences corresponding to target mRNAs but differing from the original dsRNAs [13, 14]. Genes encoding RdRPs have been identified in many organisms, but not in flies or mammals .
E. histolytica possesses the molecular machinery for RNAi. It has a gene [GenBank:XM_645408] [2, 15, 16] encoding a protein which has a single RNase III domain and possesses RNase III activity, and could perform the Dicer role as a dimer. It also has two Argonaute homologs [GenBank:XM_651344, XM_651422] [2, 15–17] and an RdRP [GenBank:XM_646217] [2, 15]. Exploitation of RNAi for knockdown of gene expression is an attractive approach for E. histolytica, as there is no evidence for meiotic division or detectable homologous recombination of genes [18–20], thus it has not been possible to generate gene knockouts [18, 21]. Multiple copies of the genome, and even nuclei, occur in the parasite due to an apparent lack of the normal cell cycle regulatory checkpoints [22, 23]. Homologous Entamoeba cell-cycle regulation genes are divergent from typical eukaryotic versions and may not have equivalent function . This presents difficulties in studying gene function or in isolating recessive mutations . The study of the function of individual genes in the past has been limited to other techniques, such as the over-expression of wild-type or mutant genes, and other methods of gene inactivation such as antisense [21, 24]. Methods of RNAi used in E. histolytica have included the use of long dsRNA expressed by an E. histolytica RNA polymerase II promoter, which was successfully used to knock down expression of the E. histolytica proteins Diaphanous, Klp5 and EhSTIRP [18, 25, 26], and the soaking of trophozoites in artificial siRNAs to knock down γ-tubulin expression . These reports of RNAi use in E. histolytica showed knockdown of a single gene or of a gene family. Here, we report in this study the success of the method of expression of short hairpin RNAs driven by the E. histolytica U6 promoter to knock down protein expression in E. histolytica of three unrelated genes.
Short hairpin RNAs (shRNAs) have a similar structure to siRNAs except the sense and antisense strands are connected at one end by a short loop, and function like siRNAs to knock down gene expression . shRNAs can be produced from an expression vector as a single transcript from a RNA polymerase III promoter. The eukaryotic U6 promoter offers two advantages over other RNA polymerase III promoters: the promoter region immediately upstream of the transcribed sequence for the U6 small nuclear RNA gene includes all the required regulatory elements [28, 29], and the termination sequence consists of 4 to 5 thymidine residues rather than a poly-A tail [28, 29]. A variety of shRNA loop and stem lengths have been tested, with the loop UUCAAGAGA  used in a number of mammalian shRNA constructs, including Gou et al (2003) , and is also used in the constructs in this study. Longer hairpins with 29-base pair stems appear to be better inhibitors of gene expression than ones with shorter 19–21 bp stems . Increased effectiveness has also been seen for similarly-sized longer artificial siRNAs, with only one siRNA apparently generated per longer shRNA or siRNA [31, 32].
Genes selected for knockdown: The three genes selected for knockdown in this study, Igl, URE3-BP, and EhC2A, are genes involved in amebic virulence under study in our laboratory; they were selected since we wanted to create an additional tool for studying the function and role of these genes in amebic virulence.
Igl, the intermediate subunit of the galactose- and N-acetyl-D-galactosamine- (Gal/GalNAc) inhibitable lectin [33, 34], is a 150 kDa protein. The Gal/GalNAc lectin, the major defined amebic adhesin, is a virulence factor mediating adherence to target cells in the first step of contact-dependent cell killing . The lectin binds to terminal galactose or GalNAc residues in glycoproteins on the surfaces of target cells, and is composed of three subunits: the heavy subunit Hgl (containing a carbohydrate-recognition domain), the intermediate subunit Igl, and the light subunit Lgl . The integral-membrane Hgl is disulfide-bonded to the GPI (glycosylphosphatidylinositol)-anchored Lgl. Igl is also GPI-anchored to the membrane . Evidence that Igl is associated non-covalently with the Hgl-Lgl heterodimer includes that Igl and the Hgl-Lgl heterodimer co-migrate in native gel electrophoresis, and affinity-purification of Igl with monoclonal antibodies results in the co-purification of the Hgl-Lgl heterodimer [3, 33, 34]. Igl is encoded by two unlinked gene copies, Igl1 [GenBank:AF337950]  and Igl2 [GenBank:XM_647302] ; [GenBank:AF337951] , producing ~1100 aa proteins that are 81% identical and contain 32 CXXC repeats. CXXC repeats are also found in a family of transmembrane kinases of E. histolytica and the Giardia lamblia variant-specific surface proteins .
URE3-BP, Upstream Regulatory Element 3-Binding Protein [GenBank:AF291721] , is a 22.6 kDa calcium-regulated transcription factor encoding two EF-hand motifs, which are associated with calcium-binding activity . URE3-BP binds to the URE3 (Upstream Regulatory Element 3) consensus motif, TATTCTATT, found in the promoter of hgl5, which is one of the genes encoding the Gal/GalNAc lectin heavy subunit, and is also present in the ferredoxin 1 (fdx1) promoter, thereby regulating the expression of these genes . The human neuronal protein DREAM (calsenilin) is the only other known example of a calcium-responsive transcription factor with EF hands .
EhC2A [GenBank:XM_650207]  is a 22 kDa calcium-binding membrane protein containing a conserved C2 domain, is associated with the ability to bind phospholipids, and has a proline-rich C-terminal tail. This protein was found to be associated to the amebic phagosome . A C2 domain, identified originally in protein kinase C, is a Ca2+-binding motif that allows calcium-dependent protein anchoring to or interaction with membranes; these domains are found in a number of signaling proteins in eukaryotes .
A gene for which we have previously shown knockdown is PATMK, Phagosome-Associated Transmembrane Kinase 96 [GenBank:XM_650501] [2, 39]. PATMK is a transmembrane kinase family member found in the early phagosome and is involved in the phagocytosis of human erythrocytes . It contains an intracellular putative kinase domain, a short membrane-spanning region, and an ectodomain containing CXXC-repeats like Igl [35, 39].
We report here the effectiveness of shRNAs in silencing genes in Entamoeba histolytica. Expression of 29-bp shRNAs driven by the E. histolytica U6 promoter was successful in knocking down protein expression of the three different and unrelated genes in E. histolytica reported in this study, and we previously showed knockdown for a fourth gene . This method of gene knockdown appears to function well for a variety of gene types, and should be a useful tool for studying gene function in this organism.
Construction of shRNA constructs
Sequences chosen to generate shRNA constructs that were successfully transfected into amebae
Location in mRNA/cDNA (bp from ATG)
Total length of target mRNA (bp)
3306 (Igl1), 3318 (Igl2)
502–530; also 406–434
EhC2A (363–391 scrambled)
Does not match any E. histolytica mRNA
Does not match any E. histolytica mRNA
Oligos used for generating shRNA constructs by PCR and transfected into amebae
U6 HindIII forward
Igl1 (272–300) R1
Igl1 (272–300) R2
Igl (1198–1226) R1
Igl (1198–1226) R2
Igl (2412–2440) R1
Igl (2412–2440) R2
Igl (2777–2805) R1
Igl (2777–2805) R2
URE3-BP (350–378) R1
URE3-BP (350–378) R2
URE3-BP (580–608) R1
URE3-BP (580–608) R2
EhC2A (363–391) R1
EhC2A (363–391) R2
EhC2A (502–530) R1
EhC2A (502–530) R2
EhC2A (363–391 scrambled) R1
EhC2A (363–391 scrambled) R2
Sequences of oligos used for amplification in qRT-PCR
mRNA/cDNA section amplified (bp from ATG)
Total length of mRNA (bp)
Igl 5' F
85–450 (Igl1), 85–459 (Igl2)
3306 (Igl1), 3318 (Igl2)
Igl 5' R
Igl 3' F
2967–3166 (Igl1), 2979–3178 (Igl2)
Igl 3' R
URE3-BP 5' F
URE3-BP 5' R
URE3-BP Middle F
URE3-BP Middle R
URE3-BP 3' F
URE3-BP 3' R
variable (detects all family members/alleles)
variable (detects all family members/alleles)
Transfectants were maintained at 15 μg/ml hygromycin. For knockdown studies, the hygromycin concentration was increased every 24 hours until the final level of selection was achieved, and was maintained for 48 hours, in order to increase the copy number of the episomal shRNA vector [41–43]. The level of hygromycin selection was increased until the desired knockdown was attained, up to 100 μg/ml. Transfected trophozoites selected with 100 μg/ml hygromycin continued to grow and divide for at least two weeks under continuous selection. A shRNA directed against green fluorescent protein (GFP) , with a sequence matching nothing in the E. histolytica genome, was utilized as a control for transfection and hygromycin selection for the Igl and URE3-BP transfectants. GFP shRNA transfectants were selected with the same level of hygromycin as other shRNA transfectants. For EhC2A, a scrambled control matching nothing in the E. histolytica genome was created, containing the same nucleotides as the EhC2A (363–391) shRNA, but in a different order. Sequences of the shRNA sense strands are shown in Table 1. Non-transfected HM1-IMSS amebae were also included, with the results for Western blotting and qRT-PCR being statistically the same as the GFP controls. Three biological replicates were grown per shRNA transfectant, and one for the nontransfected HM1:IMSS amebae. All sample trophozoites were grown in 25 cm2 tissue culture flasks, and were harvested for crude lysate and for RNA isolation on the same day from the same flask. For protein and mRNA comparison, actin was used as the "housekeeping" control gene, as a loading and normalization control.
Knockdown of Igl protein
Summary of Igl1 protein levels in Igl shRNA transfectants
shRNA Transfectant or Control Sample
% of Igl1 protein level (± SE)
100.0 ± 3.6
115.5 ± 11.8
95.3 ± 3.2
27.8 ± 1.3
42.3 ± 2.1
38.1 ± 3.1
Knockdown of Igl mRNA
Summary of Igl mRNA levels in Igl shRNA transfectants
shRNA transfectant or control sample
Igl 5' oligo pair
Igl 3' oligo pair
Igl1 oligo pair
Igl2 oligo pair
100.0 ± 4.1
100.0 ± 4.9
100.0 ± 3.0
100.0 ± 4.0
101.4 ± 4.3
96.1 ± 3.5
105.5 ± 3.1
103.9 ± 6.1
100.6 ± 5.0
103.4 ± 9.1
91.1 ± 6.9
106.0 ± 5.2
71.3 ± 2.9
67.1 ± 3.0
61.1 ± 3.2
70.2 ± 2.7
70.9 ± 2.7
62.1 ± 1.6
68.3 ± 2.5
76.8 ± 1.6
68.1 ± 3.3
62.3 ± 2.9
74.1 ± 3.3
77.8 ± 3.0
Knockdown of URE3-BP protein
Summary of URE3-BP protein levels in URE3-BP shRNA transfectants
shRNA transfectant or control sample
% of control protein level (± SE)
100 ± 9.9
111.3 ± 15.8
10.8 ± 1.0
13.8 ± 2.6
Knockdown of URE3-BP mRNA
Summary of mRNA levels in GFP shRNA control transfectants, URE3-BP shRNA transfectants, and nontransfected HM1:IMSS trophozoites
shRNA transfectant or control sample
URE3-BP 5' oligo pair
URE3-BP middle oligo pair
URE3-BP 3' oligo pair
100.0 ± 2.9
100 ± 2.8
100 ± 4.3
106.4 ± 5.8
108.9 ± 5.6
102.8 ± 5.0
67.0 ± 2.5
67.4 ± 2.0
72.2 ± 2.8
12.4 ± 0.8
13.5 ± 3.3
12.5 ± 3.8
Knockdown of EhC2A protein
Northern blots of small RNAs
We have utilized the U6 promoter to drive expression of shRNAs with a 29-bp stem and a 9-nt loop to knock down protein expression of three unrelated genes: a membrane protein, Igl, the intermediate subunit of the Gal/GalNAc lectin; URE3-BP, a calcium-regulated transcription factor, upstream regulatory element 3- binding protein; and EhC2A, a membrane-binding protein. Previously we had reported preliminary experience with this system in the near-complete knockdown of phagosome-associated transmembrane kinase 96 (PATMK) . In the work reported here, the highest level of protein knockdown for Igl was 72%, for URE3-BP 89%, and for EhC2A 97%. We concluded that this was a reliable and effective system for gene knockdown in E. histolytica. This method has advantages over other methods used for gene silencing: the U6-shRNA expression cassettes are small (420 bp), appear to be active against different types of genes, yield significant knockdown, and the expression vector, once transfected, allows continuous expression of shRNAs, thus avoiding performing multiple transfections, and the shRNAs can be easily synthesized via PCR.
Not every transfected shRNA construct was equally effective in silencing gene expression. For example, neither the EhC2A (502–530) nor the Igl (2412–2440) shRNA construct blocked gene expression. In the case of Igl (2412–2440), the run of four thymidines at positions 19–23 in the shRNA sense strand could possibly cause RNA polymerase III to terminate the transcript prematurely. In the case of the EhC2A (502–530) shRNA construct, the shRNA could bind to two locations in the mRNA, both the originally targeted region (502–530), but also about 100 nt upstream (406–434), and perhaps this had an inhibitory effect on the ability of this shRNA to allow knockdown.
Factors other than the shRNA sequence affect the ability of a shRNA to down-regulate gene expression. The secondary structure of the transcript affects the ability of the RISC to bind to its target site [44, 45], and the relative abundance and stability of an mRNA may play a significant role in determining whether a given shRNA will effectively lead to the degradation of its target message. In addition, the stability of a protein product may also be a determinant in the detection of a knockdown phenotype. The protein with the least knockdown in these studies, Igl, was the most abundant; EhC2A was the least abundant and had the most knockdown . The level of hygromycin utilized to select for transfectants was an important determinant of the extent of protein knockdown. Igl knockdown was twice as effective with 100 μg/ml as with 30 μg/ml of hygromycin selection.
The qRT-PCR data was not correlated directly with the level of protein knockdown. For the Igl transfectants, the mRNA knockdown level was not as high as the protein knockdown level, indicating the possibility that the protein could have a high turnover rate or be somewhat unstable. For URE3-BP, the URE3-BP (350–378) and (580–608) transfectants had similar levels of protein knockdown; however, the mRNA levels in the URE3-BP (350–378) transfectants were higher (67% of the control level), versus the URE3-BP (580–608) transfectants (13.5% of the control level). This difference is probably not due to partial mRNA decay, since the qRT-PCR data showed consistent URE3-BP levels among the three oligo pairs amplifying the 5', middle, and 3' sections of the transcript. One possible explanation could be that the secondary structure of the URE3-BP mRNA at the location of the URE3-BP (350–378) shRNA could interfere sufficiently with the RISC being able to cleave the mRNA but still allow RISC binding, allowing for a degree of translational inhibition in addition to some mRNA destruction.
The E. histolytica U6 promoter appears to be functional and producing shRNAs: the Northern blots of the small RNAs detected two sizes of small RNAs when probed with oligos that were complementary to the individual sense and antisense strands of the shRNAs. These may represent the unprocessed hairpin and the resulting siRNAs after Dicer processing. Surprisingly, the abundance of the small RNA was not proportional to the level of silencing. Northern blots may not be sensitive enough to identify low-level small RNA production, with low-level production adequate for protein knockdown.
We report the knockdown of three genes in this study: Igl, the intermediate subunit of the Gal/GalNAc lectin; the calcium-responsive transcription factor URE3-BP; the membrane-binding protein EhC2A, by transfecting E. histolytica with expression vectors using the E. histolytica U6 promoter to drive expression of shRNAs targeting endogenous genes. We have previously reported the knockdown of transmembrane kinase PATMK . These genes come from different families, with different functions, so this shRNA knockdown method appears robust and not specific to only one gene or gene family.
Culture of trophozoites
E. histolytica strain HM1:IMSS trophozoites were grown axenically in TYI-S-33 (Trypticase-yeast extract-iron-serum) (TYI) medium supplemented with 1× Diamond's vitamins (SAFC Biosciences, Lenexa, KS, USA), 15% heat-inactivated adult bovine serum (Gemini Bio-Products, West Sacramento, CA), 100 U of penicillin/ml and 100 μg streptomycin sulfate/ml (Gibco/Invitrogen, Carlsbad, CA, USA), at 37°C in 25 cm2 tissue culture flasks  in a volume of 50 ml, and then transfected as described below.
Transfection of amebae
Plasmid DNA was prepared using the HiSpeed Qiagen Maxi Kit (Qiagen, Valencia, CA, USA). Medium 199 (M199) (Gibco BRL/Invitrogen, Carlsbad, CA, USA) was supplemented with 5.7 mM cysteine, 25 mM HEPES, and 0.6 mM ascorbic acid , adjusted to pH 7.0 and filter-sterilized. Twenty μg plasmid DNA diluted in 100 μl supplemented M199s medium (M199S) in 2-ml microcentrifuge tubes was mixed with 15 μl of SuperFect or Attractene transfection reagent (Qiagen, Valencia, CA, USA), and incubated at room temperature to allow transfection-complex formation as per the manufacturer's instructions. Heat-inactivated bovine serum was added to the remaining M199S to a 15% concentration. Amebae were harvested by tapping the tissue culture flasks on a benchtop, were centrifuged at 200 × g for 5 min at 4°C, and suspended in M199S with serum to 2.5 × 105 amebae/ml. Tubes containing transfection complexes were filled with the suspended trophozoites, the contents mixed by inversion, and the tubes were incubated horizontally for 3 hours at 37°C. Tube contents were added to warm TYI in 25 cm2 tissue culture flasks, and incubated overnight at 37°C. 15 μg/ml hygromycin (Invitrogen, Carlsbad, CA, USA) was added for selection after the overnight incubation . After 4–5 days, 25 ml of the TYI was removed to a new 25 cm2 tissue culture flask, and 25 ml fresh TYI with hygromycin was added to each of the flasks. Transfectants were usually apparent 1–2 weeks after transfection.
E. histolyticashRNA constructs
All short hairpin RNAs used in this study were expressed by the U6 promoter [GenBank:U43841]  (Figure 1A) and cloned into the amebic expression vector pGIR310, a modification of pGIR308 [49, 50] by the addition of a short polylinker containing HindIII, SalI, and NotI restriction sites (Figure 1B). Modified pGIR310 conferred resistance to hygromycin in E. histolytica and to ampicillin in Escherichia coli (E. coli). All shRNA constructs used in these studies had the same structure: a short hairpin consisting of a 29-nucleotide sense strand, followed by the 9-nucleotide loop and the 29-nucleotide complementary antisense strand (Figure 1).
Sequence selection for shRNA constructs
The Ambion siRNA finder  was used to select possible siRNA sequences of 21 mers beginning with two adenine residues. To select sequences that would target Igl1 and Igl2 both separately and simultaneously, those portions of their coding sequences which were identical or divergent were input separately, while the entire coding sequence of URE3-BP was used to select siRNA sequences. For EhC2A the portion of the gene sequence selected for targeting was the poly-proline region (bases 301–567) since this region is least similar to the other gene family members. From the pool of selected 21 mer sequences, those with runs of more than 4 As or Ts were eliminated, and those with GC content between 30% and 50% were lengthened to 29 bp by adding the next eight bases in the genomic sequence. The TIGR E. histolytica Genome Project database  was used to check that each 29-bp sequence was unique to its gene, with non-unique ones eliminated. A minimum of four unique sequences were selected per gene. To create a scrambled control sequence, one of the selected sequences was chosen, and the bases were scrambled (each began with the AA dinucleotide); these sequences were then checked to confirm they matched nothing in the E. histolytica genome. In addition, a sequence targeted to the green fluorescent protein (GFP) was included as a control . The chosen sequences, those ultimately transfected into E. histolytica HM1:IMSS trophozoites, are shown in Table 1. Constructs that did not successfully transfect are not shown.
shRNA primer design
Primers were designed based on the method used by Gou et al (2003)  to yield PCR-generated shRNA constructs in a 2-step PCR process diagrammed in Figure 1. The final PCR product contained the E. histolytica U6 promoter followed by the sense strand of the hairpin, the 9 bp loop (TTCAAGAGA) , the antisense strand of the hairpin, and the U6 terminator sequence . An ApaI restriction site (GGGCCC) was included between the 3' end of the U6 promoter and the beginning of the shRNA sequence . To facilitate cloning of the PCR product into the expression vector, a HindIII site was added to the 5' end of the U6 promoter sequence, and a NotI site was added following the terminator sequence. The selected siRNA sequences, shown in Table 1, were used to design oligos to create shRNAs. Two rounds of PCR were employed to generate the final shRNA constructs, using one forward primer and two reverse primers, whose sequences are listed in Table 2. In the first round of PCR, the E. histolytica U6 promoter followed by the sense strand and the loop were generated using a forward primer amplifying the 5' end of the U6 promoter and a first reverse primer containing the sequence of the sense strand of the shRNA and the future loop (Figure 1A, Table 2). A second round of PCR created the completed shRNA construct: the product from the first round was used as a template, using the same forward primer as in the first round, and a second reverse primer containing the sequence of the loop, the antisense strand sequence, and the U6 termination sequence (Figure 1A, Table 2). The control GFP sequence  was used to design oligos for making a shRNA control construct. Sense strand sequences chosen to make the Igl, URE3-BP and EhC2A shRNA constructs successfully transfected into trophozoites are shown in Table 1, and PCR oligos used to amplify these sequences to generate shRNAs via PCR are shown in Table 2.
PCR conditions for generating shRNAs
Initially, E. histolytica genomic DNA was used as a template for the first round of Igl shRNA PCRs. For the URE3-BP and EhC2A shRNA PCRs, the cloned U6 promoter was used as the PCR template: the Igl shRNA plasmids were digested with HindIII and ApaI and the U6 promoter was gel-purified using the QIAquick Gel Extraction Kit (Qiagen, Valencia, CA, USA). Two rounds of PCR were used to generate the shRNA constructs.
The first PCR round generated the sense strand of the hairpin and the loop. Reaction volumes of 40 μl were set up, each consisting of 0.6 μl SAHARA™ DNA polymerase (Bioline USA Inc., Taunton, MA, USA), 4 μl 10× SAHARA™ PCR buffer, 3.2 μl 50 mM MgCl2, 2 μl dNTP mix (stock 10 mM each), 0.4 μl U6 HindIII forward oligo (100 μM stock), 0.4 μl R1 oligo (100 μM stock), 1 μl (200 ng E. histolytica genomic DNA or 25 ng gel-purified digest of HindIII/ApaI U6 promoter), and 28.4 μl sterile water. Cycling conditions were as follows: 95°C for 8 minutes, 10 cycles of 95°C 45 sec, 40°C 1 min, 68°C 1 min 30 sec; 25 cycles of 95°C 45 seconds, 52°C 1 min, 68°C 1 min 30 sec, and a 5 min final extension at 68°C. 5 μl of each PCR product was subjected to agarose gel electrophoresis to check that the products were ~380 bp.
In the second PCR round, the first round PCR product was used as a template to add the antisense strand of the hairpin, the terminator sequence and the NotI site. Each 100 μl-volume reaction contained 2 μl SAHARA™ DNA Polymerase (Bioline USA Inc., Taunton, MA, USA), 10 μl 10× SAHARA™ PCR buffer, 8 μl 50 mM MgCl2, 5 μl dNTP mix (10 mM each), 0.8 μl U6 HindIII forward oligo (100 μM), 0.8 μl R2 oligo (100 μM), 2 μl PCR product from the first PCR round, and 71.4 μl sterile water. Cycling conditions were: 95°C for 8 minutes, 10 cycles of 95°C 45 sec, 18.5°C 1 min 30 sec, 68°C 1 min 30 sec; 30 cycles of 95°C 45 seconds, 55°C 1 min, 68°C 1 min 30 sec, and a 5 min final extension at 68°C. The low annealing temperature in the early cycles of the second PCR was used since the loop is the only overlap between the first round product and the second round reverse oligo. The second round PCR products were checked by agarose gel electrophoresis for products of the correct size (~420 bp). Sometimes a smaller product was present in addition to the correct size product in the final PCR product; this was ignored since it had no effect on the subsequent cloning steps. These final PCR products were ethanol-precipitated, then they and modified pGIR310 were digested with HindIII and NotI. The digested expression vector was gel-purified using the QIAquick Gel Extraction Kit (Qiagen, Valencia, CA, USA). The PCR products were ethanol-precipitated at -20°C overnight, resuspended, ligated into modified pGIR310, transformed into E. coli, and colonies screened. Plasmids with inserts were sequenced, and those with perfect U6 promoter and hairpin sequences were cultured, plasmids were isolated using the Qiagen HiSpeed Maxiprep Kit (Qiagen, Valencia, CA, USA), and transformed into HM1:IMSS strain trophozoites as described above.
The Igl, URE3-BP, or EhC2A shRNA transfectants were grown in 25 cm2 tissue culture flasks and selected beginning with 15 μg/ml of hygromycin, with the hygromycin level increased every 24 hours until the final level of selection was reached, and this level was maintained for 48 hours before harvesting. The GFP control, all three Igl, and the URE3-BP (350–378) transfectants were selected with 100 μg/ml, the URE3-BP (580–608) shRNA transfectants with 75 μg/ml, and the EhC2A samples with 90 μg/ml hygromycin. The final concentration of hygromycin selection differs since the selection was increased until the desired level of knockdown was achieved. There were three biological replicates per shRNA transfectant, and one for the HM1:IMSS nontransfected trophozoites. Trophozoites were harvested as described above for transfection, counted, resuspended in ice cold Lysis Buffer (150 mM NaCl, 50 mM Tris, 5× Sigma protease inhibitor cocktail (P2714) (Sigma-Aldrich, St. Louis, MO, USA), 25 μg/ml E-64 (Sigma-Aldrich, St. Louis, MO, USA)) at an initial concentration of 2 × 106–5 × 106 amebae/ml, and lysed by sonication by pulsing twice for 10 seconds each with a 10 second rest on ice between pulses. Protein was quantified and sample lysates were diluted to the same protein concentration, were serially-diluted 1:2, 1:4, and 1:8 with Lysis Buffer, and were subjected to SDS-PAGE on 12% (Igl) or 15% (URE-BP and EhC2A) gels. All sample lysates and dilutions were done in triplicate (technical replicates). Gels were transferred to PVDF membrane, membranes were cut in half so each half could be probed separately, were blocked in 5% milk, and incubated with either antibodies against Igl1, URE3-BP, EhC2A, or control antibodies against actin (anti-actin from Santa Cruz Biotechnology (Santa Cruz Biotechnology, Santa Cruz, CA) or Sigma (Sigma-Aldrich, St. Louis, MO, USA)).
The ECL kit from Roche (Roche Applied Science, Indianapolis, IN, USA) was used to treat membranes after secondary antibody incubation, bands were visualized on film, film images were electronically scanned, and Scion Image Beta 4.0.3 software (Scion Corporation, Frederick, MD, USA) was used to quantify band intensity. The background value was subtracted from each band value, the ratio of Igl, EhC2A, or URE3-BP band value to the control actin band value was taken for normalization, and then the sample shRNA transfectant lines were compared to the control GFP shRNA or scrambled shRNA transfectant line. HM1:IMSS nontransfected samples were also included. Values for each shRNA transfectant were averaged, and the SE for each average was calculated using the total number of biological replicates multiplied by the number of technical replicates. Statistical analysis was performed using Student's t test (two-tailed) or ANOVA. The GraphPad QuickCalcs P-value calculator was used to calculate the P-values .
Isolation of total RNA
Igl, URE3-BP, and control GFP transfectant shRNA lines were selected with hygromycin as described above for Western blotting, and samples were collected and frozen in TRIzol reagent (Invitrogen, Carlsbad, CA, USA) at -80°C for RNA isolation at the same time as those harvested for crude lysate for protein analysis. Total RNA isolated from each shRNA transfectant and nontransfected HM1:IMSS sample using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) was treated with RNase-free recombinant DNase I (Roche, Indianapolis, IN, USA) for 30 minutes at 37°C, and purified on RNeasy columns using the RNeasy Mini kit as per the manufacturer's instructions (Qiagen, Valencia, CA, USA). Five μg RNA per sample was reverse-transcribed using SuperScriptII (Invitrogen, Carlsbad, CA, USA) and anchored oligo dT, including samples with no reverse transcriptase added (no-RT controls). To check samples for residual DNA contamination in the no-RT controls, each was screened with primers specific for the Jacob cyst-specific gene . If residual DNA contamination was observed, the RNA was treated again with DNase I as above, re-purified on RNeasy columns, and re-screened.
Quantitative reverse-transcription real-time PCR (qRT-PCR)
After the screen for residual DNA contamination was completed, the cDNA was quantified, and sample cDNAs were diluted to 100 ng/μl. HM1:IMSS cDNA was also serially-diluted for making a standard curve. All primers used for qRT-PCR in this study were selected to amplify <400 bp sections of mRNA. Amplification of actin  was performed for use as a normalization control. Oligo sequences used in qRT-PCR are shown in Table 3. Each oligo pair was checked using the E. histolytica genomic database  to validate that only the gene intended would be amplified, except for actin and Jacob, which were designed to detect all family members . An MJ Research Opticon2 DNA Engine (Bio-Rad, Hercules, CA, USA) was utilized for all qRT-PCR runs. ~200 ng of each sample or control cDNA, or serially-diluted HM1:IMSS cDNA for standard curves, was added to each sample well in a 96-well plate for each set of amplifications. cDNA from each biological replicate was run in quadruplicate (technical replicates), and there were three biological replicates per transfectant line, except for HM1:IMSS nontransfected samples, which had one biological replicate. No-RT controls were also included for each set of samples. Each well contained in addition to the cDNA: 1.25 U HotStarTaq (Qiagen, Valencia, CA, USA), 1× HotStarTaq PCR Buffer, 0.5 mM MgCl2, 200 μM each dNTP, 1 μM each oligo, 1:10,000 dilution of Sybrgreen (Invitrogen, Carlsbad, CA, USA), 0.15% Triton X-100, and water to a volume to 25 μl per well.
qRT-PCR cycling conditions were 95°C for 15 minutes, followed by 40 cycles of 95°C 30 s; 54°C 30 s; 72°C 45 s, followed by one cycle of 72°C for 3 min. At the end of amplification, a melt curve was performed from 70°C to 95°C, increasing 0.2°C every cycle with a 5-second hold. The CT values were averaged for each oligo pair for each set of technical replicates, and sample values were normalized to the housekeeping gene actin. The GFP shRNA transfectant line was used as a baseline control for comparison to the URE3-BP and Igl shRNA transfectant lines; HM1:IMSS samples were included as a secondary control. The differences in gene expression for the URE3-BP and Igl transfectant lines as compared to the GFP transfectant line were calculated by using both the relative standard curve and the comparative C(t) method (ΔΔ C(t) method) [54, 55]. Statistical analysis was performed using Student's t test (two-tailed), groups were also compared using ANOVA, and the GraphPad QuickCalcs P-value calculator  was used to calculate P-values.
Isolation of small RNAs
Three of the Igl shRNA transfectant lines, Igl (1198–1226), Igl (2412–2440), and Igl (2777–2805), as well as the two PATMK knockdown shRNA lines, PATMK (2273–2301) and PATMK (3552–3580), and the PATMK scrambled control , were grown in 25 cm2 tissue culture flasks, and selected with 30 μg/ml hygromycin, since this level of selection had yielded substantial knockdown of PATMK . Small RNAs were isolated from each sample as well as control nontransfected HM1:IMSS trophozoites using Ambion's mirVana™ miRNA Isolation Kit (Applied Biosystems/Ambion, Austin, TX, USA) as per the manufacturer's instructions.
Northern blotting of small RNAs
Oligo probes were designed to match the sense or antisense strands of each hairpin. Fifty μg of small RNAs were loaded per lane on a 12% denaturing acrylamide gel and transferred to Hybond™-N+ nylon membrane (Amersham Biosciences/GE Healthcare Biosciences Corp, Piscataway, NJ, USA) as per the manufacturer's instructions. rRNA bands were analyzed to insure equal RNA loading. Oligo probes matching to the sense or antisense strands of the hairpins were end-labelled with 32P and were hybridized with each corresponding sample blot strip overnight at 37°C overnight, washed with low and medium stringency conditions, and exposed overnight to film.
This work was supported by NIH grant AI 37941 to WAP. We thank Anindya Dutta for the suggestion to use the U6-driven shRNA system in E. histolytica. Girija Ramakrishnan provided the pGIR310 vector and designed the modifying polylinker. Carol Gilchrist provided the microarray data. Anuradha Lohia and Douglas Boettner were helpful with advice and useful discussions. Ellyn Moore and Christina Bousquet prepared the TYI medium and maintained the stock HM1:IMSS trophozoites.
- Eichinger L, Noegel AA: Comparative genomics of Dictyostelium discoideum and Entamoeba histolytica. Curr Opin Microbiol. 2005, 8 (5): 606-611. 10.1016/j.mib.2005.08.009.PubMedView Article
- Loftus B, Anderson I, Davies R, Alsmark UC, Samuelson J, Amedeo P, Roncaglia P, Berriman M, Hirt RP, Mann BJ, et al: The genome of the protist parasite Entamoeba histolytica. Nature. 2005, 433 (7028): 865-868. 10.1038/nature03291.PubMedView Article
- Petri WA, Haque R, Mann BJ: The bittersweet interface of parasite and host: lectin-carbohydrate interactions during human invasion by the parasite Entamoeba histolytica. Annu Rev Microbiol. 2002, 56: 39-64. 10.1146/annurev.micro.56.012302.160959.PubMedView Article
- Amoebiasis. Wkly Epidemiol Rec. 1997, 72 (14): 97-99.
- Fire A, Xu S, Montgomery MK, Kostas SA, Driver SE, Mello CC: Potent and specific genetic interference by double-stranded RNA in Caenorhabditis elegans. Nature. 1998, 391 (6669): 806-811. 10.1038/35888.PubMedView Article
- Novina CD, Sharp PA: The RNAi revolution. Nature. 2004, 430 (6996): 161-164. 10.1038/430161a.PubMedView Article
- Plasterk RH: RNA silencing: the genome's immune system. Science. 2002, 296 (5571): 1263-1265. 10.1126/science.1072148.PubMedView Article
- Zamore PD, Tuschl T, Sharp PA, Bartel DP: RNAi: double-stranded RNA directs the ATP-dependent cleavage of mRNA at 21 to 23 nucleotide intervals. Cell. 2000, 101 (1): 25-33. 10.1016/S0092-8674(00)80620-0.PubMedView Article
- Bernstein E, Caudy AA, Hammond SM, Hannon GJ: Role for a bidentate ribonuclease in the initiation step of RNA interference. Nature. 2001, 409 (6818): 363-366. 10.1038/35053110.PubMedView Article
- Elbashir SM, Lendeckel W, Tuschl T: RNA interference is mediated by 21- and 22-nucleotide RNAs. Genes Dev. 2001, 15 (2): 188-200. 10.1101/gad.862301.PubMed CentralPubMedView Article
- Agrawal N, Dasaradhi PV, Mohmmed A, Malhotra P, Bhatnagar RK, Mukherjee SK: RNA interference: biology, mechanism, and applications. Microbiol Mol Biol Rev. 2003, 67 (4): 657-685. 10.1128/MMBR.67.4.657-685.2003.PubMed CentralPubMedView Article
- Dykxhoorn DM, Novina CD, Sharp PA: Killing the messenger: short RNAs that silence gene expression. Nat Rev Mol Cell Biol. 2003, 4 (6): 457-467. 10.1038/nrm1129.PubMedView Article
- Pak J, Fire A: Distinct populations of primary and secondary effectors during RNAi in C. elegans. Science. 2007, 315 (5809): 241-244. 10.1126/science.1132839.PubMedView Article
- Sijen T, Fleenor J, Simmer F, Thijssen KL, Parrish S, Timmons L, Plasterk RH, Fire A: On the role of RNA amplification in dsRNA-triggered gene silencing. Cell. 2001, 107 (4): 465-476. 10.1016/S0092-8674(01)00576-1.PubMedView Article
- De S, Pal D, Ghosh SK: Entamoeba histolytica : computational identification of putative microRNA candidates. Exp Parasitol. 2006, 113 (4): 239-243. 10.1016/j.exppara.2006.01.009.PubMedView Article
- Abed M, Ankri S: Molecular characterization of Entamoeba histolytica RNase III and AGO2, two RNA interference hallmark proteins. Exp Parasitol. 2005, 110 (3): 265-269. 10.1016/j.exppara.2005.02.023.PubMedView Article
- MacFarlane RC, Shah PH, Singh U: Transcriptional profiling of Entamoeba histolytica trophozoites. Int J Parasitol. 2005, 35 (5): 533-542. 10.1016/j.ijpara.2005.02.006.PubMedView Article
- Kaur G, Lohia A: Inhibition of gene expression with double strand RNA interference in Entamoeba histolytica. Biochem Biophys Res Commun. 2004, 320 (4): 1118-1122. 10.1016/j.bbrc.2004.06.064.PubMedView Article
- Lohia A: The cell cycle of Entamoeba histolytica. Mol Cell Biochem. 2003, 253: 217-222. 10.1023/A:1026055631421.PubMedView Article
- Vayssie L, Vargas M, Weber C, Guillen N: Double-stranded RNA mediates homology-dependent gene silencing of gamma-tubulin in the human parasite Entamoeba histolytic a. Mol Biochem Parasitol. 2004, 138 (1): 21-28. 10.1016/j.molbiopara.2004.07.005.PubMedView Article
- Petri WA, Ramakrishnan G: Applying antisense technology to the study of Entamoeba histolytica pathogenesis. Trends Microbiol. 1999, 7 (12): 471-474. 10.1016/S0966-842X(99)01626-1.PubMedView Article
- Das S, Lohia A: Delinking of S phase and cytokinesis in the protozoan parasite Entamoeba histolytica. Cell Microbiol. 2002, 4 (1): 55-60. 10.1046/j.1462-5822.2002.00165.x.PubMedView Article
- Gangopadhyay SS, Ray SS, Kennady K, Pande G, Lohia A: Heterogeneity of DNA content and expression of cell cycle genes in axenically growing Entamoeba histolytica HM1:IMSS clone A. Mol Biochem Parasitol. 1997, 90 (1): 9-20. 10.1016/S0166-6851(97)00156-4.PubMedView Article
- Bracha R, Nuchamowitz Y, Mirelman D: Inhibition of gene expression in Entamoeba by the transcription of antisense RNA: effect of 5' and 3' regulatory elements. Mol Biochem Parasitol. 2000, 107 (1): 81-90. 10.1016/S0166-6851(00)00181-X.PubMedView Article
- Dastidar PG, Majumder S, Lohia A: Eh Klp5 is a divergent member of the kinesin 5 family that regulates genome content and microtubular assembly in Entamoeba histolytica. Cell Microbiol. 2007, 9 (2): 316-328. 10.1111/j.1462-5822.2006.00788.x.PubMedView Article
- MacFarlane RC, Singh U: Identification of an Entamoeba histolytica serine-, threonine-, and isoleucine-rich protein with roles in adhesion and cytotoxicity. Eukaryot Cell. 2007, 6 (11): 2139-2146. 10.1128/EC.00174-07.PubMed CentralPubMedView Article
- Yu JY, DeRuiter SL, Turner DL: RNA interference by expression of short-interfering RNAs and hairpin RNAs in mammalian cells. Proc Natl Acad Sci USA. 2002, 99 (9): 6047-6052. 10.1073/pnas.092143499.PubMed CentralPubMedView Article
- Brummelkamp TR, Bernards R, Agami R: A system for stable expression of short interfering RNAs in mammalian cells. Science. 2002, 296 (5567): 550-553. 10.1126/science.1068999.PubMedView Article
- Das G, Henning D, Wright D, Reddy R: Upstream regulatory elements are necessary and sufficient for transcription of a U6 RNA gene by RNA polymerase III. EMBO J. 1988, 7 (2): 503-512.PubMed CentralPubMed
- Gou D, Jin N, Liu L: Gene silencing in mammalian cells by PCR-based short hairpin RNA. FEBS Lett. 2003, 548 (1–3): 113-118. 10.1016/S0014-5793(03)00630-6.PubMedView Article
- Silva JM, Li MZ, Chang K, Ge W, Golding MC, Rickles RJ, Siolas D, Hu G, Paddison PJ, Schlabach MR, et al: Second-generation shRNA libraries covering the mouse and human genomes. Nat Genet. 2005, 37 (11): 1281-1288.PubMed
- Kim DH, Behlke MA, Rose SD, Chang MS, Choi S, Rossi JJ: Synthetic dsRNA Dicer substrates enhance RNAi potency and efficacy. Nat Biotechnol. 2005, 23 (2): 222-226. 10.1038/nbt1051.PubMedView Article
- Cheng XJ, Tsukamoto H, Kaneda Y, Tachibana H: Identification of the 150-kDa surface antigen of Entamoeba histolytica as a galactose- and N-acetyl-D-galactosamine-inhibitable lectin. Parasitol Res. 1998, 84 (8): 632-639. 10.1007/s004360050462.PubMedView Article
- Cheng XJ, Hughes MA, Huston CD, Loftus B, Gilchrist CA, Lockhart LA, Ghosh S, Miller-Sims V, Mann BJ, Petri WA, et al: Intermediate subunit of the Gal/GalNAc lectin of Entamoeba histolytica is a member of a gene family containing multiple CXXC sequence motifs. Infect Immun. 2001, 69 (9): 5892-5898. 10.1128/IAI.69.9.5892-5898.2001.PubMed CentralPubMedView Article
- Beck DL, Boettner DR, Dragulev B, Ready K, Nozaki T, Petri WA: Identification and gene expression analysis of a large family of transmembrane kinases related to the Gal/GalNAc lectin in Entamoeba histolytica. Eukaryot Cell. 2005, 4 (4): 722-732. 10.1128/EC.4.4.722-732.2005.PubMed CentralPubMedView Article
- Gilchrist CA, Leo M, Line CG, Mann BJ, Petri WA: Calcium modulates promoter occupancy by the Entamoeba histolytica Ca2+-binding transcription factor URE3-BP. J Biol Chem. 2003, 278 (7): 4646-4653. 10.1074/jbc.M211271200.PubMedView Article
- Okada M, Huston CD, Oue M, Mann BJ, Petri WA, Kita K, Nozaki T: Kinetics and strain variation of phagosome proteins of Entamoeba histolytica by proteomic analysis. Mol Biochem Parasitol. 2006, 145 (2): 171-183. 10.1016/j.molbiopara.2005.10.001.PubMedView Article
- Nalefski EA, Falke JJ: The C2 domain calcium-binding motif: structural and functional diversity. Protein Science. 1996, 5: 2375-2390. 10.1002/pro.5560051201.PubMed CentralPubMedView Article
- Boettner DR, Huston CD, Linford AS, Buss SN, Houpt E, Sherman NE, Petri WA: Entamoeba histolytica phagocytosis of human erythrocytes involves PATMK, a member of the transmembrane kinase family. PLoS Pathog. 2008, 4 (1): e8-10.1371/journal.ppat.0040008.PubMed CentralPubMedView Article
- Miranda R, Salgado LM, Sanchez-Lopez R, Alagon A, Lizardi PM: Identification and analysis of the u6 small nuclear RNA gene from Entamoeba histolytica. Gene. 1996, 180 (1–2): 37-42. 10.1016/S0378-1119(96)00397-6.PubMedView Article
- Vines RR, Purdy JE, Ragland BD, Samuelson J, Mann BJ, Petri WA: Stable episomal transfection of Entamoeba histolytica. Mol Biochem Parasitol. 1995, 71 (2): 265-267. 10.1016/0166-6851(95)00057-8.PubMedView Article
- Hamann L, Buss H, Tannich E: Tetracycline-controlled gene expression in Entamoeba histolytica. Mol Biochem Parasitol. 1997, 84 (1): 83-91. 10.1016/S0166-6851(96)02771-5.PubMedView Article
- Hamann L, Nickel R, Tannich E: Transfection and continuous expression of heterologous genes in the protozoan parasite Entamoeba histolytica. Proc Natl Acad Sci USA. 1995, 92 (19): 8975-8979. 10.1073/pnas.92.19.8975.PubMed CentralPubMedView Article
- Shao Y, Chan CY, Maliyekkel A, Lawrence CE, Roninson IB, Ding Y: Effect of target secondary structure on RNAi efficiency. RNA. 2007, 13: 1631-1640. 10.1261/rna.546207.PubMed CentralPubMedView Article
- Gredell JA, Berger AK, Walton SP: Impact of target mRNA structure on siRNA silencing efficiency: a large-scale study. Biotechnol Bioeng. 2008, 100: 744-755. 10.1002/bit.21798.PubMed CentralPubMedView Article
- Gilchrist CA, Houpt E, Trapaidze N, Fei Z, Crasta O, Asgharpour A, Evans C, Martino-Catt S, Baba DJ, Stroup S, et al: Impact of intestinal colonization and invasion on the Entamoeba histolytica transcriptome. Mol Biochem Parasitol. 2006, 147 (2): 163-176. 10.1016/j.molbiopara.2006.02.007.PubMedView Article
- Diamond LS, Harlow DR, Cunnick CC: A new medium for the axenic cultivation of Entamoeba histolytica and other Entamoeba. Trans R Soc Trop Med Hyg. 1978, 72 (4): 431-432. 10.1016/0035-9203(78)90144-X.PubMedView Article
- Huston CD, Boettner DR, Miller-Sims V, Petri WA: Apoptotic killing and phagocytosis of host cells by the parasite Entamoeba histolytica. Infect Immun. 2003, 71 (2): 964-972. 10.1128/IAI.71.2.964-972.2003.PubMed CentralPubMedView Article
- Ramakrishnan G, Vines RR, Mann BJ, Petri WA: A tetracycline-inducible gene expression system in Entamoeba histolytica. Mol Biochem Parasitol. 1997, 84 (1): 93-100. 10.1016/S0166-6851(96)02784-3.PubMedView Article
- Katz U, Bracha R, Nuchamowitz Y, Milstein O, Mirelman D: Comparison between constitutive and inducible plasmid vectors used for gene expression in Entamoeba histolytica. Mol Biochem Parasitol. 2003, 128 (2): 229-233. 10.1016/S0166-6851(03)00076-8.PubMedView Article
- The Ambion/Applied Biosystems siRNA Target Finder. [http://www.ambion.com/techlib/misc/siRNA_finder.html]
- TIGR Database Entamoeba histolytica Genome Project. [http://www.tigr.org/tdb/e2k1/eha1/]
- GraphPad QuickCalcs. [http://www.graphpad.com/quickcalcs]
- Cikos S, Bukovska A, Koppel J: Relative quantification of mRNA: comparison of methods currently used for real-time PCR data analysis. BMC Mol Biol. 2007, 8: 113-10.1186/1471-2199-8-113.PubMed CentralPubMedView Article
- Real-Time PCR: M. Teyfik Dorak, MD, PhD. [http://www.dorak.info/genetics/realtime.html]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.