基于ITS 序列和matK 基因序列的中国药用石斛及其混伪品的分子鉴定

为了从分子水平对中国药用石斛及其混伪品进行鉴定，本文选取了核rDNA ITS 序列和叶绿体DNA 的matK 基因序列进行研究。采用改良的CTAB 法提取石斛的基因组DNA，PCR 产物直接测序法对17 种（共32 份）药用石斛的核糖体内转录间隔区ITS 全序列进行测定，克隆测序法对12 种（共22 份）药用石斛的叶绿体的matK 基因序列进行测定，运用BioEd it，MEGA4.0 等生物软件分析了石斛属植物的rDNA ITS 序列及叶绿体的matK 基因序列的特征，比较了石斛属间、种间、种内不同居群（品种）间的序列碱基差异及遗传距离，应用邻接法构建分子系统树。主要研究结果如下： （1）建立了17 种（共32 份）药用石斛rDNA ITS 区碱基全序列数据库，其中，ITS1 的长度为228~234 bp，GC 含量为45.7%~53.0%，变异位点167 个，占总位点67.34%，信息位点106 个，占总位点42.74%，ITS2 长度为241~247 bp，GC含量为44.8%~55.7%，变异位点165 个，占总位点66.27%，信息位点115 个，占总位点46.18%。 （2）建立了12 种（共22 份）药用石斛的叶绿体matK 基因全序列数据库，叶绿体matK 基因长1410 bp，变异位点51 个，信息位点11 个。除了存在碱基替换的遗传变异外，还存在碱基的插入和缺失。 （3）通过ITS 序列比较分析了各材料间的遗传距离和碱基差异，属间的遗传距离为0.295，石斛种间的平均遗传距离为0.142，碱基相差2~156 个，种内各居群间的平均遗传距离为0.002，碱基相差1~2 个。属间的遗传距离大于种间的遗传距离，种间的遗传距离大于种内不同居群（品种）间的遗传距离。 （4）根据分析石斛叶绿体的matK 基因序列得到，外类群（密花石豆兰）与石斛属间最小遗传距离为0.027，石斛种间的平均遗传距离为0.008，种间最大的遗传距离0.014， 最小的遗传距离为0.003，碱基相差8~20 个。种内不同居群（品种）遗传距离为0.001，相差1~5 个碱基。 （5）利用17 种石斛的全序列数据库及遗传分析软件，通过对待检种rDNA I TS区进行序列测定，成功地对10 个待检种进行了鉴定，并且在原植物开花后得到了验证。 （6）运用12 种石斛的matK 基因全序列数据库及遗传分析软件，成功地对4个待检种进行了鉴定，同样在原植物开花后得到了验证。 （7）本文利用石斛的核糖体内转录间隔区ITS 序列和叶绿体的matK 基因序列数据库分别构建了NJ 树，外类群与石斛属间石斛种间以及种内不同居群（品种）间均能在NJ 树中明显分化开来，二者构建的分子系统树一致，为石斛的分子鉴定提供了依据。 In order to identify Chinese Herba Dendrobii and its adulterant species on molecular level, we studied the sequences of rDNA ITS and chloroplast matK gene. Genomic DNA of Dendrobium was extracted using the modified cetyltrimethyl ammonium bromide (CTAB) method. The PCR products of the rDNA ITS sequences of Dendrobium (32 materia ls) were purified and then sequenced. The PCR products of chloroplast matK gene of Dendrobium (22 materia ls) were purified, cloned and then sequenced. The characteristic of the sequences and the genetic dista nce were compared between Bulbophyllum odoratissimum and Dendrobium, Dendrobium interspecies, and different populations. Phylogenetic trees were constructed using the NJ method by the biology softwares including BioEd it, MEGA4.0 etc. The ma in results as follows: (1) It was built up that the database of rDNA ITS sequences of 17 species of Herba Dendrobii (32 materia ls). The ITS1 was 228~234 bp, the GC content accounting for 45.7%~53.0%. Its variable sites were 167, accounting for 67.34%. The Parsim-Informative positions were 106, accounting for 42.74%. The ITS2 was 241~247 bp, the GC accounting for 44.8%~55.7%. The variable sites were 165, accounting for 66.27%. The Parsim-Informative positions were 115, accounting for 46.18%. (2) The database of the chloroplast matK gene sequences was built up, which contained 12 species of Herba Dendrobii (22 materia ls). The matK gene sequences were about 1410bp in length. There were 51 variable sites and 11 Parsim-Informative sites. And there were nucleotides insertions and deletions in some species , in addition to the nucleotides substitutions. (3) The rDNA ITS sequences were compared and analyzed by the biology softwares. The genetic dista nce between Bulbophyllum odoratissimum and Dendrobium was 0.295. The avera ge genetic dista nce was 0.142 between Dendrobium species, and there were 2~156 variable nucleotides. The avera ge genetic dista nce between different populations was 0.002, and there were 2~156 variable nucleotides. The genetic dista nce between Bulbophyllum odoratissimum and Dendrobium was greater tha n that of Denrobium interspecies. Meanwhile, the genetic dista nce between Denrobium species was also greater tha n that of different populations (variaties). (4) The characteristics of the chloroplast matK gene sequences were obtained after analyzing by the biology softwares. The minima l genetic dista nce was 0.027 between Bulbophyllum odoratissimum and Dendrobium . The ma xima l genetic dista nce was 0.014 between Dendrobium species, and there were 20 variable nucleotides. The minima l genetic dista nce between populations was 0.003, and there were 8 variable nucleotides.The genetic dista nce between populations was 0.001, and there were 1~5 variable nucleotides. (5) The molecular Phylogeny tree was constructed on the database of rDNA ITS the sequences of 17 species of Herba Dendrobii using the biology softwares. Then we authenticated 10 materia ls on molecular level. What’s more, they had been proved when these pla nts flowered. (6) The molecular Phylogeny tree was built up on the database of chloroplast matK gene sequences of 12 species of Herba Dendrobii with the biology softwares.Then 4 materia ls were authenticated on molecular level. Moreover, they had also been proved when these pla nts were in flower. (7) The Phylogenetic trees were separately constructed on the sequences of rDNA ITS and chloroplast matK gene B. odoratissimum and Dendrobium all could be distinguished on the Phylogenetic trees. Meanwhile, the Phylogenetic trees based on two groups of sequences were coincident. rDNA ITS and matK gene sequence could be used as molecular markers for authentication of Herba Dendrobii.

