Production of infectious bovine papillomavirus from cloned viral DNA by using an organotypic raft/xenograft technique

Bovine papillomavirus type 1 (BPV-1) induces fibropapillomas in its natural host and can transform fibroblasts in culture. The viral genome is maintained as an episome within fibroblasts, which has allowed extensive genetic analyses of the viral functions required for DNA replication, gene expression, and transformation. Much less is known about BPV-1 gene expression and replication in bovine epithelial cells because the study of the complete viral life cycle requires an experimental system capable of generating a fully differentiated stratified bovine epithelium. Using a combination of organotypic raft cultures and xenografts on nude mice, we have developed a system in which BPV-1 can replicate and produce infectious viral particles. Organotypic cultures were established with bovine keratinocytes plated on a collagen raft containing BPV-1-transformed fibroblasts. These keratinocytes were infected with virus particles isolated from a bovine wart or were transfected with cloned BPV-1 DNA. Several days after the rafts were lifted to the air interface, they were grafted on nude mice. After 6–8 weeks, large xenografts were produced that exhibited a hyperplastic and hyperkeratotic epithelium overlying a large dermal fibroma. These lesions were strikingly similar to a fibropapilloma caused by BPV-1 in the natural host. Amplified viral DNA and capsid antigens were detected in the suprabasal cells of the epithelium. Moreover, infectious virus particles could be isolated from these lesions and quantitated by a focus formation assay on mouse cells in culture. Interestingly, analysis of grafts produced with infected and uninfected fibroblasts indicated that the fibroma component was not required for productive infection or morphological changes characteristic of papillomavirus-infected epithelium. This system will be a powerful tool for the genetic analysis of the roles of the viral gene products in the complete viral life cycle.

Whole genome amplification using a degenerate oligonucleotide primer allows hundreds of genotypes to be performed on less than one nanogram of genomic DNA

Genetic analysis of limiting quantities of genomic DNA play an important role in DNA forensics, paleoarcheology, genetic disease diagnosis, genetic linkage analysis, and genetic diversity studies. We have tested the ability of degenerate oligonucleotide primed polymerase chain reaction (DOP-PCR) to amplify picogram quantities of human genomic DNA for the purpose of increasing the amount of template for genotyping with microsatellite repeat markers. DNA was uniformly amplified at a large number of typable loci throughout the human genome with starting template DNAs from as little as 15 pg to as much as 400 ng. A much greater-fold enrichment was seen for the smaller genomic DOP-PCRs. All markers tested were amplified from starting genomic DNAs in the range of 0.6–40 ng with amplifications of 200- to 600-fold. The DOP-PCR-amplified genomic DNA was an excellent and reliable template for genotyping with microsatellites, which give distinct bands with no increase in stutter artifact on di-, tri-, and tetranucleotide repeats. There appears to be equal amplification of genomic DNA from 55 of 55 tested discrete microsatellites implying near complete coverage of the human genome. Thus, DOP-PCR appears to allow unbiased, hundreds-fold whole genome amplification of human genomic DNA for genotypic analysis.

Targeted development of informative microsatellite (SSR) markers

We describe a novel approach, selectively amplified microsatellite (SAM) analysis, for the targeted development of informative simple sequence repeat (SSR) markers. A modified selectively amplified microsatellite polymorphic loci assay is used to generate multi-locus SSR fingerprints that provide a source of polymorphic DNA markers (SAMs) for use in genetic studies. These polymorphisms capture the repeat length variation associated with SSRs and allow their chromosomal location to be determined prior to the expense of isolating and characterising individual loci. SAMs can then be converted to locus-specific SSR markers with the design and synthesis of a single primer specific to the conserved region flanking the repeat. This approach offers a cost-efficient and rapid method for developing SSR markers for predetermined chromosomal locations and of potential informativeness. The high recovery rate of useful SSR markers makes this strategy a valuable tool for population and genetic mapping studies. The utility of SAM analysis was demonstrated by the development of SSR markers in bread wheat.

Direct detection and isolation of restriction landmark genomic scanning (RLGS) spot DNA markers tightly linked to a specific trait by using the RLGS spot-bombing method.

