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

    Resistance to powdery mildew in the pea cultivar Xucai 1 is conferred by the gene er1

    2015-12-21 07:47:35SuliSunZhongyiWangHainingFuCanxingDuanXiaomingWangZhendongZhu
    The Crop Journal 2015年6期

    Suli Sun,Zhongyi Wang,Haining Fu,Canxing Duan,Xiaoming Wang,Zhendong Zhu*

    National Key Facility for Crop Gene Resources and Genetic Improvement,Institute of Crop Science,Chinese Academy of Agricultural Sciences, Beijing 100081,China

    Resistance to powdery mildew in the pea cultivar Xucai 1 is conferred by the gene er1

    Suli Sun1,Zhongyi Wang1,Haining Fu,Canxing Duan,Xiaoming Wang,Zhendong Zhu*

    National Key Facility for Crop Gene Resources and Genetic Improvement,Institute of Crop Science,Chinese Academy of Agricultural Sciences, Beijing 100081,China

    A R T I C L E I N F O

    Article history:

    Accepted 6 August 2015

    Available online 17 August 2015

    Erysiphe pisi er1-2 Pea Powdery mildew Xucai1

    Powdery mildew,caused by Erysiphe pisi D.C.,is a major constraint to pea production worldwide.The pea cultivar Xucai1 has shown high resistance to E.pisi under greenhouse and field conditions.The objectives of this study were to identify and characterize genes conferring resistance to powdery mildew in Xucai 1.Three crosses,Qizhen 76×Xucai 1, Bawan 6×Xucai 1,and Xucai 1×Bawan 6,were made to generate populations for genetic analysis.The resistance to E.pisi and segregation ratios in the F1,F2,and F2:3populations suggested a single recessive gene conferring the resistance of Xucai 1.Bulked segregant analysis was used to map the resistance gene using two F2populations.The resistance gene was close to markers AD60 and c5DNAmet on linkage group VI with genetic distances of 9.9 cM and 15.4 cM in the Xucai 1×Bawan 6 F2population and 8.7 cM and 8.1 cM in the Qizhen 76×Xucai1 F2population,respectively,suggesting that the resistance gene was aner1 allele.This hypothesis was confirmed by comparison of the cDNA sequences of the

    PsMLO1 gene between the parents and the PsMLO1 wild type.Three distinct types of transcripts in Xucai 1,characterized by a 129-bp deletion and 155-and 220-bp insertions, were detected,consistent with the structure of the er1-2 allele.We concluded that resistance in Xucai 1 was conferred by er1-2 and that its linked markers will be useful in pea breeding programs.

    ?2015 Crop Science Society of China and Institute of Crop Science,CAAS.Production and hosting by Elsevier B.V.This is an open access article under the CC BY-NC-ND license

    (http://creativecommons.org/licenses/by-nc-nd/4.0/).

    1.Introduction

    Pea(Pisum sativum L.)is widely distributed and has been cultivated in China for more than two thousand years[1]. According to the latest report,China ranks second only to Canada in dry pea production[2].Pea frequently suffers from biotic and abiotic stresses throughout the lifecycle of the plant.Powdery mildew,caused by Erysiphe pisi D.C.,is one of the most important pea diseases worldwide.The disease severely affects yield and quality of peas[3,4].It is of serious concern in temperate and subtropicalregions with warm,dry days and cool nights[4].Powdery mildew caused by E.pisi results in 25-50%yield losses in pea production[5-7].Severe infection by E.pisi can result in yield loss of up to 80%in susceptible cultivars[8].In pea production,agronomic practices,chemical prevention,and genetic resistance strategieshave been applied to controlthe disease.However,the control efficacy of chemical and agronomic practices is limited by many factors.Genetic resistance has become the first choice not only for its efficiency,low cost,and environmental friendliness[8,9]but also for qualitative disease resistance that are relatively easy to use in breeding programs.

    As early as 1948,Harland[10]identified powdery mildew resistance in the pea landrace Huancabamba and showed that it carried a single recessive gene.Since then,screening and genetic analysis of resistance to pea powdery mildew have been performed for more than 60 years[9].Many resistant pea accessions have been identified and characterized their gene(s)for resistance to E.pisi.Two recessive genes(er1 and er2)and one dominant gene(Er3)have been identified for powdery mildew resistance in pea germplasm[10-12].Genetic analyses of resistance to E.pisi indicated that the vast majority of resistant pea accessions carried the resistance gene er1[13].The resistance gene er2 is harbored by only a few resistant pea accessions[13].Er3 is a newly identified dominant gene from a wild relative of pea(P.fulvum)that has recently been successfully introduced into cultivated pea (P.sativum)[7,12].

    With respect to their resistance mechanisms,the modes of action at the cellular level of the two single recessive genes,er1 and er2,have been revealed by histological studies [14,15].er1 confers complete immunity or high levels of durable broad-spectrum resistance to powdery mildew by inhibiting E. pisi invasionofpea epidermalcells.The expression ofresistance gene er2 is influenced primarily by temperature and leaf age. The complete resistance conferred by er2 is revealed only at high temperature(25°C)or in mature leaves and is based mainly on post-penetration cell death,mediated by a hypersensitive response[14].The Er3 gene confers effective resistance to all E.pisi isolates tested to date[7].On Er3 plants,most E.pisi conidia are able to penetrate pea epidermalcells and form secondary hyphae,but growth of these established colonies is prevented by a strong hypersensitive response[12,16].

    DNA marker technologies and their application in genetic analysis have developed rapidly over the last decade.In peas, molecular markers linked to the three resistance genes(er1, er2,and Er3)for powdery mildew have been developed by several groups[17-25].The single recessive genes er1 and er2 have been located on pea linkage groups(LGs)VI and III, respectively[18,22,26],whereas Er3 has been mapped between the sequenced characterized amplified region(SCAR)marker Scw4637and random amplified polymorphic DNA(RAPD) marker OPAG05_1240,located on an uncertain pea LG[27].

    The recessively inherited resistance gene er1 has been used successfully for decades in pea breeding programs[13]. Recently,the resistance of er1 in pea has been found to be similar to that of MLO(mildew resistance locus O)genes in barley(Hordeum vulgare L.),Arabidopsis thaliana,and tomato (Lycopersicon esculentum Mill.)[28,29].The er1 resistant phenotype is caused by loss-of-function mutations of a pea MLO homolog,known as PsMLO1[28,29].Thus far,five er1 alleles (er1-1 to er1-5)have been characterized in pea resistant accessions,according to the mutation site and pattern of PsMLO1,each corresponding to a different PsMLO1 mutation [28-30].Among the reported five er1 alleles,er1-1 and er1-2 are used commonly in pea breeding programs[28,29].

    In China,powdery mildew is one of the most important diseases among more than 20 diseases infecting peas and severely affects the yield and quality ofpeas[31].The incidence of pea powdery mildew has reached 100%in some regions of Sichuan province[32].To date,the pathogen E.pisi is thought to be the only causal agent of pea powdery mildew in China [33-35],although two other species,E.baeumleri and E.trifolii, have been reported to cause powdery mildew on peas and showed the same symptoms as E.pisi in the Czech Republic[36] and the USA[37].Ond?ej et al.[36]and Attanayake et al.[37] showed that E.baeumleri and E.trifolii were able to overcome the er1 gene conferring resistance to E.pisi.Fondevilla et al.[38] further confirmed that E.trifolii had also overcome resistance gene Er3,but not er2.In China,screening for resistance to E. pisi in pea germplasms has been performed in greenhouses and in fields in recent years[32-35].Several pea cultivars/ landraces were highly resistant to E.pisi isolates EPBJ(NCBI, accession number KR912079)and EPYN(NCBI,accession number KR957355)[34,35].Both E.pisi isolates were identified by morphological and molecular characterization[39].However,the genetic basis of their resistance and the characteristics of their resistance gene(s)are unknown.

    Xucai 1,derived from a cross between Cuipimi and Baofeng 2,is the firstsemi-leafless pea cultivar(cv.)in China.Xucai1 has desirable traits including high yield and resistance to lodging and drought,root rot,and powdery mildew[40].Recently,Zeng et al.[33]showed that Xucai 1 was stably resistant to pea powdery mildew both in the greenhouse and in fields in 3 consecutive years.The present study was designed to clarify the genetic basis of this resistance and identify the resistance gene for powdery mildew present in Xucai1.With this aim,we performed an inheritance study using three crosses between resistant cv.Xucai1 and two susceptible cultivars(cvs.),Qizhen 76 and Bawan 6.We further characterized the resistance gene as e r1 by sequencing cDNA sequences of the PsMLO1 gene in the parents of the crosses and the pea cv.Sprinter[28],which is susceptible to E.pisi and carries the PsMLO1 wild-type allele.

    2.Materials and methods

    2.1.Plant materials and inoculum

    Three pea cultivars,Xucai 1(resistant to E.pisi),Qizhen 76 (susceptible to E.pisi),and Bawan 6(susceptible to E.pisi)were used in the experiment.Xucai 1 was kindly provided by Prof. Fengbao Wang of the Hebei Normal University of Science& Technology,Qinghuangdao,Hebei province,China;Bawan 6 was provided by Mr.Dongxu Xu of the Zhangjiakou Academy of Agricultural Sciences,Zhangjiakou,Hebei province,China; and Qizhen 76 was from Dr.Xin Chen of Jiangsu Academy of Agricultural Sciences,Nanjing,Jiangsu province,China.Three crosses were produced including reciprocal crosses between Bawan 6 and Xucai 1 and Qizhen 76×Xucai 1.Derived F1,F2, and F2:3populations of each cross were subjected to phenotypic evaluation and inheritance analysis.The F1plants from the crosses were self-pollinated to produce F2progeny.Each F2plant produced seeds for the F2:3families by self-pollination and was threshed individually.The two F2populations derived from the crosses of Qizhen 76×Xucai 1,and Xucai1×Bawan 6 and comprising 91 and 161 individuals,respectively,were used for genotypic evaluation and genetic analysis.

    The E.pisi isolate EPBJ(NCBIaccession number KR912079),a highly virulent isolate collected from a greenhouse in Beijing, was used as the inoculum in the present study[34,35].The isolate was maintained on the seedlings of susceptible pea cv.“Longwan 1.”Reproduction of inocula was achieved by continuous transfer to new healthy seedlings of pea cv.“Longwan 1”by gently shaking plants containing masses of E.pisi conidia. The inoculated plants were incubated in a growth chamber at 10±1°C with a 12 h light regime.

    2.2.Inoculation with E.pisi and disease assessment

    2.2.1.Parentalphenotypic evaluation

    Seeds of Xucai 1,Bawan 6,and Qizhen 76 were planted in 15-cm-diameter paper pots(six seeds per pot)filled with a mixture of vermiculite and peat moss(1:1).Each parent was planted in five replications using a total of 30 seeds per cultivar.The potted plants were placed under controlled conditions in a greenhouse at 18-26°C.Ten to fourteen days after planting,the plants at the fourth-or fifth-leaf stage were inoculated with E.pisi isolate EPBJ,using conidia shed by shaking heavily infected plants of pea cv.Longwan 1 over them.Ten days after inoculation,disease severity was scored on a 0-4 scale,according to the infected foliage area, macroscopic and microscopic density of mycelia and sporulation on the lower part of the plant[41,42].Plants with scores of 0-2 were classified as resistant(R)and those with scores of 3-4 as susceptible(S)[41,42].

    2.2.2.Phenotypic evaluation of F1and F2populations

    Planting was performed in a propagation greenhouse to generate offspring seeds.Plants were grown under the same conditions as described above for the parents.Inoculation was conducted under controlled conditions using the detached leaf method,together with susceptible and resistant parents as the controls.This inoculation method is effective for the evaluation of resistance to E.pisi in peas[14,16].The third or fourth leaf along with petiole from the bottom of the plant was detached from each F1and F2seedling at the fourth-or fifth-leaf stage.The cut leaves were inserted into 9-cm-diameter Petri dishes containing 2%water agar and 50 mg L-1benzimidazole to enhance longevity[43].The detached leaves were then inoculated with E.pisi isolate EPBJ using the method described above forthe inoculation ofthe parentalplants.After inoculation, the plates containing the inoculated leaves were sealed tightly with Parafilm and placed in a growth chamber at 20°C with a 14 h photoperiod.The detached leaves of the three parents Xucai 1,Qizhen 76,and Bawan 6 were treated in the same way with and without inoculation as controls.The disease severity was determined using the 0-4 scale 10 days after inoculation.The plants scored as resistant to E.pisi were subjected to repeated inoculation.

    2.2.3.Phenotypic evaluation of F2:3populations

    Twenty-four seeds selected randomly from each F2:3family were grown as described above.Disease was scored 10 days after inoculation using the 0-4 scale.F2:3families classified as either homozygous resistant or segregating for resistance to E.pisi were subjected to repeated inoculation.A chi-squared (χ2)test was used to determine the goodness-of-fit to Mendelian segregation patterns of phenotypes of all F2and F2:3populations.

    2.3.DNA extraction and pooling for bulk segregation analysis (BSA)

    Genomic DNAwas isolated fromhealthy young leaves collected from seedlings of the three pea parents,Xucai1,Qizhen 76,and Bawan 6,and of two F2mapping populations derived from Qizhen 76×Xucai 1,and Xucai 1×Bawan 6 using the cetyl trimethylammonium bromide(CTAB)extraction method with minor modifications[44].The DNA concentration was determined using an ND-1000 Nano-Drop instrument(NanoDrop Technologies Inc.,Wilmington,DE,USA).The DNAsolution was diluted to a working concentration of 50 ngμL-1with Tris-EDTA buffer(pH 8.0)and maintained at-20°C.For the bulked segregant analysis(BSA),resistant and susceptible bulks were prepared from equalamounts of DNA from 10 each of resistant and susceptible individual F2plants[45].

    2.4.Molecular marker analysis

    Initially 148 simple sequence repeat(SSR)markers that were approximately evenly distributed over a high-density pea genetic linkage map[46]were selected to screen for polymorphisms between the resistant and susceptible parents and pooled bulks. PCR amplification of SSR markers was performed in a total volume of 20μL,containing 50 ng of genomic DNA,2.5μL 10× PCR reaction buffer(20 mmol L-1MgCl2),0.2 mmol L-1of each dNTP,1.5 U of Taq DNA polymerase,and 0.2μmol L-1of primer mixture.PCR reactions were performed at 94°C for 5 min; followed by 35 cycles of 94°C for 45 s,49-60°C(depending on the primer-specific annealing temperature according to Loridon et al.[46])for 45 s and 72°C for 1 min;with a final extension at 72°C for 10 min,using a thermal cycler(Biometra,Goettingen, Germany).The amplified PCR products were separated on a 6% polyacrylamide sequencing gel after mixing with 4μL of 6× loading buffer(0.25%bromophenolblue,0.25%xylene cyanol FF, and 40%sucrose).

    SSR markers that were polymorphic between the contrasting bulks were further tested in the full F2populations derived from the crosses Qizhen 76×Xucai 1,and Xucai 1×Bawan 6. Based on the results of preliminary genetic linkage analysis, markers near the linked markers were further selected from a consensus functionalmap[47]to identify polymorphisms and analyze genetic linkage and were amplified using the method described above for SSR markers.

    Four SCAR markers,ScOPD10-650[22],ScOPE16-1600[23], ScOPO18-1200[23],and ScOPX04-880[25]linked to resistance gene er1,were also used to identify polymorphisms between the two parents and contrasting bulks and to analyze genetic linkage.PCR amplification and reaction were performed as described previously[34].

    2.5.Data analysis and genetic linkage map construction

    Linkage analysis was performed with MAPMAKER/EXP version 3.0b[48].Genetic distances were computed using the Kosambi mapping function[49].Linkage groups were determined usinga logarithm of odds(LOD)threshold of 3.0,with a maximum distance of 50 cM.The segregation of molecular markers in the F2populations was evaluated for goodness-of-fit to Mendelian segregation ratios byχ2test.The genetic linkage map was drawn with molecular markers linked to the resistance gene using MapDraw[50].

    2.6.Sequence analysis of PsMLO1 gene

    Total RNA was extracted from the young leaves of three parents Xucai 1,Bawan 6,and Qizhen 76 using an RNAprep pure Plant kit(Tiangen Biotech,Co.,Ltd.,Beijing,China), according to the manufacturer's instructions.

    First-strand cDNAs were synthesized using a BioRT Two Step RT-PCR Kit(Hangzhou Bioer Technology Co.,Ltd,Hangzhou, China)including an oligo(dT)primer.The PCRs were performed to amplify the full-length cDNAofthe target gene,PsMLO1,using the PsMLO1-specific primer pair PsMLO1F/PsMLO1R(F:5′-AAAAT GGCTGAAGAGGGAGTT-3′;R:5′-TCCACAAATCAAGCTGCTACC-3′)[29].The amplification reaction was performed under the following conditions:5 minat 95°C;35 cycles of30 s at94°C,45 s at 58°C(annealing temperature),and 70 s at 72°C;with a final extension of 10 min at72°C.Amplicons were purified with a PCR Purification Kit(Qiagen)and ligated into the pEasy_T5 vector (TransGen Biotech,Beijing,China).Sequencing reactions were performed using universal M13F or M13R primers at the Beijing Genomics Institute(BGI).The resulting cDNA sequences of the PsMLO1 gene in the parents were compared with the reported cDNA sequence of the PsMLO1 gene in the susceptible pea cv. Sprinter,having the wild-type PsMLO1 sequence(NCBI,accession number FJ463618.1)[28],using ClustalX2[51].

    2.7.Amplification of the molecular marker er1-2/MGB

    PCR amplification of marker er1-2/MGB linked to er1-2 in the parents was performed in a total volume of 20μL reaction, containing 50 ng of genomic DNA,2.5μL of 10×PCR reaction buffer(containing 20 mmol L-1MgCl2),0.8 mmol L-1of each dNTP,2.5 U oflong-range TaKaRa LA Taq(TaKaRa Bio Inc.,Otsu, Japan),and 0.2 mmolL-1of primer mixture.PCRreactions were performed in a thermalcycler(Biometra,Goettingen,Germany), at 94°C for 5 min;followed by 35 cycles of 94°C for 45 s,54°C for 3 min,and 72°Cfor 1 min;with a finalextension at 72°Cfor 15 min.The amplified PCR products were separated on a 6% polyacrylamide gel.

    3.Results

    3.1.Phenotypic evaluation of resistance to E.pisi

    3.1.1.Parentalphenotypic evaluation

    The two susceptible parents Qizhen 76 and Bawan 6 showed severe infection by E.pisi isolate EPBJ10 days after inoculation under greenhouse conditions.All plants of Qizhen 76 and Bawan 6 were covered by masses of conidia and mycelium of E.pisi,showing a disease severity of 4 and confirming their susceptibility(Fig.1;Table 1).By contrast,there were no symptoms on the resistant parent Xucai1,which appeared to be completely resistant to E.pisi isolate EPBJ(Fig.1;Table 1). These results are consistent with previous observations [31,33].

    3.1.2.Phenotypic evaluations of F1,F2,and F2:3populations

    The segregation pattern of pea resistance to powdery mildew in the F1,F2,and F2:3populations derived from the three crosses of Qizhen 76×Xucai 1,Bawan 6×Xucai 1,and Xucai 1×Bawan 6 are presented in Table 1.

    The cross of Qizhen 76×Xucai 1 generated four F1plants that were susceptible to E.pisi(Table 1).One of four F1plants produced 91 F2plants,of which 17 were resistant and 74 susceptible.Thus,the segregation of susceptibility and resistance in F2population fitted a ratio of 3:1(χ2=1.94;P=0.16), indicating monogenic recessive inheritance.Moreover,the segregation ratio in the F2:3population gave a good fit to a 1:2:1 ratio(χ2=2.41;P=0.30)(Table 1),confirming that the resistance was controlled by a single recessive gene in Xucai1.

    The cross between the susceptible parent Bawan 6 and the resistant parent Xucai 1 also generated four F1plants susceptible to E.pisi.The F2population comprised 58 plantsthat segregated into 47 susceptible and 11 resistant,with a good fit to a 3:1 ratio(χ2=1.13;P=0.28),as expected for single recessive gene inheritance(Table 1).

    Fig.1–Phenotypic evaluation of the resistant pea cultivar Xucai 1(A)and susceptible cultivars Bawan 6(B)and Qizhen 76 (C)inoculated with E.pisi isolate EPBJ.

    Table 1–The segregation patterns of pea resistance to powdery mildew in the F1,F2,and F2:3populations derived from the crosses of Qizhen 76×Xucai 1,Bawan 6×Xucai 1,and Xucai 1×Bawan 6.

    The cross between the resistant parent Xucai 1 and the susceptible parent Bawan 6 generated eightsusceptible F1plants. In the F2population of 161 plants,35 were resistant and 126 susceptible,fitting a 3:1 ratio(χ2=0.91;P=0.34).Moreover,in the F2:3population(161 families),35 were homozygous resistant,48 were homozygous susceptible,and 78 segregated in response to isolate E.pisi EPBJ(Table 1),fitting the genetic model ratio 1:2:1 (χ2=2.26;P=0.32)for a single recessive gene.

    These results indicated that the observed segregation patterns fitted 3(susceptible):1(resistant)and 1(homozygous susceptible):2(segregating):1(homozygous resistant)ratios in all F2and F2:3populations derived from the three crosses, respectively,suggesting that the resistance to E.pisi in Xucai1 was controlled by a single recessive gene.

    3.2.Molecular marker analysis

    Among the 148 SSR markers,59 were polymorphic between the parents Qizhen 76 and Xucai1 and 74 between Xucai1 and Bawan 6.These polymorphic SSR markers were further used to screen for polymorphism between the resistant and susceptible DNA pools.Only SSR marker AD60 in pea LG VI was polymorphic between the two contrasting DNA bulks differing in their response to E.pisi[46]and was thus likely to be linked to the resistance gene.SSR marker AD60 was verified as a marker linked to the er1 locus by Ek et al.[19]. This marker was accordingly genotyped in both F2populations.The marker was codominant and its segregation in both F2populations fitted the expected Mendelian ratio of 1:2:1. Further linkage analysis confirmed that the SSR marker AD60 was linked to the resistance gene(Fig.2-A),showing that the gene is located on LG VI,which harbors the er1 locus[19].

    In addition,according to this primary mapping result, molecular markers on LG VI adjacent to SSR AD60 were selected from a genetic map[47]and analyzed in the contrasting DNA bulks.The codominant marker c5DNAmet (cytosine-5 DNA-methyltransferase)was polymorphic between the bulks and thus likely to be linked to the resistance gene.Its segregation in both F2populations fitted the expected Mendelian ratio of 1:2:1.Further linkage analysis confirmed that c5DNAmet was linked to the resistance gene(Fig.2-B).

    Fig.2–The segregation patterns of molecular markers AD60(A)and c5DNAmet(B)linked to the resistance gene in some of the F2progeny derived from the cross of Qizhen 76×Xucai1.M:100 bp marker ladder;PR:resistant parent Xucai 1;PS:susceptible parent Qizhen 76;R:resistant families;H:families segregating for resistance and susceptibility;S:susceptible families.Arrows indicate recombinants between the markers and the resistance gene.

    Table 2–Amplification of four SCAR markers known to be linked to resistance gene er1 in resistant parent Xucai 1 and susceptible parents Bawan 6 and Qizhen 76(“–”indicates no amplification).

    The four known SCAR markers Sc OPD10-650[22], ScOPE16-1600[23],Sc OPO18-1200[23],and ScOPX04-880 [25]linked to er1 were used to identify polymorphism between the parents Xucai 1 and Bawan 6 and between Xucai 1 and Qizhen 76.The results of amplification by SCAR markers in the three parents are shown in Table 2.Marker ScOPD10-650 amplified a fragment of expected size(650-bp)in the resistant Xucai 1 and the susceptible Qizhen 76,but there was no amplification in the other susceptible cultivar Bawan 6.Marker ScOPE16-1600 amplified a 1600-bp fragment in the susceptible Qizhen 76,but there was no amplification in either resistant Xucai1 or susceptible Bawan 6.Marker ScOPX04-880 amplified an 880-bp fragment in both resistant Xucai 1 and susceptible Bawan 6,but produced no fragment in susceptible Qizhen 76. Marker ScOPO18-1200 amplified a 1200-bp fragment only in susceptible Bawan 6.Amplifications in the contrasting bulks pooled from F2plants were then performed using the SCAR markers polymorphic in the parents.However,these markers did not form different fragment patterns in the contrasting bulks.

    3.3.Genetic mapping of the resistance gene

    Linkage analysis of segregation data allowed construction of genetic maps for both F2mapping populations derived from the crosses of Xucai 1×Bawan 6 and Qizhen 76×Xucai 1,respectively.Maps comprising two SSRmarkers and the resistance gene for both F2populations were constructed with MapDraw[50]. Linked markers AD60 and c5DNAmet were located on the same side of the resistance gene.For the F2population derived from Qizhen 76×Xucai1,markers c5DNAmet and AD6 were linked to the resistance gene at genetic distances of 8.1 and 8.7 cM, respectively(Fig.3-A).For the F2population derived from Xucai 1×Bawan 6,the resistance gene was linked to markers AD60 and c5DNAmet at genetic distances of 9.9 and 15.4 cM,respectively (Fig.3-B).These results indicated that the powdery mildew resistance gene in Xucai1 is an allele of er1,in agreement with a previous study[19].

    Fig.3–Genetic linkage maps of F2populations derived from Qizhen 76×Xucai1(A)and Xucai1×Bawan 6(B)showing the position ofthe gene conferring resistance to powdery mildew and the molecular markers AD60 and c5DNAmet.Map distances and locus orders were determined with MAPMAKER 3.0(Lander et al.1993).Estimated genetic distances between loci are shown on the left in centiMorgans(cM).

    3.4.PsMLO1 cDNA sequence analysis

    We determined the nucleotide sequences of homologous PsMLO1 cDNAs of the contrasting parents to confirm which allele of er1 is present in Xucai 1.PsMLO1 full-length coding sequences of the three parents were obtained by cloning and sequencing of PCR products amplified by the specific primers PsMLO1F and PsMLO1R[29].On the basis of 10 clones,the homologous cDNA sequences of the PsMLO1 gene from Bawan 6 and Qizhen 76 were identical to that of PsMLO1 from the wild-type pea cv.Sprinter(susceptible to E.pisi;NCBI accession number FJ463618).Thus,Bawan 6 and Qizhen 76 carry the susceptibility gene Er1(Fig.4).In contrast,the cDNA sequences of PsMLO1 from Xucai 1 were markedly different from that ofwild-type pea cv.Sprinter and the two susceptible parents.The cDNAs from Xucai 1 comprised three distinct types of mis-spliced transcripts that were not present in the wild-type pea or the two susceptible parents.On the basis of 10 clones from Xucai 1,we obtained three PsMLO1 transcripts that were characterized by a 129-bp deletion and 155-and 220-bp insertions,respectively,in the PsMLO1 gene of wild-type Sprinter (Fig.4).Compared to the PsMLO1 transcript of the susceptible cv. Sprinter,Bawan 6,and Qizhen 76,cDNA1 from Xucai1 showed a 129-bp deletion between positions 1171 and 1299 of the PsMLO1 coding sequence,and the cDNA2 and cDNA3 from Xucai 1 had larger insertions of 155-and 220-bp at position 1263 in PsMLO1 coding sequence,respectively(Fig.4).These distinct transcripts were attributed to alternative splicing during transcription.The cDNA1 sequence also showeda 5-bp deletion frompositions 770-774 of the PsMLO1 coding sequence(Fig.5),which was attributed to an aberrant splicing event.The transcripts ofthe PsMLO1 gene in Xucai1 are consistent with those in pea cv.Stratagem(JI 2302) and Franklin which harbor the er1-2 alleles of resistance gene er1, indicating that Xucai 1 carries the er1-2 allele for resistance to powdery mildew[28,29].

    4.Discussion

    Fig.4–Comparisons of cDNA sequence of the PsMLO1 gene from the wild-type pea cultivar Sprinter(Wild_type_PsMLO1)with those identified in three parents,including two susceptible cultivars,Qizhen 76(QZ76_PsMLO1_cDNA)and Bawan 6 (BW6_PsMLO1_cDNA),and a resistant cultivar,Xucai 1(XC1_PsMLO1_cDNA1,2 and 3).Mutation sites with the deletion or insertion fragments in PsMLO1 coding sequence are boxed in black.

    In this study,the pea cv.Xucai 1 was confirmed to have high resistance(immunity)to powdery mildew,and no colonies developed on the whole plants after artificialinoculation with E.pisi(Fig.1,Table 1).The genetic basis of resistance to E.pisi in Xucai 1 was investigated using three F2and two F2:3populations derived from three crosses.Phenotypic observations of the populations showed that resistance in Xucai 1 is conferred by a single recessive gene.

    To date,only two recessive genes(er1 and er2)and one dominant gene(Er3)have been identified as conferring resistance to E.pisi[10-12].A majority ofpea germplasm resistant to E.pisi carries the recessive gene er1,which provides durable broad-spectrum resistance to powdery mildew disease caused by E.pisi[13,52].

    Fig.5–Alignment showing a 5-bp deletion from 770 to 774 bp,compared to the wild type(Wild_type_PsMLO1),in one of the transcripts derived from Xucai1(XC1_PsMLO1_cDNA1).Transcripts“QZ76_PsMLO1_cDNA”and“BW6_PsMLO1_cDNA”derived from the two susceptible cultivars Qizhen 76 and Bawan 6,respectively.

    To aid marker-assisted selection(MAS)in pea breeding programs,studies have been performed to identify different types of molecular markers,including RFLP(restriction fragment length polymorphism),RAPD,SCAR,and SSR, linked to the er1 locus[17-25,53].In these studies,various markers linked to the resistance gene er1 were reported by different researchers;however,the markers were located at different genetic distances in different mapping populations. First,Sarala[17]and Timmerman et al.[22]placed the er1 gene on pea LG VI.Dirlewanger et al.[18]found the er1 gene 9.8 cM distant from a RFLP marker,p236,whereas Timmerman et al. [22]placed the RAPD marker OPD10650at a distance of 2.1 cM from er1.This linked RAPD marker(OPD10650)was converted to a SCARmarker,ScOPD10650,which was mapped at a distance of 3.4 cM from er1[20].However,Tiwari et al.[23]found that ScOPD10650was not useful to identify Canadian pea resistant germplasm and related F3progeny derived from“Highlight (er1)×Radley.”Tiwari et al.[23]obtained three additional linked RAPD markers:OPO181200,OPE161600,and OPL61900.The first two ofthese were converted into SCARmarkers(ScOPE161200and ScOPL61600)and used successfully to identify resistant pea germplasm.Janila and Sharma[20]found that of the three SCAR markers,only ScOPD10650was linked to er1,at a genetic distance of 3.4 cM.Later,Ek et al.[19]discovered five SSR markers(PSMPSAD51,PSMPSA5,PSMPSAD60,PSMPSAA374e, and PSMPSAA369)linked to er1,and PSMPSAD60(named AD60 in this study)was identified as the most closely linked marker,at 10.4 cM distance from the er1 locus.Pereira et al.[21]identified two mutated genes,er1mut1 and er1mut2,that induced powdery mildew resistant mutants in pea,via ethylnitrosourea mutagenesis.The er1mut1 and er1mut2 loci were mapped to a pea LG containing 16 DNA markers including six linked markers.Pereira et al.[21]also attempted to remap all previously reported SCARs and five SSR markers linked to er1.Only the reported SCAR marker ScOPL61600and SSR marker PSMPSA5 were linked to the er1mut2 locus,at a longer distance.Recently, Santo et al.[54]developed specific molecular markers SCAR and a sequence-tagged site(STS)for the mutation er1mut1. Srivastava et al.[25]also developed another coupling-phase SCARmarker,ScOPX04880,linked to er1.In their study,a previous repulsion-phase marker,ScOPD10650,and a coupling-phase marker,ScOPX04880,were mapped flanking er1 at 2.2 and 0.6 cM distance,respectively.Taken together,these studies identified several markers linked to the er1 locus,with various degrees of utility,in different plant materials.In our study,in addition to the previous SCAR marker AD60,gene marker c5DNAmet was newly identified as linked to resistance gene er1.This novel marker can be used in MAS for powdery mildew resistance breeding in pea.

    In this study,the resistance gene in Xucai1 was mapped at distances of 8.7 and 9.9 cM from the linked marker AD60 in the respective F2populations(Fig.3-A and 3-B).Our result is comparable with that reported by Ek et al.[19].The parents were then tested with four previous SCAR markers linked to the er1 locus,ScOPD10650,ScOPO181200,ScOPE161600,and ScOPX04880using the methods described by Wang et al.[34]. However,the four markers were useless among the three tested parents because they were unable to distinguish the resistant parent from the susceptible ones and showed the same single band in the contrasting bulks(Table 2).This result was consistent with those of previous studies[20,21,23].We then selected molecular markers distributed in the vicinity of the linked marker AD60 on a pea linkage map[47]to identify additional markers linked to the resistance gene in our populations,and found that gene marker c5DNAmet was linked to the resistance gene in Xucai 1.Thus,we further confirmed that powdery mildew resistance gene in Xucai 1 resides on LG VI of the pea genetic map.The marker c5DNAmet mapped at 8.1 cM and 15.4 cM distances from the resistance gene in the two F2populations derived from Xucai 1×Bawan 6 and Qizhen 76×Xucai 1,respectively(Fig.3). Markers AD60 and c5DNAmet were mapped to the same side of the resistance gene on LG VI in the two F2populations (Fig.3).The mapping of the resistance gene in LG VIlinked to AD60 and c5DNAmet markers suggested that the resistance gene is an er1 allele.

    Recently,er1 resistance was demonstrated to be caused by a loss-of-function mutations of a powdery mildew susceptibility gene,PsMLO1,belonging to the MLO gene family[28,29].To date,five er1 alleles,each corresponding to a different PsMLO1 mutation,have been reported in pea accessions resistant to E.pisi [30].Humphry etal.[28]identified four er1 resistance alleles,er1-1 (JI1559),er1-2(JI2302),er1-3(JI210),and er1-4(JI1951),in naturally occurring PsMLO1 loss-of function mutations in resistant pea accessions.Pavan et al.[29]identified another er1 allele,er1-5,in pea line ROI3/02 following a mutagenesis program.Four of the five er1 alleles(excluding er1-2)are caused by point mutations in the PsMLO1 gene.Allele er1-2 comprises a large DNA insertion of unknown size and identities in the PsMLO1 gene,leading to aberrant PsMLO1 transcription[28,29].Allele er1-2 has been identified in a few pea powdery mildew resistant cultivars/lines, including Stratagem(JI2302),Franklin,Dorian,and Nadir[28,29]. Stratagem(JI 2302)is the earliest identified pea cultivar resistantto E.pisi and has been widely used in pea breeding,owing to its stable resistance[52].

    In the present study,the mutation mode that occurred in Xucai 1 is consistent with that in the pea cv.Stratagem(JI 2302)carrying the er1 resistance allele,er1-2[28,29].Thus, the resistance gene in cv.Xucai 1 is an er1 allele,er1-2.The inserted fragments of 155-and 220-bp were compared against the NCBI database by BLASTN.The result revealed that,as reported for er1-2 allele by Humphry et al.[28],the 155-and 220-bp insertions were very similar(96%and 95%,respectively)to a sequence repeated five times in the genomic BAC clone(GenBank accession number,CU655882).Interestingly,the 220-bp insertion is also highly similar(~87%identity)to part of the sequence of the pea Ogre retrotransposon(GenBank accession number AY299395).The Ogre is a giant retroelement of more than 22 kb that makes up at least 5%of the pea genome [54,55].

    Recently,Pavan et al.[30]developed different types of functionalmarkers corresponding to five er1 alleles,including an STS dominant marker er1-2/MGB targeting the er1-2 allele, according to the sequence difference.er1-2/MGB was used to test our parents in the present study,but showed no amplification in either resistant or susceptible parents.To date,er1-2/MGB has not been tested in other resistance cultivars or in derived populations carrying er1-2,except Franklin[30].Our finding may have been due to the large difference in genetic background between Xucai 1 and Franklin.Such difference can be dissected with obtaintion of complicated genome sequence of pea carrying resistance gene er1-2 after developing a codominant marker for er1-2.Four foreign pea cultivars or lines,including Stratagem (JI2302),Franklin,Dorian,and Nadir,have been shown to harbor the resistance gene er1-2[28,29].In our study,we identified another recessive resistance source harboring er1-2:the Chinese pea cultivar Xucai1.

    The desirable agronomic traits of Xucai 1,including high yield and tolerance to lodging and drought[31],would be an advantage for the use of this cultivar in breeding programs.Our study revealed that resistance to E.pisi in the Xucai 1 is conferred by the gene er1-2.This result provides valuable genetic information about Xucai 1 that will improve its use in pea breeding for resistance to E.pisi.

    5.Conclusions

    Powdery mildew resistance in Xucai 1 was governed by a recessive gene that was mapped in pea LG VI,close to markers AD60 and c5DNAmet.These results suggested that the resistance gene is an allele of er1.The analysis of PsMLO1 coding sequence in Xucai 1 revealed that the resistance gene was generated by a large insertion or deletion in the PsMLO1 gene. Thus,the resistance gene in Xucai 1 was identified as an er1 allele,specifically the er1-2 allele.

    Acknowledgements

    This study was supported by the Modern Agro-industry Technology Research System(CARS-09),the Crop Germplasm Conservation and Utilization Program(2014NWB030-14)fromthe Ministry of Agriculture of China,and the Scientific Innovation Program of Chinese Academy of Agricultural Sciences.We thank Prof.Fengbao Wang(Hebei Normal University of Science& Technology,Qinghuangdao,Hebeiprovince,China),Dr.Xin Chen (Jiangsu Academy of Agricultural Sciences,Nanjing,Jiangsu province,China),and Mr.Dongxu Xu(Zhangjiakou Academy of Agricultural Sciences,Zhangjiakou,Hebei province,China)for providing the pea cultivars.

    R E F E R E N C E S

    [1]L.J.Lang,Z.J.Zheng,Production of faba bean and pea in China, in:R.J.Summerfield(Ed.),World Crops:Cool Season Food Legumes,Current Plant Science and Biotechnology in Agriculture,Kluwer Academic Publisher,Dordrecht,Netherlands 1988,pp.1135-1152.

    [2]F.A.O.Statistical Database,Food and Agriculture Organization(FAO)of the United Nations,Rome,2014 (http://faostat3.fao.org/).

    [3]E.T.Gritton,R.D.Ebert,Interaction of planting date and powdery mildew on pea plant performance,Am.Soc. Horticult.Sci.100(1975)137-142.

    [4]P.H.Smith,E.M.Foster,L.A.Boyd,J.K.M.Brown,The early development of Erysiphe pisi on Pisum sativum L,Plant Pathol. 45(1996)302-309.

    [5]R.L.Munjal,V.V.Chenulu,T.S.Hora,Assessment of losses due to powdery mildew(Erisiphe polygoni)on pea,Indian Phytopathol.19(1963)260-267.

    [6]T.D.Warkentin,K.Y.Rashid,A.G.Xue,Fungicidal control of powdery mildew in field pea,Can.J.Plant Sci.76(1996) 933-935.

    [7]S.Fondevilla,J.I.Cubero,D Rubiales,Confirmation that the Er3 gene,conferring resistance to Erysiphe pisi in pea,is a different gene from er1 and er2 genes,Plant Breed.130(2011) 281-282.

    [8]A.Ghafoor,K.McPhee,Marker assisted selection(MAS)for developing powdery mildew resistant pea cultivars, Euphytica 186(2012)593-607.

    [9]S.Fondevilla,D.Rubiales,Powdery mildew control in pea:a review,Agron.Sustain.Dev.32(2012)401-409.

    [10]S.C.Harland,Inheritance of immunity to mildew in Peruvian forms of Pisum sativum,Heredity 2(1948)263-269.

    [11]R.J.Heringa,A.van Norel,M.F.Tazelaar,Resistance to powdery mildew(Erysiphe polygoni D.C.)in peas(Pisum sativum L.),Euphytica 18(1969)163-169.

    [12]S.Fondevilla,A.M.Torres,M.T.Moreno,D.Rubiales, Identification of a new gene for resistance to powdery mildew in Pisum fulvum,a wild relative of pea,Breed.Sci.57 (2007)181-184.

    [13]K.R.Tiwari,G.A.Penner,T.D.Warkentin,Inheritance of powdery mildew resistance in pea,Can.J.Plant Sci.77 (1997)307-310.

    [14]S.Fondevilla,T.L.W.Carver,M.T.Moreno,D.Rubiales, Macroscopic and histological characterisation of genes er1 and er2 for powdery mildew resistance in pea,Eur.J.Plant Pathol.115(2006)309-321.

    [15]R.Iglesias-Garc?a,D.Rubiales,S.Fondevilla,Penetration resistance to Erysiphe pisi in pea mediated by er1 gene is associated with protein cross-linking but not with callose apposition or hypersensitive response,Euphytica 201(2015) 381-387.

    [16]S.Fondevilla,T.L.W.Carver,M.T.Moreno,D.Rubiales, Identification and characterisation of sources of resistance to Erysiphe pisi Syd.in Pisum spp.Plant Breed.126(2007) 113-119.

    [17]K.Sarala,Linkage studies in pea(Pisum sativum L.)with reference to er gene for powdery mildew resistance and other genesPh.D.thesis Indian Agricultural Research Institute,New Delhi,India,1993.

    [18]E.Dirlewanger,P.G.Isaac,S.Ranade,M.Belajouza,R.Cousin, D.Vienne,Restriction fragment length polymorphism analysis of loci associated with disease resistance genes and developmental traits in Pisum sativum L,Theor.Appl.Genet. 88(1994)17-27.

    [19]M.Ek,M.Eklund,R.von Post,C.Dayteg,T.Henriksson,P. Weibull,A.Ceptilis,P.Issac,S.Tuvesson,Microsatellite markers for powdery mildew resistance in pea(Pisum sativum L.),Hereditas 142(2005)86-91.

    [20]P.Janila,B.Sharma,RAPD and SCAR markers for powdery mildew resistance gene er in pea,Plant Breed.123(2004) 271-274.

    [21]G.Pereira,C.Marques,R.Ribeiro,S.Formiga,M.Damaso,T. Sousa,M.Farinhó,J.M.Leit?o,Identification of DNA markers linked to an induced mutated gene conferring resistance to powdery mildew in pea(Pisum sativum L.),Euphytica 171 (2010)327-335.

    [22]G.M.Timmerman,T.J.Frew,N.F.Weeden,Linkage analysis of er1,a recessive Pisum sativum gene for resistance to powdery mildew fungus(Erysiphe pisi D.C.),Theor.Appl.Genet.88 (1994)1050-1055.

    [23]K.R.Tiwari,G.A.Penner,T.D.Warkentin,Identification of coupling and repulsion phase RAPD markers for powdery mildew resistance gene er1 in pea,Genome 41(1998)440-444.

    [24]M.Tonguc,N.F.Weeden,Identification and mapping of molecular markers linked to er1 gene in pea,J.Plant Mol.Biol. Biotechnol.1(2010)1-5.

    [25]R.K.Srivastava,S.K.Mishra,K.Singh,T.Mohapatra, Development of a coupling-phase SCAR marker linked to the powdery mildew resistance gene er1 in pea(Pisum sativum L.), Euphytica 186(2012)855-866.

    [26]V.Katoch,S.Sharma,S.Pathania,D.K.Banayal,S.K.Sharma, R.Rathour,Molecular mapping of pea powdery mildew resistance gene er2 to pea linkage group III,Mol.Breed.25 (2010)229-237.

    [27]S.Fondevilla,D.Rubiales,M.T.Moreno,A.M.Torres, Identification and validation of RAPD and SCAR markers linked to the gene Er3 conferring resistance to Erysiphe pisi D.C.in pea,Mol.Breed.22(2008)193-200.

    [28]M.Humphry,A.Reinst?dler,S.Ivanov,T.Bisseling,R. Panstruga,Durable broad-spectrum powdery mildew resistance in pea er1 plants is conferred by natural loss-of-function mutations in PsMLO1,Mol.Plant Pathol.12 (2011)866-878.

    [29]S.Pavan,A.Schiavulli,M.Appiano,A.R.Marcotrigiano,F. Cillo,R.G.F.Visser,Y.Bai,C.Lotti,L.Ricciardi,Pea powdery mildew er1 resistance is associated to loss-of-function mutations at a MLO homologous locus,Theor.Appl.Genet. 123(2011)1425-1431.

    [30]S.Pavan,A.Schiavulli,M.Appiano,C.Miacola,R.G.F.Visser, Y.Bai,C.Lotti,L.Ricciardi,Identification of a complete set of functional markers for the selection of er1 powdery mildew resistance in Pisum sativum L,Mol.Breed.31 (2013)247-253.

    [31]X.M.Wang,Z.D.Zhu,C.X.Duan,X.X.Zong,Identification and Control Technology of Disease and Pest on Faba Bean and Pea,Chinese Agricultural Science and Technology Press, Beijing,2007.

    [32]H.X.Peng,G.Yao,R.L.Jia,H.Y.Liang,Identification of pea germplasm resistance to powdery mildew,J.Southwest Agric.Univ.13(1991)384-386(in Chinese with English abstract).

    [33]L.Zeng,M.Q.Li,X.M.Yang,Identification of resistance of peas resources to powdery mildew,Grass&Turf 32(2012) 35-38(in Chinese with English abstract).

    [34]Z.Y.Wang,S.Y.Bao,C.X.Duan,X.X.Zong,Z.D.Zhu,Screening and molecular identification of resistance to powdery mildew in pea germplasm,Acta Agron.Sin.39(2013) 1030-1038(in Chinese with English abstract).

    [35]H.N.Fu,S.L.Sun,Z.D.Zhu,C.X.Duan,X.M.Yang,Phenotypic and genotypic identification of powdery mildew resistance in pea cultivars or lines from Canada,J.Plant Genet.Resour.15 (2014)1028-1033(in Chinese with English abstract).

    [36]M.Ond?ej,R.Dostálová,L.Odstr?ilová,Response of Pisum sativum germplasm resistant to Erysiphe pisi to inoculation with Erysiphe baeumleri,a new pathogen of peas,Plant Prot. Sci.41(2005)95-103.

    [37]R.N.Attanayakea,D.A.Glaweab,K.E.McPheec,F.M.Dugand, W.Chend,Erysiphe trifolii-a newly recognized powdery mildew pathogen of pea,Plant Pathol.59(2010)712-720.

    [38]S.Fondevilla,C.Chattopadhyay,N.Khare,D.Rubiales, Erysiphe trifolii is able to overcome er1 and Er3,but not er2, resistance genes in pea,Eur.J.Plant Pathol.136(2013) 557-563.

    [39]P.Baiswar,S.V.Ngachan,V.K.Verma,R.Kumar,A.K.Jha,S. Chandra,Molecular evidence of Erysiphe pisi on pea and E. trifoliorum on white clover in northeast India,Aust.Plant Dis. Notes 10(2015)12.

    [40]L.F.Dong,F.B.Wang,J.F.Fu,Development of a new semi-leafless sweet pea Xucai No.1,J.Changjiang Veget. (2008)45-46(in Chinese with English abstract).

    [41]A.Vaid,P.D.Tyagi,Genetics of powdery mildew resistance in pea,Euphytica 96(1997)203-206.

    [42]J.C.Rana,D.K.Banyal,K.D.Sharma,K.Sharma Manish,S.K. Gupta,S.K.Yadav,Screening ofpea germplasm for resistance to powdery mildew,Euphytica 189(2013)271-282.

    [43]D.Rubiales,J.K.M.Brown,A.Martín,Hordeum chilense resistance to powdery mildew and its potentialuse in cereal breeding,Euphytica 67(1993)215-220.

    [44]M.Shure,S.Wessler,N.Fedoroff,Molecular-identification and isolation of the waxy locus in maize,Cell 35(1983)225-233.

    [45]R.W.Michelmore,I.Paran,R.V.Kesseli,Identification of markers linked to disease-resistance genes by bulked segregant analysis:a rapid method to detect markers in specific genomic regions by using segregating populations,Proc.Natl.Acad.Sci.U.S.A.88(1991) 9828-9832.

    [46]K.Loridon,K.McPhee,J.Morin,P.Dubreuil,M.L.Pilet-Nayel, G.Aubert,C.Rameau,A.Baranger,C.Coyne,I. Lejeune-Henaut,J.Burstin,Microsatellite marker polymorphism and mapping in pea(Pisum sativum L.),Theor. Appl.Genet.111(2005)1022-1031.

    [47]A.Bordat,V.Savois,M.Nicolas,J.Salse,A.Chauveau,M. Bourgeois,J.Potier,H.Houtin,C.Rond,F.Murat,P.Marget,G. Aubert,J.Burstin,Translationalgenomics in legumes allowed placing in silico 5460 unigenes on the pea functional map and identified candidate genes in Pisum sativum L,Genes Genom.Genet.1(2011)93-103.

    [48]E.S.Lander,M.J.Daly,S.E.Lincoln,Constructing Genetic Linkage Maps with MAPMAKER/EXP Version 3.0:A TutorialAnd Reference Manual,Third EditionAWhitehead Institute for Biomedical Research Technical Report,Cambridge,MA,USA, 1993.

    [49]D.D.Kosambi,The estimation of map distances from recombination values,Ann.Eugen.12(1944)172-175.

    [50]R.H.Liu,J.L.Meng,MapDraw:a Microsoft excel macro for drawing genetic linkage maps based on given genetic linkage data,Hereditas 25(2003)317-321(in Chinese with English abstract).

    [51]M.A.Larkin,G.Blackshields,N.P.Brown,R.Chenna,P.A. McGettigan,H.McWilliam,F.Valentin,I.M.Wallace,A. Wilm,R.Lopez,J.D.Thompson,T.J.Gibson,D.G.Higgins, Clustal Wand Clustal X version 2.0,Bioinformatics 23(2007) 2947-2948.

    [52]R.Cousin,Resistance to powdery mildew in pea,Ann. Amélior.Plantes 15(1965)93-97.

    [53]M.Nisar,A.Ghafoor,Linkage of a RAPD marker with powdery mildew resistance er1 gene in Pisum sativum L,Russ.J.Genet. 47(2011)300-304.

    [54]T.Santo,M.Rashkova,C.Alabaca,J.Leitao,The ENU-induced powdery mildew resistant mutant pea(Pisum sativum L.)lines S(er1mut1)and F(er1mut2)harbour early stop codons in the PsMLO1 gene,Mol.Breed.32(2013)723-727.

    [55]P.Neumann,D.Pozarkova,J.Macas,Highly abundant pea LTR retrotransposon Ogre is constitutively transcribed and partially spliced,Plant Mol.Biol.53(2003)399-410.

    17 March 2015

    in revised form13 July 2015

    .

    E-mail address:zhuzhendong@caas.cn(Z.Zhu).

    Peer review under responsibility of Crop Science Society of China and Institute of Crop Science,CAAS.1Suli Sun and Zhongyi Wang contributed equally to this work.

    成人永久免费在线观看视频| 麻豆一二三区av精品| 国产精品爽爽va在线观看网站| 久久人人精品亚洲av| 国产免费男女视频| 欧美黄色淫秽网站| 日本a在线网址| 国产伦精品一区二区三区视频9 | 国产伦在线观看视频一区| 熟女人妻精品中文字幕| 日本免费a在线| 成人国产一区最新在线观看| 特级一级黄色大片| 99热这里只有精品一区| 亚洲五月天丁香| 久久99热这里只有精品18| 亚洲 国产 在线| 九色国产91popny在线| 免费观看的影片在线观看| 黄色女人牲交| 夜夜看夜夜爽夜夜摸| 亚洲国产中文字幕在线视频| 亚洲黑人精品在线| 国产精品香港三级国产av潘金莲| 精品人妻1区二区| 国产探花极品一区二区| 精品久久久久久成人av| 亚洲精品久久国产高清桃花| 99热只有精品国产| 国产一区二区三区视频了| 少妇的逼水好多| 欧美日韩一级在线毛片| 国产精品女同一区二区软件 | 国产精品久久久久久人妻精品电影| 亚洲美女视频黄频| 久久香蕉精品热| 成人一区二区视频在线观看| 97超级碰碰碰精品色视频在线观看| 18禁美女被吸乳视频| 欧美黑人欧美精品刺激| 日本与韩国留学比较| 精品国产亚洲在线| 国产伦在线观看视频一区| 成人18禁在线播放| 黄色片一级片一级黄色片| 天美传媒精品一区二区| 久久婷婷人人爽人人干人人爱| 国产高潮美女av| 国语自产精品视频在线第100页| 亚洲在线自拍视频| 性欧美人与动物交配| 久久午夜亚洲精品久久| 欧美日韩乱码在线| 亚洲熟妇中文字幕五十中出| 成人三级黄色视频| 国产免费一级a男人的天堂| 人人妻人人看人人澡| 日本在线视频免费播放| 在线观看免费午夜福利视频| 亚洲欧美一区二区三区黑人| 国产一区二区亚洲精品在线观看| eeuss影院久久| 欧美日韩黄片免| 亚洲精品一卡2卡三卡4卡5卡| 欧美乱码精品一区二区三区| 黄片大片在线免费观看| 欧美性猛交黑人性爽| 色综合亚洲欧美另类图片| 亚洲美女黄片视频| 亚洲一区二区三区色噜噜| 少妇高潮的动态图| 国产亚洲精品综合一区在线观看| 亚洲不卡免费看| 97超级碰碰碰精品色视频在线观看| 久久久久国内视频| 亚洲av免费在线观看| 国产毛片a区久久久久| 日日摸夜夜添夜夜添小说| 给我免费播放毛片高清在线观看| 亚洲,欧美精品.| 欧美乱色亚洲激情| 99热精品在线国产| 亚洲人成网站在线播| 在线观看免费午夜福利视频| 别揉我奶头~嗯~啊~动态视频| 欧美性感艳星| 国产精品久久久久久精品电影| 久久人人精品亚洲av| 日韩欧美国产一区二区入口| 色哟哟哟哟哟哟| 欧美最新免费一区二区三区 | 国产成人av教育| 1024手机看黄色片| 一个人观看的视频www高清免费观看| 听说在线观看完整版免费高清| 天堂√8在线中文| 丰满的人妻完整版| 国产一区二区在线观看日韩 | 宅男免费午夜| 国产视频一区二区在线看| 欧美日韩国产亚洲二区| 免费在线观看日本一区| 亚洲精品乱码久久久v下载方式 | 精品一区二区三区人妻视频| 免费av不卡在线播放| 激情在线观看视频在线高清| 婷婷精品国产亚洲av| 少妇高潮的动态图| 久久精品亚洲精品国产色婷小说| 99热精品在线国产| 黑人欧美特级aaaaaa片| 欧美日韩亚洲国产一区二区在线观看| 免费看日本二区| 国产v大片淫在线免费观看| 欧美色欧美亚洲另类二区| 国产日本99.免费观看| 女警被强在线播放| 麻豆成人午夜福利视频| tocl精华| 日韩有码中文字幕| 精品熟女少妇八av免费久了| 男人和女人高潮做爰伦理| 性色avwww在线观看| 国产精品久久久人人做人人爽| 白带黄色成豆腐渣| 精品99又大又爽又粗少妇毛片 | 男女之事视频高清在线观看| ponron亚洲| 99久久精品一区二区三区| 好男人在线观看高清免费视频| 91麻豆精品激情在线观看国产| 变态另类成人亚洲欧美熟女| aaaaa片日本免费| 国产一区二区亚洲精品在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产 一区 欧美 日韩| 久久久久国内视频| 宅男免费午夜| 中文亚洲av片在线观看爽| 国产成人啪精品午夜网站| 亚洲狠狠婷婷综合久久图片| 亚洲国产高清在线一区二区三| 婷婷亚洲欧美| 国产亚洲精品一区二区www| 波多野结衣巨乳人妻| 国产欧美日韩精品一区二区| 亚洲在线观看片| 级片在线观看| 精品久久久久久久久久久久久| 色视频www国产| 毛片女人毛片| 日韩高清综合在线| 国产精品 国内视频| 一个人观看的视频www高清免费观看| 国产伦人伦偷精品视频| 欧美性感艳星| 国产欧美日韩精品亚洲av| 日本精品一区二区三区蜜桃| 国产精品,欧美在线| www日本黄色视频网| 精品久久久久久,| 2021天堂中文幕一二区在线观| 老熟妇仑乱视频hdxx| 每晚都被弄得嗷嗷叫到高潮| 人妻夜夜爽99麻豆av| 成人国产一区最新在线观看| 一进一出抽搐gif免费好疼| 亚洲精品一卡2卡三卡4卡5卡| 白带黄色成豆腐渣| 内射极品少妇av片p| 免费观看精品视频网站| 国产精品一区二区免费欧美| 九色国产91popny在线| 级片在线观看| 日本在线视频免费播放| 国产中年淑女户外野战色| x7x7x7水蜜桃| 国产中年淑女户外野战色| 桃红色精品国产亚洲av| 日韩欧美精品免费久久 | 久久精品夜夜夜夜夜久久蜜豆| 国产午夜精品久久久久久一区二区三区 | 日韩欧美 国产精品| 国产成人av激情在线播放| 欧美又色又爽又黄视频| 久久久精品大字幕| 男人的好看免费观看在线视频| 男女下面进入的视频免费午夜| 九九在线视频观看精品| 法律面前人人平等表现在哪些方面| 最近最新中文字幕大全免费视频| 俄罗斯特黄特色一大片| 琪琪午夜伦伦电影理论片6080| 国产精品精品国产色婷婷| 超碰av人人做人人爽久久 | 91久久精品国产一区二区成人 | 午夜激情欧美在线| 国产中年淑女户外野战色| x7x7x7水蜜桃| 女人被狂操c到高潮| 性欧美人与动物交配| 少妇的丰满在线观看| 欧美另类亚洲清纯唯美| 色av中文字幕| 怎么达到女性高潮| 国产欧美日韩精品一区二区| 中文字幕av在线有码专区| 国产免费男女视频| 婷婷精品国产亚洲av| 欧美三级亚洲精品| 午夜福利在线观看免费完整高清在 | av福利片在线观看| 亚洲专区国产一区二区| 夜夜看夜夜爽夜夜摸| 国产一区二区三区视频了| 悠悠久久av| 午夜免费成人在线视频| avwww免费| 有码 亚洲区| 亚洲国产精品成人综合色| 亚洲五月婷婷丁香| 99热6这里只有精品| 亚洲久久久久久中文字幕| 九色国产91popny在线| 波多野结衣巨乳人妻| 欧美日韩瑟瑟在线播放| 精品国产超薄肉色丝袜足j| 97碰自拍视频| 日韩欧美一区二区三区在线观看| www.999成人在线观看| 免费看美女性在线毛片视频| 18美女黄网站色大片免费观看| 亚洲国产色片| 欧洲精品卡2卡3卡4卡5卡区| 精品久久久久久久久久免费视频| 中文资源天堂在线| 夜夜看夜夜爽夜夜摸| 婷婷精品国产亚洲av在线| 18禁国产床啪视频网站| 国产精品精品国产色婷婷| 欧美成人一区二区免费高清观看| 国产成人av教育| 亚洲天堂国产精品一区在线| 久久精品国产综合久久久| 九九在线视频观看精品| 91av网一区二区| 国产精品三级大全| 狠狠狠狠99中文字幕| 好男人在线观看高清免费视频| 99久国产av精品| 啦啦啦韩国在线观看视频| 亚洲国产精品合色在线| 国产乱人视频| www.色视频.com| 精品免费久久久久久久清纯| 别揉我奶头~嗯~啊~动态视频| 熟女少妇亚洲综合色aaa.| 人人妻人人看人人澡| 国产精品久久久久久久电影 | 夜夜爽天天搞| 久99久视频精品免费| 国产精品久久电影中文字幕| 国产探花在线观看一区二区| 一级黄色大片毛片| 丁香欧美五月| 九色成人免费人妻av| 老鸭窝网址在线观看| 黄色片一级片一级黄色片| 黄色丝袜av网址大全| 91九色精品人成在线观看| 国产午夜福利久久久久久| 国产精品香港三级国产av潘金莲| xxx96com| 国产高潮美女av| 亚洲精品国产精品久久久不卡| 国产美女午夜福利| 精品99又大又爽又粗少妇毛片 | 午夜精品久久久久久毛片777| 一本精品99久久精品77| 国产精品1区2区在线观看.| 丰满的人妻完整版| 国产高清有码在线观看视频| av片东京热男人的天堂| av专区在线播放| 悠悠久久av| 久久久久久国产a免费观看| tocl精华| 亚洲午夜理论影院| 一区福利在线观看| 一级作爱视频免费观看| 一区二区三区激情视频| 久久精品91无色码中文字幕| 国语自产精品视频在线第100页| 免费大片18禁| 午夜免费男女啪啪视频观看 | 香蕉久久夜色| 综合色av麻豆| 丁香六月欧美| 少妇的逼好多水| 制服人妻中文乱码| 高清日韩中文字幕在线| av黄色大香蕉| 亚洲激情在线av| 国产精品亚洲av一区麻豆| 国产精品综合久久久久久久免费| 青草久久国产| 国产一级毛片七仙女欲春2| 免费搜索国产男女视频| 757午夜福利合集在线观看| 亚洲性夜色夜夜综合| 脱女人内裤的视频| 88av欧美| 天堂网av新在线| 香蕉av资源在线| 成人亚洲精品av一区二区| 国产午夜精品久久久久久一区二区三区 | 级片在线观看| 亚洲激情在线av| 内地一区二区视频在线| 色噜噜av男人的天堂激情| 久久久国产成人精品二区| 91字幕亚洲| 精品欧美国产一区二区三| 一级黄色大片毛片| 午夜老司机福利剧场| 又黄又爽又免费观看的视频| 久久久久久久精品吃奶| 两性午夜刺激爽爽歪歪视频在线观看| 一夜夜www| 好看av亚洲va欧美ⅴa在| 90打野战视频偷拍视频| 看黄色毛片网站| 女生性感内裤真人,穿戴方法视频| 夜夜躁狠狠躁天天躁| 国产91精品成人一区二区三区| 法律面前人人平等表现在哪些方面| 亚洲精品日韩av片在线观看 | 亚洲专区中文字幕在线| bbb黄色大片| 久久久久久人人人人人| 18禁黄网站禁片午夜丰满| 国产高潮美女av| 亚洲欧美日韩高清专用| 久久久久久人人人人人| 无人区码免费观看不卡| av天堂在线播放| 99热6这里只有精品| 日韩人妻高清精品专区| 亚洲人成伊人成综合网2020| 精品人妻1区二区| 在线免费观看不下载黄p国产 | 国产亚洲精品久久久久久毛片| 久久天躁狠狠躁夜夜2o2o| 久久久久国产精品人妻aⅴ院| 免费看美女性在线毛片视频| 国产探花极品一区二区| 欧美激情在线99| 欧美日韩瑟瑟在线播放| 母亲3免费完整高清在线观看| 在线十欧美十亚洲十日本专区| 制服丝袜大香蕉在线| 国产成+人综合+亚洲专区| 深爱激情五月婷婷| 少妇熟女aⅴ在线视频| 69人妻影院| 日本熟妇午夜| 欧美黑人巨大hd| 国语自产精品视频在线第100页| 91久久精品国产一区二区成人 | 欧美一区二区亚洲| 少妇熟女aⅴ在线视频| а√天堂www在线а√下载| 欧美乱色亚洲激情| 久久6这里有精品| 国产精品乱码一区二三区的特点| av视频在线观看入口| 伊人久久精品亚洲午夜| 久久6这里有精品| 少妇高潮的动态图| 日韩欧美精品免费久久 | 深夜精品福利| 中文字幕久久专区| 综合色av麻豆| 久久精品国产亚洲av涩爱 | 首页视频小说图片口味搜索| 亚洲国产精品999在线| 91麻豆av在线| 久久亚洲精品不卡| 国产综合懂色| 少妇的逼好多水| 真人做人爱边吃奶动态| 国内毛片毛片毛片毛片毛片| 成人国产综合亚洲| a在线观看视频网站| 久久精品91蜜桃| 欧美色视频一区免费| 人妻夜夜爽99麻豆av| 国产真实乱freesex| 国产高清视频在线播放一区| 少妇人妻精品综合一区二区 | 草草在线视频免费看| 亚洲精品乱码久久久v下载方式 | 久久久久久久久中文| 国产精品精品国产色婷婷| 内射极品少妇av片p| 岛国在线免费视频观看| 波野结衣二区三区在线 | 欧美不卡视频在线免费观看| 日韩欧美精品免费久久 | 国产国拍精品亚洲av在线观看 | 色吧在线观看| 中文字幕久久专区| 天天添夜夜摸| 激情在线观看视频在线高清| 午夜免费激情av| 国产午夜福利久久久久久| 别揉我奶头~嗯~啊~动态视频| 好男人电影高清在线观看| 非洲黑人性xxxx精品又粗又长| 国产精品久久久久久亚洲av鲁大| 蜜桃亚洲精品一区二区三区| 日韩欧美一区二区三区在线观看| 国产色婷婷99| a级毛片a级免费在线| 搡老熟女国产l中国老女人| 听说在线观看完整版免费高清| 日韩国内少妇激情av| 男女视频在线观看网站免费| 一区二区三区激情视频| 51国产日韩欧美| 国产99白浆流出| 久9热在线精品视频| 国产爱豆传媒在线观看| 久久久国产成人精品二区| avwww免费| 最近最新中文字幕大全电影3| 哪里可以看免费的av片| 婷婷精品国产亚洲av| 色精品久久人妻99蜜桃| 悠悠久久av| 亚洲欧美激情综合另类| 最新中文字幕久久久久| 精品一区二区三区av网在线观看| 国产精品日韩av在线免费观看| 日本一本二区三区精品| 欧美性猛交黑人性爽| 内射极品少妇av片p| 男女视频在线观看网站免费| 久久精品91蜜桃| 一进一出抽搐动态| 国产午夜精品久久久久久一区二区三区 | 国产伦在线观看视频一区| 欧美黄色片欧美黄色片| 亚洲av不卡在线观看| 久久久久久久久大av| 99精品在免费线老司机午夜| 超碰av人人做人人爽久久 | 制服丝袜大香蕉在线| 老司机在亚洲福利影院| 毛片女人毛片| 女人高潮潮喷娇喘18禁视频| 国产成人啪精品午夜网站| 美女高潮喷水抽搐中文字幕| 麻豆一二三区av精品| 精品免费久久久久久久清纯| 中文亚洲av片在线观看爽| 亚洲人成电影免费在线| 美女免费视频网站| xxx96com| 乱人视频在线观看| av福利片在线观看| 禁无遮挡网站| 亚洲精品一区av在线观看| 精品久久久久久成人av| 亚洲熟妇中文字幕五十中出| 一本综合久久免费| 最后的刺客免费高清国语| 1000部很黄的大片| 老司机午夜福利在线观看视频| 香蕉av资源在线| 天堂av国产一区二区熟女人妻| 一夜夜www| 在线天堂最新版资源| 热99re8久久精品国产| 亚洲精品日韩av片在线观看 | 人妻夜夜爽99麻豆av| xxxwww97欧美| 久久久久久久久大av| 天堂av国产一区二区熟女人妻| 日本a在线网址| 女同久久另类99精品国产91| 女警被强在线播放| 一区福利在线观看| 成人欧美大片| 丝袜美腿在线中文| 12—13女人毛片做爰片一| 少妇的丰满在线观看| 国产私拍福利视频在线观看| netflix在线观看网站| 日本 欧美在线| 午夜福利高清视频| 99久久综合精品五月天人人| 国产高清视频在线观看网站| 欧美绝顶高潮抽搐喷水| 欧美乱妇无乱码| 村上凉子中文字幕在线| 九色国产91popny在线| 天天添夜夜摸| 99视频精品全部免费 在线| 国产精品美女特级片免费视频播放器| 色噜噜av男人的天堂激情| 中文字幕久久专区| 狠狠狠狠99中文字幕| 国内精品美女久久久久久| 欧美中文日本在线观看视频| 国产成年人精品一区二区| 亚洲一区二区三区不卡视频| 嫁个100分男人电影在线观看| 人人妻,人人澡人人爽秒播| 国产精品久久电影中文字幕| 成年女人永久免费观看视频| h日本视频在线播放| 日韩有码中文字幕| 久久亚洲精品不卡| 小蜜桃在线观看免费完整版高清| 嫩草影院入口| 国产精品影院久久| 亚洲天堂国产精品一区在线| 黄色视频,在线免费观看| 欧美区成人在线视频| aaaaa片日本免费| 亚洲激情在线av| 在线观看舔阴道视频| 性色avwww在线观看| 国产av在哪里看| 国产综合懂色| 国产伦人伦偷精品视频| av在线蜜桃| 少妇的丰满在线观看| 波多野结衣高清作品| 又黄又粗又硬又大视频| 国产亚洲av嫩草精品影院| 日韩成人在线观看一区二区三区| 久久久精品大字幕| 亚洲欧美日韩卡通动漫| 国产乱人视频| 久久久久精品国产欧美久久久| 最近视频中文字幕2019在线8| 国产视频一区二区在线看| 国产亚洲欧美98| 久久久久久久久久黄片| 欧美绝顶高潮抽搐喷水| 亚洲国产色片| 久久这里只有精品中国| 成人精品一区二区免费| 日本 欧美在线| 人人妻,人人澡人人爽秒播| 一本精品99久久精品77| 午夜亚洲福利在线播放| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成网站在线播| 丰满乱子伦码专区| 亚洲专区国产一区二区| 国产亚洲欧美98| 国产精品一区二区三区四区久久| 日韩精品中文字幕看吧| 在线观看免费午夜福利视频| 国产蜜桃级精品一区二区三区| 91在线观看av| 毛片女人毛片| 久久精品国产自在天天线| 国产激情欧美一区二区| 一个人观看的视频www高清免费观看| 国产三级中文精品| 国产色爽女视频免费观看| 免费看美女性在线毛片视频| 精品一区二区三区人妻视频| 亚洲国产中文字幕在线视频| 国产一区二区三区在线臀色熟女| 久久精品91蜜桃| 人妻久久中文字幕网| 九九热线精品视视频播放| 全区人妻精品视频| 99久久综合精品五月天人人| 久久精品综合一区二区三区| 国产亚洲精品久久久久久毛片| 狂野欧美激情性xxxx| 亚洲内射少妇av| 日本成人三级电影网站| 少妇的丰满在线观看| 有码 亚洲区| 精品不卡国产一区二区三区| 成年女人毛片免费观看观看9| 麻豆国产97在线/欧美| 婷婷精品国产亚洲av| 美女高潮的动态| 国产成人福利小说| 欧美区成人在线视频| 99久久综合精品五月天人人| 欧美日韩中文字幕国产精品一区二区三区| 欧美区成人在线视频| 99久国产av精品| 午夜福利免费观看在线| 国产午夜福利久久久久久| 亚洲国产中文字幕在线视频| 欧美性感艳星| 国产高潮美女av| 日韩欧美国产在线观看| 欧美成人性av电影在线观看| 此物有八面人人有两片| 观看免费一级毛片| 日韩av在线大香蕉| 欧美又色又爽又黄视频| 国产色婷婷99| 一边摸一边抽搐一进一小说|