大规模平行测序技术(MPSS)研究进展

大规模平行测序技术（massively parallel signature sequencing,MPSS）是以DNA测序为基础的大规模高通量基因分析新技术,通过标签库的建立、微珠与标签的连接、酶切连接反应和生物信息分析等步骤,获得基因表达序列.MPSS具有能测定表达水平较低、差异较小的基因,不必预先知道基因的序列,自动化和高通量等特点,是值得推广的技术.     

Massively parallel sequencing of Chikso (Korean brindle cattle) to discover genome-wide SNPs and InDels

Since the completion of the bovine sequencing projects, a substantial number of genetic variations such as single nucleotide polymorphisms have become available across the cattle genome. Recently, cataloguing such genetic variations has been accelerated using massively parallel sequencing technology. However, most of the recent studies have been concentrated on European Bos taurus cattle breeds, resulting in a severe lack of knowledge for valuable native cattle genetic resources worldwide. Here, we present the first whole-genome sequencing results for an endangered Korean native cattle breed, Chikso, using the Illumina HiSeq 2,000 sequencing platform. The genome of a Chikso bull was sequenced to approximately 25.3-fold coverage with 98.8% of the bovine reference genome sequence (UMD 3.1) covered. In total, 5,874,026 single nucleotide polymorphisms and 551,363 insertion/deletions were identified across all 29 autosomes and the X-chromosome, of which 45% and 75% were previously unknown, respectively. Most of the variations (92.7% of single nucleotide polymorphisms and 92.9% of insertion/deletions) were located in intergenic and intron regions. A total of 16,273 single nucleotide polymorphisms causing missense mutations were detected in 7,111 genes throughout the genome, which could potentially contribute to variation in economically important traits in Chikso. This study provides a valuable resource for further investigations of the genetic mechanisms underlying traits of interest in cattle, and for the development of improved genomics-based breeding tools.     

Genome analysis of the domestic dog (Korean Jindo) by massively parallel sequencing

Although pioneering sequencing projects have shed light on the boxer and poodle genomes, a number of challenges need to be met before the sequencing and annotation of the dog genome can be considered complete. Here, we present the DNA sequence of the Jindo dog genome, sequenced to 45-fold average coverage using Illumina massively parallel sequencing technology. A comparison of the sequence to the reference boxer genome led to the identification of 4 675 437 single nucleotide polymorphisms (SNPs, including 3 346 058 novel SNPs), 71 642 indels and 8131 structural variations. Of these, 339 non-synonymous SNPs and 3 indels are located within coding sequences (CDS). In particular, 3 non-synonymous SNPs and a 26-bp deletion occur in the TCOF1 locus, implying that the difference observed in cranial facial morphology between Jindo and boxer dogs might be influenced by those variations. Through the annotation of the Jindo olfactory receptor gene family, we found 2 unique olfactory receptor genes and 236 olfactory receptor genes harbouring non-synonymous homozygous SNPs that are likely to affect smelling capability. In addition, we determined the DNA sequence of the Jindo dog mitochondrial genome and identified Jindo dog-specific mtDNA genotypes. This Jindo genome data upgrade our understanding of dog genomic architecture and will be a very valuable resource for investigating not only dog genetics and genomics but also human and dog disease genetics and comparative genomics.     

基因表达快速分析新技术-MPSS简介

介绍了一种新的基因表达快速分析技术-MPSS技术的基本原理、技术路线及其优点和应用。     

Development and validation of a 10‐gene prognostic signature for acute myeloid leukaemia