Elusive relationships within order fabales: phylogenetic analyses using matK and rbcL sequence data

The order Fabales, including Leguminosae, Polygalaceae, Quillajaceae and Surianaceae, represents a novel hypothesis emerging from angiosperm molecular phylogenies. Despite good support for the order, molecular studies to date have suggested contradictory, poorly supported interfamilial relationships. Our reappraisal of relationships within Fabales addresses past taxon sampling deficiencies, and employs parsimony and Bayesian approaches using sequences from the plastid regions rbcL (166 spp.) and matK (78 spp.). Five alternative hypotheses for interfamilial relationships within Fabales were recovered. The Shimodaira-Hasegawa test found the likelihood of a resolved topology significantly higher than the one calculated for a polytomy, but did not favour any of the alternative hypotheses of relationship within Fabales. In the light of the morphological evidence available and the comparative behavior of rbcL and matK, the topology recovering Polygalaceae as sister to the rest of the order Fabales with Leguminosae more closely related to Quillajaceae + Surianaceae, is considered the most likely hypothesis of interfamilial relationships of the order. Dating of selected crown clades in the Fabales phylogeny using penalized likelihood suggests rapid radiation of the Leguminosae, Polygalaceae, and (Quillajaceae + Surianaceae) crown clades.

Filogenia e biogeografia de Droseraceae inferidas a partir de caracteres morfológicos e moleculares (18S, atpB, matK, rbcL e ITS)

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Amelia matką; dzieło za dalszy ciąg Pamiątki po dobréy matce służyć mogące.

t. 1. Część pierwsza, obeymuiąca ogólny obraz religii prawdziwey.

Molecular Characterization of Viruses Causing the Cassava Brown Streak Disease Epidemic in Eastern Africa

Cassava brown streak disease (CBSD) was described for the first time in Tanganyika (now Tanzania) about seven decades ago. Tanganyika (now Tanzania) about seven decades ago. It was endemic in the lowland areas of East Africa and inland parts of Malawi and caused by Cassava brown streak virus (CBSV; genus Ipomovirus; Potyviridae). However, in 1990s CBSD was observed at high altitude areas in Uganda. The causes for spread to new locations were not known.The present work was thus initiated to generate information on genetic variability, clarify the taxonomy of the virus or viruses associated with CBSD in Eastern Africa as well as to understand the evolutionary forces acting on their genes. It also sought to develop a molecular based diagnostic tool for detection of CBSD-associated virus isolates. Comparison of the CP-encoding sequences of CBSD-associated virus isolates collected from Uganda and north-western Tanzania in 2007 and the partial sequences available in Genbank revealed occurrence of two genetically distinct groups of isolates. Two isolates were selected to represent the two groups. The complete genomes of isolates MLB3 (TZ:Mlb3:07) and Kor6 (TZ:Kor6:08) obtained from North-Western (Kagera) and North-Eastern (Tanga) Tanzania, respectively, were sequenced. The genomes were 9069 and 8995 nucleotides (nt), respectively. They translated into polyproteins that were predicted to yield ten mature proteins after cleavage. Nine proteins were typical in the family Potyviridae, namely P1, P3, 6K1, CI, 6K2, VPg, NIa-Pro, NIb and CP, but the viruses did not contain HC-Pro. Interestingly, genomes of both isolates contained a Maf/HAM1-like sequence (HAM1h; 678 nucleotides, 25 kDa) recombined between the NIb and CP domains in the 3’-proximal part of the genomes. HAM1h was also identified in Euphorbia ringspot virus (EuRSV) whose sequence was in GenBank. The HAM1 gene is widely spread in both prokaryotes and eukaryotes. In yeast (Saccharomyces cerevisiae) it is known to be a nucleoside triphosphate (NTP) pyrophosphatase. Novel information was obtained on the structural variation at the N-termini of polyproteins of viruses in the genus Ipomovirus. Cucumber vein yellowing virus (CVYV) and Squash vein yellowing virus (SqVYV) contain a duplicated P1 (P1a and P1b) but lack the HC-Pro. On the other hand, Sweet potato mild mottle virus (SPMMV), has a single but large P1 and has HC-Pro. Both virus isolates (TZ:Mlb3:07 & TZ:Kor6:08) characterized in this study contained a single P1 and lacked the HC-Pro which indicates unique evolution in the family Potyviridae. Comparison of 12 complete genomes of CBSD-associated viruses which included two genomes characterized in this study, revealed genetic identity of 69.0–70.3% (nt) and amino acid (aa) identities of 73.6–74.4% at polyprotein level. Comparison was also made among 68 complete CP sequences, which indicated 69.0-70.3 and 73.6-74.4 % identity at nt and aa levels, respectively. The genetic variation was large enough for dermacation of CBSD-associated virus isolates into two distinct species. The name CBSV was retained for isolates that were related to CBSV isolates available in database whereas the new virus described for the first time in this study was named Ugandan cassava brown streak virus (UCBSV) by the International Committee on Virus Taxonomy (ICTV). The isolates TZ:Mlb3:07 and TZ:Kor6:08 belong to UCBSV and CBSV, respectively. The isolates of CBSV and UCBSV were 79.3-95.5% and 86.3-99.3 % identitical at nt level, respectively, suggesting more variation amongst CBSV isolates. The main sources of variation in plant viruses are mutations and recombination. Signals for recombination events were detected in 50% of isolates of each virus. Recombination events were detected in coding and non-coding (3’-UTR) sequences except in the 5’UTR and P3. There was no evidence for recombination between isolates of CBSV and UCBSV. The non-synonomous (dN) to synonomous (dS) nucleotide substitution ratio (ω) for the HAM1h and CP domains of both viruses were ≤ 0.184 suggesting that most sites of these proteins were evolving under strong purifying selection. However, there were individual amino acid sites that were submitted to adaptive evolution. For instance, adaptive evolution was detected in the HAM1h of UCBSV (n=15) where 12 aa sites were under positive selection (P< 0.05) but not in CBSV (n=12). The CP of CBSV (n=23) contained 12 aa sites (p<0.01) while only 5 aa sites in the CP gene of UCBSV were predicted to be submitted to positive selection pressure (p<0.01). The advantages offered by the aa sites under positive selection could not be established but occurrence of such sites in the terminal ends of UCBSV-HAMIh, for example, was interpreted as a requirement for proteolysis during polyprotein processing. Two different primer pairs that simultaneously detect UCBSV and CBSV isolates were developed in this study. They were used successfully to study distribution of CBSV, UCBSV and their mixed infections in Tanzania and Uganda. It was established that the two viruses co-infect cassava and that incidences of co-infection could be as high as 50% around Lake Victoria on the Tanzanian side. Furthermore, it was revealed for the first time that both UCBSV and CBSV were widely distributed in Eastern Africa. The primer pair was also used to confirm infection in a close relative of cassava, Manihot glaziovii (Müller Arg.) with CBSV. DNA barcoding of M. glaziovii was done by sequencing the matK gene. Two out of seven M. glaziovii from the coastal areas of Korogwe and Kibaha in north eastern Tanzania were shown to be infected by CBSV but not UCBSV isolates. Detection in M. glaziovii has an implication in control and management of CBSD as it is likely to serve as virus reservoir. This study has contributed to the understanding of evolution of CBSV and UCBSV, which cause CBSD epidemic in Eastern Africa. The detection tools developed in this work will be useful in plant breeding, verification of the phytosanitary status of materials in regional and international movement of germplasm, and in all diagnostic activities related to management of CBSD. Whereas there are still many issues to be resolved such as the function and biological significance of HAM1h and its origin, this work has laid a foundation upon which the studies on these aspects can be based.

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

