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

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

Phylogenetic relationships among two species of golden monkey and three species of leaf monkey inferred from rDNA variation

Restriction maps of rDNA repeats of five species of Colobinae and three outgroup taxa, Hylobates leucogenys, Macaca mulatta, and Macaca irus, were constructed using 15 restriction endonucleases and cloned 18S and 28S rRNA gene probes. The site variation between Rhinopithecus roxellana and Rhinopithecus bieti is comparable to that between Presbytis francoisi and Presbytis phayrei, implying that R. bieti is a valid species rather than a subspecies of R. roxellana. Phylogenetic analysis on the 47 informative sites supports the case for Rhinopithecus being an independent genus and closely related to Presbytis. Furthermore, branch lengths of the tree seem to support the hypothesis that the leaf monkeys share some ancestral traits as well as some automorphic characters.

A common SNP of MCPH1 is associated with cranial volume variation in Chinese population

Microcephaly (MCPH) genes are informative in understanding the genetics and evolution of human brain volume. MCPH1 and abnormal spindle-like MCPH associated (ASPM) are the two known MCPH causing genes that were suggested undergone recent positive selectio

NONINVASIVE GIANT PANDA PATERNITY EXCLUSION

Giant panda hair samples obtained by noninvasive methods served as a source of DNA for amplification of seven giant panda microsatellite loci utilizing the polymerase chain reaction. Thirteen giant pandas held in Chinese zoos were tested for identification of paternity. Some males listed as sires have been excluded as the biological father of captive-born giant pandas. Because of the death of some potential sires, paternity is still not assigned for some giant pandas, although there is a high likelihood that paternity assignment could be made if postmortem samples are available for genetic analysis. The DNA microsatellite variation assayed by the test we have developed provides a rapid, highly informative, and noninvasive method for paternity identification in giant pandas. (C) 1994 Wiley-Liss, Inc.

Phylogenetic relationships of five pika species (genus Ochotona) based on mitochondrial DNA restriction maps

Restriction site mapping of mitochondrial DNA (mtDNA) with 16 restriction endonucleases was used to examine the phylogenetic relationships of Ochotona cansus, O. huangensis, O. thibetana, O. curzoniae and O. erythrotis. A 1-kb length variation between O. erythrotis of subgenus Pika and other four species of subgenus Ochotona was observed, which may be a useful genetic marker for identifying the two subgenera. The phylogenetic tree constructed using PAUP based on 61 phylogenetically informative sites suggests that O. erythrotis diverged first, followed by O. cansus, while O. curzoniae and O. huangensis are sister taxa related to O. thibetana, The results indicate that both O. cansus and O. huangensis should be treated as independent species. If the base substitution rate of pikas mtDNA was 2% per million years, then the divergence time of the two subgenera, Pika and Ochotana, is about 8.8 Ma ago of late Miocence, middle Bao-dian of Chinese mammalian age, and the divergence of the four species in subgenus Ochotona would have occurred about 2.5 - 4.2 Ma ago, Yushean of Chinese mammalian age. This calculation appears to be substantiated by the fossil record.

Mitochondrial DNA sequence polymorphisms of five ethnic populations from northern China

To study the mitochondrial DNA (mtDNA) polymorphisms in a total of 232 individuals from five ethnic populations (Daur, n=45; Ewenki, n=47; Korean, n=48; Mongolian, n=48; Oroqen, n=44) in northern China, we analyzed the control region sequences and typed for a number of characteristic mutations in coding regions (especially the region 14576-16047), by direct sequencing or restriction-fragment-length-polymorphism (RFLP) analysis. With the exception of 14 individuals belonging to the European-specific haplogroups R2, H, J, and T, the mtDNAs considered could be assigned into the East Asian-specific haplogroups described recently. The polymorphisms in cytochrome b sequence were found to be very informative for defining or supporting the haplogroups status of East Asian mtDNAs in addition to the reported regions 10171-10659 and 14055-14590 in our previous study. The haplogroup distribution frequencies varied in the five ethnic populations, but in general they all harbored a large amount of north-prevalent haplogroups, such as D, G, C, and Z, and thus were in agreement with their ethnohistory of northern origin. The two populations (Ewenki and Oroqen) with small population census also show concordant features in their matrilineal genetic structures, with lower genetic diversities observed.