Acute myeloid leukaemia (AML) is the most common type of adult acute leukaemia and has a poor prognosis. Thus, optimal risk stratification is of greatest importance for reasonable choice of treatment and prognostic evaluation. For our study, a total of 1707 samples of AML patients from three public databases were divided into meta‐training, meta‐testing and validation sets. The meta‐training set was used to build risk prediction model, and the other four data sets were employed for validation. By log‐rank test and univariate COX regression analysis as well as LASSO‐COX, AML patients were divided into high‐risk and low‐risk groups based on AML risk score (AMLRS) which was constituted by 10 survival‐related genes. In meta‐training, meta‐testing and validation sets, the patient in the low‐risk group all had a significantly longer OS (overall survival) than those in the high‐risk group (P < .001), and the area under ROC curve (AUC) by time‐dependent ROC was 0.5854‐0.7905 for 1 year, 0.6652‐0.8066 for 3 years and 0.6622‐0.8034 for 5 years. Multivariate COX regression analysis indicated that AMLRS was an independent prognostic factor in four data sets. Nomogram combining the AMLRS and two clinical parameters performed well in predicting 1‐year, 3‐year and 5‐year OS. Finally, we created a web‐based prognostic model to predict the prognosis of AML patients ( https://tcgi.shinyapps.io/amlrs_nomogram/ ).     

基于荧光微球液相基因表达阵列的建立

建立基于荧光微球的液相基因表达阵列,用于特定基因组合的表达谱分析.采用带有不同强度荧光鉴别信号的羧基化微球,与氨基修饰的不同标签寡核苷酸序列化学偶联,制成微球阵列.多重连接依赖的探针扩增技术(MLPA)用于扩增靶基因核苷酸序列,即通过RNA标本六随机引物逆转成cDNA,与不同基因特异性的一对探针杂交,耐高温的连接酶联接,最后采用生物素标记的同一对引物扩增.PCR产物与微球阵列液相杂交,加入链亲和素标记的PE染料,上流式细胞仪检测.应用这一系统检测骨髓增生异常综合症中难治性贫血(RA)、难治性贫血伴原始细胞增多(RAEB)、难治性贫血伴转化中原始细胞增多(RAEBt)、急性髓细胞性白血病(AML)和其他组(包括再生障碍性贫血、血小板减少、巨幼贫、溶贫等)差异表达谱,差异表达结果用实时荧光定量PCR验证.共建立了5个基因的微球阵列,分别为Rap1GAP、RAC2、SPA1、RhoBTB3和内参GAPDH,每个基因检测的线性范围为0.002 5～0.1μmol,液相表达阵列具有良好的特异性和重复性(P<0.001).检测RA、RAEB、RAEBt、AML和其他组差异表达发现,RAC2、RhoBTB3、SPA-1和Rap1GAP各组间有显著性差异性存在(分别为P<0.000 1,P=0.049 1,P=0.020 6和P=0.004 6),其差异显著性与实时荧光定量PCR一致,泊松相关系数分别为0.930,0.946,0.945和0.921,具显著性(P<0.001).结果表明,成功建立了基于荧光微球的液相基因表达阵列,其敏感性高、特异性强、重复性好.     

Genomic DNA macroarrays as a tool for analysis of gene expression in Leishmania

Gene-array technologies have been applied in a wide number of organisms to study gene expression profiling under several physiological and experimental conditions. Gene-array implementations combined with the information arising from emerging genome sequencing projects are expected to be in the near future a major tool to characterize genes involved in different processes. So far, gene expression profile studies in trypanosomatids have been performed in microarrays that use a glass support to immobilize fragments of genomic DNA followed by fluorescent detection. Here, we wanted to test the potential of genomic DNA macroarrays of Leishmania infantum using nylon membranes and radioactive detection. Nylon macroarrays present a number of advantages since the processing of the membranes is based on standard Southern blotting protocols familiar to molecular biologists, and the data acquisition equipment is available to most research institutions. Nylon macroarrays were employed to search for genes showing increased mRNA abundance during an axenic differentiation of L. infantum promastigotes to amastigotes. Several clones were rescued and, after validation by Northern blot assays, these L. infantum sequences were used to screen the Leishmania major gene database. The L. major contigs with high homology to the L. infantum sequences allowed a consistent identification of the regulated genes.     

Cross-platform compatibility of Hi-Plex,a streamlined approach for targeted massively parallel sequencing

Although per-base sequencing costs have decreased during recent years, library preparation for targeted massively parallel sequencing remains constrained by high reagent cost, limited design flexibility, and protocol complexity. To address these limitations, we previously developed Hi-Plex, a polymerase chain reaction (PCR) massively parallel sequencing strategy for screening panels of genomic target regions. Here, we demonstrate that Hi-Plex applied with hybrid adapters can generate a library suitable for sequencing with both the Ion Torrent and the TruSeq chemistries and that adjusting primer concentrations improves coverage uniformity. These results expand Hi-Plex capabilities as an accurate, affordable, flexible, and rapid approach for various genetic screening applications.     

Monitoring genome-wide changes in gene expression in response to endogenous cytokinin reveals targets in Arabidopsis thaliana

Cytokinins have been implicated in developmental and growth processes in plants including cell division, chloroplast biogenesis, shoot meristem initiation and senescence. The regulation of these processes requires changes in cytokinin-responsive gene expression. Here, we induced the expression of a bacterial isopentenyl transferase gene, IPT, in transgenic Arabidopsis thaliana seedlings to study the regulation of genome-wide gene expression in response to endogenous cytokinin. Using MPSS (massively parallel signature sequencing) we identified 823 and 917 genes that were up- and downregulated, respectively, following 24 h of IPT induction. When comparing the response to cytokinin after 6 and 24 h, we identified different clusters of genes showing a similar course of regulation. Our study provides researchers with the opportunity to rapidly assess whether genes of interest are regulated by cytokinins.     

Gene expression profiling of the pH response in Shigella flexneri 2a

The pH response of Shigella flexneri 2a 301 was identified by gene expression profiling. Gene expression profiles of cells grown in pH 4.5 or 8.6 were compared with the profiles of cells grown at pH 7.0. Differential expression was observed for 307 genes: 97 were acid up-regulated, 102 were acid down-regulated, 91 were base up-regulated, and 86 were base down-regulated. Twenty-seven genes were found to be both acid and base up-regulated, and 29 genes were both acid and base down-regulated. This study showed that (1) the most pH-dependent genes regulate energy metabolism; (2) the RpoS-dependent acid-resistance system is induced, while the glutamate-dependent acid resistance system is not; (3) high pH up-regulates some virulence genes, while low pH down-regulates them, consistent with Shigella infection of the low gut; and (4) several cross-stress response genes are induced by pH changes. These results also illustrate that many unknown genes are significantly regulated under acid or basic conditions, providing researchers with important information to characterize their function.     

The Arabidopsis co-expression tool (ACT): a WWW-based tool and database for microarray-based gene expression analysis

We present a new WWW-based tool for plant gene analysis, the Arabidopsis Co-Expression Tool (ACT), based on a large Arabidopsis thaliana microarray data set obtained from the Nottingham Arabidopsis Stock Centre. The co-expression analysis tool allows users to identify genes whose expression patterns are correlated across selected experiments or the complete data set. Results are accompanied by estimates of the statistical significance of the correlation relationships, expressed as probability (P) and expectation (E) values. Additionally, highly ranked genes on a correlation list can be examined using the novel clique finder tool to determine the sets of genes most likely to be regulated in a similar manner. In combination, these tools offer three levels of analysis: creation of correlation lists of co-expressed genes, refinement of these lists using two-dimensional scatter plots, and dissection into cliques of co-regulated genes. We illustrate the applications of the software by analysing genes encoding functionally related proteins, as well as pathways involved in plant responses to environmental stimuli. These analyses demonstrate novel biological relationships underlying the observed gene co-expression patterns. To demonstrate the ability of the software to develop testable hypotheses on gene function within a defined biological process we have used the example of cell wall biosynthesis genes. The resource is freely available at http://www.arabidopsis.leeds.ac.uk/ACT/     

拟南芥不定芽发生早期的数字基因表达谱分析

目前,有关不定芽发生的研究主要集中在单基因的调控方面,缺乏转录组方面的系统研究.利用RNA-seq高通量测序技术在全基因组范围内检测了不定芽发生早期的基因表达谱,共检测到2457个差异表达基因.这些基因参与了激素代谢和信号转导、愈伤组织和侧根的形成、茎顶端分生组织的发育和光合作用等过程.进一步的途径富集分析表明,不定芽发生早期苯丙氨酸代谢和苯丙胺素合成等途径相关的基因显著富集.并且苯丙氨酸可以显著抑制不定芽的发生,暗示了苯丙氨酸代谢和苯丙胺素的合成可能在不定芽发生过程起着重要的作用.     

中华蜜蜂DNA甲基化转移酶Dnmt3基因克隆及表达谱分析

为探究中华蜜蜂Apis cerana cerana的DNA甲基化模式, 本研究采用RT PCR技术克隆了中华蜜蜂DNA甲基化转移酶3（Dnmt3）基因(GenBank登录号为JQ740768); 采用荧光定量PCR检测不同发育时期工蜂（4日龄蛹, 1, 7和30日龄成年蜂及产卵工蜂）和蜂王（4日龄蛹, 1日龄蜂王和产卵蜂王）头部的Dnmt3基因mRNA的表达量。结果表明： 该基因cDNA序列全长2 277 bp, 编码758个氨基酸残基, 预测的蛋白分子量为88.24 kD, 等电点为7.85。将中华蜜蜂与其他物种的Dnmt3基因的结构域进行比对, 同时将该基因推导的氨基酸序列与其他物种的Dnmt3氨基酸序列进行同源性比对和系统发育分析, 发现与西方蜜蜂的Dnmt3序列一致性高达99%。该基因在工蜂和蜂王不同发育时期均有表达, 1日龄工蜂与7日龄工蜂中没有显著差异(P>0.05), 30日龄工蜂中的表达量显著高于前两者 (P<0.05); 蜂王蛹中的表达量显著高于工蜂蛹 (P<0.05); 1日龄的蜂王中的表达量显著高于1日龄的工蜂(P<0.05); 产卵工蜂与产卵蜂王中的表达量没有差异（P>0.05）。这种表达情况提示其可能与工蜂劳动分工及蜜蜂卵巢发育有关。