Construction of Three High-Density Genetic Linkage Maps and Dynamic QTL Mapping of Growth Traits in Yellow River Carp (Cyprinus carpio haematopterus)

To provide the theoretical basis for researching growth, development, and molecular marker-assisted breeding of the economically important Yellow River carp (Cyprinus carpio haematopterus) using dynamic quantitative trait locus (QTL) mapping, we constructed three genetic linkage maps from 207 progeny using a new modified genotyping-by-sequencing method. The three maps contained 16,886, 16,548, and 7482 single nucleotide polymorphism markers, respectively, with an average interval of 0.36 cM, 0.45 cM, and 1.00 cM. We identified 148 QTLs related to four growth traits that were located on 25 chromosomes from three growth stages of Yellow River carp. A total of 32, 36, 43, and 37 QTLs were associated with body length, height, width, and weight, respectively. Among them, 47 QTLs were detected for only one growth trait in one stage, but all of the other QTLs were co-localized. Of the 14 main QTLs, 13 were located on chromosome 12, which suggests the presence of growth-related genes on this chromosome. We then detected 17 candidate genes within 50 K upstream and downstream of the 14 main QTLs. This is the first report of the dynamic QTL mapping of growth traits of Yellow River carp, and the results can be used in future studies of growth, development, and molecular-assisted breeding of this species.     

Identification of genetic loci associated with growth traits at weaning in yak through a genome-wide association study

A multilocus GWAS was performed to explore the genetic architecture of four growth traits in yak. In total, 354 female yaks for which measurements of body weight (BW), withers height (WH), body length (BL) and chest girth (CG) at weaning were available underwent genotyping with the Illumina BovineHD BeadChip (770K). After quality control, we retained 98 688 SNPs and 354 animals for GWAS analysis. We identified seven, 18, seven and nine SNPs (corresponding to seven, 17, seven and eight candidate genes) associated with BW, WH, BL and CG at weaning respectively. Interestingly, most of these candidate genes were reported to be involved in growth-related processes such as muscle formation, lipid deposition, feed efficiency, carcass composition and development of the central and peripheral nervous system. Our results offer novel insight into the molecular architecture underpinning yak growth traits. Further functional analyses are thus warranted to explore the molecular mechanisms whereby these genes affect these traits of interest.     

Genome-Scale Association Study of Abnormal Scale Pattern in Yellow River Carp Identified Previously Known Causative Gene in European Mirror Carp

Common carp (Cyprinus carpio) is one of the most widely studied fish species due to its great economic value and strong environmental adaptability. Scattered scale, a typical phenotype of the mirror carp that is derived from Europe, has never been observed in the Yellow River carp previously. We recently identified approximately one fourth of the F1 progenies displaying scattered scale in a full-sib Yellow River carp family in our breeding program, despite both parents that showed wild type with normal scale patterns. This family provides us unique materials to investigate the genetic basis underlying the abnormal scale mutant in Yellow River carp population. Genome-wide association study (GWAS) and association mapping were performed based on genome-wide single nucleotide polymorphisms (SNP) genotyped with common carp 250 K SNP genotyping array in 82 samples of the Yellow River carp family. We identified a 1.4 Mb genome region that was significantly associated with abnormal scattered scale patterns. We further identified a deletion mutation in fibroblast growth factor receptor 1 a1 (fgfr1a1) gene within this genome region. Amplification and sequencing analysis of this gene revealed a 311-bp deletion in intron 10 and exon 11, which proved that fgfr1a1 could be the causal gene responsible for abnormal scattered scale in the Yellow River carp family. Since similar fragment mutation with 306-bp and 310-bp deletions had been previously reported as causal mutation of scattered scale patterns in the mirror carp, we speculate that either the deletion mutation was introduced from Europe-derived mirror carp or the deletion independently occurred in the mutation hotspot in fgfr1a1 gene. The results provided insights into the genetic basis of scale pattern mutant in Yellow River carp population, which would help us to eliminate the recessive allele of the abnormal scale patterns in Yellow River carp population by molecular marker-assisted breeding.     

京海黄鸡体重性状全基因组关联分析