银杉属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平均值的两倍多。银杉的这一居群遗传结构与其它裸子植物的居群遗传结构截然不同。此外，我们提出了度量遗传多样性水平的分化指数概念及其计算方法，并发现遗传多样性水平的高低与生境的复杂程度有一定的相关性。 银杉的遗传多样性水平很低，其适应幅必然很窄，且由于严重的遗传漂变导致居群间的强烈分化，使基因流受阻，进而产生严重的近交，近交又使银杉的适应能力进一步下降，并使极为有限的遗传多样性进一步丧失，这种恶性循环会导致银杉濒危程度的加剧，甚至绝灭。 鉴于银杉独特的居群遗传结构，即有相当大一部分遗传变异存在于居群之间，取样保护时，不仅要在每个居群中取足够多的个体，而且要在尽可能多的居群中取样。就原地保护而言，由于每个居群和亚居群都有独特的基因型，因而都具有保护价值。另外，人工加强居群间的基因流（如人工授粉）对提高银杉各居群的适应能力肯定是有益的。

低等金缕梅类的系统学：分子和发育学证据

传统的低等金缕梅类包含7个科，即领春木科（Eupteleaceae）、昆栏树科（Trochodendraceae）、水青树科（Tetracentraceae）、连香树科（Cercidiphyllaceae）、折扇叶科（Myrothamnaceae）、悬铃木科（Platanaceae）和金缕梅科（Hamamelidaceae）。这些科之间的谱系关系以及金缕梅科内部的谱系关系是研究被子植物系统和进化的主要问题之一。本研究的目的在于通过对现有的和本研究中产生的分子和发育学资料的分析来对这些谱系问题进行探讨，希望使我们对这些类群的认识提高到了一个新的水平。 本研究测定了包括连香树科、领春木科、悬铃木科、水青树科和金缕梅科内所有亚科代表共计14个种的rbcL基因，matK基因和18S-26S核核糖体（nrDNA）内转录间隔区（ITS）的序列，并结合原有资料对所有低等金缕梅类（除非洲分布的折扇叶科）及其相近类群进行了分支分析。主要结论包括;传统低等金缕梅类是一个复系类群;金缕梅科也是一个复系类群;ITS序列分析支持将连香树科包含在金缕梅科内，但是，rbcL资料和matK资料不支持将连香树科包含在金缕梅科内。本研究支持金缕梅科内6个亚科的分类系统，即将其分为阿丁枫亚科（Altingioideae）、红花荷亚科（Rhodoleioideae）、马蹄荷亚科（Exbucklandioideae）、壳菜果亚科（Mytilarioideae）、双花木亚科（Disanthoideae）和金缕梅亚科（Hamamelidoideae），但壳菜果亚科的山桐柴属（Chunia）应置于马蹄荷亚科内。 本研究首次描述和报道了金缕梅科5个亚科6个种的花的形态发生。通过对本研究发现的新证据和文献中记载的资料分析表明金缕梅科的花发育的特征无明显的共同式样。从本研究提供的花发育资料得出的值得注意的结论是红花荷属在金缕梅科中的亚科地位得到支持。对牛鼻栓（Fortunearia sinensis Rehd. et Wils.）的花发育研究表明其花性分化在发育晚期才出现。 本文首次对显脉红荷（Rhodoleia henryi Tong）的胚胎学进行了研究。尽管红花荷属是金缕梅科内适应鸟传粉的特化类群，但其胚胎学的特征与金缕梅科其它植物并无明显区别。因此，本研究不支持红花荷属独立成科而支持其在金缕梅科内的亚科地位。

药用野生稻复合体（稻属）的分子系统发育与进化研究