Phylogeny and zoogeography of the cyprinid genus Epalzeorhynchos Bleeker (Cyprinidae : Ostariophysi)

A total of 62 variable osteological and external characters was found among the five currently recognized species of Epalzeorhynchos. When the genera Crossocheilus and Paracrossocheilus are combined as the outgroup, only 30 of these characters can be polarized. This includes six autapomorphies. The remaining 24 polarized characters form a data matrix which yields a single, 26-step tree with a Consistency Index (CI) of 1 and a Retention Index (RI) of 1, The analysis was also performed on a combined dataset in which the 32 unpolarized characters (characters for which the combined outgroup was dimorphic) were added and coded as missing data (i.e., "?"), Analyzing this data matrix with all multistate characters ordered generates the same single most-parsimonious tree with a length of 63 steps, a CI of 0.98 and a RI of 0.97, When either Crossocheilus or Paracrossocheilus is used as the sole outgroup, the same single most-parsimonious tree is produced although the numbers of informative characters and some of the polarities differ. Evidence is presented to support the following hypotheses: (1) E, kalopterus + E, frenatus + E. bicolor + E. munensis form a monophyletic group; (2) E. frenatus + E, bicolor + E, munensis form a monophyletic group with E, kalopterus as its sister group; this speciation event is congruent with the predictions of vicariant speciation mode I; and (3) E. bicolor and E. munensis are sister groups, again congruent with vicariant speciation mode I, Evidence presented here also supports the zoogeographical hypothesis that the faunas of the Indochinese region and the Greater Sundas are more closely related to each other than either of them is to the lower Salween basin fauna and that the lower Mekong, Chao Phraya, and Mac Khlong basin faunas are more closely related to each other than any of them is to the Greater Sundas, In addition, the monophyly of Epalzeorhynchos is also preliminarily discussed by including either Paracrossocheilus or Crossocheilus in the ingroup. It is demonstrated that E. bicornis clusters with Paracrossocheilus when Paracrossocheilus is included in the ingroup. It seems likely that the taxonomic position of E. bicornis will be resolved as more fishes of the Crossocheilus group are included in future studies.

Variation in mitochondrial DNA and phylogeny of six species of pikas (Ochotona)

Phylogenetic relationships of six species of Ochotona were investigated using mitochondrial DNA (mtDNA) restriction-site analysis. The phylogenetic tree constructed using PAUP was based on 62 phylogenetically informative sites with O. erythrotis designated as an outgroup. Two clades were evident. One contained O. curzoniae, O. huangensis, and O. thibetana. in the second, O. daurica was a sister taxon of O. cansus. The five species appear to come from different maternal lineages. The branching structure of the tree and sequence divergence confirm that huangensis and cansus are distinct species, and that mol-osa is a synonym of O. cansus rather than O. thibetana. Divergence time, estimated from genetic distances, indicates that the ancestors of the two maternal lineages diverged ca. 6.5 x 10(6) years ago. O. curzoniae diverged from O. huangensis, and O. daurica diverged from O. cansus, at about the same time (ca. 3.4 x 10(6) years ago). These data suggest a period of rapid radiation of the genus Ochotona in China, perhaps during the late Pliocene. These calculations correspond roughly to tectonic events and environmental changes in China throughout this period, and appear to be substantiated by the fossil record.

Identification and characterization of a new class of human immunodeficiency virus type 1 recombinants comprised of two circulating recombinant forms, CRF07_BC and CRF08_BC, in China

We identified a new class of human immunodeficiency virus type 1 (HIV-1) recombinants (00CN-HH069 and 00CN-HH086) in which further recombination occurred between two established circulating recombinant forms (CRFs). These two isolates were found among 57 HIV-1 samples from a cohort of injecting drug users in eastern Yunnan Province of China. Informative-site analysis in conjunction with bootscanning plots and exploratory tree analysis revealed that these two strains were closely related mosaics comprised of CRF07_BC and CRF08_BC, which are found in China. The genotype screening based on gag-reverse transcriptase sequences if 57 samples from eastern Yunnan identified 47 CRF08_BC specimens (82.5%), 5 CRF07_BC specimens (8.8%), and 3 additional specimens with the novel recombinant structure. These new "second-generation" recombinants thus constitute a substantial proportion (5 of 57; 8.8%) of HIV-1 strains in this population and may belong to a new but yet-undefined class of CRF. This might be the first example of CRFs recombining with each other, leading to the evolution of second-generation inter-CRF recombinants.

