松属植物rRNA基因的变异模式及其进化生物学意义

被子植物的rRNA基因已经得到深入研究。二倍体被子植物一般拥有1-4对18S-5.8S-26S rDNA位点和1-2对5S rDNA位点。作为特殊的多基因家族成员，rDNA会受均一化力 (homogenizing forces) 的作用，通过基因转换、不等交换等机制，形成基因的致同进化 (concerted evolution)。长期以来，我们一直认为动植物rDNA致同进化水平很高，各种拷贝的序列几乎完全一致，因此可以直接应用PCR测序的方法进行分子系统学研究。但是在裸子植物中由于研究资料的匮乏，使我们对裸子植物rDNA的变异模式了解甚少。松属植物作为裸子植物的最大类群，它的rDNA变异和进化有何特点、与被子植物是否相同，是这个重要类群的进化研究中目前尚未解决的问题。本文的研究内容从三个方面进行： (1)rDNA的染色体定位 目前，松属的18S-5.8S-26S rDNA的染色体定位研究只包括5种植物，其中的3种同时涉及到5S rDNA定位。这些研究结果表明，不同种存在相异的rDNA位点数目，甚至不同的个体的rDNA位点均有变化。其共同点是，18S-5.8S-26S rDNA位点数平均较被子植物多，5S rDNA除Pinus radiata外，在其它种里则与被子植物相似。这种现象是松属或裸子植物的共同特征，亦或是特例呢？有限的研究限制了对裸子植物rDNA的了解。本研究的目的之一就是研究松属植物rDNA的染色体空间分布特征，希望借此了解松属植物间的关系，比较裸子植物和被子植物rDNA在染色体组水平的差异。 (2)5S rDNA的分子进化 5S rDNA的序列水平的进化研究在松属中尚属空白。5S rDNA在染色体数目上没有显示裸子植物与被子植物的差异，是否意味着松属乃至裸子植物的5S rDNA也同被子植物一样——致同进化完全，序列高度一致呢？利用克隆测序方法对松属植物5S rDNA的研究无疑是有开创性的工作，可以探讨裸子植物的5S rDNA的进化机制和种间关系。 (3)杂种基因组研究 杂交物种的起源演化是当前生物学研究的热点，通过杂种基因组的研究，可以了解杂种的的基因组构成，组织方式和进化历史，探讨杂交事件对成种过程的影响及意义。这项研究涉及到高山松、云南松和油松。之所以采用这三种植物，因为等位酶、cpDNA和mtDNA证据证明高山松为油松和云南松的自然杂交种。但这些证据不足以反映杂种核基因组的重组特征和构成及其进化规律。我们利用rDNA-FISH、5S rDNA和基因组原位杂交分析三种松树间的基因组关系，为揭示高山松的进化机制和历史提供新的依据。 本项研究得到以下结果： 一. rDNA荧光原位杂交 (FISH) 通过对华山松和白皮松两种单维管束亚属植物及油松、云南松、高山松、马尾松和南亚松等五种双维管束亚属植物的18S rDNA与5S rDNA的荧光原位杂交，结果表明： ⑴ 裸子植物的18S rDNA位点数目明显多于二倍体被子植物。其中主要位点数目，油松有7对，高山松5对，云南松8对，马尾松10对，南亚松6对，白皮松3对，华山松10对，平均在7对;另外，部分松树还存在弱位点。无论强弱位点都有部分存在于染色体的着丝粒区，除了赤松 (Pinus densiflora)，在其它松科植物中并没有发现这种现象。究竟是基因转移的结果或该位点是18S rDNA的原始起源位置还有待确证。 ⑵ 5S rDNA位点相对变异较小，与被子植物相当。除了华山松5S rDNA有4对位点，马尾松只有1对位点外，其它松树的5S rDNA位点数目均为2对，并且在双维管束亚属植物中有一对属于弱位点。 ⑶ 两种rDNA存在不同连锁模式。双维管束亚属植物中，5S与18S rDNA连锁在同一染色体的同一臂或两条臂上。在同一染色体臂时，18S rDNA在臂的远端。单维管束亚属植物的5S与18S rDNA或连锁于同一染色体的同一臂上，或分别处于不同染色体。前一情况，5S rDNA位于臂的远端。据此可以说明两个亚属的rDNA结构在染色体组水平的很大分化。 ⑷ 松属植物的关系及高山松核型特征。由于5S与18S rDNA连锁关系的不同，可以将单维管束亚属和双维管束亚属分开。各亚属的不同物种可以依据杂交位点的多少、位置、信号强弱构成的核型图加以区分，并且构成一定的系统关系。杂交起源的高山松在染色体组上，表现出对油松和云南松两亲本不同染色体特征的分别继承与重组，并产生独有的特征。其II同源染色体之一18S rDNA位点的缺失，可能是染色体重组的痕迹。 二. 5S rDNA的序列变异与分子进化 利用分子克隆和DNA测序分析了油松、云南松、马尾松、白皮松和不同遗传背景的高山松居群的5S rRNA基因序列变异及基因进化规律，得到以下主要结果： ⑴ 5S rDNA的结构特征。双维管束亚属植物长度在658-728 bp，白皮松则为499-521 bp。长度差异体现在基因间隔区，而基因区极端保守，基本为120 bp。基因转录区内部存在着转录控制区，决定了5S rRNA的转录起始与转录效率。5S rRNA基因能够折叠成正常的二级结构，其中，相对于干区来说，环区要保守，但环E却表现出异乎寻常的变异，转换/颠换比值高达7.1，这种突变可能是假基因的产物。基因间隔区存在一定的保守单元，其中一些与转录的起始和终止调控相关，有些是裸子植物未知功能的特异保守区。 ⑵ 松属植物5S rDNA存在着基因组内与种间的异质性。基因组内的各个克隆中有超过80%的特异的，彼此不相同。整个5S rDNA分化距离为0.042 - 0.051，其中，间隔区的分化比基因区高，其速度约是基因区的3-7倍。比较种间5S rDNA序列发现：在122个克隆中，基因区只有50个特异的序列。基因组间的序列变异度与基因组内 (个体内) 没有明显差别。白皮松的间隔区与双维管束亚属松树的5S rDNA间隔区差异极大，几乎不能排序，而四种双维管束亚属植物的5S rDNA间隔区种间种内差异不大。 ⑶ 松属植物5S rDNA进化。PAUP分析建立的5S rRNA基因树显示，5S rRNA基因在基因组内是多系的 (polyphyletic)，表明成种事件以前，祖先种就已经存在序列的分化。观测到的5S rRNA基因序列变异状况，并非完全是致同进化或独立进化的单一因素造成的，而是二者的相互作用的结果。致同进化确实存在，只是速度较慢而已。 ⑷ 高山松5S rDNA 组成。高山松拥有最高的基因组内的序列多样性，高山松的5S rDNA拷贝既有亲本类型，又有重组类型，并且不同地理及遗传来源的高山松显示一定的分化趋势，有更多的拷贝来自母系亲本。 三. 基因组原位杂交 以油松和云南松总DNA作为探针，相互进行基因组原位杂交，结果显示云南松和油松的染色体组可以完全被对方探针标记，在现有基因组原位杂交的分辨率下不能将两个基因组区分开。说明云南松和油松基因组之间存在高比例的同源序列，两种松树的基因组组成十分相似。利用油松和云南松总DNA作为探针，对高山松的染色体组进行双探针基因组原位杂交。结果表明，高山松全部基因组都能与两亲本探针完全杂交，说明三者间有着异乎寻常的亲缘关系。但在PH失调影响下，高山松只有部分基因组被杂交，并且两种探针的杂交信号有轻微差异。这可能是高度重复序列优先杂交的结果。这些情况表明，高山松虽然在基因组构成上与两个亲本基本一致，但基因在染色体组的空间排布上是存在差异的，这一点可以从rDNA-FISH中证明。

