• <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.

    av天堂在线播放| 又大又爽又粗| 99riav亚洲国产免费| 国产单亲对白刺激| 亚洲自偷自拍图片 自拍| 一进一出抽搐gif免费好疼| 国产一区二区在线观看日韩 | 国产精品亚洲美女久久久| 999久久久精品免费观看国产| 国产精品香港三级国产av潘金莲| 在线观看美女被高潮喷水网站 | 欧美不卡视频在线免费观看 | 99热这里只有精品一区 | 国产亚洲av高清不卡| 欧美日本亚洲视频在线播放| 国产一区二区三区在线臀色熟女| 国产97色在线日韩免费| 亚洲自偷自拍图片 自拍| av中文乱码字幕在线| 国产亚洲欧美在线一区二区| 可以在线观看的亚洲视频| 淫妇啪啪啪对白视频| 久久国产精品影院| 精品国产美女av久久久久小说| 人成视频在线观看免费观看| 久久久久精品国产欧美久久久| 又黄又爽又免费观看的视频| 久久这里只有精品中国| 又大又爽又粗| 久久午夜综合久久蜜桃| 色哟哟哟哟哟哟| 美女黄网站色视频| 国产真实乱freesex| 99热这里只有是精品50| ponron亚洲| 变态另类丝袜制服| 99久久国产精品久久久| 香蕉久久夜色| 啦啦啦免费观看视频1| 在线观看美女被高潮喷水网站 | 夜夜看夜夜爽夜夜摸| 老汉色∧v一级毛片| 久久久久久免费高清国产稀缺| 伊人久久大香线蕉亚洲五| 999精品在线视频| 999精品在线视频| 三级国产精品欧美在线观看 | 日本黄大片高清| 久久婷婷成人综合色麻豆| 两个人的视频大全免费| 天堂av国产一区二区熟女人妻 | 手机成人av网站| 全区人妻精品视频| 久久精品91蜜桃| 亚洲一区二区三区色噜噜| 在线十欧美十亚洲十日本专区| 国产麻豆成人av免费视频| 禁无遮挡网站| 十八禁人妻一区二区| 久久这里只有精品19| 老司机福利观看| 国产99久久九九免费精品| 久久久久免费精品人妻一区二区| 久久国产精品影院| 午夜免费成人在线视频| 亚洲男人天堂网一区| 香蕉丝袜av| tocl精华| 亚洲一区高清亚洲精品| 日本撒尿小便嘘嘘汇集6| 老司机福利观看| 一级毛片精品| 亚洲av成人一区二区三| 欧美zozozo另类| 久久精品91无色码中文字幕| 精品熟女少妇八av免费久了| av超薄肉色丝袜交足视频| 亚洲黑人精品在线| www.www免费av| 精品无人区乱码1区二区| 精品第一国产精品| 国产熟女xx| 日韩中文字幕欧美一区二区| 五月伊人婷婷丁香| 男女床上黄色一级片免费看| 51午夜福利影视在线观看| 国产熟女午夜一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 欧美最黄视频在线播放免费| 激情在线观看视频在线高清| 男人舔奶头视频| 精品乱码久久久久久99久播| 最近最新中文字幕大全免费视频| 国产人伦9x9x在线观看| 国产午夜精品久久久久久| 男插女下体视频免费在线播放| 色综合欧美亚洲国产小说| 无限看片的www在线观看| 天天躁夜夜躁狠狠躁躁| 久久中文字幕人妻熟女| 亚洲精品国产一区二区精华液| 亚洲国产精品合色在线| 啦啦啦免费观看视频1| 亚洲精品中文字幕一二三四区| e午夜精品久久久久久久| 午夜日韩欧美国产| 国产在线观看jvid| 欧美zozozo另类| 亚洲第一欧美日韩一区二区三区| 熟妇人妻久久中文字幕3abv| 亚洲成人中文字幕在线播放| 精品久久久久久久末码| 国产精品一区二区免费欧美| 亚洲一区二区三区色噜噜| 国产成年人精品一区二区| 日韩大码丰满熟妇| 日韩大尺度精品在线看网址| 999久久久精品免费观看国产| 国产91精品成人一区二区三区| 国产精品久久久久久人妻精品电影| 欧美黄色淫秽网站| 757午夜福利合集在线观看| 欧美在线黄色| 免费人成视频x8x8入口观看| 两个人看的免费小视频| 999久久久国产精品视频| 琪琪午夜伦伦电影理论片6080| 欧美zozozo另类| 桃色一区二区三区在线观看| 国语自产精品视频在线第100页| 国产精品免费一区二区三区在线| 一本一本综合久久| 国产亚洲精品av在线| 999久久久国产精品视频| 欧美日韩中文字幕国产精品一区二区三区| 两人在一起打扑克的视频| 久久精品国产综合久久久| 日本五十路高清| 成在线人永久免费视频| 久久久国产成人精品二区| 国产精品国产高清国产av| 欧美午夜高清在线| 精品熟女少妇八av免费久了| 成年版毛片免费区| 女生性感内裤真人,穿戴方法视频| 99久久精品国产亚洲精品| 国产精品国产高清国产av| 美女黄网站色视频| 狂野欧美激情性xxxx| 又紧又爽又黄一区二区| 最近视频中文字幕2019在线8| 精品午夜福利视频在线观看一区| 久久婷婷成人综合色麻豆| 90打野战视频偷拍视频| 中文在线观看免费www的网站 | 十八禁人妻一区二区| 免费一级毛片在线播放高清视频| 久久亚洲精品不卡| 99久久国产精品久久久| 在线观看免费视频日本深夜| 俄罗斯特黄特色一大片| 最近最新中文字幕大全电影3| 欧美一区二区精品小视频在线| 我要搜黄色片| 啦啦啦观看免费观看视频高清| 国产av在哪里看| 欧美zozozo另类| 久久久精品大字幕| 啦啦啦免费观看视频1| 毛片女人毛片| 国产午夜精品论理片| 久久婷婷人人爽人人干人人爱| 午夜精品在线福利| 日本熟妇午夜| 99在线人妻在线中文字幕| 国产69精品久久久久777片 | 一二三四在线观看免费中文在| 99国产极品粉嫩在线观看| 久久人人精品亚洲av| 久久精品国产清高在天天线| 十八禁人妻一区二区| 美女扒开内裤让男人捅视频| 亚洲av中文字字幕乱码综合| 国产午夜福利久久久久久| 亚洲精品久久成人aⅴ小说| 国产精品,欧美在线| 日本精品一区二区三区蜜桃| 欧美成人一区二区免费高清观看 | 蜜桃久久精品国产亚洲av| 国产免费男女视频| 最好的美女福利视频网| 成人特级av手机在线观看| 国产片特级美女逼逼视频| 又黄又爽又刺激的免费视频.| 熟女电影av网| 亚洲av.av天堂| 一区二区三区免费毛片| 亚洲图色成人| 亚洲成人av在线免费| 国产精品久久久久久精品电影| 观看美女的网站| 一级黄片播放器| 中文在线观看免费www的网站| 黄色配什么色好看| 日日啪夜夜撸| 久久精品国产清高在天天线| 日本av手机在线免费观看| 午夜爱爱视频在线播放| 国产综合懂色| 国产亚洲精品久久久久久毛片| 成人三级黄色视频| 我要看日韩黄色一级片| 乱系列少妇在线播放| 亚洲人成网站在线播| 看免费成人av毛片| 黄片无遮挡物在线观看| 久久久久久国产a免费观看| 男插女下体视频免费在线播放| 蜜桃亚洲精品一区二区三区| 日韩成人伦理影院| 黄色配什么色好看| 亚洲中文字幕日韩| 99热网站在线观看| 毛片一级片免费看久久久久| 成人二区视频| 欧美潮喷喷水| 99久久九九国产精品国产免费| 两个人的视频大全免费| 国产精品一区二区性色av| 亚洲av成人av| 久久久国产成人免费| 国产精品野战在线观看| 日本免费一区二区三区高清不卡| 国产国拍精品亚洲av在线观看| 麻豆国产av国片精品| 国产伦精品一区二区三区四那| 又粗又爽又猛毛片免费看| 久久久久久大精品| 久久人人爽人人片av| 日日摸夜夜添夜夜添av毛片| 最好的美女福利视频网| 亚洲aⅴ乱码一区二区在线播放| 精品人妻一区二区三区麻豆| 天堂中文最新版在线下载 | 青春草亚洲视频在线观看| 国产91av在线免费观看| 黄片无遮挡物在线观看| 中文字幕精品亚洲无线码一区| 国产老妇女一区| 一个人免费在线观看电影| 精品久久久噜噜| 精品久久久久久久久亚洲| 欧美xxxx性猛交bbbb| 熟女电影av网| 美女xxoo啪啪120秒动态图| 久久久国产成人精品二区| 99久久久亚洲精品蜜臀av| 精品久久久久久久末码| 白带黄色成豆腐渣| 最近中文字幕高清免费大全6| 日韩大尺度精品在线看网址| 丰满人妻一区二区三区视频av| 久久草成人影院| 一进一出抽搐gif免费好疼| 欧美xxxx黑人xx丫x性爽| 国产伦理片在线播放av一区 | 日本三级黄在线观看| 神马国产精品三级电影在线观看| 亚洲成人久久爱视频| 日本成人三级电影网站| 99久久无色码亚洲精品果冻| 久久久久免费精品人妻一区二区| 变态另类丝袜制服| 日韩大尺度精品在线看网址| 亚洲精品色激情综合| 成熟少妇高潮喷水视频| 亚洲精品久久久久久婷婷小说 | 日韩一区二区三区影片| 天天躁夜夜躁狠狠久久av| 一本久久中文字幕| 精品不卡国产一区二区三区| av卡一久久| 精品久久久久久久久久免费视频| 三级国产精品欧美在线观看| 亚洲精品色激情综合| 一级黄片播放器| 亚洲av不卡在线观看| 美女脱内裤让男人舔精品视频 | 好男人在线观看高清免费视频| 国产精品电影一区二区三区| 久久热精品热| 少妇裸体淫交视频免费看高清| 欧美极品一区二区三区四区| 国产亚洲精品久久久久久毛片| 国产一区二区在线av高清观看| 久久精品久久久久久久性| 黄色一级大片看看| 免费不卡的大黄色大毛片视频在线观看 | 如何舔出高潮| 国产 一区精品| 美女 人体艺术 gogo| 日本黄大片高清| 午夜老司机福利剧场| 伦理电影大哥的女人| 在线观看美女被高潮喷水网站| 欧美变态另类bdsm刘玥| 成年女人永久免费观看视频| 国产高清三级在线| 亚洲激情五月婷婷啪啪| 综合色av麻豆| 日韩一区二区三区影片| 天堂av国产一区二区熟女人妻| 久久这里只有精品中国| 麻豆av噜噜一区二区三区| 永久网站在线| 日本av手机在线免费观看| 免费一级毛片在线播放高清视频| 欧美bdsm另类| av在线播放精品| 99精品在免费线老司机午夜| 日本与韩国留学比较| 午夜精品在线福利| 一个人看视频在线观看www免费| 中文欧美无线码| 一区二区三区四区激情视频 | 精品人妻视频免费看| videossex国产| 美女被艹到高潮喷水动态| 黄片wwwwww| 男人和女人高潮做爰伦理| av在线亚洲专区| 久久精品国产清高在天天线| 少妇人妻精品综合一区二区 | 国产精品,欧美在线| 亚洲人与动物交配视频| 国产精品99久久久久久久久| 色吧在线观看| 久久精品国产99精品国产亚洲性色| 人体艺术视频欧美日本| a级毛片a级免费在线| 成年av动漫网址| 国产爱豆传媒在线观看| 亚洲欧美中文字幕日韩二区| 22中文网久久字幕| 狂野欧美白嫩少妇大欣赏| 悠悠久久av| 性色avwww在线观看| 三级经典国产精品| 97热精品久久久久久| h日本视频在线播放| 最近中文字幕高清免费大全6| 麻豆成人av视频| 一级二级三级毛片免费看| 国产精品人妻久久久久久| 久久久久久久午夜电影| 久久久久久久久中文| 丰满乱子伦码专区| 中文字幕av在线有码专区| 久久久久久国产a免费观看| 国产精品电影一区二区三区| 久久久久久久午夜电影| 在线播放无遮挡| 中文字幕熟女人妻在线| 亚洲欧美成人精品一区二区| 在线观看美女被高潮喷水网站| 尾随美女入室| 嫩草影院新地址| 日本熟妇午夜| 一个人观看的视频www高清免费观看| 国产成人freesex在线| 国产黄色小视频在线观看| 欧美+日韩+精品| 日本免费a在线| av在线亚洲专区| 免费大片18禁| 看免费成人av毛片| 亚洲五月天丁香| av天堂中文字幕网| 久久久久九九精品影院| 日本-黄色视频高清免费观看| 男插女下体视频免费在线播放| 中文亚洲av片在线观看爽| 午夜免费激情av| 亚洲国产高清在线一区二区三| 午夜老司机福利剧场| 夫妻性生交免费视频一级片| 美女cb高潮喷水在线观看| 国产亚洲精品久久久久久毛片| 一本久久精品| 99久久久亚洲精品蜜臀av| 国产人妻一区二区三区在| 日韩av不卡免费在线播放| 51国产日韩欧美| 国产精品久久久久久亚洲av鲁大| 色视频www国产| 欧美日本亚洲视频在线播放| 欧美成人精品欧美一级黄| 黄片wwwwww| 亚洲美女搞黄在线观看| 亚洲av男天堂| 大又大粗又爽又黄少妇毛片口| 国产黄a三级三级三级人| a级一级毛片免费在线观看| 久久精品综合一区二区三区| 国产成人一区二区在线| 亚洲av第一区精品v没综合| 两性午夜刺激爽爽歪歪视频在线观看| 久久韩国三级中文字幕| 亚洲中文字幕日韩| 国产人妻一区二区三区在| 亚洲精品乱码久久久久久按摩| 久久人人爽人人片av| 亚洲乱码一区二区免费版| 又爽又黄无遮挡网站| 国产真实伦视频高清在线观看| 特大巨黑吊av在线直播| 亚洲欧洲国产日韩| 亚洲人成网站高清观看| 成人特级av手机在线观看| 国产精品不卡视频一区二区| 亚洲欧美日韩卡通动漫| 99视频精品全部免费 在线| 亚洲高清免费不卡视频| 天天躁日日操中文字幕| 亚洲精品久久国产高清桃花| 精品一区二区三区人妻视频| 最近的中文字幕免费完整| 2021天堂中文幕一二区在线观| 又爽又黄无遮挡网站| 成人国产麻豆网| 最近最新中文字幕大全电影3| 天堂√8在线中文| 国产成人午夜福利电影在线观看| 91在线精品国自产拍蜜月| 国产精品一区二区在线观看99 | 日产精品乱码卡一卡2卡三| 蜜桃亚洲精品一区二区三区| 亚洲精品久久国产高清桃花| 免费av不卡在线播放| 亚洲四区av| 国产成人精品久久久久久| 久久久久九九精品影院| 国产精品一区二区在线观看99 | 18+在线观看网站| 亚洲一级一片aⅴ在线观看| a级毛色黄片| 啦啦啦韩国在线观看视频| 国产黄片美女视频| 成人漫画全彩无遮挡| 免费人成在线观看视频色| 日韩欧美一区二区三区在线观看| 给我免费播放毛片高清在线观看| 91狼人影院| 亚洲自偷自拍三级| 内地一区二区视频在线| 深夜a级毛片| 一级毛片我不卡| h日本视频在线播放| 国产黄片视频在线免费观看| 欧美日本视频| 久久久久久大精品| 国产av麻豆久久久久久久| 五月玫瑰六月丁香| 亚洲久久久久久中文字幕| 日本免费一区二区三区高清不卡| 日本欧美国产在线视频| 网址你懂的国产日韩在线| 国产高清有码在线观看视频| 久久久久免费精品人妻一区二区| 91久久精品国产一区二区三区| 啦啦啦韩国在线观看视频| 国产精品一区二区三区四区久久| 精品久久久久久久久久久久久| 成人国产麻豆网| 中文字幕免费在线视频6| 婷婷色综合大香蕉| 中出人妻视频一区二区| 久久久久久伊人网av| 国产激情偷乱视频一区二区| 成人午夜高清在线视频| 免费观看在线日韩| 国产视频内射| 日日撸夜夜添| 亚洲av中文字字幕乱码综合| 亚洲aⅴ乱码一区二区在线播放| 精品熟女少妇av免费看| 麻豆乱淫一区二区| 国产蜜桃级精品一区二区三区| 亚洲丝袜综合中文字幕| 一区二区三区四区激情视频 | 久久久国产成人免费| 黄片无遮挡物在线观看| 老女人水多毛片| 淫秽高清视频在线观看| 婷婷精品国产亚洲av| 在线免费十八禁| 黄片无遮挡物在线观看| 老女人水多毛片| or卡值多少钱| 国产黄片美女视频| 一级毛片我不卡| 国产三级中文精品| 别揉我奶头 嗯啊视频| 亚洲精品粉嫩美女一区| 国内精品一区二区在线观看| 好男人在线观看高清免费视频| 22中文网久久字幕| 亚洲久久久久久中文字幕| 变态另类丝袜制服| 大香蕉久久网| 久99久视频精品免费| 欧美bdsm另类| 美女被艹到高潮喷水动态| av视频在线观看入口| 人妻夜夜爽99麻豆av| 久久精品国产亚洲网站| 国产精品久久久久久精品电影| 赤兔流量卡办理| av在线播放精品| 免费观看在线日韩| 久久精品国产99精品国产亚洲性色| 免费电影在线观看免费观看| 久久久久免费精品人妻一区二区| 一级毛片我不卡| 亚洲真实伦在线观看| 级片在线观看| 激情 狠狠 欧美| 国产精品一区二区在线观看99 | 99九九线精品视频在线观看视频| 九九在线视频观看精品| 少妇熟女欧美另类| av免费在线看不卡| 麻豆成人av视频| 毛片女人毛片| 哪个播放器可以免费观看大片| 成人美女网站在线观看视频| 亚洲av第一区精品v没综合| 国产成人精品久久久久久| 免费看日本二区| 中文精品一卡2卡3卡4更新| 日韩国内少妇激情av| 99久国产av精品国产电影| 国产一区二区三区av在线 | 国产一区二区亚洲精品在线观看| 日韩国内少妇激情av| 久久久久网色| 婷婷色综合大香蕉| 99久久精品一区二区三区| 欧美zozozo另类| 久久热精品热| 在线播放国产精品三级| 国产精品爽爽va在线观看网站| 日本成人三级电影网站| 久久精品夜色国产| 99热网站在线观看| 国产伦在线观看视频一区| 午夜精品在线福利| 欧美bdsm另类| 国产在线男女| 日韩高清综合在线| 久久99精品国语久久久| 国产 一区 欧美 日韩| 国产视频首页在线观看| 色综合色国产| 非洲黑人性xxxx精品又粗又长| 久久精品国产清高在天天线| 亚洲国产日韩欧美精品在线观看| 最近手机中文字幕大全| av天堂在线播放| 男女那种视频在线观看| 我要看日韩黄色一级片| 国产熟女欧美一区二区| 欧美色欧美亚洲另类二区| 国产精品国产三级国产av玫瑰| 两性午夜刺激爽爽歪歪视频在线观看| 一本一本综合久久| avwww免费| 亚洲人成网站在线观看播放| 亚洲精品国产av成人精品| 国产午夜精品久久久久久一区二区三区| 亚洲在线观看片| 免费人成视频x8x8入口观看| 欧美又色又爽又黄视频| 日韩欧美国产在线观看| 免费大片18禁| 男的添女的下面高潮视频| 国产av一区在线观看免费| 亚洲av一区综合| 国产免费一级a男人的天堂| 亚洲国产精品久久男人天堂| 18禁在线播放成人免费| 12—13女人毛片做爰片一| 国内久久婷婷六月综合欲色啪| 免费看美女性在线毛片视频| 国产高潮美女av| 99久久成人亚洲精品观看| 深爱激情五月婷婷| 国产高潮美女av| av黄色大香蕉| 91在线精品国自产拍蜜月| 亚洲精品自拍成人| 亚洲人成网站高清观看| 久久韩国三级中文字幕| 最好的美女福利视频网| 午夜精品一区二区三区免费看| av卡一久久| 亚洲欧美精品专区久久| 啦啦啦观看免费观看视频高清| 精品久久久久久久末码| 亚洲乱码一区二区免费版| 99精品在免费线老司机午夜| 最后的刺客免费高清国语| 精品久久久久久久久亚洲| 人体艺术视频欧美日本| 国产精品蜜桃在线观看 |