Molecular phylogenetics of the family Cyprinidae (Actinopterygii : Cypriniformes) as evidenced by sequence variation in the first intron of S7 ribosomal protein-coding gene: Further evidence from a nuclear gene of the systematic chaos in the family

The family Cyprinidae is the largest freshwater fish group in the world, including over 200 genera and 2100 species. The phylogenetic relationships of major clades within this family are simply poorly understood, largely because of the overwhelming diversity of the group; however, several investigators have advanced different hypotheses of relationships that pre- and post-date the use of shared-derived characters as advocated through phylogenetic systematics. As expected, most previous investigations used morphological characters. Recently, mitochondrial DNA (mtDNA) sequences and combined morphological and mtDNA investigations have been used to explore and advance our understanding of species relationships and test monophyletic groupings. Limitations of these studies include limited taxon sampling and a strict reliance upon maternally inherited mtDNA variation. The present study is the first endeavor to recover the phylogenetic relationships of the 12 previously recognized monophyletic subfamilies within the Cyprinidae using newly sequenced nuclear DNA (nDNA) for over 50 species representing members of the different previously hypothesized subfamily and family groupings within the Cyprinidae and from other cypriniform families as outgroup taxa. Hypothesized phylogenetic relationships are constructed using maximum parsimony and Basyesian analyses of 1042 sites, of which 971 sites were variable and 790 were phylogenetically informative. Using other appropriate cypriniform taxa of the families Catostomidae (Myxocyprinus asiaticus), Gyrinocheilidae (Gyrinocheilus aymonieri), and Balitoridae (Nemacheilus sp. and Beaufortia kweichotvensis) as outgroups, the Cyprinidae is resolved as a monophyletic group. Within the family the genera Raiamas, Barilius, Danio, and Rasbora, representing many of the tropical cyprinids, represent basal members of the family. All other species can be classified into variably supported and resolved monophyletic lineages, depending upon analysis, that are consistent with or correspond to Barbini and Leuciscini. The Barbini includes taxa traditionally aligned with the subfamily Cyprininae sensu previous morphological revisionary studies by Howes (Barbinae, Labeoninae, Cyprininae and Schizothoracinae). The Leuciscini includes six other subfamilies that are mainly divided into three separate lineages. The relationships among genera and subfamilies are discussed as well as the possible origins of major lineages. (c) 2008 Published by Elsevier Inc.

Phylogeny of the East Asian botiine loaches (Cypriniformes, Botiidae) inferred from mitochondrial cytochrome b gene sequences

The Botiinae have traditionally represented a subfamily of the Cobitidae. At present, the classification and phylogenetic relationships of the Botiinae are controversial. To address systematic and phylogenetic questions concerning this group, we sequenced the complete cytochrome b gene from 34 samples, of which 24 represented 13 species of the East Asian botiine fishes, while the other 10 were non-botiine loach species. For the 1140 bp sequences determined, 494 sites were variable ones, of which 424 were parsimony informative. With Myxocyprinus asiaticus as an outgroup, molecular phylogenetic trees were constructed using the neighbor-joining, maximum parsimony, maximum likelihood and Bayesian methods. All molecular phylogenetic trees revealed that botiine fishes form a monophyletic group and are distantly related to other loaches, suggesting that the Botiinae should be placed in their own family. Within the Botiinae, there are three genera; Botia, Parabotia, and Leptobotia, each genus forming a monophyletic group, with the genus Botia as the most ancestral split. Our molecular results are in agreement with morphological analyses of botiines, suggesting that Botia is the ancestral genus, while Leptobotia and Parabotia were resolved as more derived sister groups.

Amplified fragment length polymorphism analysis to identify the genetic structure of the Gymnocypris przewalskii (Kessler, 1876) population from the Qinghai Basin, China