We have developed a technique for isolating DNA markers tightly linked to a target region that is based on RLGS, named RLGS spot-bombing (RLGS-SB). RLGS-SB allows us to scan the genome of higher organisms quickly and efficiently to identify loci that are linked to either a target region or gene of interest. The method was initially tested by analyzing a C57BL/6-GusS mouse congenic strain. We identified 33 variant markers out of 10,565 total loci in a 4.2-centimorgan (cM) interval surrounding the Gus locus in 4 days of laboratory work. The validity of RLGS-SB to find DNA markers linked to a target locus was also tested on pooled DNA from segregating backcross progeny by analyzing the spot intensity of already mapped RLGS loci. Finally, we used RLGS-SB to identify DNA markers closely linked to the mouse reeler (rl) locus on chromosome 5 by phenotypic pooling. A total of 31 RLGS loci were identified and mapped to the target region after screening 8856 loci. These 31 loci were mapped within 11.7 cM surrounding rl. The average density of RLGS loci located in the rl region was 0.38 cM. Three loci were closely linked to rl showing a recombination frequency of 0/340, which is < 1 cM from rl. Thus, RLGS-SB provides an efficient and rapid method for the detection and isolation of polymorphic DNA markers linked to a trait or gene of interest.

Fine mapping of the tomato I-3 gene for fusarium wilt resistance and elimination of a co-segregating resistance gene analogue as a candidate for I-3

The I-3 gene from the wild tomato species Lycopersicon pennellii confers resistance to race 3 of the devastating vascular wilt pathogen Fusarium oxysporum f. sp. lycopersici. As an initial step in a positional cloning strategy for the isolation of I-3, we converted restriction fragment length polymorphism and conserved orthologue set markers, known genes and a resistance gene analogue (RGA) mapping to the I-3 region into PCR-based sequence characterised amplified region (SCAR) and cleaved amplified polymorphic sequence (CAPS) markers. Additional PCR-based markers in the I-3 region were generated using the randomly amplified DNA fingerprinting (RAF) technique. SCAR, CAPS and RAF markers were used for high-resolution mapping around the I-3 locus. The I-3 gene was localised to a 0.3-cM region containing a RAF marker, eO6, and an RGA, RGA332. RGA332 was cloned and found to correspond to a putative pseudogene with at least two loss-of-function mutations. The predicted pseudogene belongs to the Toll interleukin-1 receptor-nucleotide-binding site-leucine-rich-repeat sub-class of plant disease resistance genes. Despite the presence of two RGA332 homologues in L. esculentum, DNA gel blot and PCR analysis suggests that no other homologues are present in lines carrying I-3 that could be alternative candidates for the gene.

A novel multilocus sequence typing scheme for the opportunistic pathogen Propionibacterium acnes and characterisation of type 1 cell surface-associated antigens