松属单维管束亚属植物rDNA的染色体定位研究

松属植物的基因组十分庞大（大于20000Mbp），其中约90%是由重复序列组成的，我们对其结构和组成仍知之甚少。松属在系统分类上分为两个亚属:单维管束亚属和双维管束亚属。基因组大小研究发现单维管束亚属植物的基因组更大。rDNA作为一类有功能的多基因家族重复序列，其自身特性决定了它在基因组研究中的重要性。FISH技术为rDNA在染色体上物理定位提供了有力的工具。尽管现在对松属rDNA FISH已有不少报道，但主要集中在双维管束亚属，对单维管束亚属的研究几乎是空白。本研究选择5个单维管束亚属松属植物P. bungeana, P. koraiensis, P. armandii, P. wallichiana, P. strobus，进行rDNA FISH研究。旨在弄清18S-25S rDNA和5S rDNA在单维管束亚属植物染色体上的位点数目和分布模式。结合前人对松属双维管束亚属植物的工作，对单、双维管束亚属植物之间rDNA FISH结果进行比较，从而可以从整体上认识松属植物的18S-25S rDNA和5S rDNA在染色体上的分布式样。在此基础上进一步探讨18S-25S rDNA和5S rDNA这些重复序列在松属植物基因组结构和组成中的地位和作用。本研究主要结果如下： 1.rDNA FISH在松属染色体核型分析中的作用 本研究中5种松属单维管束亚属植物染色体数目均为2n=24，除最短一条染色体为亚中部着丝粒染色体外，其余11条均为中部着丝粒染色体，长度和臂比也十分接近，同源染色体的不容易鉴定，很难排出精确的核型。在我们的研究结果中，5个松属植物中，除了白皮松外，18S-25S rDNA和5S rDNA分布在12对染色体中的10对染色体上，这些位点可作为染色体标记，大大提高了同源染色体鉴定的准确度，但是染色体之间排序问题依然没有很好地解决。核型比较认为种间是否存在部分同源染色体关系也不是十分明确，仅Ⅺ号和Ⅻ号染色体有这种关系，这主要由于Ⅺ号和Ⅻ号染色体容易准确地鉴别出来。核型分析的精确仍有待增加标记来提高。 2．rDNA位点数目在松属两个亚属间的比较及其与基因组大小的关系 松属植物18S-25S rDNA位点通常为5-10个，5S rDNA位点为1-4个。其中单维管束亚属18S-25S rDNA位点通常为9-10个（除白皮松为4个外），5S rDNA位点为2- 4个;双维管束亚属为18S-25S rDNA位点通常为5-10个，5S rDNA位点通常为1-2个。而二倍体被子植物18S-25S rDNA位点通常为1-5个5S rDNA位点为1-3个。暗示18S-25S rDNA和5S rDNA位点数目多少和基因组大小还是有一定的相关性。因为松属植物的基因组比典型的二倍体被子植物大得多，单维管束亚属植物的的C-值又普遍比双维管束亚属植物的高。白皮松虽有些例外，18S-25S rDNA位点数目少，但信号强度大得多，代表拷贝数高，因此其基因组大小可以从rDNA拷贝数上得到解释。 3．18S-25S rDNA和5S rDNA位点在松属两个亚属之间的分布模式比较 18S-25S rDNA和5S rDNA位点在松属两个亚属染色体上的分布方式有明显不同，每个亚属均有两种分布形式，并形成各自稳定的分布模式。在单维管束植物中，18S-25S rDNA和5S rDNA位点或相邻分布于同一染色体同一臂上，5S rDNA位于臂的远端;或两位点分布于不同的染色体。而在双维管束植物中18S-25S rDNA和5S rDNA或相邻分布于同一染色体同一臂上，18S-25S rDNA在臂的远端;或两位点分布于同一染色体两条臂上。在两个亚属中，当18S-25S rDNA和5S rDNA位点位于同一条染色体臂上时，相对位置正好相反。这完全不同的rDNA分布模式的形成，可能与松属这两个亚属植物的物种形成和分化过程中染色体发生倒位或易位有关，暗示这两个亚属的基因组结构存在分化。但这各自的分布模式是否可以作为判断亚属的特征依据仍有待加大样本量证实。 4．rDNA 位点分布及变异具有系统学意义 rDNA FISH 结果符合分类中亲缘关系越近，分布模式越相似的原则，因而认为rDNA 位点在染色体上的分布模式，具有系统学意义。基于已知的松属植物rDNA FISH结果构建的系统关系，符合传统分类系统中对亚组划分。rDNA FISH结果与分子系统学的研究结果相比较认为，松属单维管亚属5种松中，以乔松和北美乔松关系最近，与同一个亚组的华山松稍远，与另一个亚组的红松更远。而白皮松作为一个特有的孑遗类群，系统位置比较特殊，分子系统学研究认为其处于基部的位置，本研究表明其rDNA位点有明显的特点：位点数目少，但信号强，反映了拷贝数多。那是否它就代表了祖先类群的位点分布模式，需要更多的基部类群的rDNA FISH结果支持。

