整合数字基因表达谱与全基因组关联分析鉴定猪血液性状候选基因  

Identification of Candidate Genes for Hematological Traits by Integrating Gene Expression Profiling and Genome-Wide Association Study in a Porcine Model

作　　者：徐盼[1] 张震[1] 章峰[1] 杨斌[1] 段艳宇[1] 

机构地区：[1]江西农业大学省部共建猪遗传改良与养殖技术国家重点实验室,南昌330045

出　　处：《中国农业科学》2016年第2期348-360,共13页Scientia Agricultura Sinica

基　　金：国家基金(青年项目)(31301950);江西省自然科学基金(2010GQN0045)

摘　　要：【目的】整合数字基因表达谱与全基因组关联分析鉴定白色杜洛克×二花脸F_2资源群体的血液性状候选基因。【方法】白色杜洛克×二花脸F_2资源群体在(240±3)d屠宰,收集血液于抗凝管中进行血常规检测。利用Illumina 60K SNP芯片对1 020头F_2资源群体进行基因分型。剔除基因型检出率<90%和孟德尔错误检出率>5%的个体。检出率<95%、次等位基因频率<5%、哈代-温伯格检验(HWE)P<5×10^(-6)、与性染色体连锁疑似常染色体的SNP被筛除。利用Illumina GA II测序仪测序对502头F_2资源群体的肝脏进行数字基因表达谱测序。测序得到的原始数据经过滤获得清洁标签后与参考标签数据库比对,将能唯一比对到参考基因序列的清洁标签数量进行标准化处理以获得标准化的基因表达量。每个转录本的表达水平进一步转化为lg2值。在少于20%的个体中表达的转录本被滤去。表型性状和基因表达性状使用R程序包中Gen ABEL内polygentic功能进行性别、批次和亲缘关系的校正。其残差使用R程序包中斯皮尔曼系数评估基因表达水平与表型数据的关联性,设定保守阈值P<5×10^(-4)时调整多重检验。将检测到的表达数量性状位点(eQTL)及其对应基因根据其位置相对照的关系进行绘图。搜寻前期GWAS最高点5.0 Mb区域内eQTL结合GWAS结果进行综合分析。Gene Ontology&KEGG pathway富集分析使用在线工具DAVID。基因共表达网络使用在线工具Gene MANIA进行构建。【结果】白色杜洛克×二花脸F_2资源群体中502个个体的20 108个肝脏转录本通过了质检。当P<5×10^(-4)时鉴别到与血红蛋白(HGB)、红细胞数目(RBC)、红细胞压积(HCT)、平均红细胞体积(MCV)、平均红细胞血红蛋白含量(MCH)和白细胞数目(WBC)关联的转录本共259个。有34个转录本与一个以上表型关联。用上述血细胞性状关联的转录本进行eQTL定位,当阈值P<10^(-5)时得到304个eQTL,【Objective】We herein integrated digital gene expression profiling and genome-wide association study in a White Duroc × Erhualian F_2 resource population to identify candidate genes for hematological traits.【Method】The White Duroc × Erhualian F_2 resource population were slaughtered at 240 ± 3 days. Blood was collected in anticoagulation tubes. A set of hematological parameters were measured using a whole blood analyzer. The 1 020 F_2 pigs were genotyped using an Illumina porcine 60 K SNP chip. Individuals with genotype-missing rates 10% and Mendellian errors 5% were removed. SNPs with a call rate 95%, minor allele frequency 5%, P value 5 × 10^（-6） for Hardy Weinberg equilibrium（HWE）, and the autosomal SNPs that were linked to sex chromosome were excluded. The liver samples of the 502 F_2 pigs were then performed digital gene expression profiles sequencing on Illumina GA II. The raw tags were filtered to obtain clean tags. The clean tags were uniquely mapped to the reference gene sequences and were defined as unambiguous clean tags. The number of unambiguous clean tags was normalized to represent the expression level of each transcript. The expression level of each transcript was further transformed to lg2 value. The transcripts that expressed less than 20% of individuals were rejected. The traits of phenotype and gene expression were adjusted for sex, batch and kinship using polygentic function of Gen ABEL in R package. The correlations between gene expressions and phenotypic traits were evaluated using the residuals by Spearman＇s correlation coefficient with a conservative threshold P 0.0005. Positions of detected eQTL were plotted against the positions of the genes for which that eQTL were found. We also searched the eQTL within the 5.0 Mb region of the peak SNP of GWAS and performed an integrated analysis of eQTL and GWAS. Gene Ontology KEGG pathway enrichment analysis was implemented by DAVID online tools and gene co-expression network was constructed by Gene MANIA online tools. 【Resul

关 键 词：猪 血液性状 数字基因表达谱 表达数量性状位点 全基因组关联分析 候选基因 

分 类 号：S828[农业科学—畜牧学]