[Frontiers in Bioscience E4, 2579-2585, June 1, 2012]

MicroRNA-27a activity is not suppressed in porcine oocytes

Li Chen1, Xiaoxiang Hu1, Yunping Dai2, Qinghe Li1, Xiaobo Wang3, Qiuyan Li2, Kai Xue4,Yan Li2, Jingwen Liang1, Yuanyuan Wang1, Xiaojuan Liu1, Ning Li1

1State Key Laboratory for AgroBiotechnology, College of Biological Sciences, China Agricultural University, Beijing 100193, China, 2Ji Pulin biotech Ltd, Beijing 100193, China, 3Center for Cardiovascular Sciences, Albany Medical College, Albany, NY 12208, 4Dalian Medical University, Dalian 116044, China

TABLE OF CONTENTS

1. Abstract
2. Introduction
3. Materials and methods
3.1. Oocyte collection and culture
3.2. Evaluation of meiotic stages of porcine oocytes
3.3. MicroRNA expression profiling assay and analysis
3.4. Taqman real-time PCR of miRNAs
3.5. Bioinformatics method
3.6. Microinjection
3.7. Western blot analysis
3.8. Transfection and luciferase assay
3.9. Statistical analysis
4. Results
4.1. Meiotic stages of porcine oocytes during in vitro maturation
4.2. miRNA expression profiling during porcine oocyte maturation
4.3. miR-27a facilitates the first cell cleavage
4.4. miR-27a represses translation of MAP2K4
5. Discussion
6. Acknowledgment
7. References

1. ABSTRACT

MicroRNAs (miRNAs) are a class of small noncoding RNAs involved in multiple cellular processes. Recent findings indicate that miRNA activity is globally suppressed in mouse oocytes. However, whether miRNAs are expressed and function in porcine oocytes remains unknown. In this study, our aims were to ascertain if miRNA biogenesis occurs and whether miRNA activity is globally suppressed during porcine oocyte maturation. First, to identify if miRNA biogenesis occurred, a TaqMan low-density array containing 365 mature human miRNAs was used to examine miRNA expression. This analysis revealed dynamic changes in miRNA expression, suggesting miRNA biogenesis during porcine oocyte maturation. Then, to identify if miRNA activity was globally suppressed in porcine oocytes, we focused on miR-27a, which functions in the cell cycle. miR-27a was found to facilitate the first cleavage and repress the translation of its messenger RNA target (MAP2K4), suggesting that miR-27a activity was not suppressed in porcine oocytes. Taken together, our results suggest that miRNA activity was not globally suppressed in porcine oocytes, at least for miR-27a. However, because we only investigated the activity of miR-27a, further experiments are definitely required to ascertain this point.

2. INTRODUCTION

In mammals, oocyte maturation is crucial for species survival, but the process of oocyte maturation is complex. In most mammals, fully grown oocytes are arrested in prophase of the first meiosis, which is referred to as the germinal vesicle (GV) stage and is characterized by a large nucleus (1). In response to specific signals such as gonadotropins (2), the immature oocytes resume meiosis, defined by germinal vesicle breakdown (GVBD). Then, the first polar body is extruded, and oocytes enter and arrest the cell cycle at meiosis II (MII) until they are stimulated by a sperm (3). During these processes, oocytes are transcriptionally quiescent (4); thus, post-transcriptional regulation may be important for oocyte maturation.

MicroRNAs (miRNAs) are a class of small noncoding RNAs (approximately 22 nucleotides) that regulate gene expression at a post-transcription level by translation repression or target degradation (5, 6). As newly identified regulators, miRNAs are critical for oogenesis and early embryo development. For example, miRNA expression has been shown to change dynamically during mouse oogenesis and bovine oocyte maturation, suggesting potential roles of miRNAs in these processes (7, 8). The miRNA pathway partially inhibits accumulation of the rp-S2 protein during Drosophila oocyte maturation (9). However, different studies found that miRNA function was lost in mouse oocytes and early embryos. miRNAs were shown to poorly repress their targets in mouse oocytes (10). Dgcr8-deficient oocytes, which lost canonical miRNAs, matured and fertilized normally, producing seemingly healthy offspring (11). Despite having been studied in many organisms, whether miRNAs express and function during porcine oocyte maturation remains unknown. The current study was therefore designed to investigate miRNAs in porcine oocytes.

