生物模式的元胞自动机模型(I):盘基网柄菌的聚集模式

对生物模式的形成机制的探讨一直是生命科学特别是发育生物学的重要课题.目前已经积累了大量的多学科的研究数据并提出了一些的理论,但生物模式形成的真正机制仍然很不清楚而需更深入的探索.本文试图运用元胞自动机方法建立一个从单细胞及其行为到细胞与细胞、细胞与胞外环境相互作用下生物模式形成的模型.并应用此模型,基于"诱导开关"概念,提出一种新的离散模型来模拟盘基网柄菌(Dictyostelium discoideum)的聚集模式.

基于细胞行为的生物学图式形成的离散模型研究

生物学图式及其形成规律一直是生命科学特别是发育生物学的重要课题；同时也是组织工程中实现体外组织构建的核心科学问题之一。长期以来，对生物图式形成的模型研究的根本不足之处是以数理方法为基础的动力学模型研究和生物学背景的结合不够。因此，本文试图遵循生物图式本身的形成过程，寻求一条与生物学相适配的途径，即以哺乳动物组织发育／活组织工程化构建为目标，以细胞行为为基点，以力学一化学藕合作用为介导，以元胞自动机方法为基础，建立生物学图式形成的一个细胞一环境整体离散模型。应用这一整体离散模型，在不同的控制参数下，对盘基网柄菌的聚集图式和杆菌的生长图式进行了系统的分析，对血管发生（vasculogenesis）的自组装图式进行了初步的新的探索，得到了与实验研究定性上一致的结果。提出了“诱导开关”概念，对盘基网柄菌（Dictyostelium discoideu），杆菌（Bacillus）和血管内皮祖细胞（Endothelial Precursor cells，EPC）三种模式生物，分别以cAMP的信号波前，营养微粒，VEGF的浓度梯度等为诱导开关量。在对盘基网柄菌细胞接收到cAMP后分泌和定向迁移形成的聚集图式的模拟中，系统地考察了影响聚集图式的各种控制参数；一个重要的结果表明细胞初始响应间期对形成的聚集模式有十分显著的影响；引入聚集速度、回转半径、盒质量分布系数等概念对盘基网柄菌的聚集图式进行了一些定量描述的探索。在对杆菌因代谢、增殖、凋亡／衰亡而形成的生长图式的模拟中，系统地定量地分析了在初始营养浓度、营养／代谢物扩散快慢、代谢抑制三者藕合作用下的生长图式；引入定向流动边界，考察了杆菌向营养入口方向的优势生长。在对血管内皮祖细胞经vEGF分子浓度梯度场的诱导进行定向迁移，分化为血管内皮细胞，并自组装形成网状的原初毛细血管丛的模拟中，建立了一个微血管发生自组装图式的新的离散模型，为以后加入力与内皮细胞的相互作用以及血管再生等构建了一个前期模型框架；初步考察了细胞的浓度，细胞分泌vEGF分子的周期，vEGF分子的扩散时间尺度等对血管发生图式的影响因素。

蔷薇类的分子系统学研究—兼论系统发育基因组学的原理和方法(以SET基因超家族为例)