稻属（Oryza L.）隶属于禾本科（Gramineae）Ehrhartoideae亚科的稻族（Oryzeae），包括两个栽培种（亚洲栽培稻O. sativa和非洲栽培稻O. glaberrima）和大约20多个野生种，广布于热带亚洲、非洲、大洋洲、中美洲和南美洲。药用野生稻复合体（O. officinalis complex）是稻属中最大、也是最复杂的一个复合体，共包括9个种，含有5种染色体组类型（B、C、BC、CD 和E）。作为栽培稻品质改良的重要基因库，药用野生稻复合体在稻属中具有重要的地位。但是，相似的形态和重叠的地理分布使部分物种的分类和鉴定一直较为困难;种内染色体组构成和倍性的不同更增加了分类鉴定的复杂性。这种情况阻碍了对这些野生稻遗传优势的有效利用。另外，由于物种间断分布和缺乏明确的二倍体亲本等原因，药用野生稻复合体内的异源多倍体起源一直存在争议。本文通过细胞核乙醇脱氢酶基因（Adh）和nrDNA的内转录间隔区（ITS）限制性片段长度多态性（RFLP）分析;叶绿体matK 基因、trnL 内含子和trnL-trnF 基因间隔区、核基因Adh和GPA1以及核糖体DNA ITS片段等序列比较的方法，对药用野生稻复合体中染色体组和物种的鉴定、种间系统发育关系，以及异源多倍体CCDD物种的起源和多倍体ITS的分子进化等进行了研究。主要研究结果如下： 1． 利用核Adh 基因限制性片段长度多态性，检测了来自国际水稻研究所基因库的64份药用野生稻复合体的样品。结果证明，所有O. rhizomatis样品都是含C染色体组的二倍体，所有O. minuta样品都是含BC 染色体组的四倍体。但是，种子库中鉴定为O. officinalis、O. punctata和O. eichingeri的样品中，同时都发现了含C染色体组的二倍体和含BC染色体组的四倍体。四倍体的O. officinalis只在印度分布，而且曾被描述为另一个种O. malampuzhaensis。 四倍体的O. punctata，也被一些学者称为O. schweinfurthiana，被发现和其二倍体一样分布广泛。值得注意的是，有两个曾被作为O. officinalis 四倍体的样品实际上是含有CD染色体组的物种O. latifolia。我们的结果增进了对国际水稻研究所种子库中部分野生稻样品染色体组构成的理解, 纠正了以往对药用野生稻复合体样品的错误鉴定，为今后进一步研究和利用这部分资源提供了种质编目的重要基础。 2． 对稻属中代表不同地理分布区的、含CD染色体组的11个样品（包括77个克隆）的ITS片段进行了测序。基于这些ITS序列的限制性片段长度多态性，提出一个快速而可靠的区分稻属CD 染色体组物种的方法。这个方法的具体步骤是：（1）利用通用引物扩增ITS 片段;（2）利用限制性内切酶FokI和/或DraⅢ消化PCR扩增产物;（3）用1%的琼脂糖胶电泳并根据消化产物的片段长度多态性来区分不同物种。 3． 利用包括两个叶绿体片段（matK和 trnL-trnF）、nrDNA内转录间隔区（ITS）和三个核基因（Adh1、Adh2和GPA1）的同源序列分析，探讨了药用野生稻复合体中二倍体物种和它们所代表的染色体组之间的系统发育关系。独立和合并的基因系统发育树都显示了一致的结果，即C染色体组和B染色体组的亲缘关系要比它们和E染色体组的近。三个含C染色体组的二倍体中，O. officinalis 和O. rhizomatis表现出较近的亲缘关系。值得注意的是，在O. eichingeri种内，尽管基于多基因的数据支持来自斯里兰卡的样品和来自非洲的样品聚成一个分支，但是较低的支持率表明, 两个地区的样品之间存在着较高的遗传分化。 4． 稻属中含CD染色体组的物种特产于拉丁美洲，包括O. alta、 O. grandiglumis 和O. latifolia。由于具有相同的染色体组类型、相似的形态特征和重叠的地理分布，这3个物种间的系统发育关系一直存在争论。另外，因为美洲大陆上没有含C和D染色体组的二倍体物种存在，对这些含CD染色体组物种的可能起源也有不同的假设被提出。使这个问题更具挑战性的是，尽管开展了世界范围的收集，至今仍没有找到含D 染色体组的二倍体物种。在本研究中，代表含C、CD和E染色体组以及含G染色体组的外类群共7个物种，共15份样品的2个叶绿体片段（matK和trnL-trnF）和3个核基因（Adh1，Adh2 和 GPA1）部分片段被测序。基于简约法、距离法和最大似然法的系统发育分析都充分支持含CD染色体组的物种起源于一次杂交事件的推论，并且显示，在物种形成时，含C染色体组的物种（O. officinalis 或O. rhizomatis 而非O. eichingeri）可能承担了母本，而含E染色体组的物种（O. australiensis）则可能承担了父本。另外，CCDD物种间非常一致的系统发育关系表明，非常大的分歧存在于 O. latifolia 和其它两个种（O. alta和O. grandiglumis）之间，这个结果倾向于将后两个种处理为同种或同种下不同分类群。 5． 基于178个克隆序列比较，探讨了ITS在稻属多倍体中的致同进化及其系统学意义。研究发现稻属异源四倍体的ITS存在不同形式的进化方式：首先，非洲BBCC四倍体O. eichingeri和O. punctata 的ITS片段同时保留了双亲拷贝，而且系统发育研究表明，二倍体的O. eichingeri和O. punctata 可能是这两个四倍体的直接祖先;其次，亚洲四倍体O. malampuzhaensis和O. minuta 的ITS仅定向保留母本ITS拷贝;另外，美洲CCDD四倍体的ITS序列发生了双向致同进化，即O. alta和O. grandiglumis的ITS位点一致化成C染色体组类型，而O. latifolia一致化成 D/E 染色体组类型。我们的研究进一步表明在利用ITS片段进行系统发育分析时，特别是涉及异源多倍体时必须慎重。 6． 利用栽培稻的微卫星引物，对含B/C染色体组的6个物种，157个体的SSR位点进行扩增。结果在这些亲缘关系稍远的野生稻中得到7个SSR位点，其中5个位点表现出多态性。比较BB、CC和BBCC物种SSR位点的每位点平均等位基因数A、多态位点百分率P和期望杂合度He ，3项指标发现，四倍体物种的遗传多样性，总体上要高于二倍体物种;二倍体物种内部，O. officnalis的遗传变异最大。另外，以遗传相关性为标准，讨论了B/C染色体组物种间的系统发育关系，同时推测了现存二倍体物种和4个BBCC四倍体物种的遗传关系。