体重性状是肉鸡重要的经济性状。为了寻找可用于京海黄鸡体重性状遗传改良的分子标记及候选基因,本文以400只京海黄鸡核心群母鸡为基础,测定了0~14周龄体重,利用简化基因组测序技术(Specific-locus amplified fragment sequencing, SLAF-seq)对京海黄鸡体重性状进行全基因组关联研究(Genome-wide association stndy, GWAS),筛选与京海黄鸡体重性状相关的SNPs位点。结果共检测到100个与京海黄鸡体重相关的SNPs位点,其中15个位点效应达到全基因组显著水平(P<1.87E-06),85个位点效应达到全基因组潜在显著水平(P<3.73E-05)。通过筛选每个显著SNP周围1 Mb区域内的基因,共找到9个可能的候选基因,其中FAM124A(Family with sequence similarity 124A)、QDPR(Quinoid dihydropteridine reductase)、WDR1(WD repeat domain 1)和SLC2A9(Solute carrier family 2 (facilitated glucose transporter), member 9) 4个基因可能是影响体重性状的重要候选基因。同时还发现,4号染色体75.6~80.7 Mb区域集中了大部分与京海黄鸡中后期体重性状显著相关的SNPs位点,该区域可能是影响京海黄鸡中后期生长体重的重要候选区域。     

小麦耐低磷相关性状的全基因组关联分析

通过对小麦耐低磷相关性状进行全基因组关联分析(GWAS,genome-wide association study),挖掘与小麦耐低磷性显著相关的单核苷酸多态性标记(SNP,single nucleotide polymorphism)位点及候选基因,为小麦耐低磷性状的遗传基础和分子机制研究提供理论参考。本试验以198份黄淮麦区小麦品种(系)为试验材料,设置低磷和正常磷营养液水培试验,利用小麦35K芯片对分布于小麦全基因组的11896个SNP,采用Q+K关联模型对小麦耐低磷性相关性状进行关联分析。结果表明,小麦耐低磷性状表现出广泛的表型变异,变异系数为15.65%~26.59%,多态性信息含量(PIC,polymorphic information content)为0.095~0.500。群体结构分析表明,试验所用自然群体可分为2个亚群,GWAS共检测到67个与小麦耐低磷相关性状显著关联的SNP位点(P≤0.001),这些位点分布在除3A、3B和3D以外的18条染色体上,单个SNP位点可解释5.826%~9.552%的表型变异。在这些显著位点中有4个SNP位点同时关联到了2个不同的耐低磷性状。对67个SNP位点进行发掘,筛选到7个可能与小麦耐低磷性有关的候选基因。TraesCS6A02G001000和TraesCS6A02G001100在锌指合成中有重要作用;TraesCS6A02G118100可能为低磷胁迫诱导基因;TraesCS5D02G536400、TraesCS1B02G154200和TraesCS5D02G536500与低磷胁迫相关酶类基因家族有关;TraesCS1D02G231200与植物DUF 538结构域蛋白有关,是植物胁迫相关调控蛋白候选基因。     

Genome-wide association study reveals the genetic basis of growth trait in yellow catfish with sexual size dimorphism

Sexual size dimorphism has been widely observed in a large number of animals including fish species. Genome-wide association study (GWAS) is a powerful tool to dissect the genetic basis of complex traits, whereas the sex-differences in the genomics of animal complex traits have been ignored in the GWAS analysis. Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish in China with significant sexual size dimorphism. In this study, GWAS was conducted to identify candidate SNPs and genes related to body length (BL) and body weight (BW) in 125 female yellow catfish from a breeding population. In total, one BL-related SNP and three BW-related SNPs were identified to be significantly associated with the traits. Besides, one of these SNPs (Chr15:19195072) was shared in both the BW and BL traits in female yellow catfish, which was further validated in 185 male individuals and located on the exon of stat5b gene. Transgenic yellow catfish and zebrafish that expressed yellow catfish stat5b showed increased growth rate and reduction of sexual size dimorphism. These results not only reveal the genetic basis of growth trait and sexual size dimorphism in fish species, but also provide useful information for the marker-assisted breeding in yellow catfish.     

生物信息学分析方法I: 全基因组关联分析概述

