- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Journal of Biomedicine and Biotechnology
Volume 2010 (2010), Article ID 915980, 8 pages
High-Throughput Sequencing of MicroRNAs in Adenovirus Type 3 Infected Human Laryngeal Epithelial Cells
1State Key Laboratory of Bioelectronics, School of Biological Science and Medical Engineering, Southeast University, SiPaiLou2#, Nanjing 210096, China
2Key Laboratory of Enteric Pathogenic Microbiology, Ministry of Health, Microbiological Laboratory, Jiangsu Center for Disease Prevention and Control (CDC), 172 Jiangsu Rd, Nanjing 210009, China
Received 11 February 2010; Revised 27 March 2010; Accepted 31 March 2010
Academic Editor: Momiao Xiong
Copyright © 2010 Yuhua Qi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Adenovirus infection can cause various illnesses depending on the infecting serotype, such as gastroenteritis, conjunctivitis, cystitis, and rash illness, but the infection mechanism is still unknown. MicroRNAs (miRNA) have been reported to play essential roles in cell proliferation, cell differentiation, and pathogenesis of human diseases including viral infections. We analyzed the miRNA expression profiles from adenovirus type 3 (AD3) infected Human laryngeal epithelial (Hep2) cells using a SOLiD deep sequencing. 492 precursor miRNAs were identified in the AD3 infected Hep2 cells, and 540 precursor miRNAs were identified in the control. A total of 44 miRNAs demonstrated high expression and 36 miRNAs showed lower expression in the AD3 infected cells than control. The biogenesis of miRNAs has been analyzed, and some of the SOLiD results were confirmed by Quantitative PCR analysis. The present studies may provide a useful clue for the biological function research into AD3 infection.
Adenovirus (AD) belongs to a family of nonenveloped icosahedral viruses containing a double-stranded DNA genome, Adenoviridae, and genus Mastadenovirus. There are 53 human AD serotypes divided into seven subgroups (A to G) based on immunologic, biologic, and biochemical characteristics . Among these, viruses in groups B1, C, and E cause respiratory infections, group B2 viruses infect the kidney and urinary tract; group F viruses cause gastroenteritis; and several group D serotypes are associated with epidemic keratoconjunctivitis . Most people recover from AD infections by themselves, but people with immunodeficiency and small children are at higher risk of fatal adenovirus infection. Studies showed that AD can evade host immune responses, such as inhibition of interferon functions by RNA and E1A, and inhibition of intrinsic cellular apoptosis in infected cells . The interaction between the virus and its host is still not clear. Recently, several studies evaluated human adenoviruses which play an important role in the study of tumorigenesis, including oncogenic transformation of cells, immune evasion, angiogenesis, and metastasis [3, 4]. There are no antiviral drugs to treat adenoviral infections, so treatment is largely directed at the symptoms. The discovery of the AD yielded further evidence that viruses can sometimes interfere with each other during the growth and replication within a host animal. Other researchers found that this interference can be mediated by interferon produced by the host animal. So interferon has become prominent in the treatment of a variety of cancers and infectious diseases, such as hepatitis C . Recently, several studies demonstrated that RNA interference (RNAi) can also be used for treatment with viral infection in vitro, such as hepatitis B virus, coxsackievirus B3, and Coxsackievirus A16 [6–8]. These findings indicate that those specific siRNAs targeted viruses are not as effective as expected.
RNAi can also be triggered by microRNAs (miRNAs), which are short ~22 nucleotide RNA sequences that bind complementary sequences in the UTR of multiple target mRNAs, playing an important role in the development, proliferation, differentiation and apoptosis of organisms. Growing evidences have indicated a strong association between microRNAs and various viral infections. For example, miR-125b, mirR-150, miR-382, and miR-223 have been reported to contribute to the maintenance of human immunodeficiency virus-1 (HIV-1) latency in resting CD4+ T cells . Epstein-Barr virus and Kaposi sarcoma herpes viruses (KSHVs) have been reported to encode miRNAs , but the functions of most of them are not known. miRNAs have been reported to be used as a tool for gene specific therapeutics against viral infections . The potential use of miRNAs as novel noninvasive biomarkers for diagnosis has been suggested in other studies [12–14].
In this study, we used a strategy of initial screening by Applied Biosystems’ next generation sequencing system which is Sequencing by Oligonucleotide Ligation and Detection (SOLiD) and validation by quantitative RT-PCR (qRT-PCR), to analyze the different miRNA expression profiling in AD3 infected Hep2 cells.
2. Materials and Methods
2.1. Cell Culture and Virus Infection
Human laryngeal epithelial (Hep2) cells were cultivated in complete RPMI medium 1640 (Invitrogen) supplemented with 10% fetal calf serum (FCS) at in an atmosphere of 5% CO2. When the Hep2 cells were grown in 25-cm2 flasks reached to 70% confluence, they were infected with AD3 at a multiplicity of infection (m.o.i.) of 0.4 under 50% tissue culture infectious doses (TCID50), and maintained after infection at in RPMI 1640 medium with 2% FCS.
2.2. RNA Isolation
Hep2 cells were infected with AD3 as described above. At 6 hours, 24 hours, 48 hours, and 72 hours post infection (p.i.), RNA samples were prepared from AD3 infected Hep2 cells and controlled Hep2 cells using a mirVana miRNA Isolation Kit (Ambion): 102–107 cultured cells were washed in 1% PBS for three times, washed cells were mixed with 600 l Lysis/Binding Solution and 1/10 volume of miRNA Homogenate Additive thoroughly by vortexing and incubated on ice for 10 minutes An equal volume of acid/phenol/chloroform (Ambion) was then added to the each aliquot. The solution was centrifuged for 5 minutes at 10,000 g at room temperature. For miRNA isolation, the extraction was mixed thoroughly with 1/3 volume of 100% ethanol and passed through a column (Ambion). The flow-through was contained the small RNAs. The filtrate was then mixed thoroughly with 2/3 volume 100% ethanol and passed through the second column (Ambion). The column was washed following the protocol, and miRNAs were eluted in 100 l of elution buffer (). For total RNA isolation, the extraction was mixed thoroughly with 1.25 volume of 100% ethanol and passed through a column (Ambion). The column was washed following the protocol and total RNA were eluted in 100 l of elution buffer (). A total of 1 l of the eluate was quantified by Nanodrop.
2.3. SOLiD Sequencing and Sequence Analysis
MiRNA samples (100 ng) isolated from AD3 infected Hep2 cells and control Hep2 cells were processed into sequencing libraries using the Small RNA Expression Kit (Applied Biosystems). Briefly, RNA was ligated overnight with the adapters from the kit, reverse transcribed, RNAse H-treated, and PCR-amplified before size selection on Agarose gels to contain 16–61 nt of inserted sequences. Libraries were amplified onto beads using emulsion PCR, deposited on slides, and sequenced using the SOLiD v 2 sequencing system (Applied Biosystems)  at the State Key Lab of Bioelectronics Laboratory, Southeast University of China. Finally, SOLiD data was first analyzed by SOLiD System Small RNA Analysis Pipeline Tool (RNA2MAP). Firstly, the parameters set for alignments were the maximum length (18 nt) and the tolerance of 3 mismatches when generating initial seeds locations. In the extension step, less than 6 mismatches were allowed in full length mapping. The miRBase sequences (Sanger) of human being was download from miRBase (http://www.mirbase.org/). Acceptable sequences were compared to the miRBase database (release 14.0). To increase signal above noise, we conservatively selected only those alignments corresponding to beads sampled a minimum of 10 times in any of the libraries.
2.4. Quantitative RT-PCR (qRT-PCR)
To confirm the expression of miRNAs by deep sequencing approach, stem-loop quantitative RT-PCR (qRT-PCR) was performed. Briefly, cDNA was synthesized from total RNA by using AMV reverse transcriptase (TaKaRa). The 20 L reactions were incubated for 15 minutes at , 30 minutes at and 5 minutes at , and then held at . Subsequently, the 20 L PCR reactions included 1 l RT-PCR product, 10 l Premix Ex Taq (TaKaRa), and 1 l SYBR green (Invitrogen). The reactions were incubated in a 96-well optical plate at for 10 minutes, followed by 40 cycles of for 15 seconds, and for 1 minute. Q-PCR assay was performed in triplicate using an ABI 7500 machine (Applied Biosystems). The expression levels of candidate miRNAs were measured in terms of threshold cycle value () and normalized to Human small nuclear RNA U6. The ratio of two groups of miRNAs was calculated by using the equation , in which , and .
2.5. Statistical Analysis
We firstly used the Z-test to determine the statistical significance of the differences between the two libraries. This approach is to look at the number of copies of a speciﬁc miRNA per cell as a fraction or proportion of the total number of miRNA molecules in that cell. The same proportion (p) of speciﬁc tags should be present in the miRNA library of all sequenced tags  So the Z-value for a specific miRNA could be calculated as following: The proportion , the expected proportion when the null hypothesis is true, is calculated as . The Z statistic is approximately normally distributed and can be compared with critical Z value for the two-sided significance level .
A Z-test was used to deal with the raw data by considering both the copies and fold change of miRNAs. Candidate miRNAs should be fulﬁl three criteria: (1) mean fold change 2, (2) having at least 10 copies by SOLiD sequencing, and (3) Z score 1.96 or −1.96 for .
3.1. Replication of AD3 in Hep2 Cells
Cytopathic effect (CPE) was first observed at 24 hours post infection and progressed to moderate and severe CPE at 72 hours and 96 hours, respectively. miRNA extracted from cells after 72 hours of infection was used for SOLiD deep sequencing.
3.2. SOLiD Sequencing of miRNAs from Hep2 Cells Infected by AD3 and Control Hep2 Cells
To search for the miRNAs expression in AD3 infected Hep2 cells and the control, we prepared these two miRNA samples for sequencing by the SOLiD system. 358,297 effective reads were obtained in AD3 infected cells and 479,414 effective reads in control cells. After filtrating reads contaminated by rRNAetc (rRNA, tRNA, snRNA, and snoRNA), repeat fractions, and degraded fractions, the remaining effective reads were mapped to the human precursor libraries. 37,800 and 78,143 reads were obtained from infected and control cells, respectively. In these reads, the most abundant size class was 22 nt in both of AD3 infected Hep2 cells and the control. The percentage was 29.3% and 34.6% respectively. Following was 23 nt (27.7% in AD3 infected Hep2 cells and 26.3% in the control) (Figure 1). Compared to the miRBase (14.0), 492 precursor, miRNA out of 734 known precursor miRNAs were identified among the AD3 infected Hep2 cells, and 540 precursor miRNAs were identified among the control Hep2 cells. All of the known miRNAs have been characterized in the Sanger miRBase sequence database. Figure 2 showed the distribution of genomic loci in the chromosome. All of the human chromosome were located in the precursor miRNAs. Chromosomes 19 and X contain more precursor miRNAs (12%, and 11%, resp.). 314 miRNAs (51%) were located in introns of protein-coding genes, 231 miRNAs (38%) were located in intergenic regions, while only 43 miRNAs (7%) were located in exons of noncoding RNAs or UTR of protein coding genes. Interestedly, 23 miRNAs (4%) were located in either an exon or an intron depending on alternative splicing of the host transcript (Table 1).
3.3. Expression Profiles of miRNAs in AD3 Infected Cells
To study the differential expression profile of known miRNAs between AD3 infected Hep2 cells and the control, both the Z-test and fold-change analysis were made based on the sequencing results. There are 80 precursor miRNAs to meet the three criteria described above, and could be selected as candidates. Among these candidates, 44 precursor miRNAs were shown upregulated expression in AD3 infected cells (Table 2). Z-test value was from 1.99 to 49.33. 19 out of 44 upregulated precursor miRNAs showed higher values (over 10). Among this, miR-1975 and miR-1302-2 have the highest Z values with 49.33 and 31.13, respectively. Table 2 showed the ratio of miRNAs expression level of AD3 infected cells versus control cells. Among this, miR-1975 and miR-1302-2 have the higher ratio (7.47-fold and 7.31-fold). 36 precursor miRNAs were shown downregulated expression in AD3 infected cells (Table 3). Among this there are 11 miRNAs that have the higher Z values (over 10). Table 3 showed the ratio of AD3 infected cells/control cells (expression level for precursor miRNAs). Among this miR-181b-2 and miR-1274b showed lowest change in AD3 infected cells (0.19 fold) compared with control cells.
3.4. Q-PCR Detection for miRNAs Profiles in AD3 Infected Cells
To confirm the differential expression of miRNAs in AD3 infected Hep2 cells, we performed some miRNAs tests with stem-loop Q-PCR assays in total RNA isolated from AD3 infected Hep2 cells and the control, respectively. The expression levels of candidate miRNAs were measured in terms of threshold cycle value () and normalized to Human small nuclear RNA U6. The ratio of miRNAs in two groups was calculated by using the equation . The results showed that there was general consistency between Q-PCR assay and SOLiD sequencing analysis. The seven upregulated miRNAs in SOLiD results showed increased ratios measured in Q-PCR assays (Figure 3). For example, miR-1974, miR-1975, and miR-7 showed increasing values of 2.51, 4.87 and 2.44 respectively measured in Q-PCR assay, and increasing values of 2.95, 7.47 and 6.34 were obtained in SOLiD sequencing (Table 2). The seven downregulated miRNAs in SOLiD results showed decrease ratios measured in Q-PCR assays (Figure 3). For example, the results showed decreased expression levels in miR-27b, miR-125b, and miR-27a in AD3 infected Hep2 cells with fold change of 0.74, 0.75, and 0.77 respectively (Figure 3). While in SOLiD sequencing analysis, these miRNAs have shown decreased expression with fold change 0.43, 0.41 and 0.38, respectively (Table 3). To further understand the miRNAs expression in AD3 infected Hep2 cells at 6 hours, 24 hours, 48 hours, and 72 hours post infection, we performed Q-PCR experiments on miRNAs chosen from Figure 3. The ratio of miRNAs in AD3 infected cells versus controls was calculated by using the equation . Surprisingly, all of the miRNAs shown in Table 4 exhibited increased expression in AD3 infected cells at 6 hours post infection. However, at 48 hours p.i. and 72 hours p.i., miR-27b, miR-125b, and miR-181b showed downregulated expression in AD3 infected cells comparing controls (Table 4). At 48 hours p.i. and 72 hours p.i., miR-1974 showed increased expression in AD3 infected cells with fold change 1.36 and 1.26 respectively. However, at 24 hours p.i., miR-1974 showed decreased expression in AD3 infected cells. Table 4 showed that the miR-17 expression was upregulated at 6, 24, and 72 hours p.i. in AD3 infected cells. But at 48 hours p.i., miR-17 showed downregulated expression in AD3 infected cells.
As second-generation sequencing platforms have matured, all three major high-throughput sequencing systems—Illumina's Genome Analyzer (Solexa), Life Technologies' ABI SOLiD, and Roche's 454 GS FLX have been used to study miRNA expression profiles associated with disease such as cancer, viral and metabolic diseases. Compared with the microarray technology, deep sequencing can generate millions of small RNA sequence reads from a given sample, giving possible to discover the novel miRNAs or study the function of virus encoded miRNAs by using bioinformatics and molecular biology tools. miRNAs have recently been recognized as major regulators of gene expression in various viral infections. Virus-encoded miRNAs seem to evolve rapidly and regulate both the viral life cycle and the interaction between viruses and their hosts. In the present study, we have screened out some miRNA candidates that have obvious differential expression from AD3 infected Hep2 cells by SOLiD deep sequencing system.
This study identified 492 precursor miRNAs in the AD3 infected cells and more precursor miRNAs in control cells by mapping the human miRBase database. Out of these almost 38% were located in intergenic regions, which usually contain their own miRNA gene promoter and regulatory units . However, as much as 51% miRNAs were located in introns, and only 7% were located in exons of noncoding RNAs or UTR of protein coding genes. These usually show a concurrent transcription and regulation expression profile originating from a common promoter with their host genes [19–21]. Some pre-miRNAs existing in drosphila and c elegants have been reported to be spliced from the introns in which they reside without having to undergo the microprocessor machinery. These miRNAs are called mirtrons [22, 23]. Since our study identified 51% miRNAs were located in introns, we think more mirtons could be discovered through bioinformatics and molecular biology tools.
The interaction between the AD and its host is still unknown. One of the mechanisms is that the virus must interfere with the host cell antiviral defense mechanisms to maximize escape and spread of the progeny virus during a virus infection . Previous studies determine the changes in the host cell gene expression profiles upon infection with AD using DNA microarray technique [24, 25]. For example, many transcriptions factor E2F-dependent host cell genes showed different expression during AD infection. 45% of the upregulated genes and 25% of the downregulated genes contained potential E2F binding sites . It is showen that E2F-dependent host cell genes are subjected to a selective regulation. However, the mechanism of AD infection after transcription is still unknown. miRNA can regulate gene expression through translational repression or mRNA cleavage. Previous research has reported that the expression of E2F1 is negatively regulated by two c-myc-regulated miRNAs, miR-17-5p, and miR-20a . In our study, upregulated miR-17 may target with transcription factor E2F1 , playing important role in cell differentiation, proliferation, and apoptosis. Out of these, miR-17 has been reported to be overexpression in B-cell lymphoma samples . The high miR-17 cluster miRNAs suppress cell death . This suggests that during a virus infection, cellular apoptosis must delay long enough and the virus needs to establish optimal conditions for replication to ensure efficient production of progeny virus. Our study may provide a useful clue for the biological function research into AD3 infected host cells. Further investigation needs to clarify the roles of identified miRNAs in the pathogenesis of AD3 infected host cells.
SOLiD sequencing provides a useful method for identification of the miRNAs profiles in AD3 infected Hep2 cells. 492 precursor miRNAs were identified in the AD3 infected Hep2 cells, and 540 precursor miRNAs were identified in the control. A total of 44 miRNAs demonstrated high expression and 36 miRNAs showed lower expression in the AD3 infected cells than the control. Further studies are required to identify the miRNA target genes and the functions of the miRNAs in the complex molecular network regulation during the virus infection host cells using bioinformatics tools.
The paper was supported by the National Natural Science Foundation of China (30901285), the Natural Science Foundation of Jiangsu Province of China (SBK200922783), and the Open Research Fund of State Key Laboratory of Bioelectronics, Southeast University of China. Both Yuhua Qi and Jin Tu contributed equally to this work.
- M. Echavarría, “Adenoviruses in immunocompromised hosts,” Clinical Microbiology Reviews, vol. 21, no. 4, pp. 704–715, 2008.
- J. A. Mahr and L. R. Gooding, “Immune evasion by adenoviruses,” Immunological Reviews, vol. 168, pp. 121–130, 1999.
- H. Guan, J. F. Williams, and R. P. Ricciardi, “Induction of neuronal and tumor-related genes by adenovirus type 12 E1A,” Journal of Virology, vol. 83, no. 2, pp. 651–661, 2009.
- G. Chinnadurai, “Modulation of oncogenic transformation by the human adenovirus E1A C-terminal region,” Current Topics in Microbiology and Immunology, vol. 273, pp. 139–161, 2004.
- K. A. Sharieff, D. Duncan, and Z. Younossi, “Advances in treatment of chronic hepatitis C: ‘pegylated’ interferons,” Cleveland Clinic Journal of Medicine, vol. 69, no. 2, pp. 155–159, 2002.
- Z. Wu, Y. Gao, L. Sun, P. Tien, and Q. Jin, “Quick identification of effective small interfering RNAs that inhibit the replication of coxsackievirus A16,” Antiviral Research, vol. 80, no. 3, pp. 295–301, 2008.
- K.-L. Wu, X. Zhang, J. Zhang, et al., “Inhibition of hepatitis B virus gene expression by single and dual small interfering RNA treatment,” Virus Research, vol. 112, no. 1-2, pp. 100–107, 2005.
- J. Yuan, P. K. M. Cheung, H. M. Zhang, D. Chau, and D. Yang, “Inhibition of coxsackievirus B3 replication by small interfering RNAs requires perfect sequence match in the central region of the viral positive strand,” Journal of Virology, vol. 79, no. 4, pp. 2151–2159, 2005.
- J. Huang, F. Wang, E. Argyris, et al., “Cellular microRNAs contribute to HIV-1 latency in resting primary CD4+ T lymphocytes,” Nature Medicine, vol. 13, no. 10, pp. 1241–1247, 2007.
- V. Nair and M. Zavolan, “Virus-encoded microRNAs: novel regulators of gene expression,” Trends in Microbiology, vol. 14, no. 4, pp. 169–175, 2006.
- S.-Y. Ying and S.-L. Lin, “Current perspectives in intronic micro RNAs (miRNAs),” Journal of Biomedical Science, vol. 13, no. 1, pp. 5–15, 2006.
- J. Skog, T. Wurdinger, S. van Rijn, et al., “Glioblastoma microvesicles transport RNA and proteins that promote tumour growth and provide diagnostic biomarkers,” Nature Cell Biology, vol. 10, no. 12, pp. 1470–1476, 2008.
- K. E. Resnick, H. Alder, J. P. Hagan, D. L. Richardson, C. M. Croce, and D. E. Cohn, “The detection of differentially expressed microRNAs from the serum of ovarian cancer patients using a novel real-time PCR platform,” Gynecologic Oncology, vol. 112, no. 1, pp. 55–59, 2009.
- S. Gilad, E. Meiri, Y. Yogev, et al., “Serum microRNAs are promising novel biomarkers,” PLoS ONE, vol. 3, no. 9, article e3148, 2008.
- L. A. Goff, J. Davila, M. R. Swerdel, et al., “Ago2 immunoprecipitation identifies predicted MicroRNAs in human embryonic stem cells and neural precursors,” PLoS ONE, vol. 4, no. 9, article e7192, 2009.
- T. D. Schmittgen and K. J. Livak, “Analyzing real-time PCR data by the comparative CT method,” Nature Protocols, vol. 3, no. 6, pp. 1101–1108, 2008.
- A. J. Kal, A. J. van Zonneveld, V. Benes, et al., “Dynamics of gene expression revealed by comparison of serial analysis of gene expression transcript profiles from yeast grown on two different carbon sources,” Molecular Biology of the Cell, vol. 10, no. 6, pp. 1859–1872, 1999.
- N. C. Lau, L. P. Lim, E. G. Weinstein, and D. P. Bartel, “An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans,” Science, vol. 294, no. 5543, pp. 858–862, 2001.
- Y.-K. Kim and V. N. Kim, “Processing of intronic microRNAs,” The EMBO Journal, vol. 26, no. 3, pp. 775–783, 2007.
- S. Baskerville and D. P. Bartel, “Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes,” RNA, vol. 11, no. 3, pp. 241–247, 2005.
- A. Rodriguez, S. Griffiths-Jones, J. L. Ashurst, and A. Bradley, “Identification of mammalian microRNA host genes and transcription units,” Genome Research, vol. 14, no. 10, pp. 1902–1910, 2004.
- K. Okamura, J. W. Hagen, H. Duan, D. M. Tyler, and E. C. Lai, “The mirtron pathway generates microRNA-class regulatory RNAs in Drosophila,” Cell, vol. 130, no. 1, pp. 89–100, 2007.
- J. G. Ruby, C. H. Jan, and D. P. Bartel, “Intronic microRNA precursors that bypass Drosha processing,” Nature, vol. 448, no. 7149, pp. 83–86, 2007.
- H. Zhao, F. Granberg, L. Elfineh, U. Pettersson, and C. Svensson, “Strategic attack on host cell gene expression during adenovirus infection,” Journal of Virology, vol. 77, no. 20, pp. 11006–11015, 2003.
- A. Dorn, H. Zhao, F. Granberg, et al., “Identification of specific cellular genes up-regulated late in adenovirus type 12 infection,” Journal of Virology, vol. 79, no. 4, pp. 2404–2412, 2005.
- K. A. O'Donnell, E. A. Wentzel, K. I. Zeller, C. V. Dang, and J. T. Mendell, “c-Myc-regulated microRNAs modulate E2F1 expression,” Nature, vol. 435, no. 7043, pp. 839–843, 2005.
- L. He, J. M. Thomson, M. T. Hemann, et al., “A microRNA polycistron as a potential human oncogene,” Nature, vol. 435, no. 7043, pp. 828–833, 2005.