We have developed a novel multilocus sequence typing (MLST) scheme and database (http://pubmlst.org/pacnes/) for Propionibacterium acnes based on the analysis of seven core housekeeping genes. The scheme, which was validated against previously described antibody, single locus and random amplification of polymorphic DNA typing methods, displayed excellent resolution and differentiated 123 isolates into 37 sequence types (STs). An overall clonal population structure was detected with six eBURST groups representing the major clades I, II and III, along with two singletons. Two highly successful and global clonal lineages, ST6 (type IA) and ST10 (type IB1), representing 64?% of this current MLST isolate collection were identified. The ST6 clone and closely related single locus variants, which comprise a large clonal complex CC6, dominated isolates from patients with acne, and were also significantly associated with ophthalmic infections. Our data therefore support an association between acne and P. acnes strains from the type IA cluster and highlight the role of a widely disseminated clonal genotype in this condition. Characterization of type I cell surface-associated antigens that are not detected in ST10 or strains of type II and III identified two dermatan-sulphate-binding proteins with putative phase/antigenic variation signatures. We propose that the expression of these proteins by type IA organisms contributes to their role in the pathophysiology of acne and helps explain the recurrent nature of the disease. The MLST scheme and database described in this study should provide a valuable platform for future epidemiological and evolutionary studies of P. acnes.

The use of high resolution melting (HRM) to map single nucleotide polymorphism markers linked to a covered smut resistance gene in barley

Using an established genetic map, a single gene conditioning covered smut resistance, Ruh.7H, was mapped to the telomere region of chromosome 7HS in an Alexis/Sloop doubled haploid barley population. The closest marker to Ruh.7H, abg704 was 7.5 cM away. Thirteen loci on the distal end of 7HS with potential to contain single nucleotide polymorphisms (SNPs) were identified by applying a comparative genomics approach using rice sequence data. Of these, one locus produced polymorphic co-dominant bands of different size while two further loci contained SNPs that were identified using the recently developed high resolution melting (HRM) technique. Two of these markers flanked Ruh.7H with the proximal marker located 3.8 cM and the distal marker 2.7 cM away. This is the first report on the application of the HRM technique to SNP detection and to rapid scoring of known cleaved amplified polymorphic sequence (CAPS) markers in plants. This simple, precise post-PCR technique should find widespread use in the fine-mapping of genetic regions of interest in complex cereal and other plant genomes.

The development of an improved PCR-based marker system for Sw-5, an important TSWV resistance gene of tomato

Sw-5 is an important disease resistance gene of tomato, providing broad resistance to Tomato spotted wilt virus (TSWV). A cleaved amplified polymorphic sequence (CAPS) marker, closely linked to the gene, has been reported. Although the Sw-5 locus has been characterised, a gene-specific marker has not been developed. This paper presents a PCR-based marker-system that consists of the co-amplification of a dominant marker representing the Sw-5 gene sequence, and the modified CAPS marker as a positive control and indicator of genotype.

A consensus genetic map of sorghum that integrates multiple component maps and high-throughput Diversity Array Technology (DArT) markers

Background: Sorghum genome mapping based on DNA markers began in the early 1990s and numerous genetic linkage maps of sorghum have been published in the last decade, based initially on RFLP markers with more recent maps including AFLPs and SSRs and very recently, Diversity Array Technology (DArT) markers. It is essential to integrate the rapidly growing body of genetic linkage data produced through DArT with the multiple genetic linkage maps for sorghum generated through other marker technologies. Here, we report on the colinearity of six independent sorghum component maps and on the integration of these component maps into a single reference resource that contains commonly utilized SSRs, AFLPs, and high-throughput DArT markers. Results: The six component maps were constructed using the MultiPoint software. The lengths of the resulting maps varied between 910 and 1528 cM. The order of the 498 markers that segregated in more than one population was highly consistent between the six individual mapping data sets. The framework consensus map was constructed using a "Neighbours" approach and contained 251 integrated bridge markers on the 10 sorghum chromosomes spanning 1355.4 cM with an average density of one marker every 5.4 cM, and were used for the projection of the remaining markers. In total, the sorghum consensus map consisted of a total of 1997 markers mapped to 2029 unique loci ( 1190 DArT loci and 839 other loci) spanning 1603.5 cM and with an average marker density of 1 marker/0.79 cM. In addition, 35 multicopy markers were identified. On average, each chromosome on the consensus map contained 203 markers of which 58.6% were DArT markers. Non-random patterns of DNA marker distribution were observed, with some clear marker-dense regions and some marker-rare regions. Conclusion: The final consensus map has allowed us to map a larger number of markers than possible in any individual map, to obtain a more complete coverage of the sorghum genome and to fill a number of gaps on individual maps. In addition to overall general consistency of marker order across individual component maps, good agreement in overall distances between common marker pairs across the component maps used in this study was determined, using a difference ratio calculation. The obtained consensus map can be used as a reference resource for genetic studies in different genetic backgrounds, in addition to providing a framework for transferring genetic information between different marker technologies and for integrating DArT markers with other genomic resources. DArT markers represent an affordable, high throughput marker system with great utility in molecular breeding programs, especially in crops such as sorghum where SNP arrays are not publicly available.

Diagnostic Molecular Markers for Phosphine Resistance in U.S. Populations of Tribolium castaneum and Rhyzopertha dominica

Stored product beetles that are resistant to the fumigant pesticide phosphine (hydrogen phosphide) gas have been reported for more than 40 years in many places worldwide. Traditionally, determination of phosphine resistance in stored product beetles is based on a discriminating dose bioassay that can take up to two weeks to evaluate. We developed a diagnostic cleaved amplified polymorphic sequence method, CAPS, to detect individuals with alleles for strong resistance to phosphine in populations of the red flour beetle, Tribolium castaneum, and the lesser grain borer, Rhyzopertha dominica, according to a single nucleotide mutation in the dihydrolipoamide dehydrogenase (DLD) gene. We initially isolated and sequenced the DLD genes from susceptible and strongly resistant populations of both species. The corresponding amino acid sequences were then deduced. A single amino acid mutation in DLD in populations of T.castaneum and R.dominica with strong resistance was identified as P45S in T.castaneum and P49S in R.dominica, both collected from northern Oklahoma, USA. PCR products containing these mutations were digested by the restriction enzymes MboI and BstNI, which revealed presence or absence, respectively of the resistant (R) allele and allowed inference of genotypes with that allele. Seven populations of T.castaneum from Kansas were subjected to discriminating dose bioassays for the weak and strong resistance phenotypes. Application of CAPS to these seven populations confirmed the R allele was in high frequency in the strongly resistant populations, and was absent or at a lower frequency in populations with weak resistance, which suggests that these populations with a low frequency of the R allele have the potential for selection of the strong resistance phenotype. CAPS markers for strong phosphine resistance will help to detect and confirm resistant beetles and can facilitate resistance management actions against a given pest population.

Population connectivity among Dry Tortugas, Florida, and Caribbean populations of mutton snapper (Lutjanus analis), inferred from multiple microsatellite loci

Determining patterns of population connectivity is critical to the evaluation of marine reserves as recruitment sources for harvested populations. Mutton snapper (Lutjanus analis) is a good test case because the last known major spawning aggregation in U.S. waters was granted no-take status in the Tortugas South Ecological Reserve (TSER) in 2001. To evaluate the TSER population as a recruitment source, we genotyped mutton snapper from the Dry Tortugas, southeast Florida, and from three locations across the Caribbean at eight microsatellite loci. Both Fstatistics and individual-based Bayesian analyses indicated that genetic substructure was absent across the five populations. Genetic homogeneity of mutton snapper populations is consistent with its pelagic larval duration of 27 to 37 days and adult behavior of annual migrations to large spawning aggregations. Statistical power of future genetic assessments of mutton snapper population connectivity may benefit from more comprehensive geographic sampling, and perhaps from the development of less polymorphic DNA microsatellite loci. Research where alternative methods are used, such as the transgenerational marking of embryonic otoliths with barium stable isotopes, is also needed on this and other species with diverse life history characteristics to further evaluate the TSER as a recruitment source and to define corridors of population connectivity across the Caribbean and Florida.

银杉的遗传多样性及系统位置的研究一兼论松科的分子系统学

银杉属Cathaya是我国著名植物学家陈焕镛和匡可任建立的松科单型属，该属是否成立及该属在松科中的系统位置至今仍众说不一。银杉Cathazyaargyro phylla是银杉属的唯一幸存种，被誉为第三纪的活化石，并被列为我国一级保护的濒危植物，也是世界植物红皮书所列少数几种濒危裸子植物之一，保护区虽已设立数个，但该物种的濒危原因仍为未解之谜。本文针对上述问题开展了松科植物的分子系统学研究，探讨了银杉属的系统位置及松科的属间关系，并通过遗传多样性的研究探讨了银杉的濒危机制及保护策略。 1．松科的分子系统学研究 为选取合适的分子标记，本文先探讨了ITS-1在松科等裸子植物中的系统学价值。对1 6种裸子植物（其中含松科10属12种）的ITS-1的扩增结果如下：ITS-1在裸子植物属闻存在极大的长度变异(约600-2600bp)，因而不适用于属级以上的系统学研究;在同属种中，ITs-1的长度虽较一致，但与被子植物的ITS-1相比，裸子植物的ITS-1太长(如在松属中超过2600bp)，难以进行大规模研究。银杉属ITS-1的扩增产物中，一个片段的长度与松属ITS-1的长度相当，另一个片段的长度与云杉属ITS-l的长度相当，这三属间的关系值得进一步研究。此外，在松科植物ITS-1的研究中必须注意真菌“污染”这一问题。 由于ITS不适用于裸子植物属间关系的研究，因而本文运用rbcL-accD基因片段及trnK基因的PCR-RFLP分析探讨银杉属的系统位置及松科的属间关系。rbcL和trnK虽然都是叶绿体基因，但位于trnK内含子区的natK基因与rbcL基因在所受的选择压及进化速率上相差甚远，通过这两个基因研究结果的比较有助于评价系统树的可靠性。此外，考虑到松科的10个属中有9个属的rbcL基因序列已被测定，我们补测了银杉属的rbcL基因序列，进而对松科10个属的rbcL基因序列进行了全面分析。PCR -RFLP分析及rbcL基因的序列分析结果如下： A. PCR扩增的松科rbcL-accD基因片段均长约2550bp，属间没有明显的长度变异，1 8种限制性内切酶的单酶切及7种组合双酶切共获得86个酶切位点，其中54个为变异位点。PAUP和MEGA软件分析均表明：银杉属、松属、黄杉属和落叶松属构成一个单系群，且银杉属与松属的关系似乎近于它与另外二属的关系，但这一结果来得到bootstrap分析的较强支持;黄杉属近缘于落叶松属;冷杉属近缘于油杉属。 B．PCR扩增出的松科trnK基因均长约2557bp，在属间也没有明显的长度变异。16种限制性内切酶的单酶切及7种组合双酶切共获得92个酶切位点，其中68个为变异位点。变异位点中的35个具有系统发育信息，其中28个位于matK基因上。分布于整个trnK基因上的informative位点和仅分布于matK基因上的informative位点分别用于PAUP(version 3.1.1)分析，且PAUP分析中分别用Wagner简约法和Dollo筒约法构建系统树。总体看来，Wager简约树和Dono简约树的结构基本一致（Cedrus的位置除外）：Abies、Keteleeria、Tsuga和Pseudolarix分为一支，且这一支在所有的简约树中均得到分辨;Pseudotsuga近缘于Larix,而Abkes相对近缘于Keteleeria; Cathaya与Ables-Keteleeria-Pseudolarix-Tsuga这一支的关系较远，它与松科的其它属聚在一起，但b∞tstrap分析也不支持它与任何一个属有很近的关 系，说明Cathaya是比较孤立的一个属。Cedrus的位置较特殊，在Wagner简约树中它位于松属以外的其它8个属的基部（因松属为外类群），而在Dollo筒约树中它位于Al}ies-Kete~eer uz-Pseudo/zr/x- Tsuga达一支的基部，并且它与这一支的关系得到了bootstrap分析的较强支持。 C．银杉的rbeL基因长1425bp，其序列与松科其它9个属的rbcL基因序列均有较大差异，从遗传距离(Mega软件中的P-distance)上看，银杉属作为一个属是无可非议的。PAUp和Mega软件分析均表现：Abies、Keteleeria、Pseudolter和Tsuga分为一支，且Cathaya与这一支的关系较远;Pseuaiotsuga近缘于Larix,且这二属间的关系得到了bootstrap分析的强烈支持;Abies近缘于Keteleeria，Pseudolarix近缘于Tsuga，且这两支也得到了bootstrap分析的较强支持;Cathaya与Pinus聚为一支，但支持强度不高。Cetirus的位置与trnK基因的PCR-RFLP分析结果极为相似。 根据rbcL-accD基因片段及trnK基因的PCR-RFLP分析以及rbcL基因的序列分析可以得出如下结论：<1> Cathaya作为一个属不仅是成立的，而且是一个较为孤立的属，将它置于其它任何属内都是不正确的，该属与Abies-Keteleer ia-Pseudolarx-Tsuga这一支的关系较远，相对而言Cathaya与Pinus的关系可能稍近一些。<2>松科主要分为两大支：一支含Abies、KeteZeera、Pseudolarx和Tsuga，该支中Abies近缘于KetelaeriaPseudolarix近缘于Tsuga;另一支含Pseudotsuga、Larix、Cathaya、Pinus和Picete，该支中Pseudotsuga明显近缘于Larix;Cedrus可能更近于Abies-Keteleeria-Pseudolarix-Tsuga这一支。 2．银杉的遗传多样性研究 本研究运用随机扩增多态DNA(RAPD)方法对银杉的1 1 3个个体（分别采自现存全部4个居群：金佛山、花坪、大瑶山和八西山)进行了遗传多样性检测。21个1O-mer的寡核苷酸引物共检测11个位点，其中42个位点是多态的，占37%，这一多态位点百分率仅为裸子植物多态位点百分率平均值的一半左右，说明银杉的遗传多样性水平很低。对多态位点的分析发现银杉屠群内发生了十分严重的遗传漂变，因而导致屠群间的强烈分化。AMOVA分析表明：银杉的遗传变异中，34.7%的变异存在于居群闻，这一数值是裸子植物GST平均值(6。8%)的5倍多;在金佛山屠群中，17%的遗传变异存在于亚居群间，这一数值也是裸子植物GsT平均值的两倍多。银杉的这一居群遗传结构与其它裸子植物的居群遗传结构截然不同。此外，我们提出了度量遗传多样性水平的分化指数概念及其计算方法，并发现遗传多样性水平的高低与生境的复杂程度有一定的相关性。 银杉的遗传多样性水平很低，其适应幅必然很窄，且由于严重的遗传漂变导致居群间的强烈分化，使基因流受阻，进而产生严重的近交，近交又使银杉的适应能力进一步下降，并使极为有限的遗传多样性进一步丧失，这种恶性循环会导致银杉濒危程度的加剧，甚至绝灭。 鉴于银杉独特的居群遗传结构，即有相当大一部分遗传变异存在于居群之间，取样保护时，不仅要在每个居群中取足够多的个体，而且要在尽可能多的居群中取样。就原地保护而言，由于每个居群和亚居群都有独特的基因型，因而都具有保护价值。另外，人工加强居群间的基因流（如人工授粉）对提高银杉各居群的适应能力肯定是有益的。