全基因组关联分析(GWAS)是动植物复杂性状相关基因定位的常用手段。高通量基因分型技术的应用极大地推动了GWAS的发展。在植物中, 利用GWAS不仅能够以较高的分辨率在全基因组水平鉴定出各种自然群体特定性状相关的基因或区间, 而且可揭示表型变异的遗传架构全景图。目前, 人们利用GWAS分析方法已在拟南芥(Arabidopsis thaliana)、水稻(Oryza sativa)、小麦(Triticum aestivum)、玉米(Zea mays)和大豆(Glycine max)等模式植物和重要农作物品系中发掘出与各种性状显著相关的数量性状座位(QTL)及其候选基因位点, 阐明了这些性状的遗传基础, 并为揭示这些性状背后的分子机理提供候选基因, 也为作物高产优质品种的选育提供了理论依据。该文对GWAS的方法、影响因素及数据分析流程进行了详细描述, 以期为相关研究提供参考。     

银川番茄斑萎病毒的分子鉴定

为明确银川番茄(Lycopersicon esculentum)是否遭受了番茄斑萎病毒(TSWV)的危害, 采用国家标准TSWV RT- PCR检测技术对银川番茄上采集的14份疑似感染TSWV病叶样本进行分子鉴定, 对克隆得到的核衣壳蛋白基因N (Nucleocapsid)序列进行多序列比对和系统进化树分析, 随后对PCR阳性样本进行蛋白检测。结果表明, 14份病叶样本中有8份扩增出长度为394 bp的TSWV N基因序列, 且8条序列完全一致; 获得的银川番茄TSWV分离物与云南番茄、中国莴苣(Lactuca sativa)、中国鸢尾(Iris tectorum)和重庆辣椒(Capsicum annuum) TSWV分离物相对近缘, 与山东、黑龙江和北京等地及国外TSWV分离物相对远缘; 利用TSWV的抗体通过Western blot对8个PCR阳性样本进一步检测, 结果也证实8个阳性样本中存在TSWV感染。该研究首次通过分子鉴定及蛋白检测证明银川番茄上存在TSWV感染, 需要加快抗TSWV番茄品种的选育工作。     

鲢EST-SSR标记的开发及其与生长性状关联性分析

研究采用Illumina高通量测序技术对鲢(Hypophthalmichthys molitrix)卵组织进行转录组测序分析, 筛选多态EST-SSR标记并与生长性状进行关联性分析。研究结果显示: (1)测序组装总Unigene 77129个, SSR共24458个, 占所测EST数据的31.7%。鲢包含SSR的EST序列主要有二碱基、三碱基、四碱基和五碱基组成。二碱基中, AC含量最丰富; (2)随机挑选175条EST序列设计SSR引物, 筛选出多态EST-SSR标记36个, Blast结果显示有10条已知基因功能; (3)全同胞家系检测结果显示有两个标记与鲢的生长性状相关。SCE26与体长和体重显著相关(P<0.05), SCE65与肥满度显著相关(P<0.05)。研究结果将为鲢的分子标记辅助育种的应用提供基础理论指导。     

Identification of pleiotropic genes and gene sets underlying growth and immunity traits: a case study on Meishan pigs

Both growth and immune capacity are important traits in animal breeding. The animal quantitative trait loci (QTL) database is a valuable resource and can be used for interpreting the genetic mechanisms that underlie growth and immune traits. However, QTL intervals often involve too many candidate genes to find the true causal genes. Therefore, the aim of this study was to provide an effective annotation pipeline that can make full use of the information of Gene Ontology terms annotation, linkage gene blocks and pathways to further identify pleiotropic genes and gene sets in the overlapping intervals of growth-related and immunity-related QTLs. In total, 55 non-redundant QTL overlapping intervals were identified, 1893 growth-related genes and 713 immunity-related genes were further classified into overlapping intervals and 405 pleiotropic genes shared by the two gene sets were determined. In addition, 19 pleiotropic gene linkage blocks and 67 pathways related to immunity and growth traits were discovered. A total of 343 growth-related genes and 144 immunity-related genes involved in pleiotropic pathways were also identified, respectively. We also sequenced and genotyped 284 individuals from Chinese Meishan pigs and European pigs and mapped the single nucleotide polymorphisms (SNPs) to the pleiotropic genes and gene sets that we identified. A total of 971 high-confidence SNPs were mapped to the pleiotropic genes and gene sets that we identified, and among them 743 SNPs were statistically significant in allele frequency between Meishan and European pigs. This study explores the relationship between growth and immunity traits from the view of QTL overlapping intervals and can be generalized to explore the relationships between other traits.     

Indel‐seq: a fast‐forward genetics approach for identification of trait‐associated putative candidate genomic regions and its application in pigeonpea (Cajanus cajan)