苏铁研究与园林规划

首先回顾了以往苏铁研究，在此基础上，对苏铁的区系地理、苏铁中种皮形态、台湾苏铁的模式产地与海南苏铁的名实问题、从叶绿体DNA matK序列及应用Trn L(UAA)序列资料分析苏铁目系统发育关系等方面进行了深入的研究。 提出了我国西南及印度支那地区和澳大利亚及太平洋岛屿为现代苏铁两大分布中心，前者是苏铁科的演化中心，而后者是次生的分化中心。首次提出钉羊齿（Rhaphidopteris）可能是现代苏铁科近祖的观点。 首次系统研究了27种苏铁属植物中种皮的形态，据此把苏铁属种子分为苏铁型、台湾苏铁型、韦德苏铁型、拳叶苏铁型、斯曼苏铁型及粉种苏铁型等6个类型。对苏铁亚属中种皮纹饰的演化趋势以及其与种间亲缘关系进行了探讨。 对台湾苏铁（Cycas taiwaniana）的模式产地进行了深入的考证，认为其来自福建厦门一带的栽培植株，而非其他学者认为的广东汕头或台湾高雄。海南苏铁（C.hainanensis）应作为台湾苏铁的异名处理。 描述了苏铁科一新亚属奥苏铁亚属Subgen. Media，二新组台湾苏铁组Sect. Taiwanianae及韦德苏铁组Sect. Wadeanae。 分子研究表明托叶铁属（Stangeria）与波温铁属（Bowenia）并不构成一个单系类群;大泽米铁属（Macroz以mia）与非洲铁属（Encephalartos）关系最近，而与鳞叶铁属（Lepidozamia）稍远。 其次，在园林规划方面，对广东省江门市白水带风景区的森林植被进行深入的研究，在此基础上开展了白水带植物景观规划与植被改造规划。在跨越林学与园林两大学科方面做了开拓性的尝试，大气势、大手笔营造森林植物景观。以小班为单位进行规划，植被改造依据植物群落演替原理，模拟南亚热带季风常绿阔叶林，将种类贫乏、结构简单的人工林逐步改造为种类丰富、结构复杂的异龄复层混交林。首次以四季不同种类的花色变化的大范围森林景观来营造南亚热带平常并不明显的季相变化;而专类园的规划注意生态与景观的结合，从植物的生态学特性、群落学特性和景观要求出发进行配置，从而克服一般北京植物园注重品种收集，按分类系统或地理分布来布置而轻视植物生态习性与景观效果的通病。对现有的海金沙(Lygodium Japonicum)等层片或单优群落加以整理形成自然且景观优美的特色层片以及以江门的乡土植物与优秀外引植物营造独木成林、万雀迎宾、层峦叠嶂、孔雀开屏、仙鹤望天等景观的规划思想为国内首创;在国内首先提出规划建设热带雨林、冷色北京植物园、佛教北京植物园、国树国花园等专类园。 在白水带规划的同时，就江门市的人居环境进行了深入的调研。对国内城市生态环境建设普遍存在的一些倾向与问题进行了反思，在总结江门城市环境建设得失的基础上提出了营造可持续发展的人居环境七个方面的战略性意见。

壳斗目的系统发育研究