3. MATERIALS AND METHODS

3.1. Oocyte collection and culture

All the ovaries were collected from prepubertal gilts at a local slaughterhouse and transported to the laboratory within 2 h at 37�C. The oocytes were aspirated from the antral follicles (3-6 mm in diameter) using an 18-gauge needle fixed to a 10-mL disposable syringe. Compact cumulus-oocyte complexes (COCs) were selected, washed, and cultured for up to 44 h at 39�C and 100% humidity under 5% CO2. For the maturation culture, improved TCM-199 (Gibco, Grand Island, NY, USA) supplemented with 10% (v/v) fetal bovine serum (FBS), 10% (v/v) porcine follicular fluid, 0.57 mmol/L cysteine, 10 ng/mL epidermal growth factor, 10 IU/mL pregnant mare serum gonadotropin, and 10 IU/mL human chorionic gonadotropin was used.

3.2. Evaluation of meiotic stages of porcine oocytes

Oocytes were evaluated for their nuclear status according to the procedures described by Sun et al. (12).

3.3. miRNA expression profiling assay and analysis

Seventy oocytes were lysed at 95�C for 5 min and reverse transcribed using a Multiplex RT set pool (Applied Biosystems, CA, USA). Complemetary DNA (cDNA) was then loaded onto a TaqMan low-density array (TLDA; Applied Biosystems) containing 365 miRNAs in addition to 3 small nucleolar RNAs that function as endogenous controls (RNU44, RNU48, and RNU6B) according to the manufacturer's protocol. All the miRNAs were then assayed using an ABI 7900 HT Fast Real-Time Polymerase Chain Reaction (PCR) system (Applied Biosystems). Relative miRNA expression data were analyzed using the ΔΔ Ct method. RNU6B was used as the internal control, and oocytes at 0 h were used as a calibrator. The normalized data from the assays were then subjected to hierarchical clustering using 2 methods: (1) mean-centering data for each miRNA (but not the samples) followed by correlation similarity metrics for miRNA and sample clustering, and (2) euclidean and correlation similarity metrics to cluster miRNAs and samples, respectively, without centering the data.

3.4. TaqMan real-time PCR of miRNAs

The oocytes were treated as described above. The TaqMan miRNA kit (Applied Biosystems) was used to quantify the miRNA (RNU6B, let-7g, miR-21, miR-19a, and miR-27a) expression according to the manufacturer's protocol. Relative miRNA expression data were analyzed using the ΔΔ Ct method. RNU6B was used as the internal control to determine relative miRNA expression.

3.5. Bioinformatics method

