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

    Nascent Polypeptide-Associated Complex Involved in the Development and Pathogenesis of Fusarium graminearum on Wheat

    2020-09-12 03:22:50XuliWngXinXiebJinLiuGuoLingWngcDewenQiu
    Engineering 2020年5期

    Xuli Wng*, Xin Xieb, Jin Liu Guo-Ling Wngc, Dewen Qiu*

    a State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China

    b College of Agriculture, Guizhou University, Guiyang 550025, China

    c Department of Plant Pathology, The Ohio State University, Columbus, OH 43210, USA

    Keywords:Fusarium head blight Nascent polypeptide-associated complex Gene knockout Pathogenicity Subcellular localization

    A B S T R A C T Reliable knowledge on pathogenic agents contributes to effective plant protection.For most plant pathogens, maintaining protein homeostasis (proteostasis) is essential for unfolding the cellular functions to survive and thrive. However, the fungal proteins involved in proteostasis remain poorly characterized in the process of pathogenesis. In this study, we characterized the function of the nascent polypeptideassociated complex (NAC) in Fusarium graminearum (F. graminearum) (FgNAC), one of the top 10 fungal pathogens with predominant scientific/economic importance. We found that FgNACα, a subunit of FgNAC, manifests high structural and functional similarity to its homologous counterparts in yeast and other species.The mutants of F.graminearum lacking NACα are viable but suffer significant defects in vegetative growth,conidial production,and pathogenesis.In addition,we show here that FgNACα can interact with another subunit of NAC (FgNACβ) in a yeast-two-hybrid assay. The subcellular localization results show that FgNACα and FgNACβ are predominantly localized in the cytoplasm. Future studies should focus on deciphering the mechanism by which NAC orchestrates protein biogenesis and consequentially modulates development and pathogenesis.

    1. Introduction

    Plants in general, and crops in particular, are inevitably and continually exposed to attack from a diverse range of phytopathogenic organisms,each with a unique mode of pathogenesis.Plant diseases have a negative impact on human health directly through agricultural and economic loss,and indirectly through biodiversity collapse and environmental change [1]. In particular, it has been estimated that the losses due to persistent and epidemic outbreaks of fungal and oomycete infection in the five most important crops (i.e., corn, wheat, rice, potatoes, and cassava) would be sufficient to feed 8.5% of the world’s population [2]. A reasonable strategy to cope with long-lasting and recurrent plant diseases should be effective protection and timely control.Characterization of the genes involved in pathogen development and pathogenesis contribute to novel strategies for plant disease control [3].

    The nascent polypeptide-associated complex (NAC) is a functionally versatile protein complex that is involved in protein biogenesis, assembly, and transportation, and that plays a key role in maintaining protein homeostasis (proteostasis) [4]. In all eukaryotes, NAC is composed of a heterodimer consisting of an α- and a β-subunit[5-7].The heterodimeric NAC associates with ribosomes in a 1:1 stoichiometry[8]and dynamically interacts with both ribosomes and nascent polypeptides in a reversible fashion[9],thereby performing a chaperone-like function [10-12]. In yeast Saccharomyces cerevisiae,NAC assists in the folding and maturation of newly synthesized proteins,prevents nascent polypeptides from inappropriate targeting to the endoplasmic reticulum [13,14],initiates protein targeting to mitochondria in vivo [15], and functionally connects to other chaperone networks such as the Hsp70 system[16]. In nematode Caenorhabditis elegans,NAC has a primary role—similar to its role in yeast—as a ribosomeassociated chaperone to regulate translation and to assist in the folding of nascent polypeptides;it is also an essential and beneficial proteostasis sensor that not only detects proteotoxic stress,but also mediates translation upon proteotoxic stress, thus providing the cell with a regulatory feedback mechanism to maintain proteostasis[9]. These unique features of the NACs and their potential roles in cellular operation make them fascinating subjects of studies in a diverse range of organisms from various perspectives.

    In the present study, we characterized the role of NAC in the ascomycete fungus Fusarium graminearum (F. graminearum; teleomorph:Gibberella zeae),with a focus on its pathogenesis on wheat.This fungus is a highly destructive plant pathogen that causes Fusarium head blight(FHB)in wheat and other small-grain cereals[17,18].Infection by F.graminearum not only results in direct yield loss in the field,but also produces mycotoxins within grains during storage, which are detrimental to human and animal health; thus,this fungus poses a threat to food safety that is of global concern[19-22]. Investigation of the genes involved in pathogenesis may provide insights into the mechanism by which F. graminearum interacts with its host plants, thus promoting the development of new strategies for the effective control of FHB. Here, through targeted gene deletion, yeast two-hybrid assays, and subcellular localization combined with pathogenicity testing,we demonstrate that the NACα subunit plays an important role in regulating the development and pathogenesis of F. graminearum in wheat.

    2. Materials and methods

    2.1. Fungal strains and culture conditions

    The wild-type strain PH-1 of F. graminearum was used for the construction of gene deletion mutants. Both the wild-type strain and its derived mutants were routinely cultured on potato dextrose agar (PDA) and V8 agar medium at 25 °C with a 12 h:12 h light:dark cycle for mycelial growth assays.The experiments were performed with three replicate plates. Asexual reproduction was induced in carboxyl methyl cellulose (CMC) [23]. All strains were preserved as conidial suspensions in 20% glycerol at -70 °C.

    2.2. Bioinformatic analysis

    The genome of F. graminearum PH-1 [24] was searched with protein basic local alignment search tool (BLASTp) using the sequence of the NACα subunit(NP_596361,EGD2)of fission yeast,

    Schizosaccharomyces pombe (S. pombe), as the query [25]. The whole-length sequencing of the NACα of the PH-1 was downloaded and then used as the query to search GenBank for homologous sequences from other species (Table S1 in Supplementary data).Homologous amino acid(aa)sequences were aligned using Clustal X version 2.0 [26]. A neighbor-joining phylogenetic tree was constructed using molecular evolutionary genetics analysis version 7.0 (MEGA7),with the node support being assessed through bootstrapping for 1000 replicates [27].

    2.3. Generation of deletion mutants and complementary strains

    The NACα gene deletion constructs were generated from the F. graminearum wild-type strain PH-1 using the split-marker approach[28,29].The upstream 824 base pairs(bp)DNA fragment of FgNACα gene was cloned and inserted into the pUC19-G418Rvector to generate pUC19-A-G418, and the downstream 961 bp DNA fragment was cloned and sequentially inserted into pUC19-AG418 to generate pUC19-A-G418-B(Fig.S1 in Supplementary data).The whole A-G418-B fragment was then amplified and transformed into F.graminearum,as described by Yuan et al.[30].After checking for orientation,the correct gene replacement vector was linearized by digesting with HindIII and was transformed into PH-1 using the polyethylene glycol(PEG)-CaCl2method.The G418-resistant transformants were isolated and purified by means of the single-spore method,and their genotypes were examined by sequencing.

    For complementation assays, a fragment containing the entire gene—including the open reading frame(ORF),the promoter region,and the terminator region—was amplified and cloned to the pMD18-T to generate the complementation construct pMD18-FgNACα. The insert sequence in the complementation construct was confirmed by DNA sequencing. The correct construct was linearized with EcoRI and introduced into the deletion mutants by cotransformation with the HindIII-linearized pUCATPH plasmid to generate the complementary strains. The G418-sensitive and hygromycin-resistant transformants were isolated and their genotypes were confirmed by polymerase chain reaction (PCR) using the primer pair ORFα-F/ORFα-R(Table S2 in Supplementary data).

    2.4. Phenotypic assays

    To measure the growth rate,5 mm mycelium plugs cut from the edge of a 3 d-old colony were inoculated on PDA plates and cultivated at 25°C. The diameter of the colony was measured daily for 4 d. To quantitatively analyze conidiation, four 5 mm mycelium plugs of the aerial hyphae of 2 d-old cultures were introduced into 100 mL of liquid CMC medium and cultured at 25 °C for 5 d.Conidia were harvested by centrifugation at 5000 r·min-1for 5 min and washed twice with sterile distilled water. The concentration of conidia was determined with a hemocytometer in the microscope, as described by Hou et al. [31]. Morphological observation and measurement were conducted under an Olympus BH-2 microscope.The width of the middle intercalary cell and the length from the elongated apical cell to the pedicellate foot cell of 100 conidia were measured for each strain examined.

    2.5. Plant infection assays

    The wheat cv. Minxian169 was used for the infection assay on both spikelets and coleoptiles. Ten microliters of the conidia suspension (4 × 105conidia·mL-1) collected from 5 d-old CMC medium was injected into the third full-sized spikelet from the base of the inflorescence of 35 to 42 d-old plants.The inoculated wheat heads were enclosed in small plastic bags for 48 h to maintain humidity. Symptomatic spikelets in each plant were counted 12 d after inoculation. The disease index scores were calculated as previous reported [32]. The number of inoculated wheat heads per treatment was 9 or 10 for each test,and all tests were repeated at least four times. Wheat coleoptiles were inoculated following the method described by Liu et al.[33].Fifty coleoptiles were inoculated for each strain and then kept at room temperature in a growth chamber. The length of brown lesions was measured 7 d after inoculation.

    2.6. Yeast two-hybrid assays and subcellular localization

    The Matchmaker GAL4 two-hybrid system 3 (Clontech, USA)was used to test interaction between the two subunits of FgNAC.The FgNACα ORF was cloned into pGADT7 to create the vector AD-FgNACα. FgNACβ ORF was amplified from the complementary DNA (cDNA) of PH-1 and cloned into the pGBKT7 vector (BDFgNACβ),which was used as the prey in the yeast two-hybrid assay.After confirmation by means of sequencing,the bait and prey vectors were co-transformed in pairs into the yeast strain AH109 with the Alkali-Cation yeast transformation kit (MP Biomedicals, USA).The Leu+and Trp+transformants were observed for growth on synthetic dextrose medium (SD-Leu-Trp-His, SD-LWH), and the expression of the LacZ reporter gene was measured as described by Zhang et al. [34].

    To test the subcellular localization of FgNACα and FgNACβ in F. graminearum, the coding regions of FgNACα and FgNACβ were amplified separately. The vector pDL2 was used to generate the C-terminal enhanced green fluorescent protein (eGFP) fusion constructs of FgNACα and FgNACβ that are under the control of the native promoter. The fusion constructs were directly used for F.graminearum protoplasts transformation with the general procedures described by Sweigard et al. [35]. All the samples were observed under an Olympus BX61 fluorescence microscope.

    2.7. Statistical analysis

    Statistical analyses were carried out using SPSS 19 software,and the significant difference between treatments was analyzed using one-way analysis of variance (ANOVA). Statistically significant differences p values of less than 0.05 were considered to be statistically significant.

    3. Results

    3.1. FgNACα is orthologous to the yeast NACα subunit

    Only a single NACα subunit gene (FGRAMPH1_01G10263) was found in the F.graminearum PH-1 genome in a BLASTp search with the fission yeast’s NACα subunit EGD2 as the query.We designated this NACα subunit gene of F. graminearum as FgNACα. FgNACα contains 1166 nucleic acids with an ORF of 627 bp,which encodes 209 aas. Similar to EGD2, FgNACα consists of a conserved NAC domain (57 aas) at the N terminal and an ubiquitin-associated(UBA) domain (40 aas) at the C terminal (Figs. 1(a) and (b)). The neighbor-joining tree constructed from the NACα homologous sequences was completely consistent with the phylogeny of the analyzed species[36].The aa sequence of FgNACα has greater than 70% similarity to that of Magnaporthe grisea (XP_366584.1), and together they formed a well-supported (bootstrap support 79%)subclade in the strongly supported (99%) clade including aa sequences of NACα or putative NACα, all from fungi (Fig. 1(c)).These results indicate that the gene FgNACα is the orthologous of the fungal NACα.

    3.2.FgNACα gene disruption and complementation in F. graminearum

    Transformants of the FgNACα deletion mutant (fgnacα) of F. graminearum that could grow in G418 medium were picked and identified by PCR, with the wild-type PH-1 used as the negative control (Fig. S1(b)). No FgNACα band (755 bp) was detected by primer pair ORFα-F/ORFα-R (Fig. S1(b), lane 3) compared with the bands amplified by primer pairs FgNACα-Ft1/Neo-F(2781 bp),indicating that FgNACα was successfully replaced by the neo gene(Fig.S1(b),lanes 1 and 2).To produce the complementary construct of the fgnacα, fgnacα::FgNACα, the 2993 bp of the full-length FgNACα gene, including the promoter and terminator, were cloned and transformed together with HindIII-linearized pUCATPH plasmid into deletion mutant strain fgnacα protoplast(Fig. S1(c)). Transformants that could grow in both G418 and hygromycin medium were selected and identified by PCR (Fig. S1(d)). The FgNACα band (755 bp) was detected by the primer pair ORFα-F/ORFα-R(Fig.S1(d),lane 4),but there was no corresponding band in the deletion mutant strain fgnacα (Fig. S1(d),lane 3),indicating that FgNACα was successfully restored in the deletion mutant strain fgnacα.

    3.3. FgNACα regulates mycelia growth

    In comparison with the wild-type strain PH-1, mycelia growth of the deletion mutant strain fgnacα was significantly reduced, as measured by colony diameter on both potato dextrose agar(PDA)and V8 medium(Fig.2(a)).The deletion mutant strain fgnacα grows slowly during the cultural periods,and the colony diameter was only about 69%of that of the wild-type strain on the third day of inoculation on the PDA plates (Fig. 2(b)). In addition, the deletion mutant strain fgnacα produced fewer aerial hyphae with irregular colony margins (Fig. 2(a)). The phenotype was hardly observable between the wild-type strain PH-1 and the complementary strain fgnacα::FgNACα (Fig. 2). Similarly, in liquid CMC medium, the deletion strain produced much less mycelia biomass in comparison with the wild-type strain and the complementary strain throughout the whole of cultivation, and similar mycelia biomass was harvested from the wild-type and complementary strains (Fig. 2(c)). Taken together, disruption of the FgNACα gene dramatically restricted the vegetation growth of F. graminearum.

    3.4. Disruption of FgNACα causes abnormal conidia production

    Fig. 1. Identification of the FgNACα orthologous gene. (a) Aa sequences comparison of NACα in S. pombe (NP_596361, EGD2) and F. graminearum (FGSG_08560.3, FgNACα);(b) functional domains of NACα in S.pombe and F.graminearum;(c) phylogenetic tree based on aa sequences of NACα. Numbers near the node represent the bootstrapping support value. The FgNACα from F. graminearum is indicated by the red diamond.

    Fig. 2. Deletion of FgNACα causes deficit in mycelia growth. (a) Mycelium growth of the wild-type strain PH-1, the deletion mutant strain fgnacα, and the complementary strain fgnacα::FgNACα on PDA and V8 medium plates at 25°C for 3 d;(b)the colony diameter of each strain was measured after incubation for 1,2,3,and 4 d in PDA medium at 25 °C; (c) the weight of the mycelia was calculated after culturing in liquid CMC medium at 25 °C for 24, 48, and 72 h. The symbol * indicates statistically significant differences p <0.05.

    To examine the role of the FgNACα gene in reproductive development, we compared conidia production and conidial germination of the mutant strains and wild-type strain. We did not find a clear difference in conidial morphology between the deletion mutant strain fgnacα and the wild-type strain or between the deletion mutant strain and the complementary strain.While the number of spores produced by the deletion mutant strain fgnacα was significantly reduced (p <0.05), only about 40% of those produced by the wild-type strain or about 42%of those produced by the complementary strain (Fig. 3(a)). In addition, the conidial germination of the deletion mutant strain fgnacα was significantly delayed,compared with the wild-type and complementary strains;the germination rate of the deletion mutant strain was only 32%of that of the wild-type strain PH-1 after incubating for 4 h in the CMC liquid medium(Fig.3(b)).Interestingly,the germination rate of the deletion mutant strain caught up quickly, accounting for about 86% of that of the wild-type strain after 8 h of cultivation and more than 92%after 12 h of cultivation(Fig.3(b)).The conidia development of the complementary stain was similar to that of the wild-type strain(Fig. 4). These results suggest that the FgNACα gene might be involved in the earlier stage of the conidiation process, and that its function might be compensated by those of other genes.

    3.5. FgNACα is required for the infection of wheat

    The results of the infection assays on wheat cv.Minxian169 are shown in Fig. 4. When the conidia were inoculated on the wheat heads for 12 d, the deletion mutant strain fgnacα only infected about 20%of the inoculated spikelet,whereas the wild-type strain PH-1 and the complementary strain caused over 47% infection rates (Fig. 4(b)). The disease symptom was much more serious when caused by the wild-type strain PH-1 and the complementary strain fgnacα::FgNACα (Fig. 4(a)). Similar patterns were also found when the conidia was inoculated on the coleoptiles(Figs. 4(c) and(d)).The brown lesions caused by the deletion mutant strains were much shorter(~0.54 cm)than those caused by the wild-type strain(1.92 cm) and the complementary strain (1.88 cm). These results suggest that FgNACα manipulates the virulence of F. graminearum on wheat.

    3.6. Interaction between FgNACα and FgNACβ and their subcellular localization

    It has been reported that the α and β subunits of NAC interact with each other and function as a heterodimer in order to facilitate the correct accumulation of the precursors of mitochondrial protein in the cytosol in yeast [37]. To test whether the two subunits functioned in the same way as they do in yeast, we investigated their interaction through yeast two-hybrid assays(Fig.5).We found that FgNACα could interact with FgNACβ in yeast (Fig. 5(a)).This means that FgNACα is physically associated with and functionally related to FgNACβ. They may form a heterodimer and protect the nascent polypeptide from proteolysis.

    Fig. 3. Comparative difference of wild-type strain PH-1, deletion mutant strain fgnacα, and complementary strain fgnacα::FgNACα in spore number and germination rate.(a)The production of conidia in PH-1,fgnacα,and fgnacα::FgNACα strains.(b)Germination was measured by the percentage of germinated conidia in CMC liquid culture after incubating for 4, 8, 12, and 24 h. The symbol * indicates statistically significant differences p <0.05. NS: no significant difference.

    Fig. 4. FgNACα is required for F.graminearum pathogenicity. (a,b) Infection assays with flowering wheat heads. Wheat heads were inoculated with conidia of the wild-type strain PH-1, the deletion mutant strain fgnacα, and the complementary strain fgnacα::FgNACα. The pictures were taken 12 d post inoculation. (c,d) Deletion of FgNACα attenuates the virulence of F. graminearum on wheat coleoptiles. The representative brown lesions on the coleoptiles infected with PH-1, fgnacα, and fgnacα::FgNACα are indicated with arrows. The pictures were taken 7 d post inoculation. The symbol * indicates statistically significant differences p <0.05.

    Fig. 5. FgNACα interacts with FgNACβ. (a) Yeast two-hybrid assays of FgNACα and FgNACβ. FgNACα was fused with the activation domain (AD) of GAL4; FgNACβ was fused with the binding domain(BD)of GAL4. The transformants of AH109 yeast were diluted 10,100,and 1000 times,and grown on synthetic dextrose medium without Leu and Trp (SD-LW), and on the plate SD-LWH. (b) Subcellular localization of FgNACα and FgNACβ in F. graminearum. Conidia expressing the FgNACα-GFP or FgNACβ-GFP fusion construct were examined. All the samples were observed under fluorescence microscope.

    In order to investigate the subcellular localization of FgNACα and FgNACβ, we generated the C-terminal GFP fusion constructs of the FgNACα and FgNACβ coding regions with the native promoters and transformed them to the wild-type strain PH-1 protoplasts. Under the fluorescence microscope, the green signals were detected in the cytoplasm of the conidia in transformation strains with FgNACα-GFP and FgNACβ-GFP (Fig. 5(b)). The results suggest that FgNACα and FgNACβ can have the same localization in the F. graminearum cells, which may provide chances for their physical interaction.

    4. Discussion

    NAC is the first discovered nonribosomal factor that cotranslationally interacts with newly synthesized polypeptides [14].Together with signal recognition particles,NAC promotes the fidelity of protein targeting to the endoplasmic reticulum [38]. Since its discovery more than 20 years ago, various potential roles for NAC have been demonstrated in several evolutionarily dispersed model organisms, such as yeast [15,39,40], tobacco [41], and nematode [9]. From these studies, diverse roles in maintaining proteostasis have been discovered[4,42].In most cases,such functional multiplication was often companied by evolutionary diversification. However, we found that NACα is evolutionarily highly conserved in terms of both the primary sequence structure and the functional domains (Fig. 1). From animals to plants and to fungi, the aa sequences of NACα are very similar and harbor an NAC domain and a UBA domain(Figs.1(a)and(b)).The phylogeny constructed from NACα sequences is highly congruent with the species evolutionary history (Fig. 1(c)). Given the fact that the αand β-subunits of NAC form a heterodimer, which was found here(Fig. 5) and in other species [5,6], it is reasonable to predict that NACβ is also evolutionarily conserved. Nevertheless, this evolution-function paradox can be explained by the fact that NAC is a key factor in protein biogenesis and functionally connects to other metabolic networks, such as the Hsp70 system [16]. Such functional versatility and evolutionary conservatism imply that the results found in F.graminearum in the present study are likely to be transferable to other plant pathogens and vice versa, and that the disease management methods developed by targeting NAC will have trans-disease applicability.

    Although the biological function of NACα has been extensively characterized in model organisms,its role in non-model organisms is still poorly understood. In the present study, we found that the NACα subunit gene of F. graminearum is not only involved in vegetative growth, but also regulates its pathogenicity positively.In other plant pathogens, it has been reported that NACα in Sclerotinia sclerotiorum regulates its pathogenicity negatively [43].Such a difference may be explained by the different lifestyle of the two pathogens. F. graminearum is a hemi-biotrophic fungus, while S. sclerotiorum is a necrotrophic fungus. In animals, NACα is also involved in the innate immune response to pathogens; its expression was significantly up-regulated when a crab was exposed to Vibrio anguillarum challenges [44]. Similarly, in Japanese flounder,NACα gene expression was up-regulated in fish tissues in response to Edwardsiella tarda challenges, and overexpression of NACα enhanced resistant gene expression [45]. In plants, it has been demonstrated that silencing the NACα gene in tobacco (Nicotiana benthamiana) inhibits the cell-to-cell movement of the brome mosaic virus[41].Moreover,NACα protein shows higher abundance in tomato resistant lines compared with susceptible lines after inoculation with Ralstonia solanacearum [46]. Clearly, the function of NACα in wheat or host plants in general and the molecular interaction between the pathogen-host systems need to be fully characterized before we can precisely determine the role of NACα in pathogenesis.

    For fungal pathogens,the ability to invade the plant cell is a prerequisite for successful colonization and infection,while the ability to maintain proteostasis is essential for unfolding the cellular functions to survive and thrive. The fungal molecules involved in proteostasis remain poorly characterized in phytopathogens. In this study,we characterized FgNACα,an evolutionarily highly conserved regulator in maintaining proteostasis,which participates in fungal development and pathogenicity. In addition, our results suggest that the function of FgNACα is likely unfolded by forming a heterodimer with FgNACβ (Fig. 5). However, pathogenesis is unfolded through molecular interactions not only within pathogens and hosts,but also between factors from pathogens and from hosts.It will be very interesting to figure out the molecular mechanism of how NAC performs different functions in pathogen-host interaction systems.

    5. Conclusion

    Effective plant protection technologies often stem from a reliable understanding of pathogenic agents. Such basic knowledge will undoubtedly help to produce novel insight into pathogenesis,and will lead to innovative plant protection strategies. In the present study,we found that the NACα of the plant pathogenic fungus F. graminearum manifests high structural and functional similarity to its homologous counterparts from yeast and other fungal species. The mutants of F. graminearum that lack NACα are viable,but suffer significant defects in vegetative growth,conidial production,and pathogenesis on wheat.These functional consequences of NACα deletion might stem from the important role that NAC plays in proteostasis,as found in yeast and other model organisms.From the perspective of plant protection, NAC deserves recognition in the development of innovative control methods for the destructive plant disease, FHB.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China (31471737, 31671984, and 31801691) and the Advanced Programs of Guizhou Province for the Returned Overseas Scholars ([2018]02).

    Compliance with ethics guidelines

    Xuli Wang, Xin Xie, Jin Liu, Guo-Liang Wang, and Dewen Qiu declare that they have no conflict of interest or financial conflicts to disclose.

    Appendix A. Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.eng.2019.07.025.

    在线观看美女被高潮喷水网站| 2021天堂中文幕一二区在线观| 精品午夜福利视频在线观看一区| 久久久久九九精品影院| avwww免费| 国产精品久久久久久久电影| 高清在线国产一区| 午夜激情欧美在线| 亚洲精品影视一区二区三区av| 精品久久国产蜜桃| 亚洲精品成人久久久久久| 欧美区成人在线视频| 真实男女啪啪啪动态图| 久久草成人影院| 午夜福利视频1000在线观看| 18禁在线播放成人免费| 亚洲精品成人久久久久久| 婷婷精品国产亚洲av| 国产精品伦人一区二区| 亚洲国产精品成人综合色| 欧美zozozo另类| 露出奶头的视频| 国产av不卡久久| 国产精品人妻久久久影院| 亚州av有码| 啪啪无遮挡十八禁网站| 亚洲精华国产精华精| 高清毛片免费观看视频网站| 久久九九热精品免费| 日本撒尿小便嘘嘘汇集6| 又爽又黄无遮挡网站| 亚洲四区av| 三级毛片av免费| 三级毛片av免费| 日本撒尿小便嘘嘘汇集6| 三级毛片av免费| 亚洲成人久久爱视频| www.www免费av| 免费人成视频x8x8入口观看| 亚洲精品国产成人久久av| 亚洲国产欧美人成| 精华霜和精华液先用哪个| 午夜视频国产福利| 十八禁国产超污无遮挡网站| 欧美国产日韩亚洲一区| 女人被狂操c到高潮| bbb黄色大片| 99久久无色码亚洲精品果冻| 久久精品国产自在天天线| 91久久精品电影网| av在线蜜桃| 大型黄色视频在线免费观看| 麻豆久久精品国产亚洲av| 老熟妇乱子伦视频在线观看| 欧美高清成人免费视频www| 亚洲欧美日韩高清在线视频| 亚洲美女搞黄在线观看 | 久久99热这里只有精品18| av在线蜜桃| 久久精品国产亚洲av涩爱 | 日本黄大片高清| 又紧又爽又黄一区二区| 很黄的视频免费| 男人和女人高潮做爰伦理| 少妇裸体淫交视频免费看高清| 亚洲av日韩精品久久久久久密| 村上凉子中文字幕在线| 91午夜精品亚洲一区二区三区 | 琪琪午夜伦伦电影理论片6080| 国产精品免费一区二区三区在线| 综合色av麻豆| 在线观看午夜福利视频| 久久九九热精品免费| 亚洲人与动物交配视频| 成人性生交大片免费视频hd| 一区二区三区四区激情视频 | 两性午夜刺激爽爽歪歪视频在线观看| 人人妻人人看人人澡| av在线亚洲专区| 搞女人的毛片| 欧美人与善性xxx| 中国美白少妇内射xxxbb| 尾随美女入室| 久久久久久久久久成人| 大又大粗又爽又黄少妇毛片口| 中文字幕精品亚洲无线码一区| 久久久色成人| 久久天躁狠狠躁夜夜2o2o| 亚洲性夜色夜夜综合| 美女xxoo啪啪120秒动态图| 哪里可以看免费的av片| 91麻豆精品激情在线观看国产| 亚洲美女视频黄频| 国产精品国产三级国产av玫瑰| 香蕉av资源在线| 精华霜和精华液先用哪个| 成人国产一区最新在线观看| 成人国产一区最新在线观看| 国产精品无大码| 人人妻人人澡欧美一区二区| 亚洲精品成人久久久久久| 国产精品福利在线免费观看| 亚洲av中文av极速乱 | 国产成年人精品一区二区| 少妇丰满av| 深夜a级毛片| 久久热精品热| 精品不卡国产一区二区三区| 中国美白少妇内射xxxbb| 99精品在免费线老司机午夜| 亚洲男人的天堂狠狠| 又爽又黄a免费视频| 啦啦啦观看免费观看视频高清| 内射极品少妇av片p| 欧美激情在线99| 淫秽高清视频在线观看| 床上黄色一级片| 我要看日韩黄色一级片| 成人午夜高清在线视频| 91久久精品国产一区二区成人| 国产精品久久久久久av不卡| 最后的刺客免费高清国语| 日韩欧美精品免费久久| 在线观看午夜福利视频| 亚洲欧美日韩高清专用| 色精品久久人妻99蜜桃| 亚洲av熟女| 欧美成人a在线观看| 欧美黑人巨大hd| 亚洲欧美日韩卡通动漫| 亚洲最大成人手机在线| 久久久色成人| 亚洲精品456在线播放app | 欧美人与善性xxx| 亚洲成av人片在线播放无| 国产色婷婷99| 亚洲在线观看片| 国产大屁股一区二区在线视频| 一本精品99久久精品77| 国产在视频线在精品| 尤物成人国产欧美一区二区三区| 直男gayav资源| 亚洲人成网站在线播| 老司机福利观看| 老司机福利观看| 如何舔出高潮| 中出人妻视频一区二区| 亚洲精品456在线播放app | 国产女主播在线喷水免费视频网站 | 免费看光身美女| 韩国av一区二区三区四区| 亚洲国产日韩欧美精品在线观看| 大又大粗又爽又黄少妇毛片口| 老司机午夜福利在线观看视频| 成人二区视频| 特大巨黑吊av在线直播| 夜夜夜夜夜久久久久| 国产麻豆成人av免费视频| 国产日本99.免费观看| 成人无遮挡网站| 成人无遮挡网站| 99国产精品一区二区蜜桃av| 淫秽高清视频在线观看| 色综合站精品国产| 亚洲av日韩精品久久久久久密| 久久精品综合一区二区三区| 床上黄色一级片| 精品久久国产蜜桃| 国产高清视频在线观看网站| 亚洲精品粉嫩美女一区| 女生性感内裤真人,穿戴方法视频| 国产精品久久久久久精品电影| 床上黄色一级片| 观看美女的网站| 欧美xxxx黑人xx丫x性爽| 99久久无色码亚洲精品果冻| 一个人观看的视频www高清免费观看| 国产精品久久视频播放| 99riav亚洲国产免费| 亚洲国产精品成人综合色| 亚洲欧美日韩卡通动漫| 亚洲第一电影网av| 成人性生交大片免费视频hd| 51国产日韩欧美| 特大巨黑吊av在线直播| av在线天堂中文字幕| 我要搜黄色片| 日韩亚洲欧美综合| 麻豆成人午夜福利视频| 欧美xxxx性猛交bbbb| 黄色丝袜av网址大全| 淫妇啪啪啪对白视频| 亚洲无线观看免费| 全区人妻精品视频| 嫩草影视91久久| 久久亚洲精品不卡| 久久久精品欧美日韩精品| 色播亚洲综合网| 成年人黄色毛片网站| 成年免费大片在线观看| 黄色丝袜av网址大全| 国产在线男女| 搡老熟女国产l中国老女人| 一个人看视频在线观看www免费| 人妻夜夜爽99麻豆av| avwww免费| 国产国拍精品亚洲av在线观看| 国产亚洲91精品色在线| 乱人视频在线观看| 色精品久久人妻99蜜桃| 午夜精品在线福利| 一卡2卡三卡四卡精品乱码亚洲| 尤物成人国产欧美一区二区三区| 亚洲无线观看免费| 91午夜精品亚洲一区二区三区 | 免费不卡的大黄色大毛片视频在线观看 | 久久国产精品人妻蜜桃| 99久久精品一区二区三区| 国产亚洲精品av在线| 亚洲中文字幕一区二区三区有码在线看| 婷婷精品国产亚洲av在线| 精品久久久久久,| 哪里可以看免费的av片| 国模一区二区三区四区视频| 少妇人妻精品综合一区二区 | 看免费成人av毛片| 国产午夜福利久久久久久| 日日啪夜夜撸| 久久精品夜夜夜夜夜久久蜜豆| 欧美性猛交黑人性爽| 欧美色视频一区免费| 亚洲图色成人| 黄色一级大片看看| 此物有八面人人有两片| 亚洲av免费在线观看| 午夜亚洲福利在线播放| 国内毛片毛片毛片毛片毛片| 久久午夜福利片| 日本撒尿小便嘘嘘汇集6| 精品国内亚洲2022精品成人| 久久久久国产精品人妻aⅴ院| 免费av不卡在线播放| 日韩欧美精品v在线| 亚洲国产欧洲综合997久久,| 性色avwww在线观看| 亚洲久久久久久中文字幕| 伦理电影大哥的女人| www日本黄色视频网| 一个人免费在线观看电影| 精品一区二区三区视频在线| 精品99又大又爽又粗少妇毛片 | 国产男人的电影天堂91| 一区二区三区免费毛片| 啦啦啦啦在线视频资源| 免费av观看视频| 婷婷亚洲欧美| 欧美一区二区精品小视频在线| 久久人人爽人人爽人人片va| 亚洲国产精品成人综合色| 最近视频中文字幕2019在线8| 12—13女人毛片做爰片一| 亚洲 国产 在线| 日本精品一区二区三区蜜桃| 超碰av人人做人人爽久久| 1024手机看黄色片| 蜜桃亚洲精品一区二区三区| 欧美日韩瑟瑟在线播放| 色噜噜av男人的天堂激情| 婷婷色综合大香蕉| 免费大片18禁| 很黄的视频免费| 久久久久国产精品人妻aⅴ院| 精品一区二区免费观看| 男人的好看免费观看在线视频| 中文在线观看免费www的网站| h日本视频在线播放| 久久久久久久久大av| 午夜精品在线福利| 联通29元200g的流量卡| 两性午夜刺激爽爽歪歪视频在线观看| 色5月婷婷丁香| 欧美高清成人免费视频www| 日韩大尺度精品在线看网址| 看黄色毛片网站| 亚洲色图av天堂| 午夜福利欧美成人| 国产精品精品国产色婷婷| 男插女下体视频免费在线播放| 校园春色视频在线观看| 91狼人影院| 日韩精品青青久久久久久| 亚洲av日韩精品久久久久久密| 久久午夜亚洲精品久久| 精品一区二区三区av网在线观看| 国产精品一区二区性色av| 最近视频中文字幕2019在线8| 蜜桃亚洲精品一区二区三区| 成人国产综合亚洲| 国产精品,欧美在线| 能在线免费观看的黄片| 两性午夜刺激爽爽歪歪视频在线观看| 99久久精品热视频| 一边摸一边抽搐一进一小说| 级片在线观看| 久久精品91蜜桃| 亚洲av美国av| 国产成人aa在线观看| 亚洲av电影不卡..在线观看| 亚洲性久久影院| 久久欧美精品欧美久久欧美| 最近视频中文字幕2019在线8| 麻豆av噜噜一区二区三区| 亚洲无线观看免费| 波多野结衣高清无吗| 国产一区二区激情短视频| 国产精品伦人一区二区| 国产乱人伦免费视频| 免费不卡的大黄色大毛片视频在线观看 | 日韩欧美免费精品| 很黄的视频免费| 国产伦在线观看视频一区| 精品人妻视频免费看| 日本爱情动作片www.在线观看 | 一本一本综合久久| 别揉我奶头 嗯啊视频| 日本熟妇午夜| 亚洲乱码一区二区免费版| 九色成人免费人妻av| 高清毛片免费观看视频网站| 九九爱精品视频在线观看| 日本熟妇午夜| 国产精品无大码| 嫩草影院精品99| 国产成人aa在线观看| 国产欧美日韩精品一区二区| 国产午夜福利久久久久久| 他把我摸到了高潮在线观看| 一个人观看的视频www高清免费观看| 亚洲天堂国产精品一区在线| 欧美又色又爽又黄视频| 啪啪无遮挡十八禁网站| 国产精品久久电影中文字幕| 国产伦在线观看视频一区| 国产精品久久久久久亚洲av鲁大| 国产国拍精品亚洲av在线观看| 男插女下体视频免费在线播放| 欧洲精品卡2卡3卡4卡5卡区| 中国美白少妇内射xxxbb| 此物有八面人人有两片| 91久久精品国产一区二区三区| 亚洲七黄色美女视频| 春色校园在线视频观看| 日本三级黄在线观看| 久久久久久久久久久丰满 | 悠悠久久av| 欧美日韩乱码在线| 久久精品国产99精品国产亚洲性色| 亚洲性久久影院| 婷婷精品国产亚洲av在线| 最近视频中文字幕2019在线8| 少妇人妻精品综合一区二区 | 麻豆成人av在线观看| 久久精品国产亚洲av香蕉五月| а√天堂www在线а√下载| 免费看av在线观看网站| 中文字幕av在线有码专区| 三级国产精品欧美在线观看| 色吧在线观看| 国内精品宾馆在线| 午夜精品一区二区三区免费看| 免费观看在线日韩| 亚洲18禁久久av| 精品欧美国产一区二区三| 一本精品99久久精品77| 熟女电影av网| 国产精品,欧美在线| 国产亚洲av嫩草精品影院| 精品一区二区三区av网在线观看| 成年女人看的毛片在线观看| АⅤ资源中文在线天堂| 国产精品嫩草影院av在线观看 | 成人特级黄色片久久久久久久| 婷婷精品国产亚洲av在线| 乱人视频在线观看| 国产成人av教育| 中文字幕av在线有码专区| 自拍偷自拍亚洲精品老妇| 我要搜黄色片| 欧美成人免费av一区二区三区| 看免费成人av毛片| 久久精品夜夜夜夜夜久久蜜豆| 色精品久久人妻99蜜桃| 国产aⅴ精品一区二区三区波| 永久网站在线| 国产91精品成人一区二区三区| 免费观看的影片在线观看| 成人综合一区亚洲| 大型黄色视频在线免费观看| 狂野欧美激情性xxxx在线观看| 91麻豆av在线| 亚洲精品一卡2卡三卡4卡5卡| 久久久成人免费电影| 国产精品无大码| 一个人免费在线观看电影| 亚洲久久久久久中文字幕| 成年女人看的毛片在线观看| 亚洲精品粉嫩美女一区| 久久久久性生活片| 亚洲午夜理论影院| 美女免费视频网站| 岛国在线免费视频观看| 欧美高清性xxxxhd video| 又紧又爽又黄一区二区| 亚洲熟妇中文字幕五十中出| 午夜激情福利司机影院| 一个人看视频在线观看www免费| 男人舔女人下体高潮全视频| 精品久久久久久久末码| 老师上课跳d突然被开到最大视频| 人妻夜夜爽99麻豆av| 国产精品自产拍在线观看55亚洲| 免费看日本二区| 国产白丝娇喘喷水9色精品| 国产高清激情床上av| 最后的刺客免费高清国语| 精品国内亚洲2022精品成人| 日日啪夜夜撸| 99国产精品一区二区蜜桃av| av黄色大香蕉| 久久久久久九九精品二区国产| 美女xxoo啪啪120秒动态图| 欧美极品一区二区三区四区| 国产毛片a区久久久久| 美女高潮喷水抽搐中文字幕| 亚洲人与动物交配视频| 国产精品乱码一区二三区的特点| 丰满的人妻完整版| 精品人妻熟女av久视频| 老司机深夜福利视频在线观看| 男女做爰动态图高潮gif福利片| 一本一本综合久久| 欧美日本视频| 久久精品国产亚洲av香蕉五月| 好男人在线观看高清免费视频| 看黄色毛片网站| 一a级毛片在线观看| 99久久精品国产国产毛片| 午夜日韩欧美国产| 午夜a级毛片| 久9热在线精品视频| 内射极品少妇av片p| 国产91精品成人一区二区三区| 一进一出抽搐动态| 亚洲国产精品合色在线| 久久人人爽人人爽人人片va| 国产主播在线观看一区二区| 免费观看人在逋| 日本一本二区三区精品| 2021天堂中文幕一二区在线观| 日韩高清综合在线| 在线国产一区二区在线| 黄色女人牲交| 亚洲人成伊人成综合网2020| 国产老妇女一区| 国国产精品蜜臀av免费| 亚洲男人的天堂狠狠| 成人毛片a级毛片在线播放| 亚洲av第一区精品v没综合| 黄色丝袜av网址大全| 日韩在线高清观看一区二区三区 | 国产精品精品国产色婷婷| 又爽又黄a免费视频| 色在线成人网| 18禁裸乳无遮挡免费网站照片| 国产精品人妻久久久久久| 免费电影在线观看免费观看| 美女高潮的动态| 亚洲av五月六月丁香网| 欧美性猛交黑人性爽| 国产av一区在线观看免费| 亚洲真实伦在线观看| 18禁裸乳无遮挡免费网站照片| 国产成人一区二区在线| 丰满乱子伦码专区| av在线老鸭窝| 亚洲自拍偷在线| 成年人黄色毛片网站| 亚洲av中文字字幕乱码综合| 他把我摸到了高潮在线观看| 一级av片app| 国产一区二区三区视频了| 一区二区三区四区激情视频 | 国产精品国产三级国产av玫瑰| 美女高潮的动态| 日韩中文字幕欧美一区二区| 黄色视频,在线免费观看| av在线老鸭窝| eeuss影院久久| 欧美xxxx性猛交bbbb| 国产精品98久久久久久宅男小说| 亚洲狠狠婷婷综合久久图片| 国产av一区在线观看免费| 国产男靠女视频免费网站| 欧美在线一区亚洲| 小蜜桃在线观看免费完整版高清| 香蕉av资源在线| 亚洲精品456在线播放app | 性插视频无遮挡在线免费观看| 黄色欧美视频在线观看| 97人妻精品一区二区三区麻豆| 久久天躁狠狠躁夜夜2o2o| 午夜免费男女啪啪视频观看 | 最好的美女福利视频网| 久久午夜亚洲精品久久| 国内精品一区二区在线观看| 色哟哟·www| 久久精品国产鲁丝片午夜精品 | av女优亚洲男人天堂| 三级男女做爰猛烈吃奶摸视频| 露出奶头的视频| 啦啦啦观看免费观看视频高清| a级毛片免费高清观看在线播放| 日韩欧美在线二视频| 亚洲自拍偷在线| 俺也久久电影网| 亚洲经典国产精华液单| 欧美日韩乱码在线| 色综合婷婷激情| 天堂动漫精品| 精品人妻一区二区三区麻豆 | 国产激情偷乱视频一区二区| 婷婷六月久久综合丁香| 性欧美人与动物交配| 动漫黄色视频在线观看| 欧美黑人巨大hd| 国产伦在线观看视频一区| 亚洲内射少妇av| 国产大屁股一区二区在线视频| 亚洲国产精品sss在线观看| 国产成人aa在线观看| av.在线天堂| 国产精品久久视频播放| 看免费成人av毛片| 亚洲狠狠婷婷综合久久图片| 老师上课跳d突然被开到最大视频| 欧美日本亚洲视频在线播放| 国产成人aa在线观看| 国产精品亚洲一级av第二区| 欧美激情在线99| 尾随美女入室| 国产精品爽爽va在线观看网站| 啦啦啦韩国在线观看视频| 久久精品国产清高在天天线| 欧美成人免费av一区二区三区| 免费看a级黄色片| 很黄的视频免费| 欧美一区二区精品小视频在线| 亚洲va在线va天堂va国产| 成年版毛片免费区| 亚洲七黄色美女视频| 熟妇人妻久久中文字幕3abv| 亚洲性夜色夜夜综合| 久久久久免费精品人妻一区二区| 成年女人永久免费观看视频| 亚洲av成人精品一区久久| 久久精品国产亚洲av香蕉五月| 欧美另类亚洲清纯唯美| 国产一区二区三区av在线 | 神马国产精品三级电影在线观看| 国产人妻一区二区三区在| 99久久久亚洲精品蜜臀av| 亚洲欧美日韩高清专用| 精品一区二区三区视频在线观看免费| 国模一区二区三区四区视频| 亚洲色图av天堂| 欧美黑人巨大hd| 国产白丝娇喘喷水9色精品| 欧美在线一区亚洲| 黄色配什么色好看| 91在线精品国自产拍蜜月| 免费看光身美女| 亚洲av第一区精品v没综合| 人妻制服诱惑在线中文字幕| 岛国在线免费视频观看| 看免费成人av毛片| 女人十人毛片免费观看3o分钟| 看免费成人av毛片| 女人十人毛片免费观看3o分钟| 又黄又爽又免费观看的视频| 两个人视频免费观看高清| 人妻制服诱惑在线中文字幕| 岛国在线免费视频观看| 国产精品国产高清国产av| 欧美一区二区亚洲| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av第一区精品v没综合| 国产精品乱码一区二三区的特点| 亚洲七黄色美女视频| 黄色日韩在线| 日本一二三区视频观看| 国产精品一区二区性色av| 精品人妻熟女av久视频| 日韩精品青青久久久久久| 国产日本99.免费观看| 欧美日本亚洲视频在线播放| 欧美3d第一页| 男人舔女人下体高潮全视频| 亚洲不卡免费看| 国产精品一及| 99热精品在线国产| 久久精品综合一区二区三区| 国产精品人妻久久久久久| 18禁在线播放成人免费| 伊人久久精品亚洲午夜| 99国产精品一区二区蜜桃av|