分子系统发育分析的主要任务包括：（1）帮助建立生命之树（tree of life）;（2）追踪基因和基因家族（gene family）的起源和进化, 以获知基因在进化过程中的功能分化和伴随发生的重要分子事件（key molecular events）和形态性状的关键创新（key innovation）。这两个方面在本研究中都有所涉及。对于前者，选用植物线粒体matR基因重建被子植物蔷薇类群的系统发育关系;对于后者，则以SET基因超家族为例，探讨其在真核生物中的进化分类以及与功能多样性的关系。 I 蔷薇类的分子系统学 蔷薇类（rosids）是基于分子数据建立的被子植物的主要分支之一，包含13个目，大约三分之一的被子植物物种。两个主要蔷薇类内部分支是豆类fabids（包含7个目）和锦葵类malvids（包含3个目）。然而，这两个分支内部，以及这两个分支与蔷薇类基部类群，包括牻牛儿苗目（Geraniales）、桃金娘目（Myrtales）和流苏子目（Crossosomatales）之间的关系大多是不清楚的。本研究中，我们选取174个物种来代表72个蔷薇类（rosids）的科，利用两个数据集，即线粒体matR单基因数据集和包括线粒体matR基因、两个质体基因(rbcL、 atpB)和一个核基因(18S rDNA) 的4基因数据集，重建蔷薇类在科以上分类阶元水平的系统发育关系。同时，还对线粒体matR基因的进化特征和用于大尺度系统发育分析的适合度和潜力进行了评价。 线粒体matR单基因数据支持malvids和大多数蔷薇类目的单系性质，然而，豆类（fabids）成员没有形成一个分支，其COM亚支，包括卫矛目（Celastrales）、酢浆草目（Oxalidales）、金虎尾目（Malpighiales）和蒜树科（Huaceae），分辨为锦葵类（malvids）的姐妹群。这个关系在最近根据花结构特征曾被提出过，但从未在之前的分子系统发育分析中得到分辨。4基因数据集支持首先是牻牛儿苗目（Geraniales），接着是桃金娘目（Myrtales）作为蔷薇类（rosids）的最基部的分支;流苏子目（Crossosomatales）是锦葵类（malvids）姐妹群，以及蔷薇类（rosids）的核心部分包括豆类（fabids），锦葵类（malvids）和流苏子目（Crossosomatales）。线粒体matR基因的进化特征分析显示，与两个叶绿体基因(rbcL 和atpB)比较，同义替代速率约是它们的1/4，而非同义替代速率接近于自身的同义替代速率，表明matR 基因具有松弛的选择压力。线粒体matR基因相对慢速的进化使非同源相似（homoplasious）突变减少，提高了系统发育信息的质量，同时，松弛的选择压力使非同义替代数量增加，弥补了慢速进化导致的系统发育信息数量不足的缺陷，这两个方面的结合使线粒体matR基因非常适用于被子植物在科以上水平的系统发育研究。 II SET基因超家族的系统发育基因组学分析 SET基因超家族基因编码含有SET结构域的蛋白，在真核生物中，SET-domain蛋白一般是多结构域（multi-domain）的。SET-domain蛋白具有对组蛋白H3和H4的N末端尾部进行赖氨酸残基甲基化修饰的酶活性;从异染色质形成到基因转录，甲基化的组蛋白广泛影响染色质水平的基因调控。依据SET结构域一级序列的相似性和结构域组织（domain architecture）特征，目前，SET-domain基因超家族被划分为4-7个家族。由于这些划分或者使用动物或者使用植物SET基因，只有少数其它类群的物种加入分析，因此这样的划分可能是不完整的。本研究采用系统发育基 因组学方法（phylogenomic approach），在真核生物范围内广泛取样，期望获得相对完整的SET-domain基因家族的 进化分类方案，在此基础上加深理解SET-domain基因的进化机制和功能多样性。 在提取了17个物种，代表5个真核超群的SET蛋白序列基础上，系统发育分析结合“结构域组织特征”鉴别了9个SET基因家族，其中一个是新的SET基因家族。以前的SET8和Class VI家族,及SMYD和SUV4-20家族分别合并为一个家族。大部分家族在进化过程中发生了2次以上的基因重复事件，通过获得不同的结构域产生具有不同功能的新基因。一个SET基因家族在进化过程中推测发生了从脊椎动物祖先向盘基网柄菌（Dictyostelium discoideum）的水平基因转移。

Detecting positive darwinian selection in brain-expressed genes during human evolution

To understand the genetic basis that underlies the phenotypic divergence between human and non-human primates, we screened a total of 7176 protein-coding genes expressed in the human brain and compared them with the chimpanzee orthologs to identity genes

Sox genes evolution in closely related young tetraploid cyprinid fishes and their diploid relative