本文采用来自三个基因组——叶绿体、线粒体和核基因组的六个DNA序列片段，对壳斗目的系统发育进行了重建，主要探讨了壳斗目八个科的科间系统演化关系及不同基因系统树间存在冲突的原因。在此基础上针对具体问题，进一步选用核光敏色素PHYC基因序列对杨梅科在壳斗目中的系统位置进行了更深入的探讨。在对壳斗目开展大量分子系统学研究的同时，还对胡桃科化香树属Platycarya的花形态发生过程进行了详细的扫描电镜观察，讨论了胡桃科单性花的起源和演化问题。主要内容包括： 1 基于四个叶绿体DNA序列的分析 以金缕梅属Hamamelis 和朴树属Celtis 为外类群，对壳斗目几乎所有属的代表共31种植物的trnL-F、matK、rbcL和atpB 序列进行了测定，通过四个序列单独和联合分析，得到如下结果：除rbcL基因树由于信息位点少而分辨率较低之外，壳斗目各科作为单系类群在各种分析中都得到了较强的bootstrap（BS）支持。壳斗目的八个科分为三支：南青冈科Nothofagaceae是基部分支;壳斗科Fagaceae接着分出;核心高等金缕梅类 (core “higher” hamamelids) 聚为一支, 这三支也都得到了强支持。基于不同序列构建的系统发育树，主要区别在于对核心高等金缕梅类六个科，即第三支内部分支关系分辨的不同。在trnL-F树上核心高等金缕梅类又分成两个亚支，一亚支是木麻黄科Casuarinaceae－(桦木科Betulaceae－核果桦科Ticodendraceae)，另一亚支是杨梅科Myricaceae－(胡桃科Juglandaceae－马尾树科Rhoipteleaceae);matK树上，杨梅科则与前一亚支，即木麻黄科和桦木科聚在了一起;atpB树上杨梅科又成了其他所有核心高等金缕梅类的姐妹群;基于四序列联合分析构建的系统发育树的拓扑结构基本上与matK基因树一致。胡桃科和马尾树科的亲缘关系在对不同序列的分析中都得到了较强的支持，但对杨梅科的聚类，支持率都很弱。 2 基于线粒体matR基因序列的分析 除南青冈科作为其他所有壳斗目类群的姐妹群得到强支持外，其余的壳斗目的科间系统发育关系都未得到很好的分辨。最简约树的严格一致树显示杨梅科和木麻黄科聚在一起得到的BS支持不强，另外桦木科和核果桦科，壳斗科和马尾树科分别聚在一起，但得到的支持都未超过50%。 3 基于核核糖体18S rDNA序列的分析 在18S rRNA基因最简约树的严格一致树上，壳斗目被分为两支，一支由壳斗科和南青冈科组成，另一支由核心高等金缕梅类组成，BS支持率均不高。核心高等金缕梅类又分成两个亚支，桦木科，核果桦科和木麻黄科组成的亚支得到了较强的支持，由胡桃科、马尾树科和杨梅科组成的另一亚支得到了更强的BS支持。胡桃科和桦木科的单系性都得到了分辨，虽然BS支持率不高。 4 六个DNA序列的联合分析 通过对六个DNA序列的单独和联合分析，探讨了引起基因树间拓扑结构冲突的可能原因。分别用MP法、NJ法和贝叶斯推论对壳斗目进行了系统发育重建，联合分析提供了分辨率最好、支持率最高的壳斗目谱系关系图： 1、 南青冈科是壳斗目其余类群的姐妹群;2、各科的单系性得到很好支持;3、壳斗科是核心高等金缕梅类的姐妹群;4、核心高等金缕梅类分为两个亚支，一亚支包含桦木科、核果桦科和木麻黄科;另一亚支由胡桃科、马尾树科和杨梅科组成，杨梅科是胡桃科和马尾树科的姐妹群，这一亚支的BS支持率仍然很弱。胡桃科和壳斗科的属间关系未得到很好分辨，多数分支的BS支持率和后验概率值都不高。 5 基于核PHYC基因序列证据对杨梅科系统位置的分析 用壳斗科的栎属Quercus和南青冈科的南青冈属Nothofagus做外类群，核心高等金缕梅类分为两支。桦木科和木麻黄科聚在一起，胡桃科和马尾树科为姐妹群再与杨梅科构成一支，这两个分支都得到很强的BS支持。核基因分析支持六个序列联合分析对核心高等金缕梅类各科间关系的分辨。 6 化香树的花器官发生 在扫描电镜下系统地研究了柔荑花序类植物化香树的雄花、雌花和两性花的发生和发育过程。结果表明：该植物雄花和两性花中小苞片和花被片缺乏;雄蕊轮状发生，成熟阶段雄蕊的不规则排列是由于花托的延伸且和苞片的基部融合后造成的;雌花中存在环状的花被片结构但极度退化，雌花两侧的小苞片和花被片的侧面部分构成小坚果的翅;化香树的两性花向心发生，雄蕊先发生，然后雌蕊发生。胡桃科中的单性花是由两性花退化而来。 本研究的主要发现和结论如下： 1第一次用来自不同基因组的多个DNA 序列探讨了壳斗目八科的系统发育关系;取样包括了所有科的几乎所有的属;所得到的系统树具有较高的分辨率和置信度。主要结论是壳斗目的八个科分为三支，南青冈科是最基部的分支，壳斗科做为核心高等金缕梅类的姐妹群第二个分出，核心高等金缕梅类聚在一起，并进一步分为两个亚支：第一亚支包括桦木科、核果桦科和木麻黄科;另一亚支则由胡桃科、马尾树科和杨梅科组成，杨梅科是胡桃科和马尾树科的姐妹群。 2用核光敏色素PHYC基因较好的解决了杨梅科的系统位置，尽管造成杨梅科在叶绿体基因树和核基因树上具有不同系统位置的原因尚需要进一步探讨。在PHYC基因树上桦木科和木麻黄科聚在一起，胡桃科和马尾树科为姐妹群再与杨梅科构成一支，这两个分支都得到很强的支持。 3首次在壳斗目植物中用扫描电镜观察到了雄花的发生过程和两性花的发生方式，澄清了在化香树属植物中关于雄蕊排列方式、花被式样、以及果翅来源等问题的疑惑或争论。

分子系统学分析平台的建立和应用