Identification of candidate genomic regions associated with target traits using conventional mapping methods is challenging and time‐consuming. In recent years, a number of single nucleotide polymorphism (SNP)‐based mapping approaches have been developed and used for identification of candidate/putative genomic regions. However, in the majority of these studies, insertion–deletion (Indel) were largely ignored. For efficient use of Indels in mapping target traits, we propose Indel‐seq approach, which is a combination of whole‐genome resequencing (WGRS) and bulked segregant analysis (BSA) and relies on the Indel frequencies in extreme bulks. Deployment of Indel‐seq approach for identification of candidate genomic regions associated with fusarium wilt (FW) and sterility mosaic disease (SMD) resistance in pigeonpea has identified 16 Indels affecting 26 putative candidate genes. Of these 26 affected putative candidate genes, 24 genes showed effect in the upstream/downstream of the genic region and two genes showed effect in the genes. Validation of these 16 candidate Indels in other FW‐ and SMD‐resistant and FW‐ and SMD‐susceptible genotypes revealed a significant association of five Indels (three for FW and two for SMD resistance). Comparative analysis of Indel‐seq with other genetic mapping approaches highlighted the importance of the approach in identification of significant genomic regions associated with target traits. Therefore, the Indel‐seq approach can be used for quick and precise identification of candidate genomic regions for any target traits in any crop species.     

Genome‐wide detection of genetic loci and candidate genes for teat number and body conformation traits at birth in Chinese Sushan pigs

Body conformation at birth and teat number are economically important traits in the pig industry, as these traits are usually explored to evaluate the growth and reproductive potential of piglets. To detect genetic loci and candidate genes for these traits, we performed a GWAS on 269 pigs from a recently developed Chinese breed (Sushan) using 38  128 informative SNPs on the Affymetrix Porcine SNP 55K Array. In total, we detected one genome‐wide significant (P = 1.31e‐6) SNP for teat number on chromosome X and 15 chromosome‐wide significant SNPs for teat number, body weight, body length, chest circumference and cannon circumference at birth on chromosomes 1, 3, 4, 6, 7, 9, 10, 13, 14, 15, 17 and 18. The most significant SNP had an additive effect of 0.74 × total teat number, explaining 20% of phenotypic variance. Five significant SNPs resided in the previously reported quantitative trait loci for these traits and seven significant SNPs had a pleiotropic effect on multiple traits. Intriguingly, 12 of the genes nearest to the significant SNPs are functionally related to body conformation and teat number traits, including SPRED2, MKX, TMSB4X and ESR1. GO analysis revealed that candidate genes proximal to the significant SNPs were enriched in the G‐protein coupled receptor and steroid hormone‐mediated signaling pathway. Our findings shed light on the genetic basis of the measured traits and provide molecular markers especially for the genetic improvement of teat number in Sushan and related pigs.     

Quantitative genetic estimates of growth-related traits in the common carp (Cyprinus carpio L.): A review

The common carp (Cyrpinus carpio L.) is the oldest cultured and the most domesticated fish species, as well as one of the most important freshwater fishes in the world. However, scientific studies on evaluating the growth-related quantitative traits in this fish are limited. Heritability, the most important parameter in selective breeding programs, was extensively studied for the growth-related traits. The values varied widely among the experiments and methods used because of the existence of common environmental, dominance and maternal effects. However, correlations in phenotypic and genetic levels first evaluated several years ago were limited. On the other hand, heterosis was widely reported and easily obtained for growth-related traits in the common carp. Meanwhile, genotype environment interaction and prediction of breeding values have been studied recently, and are very important in conducting selective breeding programs. The developmental quantitative genetics of growth-related traits was first analyzed in the common carp for reasonable selection during ontogeny. It is expected that genetic improvement will be achieved by carrying out direct selective breeding in the common carp.     

Quantitative genetic estimates of growth-related traits in the common carp (Cyprinus carpio L.): A review

