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

    Developing high-efficiency base editors by combining optimized synergistic core components with new types of nuclear localization signal peptide

    2020-06-22 05:46:00FeipengWangChengweiZhangWenXuShuangYuanJinlingSongLuLiJiuranZhaoJinxiaoYang
    The Crop Journal 2020年3期

    Feipeng Wang,Chengwei Zhang,Wen Xu,Shuang Yuan,Jinling Song,Lu Li,Jiuran Zhao*,Jinxiao Yang*

    Beijing Key Laboratory of Maize DNA Fingerprinting and Molecular Breeding, Beijing Academy of Agriculture & Forestry Sciences, Beijing 100097,China

    ABSTRACT The clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein(Cas)system has been widely used for genome editing.In this system,the cytosine base editor (CBE) and adenine base editor (ABE) allow generating precise and irreversible base mutations in a programmable manner and have been used in many different types of cells and organisms. However, their applications are limited by low editing efficiency at certain genomic target sites or at specific target cytosine(C)or adenine(A)residues.Using a strategy of combining optimized synergistic core components, we developed a new multiplex super-assembled ABE (sABE) in rice that showed higher base-editing efficiency than previously developed ABEs.We also designed a new type of nuclear localization signal(NLS)comprising a FLAG epitope tag with four copies of a codon-optimized NLS(F4NLSr2)to generate another ABE named F4NLS-sABE. This new NLS increased editing efficiency or edited additional A at several target sites.A new multiplex super-assembled CBE(sCBE)and F4NLSr2 involved F4NLS-sCBE were also created using the same strategy. F4NLS-sCBE was proven to be much more efficient than sCBE in rice.These optimized base editors will serve as powerful genome-editing tools for basic research or molecular breeding in rice and will provide a reference for the development of superior editing tools for other plants or animals.

    1.Introduction

    The clustered regularly interspaced short palindromic repeats(CRISPR)-CRISPR-associated protein(Cas)system has been used for genome editing in many different types of cells and organisms[1-4].Base editing is a newly developed genome-editing technology that permits the programmable generation of precise and irreversible base mutations without inducing double-strand breaks or requiring exogenous repair template DNA[5-7].Thus,base editing has been widely used to generate point mutations for basic research, gene therapy, or molecular breeding [8-17]. There are two base editors(BEs):the cytosine base editor(CBE),which converts target cytosineguanine(C-G)to thymine-adenine(T-A)base pairs,and the adenine base editor(ABE),which converts A-T to G-C.The most commonly used deaminase in CBEs in animals and plants is the rat cytidine deaminase rAPOBEC1 or Petromyzon marinus cytidine deaminase 1(PmCDA1) [5,6,8,11], while that in ABEs is a wild-type adenine deaminase ecTadA fused with its variant ecTadA* (ecTadA &ecTadA*) [7,18,19]. Both CBEs and ABEs comprise several essential components including sgRNA, Cas9 nickase (Cas9n), a deaminase,and nuclear localization signal (NLS) peptides. CBEs,but not ABEs,also require a uracil DNA glycosylase inhibitor(UGI).

    Although ABEs and CBEs represent a useful system for base editing,they still show low editing efficiency at certain genomic target sites or at specific target C or A residues. To improve editing efficiency, the core components have been optimized for use in animals and plants. For example, the scaffold of sgRNA was modified to obtain enhanced sgRNA (esgRNA) and the Cas9n DNA sequence was extensively codon-optimized[18,20-23]. More effective cytidine deaminases have also been developed and employed,such as the hyperactive human AID,Anc689 APOBEC, human APOBEC3A, evoAPOBEC1, and evoCDA1, while a single ecTadA* in ABE was shown [24-28] to be more effective than the conventional ecTadA&ecTadA*.For the NLS,both a codon-optimized bipartite NLS(bpNLS)and an additional FLAG epitope tag at the N terminus of NLS (FNLS)improved base-editing efficiency[23,25,29].Additional copies of the UGI increased the editing frequency of CBEs[30].

    In this study,we developed a multiplex super-assembled ABE system (sABE) by combining optimized synergistic core components, including esgRNA, a plant codon-optimized Cas9n, rice codon-optimized ecTadA&ecTadA*, and bpNLS.We found this new system was more efficient than previously developed ABEs.We further optimized sABE by using a new FLAG with 4*NLS(F4NLS) to generate F4NLS-sABE, which increased editing efficiency or edited additional A at several target sites. We also developed a multiplex super-assembled CBE system (sCBE) for rice with a rice codon-optimized rAPOBEC1 and double copies of UGI(2*UGI)as well as the same components of the sABE system(except for ecTadA & ecTadA*). The optimized F4NLS-sCBE also displayed greatly improved editing efficiency.

    2. Materials and methods

    2.1. Plasmid construction

    The nucleotide sequences of key components are described in the supplemental files. sABE was constructed based on our previous reported vector SpCas9n-pBE-basic [20]. The fragment 3 in SpCas9n-pBE-basic was replaced by the codonoptimized bpNLSr&edTadAr&ecTadA*r&Cas9np&bpNLSr&TNos fusion sequence with SnaB I and Avr II, synthesized by GenScript Corp.(Nanjing,Jiangsu,China).The OsU3 promoter was then substituted by the OsU6a promoter with Kpn I and BamH I to generate vector sABE-basic. Three synthesized fragments: FNLSr1& edTadAr& ecTadA*r& Cas9np& NLSr1&TNos, F4*NLSr2& edTadAr& ecTadA*r& Cas9np& 4*NLSr2&TNos, and edTadAr& ecTadA*r& Cas9np& 4*NLSr2& TNos,were used to replace bpNLSr&edTadAr&ecTadA*r&Cas9np&bpNLSr& TNos in sABE-basic to generate the vectors FNLSsABE-basic, F4NLS-sABE-basic, and 4NLS-sABE-basic, respectively. Following Ma et al. [31], targets were added to each basic vector to generate sABE, FNLS-sABE, F4NLS-sABE, and 4NLS-sABE for rice transformation. These four ABEs all included the target sites T1 to T5. Two plasmids were constructed for each ABE. The target sites T1, T2, and T3 were placed under the control of respectively the OsU6a,OsU3, and OsU6c promoters in one plasmid; and the target sites T4 and T5 were placed under the control of the OsU6a and OsU6c promoters in the other plasmid. For the CBE systems,the codon-optimized bpNLSr&rAPOBECr&Cas9np&2*UGIr&bpNLSr&T35s fusion sequence with SnaB I and Avr II was synthesized by GenScript and then used to replace fragment 3 in SpCas9n-pBE-basic to generate vector sCBEbasic. Three fragments composed of synthesized FNLSr1&rAPOBECr,a PCR fragment of Cas9np,and synthesized 2*UGIr&NLSr1&T35s were used to replace fragment 3 in SpCas9n-pBEbasic, leading to vector FNLS-sCBE-basic. Similarly, synthesized F4*NLSr2& rAPOBECr, a PCR fragment of Cas9npand synthesized 2*UGIr& 4*NLSr2 & T35s were combined to substitute fragment 3 to obtain vector F4NLS-sCBE-basic.Target sites T6,T7,T8,and T9 were cloned before the esgRNAs in each basic vector using Bsa I following Xie et al. [32] and finally three vectors, sCBE, FNLS-sCBE, and F4NLS-sCBE were constructed for rice transformation. Target site sequences and primers used are listed separately in Tables S1 and S2.

    2.2. Rice transformation

    The Agrobacterium tumefaciens strain EHA105 was transformed with the binary vectors using the freeze/thaw method.Embryogenic calli induced from mature seeds of rice variety Nipponbare (Oryza sativa L. japonica. cv. Nipponbare) were used for Agrobacterium-mediated rice transformation as previously [33] described. Three days after Agrobacterium infection, calli were cultured on selection medium (containing 50 μg mL?1hygromycin)for four weeks to obtain hygromycinresistant calli.Resistant calli were transferred to regeneration medium (not containing hygromycin) for shoot induction for one month. The shoots were transferred to rooting medium for root induction when they were 4-5 cm long. After two weeks,T0plants were harvested.

    2.3. DNA extraction and identification of transgenic T0 plants

    T0plants were subjected to genomic DNA extraction using a DNA-quick Plant System kit (TIANGEN Biotech, Beijing,China). The target locus was amplified by PCR with Cas9-specific primers (Table S2) and samples with an 853-bp nucleic acid fragment in agarose gel electrophoresis were identified as transgenic T0plants.

    2.4. Mutant identification

    Several transgenic T0plants were used to identify C-to-T or Ato-G conversions. Target loci were amplified with specific primers and the PCR products were purified using an EasyPure PCR Purification Kit (TransGen Biotech). The PCR products were subjected to Sanger sequencing (Tsingke Biological Technology,Beijing,China)to identify mutations.Base editing efficiency in T0plants was defined as the percentage of mutants with any target C-to-T or A-to-G substitution among all the transgenic samples. Editing efficiency of single C-to-T conversion was defined as the percentage of mutants with Cto-T substitution at a specific single position among all the transgenic samples. The primers used are listed in Table S2.

    3. Results and discussion

    3.1.Development of a multiplex super-assembled ABE system(sABE) in rice

    All of the core components of an ABE, including sgRNA, NLS,Cas9n, and adenine deaminase, affect its editing efficiency.Previous studies[18,20-22,25,29]have shown that esgRNA and codon-optimized bpNLS increase the base-editing efficiency of ABEs and CBEs in rice. Because the Cas9n DNA sequence optimized with favored human codons showed an improved base mutation frequency in human cells [23], we speculated that the editing efficiency of the adenine deaminase in the rice system might be increased by introduction of favored rice codons.Accordingly,to improve the base editing efficiency of ABE, we combined esgRNA, bpNLS, a plant-favored codonoptimized Cas9n and a rice-favored adenine deaminase to generate the super-assembled ABE for rice.

    We developed a codon-optimized bpNLS [25] for rice,designated as bpNLSr. For Cas9n, we chose the previously reported high-efficiency plant codon-optimized Cas9n(Cas9np),which has a higher GC content in the 5′terminal region like most genes in Gramineae genomes [31]. For the adenine deaminases ecTadA and ecTadA*,we used the same design as previously reported [34]. These were separately codonoptimized for rice and were designated as ecTadArand ecTadA*r, respectively. These codon-optimized adenine deaminases showed a higher rice-favored codon usage frequency than other previously used [18,29,35,36] adenine deaminases (Fig. S1). Finally, we combined esgRNA with bpNLSr, Cas9np, ecTadAr, and ecTadA*rto generate sABE(Fig.1).For greater efficiency,we designed sABE as a multiplex base-editing system in which at least three targets could be edited simultaneously(Fig.1).

    3.2. Improved base-editing efficiency of sABE

    To directly compare the performance of sABE with that of other ABEs, the previously reported T1, T2, T3, T4, and T5 target sites were selected from the genes OsNRT1.1B,OsSPL14,OsWRKY45, OsSLR1, and OsACC, respectively [18,19,36] (Table S1). OsNRT1.1B encodes a nitrogen transporter [37], OsSPL14 encodes a protein involved in plant architecture [38],OsWRKY45 encodes a critical regulator of the salicylic acid signaling pathway [39], OsSLR1 encodes a DELLA protein that serves as a repressor in the gibberellic acid signaling pathway[40], and OsACC encodes an acetyl-coenzyme A carboxylase that is essential for lipid biosynthesis[41].

    Fig. 1 - Diagram of the sABE, FNLS-sABE, F4NLS-sABE, and 4NLS-sABE base editors used in this study.

    All targets could be edited by previously developed ABEs(Table S3).However,the new sABE showed much higher base editing efficiencies at all tested targets in T0plants,except for one target that was edited with comparable efficiency(Tables 1, S4 and S5).At the target site T1, the base-editing efficiency of sABE (75%) was 1.3-fold that of the previously developed ABE(59.1%)(Tables 1,S3).Adenine residues at positions 4,6,8,and 12 in the protospacer region(counting from the 5′ end of the target sequence)were efficiently edited by sABE(Table 1),whereas the previously developed ABE could not edit the A residue at position 12,nor could it produce the mutation types of A4A6 > G4G6, A8A12 > G8G12, and A4A6A8 > G4G6G8 in T0transgenic lines (Table S3, Fig. 2). For target T2, sABE showed highly improved base-editing efficiency (70.8%), 2.7-fold that of the previously developed ABE (26.1%) (Tables 1, S3).Although both sABE and the previously developed ABE could edit A residues at positions 5, 7, and 10, the latter could not produce the mutation type A5A7A10>G5G7G10(Tables 1,S3,Fig. 3-A, B). The base-editing efficiency of sABE at target site T3 was 75% in T0plants (Table 1, Figs. S2-A, S3) and 80% in calli (data not shown),28% higher than that of the previously developed ABE (62.3%) in calli (Table S3). Interestingly, unlike at the T1 and T2 target sites,sABE edited only the A residue at position 7, whereas the previously developed ABE edited A residues at positions 4 and 7(Table S3).Among all five targets,T4 was the only one with the same base-editing frequency(12.5%) for both sABE and the previously developed ABE(Tables 1, S3, Figs. S2-B, S4). As the editing window of ABE ranged from positions 4 to 7 [7] and there was only one editable A residue at position 6 in this range,only A6>G6 was produced by both ABEs(Tables 1,S3).The editing efficiency of the previously developed ABE was 20.6% at T5, but that of sABE was 3.2-fold higher (66.7%) (Tables 1 and S3). The point mutation C2186R in OsACC endows rice with resistance to the herbicide haloxyfop [18]. C2186R corresponds to A7 > G7 in target T5.Given that the targeted A7>G7 mutation frequency of sABE (66.7%) was much higher than that of the previously developed ABE (1.3%) (Tables 1, S3, Fig. 3-C, D), sABE could be used for efficient generation of haloxyfop-resistant rice lines in future.

    3.3. Optimization of sABE using different NLS

    The base editor sABE showed higher base-editing efficiency than did a previously developed ABE at most tested sites.However, further optimization of sABE was desirable to improve its base-editing efficiency, especially at sites with low editing efficiency and sites where fewer A residues were edited by sABE. Given that various types, numbers, and locations of the NLS have been shown to influence the editing efficiency of ABE and CBE [18,23,25,29], we focused on optimizing sABE by changing the bpNLSr.

    We first selected the previously developed NLS with a FLAG epitope tag at the N terminus(FNLS)[23].Although this NLS resulted in increased editing efficiency of CBE in human and mouse cells, it has never been used to increase baseediting frequency in rice. The FNLS was codon-optimized for rice and then designated as FNLSr1. FNLSr1was then used to replace bpNLSrin sABE, generating FNLS-sABE (Fig. 1). The editing efficiency at the target sites T1 to T5 was tested with FNLS-sABE (Tables S6 and S7, Figs. S5-S9). Unexpectedly, in rice T0plants, FNLS-sABE displayed 2.1-fold and 1.2-fold higher editing efficiencies than sABE only at sites T4 and T5,while keeping comparable efficiency at target T1 (77.3% vs.75.0%) (Table 1). The editing efficiencies of FNLS-sABE at targets T2 (54.5%) and T3 (63.6%) were lower than those of sABE (70.8% for T2 and 75.0% for T3) (Table 1). Statistical results suggested that FNLSr1did not show better base editing efficiency improvement than bpNLSr(Fig.S10).FNLSr1did not perform as well as in mouse and human cells, where FNLS greatly increased C-to-T base editing [23]. This may be because the FNLSr1we used was not a suitable codonoptimized type for rice, because of the difference in genome environments in plant and animal cells, or because theperformance of FNLS differed between CBE and ABE. Whatever the reason, the results showed that not all NLS that perform well in animals will perform well in rice, and they should be thoroughly tested in each experimental system.

    Table 1-Summary of base editing efficiencies at several sites for sABE,FNLS-sABE,F4NLS-sABE and 4NLS-sABE.

    Next, we used our previously reported codon-optimized NLS for rice [34] and obtained FNLSr2. Considering that the number and locations of NLS affect the base-editing efficiency of ABE, and that four copies of N-terminal SV40-NLS resulted in 10-fold more efficient editing in knockout mouse cells[42],we combined F4*NLSr2and 4*NLSr2with sABE to generate F4NLS-sABE (Fig. 1). Given that in a previous study [18], the presence of three NLS at the C-terminus of Cas9n maximized the editing efficiency of ABE, we also constructed 4NLS-sABE with only one copy of 4*NLSr2at the C terminus of Cas9np(Fig.1). At target T4, F4NLS-sABE increased the base editing efficiency up to 38.1%, 3-fold that of sABE (12.5%) and 1.5-fold that of FNLS-sABE (26.1%) (Tables 1, S8, Fig. S11).Compared with sABE at T1 and T5 (75.0% and 66.7%efficiencies, respectively), F4NLS-sABE showed comparable base-editing efficiencies (83.3% and 61.9%), but did not edit any other A residues (Tables 1, S8 and S9, Figs. S12 and S13).Although F4NLS-sABE edited the targets T2 and T3 slightly less effectively than did sABE,the A residues at positions 12 at T2 and 4 at T3 were efficiently mutated by F4NLS-sABE but not by sABE (Tables 1, S9, Figs. S14 and S15). Thus, although F4NLS-sABE did not perform as well as we expected, it outperformed the previously developed ABEs and showed increased editing efficiency at many targets, especially those with low editing efficiency by sABE, or edited additional A.The base-editing efficiency of 4NLS-sABE was poor, and it edited a narrower range of A residues than did F4NLS-sABE and sABE (Tables 1, S10 and S11, Figs. S16-S20). This result was also unexpected,given that the version with three NLS at the C terminus of Cas9n performed best among the seven tested NLS combinations [18]. Given that two copies of three NLS, NLS with the FLAG epitope tag, and bpNLS were not tested in Li's study [18], it is unknown whether these would perform better.

    Fig.2-Targeted editing of rice genome at the T1 target site by sABE.(A)Schematic view of the T1 target site in OsNRT1.1B.The target sequence is highlighted in red.(B)Sequence chromatograms of targeted base editing at the T1 target site edited by sABE.Arrows point to positions with an edited base.

    The recently developed SurroGate systems [34], show greatly improved base-editing efficiencies of multiplex ABE or CBE in rice and the editing efficiency of F4NLS-sABE could be further optimized by use of these systems. Moreover, a simplified ABE system containing only ecTadA* instead of ecTadA & ecTadA* showed much higher editing efficiency at the T1 to T5 targets in rice[28].We speculate that the editing efficiency of F4NLS-sABE could be further improved by use of only ecTadA*and that the combination of SurroGate and only ecTadA* might perform best for optimizing F4NLS-sABE. In addition, considering the off-target effects of adenine base editors that usually occur at the RNA but not the DNA level[43-46], an F148A mutation to both ecTadA and ecTadA*,which was used to ensure a complete absence of RNA offtarget effects while maintaining efficient DNA on-target activity in human cells [47], could be tested in F4NLS-sABE for more precise editing.

    3.4. Higher editing efficiency of optimized F4NLS-sCBE in rice

    Considering the good performance of sABE and F4NLS-sABE,we hypothesized that similar architectures would also work well in multiplex CBE in rice. Because FNLS significantly increased C-to-T base editing in animals [23], we also tested FNLS in the CBE system, selecting the commonly used rAPOBEC1-based CBE system with UGI. We first codonoptimized rAPOBEC1 and UGI with rice-favored codons,generating rAPOBEC1rand UGIr(Figs.S21 and S22).rAPOBEC1r,Cas9np, and introduced two copies of UGIr(2*UGIr) together with bpNLSr,or FNLSr1and NLSr1or F4*NLSr2and 4*NLSr2into a tRNA-sgRNA system to generate three multiplex CBE systems:sCBE,FNLS-sCBE, and F4NLS-sCBE(Fig.4-A).

    Fig.3-Targeted editing of rice genome at the T2 and T5 target sites by sABE.(A)Schematic view of the T2 target site in OsSPL14.(B)Sequence chromatograms of targeted base editing at the T2 target site edited by sABE.(C)Schematic view of the T5 target site in OsAcc.(D)Sequence chromatograms of targeted base editing at the T5 target site edited by sABE.The target sequence is highlighted in red (A,C).Arrows point to positions with an edited base(B,D).T-to-C conversions in the opposite strand are shown(B,D).

    We selected four target sites, T6 to T9, from OsWaxy,OsCDC48, and OsSNB (Table S1) and combined them into one vector. In transgenic T0plants, we found that although sCBE functioned as a super-assembled system,there were still sites that could not be edited (Tables 2, S12). Of 32 T0plants, none showed C-to-T editing at T6 or T8 by sCBE and we observed very low editing efficiency at T7 and medium editing efficiency at T9 (Table 2, Figs. S23 and S24). However, F4NLSsCBE enabled efficient editing at T6 and showed 8-fold and 1.5-fold higher editing efficiency than sCBE at targets T7 and T9, respectively (Tables 2, S13, Figs. 4-B, C, S2-C, D, S25 and S26). F4NLS-sCBE showed increased editing efficiency for all the edited single C residues (Table S14, Fig. 4-D). FNLS-sCBE also performed better than sCBE (Tables 2, S15, Figs. S27 and S28),which is consistent with those obtained in animals[23].

    The editing efficiency of F4NLS-sCBE was thus superior to those of sCBE and FNLS-sCBE. The finding that F4NLS-sCBE performed better than FNLS-sCBE suggests that FNLSr2is more suitable than FNLSr1for CBE,or that additional copies of NLS can increase editing efficiency in the CBE system. In contrast, F4NLS-sABE did not show increased editing efficiency at all tested targets.This may be because of differences between the ABE and CBE systems, so that FNLSr2may be more suitable for CBE than for ABE,or additional copies of NLS might work well in CBE but not in ABE.To further improve the editing efficiency of F4NLS-sCBE,further studies should focus on using more effective cytidine deaminases, such as PmCDA1, which showed better base editing performance than rAPOBEC1 [48,49], other evolved deaminases(evoAPOBEC1, evoCDA1, Anc689 APOBEC) [25], and the Surro-Gate system together or separately.The third-generation base editor(BE3)-type CBEs induces genome-wide off-target mutations not only at the DNA but also at the RNA level in human and rice cells [43-46]. To avoid off-target effects at the RNA level, F4NLS-sCBE could be further optimized by use of engineered cytosine deaminases, such as rAPOBEC1 with W90Y + R126E point mutations and human APOBEC3A with R128A orY130F mutations [47]. CBEs' off-target effects at the DNA level have remained a challenge that invites further study.

    Fig.4-Cytosine base editing by multiplex super-assembled CBEs in rice.(A)Diagram of the sCBE,FNLS-sCBE and F4NLS-sABE base editors used in this study.(B)Schematic view of the T6 target site in OsWaxy.The target sequence is highlighted in red.(C)Sequence chromatograms of Line 22 at the OsWaxy target site edited by F4NLS-sCBE.The arrow points to the position with an edited base.G-to-A conversions in the opposite strand are shown.(D)Frequencies of single C-to-T conversions at the target site of T6,T7,T9 edited by sCBE and F4NLS-sCBE.

    Table 2-Base editing efficiencies of sCBE,FNLS-sCBE,and F4NLS-sCBE at several sites.

    4. Conclusions

    We used a strategy of combining many synergistic core components to develop new multiplex super-assembled ABE and CBE systems in rice, sABE and sCBE. The base-editing efficiencies of sABE were higher than those of previously developed ABEs.On this basis,we also designed a new type of NLS, F4NLSr2, to generate F4NLS-sABE, which enhanced editing efficiency or edited additional target A at several target sites. We further applied F4NLSr2to develop a more efficient multiplex super-assembled CBE system in rice,F4NLS-sCBE. These improved base editors will serve as powerful tools for basic research or molecular breeding in rice,and provide a reference for the development of improved genome-editing tools for other plants and animals.

    Declaration of competing interest

    The authors have submitted patent applications based on the results reported in this paper.

    Acknowledgments

    This work was supported by the Beijing Scholars Program[BSP041]. We thank Jennifer Smith, PhD, from Liwen Bianji,Edanz Group China (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.

    Author contributions

    Jinxiao Yang and Jiuran Zhao designed the experiments.Feipeng Wang constructed all vectors. Jinling Song and Lu Li performed rice transformation.Shuang Yuan performed plant DNA extraction and target site amplification. Feipeng Wang,Chengwei Zhang, and Wen Xu analyzed the results. Jinxiao Yang and Jiuran Zhao supervised the project. Jinxiao Yang,Jiuran Zhao, Chengwei Zhang, and Feipeng Wang wrote the manuscript.

    Appendix A.Supplementary data

    Supplementary data for this article can be found online at https://doi.org/10.1016/j.cj.2020.01.003.

    亚洲三级黄色毛片| 最新中文字幕久久久久| 国产精品熟女久久久久浪| 哪个播放器可以免费观看大片| 赤兔流量卡办理| 青春草国产在线视频| 丰满乱子伦码专区| 男的添女的下面高潮视频| 亚洲三级黄色毛片| 1024手机看黄色片| 波多野结衣高清无吗| 国产精品人妻久久久久久| 成人av在线播放网站| www.色视频.com| 特级一级黄色大片| 一区二区三区免费毛片| 亚洲电影在线观看av| 亚洲av不卡在线观看| 中文字幕熟女人妻在线| 亚洲av中文av极速乱| 久久久久久久国产电影| 久久6这里有精品| 国语自产精品视频在线第100页| 国产精品久久久久久精品电影小说 | 两个人的视频大全免费| 伊人久久精品亚洲午夜| 国产精品国产三级国产专区5o | 国产精品福利在线免费观看| 国产伦精品一区二区三区四那| 亚洲av电影不卡..在线观看| 中文乱码字字幕精品一区二区三区 | 国产黄色小视频在线观看| 只有这里有精品99| 欧美日韩综合久久久久久| 欧美高清性xxxxhd video| 我要看日韩黄色一级片| 汤姆久久久久久久影院中文字幕 | 1000部很黄的大片| 99久国产av精品| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品专区欧美| 国产精品国产高清国产av| 亚洲在线观看片| 成年免费大片在线观看| 亚洲真实伦在线观看| 联通29元200g的流量卡| av视频在线观看入口| 国产私拍福利视频在线观看| 亚洲真实伦在线观看| 国产v大片淫在线免费观看| 日本午夜av视频| 18+在线观看网站| 秋霞伦理黄片| 在线免费十八禁| 亚洲av一区综合| 在线免费观看的www视频| 国产午夜精品久久久久久一区二区三区| 国产老妇女一区| 特级一级黄色大片| 亚洲在线自拍视频| 超碰97精品在线观看| 毛片女人毛片| 我要搜黄色片| 噜噜噜噜噜久久久久久91| 久久久午夜欧美精品| 女人久久www免费人成看片 | 午夜福利在线观看免费完整高清在| 午夜免费男女啪啪视频观看| 色尼玛亚洲综合影院| 国产美女午夜福利| 毛片一级片免费看久久久久| 99在线人妻在线中文字幕| 日韩精品有码人妻一区| 中文字幕av成人在线电影| 日韩欧美精品v在线| 最近中文字幕2019免费版| 国产成人午夜福利电影在线观看| 久久草成人影院| 成人av在线播放网站| 六月丁香七月| 欧美性猛交╳xxx乱大交人| 国产高清视频在线观看网站| 男人舔女人下体高潮全视频| 国产成人精品久久久久久| 一个人看视频在线观看www免费| 国产探花在线观看一区二区| 男人狂女人下面高潮的视频| 国产亚洲精品久久久com| 天美传媒精品一区二区| 国产成人福利小说| 九九热线精品视视频播放| 日韩精品有码人妻一区| 国产女主播在线喷水免费视频网站 | 免费搜索国产男女视频| 日本av手机在线免费观看| 欧美成人一区二区免费高清观看| 18禁在线播放成人免费| 国产精品日韩av在线免费观看| 久久精品影院6| 日韩三级伦理在线观看| 2021天堂中文幕一二区在线观| 国产精品一区二区三区四区免费观看| 一级毛片aaaaaa免费看小| videos熟女内射| 亚洲av中文字字幕乱码综合| 黄色欧美视频在线观看| 欧美一级a爱片免费观看看| 国产成人a区在线观看| 亚洲国产成人一精品久久久| 国产人妻一区二区三区在| 毛片女人毛片| 老司机影院成人| 成人无遮挡网站| 午夜老司机福利剧场| 男插女下体视频免费在线播放| 日本免费a在线| 久久6这里有精品| 18禁在线无遮挡免费观看视频| 亚洲av一区综合| 一级毛片久久久久久久久女| 欧美性猛交╳xxx乱大交人| 超碰av人人做人人爽久久| АⅤ资源中文在线天堂| 日韩大片免费观看网站 | 亚洲av男天堂| 天堂√8在线中文| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品久久久久久精品电影| 亚洲精品,欧美精品| 天天躁夜夜躁狠狠久久av| 国产女主播在线喷水免费视频网站 | 精品无人区乱码1区二区| 一本一本综合久久| 在线免费十八禁| av女优亚洲男人天堂| 丝袜美腿在线中文| 国产精品女同一区二区软件| 精品久久国产蜜桃| 亚洲自偷自拍三级| 色综合亚洲欧美另类图片| 欧美日韩精品成人综合77777| 欧美日本亚洲视频在线播放| 久久欧美精品欧美久久欧美| 午夜福利网站1000一区二区三区| 日本免费在线观看一区| 亚洲av二区三区四区| 久久久久国产网址| 国产成人freesex在线| 色5月婷婷丁香| 日韩欧美在线乱码| 欧美日韩国产亚洲二区| 黄色一级大片看看| 亚洲av二区三区四区| 国产免费一级a男人的天堂| 久久久午夜欧美精品| 国产爱豆传媒在线观看| 久久久久久九九精品二区国产| 日本与韩国留学比较| 国产一区二区在线观看日韩| av国产免费在线观看| 热99在线观看视频| 欧美bdsm另类| 亚洲国产欧美在线一区| 国产国拍精品亚洲av在线观看| 成人特级av手机在线观看| 日本免费a在线| 国产真实乱freesex| 亚洲aⅴ乱码一区二区在线播放| 丝袜美腿在线中文| 99久国产av精品国产电影| 99在线视频只有这里精品首页| 久久精品夜色国产| 久热久热在线精品观看| 久久久国产成人精品二区| 天堂影院成人在线观看| 久久久精品大字幕| 久久久久久久久久久丰满| 国产av一区在线观看免费| 成人综合一区亚洲| 亚洲在久久综合| 熟女电影av网| 国产高潮美女av| 国产高潮美女av| 亚洲精品久久久久久婷婷小说 | 欧美又色又爽又黄视频| 人妻夜夜爽99麻豆av| 在线a可以看的网站| 国产成人freesex在线| 国产成人freesex在线| 欧美激情久久久久久爽电影| 亚洲最大成人中文| 亚洲国产精品合色在线| av福利片在线观看| 最近中文字幕2019免费版| 大话2 男鬼变身卡| 91午夜精品亚洲一区二区三区| 美女大奶头视频| 婷婷六月久久综合丁香| 精品人妻视频免费看| 91aial.com中文字幕在线观看| 国产精品一区二区三区四区久久| 国产黄色小视频在线观看| 国产精品精品国产色婷婷| 免费大片18禁| 91久久精品国产一区二区成人| 91精品伊人久久大香线蕉| 久久久久九九精品影院| 插阴视频在线观看视频| 禁无遮挡网站| 亚洲人与动物交配视频| 色视频www国产| 亚洲国产成人一精品久久久| 亚洲精品一区蜜桃| 国产精品99久久久久久久久| 久久久久性生活片| 精品久久久久久成人av| 国产精品电影一区二区三区| 岛国毛片在线播放| 欧美三级亚洲精品| 亚洲成人精品中文字幕电影| 26uuu在线亚洲综合色| 亚洲自偷自拍三级| 青青草视频在线视频观看| 成人无遮挡网站| 美女高潮的动态| 久久久a久久爽久久v久久| 一级毛片aaaaaa免费看小| 日韩国内少妇激情av| 欧美日韩一区二区视频在线观看视频在线 | 99久久人妻综合| 成人午夜精彩视频在线观看| 一个人看视频在线观看www免费| 你懂的网址亚洲精品在线观看 | 偷拍熟女少妇极品色| 亚洲综合精品二区| 国产精品国产三级国产专区5o | 2021天堂中文幕一二区在线观| 国产 一区精品| 久久亚洲精品不卡| 亚洲自偷自拍三级| 99九九线精品视频在线观看视频| av国产久精品久网站免费入址| 桃色一区二区三区在线观看| 蜜桃亚洲精品一区二区三区| 日本三级黄在线观看| 成人欧美大片| 人妻制服诱惑在线中文字幕| 国产成人精品久久久久久| av在线观看视频网站免费| 日韩制服骚丝袜av| 欧美性猛交╳xxx乱大交人| 99久久精品国产国产毛片| 在线免费观看不下载黄p国产| 特大巨黑吊av在线直播| 国产大屁股一区二区在线视频| 天堂av国产一区二区熟女人妻| 国内精品一区二区在线观看| 欧美三级亚洲精品| 2021天堂中文幕一二区在线观| 又爽又黄无遮挡网站| 成人亚洲精品av一区二区| 欧美一级a爱片免费观看看| 国产伦在线观看视频一区| 国产精品1区2区在线观看.| 一区二区三区高清视频在线| 亚洲成人中文字幕在线播放| 我的女老师完整版在线观看| 日韩强制内射视频| videossex国产| 韩国av在线不卡| 直男gayav资源| 亚洲精华国产精华液的使用体验| 欧美激情久久久久久爽电影| 纵有疾风起免费观看全集完整版 | 天堂网av新在线| 亚洲性久久影院| 成年女人看的毛片在线观看| av国产免费在线观看| 成人漫画全彩无遮挡| 亚洲av免费高清在线观看| 99久久成人亚洲精品观看| 日韩在线高清观看一区二区三区| 亚洲av日韩在线播放| 美女黄网站色视频| 六月丁香七月| 午夜福利成人在线免费观看| 91狼人影院| 欧美+日韩+精品| av播播在线观看一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 中文字幕熟女人妻在线| 亚洲婷婷狠狠爱综合网| 中文字幕制服av| 久久这里有精品视频免费| 亚洲中文字幕一区二区三区有码在线看| 日本wwww免费看| 久久综合国产亚洲精品| 国产精品熟女久久久久浪| 欧美成人午夜免费资源| 欧美人与善性xxx| 天天躁夜夜躁狠狠久久av| 少妇的逼好多水| 狠狠狠狠99中文字幕| 欧美xxxx性猛交bbbb| 国产免费又黄又爽又色| 在线a可以看的网站| 哪个播放器可以免费观看大片| 国产国拍精品亚洲av在线观看| 美女cb高潮喷水在线观看| 麻豆av噜噜一区二区三区| a级一级毛片免费在线观看| 日本黄色视频三级网站网址| 久久婷婷人人爽人人干人人爱| 国产成人a∨麻豆精品| 超碰97精品在线观看| 色尼玛亚洲综合影院| 国产伦精品一区二区三区四那| 成人毛片60女人毛片免费| 日韩制服骚丝袜av| av卡一久久| 国产色爽女视频免费观看| 亚洲经典国产精华液单| 男的添女的下面高潮视频| 精品国产三级普通话版| 一级毛片久久久久久久久女| 日本欧美国产在线视频| 春色校园在线视频观看| 美女高潮的动态| 干丝袜人妻中文字幕| 亚洲三级黄色毛片| 国产熟女欧美一区二区| 国产真实伦视频高清在线观看| 国产精品日韩av在线免费观看| 亚洲av成人精品一二三区| 免费大片18禁| 国产中年淑女户外野战色| 久久久久久久亚洲中文字幕| 国产黄片美女视频| 亚洲精品日韩av片在线观看| 综合色av麻豆| 2021少妇久久久久久久久久久| 欧美成人免费av一区二区三区| 人人妻人人看人人澡| 成人毛片60女人毛片免费| 美女黄网站色视频| 18禁在线无遮挡免费观看视频| 欧美成人免费av一区二区三区| 尾随美女入室| 亚洲av中文av极速乱| 青春草亚洲视频在线观看| 久99久视频精品免费| 亚洲av成人精品一区久久| 人体艺术视频欧美日本| 97超视频在线观看视频| 国内精品宾馆在线| 啦啦啦韩国在线观看视频| 亚洲精品aⅴ在线观看| 淫秽高清视频在线观看| 国产精品国产三级国产专区5o | 波多野结衣高清无吗| 国产精品一区二区三区四区久久| 国产亚洲精品久久久com| 免费观看精品视频网站| 日韩人妻高清精品专区| 国产精品三级大全| 国产老妇伦熟女老妇高清| av在线老鸭窝| 人体艺术视频欧美日本| 91狼人影院| av免费观看日本| 亚洲欧美日韩卡通动漫| av黄色大香蕉| 夜夜看夜夜爽夜夜摸| 亚洲自偷自拍三级| 午夜精品在线福利| 可以在线观看毛片的网站| av在线蜜桃| 国产一区二区在线av高清观看| 成人性生交大片免费视频hd| 国产又色又爽无遮挡免| 中文欧美无线码| 91在线精品国自产拍蜜月| 蜜桃亚洲精品一区二区三区| 岛国毛片在线播放| 亚洲国产精品专区欧美| 国内精品宾馆在线| 国产成人a区在线观看| 国产免费福利视频在线观看| 99在线人妻在线中文字幕| 免费在线观看成人毛片| 一个人看视频在线观看www免费| 亚洲美女视频黄频| 国产成人a∨麻豆精品| 午夜久久久久精精品| 夜夜看夜夜爽夜夜摸| 亚洲精品国产成人久久av| 国产69精品久久久久777片| 亚洲性久久影院| 久久这里只有精品中国| 亚洲精品影视一区二区三区av| 国产久久久一区二区三区| 三级男女做爰猛烈吃奶摸视频| 国产成人91sexporn| 大香蕉97超碰在线| 男女那种视频在线观看| 人妻夜夜爽99麻豆av| 国国产精品蜜臀av免费| 一级黄片播放器| 内射极品少妇av片p| av免费观看日本| 亚洲高清免费不卡视频| 久久人妻av系列| 成人毛片a级毛片在线播放| 免费大片18禁| 国产精品永久免费网站| 亚洲av日韩在线播放| 免费观看性生交大片5| 日韩视频在线欧美| 99热网站在线观看| 91精品国产九色| 亚洲电影在线观看av| 国产成年人精品一区二区| 三级毛片av免费| 亚洲,欧美,日韩| 国内揄拍国产精品人妻在线| 一本久久精品| 欧美潮喷喷水| 人人妻人人澡欧美一区二区| 成人午夜高清在线视频| 亚洲精品一区蜜桃| 身体一侧抽搐| 国产亚洲最大av| 国产精品蜜桃在线观看| 精品国产一区二区三区久久久樱花 | 99久久精品国产国产毛片| 亚洲丝袜综合中文字幕| 久久99精品国语久久久| 日日摸夜夜添夜夜爱| 国产综合懂色| 午夜福利高清视频| 成人漫画全彩无遮挡| 中文字幕亚洲精品专区| 插逼视频在线观看| 成年女人永久免费观看视频| 亚洲一区高清亚洲精品| 看免费成人av毛片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产伦在线观看视频一区| 插逼视频在线观看| 噜噜噜噜噜久久久久久91| 国产乱人偷精品视频| 久久久午夜欧美精品| 日韩欧美精品免费久久| 中文资源天堂在线| 高清日韩中文字幕在线| 青春草国产在线视频| 波野结衣二区三区在线| 亚洲国产欧洲综合997久久,| 亚洲精品456在线播放app| av.在线天堂| 中国国产av一级| 女人被狂操c到高潮| 久久久久久久久久久免费av| 亚洲av.av天堂| 在线免费观看不下载黄p国产| 网址你懂的国产日韩在线| 少妇的逼好多水| 九九爱精品视频在线观看| 三级国产精品欧美在线观看| 日韩强制内射视频| 精品99又大又爽又粗少妇毛片| 寂寞人妻少妇视频99o| 一夜夜www| 久久久精品大字幕| eeuss影院久久| 男人的好看免费观看在线视频| 日韩成人伦理影院| 国产精品一区二区三区四区久久| 久热久热在线精品观看| 一区二区三区高清视频在线| a级毛片免费高清观看在线播放| 嘟嘟电影网在线观看| 免费av不卡在线播放| or卡值多少钱| 别揉我奶头 嗯啊视频| 成人毛片a级毛片在线播放| 成人午夜高清在线视频| 亚洲在线自拍视频| 国产不卡一卡二| 51国产日韩欧美| 日本免费在线观看一区| 韩国高清视频一区二区三区| 欧美不卡视频在线免费观看| 国国产精品蜜臀av免费| 成人亚洲精品av一区二区| 男女啪啪激烈高潮av片| 丰满少妇做爰视频| 一个人看视频在线观看www免费| 精品久久国产蜜桃| 久久久久久久国产电影| 亚洲激情五月婷婷啪啪| 99久久无色码亚洲精品果冻| videos熟女内射| 国产 一区精品| 在线免费十八禁| 黄色日韩在线| 青春草亚洲视频在线观看| 国产精品麻豆人妻色哟哟久久 | 国产老妇女一区| 日韩成人av中文字幕在线观看| 国产免费又黄又爽又色| 国产精品av视频在线免费观看| 麻豆精品久久久久久蜜桃| 欧美精品一区二区大全| 欧美成人a在线观看| 搞女人的毛片| 国产亚洲午夜精品一区二区久久 | 2022亚洲国产成人精品| 综合色丁香网| 国产大屁股一区二区在线视频| 亚洲精品,欧美精品| 18禁裸乳无遮挡免费网站照片| 伦精品一区二区三区| 国产国拍精品亚洲av在线观看| 亚洲精品乱码久久久v下载方式| 国产极品天堂在线| 老司机影院成人| 全区人妻精品视频| 性插视频无遮挡在线免费观看| 看黄色毛片网站| 亚洲怡红院男人天堂| 啦啦啦啦在线视频资源| 两个人的视频大全免费| 日韩高清综合在线| 99久久九九国产精品国产免费| 少妇裸体淫交视频免费看高清| 麻豆一二三区av精品| 国产成年人精品一区二区| 亚洲图色成人| 男女下面进入的视频免费午夜| 人妻系列 视频| ponron亚洲| 春色校园在线视频观看| 内射极品少妇av片p| 国产成人精品婷婷| 国产色爽女视频免费观看| 精品一区二区免费观看| 中文字幕亚洲精品专区| 国产亚洲5aaaaa淫片| 亚洲成人久久爱视频| 人体艺术视频欧美日本| 亚洲国产精品成人综合色| 国产成人一区二区在线| 九色成人免费人妻av| 久久久久久久久久久丰满| 国产精品福利在线免费观看| 乱系列少妇在线播放| 你懂的网址亚洲精品在线观看 | 日韩欧美国产在线观看| 男人舔女人下体高潮全视频| 免费看日本二区| 亚洲精品日韩在线中文字幕| 亚洲国产精品成人综合色| 欧美激情久久久久久爽电影| 99九九线精品视频在线观看视频| 波多野结衣巨乳人妻| 长腿黑丝高跟| 欧美97在线视频| 亚洲国产欧洲综合997久久,| 国产在视频线在精品| 亚洲真实伦在线观看| 国产真实乱freesex| 两个人的视频大全免费| 大香蕉97超碰在线| 午夜福利在线观看免费完整高清在| 日韩欧美精品免费久久| 成年版毛片免费区| 日日摸夜夜添夜夜爱| 不卡视频在线观看欧美| 一级毛片久久久久久久久女| 日日摸夜夜添夜夜添av毛片| 国产精品爽爽va在线观看网站| a级一级毛片免费在线观看| 国产一区二区在线av高清观看| 欧美一级a爱片免费观看看| 国产一区亚洲一区在线观看| 国产精品麻豆人妻色哟哟久久 | 国产 一区精品| 啦啦啦啦在线视频资源| 91午夜精品亚洲一区二区三区| 久久久国产成人免费| 成年女人看的毛片在线观看| 精品欧美国产一区二区三| 人妻夜夜爽99麻豆av| 国产淫片久久久久久久久| 午夜福利在线观看免费完整高清在| 99久久精品国产国产毛片| 亚洲欧美精品自产自拍| 毛片女人毛片| 免费观看在线日韩| 亚洲四区av| 91久久精品电影网| 国产精品久久久久久精品电影小说 | 久久草成人影院| 婷婷色综合大香蕉| 成年av动漫网址| 亚洲国产精品专区欧美| 观看免费一级毛片| 国产男人的电影天堂91| 国产精品电影一区二区三区| 国产一区二区亚洲精品在线观看| 人人妻人人看人人澡| 国产精品女同一区二区软件| 少妇熟女aⅴ在线视频| 免费搜索国产男女视频|