中国普通野生稻（Oryza rufipogon Griff.）群体遗传学及保护策略的研究

普通野生稻的基因资源对水稻的育种已经起过并将继续起至关重要的作用。所以，研究和保存其遗传多样性具有重要的意义。首先对普通野生稻在我国的分布现状作了较为全面的野外调查。采集了69个群体，近2000个个体的硅胶干燥叶片。研究了从硅胶干燥的小量叶片中高效率、高质量制备DNA的方法。并用此方法制备了44个代表群体，1168个个体的总DNA，建立了中国普通野生稻的总DNA库。随机选取19个引物对29个群体、651个个体的总DNA进行了RAPD扩增。统计出群体和地区内的多态性片段百分比（PPB），RAPDistance和WINAMOVA程序联用分析遗传变异在地区间、地区内群体间和群体内的组成。广西濠江流域5个野生群体的分析结果表明了中下游群体的PPB值显著高于中上游群体。从而从分子水平上证实了水流对无性繁殖体的传播影响了沿江分布群体的遗传多样性和遗传结构。6个现有原位保护群体的分析结果说明了5个群体值得进一步原位保护，而另一个群体适合异位保护。中国境内29个群体的统计分析揭示出群体间的遗传多样性差异很大（PPB值从23.29%到46.18%）。广东具有的PPB值最高（79.12%），其次是海南（77.11%）和广西（68.67%）。因而，华南是我国普通野生稻的遗传多样性中心。18个群体AMOVA分析结果均表明了遗传变异主要存在于群体内，但群体间业已有较大的遗传分化。从175个随机引物中筛选出9个引物，成功地鉴定出5个群体的克隆多样性和克隆结构。这二者与群体的生态环境，特别是水分状况密切相关。也同时受人类干扰强度的影响。而与群体所处的纬度没有明显的相关性。 根据以上结果，普通野生稻原位保护的策略是选取分布在不同地理区域、遗传多样性水平高的群体。分布在广东、海南和广西的6个群体，具有很高的遗传多样性，建议加以原位保护。此外，孤岛状隔离在我国普通野生稻分布区边缘的3个小群体也应立即进行原位保护。此外，原位保护的群体应给予适度的干扰以维系较高的克隆多样性和遗传多样性。在进行遗传多样性研究和异位保护时，每群体的取样个体数应不少于25株，相邻取样个体间的距离应大于12米。