The common carp (Cyrpinus carpio L.) is the oldest cultured and the most domesticated fish species, as well as one of the most important freshwater fishes in the world. However, scientific studies on evaluating the growth-related quantitative traits in this fish are limited. Heritability, the most important parameter in selective breeding programs, was extensively studied for the growth-related traits. The values varied widely among the experiments and methods used because of the existence of common environmental, dominance and maternal effects. However, correlations in phenotypic and genetic levels first evaluated several years ago were limited. On the other hand, heterosis was widely reported and easily obtained for growth-related traits in the common carp. Meanwhile, genotype environment interaction and predic-tion of breeding values have been studied recently, and are very important in conducting selective breeding programs. The developmental quantitative genetics of growth-related traits was first analyzed in the common carp for reasonable selection during ontogeny. It is expected that genetic improvement will be achieved by carrying out direct selective breeding in the common carp.     

羊重要性状全基因组关联分析研究进展

全基因组关联分析(genome-wide association study, GWAS)是一种复杂性状功能基因鉴定的分析策略,已成为挖掘畜禽重要经济性状候选基因的重要手段。随着绵羊和山羊基因组完成和公布,以及不同密度的SNP (single nucleotide polymorphism)芯片的推出并进行商业化推广,不仅大大丰富了羊标记辅助选择可利用的分子标记,而且还为开展重要性状的分子机理的探索提供了重要技术支撑。本文主要针对羊角、羊毛、羊奶、生长发育、肉质、繁殖和疾病等重要性状的GWAS研究所用的群体、主要研究方法和研究结果进行了综述,并对GWAS方法研究现状进行了归纳,以期为进一步利用GWAS进行羊的各种性状的遗传基础研究提供参考。     

Robust Gene-Gene Interaction Analysis in Genome Wide Association Studies

Genome-wide association studies (GWAS) have successfully discovered hundreds of associations between genetic variants and complex traits. Most GWAS have focused on the identification of single variants. It has been shown that most of the variants that were discovered by GWAS could only partially explain disease heritability. The explanation for this missing heritability is generally believed to be gene-gene (GG) or gene-environment (GE) interactions and other structural variants. Generalized multifactor dimensionality reduction (GMDR) has been proven to be reasonably powerful in detecting GG and GE interactions; however, its performance has been found to decline when outlying quantitative traits are present. This paper proposes a robust GMDR estimation method (based on the L-estimator and M-estimator estimation methods) in an attempt to reduce the effects caused by outlying traits. A comparison of robust GMDR with the original MDR based on simulation studies showed the former method to outperform the latter. The performance of robust GMDR is illustrated through a real GWA example consisting of 8,577 samples from the Korean population using the Homeostasis Model Assessment of Insulin Resistance (HOMA-IR) level as a phenotype. Robust GMDR identified the KCNH1 gene to have strong interaction effects with other genes on the function of insulin secretion.     

Validation of growth-related quantitative trait loci markers in different Exopalaemon carinicauda families for marker-assisted selection

We detected growth-related QTL and associated markers from the backcross population of Exopalaemon carinicauda in the previous study. Based on our previous study, the 47 SNP markers associated with candidate growth trait QTL were selected to analyze the association between these markers and body weight (BW), body length and abdominal segment length traits in four different populations including wild population, a full-sib family, a half-sib family and a backcross population for evaluating their potential application of marker-assisted selection in E. carinicauda. The general linear model (GLM) and mixed linear model were applied and the associations between SNP loci and three growth-related traits verified. The results showed that the Marker79268 and Marker100644 were significantly associated with the BW trait in more than three populations by the GLM method. The Marker100644 was significantly associated with BW in the full-sib family, half-sib family and backcross populations by the GLM and mixed linear model methods. Our findings will provide useful SNP markers to go forward to improve growth performance through more refined marker-assisted selection in E. carinicauda.     

CRISPR/Cas9 genome editing demonstrates functionality of the autoimmunity-associated SNP rs12946510

Genome-wide association studies (GWAS) map genetic associations of complex traits with precision limited to a linkage disequilibrium group. To translate GWAS results into new understanding of disease mechanisms, individual causative polymorphisms and their target genes should be identified. CRISPR/Cas9 genome editing can be used to create isogenic cell lines bearing alternative genotypes of candidate single-nucleotide polymorphisms to test their causality and to reveal gene targets. An intergenic polymorphism rs12946510 is associated with multiple sclerosis, inflammatory bowel disease and asthma. We created sublines of the T-helper cell line bearing alternative genotypes of rs12946510 and showed that its risk (“T”) allele is associated with lower expression of IKZF3 and ORMDL3 genes and reduced cell activation. Our editing procedure can become an effective tool for discovering new genes involved in pathogenesis of complex diseases.