Previous study and analysis of cytochrome b suggested that polyploidization event in the genus Tor occurred about 10 Mya ago. In order to understand evolutionary fates of Sox gene in the early stage of genome duplication at the nucleotide level, PCR surveys for Sox genes in three closely related cyprinid fishes T douronensis (2n = 100), T qiaojiensis (2n = ?), T sinensis (2n = 100) and their relative T brevifilis (2n = 50) were performed. Totally, 52 distinct Sox genes were obtained in these four species, representing SoxB, SoxC, and SoxE group. As expected, isoforms of some Sox genes correspond with the ploidy of species, such as two copies of Sox9a exist in tetraploid species. Analysis indicated that duplicated Sox gene pairs caused by polyploidization evolved independently of each other within polyploid species. Results of substitution rate showed nearly equal rate of nonsynonymous substitution of duplicated Sox orthologs among different polyploid species and their diploid relative orthologs, suggesting at the early stage of genome duplicated Sox orthologs are under similar selective constraints in different polyploidy species and their diploid relative at the amino acid level. All PCR fragments of Sox genes obtained in this study are not accompanied by obvious increase in mutations and pseudogene formation which means that they are under strong purifying selection, suggesting that they are functional at the DNA level. Cenealogical analysis revealed that T qiaojiensis was tetraploid, and T douronensis, T qiaojiensis as well as T sinensis had an allotetraploid ancestor. (C) 2009 Elsevier B.V. All rights reserved.

Characterization of C-C chemokine receptor subfamily in teleost fish

Chemokines and their receptors play important roles in nervous and immune systems. Little information, however, exists concerning this gene family in teleost fish. In the present Study, 17 C-C chemokine receptors genes were identified from Danio rerio, 9 from Gasterosteus aculeatus, 10 from Oryzius latipes, 8 from Takifugu rubripes and 5 from Tetraodon nigroviridis. Phylogenetic analysis showed that the orthologs to mammalian CCR6, 7, 8, 9 and CCRL1 receptors were evident in zebrafish, but the clear orthologs to mammalian CCR1, 2, 3, 4, 5 and 10 were not found in zebrafish. The gene structure of zebrafish CCR (zfCCR) was further analyzed. The open reading frame of zfCCR3-1, zfCCR3-3, zfCCR6-1, zfCCR6-2, zfCCR8-2 contain one exon, and two exons were identified for zfCCR2-1, zfCCR2-2, zfCCR4 and zfCCRLI-1, three exons for zfCCR3-2, zfCCR5 and zfCCR7, four exons for zfCCR8-1 and zfCCR9-1. The expression analyses showed that in zebrafish, most C-C chemokine receptor genes Were expressed in fertilized eggs and oocytes, and all the receptor genes were expressed in larval stages. The zfCCR2-2, zfCCR3-1, zfCCR4 and zfCCR6-2 genes were expressed in all normal organs examined, whereas not for zfCCR2-1, zfCCR3-3, zfCCR6-1, zfCCR8-1, zfCCR9-2 and zfCCRL1-2. The expression of zfCCR3-2, zfCCR5, zfCCR7, zfCCR9-1 and zfCCRLI-1 were detected in the majority organs. and zfCCR8-2 and zfCCR8-3 detected only in brain. The differential expression pattern of different paralogues in organs may indicate their difference in function, which requires further investigation. (C) 2008 Elsevier Ltd. All rights reserved.

Sox genes in grass carp (Ctenopharyngodon idella) with their implications for genome duplication and evolution

The Sox gene family is found in a broad range of animal taxa and encodes important gene regulatory proteins involved in a variety of developmental processes. We have obtained clones representing the HMG boxes of twelve Sox genes from grass carp (Ctenopharyngodon idella), one of the four major domestic carps in China. The cloned Sox genes belong to group B1, B2 and C. Our analyses show that whereas the human genome contains a single copy of Sox4, Sox11 and Sox14, each of these genes has two co-orthologs in grass carp, and the duplication of Sox4 and Sox11 occurred before the divergence of grass carp and zebrafish, which support the "fish-specific whole-genome duplication" theory. An estimation for the origin of grass carp based on the molecular clock using Sox1, Sox3 and Sox11 genes as markers indicates that grass carp (subfamily Leuciscinae) and zebrafish (subfamily Danioninae) diverged approximately 60 million years ago. The potential uses of Sox genes as markers in revealing the evolutionary history of grass carp are discussed.

Identification of an Atlantic salmon IFN multigene cluster encoding three IFN subtypes with very different expression properties

