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

    Rapid Detection of Wheat Blast Pathogen Magnaporthe oryzae Triticum Pathotype Using Genome-Specific Primers and Cas12a-mediated Technology

    2021-03-22 07:37:50HouxingKngPngKngyuHuYufiDngMriBllizziDipliRniGuptNurUinMhmuAlfroUrshimSnjoyKumrPulGryPtrsonYilinZhouXupingZhouTofzzlIslmGuoLingWng
    Engineering 2021年9期

    Houxing Kng*, Y Png#, Kngyu Hu, Yufi Dng Mri Bllizzi, Dipli Rni Gupt,Nur Uin Mhmu, Alfro S. Urshim, Snjoy Kumr Pul, Gry Ptrson, Yilin ZhouXuping Zhou M Tofzzl Islm,*, Guo-Ling Wng,*

    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 Department of Plant Pathology, The Ohio State University, Columbus, OH 43210, USA

    c Institute of Biotechnology and Genetic Engineering (IBGE), Bangabandhu Sheikh Mujibur Rahman Agricultural University, Dhaka 1706, Bangladesh

    d Centro de Ciências Agrárias, Universidade Federal de S?o Carlos, Araras, SP CEP 13600-000, Brazil

    e Foreign Disease-Weed Science Research Unit (FDWSRU), United States Department of Agriculture-Agricultural Research Service (USDA-ARS), Ft. Detrick, MD 21702, USA

    Keywords:Wheat blast Magnaporthe oryzae Triticum Cas12a Nucleic acid rapid lateral flow immunoassay Field detection

    ABSTRACT Wheat blast,caused by the fungus Magnaporthe oryzae Triticum(MoT)pathotype,is a devastating disease persistent in South America and Bangladesh.Since MoT generally fails to cause visual symptoms in wheat until the heading stage when the infection would have advanced,disease control by fungicide application solely based on the detection of visual symptoms is ineffective. To develop an accurate and sensitive method to detect MoT at the seedling and vegetative stages for disease control, we sequenced the genomes of two MoT isolates from Brazil and identified two DNA fragments, MoT-6098 and MoT-6099, that are present in the MoT genome but not in the genome of the rice-infecting Magnaporthe oryzae Oryzae(MoO)pathotype. Using polymerase chain reaction(PCR),we confirmed the specificity of the two markers in 53 MoT and MoO isolates from South America and Bangladesh. To test the efficiency of the two markers, we first established a loop-mediated isothermal amplification (LAMP) method to detect MoT at isothermal conditions, without the use of a PCR machine. Following this, we used the Cas12a protein and guide RNAs(gRNAs)to target the MoT-6098 and MoT-6099 sequences.The activated Cas12a showed indiscriminate single-stranded deoxyribonuclease (ssDNase) activity. We then combined targetdependent Cas12a ssDNase activation with recombinase polymerase amplification(RPA)and nucleic acid lateral flow immunoassay(NALFIA)to develop a method that accurately,sensitively,and cost-effectively detects MoT-specific DNA sequences in infected wheat plants.This novel technique can be easily adapted for the rapid detection of wheat blast and other important plant diseases in the field.

    1. Introduction

    The hemibiotrophic fungal pathogen, Magnaporthe oryzae(M.oryzae),affects many cereals including the two most important staple crops, rice and wheat. Rice blast, caused by M. oryzae Oryza(MoO)pathotype,was reported nearly three centuries ago in China and Japan and is now found in over 85 countries [1]. Wheat blast,caused by M.oryzae Triticum(MoT)pathotype,was initially identified in six municipalities of the Parana State in Brazil in 1985 [2],and its distribution gradually expanded in South America to about 3 × 106hm2, in the early 1990s [3]. For the first time, the disease was spotted outside of South America in eight districts in Bangladesh in 2016, where it rapidly became a serious threat to wheat production [4]. In early 2017, various newspapers reported that symptoms resembling those of wheat blast were observed for the first time in Indian wheat fields close to the India–Bangladesh border[5].Between 2016 and 2019,the disease spread to a total of 20 districts in Bangladesh,and it continues to pose a serious threat to the food nutritional security of Bangladesh [6]. Based on the agro-climatic conditions of the region where wheat blast was first identified in Bangladesh in 2016,a recent study reported that up to 7 × 106hm2of wheat fields in India,Pakistan,and Bangladesh are vulnerable to MoT infection [7]. With a conservative estimate of 5%–10% yield loss, the annual potential economic loss could be 132 million–264 million USD. Therefore, an integrated approach is urgently required to control the disease in Bangladesh and to prevent this disease from spreading to other regions of the world.

    It remains unclear how MoT infects wheat plants,the identity of its alternative hosts, and those climate factors that favor an epidemic.Although MoT is capable of infecting both leaves and spikes of the wheat plant,spike infections are of greater economic significance. Infection in the rachis or peduncle can block the translocation of photosynthates, killing the upper parts of the spike. The disease can cause yield losses of up to 100%under pathogen favorable conditions [8]. However, as spike blast is seen developing in fields without visual symptoms of leaf blast, the correlation between the two is unclear. Therefore, controlling spike blast through application of fungicides can be difficult as farmers cannot always use lesions on leaves before the heading stage to warrant fungicide application. Additionally, wheat blast is also a seedtransmittable disease [9]. Screening seeds and certifying them for MoT will reduce the risk of further spread of the disease to other countries. Therefore, development of a rapid, sensitive, and lowcost method for MoT detection in wheat seedlings in the field, in alternative hosts in nearby areas, and in wheat seeds in commercial shipments, is critical for the control of wheat blast disease in South America and Bangladesh.

    Pieck et al. [10] identified DNA markers that appeared to be unique to MoT strains. One of the markers, MoT3, showed specificity when tested on DNA from 284 M.oryzae strains from 11 host species,collected from several countries.A 360-base pair(bp)fragment of the MoT3 marker was amplified in the MoT strains but not in the MoO strains. However, the fragment is located within a single-copy gene encoding a retinol dehydrogenase, which is present in MoT, MoO, and M. oryzae Lolium (MoL) strains. A recent study analyzed the polymorphism among 81 previously assembled Magnaporthe genomes and found that the MoT3 sequence is not present in all the MoT strains [11]. A new DNA marker, C17, was identified that is specific to the Triticum lineage with high sensitivity [11]. In quantitation real-time (qRT)-polymerase chain reaction (PCR) reactions, it successfully amplified DNA from 100% of the 30 wheat-borne strains used in the study. However, it is not clear whether both MoT3 and C17 are useful for rapid MoT detection in infected wheat plants or not in the field.

    In recent years, several new technologies have been developed to rapidly detect plant diseases in both laboratory and field conditions. For example, the loop-mediated isothermal amplification(LAMP) [12,13] and recombinase polymerase amplification (RPA)[14] are being widely used for plant disease detection. It is worth mentioning that a revolutionary detection method for human diseases, which integrates the programmable nucleotide-target carried by Cas proteins, such as Cas12a and Cas13a, with the nucleic acid rapid lateral flow immunoassay (NALFIA) method, was recently developed [15–17]. This method is an accurate, rapid,specific, and cost-effective alternative that can be easily adapted for plant disease detection.

    In this study,we aimed to develop a novel method for the rapid and accurate detection of MoT on field,using RPA,Cas12a,and NALFIA. First, we de novo sequenced the genome of two wheat blast strains from Brazil using next generation sequencing(NGS)strategy.Second,we performed a pair-wise gene sequence alignment analysis between the MoO and MoT genomes using the BLAT software[18],which led to the identification of two MoT-specific fragments:MoT-6098 and MoT-6099. Third, we designed a few sets of primers based on the MoT-6098 and MoT-6099 sequences and verified the sensitivity of the primers via PCR and LAMP. Finally, we combined the Cas12a-based detection method with RPA and NALFIA technologies and developed a sensitive and low-cost method for detecting MoT in infected wheat plants in the field. Our study provides a new detection method for the devastating wheat blast disease.

    2. Materials and methods

    2.1. Fungal materials

    The MoO strains used in this study (RB22, R01-1, P131, Guy11,N60, N63, and N71) were cultured in the Institute of Plant Protection, Chinese Academy of Agricultural Science, China. The DNA samples of all of the MoT strains were prepared at the Bangabandhu Sheikh Mujibur Rahman Agricultural University; the Foreign Disease-Weed Science Research Unit of United States Department of Agriculture-Agricultural Research Service at Ft.Detrick, USA; and the Centro de Ciências Agrárias, Universidade Federal de S?o Carlos, Brazil.

    2.2. Primers used in this study

    All primers used in this study are listed in Table 1.

    2.3. Genome sequencing and comparative genomic analysis

    The genomic DNA of MoT strains PR01-37.V.1 and PR01-37.V.3 was sequenced using the Illumina sequencing platform (Hisep X_ten, BioMarker Co., Ltd., China). In brief, DNA of PR01-37.V.1 and PR01-37.V.3 were randomly sheared into ~500 bp fragments and cloned to the sequencing vector for library construction,before the library was sequenced using the Illumina system. The clean sequence read data were obtained after sequence quality control.A de novo assembly of the two genomes was performed using SPAdes-3.13.0 software (St. Petersburg State University, Russia).The Augustus.2.5.5 software (University of Gottingen, Germany)was used for gene prediction of MoO and MoT genomes. A pairwise gene alignment pipeline,which utilized BLAT[18]as the core alignment tool, was then developed. The potential MoT-specific genes were then identified by analyzing the gene presence/absence from the global pair-wise gene sequence alignment among the two MoT strains(PR01-37.V.1 and PR01-37.V.3)and three MoO strains.The potential specific genes were then aligned to the additional MoO and MoT genomes [4,19]. Specifically, a two-step approach was used to identify the MoT specific genes. First, we aligned the 378 candidate genes to National Center for Biotechnology Information Search database nr databases to remove the MoO genes. Second, we aligned the 100 candidate MoT genes to more sequenced MoT and MoO genomes to test whether they were present in all of the MoT genomes but absent in all of the MoO genomes.Finally,two MoT-specific genes were identified in the analysis. The detailed bioinformatic analysis procedure is shown in Fig. 1.

    2.4. Loop-mediated isothermal amplification

    The LAMP reaction mixtures were prepared to a final volume of 20 μL and contained 1 μL of template DNA (20 μmol·L-1),2.0 μmol·L-1of the forward and backward inner primers (FIP and BIP), 0.5 μmol·L-1F3 and B3 primers (Table 1), 1 mmol·L-1of deoxy-ribonucleoside triphosphate (dNTP), 8 U (1U = 16.67 nkat)of Bst DNA polymerase (New England Biolabs, USA), and 1× ThermoPol reaction buffer (20 mmol·L-1Tris-HCL, 10 mmol·L-1KCl,2 mmol·L-1MgSO4, 10 mmol·L-1(NH4)2SO4, and 0.1% Triton X-100, pH 8.8). The mixture was kept at 65 °C for 1 h before the reaction was terminated by moving the mixture to 85°C for 5 min.

    PrimerSequence (from 5′ to 3′)TargetReferenceTm (°C)Mo.T6098 FACCAATATCACCTGAACGCAGACATMoTFor PCR58.4 Mo.T6098 RGATTCCAGATTCACCACCAAAACAGMoTFor PCR56.4 Mo.T6099 FTCTGTATTTCACACTTGGGCTTTGGMoTFor PCR57.5 Mo.T6099 RAACGTCATGTAGTGCGTCTTGTTGAMoTFor PCR59.1 Pot2 FCGTCACACGTTCTTCAACCMoT & MoOFor PCR53.6 Pot2 RCGTTTCACGCTTCTCCGMoT & MoOFor PCR53.1 Mo.T6098-SFaTTTCTGCTCGTTGGGGAGACMoTFor RT-PCR56.4 Mo.T6098-SRaATGTACCCCCGTACCCATCAMoTFor RT-PCR56.3 Mo.T6099-SFaATGTTGGTAACCACCGGTCCMoTFor RT-PCR56.3 Mo.T6099-SRaCTTCGTATGCCTGCCGTTTGMoTFor RT-PCR56.4 Pot2-SFaACGACCCGTCTTTACTTATTTGGMoT & MoOFor RT-PCR54.8 Pot2-SRaAAGTAGCGTTGGTTTTGTTGGATMoT & MoOFor RT-PCR55.4 Fam-2229-F5′-[FAM-dT]CGTCACATAGGGTATATCTTTTCACAACAGssDNAFor RPA57.6 2229-nof-probAACCTACACTCGACTCACCGCAGTTCCCC AG (THF) CGCAGTATCGTAGCCG, 3′, C3 SpacerssDNAFor RPA74.6 5B-2229-R5′-BIOTIN-TTGTCATTCGCGGTCTTAGAA GCACTAACTAssDNAFor RPA61.3 P6098_370-gRNAAAUUUCUACUGUUGUAGAU CCGAUUUGCUGGUCGAAGCAUGUCas12a-bindingFor Cas12a—Mo.T6098F3GCCCTGCATCATGTACCAATMoTFor LAMP54.7 Mo.T6098B3ACTCGTCGGCATCTGTCAMoTFor LAMP55.1 Mo.T6098 FIPAGCCCTGTCTCCGCTTGTGTCACCTGAAC GCAGACATCGMoTFor LAMP72.4 Mo.T6098 BIPATGTCATTGACGGCGTCGCCTGACGTCGA TTCCCCAGTMoTFor LAMP72.5 Mo.T6099F3CGCATCTGGAGGGGAACTMoTFor LAMP55.4 Mo.T6099B3CTATCACCTGGCGGAGCTMoTFor LAMP55.2 Mo.T6099 FIPAATCGGTGGCTGGCAGAGACGTCTCTGCA CACGGTCGAMoTFor LAMP73.3 Mo.T6099 BIPGGTCGCTAACGTGGTGACGGTTCTCGGTG GCGAGGATGMoTFor LAMP72.8

    Fig. 1. Pipeline for the identification of two MoT-specific sequences. The ellipses at the upper-left side of the figure represent different M. oryzae genomes and Gene1, 2, 3,4,..., N represent all of the predicted genes in the sequenced genomes. The arrows between two genes in two genomes indicate the bi-directional pairwise gene sequence alignments, the arched arrows between the strains represent the sequence alignments of individual genes in the two genomes. nr: non-redundant protein sequence.

    2.5. Detection of LAMP reaction products

    A 2 μL volume of 1:10 diluted SYBR Green I (10 000×, Invitrogen, USA) was added to the LAMP products. Reactions that turned green were recorded as positive, and those that remained brown were recorded as negative.The amplified LAMP reactions were also analyzed on a 2% agarose gel that was stained with ethidium bromide.The gel was photographed under ultraviolet light using a Gel Doc XR+ Imager system (Bio-Rad, USA).

    2.6. Recombinase polymerase amplification

    RPA was performed as previously described by Piepenburg et al.[14] using the kit (#WLN8203KIT) from AMP-Future Biothech Co.Ltd.(China).For the first RPA, MoT-6098/MoT-6099-specific primer pairs were designed,and a tetrahydrofuran(THF)probe was added as the Nfoase recognition site.The amplicon should contain at least one protospacer adjacent motif (PAM) site (four-base domain of ‘‘TTTA”) for Cas12a recognition. For the second RPA,single-stranded DNA (ssDNA)-specific primer pairs were designed and labeled with 5-carboxyfluorescein (5-FAM) at the 5′end and with biotin at the 3′end. For the RPA, the instructions in the RPA kit were followed. In brief, the primers were mixed with the RPA solution containing the recombinase, polymerase, and the buffers,following which 1 μL of sample DNA (or of Cas12a-digested solution) was added; the preparation (total volume = 50 μL) was then kept at 37–39 °C for 8–12 min.

    2.7. ssDNA design

    ssDNA design has two requirements: ①The ssDNA should be short (to facilitate rapid PCR detection) and ②there should be no identical sequence between the ssDNA and the target genome sequence. To meet these requirements, a computer program was used to randomly design 10 000 120-bp DNA sequences. These 120-bp sequences were then aligned to the genomes of wheat,MoT, and MoO genomes using the BLAST software [20], and the sequences with an E-value ≤e-10or more than 10 bp identical to the target were discarded. The following sequence was randomly selected among sequences that had no detectable similarity with wheat, MoT, and MoO sequences: ‘‘ACGTCGTCACA TAGGGTATATCTTTTCACAACAGTAACCTACACTCGACTCACCGCAGTT CCCCAGA CGCAGTATCGTAGCCGCGGATAGTTAGTGCTTCTAAGACCG CGAATGACAACA.”

    2.8. Cas12a-based detection

    For the RPA with Cas12a, a previously described method was used[16].The Cas12a protein was purchased from New England Bio-Labs (NEB#M0653S, USA). In brief, Cas12a was mixed with small guide RNA (sgRNA) (equimolar volume) and reaction buffer, and the Cas12a mixture was incubated for~10 min at room temperature.The RPA-amplified target DNA solution was then incubated in the Cas12a mixture at 37 °C for about 10 min to activate the Cas12a ssDNA digestion activity.Finally,the designed ssDNA and the activated Cas12a protein were combined for ssDNA digestion.

    2.9. NALFIA detection using PCRD strips

    Abingdon Health PCRD test cassettes (#FD51673, USA) were used for DNA visualization. In brief, 5 μL of the amplicon from the RPA reaction was mixed with 70 μL of PCRD extraction buffer,and the total volume (75 μL) was added to the sample well of the PCRD text cassette. The result was assessed after 3–5 min.

    3. Results

    3.1. Identification of the MoT-specific sequences through comparative genomics

    To investigate the genomic differences between MoT and MoO,we extracted DNA from two MoT strains, PR01-37.V.1 and PR01-37.V.3, which were collected from infected wheat plants in Brazil.We then re-sequenced the two strains using the Illumina sequencing platform(USA)with a genome coverage of about 150 coverage.After genome assembly, we obtained a ~38.5 megabase (Mb) genome sequence for PR01-37.V.1 and a ~40.0 Mb genome sequence for PR01-37.V.3. We also downloaded the genome sequences of three MoO strains, Guy11, P131, and Y34, from GenBank [21].Through pair-wise gene sequence alignments, we identified 378 MoT genes that were absent in the three MoO strains. Alignment of these 378 genes with the GenBank protein database (nr) led to the identification of 289 genes with high homology with those in the databases. Among them, 189 had hits of MoO genes, and the other 100 were MoT-specific candidate genes (Table S1 in Appendix A). Except for a few genes that were highly homologous with bacterial genes,most of the genes’hits are fungal genes,and more than half (57 of 100) of those hits belonged to six fungal species(Table S2 in Appendix A). When we aligned the 57 genes to other~20 sequenced MoT and MoO genomes [4,19], we identified two genes, MoT-6098 and MoT-6099, that were specific to MoT. MoT-6098 encodes a putative acid trehalase protein, while MoT-6099 encodes a duf341 family protein. Because MoT-6098 and MoT-6099 were present in all of the sequenced MoT genomes but not in any of the sequenced MoO genomes, we chose these two MoTspecific genes to develop MoT-specific markers.

    3.2. Sensitivity of the MoT-6098 and MoT-6099 primer pairs for distinguishing MoT from MoO

    To determine the specificity of the MoT-6098 and MoT-6099 sequences in the MoT and MoO strains,we designed a set of primers(MoT-6098F and MoT-6098R, and MoT-6099F and MoT-6099R, in Table 1),and used them in PCR reactions with DNA samples isolated from 53 diverse M. oryzae strains. The 53 strains were collected from 14 grass species in nine countries (Tables 2 and 3), and consisted 15 strains from rice, 19 strains from wheat, and 19 strains from possible alternative hosts, such as barley, oat, Setaria, Brachiaria, Lolium, Bromus, Eleusine, Eragrostis, Festuca, Setaria viridis,and Stenotaphrum (Tables 2 and 3). We PCR tested the 25 strains that were collected from Bangladesh and Brazil(Table 2)in the Urashima and Islam laboratories.Among them,two DNA samples from Fusarium graminearum(F.graminearum),which causes the Fusarium head blight of wheat, isolated in the Urashima laboratory were included to test the primer specificity between fungal species. In addition,we included four DNA samples from isolates infecting rice and Lolium as the controls,which were available in the Wang laboratory.The PCR analysis indicated that the two genes were present in all eight MoT strains,one barley strain,and one US Lolium strain,but not in any of the strains from rice, Setaria, Brachiaria, oat, and two F. graminearum (Fig. 2 and Table 2). Following this, we tested the 24 M.oryzae strains that were deposited in the USDA-ARS’s Foreign Disease-Weed Science Research Unit(Table 3). Except for two strains collected from Brachiaria species,all other strains contained the two genes. Among the 22 MoT-positive strains, 11 were collected from wheat and the rest were isolated from alternative hosts of MoT,consisting Lolium,Eleusine,Bromus,Eragrostis,Stenotaphrum secundatum, and Festuca, collected in Brazil, Bolivia, Paraguay,Japan, and the USA (Fig. S1 in Appendix A and Table 3).

    We then conducted a sensitivity assay to detect the efficiency of the two primer pairs. Five concentrations (0.01, 0.1, 1, 10, and 50 ng·μL-1) of the DNA from a culture of MoT strain BO12TA301 were used in the PCR reactions with the MoT-6098F/R and MoT-6099F/R primer pairs.The Pot2R/F primer pair was used as the control. The PCR results indicated that MoT-6099F/R generated a DNA band with 0.1 ng·μL-1DNA and MoT-6098F/R generated a DNA band with 10 ng·μL-1DNA(Fig.3(a)).We also tested the primer sensitivity using real-time PCR with the same DNA dilutions that were used in the conventional gel-based PCR. Reactions containing only 0.01 ng·μL-1of template DNA yielded a positive cycle threshold(Ct)value for both primer pairs,indicating that amplification of the target sequence by these primers is very sensitive(Figs.3(b)–(d)).

    3.3. Detection of MoT in infected plants

    To determine whether the MoT-6098F/R and MoT-6099F/R primer pairs could be used to detect MoT in infected wheat plants,we isolated DNA from leaves and panicles of wheat plants infected with the MoT strain BTJP4-5, in a growth chamber in the Bangabandhu Sheikh Mujibur Rahman Agricultural University inBangladesh. As controls, we also isolated DNA from naturally infected panicles of rice and mycelia of MoT isolate BTJP4-5 and MoO isolate of RB-13b. The PCR results indicated that the MoT-6098 and MoT-6099 related bands were amplified with DNA obtained from MoT-infected wheat leaves and panicles and MoT mycelia, but not from the DNA isolated from MoO-infected rice leaves or MoO mycelia (Fig. 4). These results demonstrated that the MoT-6098 and MoT-6099 primers can effectively detect MoT in infected wheat plants.

    Well #IsolateHostOriginMo.T6098Mo.T6099Pot2 1 BTGS-4-CWheatBangladesh+a++2 BTJP 4-1WheatBangladesh+++3 BTJP 4-5WheatBangladesh+++4 BTGPWheatBangladesh+++5 BTMP 182-2WheatBangladesh+++6 BTMP 1839-2WheatBangladesh+++7 BTMP 1845-3WheatBangladesh+++8 PR01 83WheatBrazil+++9 CVPNO2-17-CBarleyBrazil+++10RBSA 18-B-1RiceBangladesh-b-+11RBCH 1814-2RiceBangladesh--+12RBME 18-16-2RiceBangladesh--+13RBME 1819-3RiceBangladesh--+14RBDI 1831-5RiceBangladesh--+15RBRA 18 34-1RiceBangladesh--+16RBRA1836-3RiceBangladesh--+17RBTA 1847-3RiceBangladesh--+18RBTA 1711-1RiceBangladesh--+19Rb-13-BRiceBangladesh--+20SgCCA2aSetaria geniculateBrazil--+21WhGi02-01Fusarium from wheatBrazil--+22FSFusarium from wheatBrazil--+23VP03-17RiceBrazil--+24BpCCA3dBrachiaria plantagineaBrazil--+25TrAv01-17OatBrazil--+A BN-0050RiceAfrica--+B KJ201RiceRepublic of Korea--+C RB22RiceRepublic of Korea--+D Cl-6RiceColombia--+E PL-2LoliumUSA+++F Water---

    Well #IsolateHostOriginMo.T6098Mo.T6099Pot2 1 T-17WheatBrazil+a++2 BZ-45WheatBrazil+++3 BZ-19WheatBrazil+++4 T-42WheatBrazil+++5 BZ-6WheatBrazil+++6 T-2WheatBrazil+++7 T-25WheatBrazil+++8 BO12TA301WheatBolivia+++9 B-2WheatBolivia+++10P-13WheatParaguay+++11P-3-1WheatParaguay+++12P-28Bromus spp.Paraguay+++13P-29Bromus spp.Paraguay+++14EC4JEleusine coracanaJapan+++15U-203Eragrostis curvulaUSA+++16U-199Festuca arundinaceaUSA+++17U-116Green foxtailUSA+++18U-69Lolium multiflorumUSA+++19S-3SetariaUSA+++20U-100Setaria faberiUSA+++21PL-2LoliumUSA+++22A-2Stenotaphrum secundatumUSA+++23P-26Brachiaria spp.Paraguay-b-+24P-30Brachiaria spp.Paraguay--+A BN-0050RiceAfrica--+B KJ201RiceRepublic of Korea--+C RB22RiceRepublic of Korea--+D Cl-6RiceColombia--+E Water---

    Fig. 2. MoT-6098 and MoT-6099 sequences are presence in all tested MoT strains but absence in all tested MoO strains detected by PCR. Photographs of gels (1%) obtained with MoT and MoO samples that were amplified with (a) the MoT-6098F/R primer pair, (b) the MoT-6099F/R primer pair, and (c) the control Pot2F/2R primer pair. DNA sample information (lane 1–25 and a–f in the gels) is listed in Table 2.

    3.4. Use of LAMP to detect MoT

    Since LAMP is a specific, rapid, and efficient method for DNA detection under isothermal conditions, we used LAMP to detect the MoT-6098 and MoT-6099 fragments in MoT. The amplified DNA was visualized in both agarose gels and 1000× SYBR Green dye solutions. In the agarose gel assay, amplified bands were present in the reactions with the DNA from three MoT strains but not from five MoO strains (Fig. 5(a)). Consistent with the gel-based results, green color was obtained with the SYBR Green dye solutions containing the DNA of three MoT strains but not with the solutions containing the DNA of five MoO strains(Fig.5(b)). These results showed that LAMP could be used to detect MoT rapidly and efficiently under isothermal conditions.

    Fig.4. Detection of MoT in infected wheat plants.(a)Early stage of leaf infection(circle indicates a water-soaked lesion in a leaf 14 days after sowing).(b)Panicle infection(circle shows a gray spot in the neck of a spike 105 d after sowing).(c)PCR detection of MoT sequences in infected young and adult plants using agarose gel electrophoresis.Lane 1 is the DNA marker lane. Lanes 2 to 6 represent DNA isolated from the following samples: mycelia of MoT strain BTJP4-5(lane 2),wheat leaves artificially inoculated with MoT strain BTJP4-5(lane 3),the panicle neck of a naturally infected wheat plant in an MoT-infected wheat field(lane 4),mycelia of MoO strain RB-13b(lane 5),and the panicle neck of a naturally infected rice plant in an MoO-infected rice field (lane 6).

    Fig. 5. LAMP detection of MoT sequences. The LAMP assay was conducted with eight samples (labeled 1 to 8) containing MoT genomic DNA (samples 1 to 3) and MoO genomic DNA(samples 4 to 8).(a)LAMP-amplified genomic DNA visualized on a 2%agarose gel;the left side and the right side of the gel indicate results obtained with LAMP primers that targeted MoT-6098 and MoT-6099,respectively. (b) LAMP-amplified genomic DNA visualized with SYBR Green dye.The left side and the right side indicate the results obtained with LAMP primers that targeted MoT-6098 and MoT-6099, respectively.

    3.5.Integrating the Cas12a protein with RPA and NALFIA technologies for the rapid detection of MoT in the field

    Fig.6. Rapid detection of MoT by integrating the Cas12a protein with RPA and NALFIA technologies.(a)The Cas12a-mediated detection method when a DNA sample contains the MoT-6098 or MoT-6099 sequence.(b)The Cas12a-mediated detection method when a DNA sample does not contain the MoT-6098 or MoT-6099 sequence.(c)Detection of the amplified DNA from the first RPA on an agarose gel using the MoT-6098 primer pair. (d) NALFIA detection of the MoT-6098 sequence in MoO and MoT strains after the second RPA on the strips. (e) NALFIA detection of the MoT-6098 sequence using DNA isolated from mycelia of the MoT strain BTJP4-5 (lane 1), a wheat leaf artificially inoculated with MoT strain BTJP4-5 (lane 2), and the panicle neck of a wheat plant that was naturally infected with MoT in the field (lane 3). CTL-: negative control;CTL+: positive control.

    The application of the LAMP technology in field is limited because it requires high temperature conditions(60–65°C).Cas12a is a recently identified nuclease that can nonspecifically digest ssDNA or single-stranded RNA (ssRNA) after sgRNA-mediated DNA binding of the target sequence [12]. Based on Cas12a’s sequence recognition specificity and ssDNA nuclease activity, we used this protein in combination with the RPA and NALFIA technologies for rapid MoT detection (Figs. 6(a) and (b)). This process first uses RPA to amplify the MoT-specific sequences in MoT, and the amplified product is then incubated with sgRNA that targets MoT-6098 and MoT-6099 sequences and the Cas12a protein. If the DNA sample of MoT contains the target sequences, the nuclease activity of Cas12a will be activated and the Cas12a protein will nonspecifically digest the randomly designed ssDNA probes in the reaction. A second RPA then uses the FAM-labeled primer (left primer)and the biotin-labeled primer(right primer)for the ssDNA fragments, and the amplicon is loaded onto lateral flow strips(PCRD)for the detection of the ssDNA fragment.If the MoT sample contains the sgRNA target sequences, the activated Cas12a will digest the ssDNA, and thus, the biotin- and FAM-labeled fragment of ssDNA will not appear on the PCRD strip (Fig. 6(a)). In contrast,the ssDNA will be amplified in the second RPA and a biotin- and FAM-labeled ssDNA band will appear on the PCRD strip (Fig.6(b)). We selected four MoT and four MoO strains for the RPA and NALFIA assays and included a 120-bp randomly designed DNA fragment (see experimental procedures in Appendix A) for Cas12a digestion. The agarose gel analysis confirmed the amplification of the MoT-6098 and MoT-6099 fragments in the first RPA reaction for the four MoT samples, but not for the four MoO samples(Fig.6(c)).After the second RPA,the reaction solutions were loaded onto the PCRD strips for NALFIA detection.The results showed that the ssDNA band(the second band from top)was clearly evident in all four MoO samples and in the positive control,but not in the four MoT samples or in the negative control(Fig.6(d)).

    Fig.7. Sensitivity assays of the Cas12a-RPA-NALFIA-mediated method and its comparison with traditional PCR.(a)Sensitivity of PCR to detect the MoT-6098 sequence using different concentration of template DNA from mycelia on an agarose gel. (b) Sensitivity of PCR to detect the MoT-6098 sequence in the inoculated leaves on an agarose gel.(c) Sensitivity of the Cas12a-RPA-NALFIA-mediated method to detect the MoT-6098 sequence from mycelia on PCRD strips. (d) Sensitivity of the Cas12a-RPA-NALFIAmediated method to detect the MoT-6098 sequence in the inoculated leaves on PCRD strips. DPI: days post inoculation.

    Furthermore, when used with the DNA isolated from MoTinfected wheat leaves,the NALFIA assay detected MoT in both artificially infected leaves and naturally infected panicles(Fig.6(e)).To further test the sensitivity of the Cas12a-based NALFIA method,we compared it with the conventional PCR method using different concentrations of DNA from fungal mycelia (Fig. 7). As shown in Figs. 7(a) and (c), a week band was observed when the fungal DNA concentration was 0.01 μg·μL-1.In contrast,the Cas12a-based NALFIA method clearly detected the target sequence in as little as 0.001 μg·μL-1of the fungal DNA. When DNA from the inoculated leaves was used, PCR detected the target sequence four days post-inoculation (Fig. 7(b)). However, the Cas12a-based NALFIA method detected the target sequence as early as two days postinoculation (Fig. 7(d)). These results demonstrate that the new detection method can rapidly, accurately, and efficiently identify MoT in wheat fields, thereby providing a simple and economical platform for the detection of the devastating wheat blast disease,in South American and South Asian countries.

    4. Discussion

    With the development of NGS technology,researchers can now rapidly and inexpensively sequence an individual genome.Although NGS technology enables researchers to identify genomic differences among closely related organisms, identifying sequences that are unique to a specific species or pathotype remain challenging. To identify sequences in MoT that are not present in MoO, we re-sequenced two MoT strains from Brazil. Using a pairwise, two-step global alignment method, we aligned the potentially different genes in the two MoT strains with sequenced MoO and MoT strains to identify MoT-specific genes. Among the 57 MoT-unique candidate genes, we found two sequences that were present in MoT and other M. oryzae strains collected from alternative hosts, but not in any MoO strains. Our PCR analysis with 53 M. oryzae strains confirmed the results and demonstrated that the combination of NGS sequencing with a two-step global alignment method swiftly enables accurate identification of unique sequences when comparing two highly similar genomes.

    Rapid and accurate diagnostic assays of emergent and reemergent diseases are required to control crop diseases. In some cases, disease control depends on the immediate identification of the causal pathogens in the field. Immunological-based detection(IBD) techniques, such as enzyme-linked immunosorbent assay,have been used for on-site pathogen detection since the 1970s.As the antibodies used in the conventional IBD method are not stable, the method now uses modified colloidal gold or other substitute particles that are activated and therefore visualized by specific antibodies.The limitation of the IBD method is that specific antibodies are unavailable for most plant pathogens. Over the last four decades,PCR and RT-PCR have been used extensively to detect pathogen DNA or RNA sequences in infected plants[22].The disadvantages of the PCR method are the requirement of a laboratory with a PCR machine and the relatively long time (4–8 h) required to obtain the diagnostic result. In contrast, NALFIA is a paperbased platform for the detection and quantification of analytes in complexes, providing results within 5–30 min [23]. The low cost and ease of the assay have resulted in the expansion of the NALFIA application in multiple fields, including plant pathogen detection.The combination of NALFIA with nucleotide detection methods,such as LAMP and RPA, has recently become a popular way to detect plant pathogens [24]. Another diagnostic method uses the Cas12a protein, which has indiscriminate ssDNA cleavage activity that completely degrades ssDNA molecules after binding to the target sequence through sgRNA [16]. By combining Cas12a singlestranded deoxyribonuclease activation with isothermal amplification, a new method called DNA endonuclease-targeted clustered regularly interspaced short palindromic repeats trans reporter was developed for the rapid and specific detection of a virus in human samples [16].

    In this study,we first determined the sensitivity of primers with PCR and LAMP at isothermal conditions. Following which, we included the Cas12a protein and a random DNA fragment in the first RPA reaction to amplify the target sequences in MoT with MoT-6098 or MoT-6099 primer pairs, at 37 °C. In the second RPA reaction, biotin- and FAM-labeled primers were used to amplify the random DNA fragment for PCRD strip detection.Positive results were obtained with DNA isolated from MoT mycelial tissues and from MoT-infected wheat leaves. The entire procedure required only about 0.5 h and was conducted at 25–40°C without the need for any major laboratory equipment.The total cost of the assay per sample was about 6 USD when commercial PCRD strips were used.The cost per sample can be reduced to less than 4 USD if the strips are manufactured on a large scale using non-patented methods.Therefore,the simple method developed in this study can be easily adapted by plant pathologists, quarantine specialists, and agriculture extension workers to rapidly detect MoT in infected plants,seed lots, and alternate hosts.

    Acknowledgements

    This research was supported by grants from the SEEDS program of the Ohio Agricultural Research and Development Center(OARDC) and the Coordinated Research Program of the International Atomic Energy Agency to Guo-Liang Wang. The financial support from the Krishi Gobeshona Foundation of Bangladesh under project KGF TF 50-C/17 to Tofazzal Islam of the Institute of Biotechnology and Genetic Engineering of BSMRAU, Bangladesh,is highly acknowledged. Assistance from Dr. Francesca Hand for the LAMP assay is appreciated. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the US Department of Agriculture. USDA is an equal opportunity provider and employer.

    Data availability statement

    The data that support the findings of this study have being submitted to GenBank and the accession numbers are JAAXMV000000000 and JAAXMU000000000.

    Compliance with ethics guidelines

    Houxiang Kang, Peng Ye, Kangyu Hua, Yufei Deng, Maria Bellizzi, Dipali Rani Gupta, Nur Uddin Mahmud, Alfredo S.Urashima, Sanjoy Kumar Paul, Gary Peterson, Yilin Zhou,Md Tofazzal Islam, Xueping Zhou, and Guo-Liang Wang 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.2020.07.016.

    av天堂在线播放| 国产亚洲欧美在线一区二区| 亚洲精品中文字幕一二三四区| 成人黄色视频免费在线看| 欧美国产精品一级二级三级| 热99国产精品久久久久久7| 18禁观看日本| 精品人妻熟女毛片av久久网站| 天天操日日干夜夜撸| 成人国语在线视频| 国产精品乱码一区二三区的特点 | 最近最新中文字幕大全电影3 | 一个人免费在线观看的高清视频| 久久午夜亚洲精品久久| 亚洲色图综合在线观看| 女人被躁到高潮嗷嗷叫费观| 欧美激情极品国产一区二区三区| 亚洲熟女精品中文字幕| 啦啦啦视频在线资源免费观看| 国产亚洲精品久久久久久毛片 | 人人澡人人妻人| 在线观看66精品国产| 久久国产精品男人的天堂亚洲| 极品教师在线免费播放| 老司机午夜福利在线观看视频| 黑人操中国人逼视频| ponron亚洲| 99re6热这里在线精品视频| 久久久久国内视频| 人人妻人人添人人爽欧美一区卜| 三上悠亚av全集在线观看| 国产精品二区激情视频| 午夜精品久久久久久毛片777| 两个人看的免费小视频| 亚洲精品国产区一区二| 国产一区二区激情短视频| 制服诱惑二区| 丁香六月欧美| 久久人人97超碰香蕉20202| 日韩三级视频一区二区三区| 久99久视频精品免费| 久久青草综合色| 国产一区二区三区视频了| 人妻一区二区av| 大香蕉久久成人网| 高潮久久久久久久久久久不卡| 黄网站色视频无遮挡免费观看| 丰满的人妻完整版| 制服诱惑二区| 亚洲色图综合在线观看| 日本五十路高清| 久久国产精品大桥未久av| 三级毛片av免费| 久久国产亚洲av麻豆专区| 久久久久久人人人人人| 人人妻人人澡人人爽人人夜夜| 国产成人av教育| 91精品国产国语对白视频| 午夜视频精品福利| 天天添夜夜摸| 国产区一区二久久| 9热在线视频观看99| 精品视频人人做人人爽| 国产精品偷伦视频观看了| 制服人妻中文乱码| 精品国产亚洲在线| 亚洲熟妇中文字幕五十中出 | 亚洲精品久久午夜乱码| 国产精品影院久久| 精品国产乱码久久久久久男人| 久久久久视频综合| 少妇 在线观看| 日韩欧美一区二区三区在线观看 | 99国产精品一区二区三区| 国产精品久久久av美女十八| 麻豆av在线久日| 久久ye,这里只有精品| 亚洲午夜精品一区,二区,三区| 美女午夜性视频免费| 51午夜福利影视在线观看| 欧美成人免费av一区二区三区 | 黑丝袜美女国产一区| bbb黄色大片| 精品午夜福利视频在线观看一区| 久久国产精品男人的天堂亚洲| 亚洲精品中文字幕在线视频| 激情视频va一区二区三区| 动漫黄色视频在线观看| 亚洲少妇的诱惑av| 国产精品国产av在线观看| 老汉色av国产亚洲站长工具| 久久久久国产精品人妻aⅴ院 | 亚洲国产毛片av蜜桃av| 99久久人妻综合| 精品国产美女av久久久久小说| 中文字幕人妻丝袜制服| 成熟少妇高潮喷水视频| 日本黄色视频三级网站网址 | 9191精品国产免费久久| 国产主播在线观看一区二区| 久久香蕉激情| 在线观看舔阴道视频| 日韩成人在线观看一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 国产精品 欧美亚洲| 亚洲中文字幕日韩| 久久精品人人爽人人爽视色| 久久久久视频综合| 一级a爱视频在线免费观看| 一进一出好大好爽视频| 99热国产这里只有精品6| 黑人操中国人逼视频| 欧美精品高潮呻吟av久久| 狠狠狠狠99中文字幕| 欧美中文综合在线视频| 黄网站色视频无遮挡免费观看| 老汉色av国产亚洲站长工具| 久久久久国产一级毛片高清牌| 成年版毛片免费区| 高清视频免费观看一区二区| 十分钟在线观看高清视频www| 国产精品乱码一区二三区的特点 | 午夜福利视频在线观看免费| 亚洲av熟女| 男女高潮啪啪啪动态图| 亚洲专区中文字幕在线| 在线看a的网站| 大片电影免费在线观看免费| 黄网站色视频无遮挡免费观看| 在线天堂中文资源库| 精品无人区乱码1区二区| 高清欧美精品videossex| 美女午夜性视频免费| 这个男人来自地球电影免费观看| 亚洲成人手机| 精品一区二区三卡| 国产极品粉嫩免费观看在线| 亚洲精品一卡2卡三卡4卡5卡| 亚洲久久久国产精品| 黄色片一级片一级黄色片| 亚洲一区二区三区不卡视频| 亚洲,欧美精品.| 别揉我奶头~嗯~啊~动态视频| 丝袜美足系列| 久久精品熟女亚洲av麻豆精品| 欧美日韩成人在线一区二区| 久久精品国产综合久久久| 国产又爽黄色视频| 如日韩欧美国产精品一区二区三区| 黄色毛片三级朝国网站| 美女午夜性视频免费| 高清毛片免费观看视频网站 | 色老头精品视频在线观看| 一进一出好大好爽视频| 夜夜夜夜夜久久久久| 国产乱人伦免费视频| 亚洲久久久国产精品| 日韩一卡2卡3卡4卡2021年| 免费在线观看视频国产中文字幕亚洲| 一区福利在线观看| 久久这里只有精品19| 在线十欧美十亚洲十日本专区| av视频免费观看在线观看| 久久精品91无色码中文字幕| 精品熟女少妇八av免费久了| 91麻豆av在线| 国产在视频线精品| 巨乳人妻的诱惑在线观看| a在线观看视频网站| 成人精品一区二区免费| 精品国产美女av久久久久小说| 在线观看免费午夜福利视频| av免费在线观看网站| 天堂√8在线中文| 欧美黄色淫秽网站| 丝袜在线中文字幕| 交换朋友夫妻互换小说| 国精品久久久久久国模美| 精品亚洲成a人片在线观看| 夜夜爽天天搞| 9191精品国产免费久久| 亚洲,欧美精品.| 老司机午夜十八禁免费视频| 亚洲中文字幕日韩| 国产成人精品久久二区二区91| 制服人妻中文乱码| 国产又爽黄色视频| 精品国产一区二区久久| 欧美日韩精品网址| 大陆偷拍与自拍| 法律面前人人平等表现在哪些方面| 黄色 视频免费看| 亚洲精品中文字幕一二三四区| 久久人人爽av亚洲精品天堂| 天天影视国产精品| 黄色片一级片一级黄色片| 在线十欧美十亚洲十日本专区| 中文字幕人妻丝袜一区二区| 熟女少妇亚洲综合色aaa.| 国产精品亚洲一级av第二区| 日本a在线网址| 淫妇啪啪啪对白视频| tube8黄色片| 久久精品91无色码中文字幕| 午夜久久久在线观看| 少妇猛男粗大的猛烈进出视频| 女人久久www免费人成看片| 黄色片一级片一级黄色片| 久久精品国产综合久久久| 午夜视频精品福利| videosex国产| 日本vs欧美在线观看视频| 午夜免费鲁丝| 日韩中文字幕欧美一区二区| 国产精品久久久人人做人人爽| 母亲3免费完整高清在线观看| av片东京热男人的天堂| 国产片内射在线| 一本大道久久a久久精品| 国产深夜福利视频在线观看| 人人妻人人澡人人爽人人夜夜| 三级毛片av免费| 黑人操中国人逼视频| 男女午夜视频在线观看| 丝袜美腿诱惑在线| 成年女人毛片免费观看观看9 | 亚洲精品成人av观看孕妇| 国产不卡av网站在线观看| 国产亚洲精品久久久久5区| 极品人妻少妇av视频| 欧美乱码精品一区二区三区| www.自偷自拍.com| 亚洲av日韩精品久久久久久密| 美女 人体艺术 gogo| 精品久久久精品久久久| 亚洲国产精品sss在线观看 | 亚洲熟妇中文字幕五十中出 | 一本综合久久免费| 中文字幕人妻丝袜一区二区| 色尼玛亚洲综合影院| 男人舔女人的私密视频| 亚洲五月婷婷丁香| 亚洲一区中文字幕在线| 欧美日韩亚洲综合一区二区三区_| 黑人欧美特级aaaaaa片| 久久亚洲精品不卡| 两人在一起打扑克的视频| 日本黄色视频三级网站网址 | 日韩中文字幕欧美一区二区| a级毛片在线看网站| 天天躁日日躁夜夜躁夜夜| 男女床上黄色一级片免费看| 老司机午夜福利在线观看视频| 一进一出抽搐gif免费好疼 | www.精华液| 老汉色∧v一级毛片| 欧美日韩乱码在线| 黄色a级毛片大全视频| 啦啦啦免费观看视频1| 国产国语露脸激情在线看| 在线观看66精品国产| 少妇裸体淫交视频免费看高清 | 日韩免费高清中文字幕av| 天堂√8在线中文| 一级片'在线观看视频| 嫩草影视91久久| 91国产中文字幕| 国产精品乱码一区二三区的特点 | 欧美黑人欧美精品刺激| 国产男女内射视频| 超碰97精品在线观看| 免费人成视频x8x8入口观看| 国产免费男女视频| 亚洲熟妇熟女久久| 久久香蕉激情| 国产精品自产拍在线观看55亚洲 | 精品久久久久久电影网| 久久国产亚洲av麻豆专区| 久久久久久久国产电影| 一夜夜www| 人人澡人人妻人| 欧美日韩亚洲高清精品| 日本a在线网址| 丰满人妻熟妇乱又伦精品不卡| 久久久久久人人人人人| 黄频高清免费视频| 欧美国产精品一级二级三级| 午夜视频精品福利| 51午夜福利影视在线观看| 脱女人内裤的视频| 免费看a级黄色片| 亚洲av日韩在线播放| av国产精品久久久久影院| 老熟女久久久| 国产午夜精品久久久久久| 大香蕉久久成人网| 亚洲欧美一区二区三区久久| 午夜福利乱码中文字幕| 操出白浆在线播放| 亚洲一区中文字幕在线| 国产三级黄色录像| 国产精品亚洲av一区麻豆| 亚洲精品在线美女| 亚洲国产精品一区二区三区在线| 成人18禁在线播放| 手机成人av网站| 在线观看免费视频网站a站| 男人的好看免费观看在线视频 | 一a级毛片在线观看| av超薄肉色丝袜交足视频| 亚洲熟女精品中文字幕| 久久热在线av| 国产视频一区二区在线看| 老熟女久久久| 夜夜爽天天搞| 正在播放国产对白刺激| 动漫黄色视频在线观看| 亚洲av日韩在线播放| 久久国产亚洲av麻豆专区| 欧美激情 高清一区二区三区| 中文字幕色久视频| 日韩免费av在线播放| 免费看十八禁软件| 色播在线永久视频| 精品一品国产午夜福利视频| av片东京热男人的天堂| 妹子高潮喷水视频| 免费观看人在逋| 精品国内亚洲2022精品成人 | 亚洲av熟女| 精品一区二区三区av网在线观看| 久久久久久人人人人人| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲九九香蕉| 大片电影免费在线观看免费| 黄色 视频免费看| 老汉色av国产亚洲站长工具| 欧美日韩视频精品一区| 亚洲国产中文字幕在线视频| 成人av一区二区三区在线看| 国产午夜精品久久久久久| 精品欧美一区二区三区在线| 热re99久久精品国产66热6| 少妇裸体淫交视频免费看高清 | 欧美激情高清一区二区三区| 91大片在线观看| tube8黄色片| 成人黄色视频免费在线看| 免费在线观看视频国产中文字幕亚洲| 久久中文看片网| 成人国语在线视频| 国产99白浆流出| 久久精品国产亚洲av高清一级| 亚洲精品久久午夜乱码| 一级a爱片免费观看的视频| 丝袜人妻中文字幕| 又黄又粗又硬又大视频| 国产精品.久久久| 成人18禁在线播放| 一级a爱片免费观看的视频| 亚洲免费av在线视频| 国产男靠女视频免费网站| 美女 人体艺术 gogo| 亚洲专区中文字幕在线| 很黄的视频免费| 久久精品国产a三级三级三级| 美女 人体艺术 gogo| 日韩制服丝袜自拍偷拍| 精品久久久久久电影网| 精品久久久久久,| 少妇的丰满在线观看| 丁香欧美五月| www.精华液| 人妻 亚洲 视频| 午夜久久久在线观看| 久久精品aⅴ一区二区三区四区| 国产精品久久电影中文字幕 | 99久久精品国产亚洲精品| 建设人人有责人人尽责人人享有的| 午夜亚洲福利在线播放| 丰满的人妻完整版| 热99re8久久精品国产| 精品少妇久久久久久888优播| 日韩欧美国产一区二区入口| 亚洲男人天堂网一区| 国产97色在线日韩免费| 777久久人妻少妇嫩草av网站| 久久精品国产清高在天天线| 午夜精品久久久久久毛片777| av免费在线观看网站| 久久久国产精品麻豆| 欧洲精品卡2卡3卡4卡5卡区| 久久人人97超碰香蕉20202| 色精品久久人妻99蜜桃| 啦啦啦视频在线资源免费观看| 亚洲视频免费观看视频| 女警被强在线播放| 午夜激情av网站| 亚洲av第一区精品v没综合| 男男h啪啪无遮挡| 99热只有精品国产| 成熟少妇高潮喷水视频| 国产深夜福利视频在线观看| 嫁个100分男人电影在线观看| 亚洲国产欧美网| 99re6热这里在线精品视频| 首页视频小说图片口味搜索| 婷婷成人精品国产| 亚洲av第一区精品v没综合| 国产精品一区二区在线不卡| 天天躁夜夜躁狠狠躁躁| 亚洲aⅴ乱码一区二区在线播放 | 日韩熟女老妇一区二区性免费视频| www.熟女人妻精品国产| 国产成人系列免费观看| 久9热在线精品视频| 国产精品98久久久久久宅男小说| 国产av精品麻豆| 一级毛片女人18水好多| 国产av一区二区精品久久| 18禁裸乳无遮挡动漫免费视频| 三上悠亚av全集在线观看| 久久精品亚洲av国产电影网| 水蜜桃什么品种好| 精品国产乱码久久久久久男人| 丁香六月欧美| 亚洲精品美女久久久久99蜜臀| 午夜成年电影在线免费观看| 无遮挡黄片免费观看| 精品熟女少妇八av免费久了| 大型黄色视频在线免费观看| 亚洲一区高清亚洲精品| 黄色女人牲交| 亚洲精品久久午夜乱码| 精品国产一区二区三区久久久樱花| 国产精品久久久久成人av| 国产成人精品久久二区二区91| 黄色成人免费大全| 欧美色视频一区免费| 老熟女久久久| 校园春色视频在线观看| 深夜精品福利| 亚洲欧美一区二区三区久久| av网站免费在线观看视频| 乱人伦中国视频| 很黄的视频免费| 在线观看一区二区三区激情| 国产人伦9x9x在线观看| 亚洲熟妇熟女久久| 又黄又粗又硬又大视频| 一边摸一边抽搐一进一小说 | 欧美黄色淫秽网站| 黄色女人牲交| 啦啦啦在线免费观看视频4| 国产精品一区二区在线观看99| 成人永久免费在线观看视频| 久久天躁狠狠躁夜夜2o2o| 国产99白浆流出| x7x7x7水蜜桃| 久久亚洲精品不卡| 看免费av毛片| 欧美日韩乱码在线| 亚洲欧美激情在线| 国产麻豆69| 亚洲第一av免费看| 两个人看的免费小视频| 国产男女超爽视频在线观看| 欧美日韩精品网址| 一级作爱视频免费观看| 国产野战对白在线观看| 多毛熟女@视频| 高清视频免费观看一区二区| 韩国av一区二区三区四区| 69av精品久久久久久| 精品国产超薄肉色丝袜足j| 国产极品粉嫩免费观看在线| 国产一区在线观看成人免费| 我的亚洲天堂| 国产有黄有色有爽视频| 麻豆乱淫一区二区| 很黄的视频免费| av欧美777| 怎么达到女性高潮| 狠狠婷婷综合久久久久久88av| 国产精品免费视频内射| 99香蕉大伊视频| 国产一区有黄有色的免费视频| 成人18禁在线播放| 久久久久久久午夜电影 | 成年动漫av网址| 精品久久久久久久久久免费视频 | 久久天躁狠狠躁夜夜2o2o| 女人被躁到高潮嗷嗷叫费观| 亚洲成国产人片在线观看| 国产免费av片在线观看野外av| 精品国内亚洲2022精品成人 | 精品一区二区三区av网在线观看| 国产高清激情床上av| 操美女的视频在线观看| 高清av免费在线| 十八禁高潮呻吟视频| 热99re8久久精品国产| 亚洲成人免费电影在线观看| 国产av精品麻豆| xxxhd国产人妻xxx| 日韩三级视频一区二区三区| 精品少妇久久久久久888优播| 高清视频免费观看一区二区| 一级毛片女人18水好多| 19禁男女啪啪无遮挡网站| 久久精品国产综合久久久| 日韩视频一区二区在线观看| 91在线观看av| 国产精品一区二区在线观看99| 一本综合久久免费| 中文字幕另类日韩欧美亚洲嫩草| 精品国产国语对白av| 视频区欧美日本亚洲| 黄色丝袜av网址大全| 久久久精品免费免费高清| 国产精品久久电影中文字幕 | 日本精品一区二区三区蜜桃| 久久久久久久久免费视频了| 国产蜜桃级精品一区二区三区 | 欧美日韩亚洲高清精品| av网站免费在线观看视频| 欧美不卡视频在线免费观看 | 久久精品国产清高在天天线| 搡老乐熟女国产| 捣出白浆h1v1| 黄色a级毛片大全视频| 亚洲成av片中文字幕在线观看| 欧美最黄视频在线播放免费 | 亚洲精品在线观看二区| x7x7x7水蜜桃| 亚洲第一av免费看| 黑人猛操日本美女一级片| av视频免费观看在线观看| 国产成人啪精品午夜网站| 久久精品国产清高在天天线| 黄色丝袜av网址大全| 久久精品国产清高在天天线| 黄色丝袜av网址大全| 欧美老熟妇乱子伦牲交| 国产人伦9x9x在线观看| 久久国产精品男人的天堂亚洲| 大码成人一级视频| 欧美在线黄色| 色婷婷av一区二区三区视频| 免费看a级黄色片| 国产精品亚洲av一区麻豆| 精品国产一区二区三区四区第35| 国产激情久久老熟女| 热99久久久久精品小说推荐| 国产在线精品亚洲第一网站| 日本黄色日本黄色录像| 一区二区三区国产精品乱码| 俄罗斯特黄特色一大片| 亚洲国产精品合色在线| 欧美精品人与动牲交sv欧美| 国产亚洲av高清不卡| 男人操女人黄网站| 老熟妇仑乱视频hdxx| aaaaa片日本免费| 国产精品久久久av美女十八| √禁漫天堂资源中文www| 国产不卡av网站在线观看| 看免费av毛片| 大陆偷拍与自拍| 欧美 亚洲 国产 日韩一| 免费人成视频x8x8入口观看| 在线观看免费视频日本深夜| tocl精华| av天堂在线播放| 一级片'在线观看视频| 午夜福利一区二区在线看| 国产成人av激情在线播放| 欧美激情极品国产一区二区三区| 中文字幕高清在线视频| 一夜夜www| 一区二区三区激情视频| 一区在线观看完整版| 国产野战对白在线观看| 国产精品 欧美亚洲| 日本黄色日本黄色录像| 国产在视频线精品| 亚洲五月天丁香| 国产无遮挡羞羞视频在线观看| 久久这里只有精品19| 国产精品香港三级国产av潘金莲| 黄片播放在线免费| 久久国产精品大桥未久av| 亚洲欧美精品综合一区二区三区| 国产精品偷伦视频观看了| 女人高潮潮喷娇喘18禁视频| av一本久久久久| 91成年电影在线观看| 亚洲色图av天堂| 69精品国产乱码久久久| 天天躁夜夜躁狠狠躁躁| 免费久久久久久久精品成人欧美视频| 捣出白浆h1v1| 成年人午夜在线观看视频| 免费观看a级毛片全部| 巨乳人妻的诱惑在线观看| 欧美大码av| 大陆偷拍与自拍| 亚洲熟女毛片儿| 日本vs欧美在线观看视频| 久久天堂一区二区三区四区| 亚洲av成人一区二区三| 中亚洲国语对白在线视频| 国产高清视频在线播放一区| 久久草成人影院| 欧美黄色片欧美黄色片|