分子系统学建立在实验和计算的基础之上。DNA快速测序技术的普及为分子系统学家提供了大量数据，而序列分析技术则是探索数据发现知识的重要工具。在基因组时代，随着大量模式生物完整基因组序列的获得，分子系统学正面临着前所未有的机遇和挑战。一方面，生命之树计划有助于确定新的模式生物和开展相应的基因组计划;另一方面，模式生物的基因组计划有助于阐明它们之间的进化关系和基因组的进化模式。更为重要的是，分子系统学序列分析技术已经发展成为探索与整合基因组数据的强有力工具，从而在生命科学中发挥重要作用。事实上，分子系统学和基因组学的相互渗透正在形成一门崭新的交叉学科——系统发育基因组学。 为了奠定分子系统学研究中信息管理和数据分析工作的坚实基础，我们建立了分子系统发育分析平台。该平台为研究人员提供专业数据库服务和数据分析技术支持，以及相关的网络资源。 分子系统发育分析平台包括了3个专业数据库。第一个是DNA凭证标本数据库。该数据库中的记录包括了7项字段：英文科名、中文科名、物种拉丁名、采集人、采集号、采集地和采集时间。用户可以通过设定单个或多个字段的取值进行检索。截止2004年6月1日，该数据库共包括3491条标本记录。第二个是引物数据库。PCR引物是分子系统学实验的重要条件之一。该数据库中的记录包括3项字段：引物名称、序列内容和退火温度。用户可以通过设定单个或多个字段的取值进行检索。截止2004年6月1日，该数据库共包括170条用于扩增植物细胞核、叶绿体和线粒体基因组DNA序列的引物记录。第三个是生物计算数据库。该数据库为研究人员提供传输和保存序列分析数据和结果文件的服务。 为了确保数据库的安全性和使用性，我们开发了数据库的接口和检索工具，以及系统管理员和用户资格认证程序。通过前者，使用者可以进行数据的上传、下载、管理和检索等操作。而后者则是对不同使用者身份和权限进行设定。管理员的权限高于用户，主要负责本系统的日常维护和管理工作，以及对新增管理员和用户进行资格认证。 分析技术支持旨在帮助用户快速掌握常用的系统发育分析方法，进行有效的数据分析，从复杂的统计学算法和计算机程序中解放出来，将精力集中于计算结果的生物学解释。在该部分中，我们首先简要介绍了常用的分析方法，并且针对分子系统学中的不同问题提供了相应的解决方案。这些问题包括：系统发育重建、替代速率和分歧时间的估计、祖先分布区的重建、性状进化假说的检验、以及密码子水平适应性进化的检测。我们特别强调了似然比检验和贝叶斯推测作为方法论上的重要进展在分子系统学中所发挥的关键作用。本部分还包括大量常用的分子系统学程序或软件包及其快速使用说明和命令模块。下载安装之后，用户即可按照说明使用命令模块进行数据分析。 此外，该平台还提供了一些常用的网络资源地址，如生物信息中心、分子进化和系统发育实验室、专业期刊和相关数据库等。 最后还给出了4个应用实例，即针对特定分子系统学问题的解决方案和初步的分析结果。 第一个例子说明系统发育重建方法的应用。为了确定杨梅科的系统学位置，对6种DNA序列和叶绿体trnL-F区内的间隔性状进行了分析。单个分析表明这6种序列之间在系统学信息上存在显著差异。叶绿体基因组序列的合并分析强烈支持杨梅科和（木麻黄科，（桦木科，核果桦科））的姐妹群关系，而间隔性状的存在能够充分提高其分辨率和支持率。 第二个例子说明如何推测历史生物地理学过程。我们对壳斗目8科25属植物叶绿体基因组的trnL-F、matK、rbcL和atpB的合并序列进行了最大简约分析，得到唯一的最大简约树。基于该系统树和25属植物的地理分布数据，采用扩散－替代分析方法重建了系统树每个节点上的祖先分布区，推测了壳斗目的分布历史。结果表明，壳斗目的历史生物地理学过程由3次替代事件和20次扩散事件组成。其中最重要的替代事件是由于冈瓦纳大陆和劳亚大陆分离所导致的南青冈科及其姐妹群之间的分化。另外，在壳斗科和核心高等金缕梅类中多次发生从欧亚大陆到北美洲、甚至南美洲的平行扩散事件。 第三个例子说明如何估计分歧时间。我们仍然使用扩散－替代分析中所用的最大简约树作为分析的依据，并根据等级制似然比检验确定的最优替代模型对该系统树的支长进行了最大似然优化。似然比检验表明，该系统树不服从分子钟假说。我们以冈瓦纳大陆和劳亚大陆分离的地质事件和5个属的最早化石记录作为标定点，采用罚分似然法在没有分子钟的条件下估计了壳斗目的科间分歧时间。结果表明，绝大多数科间分歧事件都发生在白垩纪。 第四个例子说明如何检测密码子水平的适应性进化。分支间可变选择压力模型的似然比检验表明SARS冠状病毒的S基因在跨种传播过程中发生了正选择。

稻族的分子系统学研究

稻族（Oryzeae）是禾本科中包含多种经济植物的重要类群，现有大约12个属，广布于全球的热带和温带地区。虽然有证据表明稻族是一个单系类群，但稻族的分类处理和属间系统发育关系以及稻族的生物地理学等方面仍存在许多悬而未决的问题。本研究利用了7个DNA片段，包括2个叶绿体基因片段（trnL和matK）、1个线粒体基因片段（nad1）和4个分布在不同染色体上的核基因片段（Adh1、Adh2、GPA1和Waxy）以及形态性状（87个）重建了稻族的系统发育关系，并在此基础上利用6个DNA片段（trnL、matK、nad1、Adh1、Adh2和GPA1）探讨了稻属基部类群的系统发育关系。同时在系统发育研究的基础上，进一步探讨了稻族物种间的分歧时间及生物地理学问题。主要结果如下： 1.稻族的系统学 多基因和形态性状的分析表明：1）稻族是单系类群可分为两个主要分支，相当于传统的两个亚族，第一个亚族（Zizaniinae）包括稻属﹑假稻属和Porteresia，而第二个亚族（Oryzinae）包括其余8个属。两性花是稻族的原始状态，而单性花多次起源，共起源了3次，因此单性花和两性花的区别不宜作为划分亚族的依据;2）Zizania 与分布在南美的单型属Rhynchoryza关系最近;3）一些单型属（Hydrochloa、Porteresia和Prosphytochloa）的确立得不到分子证据的支持，特别是Hydrochloa，形态证据分析的结果明确不支持建立该单型属;4）形态分析结果支持Zizania 4个种聚为一支，而东亚的一个种（Zizania latifolia）为该属基部类群。 2.稻属基部类群的系统发育关系 多基因系统发育分析表明：1）稻属是一个单系类群;2）所有数据都支持稻属10种染色体组类型各自为单系，包括稻属中最后确定的染色体组类型HK（Oryza schlechteri和Porteresia coarctata）;3）F染色体组与HK和HJ中的H染色体组关系最近，类似于E染色体组与异源四倍体的CD染色体组中的D染色体组关系最近; 4）双亲遗传的核基因和母系遗传的线粒体和叶绿体基因的对比分析表明，HJ染色体组两个物种的母系来源是H，而HK染色体组两个物种母系来源是K。三个染色体组相比较，H分化较早，J次之，K分化最晚;5）在系统发育研究基础上，本研究认为稻属4个复合体的划分是可信的，原来未划分到任何一个复合体中的两个种O. brachyantha和O. schlechteri以及新近归入稻属的Oryza coarctata都应归入O. ridleyi复合体。 3.稻族各谱系的分化时间及其生物地理学意义 利用分子钟及其改进方法对稻族各谱系的起源和分化时间进行了估测，并在此基础上探讨了形成各类群地理分布式样的可能原因。结果表明，稻族在始新世（Eocene）晚期（3640万年前）起源于东南亚随后分为两大支。稻属基部类群G染色体组物种与稻属其他物种在1200万年前分开;稻属中A/B/C/E染色体组类群在700万年前开始分化。 在稻族的第一大支Oryzinae亚族中，稻属和假稻属1400万年前分开后，通过远程扩散及随后的漫长的历史过程形成了目前的分布格局。稻属从东南亚起源并逐渐扩散到非洲、澳洲及美洲;而假稻属则从非洲出发扩散到全球的热带和亚热带地区。这两个属的演化历史非常相似，都是全球广泛分布且起源时间较晚。 在稻族第二大支Zizaniinae亚族中，Zizania是稻族中唯一欧亚-北美间断分布属，与其分布在南美的姊妹群Rhynchoryza在2554万年前分开，这一支可能是先从南美扩散到东亚，然后再从东亚扩散到北美;Zizaniopsis和Luziola两属在南美洲和北美洲都有分布，其分歧时间为2180万年前，这两属具有类似的进化历史，即通过上新世（Pliocene）末期隆起的巴拿马陆桥形成现在的分布格局;Chikusichola与Potamophila/Prosphytochloa这一支从中新世（Miocene）早期（2270万年前）开始分化，由于Chikusichola（分布在东南亚）、Potamophila（澳大利亚）和Prosphytochloa（非洲）相距很远，这3个属的扩散必然与跨洋远距离扩散有关。