A cluster of 11 interferon (IFN) genes were identified in the Atlantic salmon genome linked to the growth hormone I gene. The genes encode three different IFN subtypes; IFNa (two genes), IFNb (four genes) and IFNc (five genes), which show 22-32% amino acid sequence identity. Expression of the fish IFNs were studied in head kidney, leukocytes or To cells after stimulation with the dsRNA poly I:C or the imidazoquinoline S-27609. In mammals, poly I:C induces IFN-beta through the RIG-I/MDA5 or the TLR3 pathway, both of which are dependent on NF-kappa B. In contrast, S-27609 induces mammalian IFN-alpha in plasmacytoid dendritic cells through the TLR7 pathway independent of NF-kappa B. The presence of an NF-kappa B site in their promoters and their strong up-regulation by poly I:C, suggest that salmon IFNa1/IFNa2 are induced through similar pathways as IFN-beta. In contrast, the apparent lack of NF-kappa B motif in the promoter and the strong upregulation by S-27609 in head kidney and leukocytes, suggest that IFNb genes are induced through a pathway similar to mammalian IFN-alpha. IFNc genes showed expression patterns different from both IFNa and IFNb. Taken together, salmon IFNa and IFNb are not orthologs of mammalian IFN-beta and IFN-alpha, respectively, but appear to utilize similar induction pathways. (C) 2008 Elsevier Ltd. All rights reserved.

The expression of AmphiTCTP, a TCTP orthologous gene in amphioxus related to the development of notochord and somites

The translationally controlled tumor protein (TCTP) is highly conserved and has been widely found in eukaryotic organisms. Here, we report the phylogenetic analysis and developmental expression of AmphiTCTP, a TCTP homologous gene in cephalochordate amphioxus. Phylogenetic analysis indicates that the putative protein of AmphiTCTP is close to its vertebrate orthologs. The mRNA of AmphiTCTP is found in fertilized eggs, early cleavage embryo and most of the early developmental stages by in situ hybridization and RT-PCR, but its expression is not detectable from late cleavage stage to mid-gastrula. The expression of AmphiTCTP in zygotes and early cleavage stages shows that AmphiTCTP may be a maternal gene. From the early neurula stage onward, AmphiTCTP transcript is localized in the presumptive notochord, presomitic mesoderm, and nascent somites. However, its expression is gradually down-regulated after the notochord and somites have been formed. The expression pattern of AmphiTCTP thus coincides with the differentiation of the notochord and somites, this suggests that AmphiTCTP may not be a housekeeping gene and may play an important role in mesoderm development. (c) 2007 Elsevier Inc. All rights reserved.

Discovery of genes involved in defense/immunity functions in a haemocytes cDNA library from Fenneropenaeus chinensis by ESTs annotation

Crustacean haemocytes play important roles in the host immune response including recognition, phagocytosis, melanization, cytotoxicity and inter-cellular signal communication. Expressed sequence tags (ESTs) analysis is proved to be an efficient approach not only for gene discovery, but also for gene expression profiles performance. In order to further understand the innate immune system and defense mechanisms of Chinese shrimp at molecular level, complementary DNA library is constructed from the haemocyte tissue of Fenneropenaeus chinensis. A total of 2371 cDNA clones are successfully sequenced and the average sequence length is 460 bp. About 50% are identified as orthologs of known genes from other organisms by BLASTx and BLASTn program. By sequences comparability and analysis, 34 important genes including 177 ESTs are identified that may be involved in defense or immune functions in shrimp based on the known knowledge. These genes are categorized into five categories according to their putative functions in shrimp immune system: 13 genes are different types of antimicrobial peptides (AMP, penaeidin, antilipopolysaccharide factor, etc.), and their proportion is about 3 8%; 11 genes belong to prophenoloxidase system (prophenoloxidase, serine proteinase, serine proteinase inhibitor, etc.), and their proportion is about 32%; five genes have high homology with clotting protein (lectin, transglutaminase, etc), and their proportion is about 15%; three genes may be involved in inter-cell signal communication (peroxinectin, integrin), and their proportion is about 9%; two genes have been identified to be chaperone proteins (Hsc70, thioredoxin peroxidase), and their proportion is about 6%. These EST sequences enrich our understanding of the immune genes of F chinensis and will help farther experimental research into immune factors and improve our knowledge of the immune mechanisms of shrimp. (c) 2007 Elsevier B.V. All rights reserved.