疣粒野生稻居群遗传结构及其形成与维持机制的研究

结合野外生态调查与分子标记检测，本文探讨了显性遗传标记应用于居群遗传学研究时实验与数据分析需要注意的问题，并在此基础上对中国分布的疣粒野生稻（Oryza granulata Nees et Am. ex Watt.）居群遗传多样性与居群遗传结构进行了研究。然后从metapopulation结构动态、无性系生长和物种形成等角度研究了遗传结构的空间格局与生态学、系统学因素之间的相互作用，并将研究的空间尺度从聚群内（colony）、居群内、地区内推广到地区间和整个物种水平，反映了在不同空间等级上疣粒野生稻相异的进化模式。最后，综合上述结果，提出了保护疣粒野生稻的原则和策略。结果如下： 1．根据对分布于中国云南和海南33个分布点疣粒野生稻居群所做的野外生态学调查，该物种目前在中国的30个县（市）有分布，比1978-1982年全国野生稻普查时增加了3个县市（海南通什、云南思茅和勐腊）。疣粒野生稻具有较强的耐荫和抗旱能力，在群落总盖度为90-210%范围内生长良好。它是一种典型的适应中度干扰的物种，生长于有一定干扰的斑块状生境中。疣粒野生稻在群落内的分布格局为聚集型，其居群密度较小（1.13-2.95株/m2），依靠重力下落和动物传播种子。由于对热区资源的掠夺性开发，总计有12.9%的居群已因人为干扰而灭绝，83.9%的居群处于严重的危胁之下，处于濒危状态。对疣粒野生稻的破坏在不同地区间程度不同，生境恶化和放牧是造成其居群灭绝的最主要原因，对其进行保护已经迫在眉睫。在调查的基础上，本研究建立了含该物种34个居群共1109份个体样品的总DNA库，作为易位保护的一种措施，主要用于居群遗传学和保护生物学研究。 2．利用上述总DNA库中的材料，首先采用随机扩增多态（RAPD）对几类显性标记的居群遗传结构参数进行了比较。在衡量遗传多样性水平时，多态位点比率（PPB）会低估变异的程度，其价值不如Shannon多样性指数和Nei多样性指数。在计算个体之间的遗传关系时，Mantel检测表明17种相似性系数之间存在极显著的相关性。同时，基于Φst。遗传距离的分子方差分析（AMOVA）和基于Hardy-Weinberg平衡假设的Nei氏遗传距离分析的结果间具有显著的相关性，它们都适用于对疣粒野生稻居群遗传结构进行研究，且在使用后者时应对数据进行Lynch-Milligan矫正，剔除隐性基因型（0）频率小于3/N（N为样本总数）的条带数据，以矫正显性遗传方式对变异估计偏低的影响。此外，各类遗传结构分析参数之间的高度相关性也与疣粒野生稻居群内遗传多样性低，杂合体比率较低有密切关系。 利用RAPD和inter-简单重复序列（ISSR）对来自中国20个居群疣粒野生稻混合样品，以及海南（M5）和云南（M27）两个居群各20个植株的遗传多样性进行了检测。ISSR的实验稳定性优于RAPD，且总的来说它能检测到更多的遗传变异。前者与它在PCR反应时退火温度较高，引物．模板复合物较稳定有关;而后者则与其引物靶序列容易在细胞分裂中产生突变有关。Mantel检测结果表明，衡量样品间的遗传关系时，这两种标记的分析结果在物种水平存在极显著的相关性（r = 0.917, t = 12789），而在居群水平不相关（r < 0.200）。这不但与它们所扩增的相应基因组片断的变异方式及在居群内分辩率下降有关，同时也反映了疣粒野生稻居群内和居群间存在着不同的进化模式。 3．利用RAPD和ISSR标记，对中国20个居群疣粒野生稻混合样品的遗传多样性进行了研究。RAPD和ISSR分别扩增出209和122条PCR条带，其中各有64.1l％（134条带）和72.95％（89条带）为多态条带。基于Jaccard系数的UPGMA分析表明，同一地区内居群的遗传变异比较小。20个居群按来源聚为云南与海南两类，其间产生了一定程度的遗传分化。这种遗传多样性分布的特点可能与其起源、分布格局、交配系统和种子散布方式有关。此外，尽管混合取样会低估一定的遗传变异能力，也不能得到关于居群内的遗传结构情况，但它仍然是一种获取遗传多样性信息的高效方法，适用于对研究材料进行日常管理和评价。 4．按照居群取样的方法，利用RAPD对来自中国云南和海南的20个疣粒野生稻居群共396个植株进行了遗传结构分析。联合ISSR，对其中5个居群初步的分析表明该物种在居群内的遗传多样性水平很低，RAPD的多态位点比率（PPB）在居群内从4.52％到13.06％;而ISSR的PPB值在居群内从7.08％-26.55％。AMOVA分析表明，对于RAPD来说，云南与海南两地区之间遗传变异的量占总变异量的73.85％，地区之内占19.45％，而居群内仅占6.70％。对于ISSR，疣粒野生稻地区间，地区内和居群内遗传变异的分布比率分别为49.26％、38.070A和12.66％。UPGMA聚类将同一居群内的个体聚为一支，并将居群按来源分为云南和海南两类。由于疣粒野生稻在群落内的分布为典型的metapopulation格局，伴随各聚群（colony）在群落次生演替过程中周转（灭绝与定植）时发生的遗传漂变、建立者效应和居群内强烈的近交是造成其居群内遗传多样性极低的主要原因。 5．利用10个ISSR标记对中国4个疣粒野生稻居群内的无性系生长与基因型遗传多样性进行了分析。在小尺度取样（个体间隔1.0-1.5m）的情况下，所有居群中均检测到明显的无性系生长现象。参与无性系生长的个体百分比在各居群中从25％-60％不等，Simpson多样性指数表明疣粒野生稻居群内的基因型多样性保持在较高水平（0.837-0.958）。尽管如此，AMOVA对居群内遗传变异进行方差剖分的结果表明参与无性生长的个体所含有的变异量平均只占总变异量的16.7%。因此，疣粒野生稻居群内遗传变异的来源主要依靠有性生殖来维持。同时，处于人为中度干扰之下的疣粒野生稻居群不但个体密度较高，其居群遗传多样性也未因此而降低。 6．在假设RAPD在同一种及其近缘种内PCR产物同源性较高的前提下，利用该技术对来自世界的23份O. granulata和O. meyeriana样品进行了遗传多样性分析和系统学研究。在物种水平，O. granulata具有非常高的遗传多样性（多态位点比率达83.54％），表明该物种进化历史中存在大规模居群瓶颈效应的可能性较小。O.gramulata与O. meyeriana各居群的遗传分化与岛屿形成导致的地理隔离之间有密切的关系。基于Nei & Li遗传相似性系数，利用Neighbor-Joining和UPGMA聚类法构建的两个系统树并不完全一致。主坐标分析（PCoA）支持NJ法的结果：来自O.meyeriana的两个样品倾向于聚为一类，并获得bootstrap分析的支持，但它们的遗传变异范围并未超出O. granulata。因此，我们的结果支持将这两个种进行归并。 7．由于疣粒野生稻在物种水平的遗传多样性非常高，变异主要存在于各地区之间，因此要最大程度地维持该物种遗传多样性，使之不发生遗传侵蚀意味着保护应该针对整个物种的分布区进行。对于分布于中国的居群来说，一方面由于变异主要存在于云南和海南两地区之间，另一方面由于地区内和居群内的遗传多样性相对较低，因此，云南和海南应做为中国保护该物种的两个中心。此外，由于一定程度的人为干扰有利于为该物种创造适宜生境，增加其居群密度，且不会致使遗传变异能力下降，因此在实施就地保护时应充分考虑将其与当地的经济开发项目相结合，达到自然保护与地方经济可持续发展的目的。

