借助斑马鱼EST数据库从鲤鱼微卫星序列中寻找蛋白编码基因(英文）

应用 in silico的方法,利用Blastn 和Blastx 搜索引擎,将鲤鱼微卫星序列与GenBank数据库进行同源序列比对。利用Blastn,将侧翼序列长度>50 bp的875个鲤鱼微卫星序列与斑马鱼的EST数据库首先进行比对,结果找到了121个同源序列。随后采用Blastx搜索蛋白质数据库,有94个微卫星位点存在同源蛋白。除了33个假定和3个未知蛋白外,剩余的58个微卫星位点被成功地进行了功能注释,而且其中的7个位点已经定位在了鲤鱼连锁图谱上。另外,通过PCR-SSCP的方法,将两个与鲤鱼微卫星侧翼序列相匹配的斑马鱼EST序列开发成鲤鱼的STS标记,并将其中的一个标记HLJZe33定位到鲤鱼连锁图谱上。以上研究结果表明,通过比较基因组研究,模式生物斑马鱼的很多遗传和基因组资源都可以被利用到鲤鱼的基因组研究中。     

鲤EST标记与耐低温性状的相关性分析及定位

利用GLM模型对12个表达序列标签(expressed sequence tag,EST)标记的基因型与鲤耐低温性状进行相关性分析,然后使用OneMap软件将这些EST标记进行连锁定位研究,并通过Blast x搜索引擎对这些候选EST进行注释。结果显示,EST标记CC009(P<0.05)和CC115(P<0.01)与鲤鱼耐低温性状显著相关;12个EST标记中有8个标记分别连锁定位到6个连锁群中,其中与耐低温相关的CC009和CC115分别定位到鲤连锁图谱的第38号连锁群和第2号连锁群;蛋白质数据库同源性比对发现,CC009与斑马鱼(Danio rerio)尿嘧啶激酶1(uridine-cytidine kinaseI)的同源性高达94%;而CC115为原绿球藻(Prochlorococcus marinus str)的假定糖基转移酶(putative glycosyl transferase)的同源性为56%。     

镜鲤体重的QTL定位

本研究利用217个微卫星标记和336个SNPs标记对德国镜鲤F2代68个个体基因组DNA进行基因型检测.其中507个标记共组成62个连锁群,覆盖基因组总长度为2 805.85 cM,标记间平均距离为6.31 cM;利用软件MapQTL 4.0采用区间作图法对体重性状进行QTL定位分析.研究结果共检测到14个与体重性状有关的QTLs,分布于9个连锁群.其中BW-5-1有最大的LOD值,为4.46;BW-1-1的LOD值最小,为2.25.单个QTL平均解释表型变异介于14.10%~45.50%之间,其中贡献率大于20%的主效QTLs有9个.通过BLASTX与斑马鱼蛋白质序列数据库进行序列比对,找到了与斑马鱼酰基辅酶A脱氢酶蛋白、胰淀粉酶α2蛋白、Apoeb protein和甘油醛-3-磷酸脱氢酶蛋白同源的分子标记.本研究结果对分子标记辅助育种具有重要应用价值.     

鲤鱼乳酸脱氢酶活性的QTL检测

利用SSR分子标记结合“拟测交”策略, 以荷包红鲤抗寒品系和柏氏鲤远缘杂交所产生F1代为亲本及其F2代为作图群体, 应用Windows Map Manager 2.0软件的标记回归法进行乳酸脱氢酶(Lactate dehydrogenase, LDH)活性数量性状基因座单标记定位分析。在乳酸脱氢酶性状的标记回归研究中, 共发现12个标记与LDH性状关联, 对性状的贡献率为4.00%~10.00%, 其中HLJE222达到极显著水平(P<0.01)。为了进一步验证, 利用已有的生物信息学工具, 将与性状连锁的EST序列与公共数据库中的核酸序列进行同源搜索分析, 发现标记HLJE222的EST序列与斑马鱼DAZ associated protein 1的mRNA序列相匹配(相似性为94%), 与公共数据库中的已知蛋白序列进行相似性比较, 同样发现HLJE222的EST与斑马鱼DAZ associated protein 1相匹配(相似性为97%)。结果表明HLJE222位点与影响乳酸脱氢酶活性相关的基因连锁。     

草鱼Ⅰ型微卫星标记的发掘及其多态性检测

表达序列标签(EST)是发掘Ⅰ型微卫星标记的重要资源。研究运用生物信息学方法,从草鱼头肾组织3027条EST序列中搜索到322个微卫星位点,占整个EST数据库的10.6%。其中,二核苷酸重复位点151个(46.9%),三核苷酸重复位点137个(42.5%),四、五、六核苷酸重复位点较少;在二核苷酸重复位点中,AC/GT重复位点最为丰富,占二核苷酸重复位点总数的50.3%,AG/CT重复次之,占二核苷酸重复位点总数的40.4%,AT和GC重复较少。10个微卫星位点的多态性检测结果显示,4个位点在草鱼测试群体中呈多态性,多态性位点的平均多态信息含量(PIC)和平均遗传杂合度(H)分别为0.5236和0.5441,其中,2个多态性位点的PIC值大于0.5,呈现高度多态性特征。Ⅰ型微卫星标记将为草鱼遗传连锁图谱构建和QTL分析提供有效的基因分子标记。       

用于目标序列染色体定位的镜鲤BAC-FISH体系构建

为了构建用于镜鲤(Cyprinus carpio var. specularis)特定基因组序列染色体定位的实验体系, 在细菌人工染色体(Bacterial Artificial Chromosome, BAC)文库筛选池中对已知短序列基因组片段进行PCR扩增, 筛选出包含目标序列的BAC克隆, 提取BAC质粒进行缺刻平移标记制备探针, 开展荧光原位杂交(Fluorescence in situ hybridization, FISH)实验。通过对染色体片前处理、BAC质粒探针制备、C0t-1 DNA封闭基因组重复序列、预杂交、荧光染料选择、信号放大等一系列实验条件和方法的探索优化, 成功实现了目标序列在镜鲤有丝分裂中期染色体上的定位。定位对象既包括在染色体上有单一位点的序列, 如斑马鱼微卫星标记Z6884和Z4268, 也包括在染色体上有多个位点的重复序列, 如黄河鲤性别相关标记CCmf1。来自斑马鱼同一条染色体上的两个微卫星标记被分别定位于镜鲤不同染色体上, 为鲤鱼染色体数目加倍的进化假设提供了一项直接实验证据, 同时将现有遗传连锁图谱与染色体对应起来, 可作为染色体识别和细胞遗传学图谱构建的依据。黄河鲤性别相关重复序列被定位于不少于四条染色体上, 为性别决定相关基因的筛查提供了研究线索。这一BAC-FISH实验体系将成为鲤细胞遗传学图谱构建、基因组进化和比较基因组学研究中的重要研究工具。       

利用鸡F2资源群体构建1号染色体遗传连锁图谱

在鸡1号染色体上选取23个微卫星标记,利用东北农业大学鸡F2资源群体构建了遗传连锁图谱。选用369只F2个体用于基因型测定。结果表明23个微卫星位点除MCW0058为低度多态,其他位点均为中高度多态。构建的连锁图谱覆盖1号染色体全长,总共637.9 cM。MCW0115和ROS0025标记顺序与EL图谱不同,但与WAU图谱一致。其他标记顺序与3大参考家系标记顺序一致,图谱总长和标记间距离大于3大参考家系。此连锁图谱的构建为数量性状位点（QTL）定位奠定了良好的基础。     

鲤鱼多态性EST标记的筛选与特性分析

表达序列标签(Expressed sequence tag,EST)标记在基因组作图和分子标记辅助育种研究中具有重要价值.筛选和开发可用于遗传作图的鲤鱼多态性EST标记,对研究鲤鱼的基因组结构、遗传多样性研究和遗传育种有重要意义.本研究根据GenBank数据库中鲤鱼EST序列设计了67对EST引物,有47对在鲤鱼基因组DNA中成功扩增得到稳定的特异性条带,经单链构象多态性(SSCP)分析,12对(25.5%)引物扩增的EST在1个鲤鱼回交家系中具有多态性,其中6个为鲤鱼功能基因,5个与斑马鱼(Danio rerio)功能基因具较高相似性,1个与斑点叉尾(Ictalurus punctatus)功能基因具较高相似性.这些多态性的EST在鲤鱼二倍体家系和单倍体群体中的等位基因分离均符合孟德尔规律(1∶1或3∶1),单倍体中的SSCP条带数目为二倍体的二分之一.选用其中4对多态性的EST引物对洞庭湖鲤鱼进行初步群体遗传分析,结果显示基因多样性(H)为0.437,远低于微卫星标记所揭示的该群体基因多样性(接近于1).结果表明,多态性的EST-SSCP标记尽管遗传变异性较低,但是作为来源于编码区的Ⅰ型遗传标记,在遗传作图和种群遗传适应性等研究中有较好的应用潜力.     

利用两个测序水稻品种构建微卫星连锁图谱

利用已完成基因组测序的两个水稻品种日本晴和931l的数据库成功开发出水稻微卫星新标记,并利用由90个单株组成的日本晴×9311 F2作图群体,构建了一张包含152个SSR标记位点、覆盖基因组总长度2 455.7 cM的连锁图谱,有46个SSR新标记为自主开发,该图谱标记间的平均遗传距离为16.16 cM;并将未能在Temnykh等人(2001)构建的图谱上定位的微卫星标记RM345和RM494定位在第6染色体上.通过与Temnykh等人(2001)和兰涛等人(2003)所构建的图谱从作图群体的类型和大小、标记的类型和数量、标记在染色体上的线性排列顺序等几个方面进行比较,所绘制的图谱其标记在染色体线性排列上与Temnykh等人绘制的图谱具有很高的一致性,达93.81%.     

微卫星位点筛选方法综述

微卫星标记因其丰富的多态性和共显性等特点,已得到了广泛的应用.应用微卫星标记首先需要获得微卫星位点的序列信息,用来设计引物.获得微卫星位点的方法有多种,本文综述了获得和富集微卫星位点的常用方法.最简便、最省时的方法是从公共数据库(如EMBL、Genbank、EST数据库等)或已发表的文献中查找到微卫星位点,但只限于已经有序列数据发布的物种.第二种方法是种间转移扩增,即从相近物种的数据库中查找微卫星位点,或使用已有数据发表的遗传距离相近物种的微卫星标记.第三种方法是从基因组DNA中筛选微卫星位点,其中用于富集微卫星的方法有引物法、磁珠杂交法、尼龙膜杂交法以及RAPD技术法.     

Putative in silico mapping of DNA sequences to livestock genome maps using SSLP flanking sequences

In this study, an in silico approach was developed to identify homologies existing between livestock microsatellite flanking sequences and GenBank nucleotide sequences. Initially, 1955 bovine, 1570 porcine and 1121 chicken microsatellites were downloaded and the flanking sequences were compared with the nr and dbEST databases of GenBank. A total of 74 bovine, 44 porcine and 37 chicken microsatellite flanking sequences passed our criteria and had at least one significant match to human genomic sequence, genes/expressed sequence tags (ESTs) or both. GenBank annotation and BLAT searches of the UCSC human genome assembly revealed that 38 bovine, 13 porcine and 17 chicken microsatellite flanking sequences were highly similar to known human genes. Map locations were available for 67 bovine, 44 porcine and 21 chicken microsatellite flanking sequences, providing useful links in the comparative maps of humans and livestock. In support of our approach, 112 alignments with both microsatellite and match mapping information were located in the expected chromosomal regions based on previously reported syntenic relationships. The development of this in silico mapping approach has significantly increased the number of genes and EST sequences anchored to the bovine, porcine and chicken genome maps and the number of links in various human-livestock comparative maps.     

RS-SNP: a random-set method for genome-wide association studies

Background

Common carp is one of the most important aquaculture teleost fish in the world. Common carp and other closely related Cyprinidae species provide over 30% aquaculture production in the world. However, common carp genomic resources are still relatively underdeveloped. BAC end sequences (BES) are important resources for genome research on BAC-anchored genetic marker development, linkage map and physical map integration, and whole genome sequence assembling and scaffolding.

Result

To develop such valuable resources in common carp (Cyprinus carpio), a total of 40,224 BAC clones were sequenced on both ends, generating 65,720 clean BES with an average read length of 647 bp after sequence processing, representing 42,522,168 bp or 2.5% of common carp genome. The first survey of common carp genome was conducted with various bioinformatics tools. The common carp genome contains over 17.3% of repetitive elements with GC content of 36.8% and 518 transposon ORFs. To identify and develop BAC-anchored microsatellite markers, a total of 13,581 microsatellites were detected from 10,355 BES. The coding region of 7,127 genes were recognized from 9,443 BES on 7,453 BACs, with 1,990 BACs have genes on both ends. To evaluate the similarity to the genome of closely related zebrafish, BES of common carp were aligned against zebrafish genome. A total of 39,335 BES of common carp have conserved homologs on zebrafish genome which demonstrated the high similarity between zebrafish and common carp genomes, indicating the feasibility of comparative mapping between zebrafish and common carp once we have physical map of common carp.

Conclusion

BAC end sequences are great resources for the first genome wide survey of common carp. The repetitive DNA was estimated to be approximate 28% of common carp genome, indicating the higher complexity of the genome. Comparative analysis had mapped around 40,000 BES to zebrafish genome and established over 3,100 microsyntenies, covering over 50% of the zebrafish genome. BES of common carp are tremendous tools for comparative mapping between the two closely related species, zebrafish and common carp, which should facilitate both structural and functional genome analysis in common carp.     

Human chromosome 3: integration of 60 NotI clones into a physical and gene map

Sequence tagged sites generated for 60 NotI clones (NotI-STSs) from human chromosome 3-specific NotI-jumping and NotI-linking libraries were physically located using PCR screening of a radiation hybrid (RH) GeneBridge4 panel. The NotI map of chromosome 3 was generated using these RH-mapping data and those obtained earlier by FISH and sequencing of the corresponding NotI clones. The sequences of the NotI clones showed significant homologies with known genes and/or ESTs for 58 NotI-STSs (97%). These 58 NotI clones displayed 91-100% identity to 54 genes and 23 cDNA/EST clones. One known and two hypothetical protein-coding genes were localized for the first time and nine cDNA clones (unknown genes) were also carefully mapped only in this work. Three newly mapped genes are histone gene H1X (NR1-BK20C) and genes for hypothetical proteins THC1032178 and THC1024604 (NL1-243).     

Bioinformatic identification of genes encoding C1q-domain containing proteins in zebrafish

C1q is the first subcomponent of classical pathway in the complement system and a major link between innate and acquired immunities. The globular (gC1q) domain similar with C1q was also found in many non-complement C1q-domain-containing (C1qDC) proteins which have similar crystal structure to that of the multifunctional tumor necrosis factor (TNF) ligand family, and also have diverse functions. In this study, we identified a total of 52 independent gene sequences encoding C1q-domain-containing proteins through comprehensive searches of zebrafish genome, cDNA and EST databases. In comparison to 31 orthologous genes in human and different numbers in other species, a significant selective pressure was suggested during vertebrate evolution. Domain organization of C1q-domain-containing (C1qDC) proteins mainly includes a leading signal peptide, a collagen-like region of variable length, and a C-terminal C1q domain. There are 11 highly conserved residues within the C1q domain, among which 2 are invariant within the zebrafish gene set. A more extensive database searches also revealed homologous C1qDC proteins in other vertebrates, invertebrates and even bacterium, but no homologous sequences for encoding C1qDC proteins were found in many species that have a more recent evolutionary history with zebrafish. Therefore, further studies on C1q-domain-containing genes among different species will help us understand evolutionary mechanism of innate and acquired immunities.     

Proteomic profiling and protein identification by MALDI-TOF mass spectrometry in unsequenced parasitic nematodes

Lack of genomic sequence data and the relatively high cost of tandem mass spectrometry have hampered proteomic investigations into helminths, such as resolving the mechanism underpinning globally reported anthelmintic resistance. Whilst detailed mechanisms of resistance remain unknown for the majority of drug-parasite interactions, gene mutations and changes in gene and protein expression are proposed key aspects of resistance. Comparative proteomic analysis of drug-resistant and -susceptible nematodes may reveal protein profiles reflecting drug-related phenotypes. Using the gastro-intestinal nematode, Haemonchus contortus as case study, we report the application of freely available expressed sequence tag (EST) datasets to support proteomic studies in unsequenced nematodes. EST datasets were translated to theoretical protein sequences to generate a searchable database. In conjunction with matrix-assisted laser desorption ionisation time-of-flight mass spectrometry (MALDI-TOF-MS), Peptide Mass Fingerprint (PMF) searching of databases enabled a cost-effective protein identification strategy. The effectiveness of this approach was verified in comparison with MS/MS de novo sequencing with searching of the same EST protein database and subsequent searches of the NCBInr protein database using the Basic Local Alignment Search Tool (BLAST) to provide protein annotation. Of 100 proteins from 2-DE gel spots, 62 were identified by MALDI-TOF-MS and PMF searching of the EST database. Twenty randomly selected spots were analysed by electrospray MS/MS and MASCOT Ion Searches of the same database. The resulting sequences were subjected to BLAST searches of the NCBI protein database to provide annotation of the proteins and confirm concordance in protein identity from both approaches. Further confirmation of protein identifications from the MS/MS data were obtained by de novo sequencing of peptides, followed by FASTS algorithm searches of the EST putative protein database. This study demonstrates the cost-effective use of available EST databases and inexpensive, accessible MALDI-TOF MS in conjunction with PMF for reliable protein identification in unsequenced organisms.     

Localization of retina/pineal-expressed sequences: identification of novel candidate genes for inherited retinal disorders.

More than 100 genes causing inherited retinal diseases have been mapped to chromosomal locations, but less than half of these genes have been cloned. Mutations in many retina/pineal-specific genes are known to cause inherited retinal diseases. Examples include mutations in arrestin, rhodopsin kinase, and the cone-rod homeobox gene, CRX. To identify additional candidate genes for inherited retinal disorders, novel retina/pineal-expressed EST clusters were identified from the TIGR Human Gene Index database and mapped to specific chromosomal sites. After known human gene sequences were excluded, and repeat sequences were masked, 26 novel retina and pineal gland cDNA clusters were identified. The retinal expression of each novel EST cluster was confirmed by PCR assay of a retinal cDNA library, and each cluster was localized in the genome using the GeneBridge 4.0 radiation hybrid panel. In silico expression data from the TIGR database suggest that these EST clusters are retina/pineal-specific or predominantly expressed in these tissues. This combination of database analysis and laboratory investigation has localized several EST clusters that are potential candidates for genes causing inherited retinopathy.     

Bioinformatic identification of genes encoding Clq-domaincontaining proteins in zebrafish

C1q is the first subcomponent of classical pathway in the complement system and a major link between innate and acquired immunities. The globular （gC1q） domain similar with C1q was also found in many non-complement C1q-domain-containing （C1qDC） proteins which have similar crystal structure to that of the multifunctional tumor necrosis factor （TNF） ligand family, and also have diverse functions. In this study, we identified a total of 52 independent gene sequences encoding C1q-domain-containing proteins through comprehensive searches of zebrafish genome, cDNA and EST databases. In comparison to 31 orthologous genes in human and different numbers in other species, a significant selective pressure was suggested during vertebrate evolution. Domain organization of C1q-domain-containing （C1qDC） proteins mainly includes a leading signal peptide, a collagen-like region of variable length, and a C-terminal C1q domain. There are 11 highly conserved residues within the C1q domain, among which 2 are invariant within the zebrafish gene set. A more extensive database searches also revealed homologous C1qDC proteins in other vertebrates, invertebrates and even bacterium, but no homologous sequences for encoding C1qDC proteins were found in many species that have a more recent evolutionary history with zebrafish. Therefore, further studies on C1q-domain-containing genes among different species will help us understand evolutionary mechanism of innate and acquired immunities.