• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Differential RNA Editing of Mitochondrial Genes in WA-Cytoplasmic Based Male Sterile Line Pusa 6A, and Its Maintainer and Restorer Lines

    2019-10-08 06:44:56UmakantaNgangkhamSwarupKumarParidaAshokKumarSinghTrilochanMohapatra
    Rice Science 2019年5期

    Umakanta Ngangkham, Swarup Kumar Parida, Ashok Kumar Singh, Trilochan Mohapatra

    Differential RNA Editing of Mitochondrial Genes in WA-Cytoplasmic Based Male Sterile Line Pusa 6A, and Its Maintainer and Restorer Lines

    Umakanta Ngangkham1, 2, Swarup Kumar Parida1, 3, Ashok Kumar Singh4, Trilochan Mohapatra1, 5

    ()

    RNA editing changes the nucleotides at the transcript level of mitochondrial genes which results in synthesis of functional proteins. This study was designed to find the editing sites which could be implicated in male fertility restoration and to develop editing based markers for differentiation of cytoplasmic male sterility and maintainer lines from each other. DNA and RNA from young panicles were isolated from three-line system of hybrid rice PRH10,wild abortive (WA) cytoplasm based male sterile (A linePusa6A), maintainer (B line Pusa6B) and restorer (R line PRR78) lines. Pusa6A and PRR78 having the same WA cytoplasm are allo-nuclear and iso-cytpolasmic lines. The genomic and cDNA amplicons for eight mitochondrial genes (,,,,,,and) were sequenced and compared. Differences in genomic and cDNA sequences were considered as editing. Two hundred and thirty editing sites having base substitution or insertion/deletion were identified with the highest in(5.74%) and the lowest in(0.60%). The highest editing sites were observed in fertile maintainer Pusa6B followed by PRR78 and Pusa6A,of which random five editing sites in five different rice mitochondrial transcripts namely,,,andwere chosen and validated through cleaved amplified polymorphism sequence (CAPS) analysis and found to be partially edited in four genes. The identical editing sites of different mitochondrial genes from maintainer and restorer lines might reflect their possible contribution to fertility restoration of sterile WA cytoplasm.

    cleaved amplified polymorphism sequence; cytoplasmic male sterility; hybrid rice; RNA editing; mitochondrial gene; wild abortive; nucleotide

    RNA editing is an important posttranscriptional modification process that increases the diversities of transcriptomes and proteomes in eukaryotic organisms(Wu et al, 2017). In the mitochondria of flowering plants, the primary transcripts undergo ‘editing’ that involves predominantly C to U substitutions(Covello and Gray, 1989). Besides substitutions, insertions and/or deletions of other nucleotides have also been reported (Hanson et al, 1996; Smith et al, 1997; Stuart et al, 1997). The consequences of RNA editing are either modification of the coded information for some amino acids or generation of new initiation and/or termination codons of the proteins. The first report of RNA editing was described in 1986 which involved addition or deletion of uridine residues in the RNAs of kinetoplastid, the trypanosomal mitochondria (Benne et al, 1986). Since then, it has been reported as a frequent phenomenon of both monocotyledonous and dicotyledonous higher plant species (Hiesel et al, 1989; Gray et al, 1992), abundantly in mitochondria as compared to chloroplast (Chateigner-Boutinet al, 2008). The first evidence of RNA editing in the plant mitochondria comes from comparison of the protein sequence of subunit 9 () of wheat mitochondrial ATP synthase (ATPase) and its corresponding gene sequence (Graves et al, 1990). Specific RNA editing data for plant mitochondrialmRNA are able to illustrate the basic features of RNA editing in plant mitochondria (Cattaneo, 1991). Higher degree of conservation editing sites in different plant species has indicated sharing of editing sites by both monocots and dicots. However, species specific RNA editing patterns have also been reported in several plant mitochondrial genes (Covello and Gray, 1990).

    Cytoplasmic male sterility (CMS), a maternally inherited trait, is widely observed in more than 150 plant species (Liu et al, 2007), where the male organ is sterile/non-functional but having functional female organ leading to compulsory cross-pollination. In many crops, the fertility of male sterile line is restored by fertility restorer () gene(s), which are encoded by nuclear genome. Such systems have been exploited in many crops for commercial production of hybrid seeds and to understand the nuclear-cytoplasmic organelles interaction. The best studied case of a nucleo-cytoplasmic incompatibility leading to the expression of CMS phenotype is the CMS-T of maize. In this case, mitochondrial dysfunction and male-sterility are associated with mitochondrial DNA rearrangements creating a new ORF (), which is transcribed and translated into a chimeric polypeptide (Dewey et al, 1986). Other examples are the chimeric genes-of petunia (Young and Hanson, 1987), andandof the ‘Polima’ CMS in rapeseed (Handa et al, 1995). The transcripts of all these chimeric genes, with the exception of, were found to be edited. Correlation of RNA editing in the plant mitochondrial gene with CMS has also been observed forin transgenic tobacco (Hernould et al, 1993). Transgenic tobacco plants showed male sterility when the unedited(non-functional) was introduced, while all the transgenic plants carrying the edited (functional) form ofwere found to be fertile. Cell type specific RNA editing inmitochondrial gene is involved in CMS in(Howad and Kempken, 1997). The complete editing of angene that might restore the fertility of BT-type CMS has also been reported in rice (Iwabuchi et al, 1993). Das et al(2010) and Chakraborty et al (2015) have reported that the sterility of wild abortive (WA) cytoplasm is due to competition between translated products of edited and uneditedtranscripts, which leads to impaired biogenesis of ATPase complex in rice mitochondria. Most of the mitochondrial DNA rearrangements in CMS identified in plants are associated with the genes encoding subunits of ATP synthase(Ji et al, 2014).

    The present study was carried out to gain insight into RNA editing taking place in eight mitochondrial genes from WA-cytoplasm which involves in vital mitochondrial processes and are widely used in rice hybrid breeding, and to examine possible relationship of any editing sites with expression of CMS.

    Materials and Methods

    Rice materials

    Rice materials included three lines of hybrid rice (L.), Pusa Rice Hybrid10 (PRH10) seed production viz. CMS line Pusa6A, maintainer line Pusa6B and restorer line Pusa Rice Restorer78 (PRR78). These lines were grown in the research field of the Indian Agricultural Research Institute, New Delhi, India, following standard agronomic practices.

    Isolation of total genomic DNA and cellular RNA

    Total cellular DNA was isolated from young panicles of Pusa6A, Pusa6B and PRR78 following Doyle and Doyle (1990). DNA was purified, quantified and diluted to 50ng/μL for PCR reaction. Total cellular RNA was isolated from young panicles using SV Total RNA Isolation System (www.promega.com) following the manufacturer’s protocol.

    Primer designing, PCR and RT-PCR amplification

    The primer sequences of eight rice mitochondrial genes namely,,,,,,andwere taken from Ngangkham et al (2010) (Table 1). The PCR constituents and cycling conditions with varying annealing temperature were carried out according to Ngangkham et al (2010). For reverse transcript PCR (RT-PCR), RNA was convertedto cDNA using ImProm-IITMReverse Transcription System (www.promega.com) following the manufacturer’s protocol and cDNAs were directly used for gene specific PCR amplification. The amplicons were purified using E.Z.N.A. Gel Extraction Kit (www.Omegabiotek.com).

    Cloning and sequencing

    Purified PCR and RT-PCR products were cloned into pGEM-T Easy Vector System (www.promega.com),following the manufacturer’s instruction. Recombinant white colonies were selected and used for plasmid isolation using AxyPrep Easy-96 plasmid DNA Kit (Axygen Biosciences, USA). Purified recombinant plasmids were sequenced with universal primers in both forward and reverse directions using a capillary-based Automated DNA Sequencer (MegaBACE 4000, Amersham Biosciences, USA). The vector sequences from the sequencing data were cleaned using VecScreen of NCBI online tools (http://www.ncbi.nlm.nih.gov/VecScreen/VecScreen.html).

    Detection of RNA editing sites through multiple sequence alignment and cleaved amplified polymorphism sequence (CAPS) analysis

    For discovering of RNA editing sites, the high quality sequences obtained from genomic and cDNA of each rice mitochondrial gene were aligned and compared using BioEdit Sequence Alignment Editor (http://www.mbio.ncsu.edu/bioedit/bioedit.html).RNA editing site is considered only when a specific nucleotide base alteration was detected in the cDNA sequences, as compared to the corresponding genomic sequences of a rice line. Frequency of RNA editing for each rice mitochondrial genes was also estimated. To validate the presence of editing sites identified through sequencing allignment, SNP2CAPS software tool (Thiel et al, 2004) was used for identification of appropriate restriction enzymes whose target sites were altered or introduced due to editing. The restriction enzymes predicted by SNP2CAPS were used to validate the SNP in the enzyme target sites in the mitochondrial genes. PCR amplicon of each gene was quantified on 1% agarose gel along with different known quantity of uncut lambda DNA as standard for comparison. Digestion of amplified products was carried out according to the manufacturer’s protocol and the digested products were resolved on 1.5% agarose gel using 0.5× TBE buffer. Fragment size of the digested products was determined using 100bp DNA ladder plus (Fermentas, Lithuania) and further compared with that of the undigested control to infer editing site.

    Results

    Editing in Pusa 6A, Pusa 6B and PRR78

    Comparison of the cDNA sequence of eight rice mitochondrial genes for sterile (Pusa6A), maintainer (Pusa6B) and restorer (PRR78) lines individually, with their corresponding genomic DNA sequences, led to identification of 230 nucleotides being edited, which consisted of substitution and insertion/deletion of nucleotides, from a total of 6765 nucleotides, with an overall editing frequency of 3.39% (Table 2). Maximum edited nucleotides (52 sites) were observed in, followed by(41),(36),(35),(29),(23),(8) and(6). The frequency of RNA editing ranged from 0.60% () to 5.74% () among the eight mitochondrial genes studied. Overall, there were 196 substitutions (74.24%), 37 insertions (14.02%) and 31 deletions (11.74%) in the three rice lines(Table 2).

    In the CMS line Pusa6A, there were 42 nucleotides edited (Table 2), which was the lowest among the three lines studied. Of these, 1 was deletion, 13 were insertions and 28 were substitutions. The maximum of 11 edited nucleotides were found intranscripts (1.36%), of which 3 were insertions, 1 was deletion and 7 were substitutions. This was followed by 10 edited nucleotides intranscript, 9 in, 8 in,2 inand 1 each inand, while in, there was no editing.

    Table 2. Statistics of eight mitochondrial genes in Pusa 6A, Pusa 6B and PRR78.

    NS, Number of substitution sites; NI, Number ofinsertion sites; ND, Number ofdeletion sites; NT, Number of total sites edited; F, Frequency; TPPE, Total physical position edited; TF, Total editing frequency; ‘–’, None.

    In the case of PRR78, there were 109 edited bases which included 82 substitutions, 14 insertions and 13 deletions. The maximum edited nucleotides were observed in(36), followed by(32),(14),(12),(11),(3) and(1) with no editing intranscript. When compared to iso-cytoplasmic Pusa6A line, there were 67 additional edited nucleotides in four genes viz.,,andin PRR78 (Table 2).

    In contrast to the above two parental lines carrying WA cytoplasm that has male sterilization effect, the frequency of editing was much higher in the maintainer line, Pusa6B carrying normal fertile cytoplasm. There were a total of 113 edited nucleotides distributed across all the eight genes studied. Of these, 86 were base substitutions, 10 were insertions and 17 were deletions. The majority of these edited bases were different from those in the other two rice lines. The maximum editing took place intranscript (36), followed by(24),(19),(15),(8),(6),(3) and(2). As compared to isonuclear CMS line Pusa6A, Pusa6B had 71 new edited sites distributed across the eight mitochondrial genes.

    CAPS analysis of editing sites

    The RNA editing sites detected in the transcripts of eight rice mitochondrial genes were validated through CAPS analysis using appropriate restriction enzymes as predicted by SNP2CAPS. The editing sites in eight rice mitochondrial genes were predicted to have target sites for 18 different tetra-cutter restriction enzymes. Of these, five different rice mitochondrial transcripts (,,,and) were chosen and examined through CAPS analysis.

    In case of, a 406 bp amplicon was obtained from both mitochondrial DNA and cDNA in all the three lines. When digested withI, the amplicons from genomic DNA were not affected in any of the three lines, while the cDNA got digested only in Pusa6B giving a fragment of 390 bp, but not in the other two lines. This validated the T/G substitution in theI site GATC in Pusa6B (Fig. 1-A).

    PCR amplification using primers specific togene gave 810 bp fragments for the mitochondrial DNA templates of Pusa6A, Pusa6B and PRR78, which weren’t digested byCH4IV. The same 810 bp RT-PCR product forcould not be digested in maintainer line Pusa6B. In contrast, cDNAs from Pusa6A and PRR78 were digested into fragments of 350 bp and 300 bp, and even lesser size (Fig.1-B), revealing post transcriptional creation of two or more ACGT target sites for the enzymeCH4IV including a C/A substitution in the sequenced region. This result again validated editing intranscript in Pusa6A and PRR78 but not in Pusa6B.

    In case ofgene, a 1001 bp fragment was amplified in all the three lines when mitochondrial DNA was used as a template, which could not be digested by509I. The 1001 bp fragments from cDNA of Pusa6A and PRR78 were digested in where there was no digestion for cDNA of Pusa6B (Fig.1-C). In this digestion, cDNA amplicon from Pusa6A and PRR78 gave a multiple fragments including four distinct and detectable fragments with different sizes of 701, 551, 450 and 300 bp. The combination of two digested fragment pairs such as 700 with 300 bp and 550 with 450 bp might be constituted the pre-digested amplicon. This result indicated the occurrence of two sites edited in sterile cytoplasms which wasn’t edited completely and cDNAs of Pusa6A and PRR78 consistedof both the partial edited transcripts.

    Gene specific primers forgene gave a PCR amplicon of size 1002bp in all the three lines, which was cleaved by enzymeI only in Pusa6B to give two fragments (Fig.1-D), whereas, the amplicons in Pusa6A and PRR78 remained uncut. This revealed a target site for the enzyme existing in the mtDNA only in Pusa6B. When the RT-PCR products of 1002 bp were digested using the same enzyme, digestion pattern in Pusa6B was identical to that observed for PCR product suggesting that no further change had occurred post transcriptionally. However, in case of Pusa6A and PRR78, the RT-PCR products for thegot digested to produce multiple fragments in contrast to PCR product. This indicated that multiple editing led to creation of target site CTAG for the enzymeI.

    Fig. 1. Validation of RNA editing in mitochondrialgenes by cleaved amplified polymorphism sequence analysis.

    M, 100 bp DNA ladder plus; 6A, CMS line Pusa6A; 6B, Maintainer line Pusa6B; PRR78, Restorer line PRR78.

    The recognition sites ofI (GATC) for(A),CH4IV (ACGT) for(B),509I (AATT) for(C),I (CTAG) for(D) andI (AGCT) for(E) were affected by SNPs T/G, C/A, T/A, G/C and A/T as predicted by SNP2CAPS, respectively. The digested PCR fragments were separated in 1.2% agarose gel.

    Similar to thegene patterns, the amplicon of 749 bp obtained from thegene was digested byI in Pusa6B giving two fragments of 400 bp and 349bp but not in Pusa6A and PRR78, revealing thereby presence of one target site AGCT for the enzyme in Pusa6B only. In RT-PCR using RNA templates, a 749bp amplicon was obtained in all the three lines. This fragment got digested to give multiple fragments in Pusa6A and PRR78 unlike to the PCR products, while the restriction fragment size remained identical to that from mtDNA template in Pusa6B (Fig.1-E). This result also revealed that editing took place at multiple sites in thetranscript of Pusa6A and PRR78.

    Discussion

    RNA editing, a post transcriptional process, and CMS are two important phenomena associated with higher plant mitochondria where inadequate or defective RNA editing in specific mitochondrial genes causes male sterility (Hu et al, 2013). RNA editing is reported to regulate expression of several CMS-related mitochondrial genes (Howad and Kempken, 1997). Presence of edited and unedited forms of mRNA may give rise to biologically active proteins having different physical and/or functional properties (Gray and Covello, 1993). In some cases, RNA editing is an absolute requirement, which restores critical amino acid residue(s), in absence of which non-functional proteins might get synthesized as in case of wheatmRNA (Covello and Gray, 1989).

    Effect of isonuclear-alloplasmic interaction on RNA editing

    Alloplasmic lines are created by exchanging the cytoplasm of one species line with another through substitution backcrossing while maintaining the nuclear genome, which leads to disruption of the interactions between the nucleus and cytoplasm that have co-evolved (Kianian and Kianian, 2014). It provides an opportunity to understand the interactions between the cytoplasmic organelles and the nucleus. The fertile maintainer line Pusa6B with totally different mitochondrial genome from the sterile line Pusa6A having WA sterile cytoplasm had identified 113 editing sites in eight mitochondrial genes, which is expected according to the hypothesis that decreased RNA editing frequency contributes to male sterility (Howad and Kempken, 1997; Yu et al, 2013). Presence of additional 71 (63.71%) edited nucleotides in Pusa6B as compared to Pusa6A covering all the eight mitochondrial genes might indicate occurrences of higher and different mode of post-transcriptional editing in the mitochondrial genome of Pusa6B, since it has male fertile cytoplasm compared to male sterile cytoplasm Pusa6A (Table 2). While compared to Pusa6A, the maximum additional edited nucleotides in Pusa6B were observed in, followed bywith the minimum inandalthough nuclear background were identical, and therefore the same trans-factors for editing process. Since most of the maximum additional edited nucleotides were observed in genes involved in ATP synthesis, this shows the important role of editing in vital mitochondrial processes. These results also showed that the cytoplasm may play a vital role in the editing process. Similar results are also reported by Hu et al (2013), in which five isonuclear-alloplasmic male sterile lines, when compared to their corresponding CMS lines of rice, display different editing efficiencies at the conserved sites of,andgenes.

    Effect of allonuclear-isoplasmic interaction on RNA editing

    The male sterile line Pusa6A and the restorer line PRR78 are allonuclear-isoplasmic parental lines since both having the same WA based sterile cytoplasm and thus the same mitochondrial genome constitution. Unlike in Pusa6A line, there were 67 additional edited nucleotides in four genes viz.,,andgenes in PRR78. However, there was a decrease in number of total sites edited for. Out of the total 12 edited nucleotides in the remaining three genes viz.,and, six editing sites had similar editing patterns in both iso-cytoplasmic Pusa6A and PRR78 lines, which couldn’t be correlated to male sterility. However, a total of 35 editing sites in six genes (,,,,and) observed in Pusa6A but not in PRR78, indicating that RNA editing isn’t specific to the restorer lines. However, 67 new editing sites in PRR78 compared to its iso-cytoplasmic Pusa6A, pointed to difference in nuclear-cytoplasmic interaction, which can be ascribedto the presence of fertility restorer () gene in PRR78.

    Varying degree of editing in different mitochondrial transcripts in different lines

    Comparison of the transcript sequence of eight rice mitochondrial genes in the parental lines of PRH10, viz. sterile (Pusa6A), maintainer (Pusa6B) and restorer (PRR78) lines revealed a total of 3.39% edited nucleotide, ranging from 0.60% to 5.74%, which consisted of substitutions and insertions/deletions. The frequency of RNA editing ranged from 0.60% () to 5.47% () in the eight mitochondrial genes studied, which are agreeable to those reported earlier (Giege and Brennicke, 1999; Picardi et al, 2010), and similar in-consistency of editing frequency across these eight genes as a function of gene length, which is also reported in earlier studies in,and(Mower and Palmer, 2006). Abundance of base substitution (74.2%) type of post transcriptional editing at multiple positions in transcripts of the rice mitochondria is consistent with the earlier observations in several monocotyledonous and dicotyledonous plant mitochondrial genes (Schuster et al, 1990). This result has suggested the existence of apparent coding anomalies in rice mitochondria. When we compared the RNA editing sites between Pusa6B and PRR78, there were 21 editing sites in four genes namely,(11),(6),(3) and(1), which were identical in both the fertile lines but absent in sterile line Pusa6A, for which further analysis is required to validate this editing sites with the male sterility mechanism. Since, alteration of RNA editing in CMS lines compared to other normal cytoplasm is not always necessarily the mechanism of CMS in plants (Horn et al, 2014). A relevant example in rice is that the presence ofgene, whose protein interacts with the nuclear encoded mitochondrial protein COXII, triggers premature tapetal programmed cell death and consequent pollen abortion (Luo et al, 2013) and hence male sterility in the WA CMS system, which is in contradict to Das et al (2010) who had earlier showed that an uneditedtranscript as molecular cause for this CMS system in rice and co-segregated with the male sterility. RNA editing is also influenced by many factors, such as tissue type, organ type, developmental stage, nucleo-cytoplasmic interactions and nuclear background (Grosskopf and Mulligan, 1996; Zehrmann et al, 2008) as well as by environmental stress (Liew et al, 2017). Such kinds of phenomena are also reported by Hu et al (2013),which indicated that nucleo-cytoplasmic interactions affect the RNA editing of,andtranscripts. Therefore, further analysis and characterization of these potential RNA editing sites in candidate mitochondrial genes is required to identify the mitochondrial gene responsible for male sterility.

    Partial/incomplete editing in mitochondrial transcripts

    Presence of multiple digested fragments of RT-PCR products of five mitochondrial genes namely,,,andduring the CAPS analysis, showed multiple sites of RNA editing. However, there was discrepancy between the number of editing sites revealed by sequencing and that based on CAPS analysis. For instance, in case ofgene, one edited sites corresponding to restriction enzyme target site (AATT) for restriction enzyme,509I based on sequence analysis was predicted. However, digestion of the same amplicon using the same restriction enzyme509I leads to four fragments and mixed of smaller fragments as compared to two digested fragments expected. This might be due to occurrence of incomplete/partial editing of gene transcripts and sequencing only one cloned RT-PCR product that was sequenced which do not represent the whole populations of differentially edited transcripts. An editing site of a particular gene is considered as partially/incompletely edited when the site is not edited in all the transcript pool. It is well known that partially/incompletely edited transcripts do exist in plant mitochondria (Schuster et al, 1990;Mower and Palmer, 2006; Zehrmannet al, 2008),which were represented in the RT-PCR products that were digested. While sequence comparison of genomic and multiple cDNA clones ofgene inmitochondrion, Schuster et al (1990) found none of the cDNA clones were completely edited for all potential editing sites, suggesting that such completely edited transcripts might be rare as compared to abundant of mixture transcripts edited partially. Abundance of partial editing process is further shown by sequence comparison ofmitochondrion through RNAseq which detects 76% of editing sites as partial editing, of which 28% is found to be tissue specific(Picardi et al, 2010). They further suggested that partialediting was derived from incomplete editing transcripts as well as partial edited from mixed tissue sample. Therefore, the discrepancy of editing sites in these studies might be due to failure of potential partial editing site detection because of restriction of single cDNA clone sequencing thereby precluding multiple partially edited transcripts from the transcript pool and assessment of tissue specificity of editing sites.

    Benne R, van den Burg J, Brakenhoff J P, Sloof P, van Boom JH, Tromp M C. 1986. Major transcript of the frameshiftedgene from trypanosome mitochondria contains four nucleotides that are not encoded in the DNA.,46(6): 819–826.

    Cattaneo R. 1991. Different types of messenger RNA editing., 25: 71–88.

    Chakraborty A, Mitra J, Bhattacharyya J, Pradhan S, Sikdar N, Das S, Chakraborty S, Kumar S, Lakhanpaul S, Sen S K. 2015. Transgenic expression of an unedited mitochondrialgene product from wild abortive (WA) cytoplasm of rice (L.) generates male sterility in fertile rice lines.,241:1463–1479.

    Chateigner-Boutin A L, Ramos-Vega M, Guevara-Garcia A, Andres C, de la Luz Gutierrez-Nava M, Cantero A, Delannoy E, Jiménez L F, Lurin C, Small I, León P. 2008. CLB19,a pentatricopeptide repeat protein required for editing ofandchloroplast transcripts.,56(4): 590–602.

    Covello PS, Gray MW. 1989. RNA editing in plant mitochondria.,341:662–666.

    Covello P S, Gray M W. 1990. Differences in editing at homologous sites in messenger RNAs from angiosperm mitochondria.,18(17):5189–5196.

    Das S, Sen S, Chakraborty A, Chakraborti P, Maiti M K, Basu A, Basu D, Sen S K. 2010. An unedited 1.1 kb mitochondrialgene transcript in the wild abortive cytoplasmic male sterility WA-CMS system ofL. subsp,10: 39.

    Dewey R E, Levings C S, Timothy D H. 1986. Novel recombinations in the maize mitochondrial genome produce a unique transcriptional unit in themale sterile cytoplasm.,44(3):439–449.

    Doyle J J, Doyle J L. 1990. Isolation of plant DNA from fresh tissue., 12:13–15.

    Giege P, Brennicke A. 1999. RNA editing inmitochondria effects 441 C to U changes in ORFs.,96:15324–15329.

    Graves P V, Bégu D, Velours J, Neau E, Belloc F, Litvak S, Araya A. 1990. Direct protein sequencing of wheat mitochondrial ATP-synthase subunit 9 confirms RNA editing in plants.,214(1):1–6.

    Gray M W, Hanic-Joyce P J, Covello P S. 1992. Transcription, processing and editing in plant mitochondria.,43:145–175.

    Gray M W, Covello P S. 1993. RNA editing in plant mitochondria and chloroplasts.,7:64–71.

    Grosskopf D, Mulligan RM. 1996. Developmental and tissue specificity of RNA editing in mitochondria of suspensioncultured maize cells and seedlings.,29(6):556–563.

    Handa H, Gualberto J M, Grienenberger J M. 1995. Characterization of the mitochondrialgene and its derivative,, a chimeric open reading frame specific to one mitochondrial genome theCMS in rapeseedL.,28(6):546–552.

    Hanson M R, Sutton C A, Lu B W. 1996. Plant organelle gene expression: Altered by RNA editing.,1(2):57–64.

    Hernould M, Suharsono S, Litvak S, Araya A, Mouras A. 1993. Male sterility induction in transgenic tobacco plants with an uneditedmitochondrial gene from wheat.,90(6):2370–2374.

    Hiesel R, Wissinger B, Schuster W, Brennicke A. 1989. RNA editing in plant mitochondria.,246:1632–1634.

    Horn R, Gupta K J, Colombo N. 2014. Mitochondrion role in molecular basis of cytoplasmic male sterility.,19:198–205.

    Howad W, Kempken F. 1997. Cell type-specific loss ofRNA editing in cytoplasmic male sterile,94:11090–11095.

    Hu J H, Yi R, Zhang H Y, Ding Y. 2013. Nucleo-cytoplasmic interactions affect RNA editing of,, andin alloplasmic male sterile rice (L.) lines., 13(2):87–95.

    Iwabuchi M, Kyozuka J, Shimamoto K. 1993. Processing followed by complete editing of an altered mitochondrialRNA restores fertility of cytoplasmic male-sterile rice.,12(4):1437–1446.

    Ji J J, Huang W, Li D W, Yin Y X, Chai W G, Gong Z H. 2014. A CMS-related gene,, causes increased ATP hydrolysis activity of the mitochondrial F1Fo-ATP synthase and induces male sterility in pepper (L).,32(4):888–899.

    Kianian P M A, Kianian S F. 2014. Mitochondrial dynamics and the cell cycle.,5:222.

    Liew Y J, Li Y, Baumgarten S, Voolstra C R, Aranda M. 2017. Condition-specific RNA editing in the coral symbiont.,13(2):e1006619.

    Liu Z L, Xu H, Guo J X, Liu Y G. 2007. Structural and expression variations of the mitochondrial genome conferring the WA type of CMS in rice.,49(6):908–914.

    Luo D P, Xu H, Liu Z L, Guo J X, Li H Y, Chen L T, Fang C, Zhang Q Y, Bai M, Yao N, Wu H, Wu H, Ji C H, Zheng H Q, Chen Y L, Ye S, Li X Y, Zhao X C, Li R Q, Liu Y G. 2013. A detrimental mitochondrial-nuclear interaction causes cytoplasmic male sterility in rice.,45(5):573–577.

    Mower J P, Palmer J D. 2006. Patterns of partial RNA editing in mitochondrial genes of.,276(3):285–293.

    Ngangkham U, Parida S K, De S, Kumar A R, Singh A K, Singh NK, Mohapatra T. 2010.Genic markers for wild abortive WA cytoplasm based male sterility and its fertility restoration in rice.,26:275–292.

    Picardi E, Horner DS, Chiara M, Schiavon R, Valle G, Pesole G. 2010. Large-scale detection and analysis of RNA editing in grape mtDNA by RNA deep-sequencing.,38(14): 4755–4767.

    Schuster W, Hiessel R, Wissiinger B, Brennicke A. 1990.RNA editing in the cytochrome b locus of the higher plantincludes a U to C transition.,10(5):2428–2431.

    Smith H C, Gott J M, Hanson M R. 1997.A guide to RNA editing.,3(10):1105–1123.

    Stuart K, Allen T E, Heideman S, Seiwert S D. 1997. RNA editing of kinetoplast protozoa.,61(1):105–120.

    Thiel T, Kota R, Grosse I, Stein N, Graner A. 2004. SNP2CAPS: A SNP and INDEL analysis tool for CAPS marker development.,32(1):5.

    Wu B, Chen H M, Shao J J, Zhang H, Wu K, Liu C. 2017. Identification of symmetrical RNA editing events in the mitochondria ofby strand specific RNA sequencing.,7:42250.

    Young E G, Hanson M R. 1987. A fused mitochondrial gene associated with cytoplasmic male sterility is developmentally regulated.,50(1):41–49.

    Yu J H, Zhao Y X, Qin Y T, Yue B, Zheng Y L, Xiao H L. 2013. Discovery of microRNAs associated with the S type cytoplasmic male sterility in maize.,12(2):229–238.

    Zehrmann A, van der Merwe J A, Veritskiy D, Brennicke A, Takenaka M. 2008. Seven large variations in the extent of RNA editing in plant mitochondria between three ecotypes of,8(4):319–327.

    19 July 2018;

    10 October 2018

    s:Trilochan Mohapatra(tmnrcpb@gmail.com); Umakanta Ngangkham (ukbiotech@gmail.com)

    Copyright ? 2019, China National Rice Research Institute. Hosting by Elsevier B V

    This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/)

    Peer review under responsibility of China National Rice Research Institute

    http://dx.doi.org/10.1016/j.rsci.2019.08.002

    (Managing Editor: Wang Caihong)

    桃红色精品国产亚洲av| 国产毛片a区久久久久| 日韩高清综合在线| 一夜夜www| 日韩强制内射视频| 校园人妻丝袜中文字幕| 国产高清三级在线| 成年女人永久免费观看视频| 久久精品国产自在天天线| 无遮挡黄片免费观看| 22中文网久久字幕| 久久久久国产精品人妻aⅴ院| 久久久精品欧美日韩精品| 久久精品久久久久久噜噜老黄 | 免费一级毛片在线播放高清视频| 国产精品野战在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲精品日韩av片在线观看| 日本精品一区二区三区蜜桃| 亚洲av第一区精品v没综合| 午夜视频国产福利| 狂野欧美白嫩少妇大欣赏| 18禁裸乳无遮挡免费网站照片| 一个人观看的视频www高清免费观看| 亚洲欧美激情综合另类| 亚洲,欧美,日韩| 嫩草影视91久久| 九色成人免费人妻av| 毛片女人毛片| 免费在线观看日本一区| 中文字幕av成人在线电影| АⅤ资源中文在线天堂| 国产真实乱freesex| 亚洲av成人精品一区久久| 日韩中文字幕欧美一区二区| 亚洲自偷自拍三级| 国产精品一区www在线观看 | 99在线人妻在线中文字幕| a在线观看视频网站| 亚洲精品乱码久久久v下载方式| 欧美精品啪啪一区二区三区| 国产不卡一卡二| 校园人妻丝袜中文字幕| 亚洲美女搞黄在线观看 | 午夜a级毛片| 日本爱情动作片www.在线观看 | 两性午夜刺激爽爽歪歪视频在线观看| 一级a爱片免费观看的视频| 日韩欧美免费精品| 女人被狂操c到高潮| 在线播放无遮挡| 人人妻,人人澡人人爽秒播| 久久久久久久午夜电影| 国产精品美女特级片免费视频播放器| 韩国av一区二区三区四区| 久久久午夜欧美精品| 久久精品综合一区二区三区| 精品无人区乱码1区二区| av专区在线播放| 国产高潮美女av| 99在线人妻在线中文字幕| 午夜精品久久久久久毛片777| 亚洲三级黄色毛片| 亚洲专区中文字幕在线| 18禁在线播放成人免费| 久久久国产成人精品二区| 91狼人影院| 国产三级在线视频| av视频在线观看入口| 国产v大片淫在线免费观看| 熟妇人妻久久中文字幕3abv| 欧美bdsm另类| 亚洲精品影视一区二区三区av| 亚洲精品乱码久久久v下载方式| .国产精品久久| 国产蜜桃级精品一区二区三区| 嫩草影院入口| 天堂网av新在线| 亚洲成人久久性| 少妇的逼好多水| 成人性生交大片免费视频hd| a级一级毛片免费在线观看| 看免费成人av毛片| 永久网站在线| 精品乱码久久久久久99久播| av视频在线观看入口| 韩国av在线不卡| 久久精品国产亚洲av天美| 少妇的逼水好多| 国产白丝娇喘喷水9色精品| 国产精品国产高清国产av| 亚洲成人久久性| 亚洲美女视频黄频| 亚洲狠狠婷婷综合久久图片| 性色avwww在线观看| 两人在一起打扑克的视频| 99热这里只有是精品50| 人妻制服诱惑在线中文字幕| 国产精品一区二区三区四区久久| 99久久精品热视频| 丝袜美腿在线中文| 欧美日韩国产亚洲二区| 久久精品夜夜夜夜夜久久蜜豆| 男女视频在线观看网站免费| 成人综合一区亚洲| 最好的美女福利视频网| 亚洲国产精品成人综合色| 日本在线视频免费播放| 亚洲成人久久性| 久久久久久大精品| 夜夜夜夜夜久久久久| 国产精品99久久久久久久久| 五月伊人婷婷丁香| 高清毛片免费观看视频网站| 国产av麻豆久久久久久久| 男人舔奶头视频| 偷拍熟女少妇极品色| 在线看三级毛片| 特级一级黄色大片| 亚洲综合色惰| 热99re8久久精品国产| 国产免费男女视频| 久久久久久伊人网av| 成年人黄色毛片网站| 国产在线精品亚洲第一网站| 国产精品福利在线免费观看| 日韩av在线大香蕉| 成人综合一区亚洲| 日本色播在线视频| 午夜福利视频1000在线观看| netflix在线观看网站| 日本与韩国留学比较| 狠狠狠狠99中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 国语自产精品视频在线第100页| 他把我摸到了高潮在线观看| 久久九九热精品免费| 人妻少妇偷人精品九色| www.色视频.com| 亚洲精品国产成人久久av| 舔av片在线| 午夜精品在线福利| 性欧美人与动物交配| 日日夜夜操网爽| 网址你懂的国产日韩在线| 亚洲最大成人中文| 日本成人三级电影网站| 国产精品无大码| 国产一区二区激情短视频| 日韩欧美精品免费久久| 欧美日韩亚洲国产一区二区在线观看| 成人精品一区二区免费| 成年女人看的毛片在线观看| www日本黄色视频网| 中文字幕熟女人妻在线| 别揉我奶头~嗯~啊~动态视频| 国产在视频线在精品| 免费看光身美女| 亚洲美女黄片视频| 亚洲熟妇熟女久久| 国产成人影院久久av| 九九爱精品视频在线观看| 黄色一级大片看看| 亚洲精品色激情综合| 国产精品一区www在线观看 | 欧美一区二区亚洲| 亚洲人成网站高清观看| 精品人妻视频免费看| 精品久久久久久久末码| 免费av观看视频| 少妇被粗大猛烈的视频| 成人性生交大片免费视频hd| 国产精品98久久久久久宅男小说| 国产一区二区在线av高清观看| 国产亚洲欧美98| 婷婷精品国产亚洲av| 欧美精品国产亚洲| 亚洲av电影不卡..在线观看| 1000部很黄的大片| 欧美在线一区亚洲| 日韩欧美在线乱码| 亚洲av免费在线观看| 国产中年淑女户外野战色| 国产精品国产高清国产av| 亚洲av成人av| 精品久久久久久成人av| 国产精品综合久久久久久久免费| 亚洲国产精品成人综合色| 大又大粗又爽又黄少妇毛片口| 国产真实伦视频高清在线观看 | 又粗又爽又猛毛片免费看| 日日干狠狠操夜夜爽| 国产成人a区在线观看| 免费看a级黄色片| 日韩欧美精品免费久久| 尾随美女入室| 波多野结衣高清作品| 99热这里只有是精品在线观看| 三级男女做爰猛烈吃奶摸视频| 有码 亚洲区| 特大巨黑吊av在线直播| 久久精品人妻少妇| a在线观看视频网站| 91在线观看av| 亚洲国产欧洲综合997久久,| 亚洲图色成人| 日韩精品青青久久久久久| 最近最新免费中文字幕在线| 亚洲性久久影院| 久久久久久伊人网av| 在线看三级毛片| 婷婷精品国产亚洲av| 91久久精品国产一区二区三区| 亚洲avbb在线观看| 日日啪夜夜撸| 国产欧美日韩一区二区精品| 久久久久国产精品人妻aⅴ院| 亚洲内射少妇av| 内地一区二区视频在线| 日本与韩国留学比较| 99久久精品国产国产毛片| 国产精品99久久久久久久久| 国产精品日韩av在线免费观看| 蜜桃久久精品国产亚洲av| 久久精品人妻少妇| 久久久色成人| 三级男女做爰猛烈吃奶摸视频| 中文字幕熟女人妻在线| 国产精品电影一区二区三区| avwww免费| 国产人妻一区二区三区在| 国内精品一区二区在线观看| 国产一区二区三区在线臀色熟女| 淫妇啪啪啪对白视频| 美女大奶头视频| 人妻少妇偷人精品九色| 午夜福利在线观看免费完整高清在 | 久久久精品大字幕| 久久婷婷人人爽人人干人人爱| 热99re8久久精品国产| 在线a可以看的网站| 国产精品自产拍在线观看55亚洲| 欧美最黄视频在线播放免费| www.色视频.com| 国产乱人视频| 在现免费观看毛片| 在线天堂最新版资源| 超碰av人人做人人爽久久| 精品久久久久久久久久免费视频| www.色视频.com| 免费看av在线观看网站| 成人国产一区最新在线观看| 日本欧美国产在线视频| 老熟妇乱子伦视频在线观看| 搡老妇女老女人老熟妇| 久久久久免费精品人妻一区二区| 成人特级av手机在线观看| 在线播放国产精品三级| 国内揄拍国产精品人妻在线| 在线看三级毛片| 色哟哟哟哟哟哟| 亚洲精华国产精华精| 尾随美女入室| 中文字幕熟女人妻在线| 成人精品一区二区免费| 精品久久国产蜜桃| 国产探花在线观看一区二区| 亚洲自偷自拍三级| 真实男女啪啪啪动态图| 国产主播在线观看一区二区| 色av中文字幕| 少妇人妻一区二区三区视频| 欧美丝袜亚洲另类 | 成年免费大片在线观看| 我的老师免费观看完整版| 免费搜索国产男女视频| 搞女人的毛片| 在线观看美女被高潮喷水网站| 99热精品在线国产| 在线免费观看的www视频| 亚洲avbb在线观看| 欧美bdsm另类| 色av中文字幕| 久久人人精品亚洲av| av女优亚洲男人天堂| av天堂中文字幕网| 色精品久久人妻99蜜桃| 精品福利观看| 夜夜爽天天搞| 少妇熟女aⅴ在线视频| 久久久精品欧美日韩精品| 中文字幕久久专区| 十八禁网站免费在线| 欧美国产日韩亚洲一区| 国内揄拍国产精品人妻在线| 国产精品伦人一区二区| 日韩欧美国产一区二区入口| 可以在线观看的亚洲视频| 欧美三级亚洲精品| 夜夜爽天天搞| 国产欧美日韩精品一区二区| 国产三级在线视频| 69人妻影院| 97超级碰碰碰精品色视频在线观看| 国内精品久久久久精免费| 欧美又色又爽又黄视频| 国产精品三级大全| 热99re8久久精品国产| 亚洲国产日韩欧美精品在线观看| 搡老岳熟女国产| videossex国产| 国产在线男女| 搡老岳熟女国产| 亚洲欧美清纯卡通| a级一级毛片免费在线观看| 久久精品国产自在天天线| 亚洲四区av| 91在线观看av| 丰满的人妻完整版| bbb黄色大片| 不卡视频在线观看欧美| 在线免费观看的www视频| xxxwww97欧美| 在线观看免费视频日本深夜| 色哟哟·www| 亚洲在线自拍视频| 五月玫瑰六月丁香| 精品久久久久久久久久免费视频| 亚洲av美国av| 国产 一区 欧美 日韩| 国产精品av视频在线免费观看| 蜜桃亚洲精品一区二区三区| 婷婷色综合大香蕉| 精品久久久久久成人av| 婷婷六月久久综合丁香| 久久久久久久久久久丰满 | 国内精品久久久久久久电影| av天堂在线播放| 少妇人妻精品综合一区二区 | 看片在线看免费视频| 天堂动漫精品| 美女被艹到高潮喷水动态| 99热只有精品国产| 日本一本二区三区精品| 欧美精品国产亚洲| 国产欧美日韩精品一区二区| a级毛片免费高清观看在线播放| 日日夜夜操网爽| 长腿黑丝高跟| 精品国内亚洲2022精品成人| 国产国拍精品亚洲av在线观看| 国内精品久久久久精免费| 深夜精品福利| 欧美一级a爱片免费观看看| 国产私拍福利视频在线观看| 真实男女啪啪啪动态图| 亚洲av成人av| 精品久久久久久久人妻蜜臀av| 国产一区二区在线av高清观看| 亚洲三级黄色毛片| 色综合站精品国产| 亚洲国产日韩欧美精品在线观看| 制服丝袜大香蕉在线| 色5月婷婷丁香| 嫩草影院精品99| 欧美日本视频| 观看美女的网站| av在线亚洲专区| 偷拍熟女少妇极品色| 国产伦在线观看视频一区| 狂野欧美激情性xxxx在线观看| 午夜免费男女啪啪视频观看 | 97超级碰碰碰精品色视频在线观看| 99国产极品粉嫩在线观看| 亚洲欧美激情综合另类| 国产精品福利在线免费观看| 亚洲一区高清亚洲精品| 亚洲精品成人久久久久久| 91精品国产九色| 亚洲国产精品sss在线观看| 日本黄色片子视频| 午夜福利欧美成人| 尤物成人国产欧美一区二区三区| 99久久中文字幕三级久久日本| 国产又黄又爽又无遮挡在线| 亚洲av成人av| 人人妻人人看人人澡| 久久久久久大精品| 在线观看舔阴道视频| 国产精品精品国产色婷婷| 免费一级毛片在线播放高清视频| 听说在线观看完整版免费高清| 久久久久免费精品人妻一区二区| 我要搜黄色片| 国产乱人伦免费视频| 色精品久久人妻99蜜桃| 日本-黄色视频高清免费观看| 99久久精品国产国产毛片| 国内精品美女久久久久久| 午夜精品久久久久久毛片777| 午夜福利高清视频| 麻豆av噜噜一区二区三区| 无遮挡黄片免费观看| 精品午夜福利视频在线观看一区| 99久久精品热视频| 欧美一级a爱片免费观看看| 亚洲内射少妇av| 午夜免费男女啪啪视频观看 | 国产单亲对白刺激| 婷婷六月久久综合丁香| 国产成人a区在线观看| 国产国拍精品亚洲av在线观看| 精品久久久久久久人妻蜜臀av| 久久久精品欧美日韩精品| 国产精品99久久久久久久久| 国产精品久久久久久av不卡| 亚洲人成网站在线播放欧美日韩| 久久午夜亚洲精品久久| ponron亚洲| 一进一出抽搐动态| 国产高潮美女av| 在线免费观看不下载黄p国产 | 九九爱精品视频在线观看| 久久精品国产清高在天天线| 啦啦啦韩国在线观看视频| 综合色av麻豆| 国内精品一区二区在线观看| 午夜久久久久精精品| 免费电影在线观看免费观看| 欧美最黄视频在线播放免费| 深爱激情五月婷婷| 婷婷精品国产亚洲av| 国产亚洲av嫩草精品影院| 成人毛片a级毛片在线播放| 三级毛片av免费| 午夜福利高清视频| 久久精品国产99精品国产亚洲性色| 色噜噜av男人的天堂激情| 天堂影院成人在线观看| 国产欧美日韩精品亚洲av| 中文字幕高清在线视频| av视频在线观看入口| 波野结衣二区三区在线| 久久午夜亚洲精品久久| av天堂在线播放| 亚洲av不卡在线观看| 久久久久久伊人网av| 欧美一级a爱片免费观看看| a级毛片免费高清观看在线播放| 亚洲,欧美,日韩| 性插视频无遮挡在线免费观看| 国产视频一区二区在线看| 成年人黄色毛片网站| 嫩草影视91久久| 成人二区视频| www.色视频.com| 免费在线观看日本一区| 日韩欧美免费精品| 欧美日韩综合久久久久久 | 小说图片视频综合网站| 成人性生交大片免费视频hd| 国产 一区精品| 婷婷丁香在线五月| 日韩中字成人| 国产真实乱freesex| 亚洲精品456在线播放app | 久久精品久久久久久噜噜老黄 | 亚洲精品一卡2卡三卡4卡5卡| 色哟哟·www| 国内精品久久久久精免费| 日本 av在线| 桃色一区二区三区在线观看| 人妻夜夜爽99麻豆av| 嫩草影院入口| 人妻制服诱惑在线中文字幕| or卡值多少钱| 免费观看精品视频网站| 一a级毛片在线观看| 18禁黄网站禁片午夜丰满| av福利片在线观看| 一个人免费在线观看电影| 在线观看舔阴道视频| 亚洲av第一区精品v没综合| 麻豆成人av在线观看| 美女cb高潮喷水在线观看| 国产乱人视频| 校园春色视频在线观看| 少妇丰满av| 男人狂女人下面高潮的视频| 国产精品久久电影中文字幕| 中文亚洲av片在线观看爽| 欧美一区二区亚洲| 中文字幕久久专区| 亚洲精品色激情综合| 久久久久久久亚洲中文字幕| 91麻豆精品激情在线观看国产| 在线免费观看不下载黄p国产 | 国产精品久久久久久精品电影| 午夜爱爱视频在线播放| 999久久久精品免费观看国产| 亚洲欧美清纯卡通| 在线天堂最新版资源| 最后的刺客免费高清国语| 在线观看66精品国产| 午夜福利在线观看吧| 黄色日韩在线| 岛国在线免费视频观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲最大成人av| 看黄色毛片网站| 色综合站精品国产| 人人妻人人看人人澡| 我要看日韩黄色一级片| 欧美日本视频| 亚洲av一区综合| 日本熟妇午夜| 搡老妇女老女人老熟妇| 伦理电影大哥的女人| 日本欧美国产在线视频| 少妇人妻一区二区三区视频| 精品久久久久久久久久免费视频| 亚洲av一区综合| 欧美三级亚洲精品| 一级av片app| 亚洲18禁久久av| 不卡视频在线观看欧美| 99在线视频只有这里精品首页| 97人妻精品一区二区三区麻豆| 国产高清三级在线| 国产黄色小视频在线观看| 午夜福利在线观看免费完整高清在 | 亚洲自偷自拍三级| 伊人久久精品亚洲午夜| 亚洲四区av| 美女大奶头视频| 午夜精品一区二区三区免费看| 国产一区二区亚洲精品在线观看| 黄色欧美视频在线观看| 床上黄色一级片| 国产亚洲精品av在线| 91在线精品国自产拍蜜月| 精品福利观看| 日韩欧美精品免费久久| 国产精品久久久久久亚洲av鲁大| 久久精品国产鲁丝片午夜精品 | 国产精品98久久久久久宅男小说| 在线免费观看不下载黄p国产 | 无人区码免费观看不卡| 精品人妻偷拍中文字幕| 色尼玛亚洲综合影院| 亚洲狠狠婷婷综合久久图片| 久久香蕉精品热| 欧美3d第一页| 亚洲国产欧洲综合997久久,| 日韩精品中文字幕看吧| 国产麻豆成人av免费视频| 久久精品夜夜夜夜夜久久蜜豆| 日日夜夜操网爽| 天天躁日日操中文字幕| 日本五十路高清| 久久久久久久亚洲中文字幕| 国内毛片毛片毛片毛片毛片| 日本黄大片高清| 欧美另类亚洲清纯唯美| 欧美日韩中文字幕国产精品一区二区三区| 国产乱人伦免费视频| 亚洲人成伊人成综合网2020| 国产精品免费一区二区三区在线| 久久精品国产亚洲av涩爱 | 久久久久九九精品影院| 国产 一区精品| 欧美潮喷喷水| 毛片一级片免费看久久久久 | 亚洲欧美日韩高清专用| avwww免费| 免费av观看视频| 免费观看的影片在线观看| www日本黄色视频网| 免费一级毛片在线播放高清视频| 一区二区三区免费毛片| 亚洲自拍偷在线| 免费人成在线观看视频色| 嫩草影院入口| 久久人人精品亚洲av| 联通29元200g的流量卡| 最近最新免费中文字幕在线| 亚洲乱码一区二区免费版| 免费观看精品视频网站| 最近最新免费中文字幕在线| 亚洲美女黄片视频| 成人一区二区视频在线观看| 毛片一级片免费看久久久久 | 国产精品美女特级片免费视频播放器| 色综合亚洲欧美另类图片| 国产精品99久久久久久久久| 国产精品一区www在线观看 | 人人妻人人看人人澡| 久久国内精品自在自线图片| 男人舔女人下体高潮全视频| netflix在线观看网站| 给我免费播放毛片高清在线观看| 精品久久久久久久久亚洲 | 色精品久久人妻99蜜桃| 午夜福利18| 日韩一区二区视频免费看| 亚洲va在线va天堂va国产| 黄片wwwwww| 91在线观看av| 国产淫片久久久久久久久| 久久99热这里只有精品18| 欧美绝顶高潮抽搐喷水| 狂野欧美白嫩少妇大欣赏| 我的女老师完整版在线观看| 亚洲国产精品合色在线| 国产成人av教育|