Gene structure and transcription of IRF-1 and IRF-7 in the mandarin fish Siniperca chuatsi

The genes of IRF-1 and IRF-7 have been cloned from the mandarin fish (Siniperca chuatsi). The IRF-1 gene has 4919 nucleotides (nt) and contains 10exons and 9introns, with an open reading frame (ORF) of 903 ntencoding301 aa. The IRF-7 gene has 6057 nt and also contains 10exons and 9introns, with an ORF of 1308 nt encoding 436 aa. The IRF-1 and IRF-7 genes have only one copy each in the genome. The transcription of IRF-1 and IRF-7 in different organs was analyzed by real-time PCR, and both molecules were constitutively expressed. The IRF-I and IRF-7 mRNAs were abundant in gill, spleen, kidney and pronephros. The temporal transcriptional changes for IRF-1, IRF-7 and Mx were investigated within 48 h after poly I: C stimulation in liver, gill, spleen and pronephros. An increased transcription was detected for IRF-1 and IRF-7 12 h post-stimulation, being earlier than the transcription of Mx protein; however, IRF-1 and IRF-7 transcription decreased while the Mx protein was stable at 48 h post-stimulation. (c) 2007 Published by Elsevier B.V.

Complete nucleotide sequence of the S10 genome segment of grass carp reovirus (GCRV)