Target genes were predicted using the 3 most known target prediction databases: TargetScan 5.1 (http://www.targetscan.org), miRBase (http://microrna.sanger.ac.uk/targets/v5), and PicTar (http://pictar.bio.nyu.edu). Pathway analysis using the KEGG Orthology-based Annotation System software was conducted to determine the intersection of the targets obtained by the 3 algorithms.

3.6. Microinjection

The concentration of each type of RNA (miR-27a mimics, inhibitors, or scrambles; Ambion, Inc.) in the solutions was adjusted to 0.5 mg/mL. The microinjection was performed in 150 μL of the culture medium using microinjectors (IM-5B; Narishige, Tokyo, Japan) equipped with manipulators (MMN-1; Narishige) mounted on an inverted microscope (Diaphoto200; Nikon, Kawasaki, Japan). Approximately 50 pL of the RNA solution was injected into the cytoplasm of the noncultured COCs or MII oocytes by continuous pneumatic pressure. After microinjection, the COCs were cultured for 44 h and denuded for Western blot analysis. The MII oocytes were parthenogenetically activated as described by Wei et al. (13). Scrambles were used as the negative control in the experiment.

3.7. Western blot analysis

For detection of MAP2K4, oocytes or cell lines were dissolved in lysis buffer (80 mM b-glycerophosphate, 25 mM p-nitrophenylphosphate, 25 mM 3-Morpholinopropanesulfonic Acid (pH 7.2), 20 mM ethylene glycol tetraacetic acid, 20 mM MgCl2�6H2O, 1.1 mM Na3VO4, 1 mM phenylmethanesulfonyl fluoride, 2 μg/mL aprotinin, 200 μM leupeptin, 5 mM NaF, and 1.0% NP-40). The protein samples were denatured at 98�C for 10 min and then separated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis gel and transferred to a polyvinylidene fluoride membrane. After the transfer, the membranes were blocked with Tris-buffered saline containing 0.1% Tween-20 and 5% low-fat dry milk and probed using primary antibodies. Antibodies against MAP2K4 (Abcam, Inc., UK) and β-actin (Santa Cruz Biotechnology, Inc., USA) were used. Finally, the results were obtained using horseradish peroxidase-linked secondary antibodies (Sigma-Aldrich, USA) and enhanced chemiluminescence.

3.8. Transfection and luciferase assay

The 293FT cells were maintained in Dulbecco's Modified Eagle's Medium (Sigma-Aldrich, USA) supplemented with 10% FBS (Gibco, USA) and maintained at 37�C in the presence of 5% CO2. Pig MAP2K4 3′-untranslated region (UTR) was amplified from the cDNA of the pig ovary using primers 5′-CCCAGACACCATGTGCAATAAG-3′ and 5′-GGTGCGTTTTATTACAACTGACC-3′, cloned into the EZ-T Vector (GenStar BioSolutions, China), and then finally subcloned into the psiCHECK-2 Vector (Promega, USA). Once 70-80% of the cells were confluent, 200 ng of the reporter-gene constructs was cotransfected with 30 nM miR-27a mimics or inhibitors using Lipofectamine 2000 (Invitrogen, USA). The scrambles (Ambion, Inc., USA) were transfected as the negative control. After 24 h of transfection, firefly and Renilla luciferase activities were determined using Promega's dual luciferase assay (Progema, USA). The Renilla luciferase was normalized using the firefly luciferase for correction of transfection efficiency.

3.9. Statistical analysis

All of the data used for the statistical analysis were analyzed by the Student t-test; p < 0.05 was defined as the level of statistical significance.

4. RESULTS

4.1. Meiotic stages of porcine oocytes during in vitro maturation

Oocytes at different meiotic stages were prepared. As shown in Figure 1, at 0 h of culture, 97.4% of the oocytes were at the GV stage. After 24 h, 92.6% of the oocytes underwent GVBD. At the 36-h mark, 38.3% of the oocytes were at anaphase I or telophase I, whereas 46.7% were at MII. Finally, by 44 h, the proportion of MII oocytes increased to 77.6%.

4.2. miRNA expression profiling during porcine oocyte maturation

To determine whether miRNA biogenesis occurred during porcine oocyte maturation, oocytes at different meiotic stages were prepared according to the results shown in Figure 1. Moreover, the COCs cultured for 36 h were denuded and cultured for an additional 8 h (36 + 8 h) to evaluate whether cumulus cells affect miRNA expression during late-stage oocyte maturation. Next, these oocytes were examined for miRNA expression using the TLDA system containing 365 mature human miRNAs. As shown in Figure 2A, dynamic changes in miRNA abundance were observed during oocyte maturation. Notably, the total miRNA amount increased by approximately 60% at 24-36 h (p = 0.0007), but it was subsequently downregulated to very low levels, only 41% of the amount at 36 h (p = 0.00074; Figure 2B). In MII oocytes at 44 h, 17 miRNAs were downregulated, whereas only 1 miRNA was significantly upregulated compared with GV oocytes at 0 h (p < 0.05; Table 1). However, the oocytes that matured in 36 h with the cumulus cells and were cultured for an additional 8 h in the absence of cumulus cells showed miRNA expression profiles very similar to those of oocytes that matured continuously in 44 h in the presence of cumulus cells. The TLDA results were further confirmed using the TaqMan real-time PCR method (Figure 2C). Moreover, the pathway analysis for the potential target genes of the differentially expressed miRNAs in MII oocytes was conducted and found to be associated with the cellular processes, including the mitogen-activated protein kinase signaling pathway, cell cycle, transforming growth factor-β signaling pathway, and Wnt signaling pathway (Figure 2D). Together, these results suggest that miRNA biogenesis occurs during porcine oocyte maturation.

4.3. miR-27a facilitates the first cell cleavage

To investigate whether miRNA function was globally suppressed in porcine oocytes, we focused on miR-27a. In breast cancer cells, miR-27a has been reported to facilitate cells undergoing the G2/M transition and contribute to cell proliferation (14). Thus, if miR-27a activity was not lost in porcine oocytes, it might facilitate the first cleavage. The current study found that miR-27a expression was significantly increased in 2-cell embryos (p < 0.05) compared with MII oocytes (Figure 3A). Further investigation showed that the first cleavage rates were notably increased (p < 0.01) for MII oocytes microinjected with miR-27a mimics but were significantly decreased (p < 0.05) by the miR-27a inhibitor (Figure 3B), indicating that miR-27a facilitated the first cell cleavage.

4.4. miR-27a represses MAP2K4 translation

To further test if miR-27a repressed the translation of messenger RNA (mRNA) targets in porcine oocytes, a bioinformatics analysis was used to predict potential target genes of miR-27a. Among the predicted genes, Pig MAP2K4 3′-UTR showed a highly conserved binding site (Figure 4A). Further analysis showed that the endogenous MAP2K4 protein level increased in the 293FT cells transfected with the miR-27a inhibitor but decreased when transfected with the miR-27a mimic, indicating a correlation between miR-27a and MAP2K4 protein expression (Figure 4B). Next, MAP2K4 was demonstrated to be a direct target of miR-27a through cloning of the 3′-UTR into a luciferase reporter vector. The luciferase activity was significantly decreased in the 293FT cells transfected with the miR-27a mimic. In contrast, the miR-27a inhibitor remarkably increased luciferase activity (Figure 4C). To further assess whether miR-27a repressed MAP2K4 translation in porcine oocytes, the COCs at the GV stage were microinjected with an miR-27a mimic or inhibitor, cultured for 44 h, and then analyzed for the presence of the MAP2K4 protein. As shown in Figure 4D, the MAP2K4 protein level significantly increased in the oocytes microinjected with the miR-27a inhibitor but decreased in the oocytes microinjected with the miR-27a mimic, suggesting that miR-27a was able to repress target translation.

5. DISCUSSION

To our knowledge, this is the first study to examine miRNA biogenesis and activity during porcine oocyte maturation. We found that miRNA expression underwent dynamic changes during oocyte maturation. In addition, 18 miRNAs were found to be differentially expressed in MII oocytes (Table 1), and the predicted targets of these miRNAs were mostly related to development (Figure 2D). Moreover, by examining the effect of miR-27a on the first cleavage and its ability to repress MAP2K4 translation, we found that miR-27a activity was not suppressed in porcine oocytes. This finding is in contrast to those of recent studies on mouse oocytes and early embryos in which miRNA activity is globally suppressed (10, 11). This difference might be due to species differences. Thus, miRNA mechanisms underlying porcine oocyte development might differ from those of mouse oocytes. However, in the current study, we examined only the miR-27a activity; thus, we cannot exclude the possibility that miRNAs other than miR-27a lose their activities in porcine oocytes. Based on this aspect, further investigations are definitely required. However, it is a fact that miRNA biogenesis and changes occur during porcine oocyte maturation.

Oocyte maturation is accompanied by mRNA deadenylation and degradation. In mouse oocytes, approximately 85 pg of mRNA exists at the GV stage, but 30% of this amount undergoes degradation during oocyte maturation (15) and the protein synthesis rate also decreases after GVBD. These mRNAs are thought to be post-transcriptionally regulated, typically through motifs embedded in the 3′-UTR (16). Furthermore, miRNA is a post-transcription regulator that results in mRNA degradation or translation inhibition by binding the target site embedded in the 3′-UTR. It might be universal for miRNA to mediate clearance of maternal mRNA in the regulation of maternal to zygotic transition (17). For example, in zebrafish, miR-430 was observed to accelerate the deadenylation and clearance of maternal mRNA and was therefore critical for the maternal-to-zygotic transition (17). This mechanism might be suitable for oocyte maturation as well. Thus, in this study, the total amount of miRNA that increased in the oocytes at 36 h might be related with the deadenylation and clearance of mRNA accumulated during oocyte growth, thus resulting in the decreased rate of protein synthesis after GVBD. After mRNA is deadenylated and degraded, miRNAs might not be required; thus, in our study, the miRNA level was lowest at 44 h. Therefore, dynamic miRNA changes during porcine oocyte maturation might be related with deadenylation and degradation of maternal mRNA, thus enhancing the maturation process.

Cumulus cells are a group of granulose cells that surround the oocyte in the antral follicle and presumably participate in meiotic arrest, meiotic resumption, and cytoplasmic maturation (18). Removal of the cumulus cells before in vitro maturation was deleterious to oocyte maturation and decreased the sperm penetration rates in mice, rats, and cattle (19, 20). However, in pigs, after the COCs were cultured for 36 h, the miRNA expression profiles were similar between oocytes cultured with (44 h) and those cultured without cumulus cells (36 + 8 h; Figure 2A). It is likely that the junctions between the cumulus cells and the oocytes were disconnected after 36 h of in vitro maturation; thus, cumulus cells may not be required for late-stage oocyte maturation as far as miRNA expression is concerned. It has been reported that cumulus-oocyte communication changed dynamically during porcine oocyte maturation in vitro. Between the MI (0-24 h) and MII (36-44 h) stages, cumulus-oocyte transzonal projections became disconnected because of cumulus expansion (21). Thus, the similar expression profile in our study might also be the result of cumulus expansion. To our knowledge, this is the first study to examine the relationship between oocytes and cumulus cells from the viewpoint of miRNA expression.

In conclusion, by observing the miRNA biogenesis and the miR-27a function in the first cleavage and the repression of MAP2K4 translation, we speculate that miRNA activity might not be globally suppressed in porcine oocytes. However, further experiments are definitely required to confirm this point.

6. ACKNOWLEDGMENTS

This work is supported by the National Basic Research Program of China (Project Number: 2007CB947402) and National High Tech Research Program of China (Project Number: 2010AA10A103).

7. REFERENCES

1. Kume, S., T. Endo, Y. Nishimura, K. Kano & K. Naito: Porcine SPDYA2 (RINGO A2) stimulates CDC2 activity and accelerates meiotic maturation of porcine oocytes. Biology of reproduction, 76, 440 (2007)
doi:10.1095/biolreprod.106.057588
PMid:17151349

2. Ayalon, D., A. Tsafriri, H. R. Lindner, T. Cordova & A. Harell: Serum gonadotrophin levels in pro-oestrous rats in relation to the resumption of meiosis by the oocytes. Reproduction, 31, 51 (1972)
doi:10.1530/jrf.0.0310051

3. Sun, Q. Y. & T. Nagai: Molecular mechanisms underlying pig oocyte maturation and fertilization. Journal of Reproduction and Development, 49, 347-359 (2003)
doi:10.1262/jrd.49.347
PMid:14967910

4. Liang, C. G., Y. Q. Su, H. Y. Fan, H. Schatten & Q. Y. Sun: Mechanisms regulating oocyte meiotic resumption: roles of mitogen-activated protein kinase. Molecular Endocrinology, 21, 2037 (2007)
doi:10.1210/me.2006-0408
PMid:17536005

5. Bartel, D. P.: MicroRNAs: Genomics, Biogenesis, Mechanism, and Function. Cell, 116, 281-297 (2004)
doi:10.1016/S0092-8674(04)00045-5

6. Gaur, A., D. A. Jewell, Y. Liang, D. Ridzon, J. H. Moore, C. Chen, V. R. Ambros & M. A. Israel: Characterization of microRNA expression levels and their biological correlates in human cancer cell lines. Cancer research, 67, 2456 (2007)
doi:10.1158/0008-5472.CAN-06-2698
PMid:17363563

7. Tang, F., M. Kaneda, D. O'Carroll, P. Hajkova, S. C. Barton, Y. A. Sun, C. Lee, A. Tarakhovsky, K. Lao & M. A. Surani: Maternal microRNAs are essential for mouse zygotic development. Genes & development, 21, 644 (2007)
doi:10.1101/gad.418707

8. Tesfaye, D., D. Worku, F. Rings, C. Phatsara, E. Tholen, K. Schellander & M. Hoelker: Identification and expression profiling of microRNAs during bovine oocyte maturation using heterologous approach. Molecular reproduction and development, 76, 665-677 (2009)
doi:10.1002/mrd.21005
PMid:19170227

9. Nakahara, K., K. Kim, C. Sciulli, S. R. Dowd, J. S. Minden & R. W. Carthew: Targets of microRNA regulation in the Drosophila oocyte proteome. Proceedings of the National Academy of Sciences of the United States of America, 102, 12023 (2005)
doi:10.1073/pnas.0500053102
PMid:16099838    PMCid:1189302

10. Ma, J., M. Flemr, P. Stein, P. Berninger, R. Malik, M. Zavolan, P. Svoboda & R. M. Schultz: MicroRNA activity is suppressed in mouse oocytes. Current Biology, 20, 265-270
doi:10.1016/j.cub.2009.12.042
PMid:20116252    PMCid:2824427

11. Suh, N., L. Baehner, F. Moltzahn, C. Melton, A. Shenoy, J. Chen & R. Blelloch: MicroRNA function is globally suppressed in mouse oocytes and early embryos. Current Biology, 20, 271-277
doi:10.1016/j.cub.2009.12.044
PMid:20116247    PMCid:2872512

12. Sun, Q. Y., L. Lai, A. Bonk, R. S. Prather & H. Schatten: Cytoplasmic changes in relation to nuclear maturation and early embryo developmental potential of porcine oocytes: Effects of gonadotropins, cumulus cells, follicular size, and protein synthesis inhibition*. Molecular reproduction and development, 59, 192-198 (2001)
doi:10.1002/mrd.1022
PMid:11389554

13. Wei, H., Q. Li, J. Li, Y. Li, Y. Dai, Y. Ma, K. Xue & N. Li: Effect of leptin on oocyte maturation and subsequent pregnancy rate of cloned embryos reconstructed by somatic cell nuclear transfer in pigs. Progress in Natural Science, 18, 1583-1587 (2008)
doi:10.1016/j.pnsc.2008.05.018

14. Mertens-Talcott, S. U., S. Chintharlapalli, X. Li & S. Safe: The oncogenic microRNA-27a targets genes that regulate specificity protein transcription factors and the G2-M checkpoint in MDA-MB-231 breast cancer cells. Cancer research, 67, 11001 (2007)
doi:10.1158/0008-5472.CAN-07-2416
PMid:18006846

15. Paynton, B. V., R. Rempel & R. Bachvarova: Changes in state of adenylation and time course of degradation of maternal mRNAs during oocyte maturation and early embryonic development in the mouse* 1. Developmental biology, 129, 304-314 (1988)
doi:10.1016/0012-1606(88)90377-6

16. Evsikov, A. V., J. H. Graber, J. M. Brockman, A. Hampl, A. E. Holbrook, P. Singh, J. J. Eppig, D. Solter & B. B. Knowles: Cracking the egg: molecular dynamics and evolutionary aspects of the transition from the fully grown oocyte to embryo. Genes & development, 20, 2713 (2006)
doi:10.1101/gad.1471006

17. Giraldez, A. J., Y. Mishima, J. Rihel, R. J. Grocock, S. Van Dongen, K. Inoue, A. J. Enright & A. F. Schier: Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science, 312, 75 (2006)
doi:10.1126/science.1122689
PMid:16484454

18. Tanghe, S., A. Van Soom, H. Nauwynck, M. Coryn & A. De Kruif: Minireview: Functions of the cumulus oophorus during oocyte maturation, ovulation, and fertilization. Molecular reproduction and development, 61, 414-424 (2002)
doi:10.1002/mrd.10102
PMid:11835587

19. Vanderhyden, B. C. & D. T. Armstrong: Role of cumulus cells and serum on the in vitro maturation, fertilization, and subsequent development of rat oocytes. Biology of reproduction, 40, 720 (1989)
doi:10.1095/biolreprod40.4.720
PMid:2752073

20. Fukui, Y. & Y. Sakuma: Maturation of bovine oocytes cultured in vitro: relation to ovarian activity, follicular size and the presence or absence of cumulus cells. Biology of Reproduction, 22, 669 (1980) (doi not found)
PMid:7388112

21. Suzuki, H., B. S. Jeong & X. Yang: Dynamic changes of cumulus-oocyte cell communication during in vitro maturation of porcine oocytes. Biology of reproduction, 63, 723 (2000)
doi:10.1095/biolreprod63.3.723
PMid:10952913

Key Words: MicorRNA, Porcine Oocyte, Meiotic Maturation, MAP2K4, microRNA-27a

Send correspondence to: Ning Li, State Key Laboratory for Agrobiotechnology, China Agricultural University, Beijing 100193, China, Tel: 86-10-62733342, Fax: 86-10-62733342, E-mail: kuaileqianjin@gmail.com