稻属A基因组的分子系统发育--兼论亚洲栽培稻的起源

稻属Oryza隶属禾本科Poaceae，包括20多个野生种和2个栽培种(亚洲栽培稻O. sativa L和非洲栽培稻O. glaberrima Steud) ，广泛分布于全球热带和亚热带。稻属物种可划分为10个基因组(又称染色体组)类型：A, B, C, BC, CD, E, F, G, HJ 和 HK。栽培稻所属的A基因组是稻属中物种数目最多、地理分布最广的基因组类型，由8个种组成。由于栽培稻属于A基因组，故A基因组物种是栽培稻遗传改良的巨大基因源。数十年来，国际上许多学者对A基因组类群开展了大量涉及形态、细胞、同工酶和分子标记方面的研究，但由于A基因组物种间遗传关系十分接近，形态上差异小且地理分布重叠，使得A基因组物种的系统发育、物种起源和生物地理学等方面存在诸多悬而未决的问题，是稻属中分类和鉴定困难较多的类群。本文利用核基因内含子序列，结合转座子插入分析，重建了A基因组的系统发育，估测了各类群的分化时间;与此同时，基于多克隆测序和基因谱系分析，探讨了O. rufipogon和O. nivara遗传关系以及亚洲栽培稻起源。主要研究结果如下： 1. A基因组的系统发育 在水稻全基因组数据库搜索的基础上，测定了4个单拷贝核基因（Adh1 及3个未注释基因）的内含子序列，构建了稻属A基因组8个种的系统发育关系。基于最大简约法和贝叶斯法的系统发育分析表明：1）澳大利亚的O. meridionalis为A基因组的基部类群;2）亚洲栽培稻两个亚种O. sativa ssp. japonica 和 O. sativa ssp. indica分别和不同的野生类群聚为独立的两个分支，支持japonica 和 indica为多次起源;3）O. rufipogon和O. nivara在系统发育树上完全混在一起，显示出二者间不存在遗传分化;4）非洲一年生野生种O. barthii是非洲栽培稻O. glaberrima的祖先，而非洲多年生野生种O. longistaminata与O. glaberrima/O. barthii.亲缘关系较远;5）分子钟方法估测A基因组类群约在2百万年前（2.0MYA）开始分化，亚洲栽培稻和非洲栽培稻，以及亚洲栽培稻的两个亚种则分别在0.7和 0.4 MYA左右开始分化。此外，通过核基因内含子序列与其它常用片段如ITS，matK等对比分析表明，进化速率相对较快的核基因内含子序列可以有效地用于近缘类群的系统发育研究。 2. Oryza rufipogon 和O. nivara群体遗传研究及亚洲栽培稻起源 对于亚洲野生类群O. rufipogon和O. nivara是合并为一个种还是处理为两个独立的种一直存在争议。在系统发育研究基础上，我们选取4个核基因内含子或5’-UTR区（Waxy, LHS，CatA和1个未注释基因），对采自整个分布区的群体样品进行了多克隆测序，结果表明：1）检测到O. rufipogon和O. nivara均有较高的核苷酸多态性，4个位点上π值和θw值平均分别为0.011和0.014;2）且二者在遗传上没有明显分化，两个类群在4个核基因位点上均检测到大量共享多态（shared polymorphism），未发现固有差异(fixed difference)，表明它们历史上可能属于一个大群体，支持将二者作为种内不同生态型或亚种处理;3）基因谱系树表明亚洲栽培稻的两个亚种indica和japonica分别和不同的O. rufipogon (包括O. nivara)群体聚在一起，进一步从基因谱系角度支持亚洲栽培稻多次起源假说。 3．转座子在群体遗传与系统发育研究中的应用 鉴于目前植物谱系地理学研究中缺乏具有足够信息量的分子标记用于检测种内遗传变异，我们选取3个核基因中的转座子，通过对取自O. rufipogon和O. nivara整个分布区的37份样品的克隆测序，探讨了进化速率快、信息含量丰富的转座子序列在群体遗传上的应用。结果表明：1）无论在物种水平还是群体水平，转座子能检测到比包括内含子在内的其它DNA区域高得多的遗传变异;2）在物种水平上，异交多年生的O. rufipogon和自交一年生的O. nivara多样性均较高，且2个种间相差很小，二者在3个位点上平均核苷酸多样性π值均为0.013，差别主要表现在O. rufipogon杂合位点比例（46.1%）明显高于O. nivara（9.1%），说明交配系统不同并不一定和物种多样性水平相关;3）是否发生转座子序列插入是有价值的系统发育信息，发生在不同染色体上3个基因中的转座子插入进一步证实A基因组基部类群是O. meridionalis;通过叶绿体中3个转座子的插入现象推断了稻族一些四倍体物种，如稻属BC基因组的一些类群的母本来源。