Hemorrhagic disease, caused by the grass carp reovirus (GCRV), is one of the major diseases of grass carp in China. Little is known about the structure and function of the gene segments of this reovirus. The S10 genome segment of GCRV was cloned and the complete nucleotide sequence is reported here. The S10 is 909 nucleotides long and contains a large open reading frame (ORF) encoding a protein of 276 amino acids with a deduced molecular weight of approximately 29.7 kDa. Comparisons of the deduced amino acid sequence of GCRV S10 with those of other reoviruses revealed no significant homologies. However, GCRV S10 shared a putative zinc-finger sequence and a similar distribution of hydrophilic motifs with the outer capsid proteins encoded by Coho salmon aquareovirus (SCSV) S10, striped bass reovirus (SBRV) S10, and mammalian reovirus (MRV) S4. It was predicted that this segment gene encodes an outer capsid protein.

Genome-wide survey of putative serine/threonine protein kinases in cyanobacteria

Background: Serine/threonine kinases (STKs) have been found in an increasing number of prokaryotes, showing important roles in signal transduction that supplement the well known role of two-component system. Cyanobacteria are photoautotrophic prokaryotes able to grow in a wide range of ecological environments, and their signal transduction systems are important in adaptation to the environment. Sequence information from several cyanobacterial genomes offers a unique opportunity to conduct a comprehensive comparative analysis of this kinase family. In this study, we extracted information regarding Ser/Thr kinases from 21 species of sequenced cyanobacteria and investigated their diversity, conservation, domain structure, and evolution. Results: 286 putative STK homologues were identified. STKs are absent in four Prochlorococcus strains and one marine Synechococcus strain and abundant in filamentous nitrogen-fixing cyanobacteria. Motifs and invariant amino acids typical in eukaryotic STKs were conserved well in these proteins, and six more cyanobacteria- or bacteria-specific conserved residues were found. These STK proteins were classified into three major families according to their domain structures. Fourteen types and a total of 131 additional domains were identified, some of which are reported to participate in the recognition of signals or substrates. Cyanobacterial STKs show rather complicated phylogenetic relationships that correspond poorly with phylogenies based on 16S rRNA and those based on additional domains. Conclusion: The number of STK genes in different cyanobacteria is the result of the genome size, ecophysiology, and physiological properties of the organism. Similar conserved motifs and amino acids indicate that cyanobacterial STKs make use of a similar catalytic mechanism as eukaryotic STKs. Gene gain-and-loss is significant during STK evolution, along with domain shuffling and insertion. This study has established an overall framework of sequence-structure-function interactions for the STK gene family, which may facilitate further studies of the role of STKs in various organisms.