Amplified fragment length polymorphism (AFLP) was used to analyse the genetic structure of 45 individuals of Gymnocypris przewalskii (Kessler, 1876), an endangered and state-protected rare fish species, from three areas [the Heima (HM), Buha (BH) and Shaliu rivers (SL), all draining into Qinghai Lake]. A total of 563 polymorphic loci were detected. The HM, BH and SL populations have 435, 433 and 391 loci, respectively (Zhu and Wu, 1975), which account for 77.26%, 76.91% and 69.45% of the total number of polymorphic loci of each population, respectively. The Nei indices of genetic diversities (H) of the three populations were calculated to be 0.2869 (HM), 0.2884 (BH) and 0.2663 (SL), respectively. Their Shannon informative indices are 0.4244, 0.4251 and 0.3915, respectively. Research results show that the mean genetic distance between HM and BH is the smallest (0.0511), between BH and SL is the second shortest (0.0608), and between HM and SL is the largest (0.0713), with the mean genetic distance among the three populations being over 0.05. Data mentioned above indicate that the three populations have a certain genetic differentiation. The total genetic diversity (H-t = 0.3045) and the mean value of genetic diversity within the population (H-s = 0.2786) indicate that the variations have mainly come from within the population.

Sequence variations of the S7 ribosomal protein gene in primitive cyprinid fishes: Implication on phylogenetic analysis

Cyprinidae is the largest fish family in the world and contains about 210 genera and 2010 species. Appropriate DNA markers must be selected for the phylogenetic analyses of Cyprinidae. In present study, the 1st intron of the S7 ribosomal protein (r-protein) gene is first used to examine the relationships among cyprinid fishes. The length of the 1st intron obtained by PCR amplification ranges from 655 to 859 by in the 16 cyprinid species investigated, and is 602 by in Myxocyprinus asiaticus. Out of the alignment of 925 nucleotide sites obtained, the parsimony informative sites are 499 and occupy 54% of the total sites. The results indicate that the 1st intron sequences of the S7 r-protein gene in cyprinids are rich in informative sites and vary remarkably in sequence divergence from 2.3% between close species to 66.6% between distant species. The bootstrap values of the interior nodes in the NJ (neighbor-joining) and MP (most-parsimony) trees based on the present S7 r-protein gene data are higher than those based on cytochrome b and the d-loop region respectively. Therefore, the 1st intron sequences of the S7 r-protein gene in cyprinids are sensitive enough for phylogenetic analyses, and the 1st intron is an appropriate genetic marker for the phylogenetic reconstruction of the taxa in different cyprinid subfamilies. However, attempts to discuss whether the present S7 r-protein gene data can be applied to the phylogeny of the taxa at the level of the family or the higher categories in Cypriniformes need further studies.

Retroposon analysis of major cetacean lineages: The monophyly of toothed whales and the paraphyly of river dolphins

SINE (short interspersed element) insertion analysis elucidates contentious aspects in the phylogeny of toothed whales and dolphins (Odontoceti), especially river dolphins. Here, we characterize 25 informative SINEs inserted into unique genomic loci during evolution of odontocetes to construct a cladogram. and determine a total of 2.8 kb per taxon of the flanking sequences of these SINE loci to estimate divergence times among lineages. We demonstrate that: (i) Odontocetes are monophyletic; (ii) Ganges River dolphins, beaked whales, and ocean dolphins diverged (in this order) after sperm whales; (iii) three other river dolphin taxa, namely the Amazon, La Plata, and Yangtze river dolphins, form a monophyletic group with Yangtze River dolphins being the most basal; and (iv) the rapid radiation of extant cetacean lineages occurred some 28-33 million years B.P., in strong accord with the fossil record. The combination of SINE and flanking sequence analysis suggests a topology and set of divergence times for odontocete relationships, offering alternative explanations for several long-standing problems in cetacean evolution.

Structures of an asymmetrically coupled double-well superlattice by double-crystal X-ray diffraction

An asymmetrically coupled (GaAs/AlAs/GaAs/AlAs)/GaAs (001) double-well supperlattice is studied by HRDCD (high resolution double-crystal X-ray diffractometry). The intensity of satellite peaks is modulated by wave packet of different sublayers. In the course of simulation, the satellite peaks in the vicinity of the node points of wave packet are very informative for precise determination of sublayer thickness and for improving accuracy.