三种兰科植物的保护遗传学研究

本研究在野外调查的基础上，采用随机扩增多态DNA (RAPD)分析和形态学方法，研究了我国三种珍稀濒危兰科植物硬叶兜兰（Paphiopedilummicranthum）、麻栗坡兜兰（P. malipoense）和独花兰（Changnienia amoena）的遗传多样性与群体遗传结构，主要结果如下： 1．采用1 2个引物对分布于我国云贵地区的4个硬叶兜兰群体共161个体进行RAPD扩增和分析，得出物种水平的多态条带百分率（PPB）为71.6%，Nci的基因多样度(h)为0.217，Shannon多样性指数（1)为0.3301;4个群体的平均多样性水平为PPB=45.2%，h=0.1457，1= 0.2204：低于远交兰花的平均水平。分子方差分析(AMOVA)表明，在总遗传变异中，群体间遗传变异占20.31%．群体内占79.69%;POPGENE给出的基因分化系数 （Gst）为0.2958;遗传分化略高于远交物种的平均水平。空间自相关分析表明，所检测的两个群体中存在明显的空间结构，基因型在群体中以不同的小斑块存在。遗传距离和空间距离不存在相关关系。 2．用于麻栗坡兜兰的RAPD引物同上，但取样范围只有贵州的2个群体共10个个体。就所研究的个体柬看，麻栗坡兜兰的遗传多样性明显低于远交兰花物种的平均水平。物种水平上，多态条带百分率(PPB)为49.5%。Nei的基因多样度（h）为0. 1174， Shannon多样性指数(I)为0.1764：在群体水平上，上述三个指标的平均值则分别为12. 75%、0.0486和0.0712，均大大低于硬叶兜兰。然而，尽管作了种种努力，麻栗坡兜兰的取样个体数量仍很少，因此所得结果可能会有误差。 3．用16个引物对分布于河南、湖北、湖南、江西4个省11个独花兰群体共216个体进行了RA PD扩增和分析，独花兰在物种水平PPB=80. 7%，h=0.197．1=0. 3116;在群体水平，上述三个指标的平均值则分别为40. 9%、0.1247和0. 1902，均低于远交兰花的平均水平。AMOVA分析表明，11个独花兰群体间的遗传变异占43.48%，群体内的占56.52%：在神农架和新宁地区内部，群体间的遗传变异分别占13.68%和49.3g%(AMOVA)。POPGENE给出的11个群体的基因分化系数(Gst)为0.3580．神农架和新宁地区内的Gst，值分别为0.1194和0.2597。可见，群体间的遗传分化明显高于远交物种的平均水平。空间自相关分析表明，独花兰的遗传变异在群体内不存在明显的空间结构。群体之间的遗传距离和空间距离不存在相关关系。 4．对独花兰7个群体形态性状的分析发现，12个形态性状在群体内均有较高的变异性，cv值变动于0.022-0.30O。庐山群体(LS)在所有性状上的平均值均为最高。营养性状和花部性状的变异性基本一致。除花葶长和花距直径与某些花部性状之间没有显著的相关关系外，各性状之间均有显著的相关性。对XN4群体的统计没有发现假磷茎数目与其他性状之间存在显著相关性。 根据以上对硬叶兜兰、麻粟坡兜兰和独花兰遗传多样性和群体遗传结构韵研究，结合其他方面的资料;对三种兰花的濒危机制进行了初步的分析。首先，人为采挖和破坏是导致这些兰花物种濒危的直接原因，尤其是麻栗坡兜兰。其次， 适宜兰花生存的生境正在只益萎缩、退化和片段化。这两方面因素的共同作用导致上述兰花群体的数目和规模日益下降，由此引发的遗传多样性降低和遗传结构的改变进一步加剧其濒危状况。对于独花兰而言，较低的繁殖能力又使其生存状态雪上加霜。针对三个物种不同的繁殖特性和遗传学状况，提出如下保护措施。(1)硬叶兜兰由于繁殖能力较强、现存个体尚多，遗传多样性损失不甚严重，因此以保护其所在的生境为基础、实施原位保护，是比较合适的保护策略。(2)麻粟坡兜兰目前受破坏程度非常严重;所剩个体很少，遗传多样性较低，已经很难进行有效的原位保护。因此;应利用迁地保护手段抢救目前尚存的个体。(3)独花兰的繁殖能力较弱，因此在保护生境和严禁采摘的基础上，可采用人工授粉等方式，提高结实率、增加繁殖效率，促使其复壮：在进行迁地保护时，则应注意不同群体间存在较大遗传变异而群体内多样性较低这一现实。