Genome-wide survey of putative serine/threonine protein kinases in cyanobacteria

Background: Serine/threonine kinases (STKs) have been found in an increasing number of prokaryotes, showing important roles in signal transduction that supplement the well known role of two-component system. Cyanobacteria are photoautotrophic prokaryotes able to grow in a wide range of ecological environments, and their signal transduction systems are important in adaptation to the environment. Sequence information from several cyanobacterial genomes offers a unique opportunity to conduct a comprehensive comparative analysis of this kinase family. In this study, we extracted information regarding Ser/Thr kinases from 21 species of sequenced cyanobacteria and investigated their diversity, conservation, domain structure, and evolution. Results: 286 putative STK homologues were identified. STKs are absent in four Prochlorococcus strains and one marine Synechococcus strain and abundant in filamentous nitrogen-fixing cyanobacteria. Motifs and invariant amino acids typical in eukaryotic STKs were conserved well in these proteins, and six more cyanobacteria- or bacteria-specific conserved residues were found. These STK proteins were classified into three major families according to their domain structures. Fourteen types and a total of 131 additional domains were identified, some of which are reported to participate in the recognition of signals or substrates. Cyanobacterial STKs show rather complicated phylogenetic relationships that correspond poorly with phylogenies based on 16S rRNA and those based on additional domains. Conclusion: The number of STK genes in different cyanobacteria is the result of the genome size, ecophysiology, and physiological properties of the organism. Similar conserved motifs and amino acids indicate that cyanobacterial STKs make use of a similar catalytic mechanism as eukaryotic STKs. Gene gain-and-loss is significant during STK evolution, along with domain shuffling and insertion. This study has established an overall framework of sequence-structure-function interactions for the STK gene family, which may facilitate further studies of the role of STKs in various organisms.

Genome-wide analysis of restriction-modification system in unicellular and filamentous cyanobacteria

Cyanobacteria are an ancient group of gram-negative bacteria with strong genome size variation ranging from 1.6 to 9.1 Mb. Here, we first retrieved all the putative restriction-modification (RM) genes in the draft genome of Spirulina and then performed a range of comparative and bioinformatic analyses on RM genes from unicellular and filamentous cyanobacterial genomes. We have identified 6 gene clusters containing putative Type I RMs and 11 putative Type II RMs or the solitary methyltransferases (MTases). RT-PCR analysis reveals that 6 of 18 MTases are not expressed in Spirulina, whereas one hsdM gene, with a mutated cognate hsdS, was detected to be expressed. Our results indicate that the number of RM genes in filamentous cyanobacteria is significantly higher than in unicellular species, and this expansion of RM systems in filamentous cyanobacteria may be related to their wide range of ecological tolerance. Furthermore, a coevolutionary pattern is found between hsdM and hsdR, with a large number of site pairs positively or negatively correlated, indicating the functional importance of these pairing interactions between their tertiary structures. No evidence for positive selection is found for the majority of RMs, e. g., hsdM, hsdS, hsdR, and Type II restriction endonuclease gene families, while a group of MTases exhibit a remarkable signature of adaptive evolution. Sites and genes identified here to have been under positive selection would provide targets for further research on their structural and functional evaluations.

The complete mitochondrial genome of the large yellow croaker, Larimichthys crocea (Perciformes, Sciaenidae): Unusual features of its control region and the phylogenetic position of the Sciaenidae

To understand the systematic status of Larimichthys crocea in the Percoidei, we determined the complete mitochondrial (mt) genome sequence using 454 sequencing-by-synthesis technology. The complete mt genome is 16,466 bp in length including the typical structure of 22 tRNAs, 2 rRNAs, 13 protein-coding genes and the noncoding control region (CR). Further sequencing for the complete CR was performed using the primers Cyt b-F and 12S-R on six L crocea individuals and two L polyactis individuals. Interestingly, all seven CR sequences from L crocea were identical while the three sequences from L polyactis were distinct (including one from GenBank). Although the conserved blocks such as TAS and CSB-1, -2, and -3 are readily identifiable in the control regions of the two species, the typical central conserved blocks CSB-D, -E, and -F could not be detected, while they are found in Cynoscion acoupa of Sciaenidae and other Percoidei species. Phylogenetic analysis shows that L crocea is a relatively recently emerged species in Sciaenidae and this family is closely related to family Pomacanthidae within the Percoidei. L crocea, as the first species of Sciaenidae with complete mitochondrial genome available, will provide important information on the molecular evolution of the group. Moreover, the genus-specific pair of primers designed in this study for amplifying the complete mt control region will be very useful in studies on the population genetics and conservation biology of Larimichthys. (c) 2008 Elsevier B.V. All rights reserved.

The complete mitochondrial genome sequence of the cutlassfish Trichiurus japonicus (Perciformes: Trichiuridae) : Genome characterization and phylogenetic considerations

Mitochondrial genome sequence and structure analysis has become a powerful tool for studying molecular evolution and phylogenetic relationships. To understand the systematic status of Trichiurus japonicus in suborder Scombroidei, we determined the complete mitochondrial genome (mitogenome) sequence using the long-polymerase chain reaction (long-PCR) and shotgun sequencing method. The entire mitogenome is 16,796 by in length and has three unusual features, including (1) the absence of tRNA(Pro) gene, (2) the possibly nonfunctional light-strand replication origin (O-L) showing a shorter loop in secondary structure and no conserved motif (5'-GCCGG-3'), (3) two sets of the tandem repeats at the 5' and 3' ends of the control region. The three features seem common for Trichiurus mitogenomes, as we have confirmed them in other three T. japonicus individuals and in T nanhaiensis. Phylogenetic analysis does not support the monophyly of Trichiuridae, which is against the morphological result. T. japonicus is most closely related to those species of family Scombridae; they in turn have a sister relationship with Perciformes members including suborders Acanthuroidei, Caproidei, Notothenioidei, Zoarcoidei, Trachinoidei, and some species of Labroidei, based on the current dataset of complete mitogenome. T japonicus together with T. brevis, T lepturus and Aphanopus carbo form a clade distinct from Lepidopus caudatus in terms of the complete Cyt b sequences. T. japonicus mitogenome, as the first discovered complete mitogenome of Trichiuridae, should provide important information on both genomics and phylogenetics of Trichiuridae. (C) 2009 Elsevier B.V. All rights reserved.

The complete mitochondrial genome of the mantid shrimp Oratosquilla oratoria (Crustacea: Malacostraca: Stomatopoda): Novel non-coding regions features and phylogenetic implications of the Stomatopoda

The complete mitochondrial (mt) genome sequence of Oratosquilla oratoria (Crustacea: Malacostraca: Stomatopoda) was determined; a circular molecule of 15,783 bp in length. The gene content and arrangement are consistent with the pancrustacean ground pattern. The mt control region of O. oratoria is characterized by no GA-block near the 3' end and different position of [TA(A)]n-blocks compared with other reported Stomatopoda species. The sequence of the second hairpin structure is relative conserved which suggests this region may be a synapomorphic character for the Stomatopoda. In addition, a relative large intergenic spacer (101 bp) with higher A + T content than that in control region was identified between the tRNA(Glu) and tRNA(Phe) genes. Phylogenetic analyses based on the current dataset of complete mt genomes strongly support the Stomatopoda is closely related to Euphausiacea. They in turn cluster with Penaeoidea and Caridea clades while other decapods form a separate group, which rejects the monophyly of Decapoda. This challenges the suitability of Stomatopoda as an outgroup of Decapoda in phylogenetic analyses. The basal position of Stomatopoda within Eumalacostraca according to the morphological characters is also questioned. (C) 2010 Elsevier Inc. All rights reserved.