Genome-wide scan for runs of homozygosity in the composite Montana Tropical® beef cattle

The aim of this study was to assess the distribution of runs of homozygosity (ROH) and autozygosity islands in the composite Montana Tropical^® beef cattle to explore hotspot regions which could better characterize the different biological types within the composite breed. Montana animals (n = 1,436) were genotyped with the GGP-LD BeadChip (~30,000 markers). ROH was identified in every individual using the plink v1.90 software. Medium and long ROH prevailed in the genome, which accounted for approximately 74% of all ROH detected. On an average, 2.0% of the genome was within ROH, agreeing with the pedigree-based inbreeding coefficient. The Montana cattle with a higher proportion of productive breed types showed the highest number of autozygosity islands (n = 17), followed by those with a higher proportion of breeds adapted to tropical environments (n = 15). Enriched terms (p < .05) associated with the immune and inflammatory response, homeostasis, reproduction, mineral absorption, and lipid metabolism were described within the autozygosity islands. In this regard, over-represented GO terms and KEGG pathways described in this population may play a key role in providing information to explore the genetic and biological mechanisms together with the genomic regions underlying each biological type that favoured their optimal performance ability in tropical and subtropical regions.     

Relationship between inbreeding and milk production traits in two Italian dairy sheep breeds

The effects of inbreeding in livestock species breeds have been well documented and they have a negative impact on profitability. The objective of this study was to evaluate the levels of inbreeding in Sarda (SAR, n = 785) and Valle del Belice (VdB, n = 473) dairy sheep breeds and their impact on milk production traits. Two inbreeding coefficients (F) were estimated: using pedigree (F_PED), or runs of homozygosity (ROH; F_ROH) at different minimum ROH lengths and different ROH classes. After the quality control, 38,779 single nucleotide polymorphisms remained for further analyses. A mixed-linear model was used to evaluate the impact of inbreeding coefficients on production traits within each breed. VdB showed higher inbreeding coefficients compared to SAR, with both breeds showing lower estimates as the minimum ROH length increased. Significant inbreeding depression was found only for milk yield, with a loss of around 7 g/day (for SAR) and 9 g/day (VdB) for a 1% increase of F_ROH. The present study confirms how the use of genomic information can be used to manage intra-breed diversity and to calculate the effects of inbreeding on phenotypic traits.     

Genome‐wide analysis of homozygosity regions in european simmental bulls

The study of Runs of Homozygosity (ROH) is a useful approach for the characterization of the genome of livestock populations. Due to their high relationship with autozygosity, ROH allow to make inference about population genetic history, to estimate the level of inbreeding, to assess within breed heterogeneity and to detect the footprints of selection on livestock genomes. Aim of this study was to investigate the distribution of runs of homozygosity in bulls belonging to five European Simmental populations and to assess the relationship between three production traits (milk yield, fat and protein contents) and autozygosity. ROH count, distribution and ROH‐based coefficient of inbreeding (F_ROH) were calculated for 3,845 Simmental bulls of five different European countries: Austria (AT), Switzerland (CH), Czech Republic (CZ), Germany (DE) and Italy (IT). Average values of ROH number per animal, and total genome length covered by ROH were 77.8 ± 20.7 and 205 ± 74.4 Mb, respectively. Bulls from AT, DE and IT exhibited similar ROH characteristics. Swiss animals showed the highest (12.6%), while CZ the lowest (4.6%) F_ROH coefficient. The relationship between ROH occurrence and milk production traits was investigated through a genome‐wide ROH‐traits association analysis (GWRA). A total of 34 regions previously associated with milk traits (yield and/or composition) were identified by GWRA. Results of the present research highlight a mixed genetic background in the 5 European Simmental populations, with the possible presence of three subgroups. Moreover, a strong relationship between autozygosity and production traits has been detected.     

Estimates of autozygosity derived from runs of homozygosity: empirical evidence from selected cattle populations

Using genome‐wide SNP data, we calculated genomic inbreeding coefficients (F_ROH > _1 Mb, F_{ROH > 2 Mb}, F_{ROH > 8 Mb} and F_{ROH > 16 Mb}) derived from runs of homozygosity (ROH) of different lengths (>1, >2, >8 and > 16 Mb) as well as from levels of homozygosity (F_HOM). We compared these values of inbreeding coefficients with those calculated from pedigrees (F_PED) of 1422 bulls comprising Brown Swiss (304), Fleckvieh (502), Norwegian Red (499) and Tyrol Grey (117) cattle breeds. For all four breeds, population inbreeding levels estimated by the genomic inbreeding coefficients F_{ROH > 8 Mb} and F_{ROH > 16 Mb} were similar to the levels estimated from pedigrees. The lowest values were obtained for Fleckvieh (F_PED = 0.014, F_{ROH > 8 Mb} = 0.019 and F_{ROH > 16 Mb} = 0.008); the highest, for Brown Swiss (F_PED = 0.048, F_{ROH > 8 Mb} = 0.074 and F_{ROH > 16 Mb} = 0.037). In contrast, inbreeding estimates based on the genomic coefficients F_{ROH > 1 Mb} and F_{ROH > 2 Mb} were considerably higher than pedigree‐derived estimates. Standard deviations of genomic inbreeding coefficients were, on average, 1.3–1.7‐fold higher than those obtained from pedigrees. Pearson correlations between genomic and pedigree inbreeding coefficients ranged from 0.50 to 0.62 in Norwegian Red (lowest correlations) and from 0.64 to 0.72 in Tyrol Grey (highest correlations). We conclude that the proportion of the genome present in ROH provides a good indication of inbreeding levels and that analysis based on ROH length can indicate the relative amounts of autozygosity due to recent and remote ancestors.     

Genome‐wide association studies reveal additional related loci for fatty acid composition in a Duroc pig multigenerational population

The aim of the present study was to detect quantitative trait loci affecting fatty acid composition in back fat and intramuscular fat in a Duroc pig population comprising seventh‐generation pedigrees using genome‐wide association studies (GWAS). In total, 305 animals were genotyped using single nucleotide polymorphisms (SNPs) array and five selected SNPs from regions containing known candidate genes related to fatty acid synthesis or metabolism. In total, 24 genome‐wide significant SNP regions were detected in 12 traits, and 76 genome‐wide suggestive SNP regions were detected in 33 traits. The Sus scrofa chromosome (SSC) 7 at 10.3 Mb was significantly associated with C17:0 in intramuscular fat, while the SSC9 at 13.6 Mb was significantly associated with C14:0 in intramuscular fat. The SSC12 at 1.0 Mb was significantly associated with C14:0 in back fat and the SSC14 at 121.0 Mb was significantly associated with C18:0 in intramuscular fat. These regions not only replicated previously reported loci containing some candidate genes involved in fatty acid composition (fatty acid synthase and stearoyl‐CoA desaturase) but also included several additional related loci.     

Combined approaches to identify genomic regions involved in phenotypic differentiation between low divergent breeds: Application in Sardinian sheep populations

Selective breeding has led to modifications in the genome of many livestock breeds. In this study, we identified the genomic regions that may explain some of the phenotypic differences between two closely related breeds from Sardinia. A total of 44 animals, 20 Sardinian Ancestral Black (SAB) and 24 Sardinian White (SW), were genotyped using the Illumina Ovine 50K array. A total of 68, 38 and 15 significant markers were identified using the case–control genome‐wide association study (GWAS), the Bayesian population differentiation analysis (F_ST) and the Rsb metric, respectively. Comparisons among the approaches revealed a total of 22 overlapping markers between GWAS and F_ST and one marker between GWAS and Rsb. Three markers detected by Rsb were also located near (<2 Mb) to highly significant regions identified by GWAS and F_ST analyses. Moreover, one candidate marker identified by GWAS and F_ST approaches was located in a run of homozygosity island that was shared by both breeds. We identified several genes involved in many phenotypic differences (such as stature and growth, reproduction, ear size, coat colour, behaviour) between the two analysed breeds. This study shows that combining several genome‐wide approaches could improve discovery of regions involved in the variability of breeding traits and responsible for the phenotypic diversity even between closely related breeds. Overall, the combination of such genome‐wide methods can be extended to other livestock breeds that share between them a similar genetic background, to understand the process that shapes the patterns of genetic variability between closely related populations.     

肃南牦牛群体遗传结构、选择信号分析和ROH检测

【目的】试验旨在对肃南牦牛进行群体遗传结构、选择信号分析和连续纯合片段(ROH)检测,挖掘肃南牦牛的种质特性相关基因。【方法】选择肃南牦牛、巴州牦牛、斯布牦牛、九龙牦牛和天祝白牦牛5个牦牛品种共计48头牦牛进行全基因组重测序,采用基因组变异检测流程(GATK)获得高质量的单核苷酸多态性(SNP)标记进行下游分析;对5个牦牛品种的基因型数据进行主成分分析、祖先成分分析和系统进化树分析以确定其群体结构;对5个牦牛群体进行ROH检测以及近交系数计算;采用综合单倍型评分(iHS)方法筛选共有高频ROH区域内受选择位点。【结果】5个牦牛品种共鉴定出15 092 883个SNPs,主要分布于内含子区域。亲缘关系计算结果显示,所有个体均不存在三代以内亲缘关系,满足后续分析要求。主成分分析表明,5个牦牛品种可以显著地分为5个类群,其中肃南牦牛群体内遗传变异较小。群体结构分析显示,肃南牦牛包含其他4种牦牛祖先成分,其中天祝白牦牛祖先成分所占的比例最大。ROH分析显示,5个牦牛品种共检测到8 426个ROHs片段,其中高频ROH区域421个。相比于其他4个牦牛品种,肃南牦牛的ROH片段总长度和数目最多,具有较高的连锁程度,其近交系数远大于其他牦牛群体。在肃南牦牛高频ROH区域上注释得到465个候选基因,主要富集到与机体发育、大脑形态形成、脂肪氧化相关的通路,包括胚胎骨骼系统形态发生、前后模式规范、甲状腺发育通路和Hippo通路,其中与胚胎发育及组织分化过程相关的基因有同源框A3(HOXA3)、HOXA5、HOXD3,与肉品质相关的基因有花生四烯酸12B (ALOX12B)、ALOX15B、ALOXE3,与中脑发育相关的基因为成纤维细胞生长因子8(FGF8)。在7号染色体1个高频ROH区域中(Chr7:12 661 870-13 045 935)鉴定到3个共享基因(核内不均一性核糖核蛋白K (HNRNPK)、驱动蛋白27(KIF27)、G激酶锚定蛋白1(GKAP1))与牦牛共有的高原适应性有关。对共有高频ROH区域内的位点进行iHS分析,鉴定到的受选择基因主要与抗病性、内质网分泌蛋白加工及细胞周期调节相关,包括N-α乙酰转移酶25(NAA25)、内质网分子伴侣29(ERP29)、跨膜蛋白16(TMEM116)、TRAF型锌指结构域蛋白1(TRAFD1)、HECT结构域E3泛素连接酶4(HECTD4)基因。【结论】本研究从全基因组水平系统评估肃南牦牛的遗传多样性和遗传背景,鉴定了肃南牦牛ROH区域内与表型相关的基因,并重点挖掘了与其他牦牛品种共享高频ROH区域内的受选择基因,为肃南牦牛种质资源开发、利用提供了重要理论依据。     

Genome-wide identification of runs of homozygosity islands in the Gyr breed (Bos indicus)

Runs of homozygosity (ROH) are contiguous homozygous regions of the genome. These regions can be used to identify genes associated with traits of economic interest, as well as inbreeding levels. The aim of the present study was to analyse the length and distribution of ROH islands in Gyr cattle and to identify genes within these regions. A population of 173 animals selected for beef production and a population of 291 animals selected for dairy production were used. Differences in the number of short ROH (ROH_{1-2 Mb}) were observed between the two populations, while the number of long ROH (ROH_{>16 Mb}) was similar. ROH islands with the highest incidences (>0.50) overlapped in several segments of the genome in the two populations. The genes identified were associated with milk production, growth, reproduction, immune response and resistance traits. Our results contribute to the understanding of how selection can shape the distribution of ROH and ROH islands within the same breed when animals are selected for different purposes such as dairy or beef production.     

Runs of homozygosity and analysis of inbreeding depression

Pedigree information was traditionally used to assess inbreeding. The availability of high-density marker panels provides an alternative to assess inbreeding, particularly in the presence of incomplete and error-prone pedigrees. Assessment of autozygosity across chromosomal segments using runs of homozygosity (ROH) has emerged as a valuable tool to estimate inbreeding due to its general flexibility and ability to quantify the chromosomal contribution to genome-wide inbreeding. Unfortunately, the identification of ROH segments is sensitive to the parameters used during the search process. These parameters are heuristically set, leading to significant variation in the results. The minimum length required to identify an ROH segment has major effects on the estimation of inbreeding and inbreeding depression, yet it is arbitrarily set. To overcome this limitation, a search algorithm to approximate mutation enrichment was developed to determine the minimum length of ROH segments. It consists of finding genome segments with significant effect differences in trait means between animals with high and low burdens of autozygous intervals with a specific length. The minimum length could be determined heuristically as the smallest interval at which a significant signal is detected. The proposed method was tested in an inbred Hereford cattle population genotyped for 30,220 SNPs. Phenotypes recorded for six traits were used for the approximation of mutation loads. The estimated minimum length was around 1 Mb for yearling weight (YW) and average daily gain (ADG) and 4 Mb for birth weight and weaning weight. These trait-specific thresholds estimated using the proposed method could be attributed to a trait-dependent effect of homozygosity. The detection of significant inbreeding effects was well aligned with the estimated thresholds, especially for YW and ADG. Although highly deleterious alleles are expected to be more frequent in recent inbreeding (long ROH), short ROH segments (<5 Mb) could contain a large number of less deleterious mutations with substantial joint effects on some traits (YW and ADG). Our results highlight the importance of accurate estimation of the ROH-based inbreeding and the necessity to consider a trait-specific minimum length threshold for the identification of ROH segments in inbreeding depression analyses. These thresholds could be determined using the proposed method provided the availability of phenotypic information.     

Detecting selection signatures on the X chromosome of the Chinese Debao pony

The X chromosome shows a special interaction between demographic factors and genetic variation, and the analysis of X‐linked genomic variation can therefore provide insights into the unique effects of demography and selection on the horse genome that cannot be readily detected by autosomal markers. Debao (DB) ponies have experienced intense selective pressure for the development of their small stature (<106 cm at adult height). To identify selective sweeps on the X chromosome of the DB pony, we performed a genome‐wide scan of three Chinese horse breeds using an Equine SNP70 BeadChip. Using Yili and Mongolian horses (>134 cm at adult height) as reference groups, both F_ST and XP‐EHH revealed that five regions on the X chromosome were under strong selection, resulting in 95 overlapping genes. Seven of these genes, SMS, PHEX, ACSL4, CHRDL1, CACNA1F, DKC1 and CDKL5, are involved in bone development, growth hormone secretion and fat deposition. The region showing the strongest selection pressure was located at the position of 86.6–87.5 Mb. The subsequent genome‐wide association analysis of the adult height of three Chinese horse breeds detected the two most significant SNPs in the same region, and these two SNPs overlapped with the gene CHRDL1. As a member of the bone morphogenetic protein (BMP) superfamily, CHRDL1 antagonizes the function of BMP4 and plays an important role in embryonic bone formation and cartilage generation. Our results provide new insights into the X‐linked selection in Chinese Debao pony.     

Genome‐wide association study for egg production and quality in layer chickens

Discovery of genes with large effects on economically important traits has for many years been of interest to breeders. The development of SNP panels which cover the whole genome with high density and, more importantly, that can be genotyped on large numbers of individuals at relatively low cost, has opened new opportunities for genome‐wide association studies (GWAS). The objective of this study was to find genomic regions associated with egg production and quality traits in layers using analysis methods developed for the purpose of whole genome prediction. Genotypes on over 4500 birds and phenotypes on over 13 000 hens from eight generations of a brown egg layer line were used. Birds were genotyped with a custom 42K Illumina SNP chip. Recorded traits included two egg production and 11 egg quality traits (puncture score, albumen height, yolk weight and shell colour) at early and late stages of production, as well as body weight and age at first egg. Egg weight was previously analysed by Wolc et al. ( 2012 ). The Bayesian whole genome prediction model – BayesB (Meuwissen et al. 2001 ) was used to locate 1 Mb regions that were most strongly associated with each trait. The posterior probability of a 1 Mb window contributing to genetic variation was used as the criterion for suggesting the presence of a quantitative trait locus (QTL) in that window. Depending upon the trait, from 1 to 7 significant (posterior probability >0.9) 1 Mb regions were found. The largest QTL, a region explaining 32% of genetic variance, was found on chr4 at 78 Mb for body weight but had pleiotropic effects on other traits. For the other traits, the largest effects were much smaller, explaining <7% of genetic variance, with regions on chromosomes 2, 12 and 17 explaining above 5% of genetic variance for albumen height, shell colour and egg production, respectively. In total, 45 of 1043 1 Mb windows were estimated to have a non‐zero effect with posterior probability > 0.9 for one or more traits.     

Signatures of selection in Charolais beef cattle identified by genome‐wide analysis

Charolais cattle are one of the most important breeds for meat production worldwide; in México, its selection is mainly made by live weight traits. One strategy for mapping important genomic regions that might influence productive traits is the identification of signatures of selection. This type of genomic features contains loci with extended linkage disequilibrium (LD) and homozygosity patterns that are commonly associated with sites of quantitative trait locus (QTL). Therefore, the objective of this study was to identify the signatures of selection in Charolais cattle genotyped with the GeneSeek Genomic Profiler Bovine HD panel consisting of 77 K single nucleotide polymorphisms (SNPs). A total 61,311 SNPs and 819 samples were used for the analysis. Identification of signatures of selection was carried out using the integrated haplotype score (iHS) methodology implemented in the rehh R package. The top ten SNPs with the highest piHS values were located on BTA 4, 5, 6 and 14. By identifying markers in LD with top ten SNPs, the candidate regions defined were mapped to 52.8–59.3 Mb on BTA 4; 67.5–69.3 on BTA 5; 39.5–41.0 Mb on BTA 6; and 26.4–29.6 Mb on BTA 14. The comparison of these candidate regions with the bovine QTLdb effectively confirmed the association (p < 0.05) with QTL related to growth traits and other important productive traits. The genomic regions identified in this study indicated selection for growth traits on the Charolais population via the conservation of haplotypes on various chromosomes. These genomic regions and their associated genes could serve as the basis for haplotype association studies and for the identification of causal genes related to growth traits.     

Investigations on the pattern of linkage disequilibrium and selection signatures in the genomes of German Piétrain pigs

The aim of this study was to study the population structure, to characterize the LD structure and to define core regions based on low recombination rates among SNP pairs in the genome of Piétrain pigs using data from the PorcineSNP60 BeadChip. This breed is a European sire line and was strongly selected for lean meat content during the last decades. The data were used to map signatures of selection using the REHH test. In the first step, selection signatures were searched genome‐wide using only core haplotypes having a frequency above 0.25. In the second step, the results from the selection signature analysis were matched with the results from the recently conducted genome‐wide association study for economical relevant traits to investigate putative overlaps of chromosomal regions. A small subdivision of the population with regard to the geographical origin of the individuals was observed. The extent of LD was determined genome‐wide using r² values for SNP pairs with a distance ≤5 Mb and was on average 0.34. This comparable low r² value indicates a high genetic diversity in the Piétrain population. Six REHH values having a p‐value < 0.001 were genome‐wide detected. These were located on SSC1, 2, 6 and 17. Three positional candidate genes with potential biological roles were suggested, called LOC100626459, LOC100626014 and MIR1. The results imply that for genome‐wide analysis especially in this population, a higher marker density and higher sample sizes are required. For a number of nine SNPs, which were successfully annotated to core regions, the REHH test was applied. However, no selection signatures were found for those regions (p‐value < 0.1).     

Characteristic features of porcine endogenous retroviruses in Vietnamese native pigs

We aimed to clarify the genomic characteristics of porcine endogenous retroviruses (PERVs) in Vietnamese native pig (VnP) breeds. First, we investigated genetic polymorphisms in β‐ and γ‐like PERVs, and we then measured the copy numbers of infectious γ‐like PERVs (PERV‐A, B, and C). We purified genomic DNA from 15 VnP breeds from 12 regions all over the country and three Western pig breeds as controls, and investigated genetic polymorphisms in all known PERVs, including the beta (β)1–4 and gamma (γ)1–5 groups. PERVs of β1, β2, β3, and γ4 were highly polymorphic with VnP‐specific haplotypes. We did not identify genetic polymorphisms in β4, γ1, or γ2 PERVs. We then applied a real‐time polymerase chain reaction–based method to estimate copy numbers of the gag, pol, and env genes of γ1 PERVs (defined as A, B, and C). VnP breeds showed significantly lower copy number of the PERV genes compared with the Western pig breeds (on average, 16.2 and 35.7 copies, respectively, p < .05). Two VnP breeds showed significantly higher copy number compared with the other VnPs (p < .05). Our results elucidated that VnPs have specific haplotypes and a low copy number of PERV genes.     

Detection and analysis of genome-wide copy number variation in the pig genome using an 80 K SNP Beadchip

Copy number variation (CNV) is an important source of genetic variability in human or animal genomes and play key roles in phenotypic diversity and disease susceptibility. In the present study, we performed a genome-wide analysis for CNV detection using SNP genotyping data of 857 Large White pigs. A total of 312 CNV regions (CNVRs) were detected with the PennCNV algorithm, which covered 57.76 Mb of the pig genome and correspond to 2.36% of the genome sequence. The length of the CNVRs on autosomes ranged from 1.77 Kb to 1.76 Mb with an average of 185.11 Kb. Of these, 220 completely or partially overlapped with 1,092 annotated genes, which enriched a wide variety of biological processes. Comparisons with previously reported pig CNVR revealed 92 (29.49%) novel CNVRs. Experimentally, 80% of CNVRs selected randomly were validated by quantitative PCR (qPCR). We also performed an association analysis between some of the CNVRs and reproductive traits, with results demonstrating the potential importance of CNVR61 and CNVR283 associated with litter sizes. Notably, the GPER1 gene located in CNVR61 plays a key role in reproduction. Our study is an important complement to the CNV map in the pig genome and provides valuable information for investigating the association between genomic variation and economic traits.     

Genome‐wide association studies for two exterior traits in Chinese Dongxiang spotted pigs

Pigs show extensive variation in exterior appearance. This variation has been explored as one of the selection targets to form breeding features in the pig industry. In this study, a customized Affymetrix Axiom Pig1.4M array plate was used to conduct genome‐wide association studies (GWAS) for two exterior traits—spotted coat color (SC) and facial type (FT)—in Chinese Dongxiang spotted pigs. Two and 12 single nucleotide polymorphisms (SNPs) were identified to be associated with SC and FT at the 5% genome‐wide significant level, respectively. Of these SNPs, two associated with SC were located around the DACH1 gene on Sus Scrofa chromosome (SSC)11. Eleven SNPs associated with FT were located within a 7.68 Mb region (29.03–36.71 Mb) on SSC7. Based on the GWAS results and the biological functions of genes, we highlight EDNRB as a candidate gene for SC and HMGA1 and RPS10 as the potential genes affecting facial variation. The findings contribute to the final characterization of causative genes and mutations underlying the effects of the SSC7 and SSC11 loci and improve our understanding of the genetic basis of phenotypic variation in Chinese indigenous pigs.     

Use of high‐density SNP data to identify patterns of diversity and signatures of selection in broiler chickens

The development of broiler chickens over the last 70 years has been accompanied by large phenotypic changes, so that the resulting genomic signatures of selection should be detectable by current statistical techniques with sufficiently dense genetic markers. Using two approaches, this study analysed high‐density SNP data from a broiler chicken line to detect low‐diversity genomic regions characteristic of past selection. Seven regions with zero diversity were identified across the genome. Most of these were very small and did not contain many genes. In addition, fifteen regions were identified with diversity increasing asymptotically from a low level. These regions were larger and thus generally included more genes. Several candidate genes for broiler traits were found within these ‘regression regions’, including IGF1, GPD2 and MTNR1AI. The results suggest that the identification of zero‐diversity regions is too restrictive for characterizing regions under selection, but that regions showing patterns of diversity along the chromosome that are consistent with selective sweeps contain a number of genes that are functional candidates for involvement in broiler development. Many regions identified in this study overlap or are close to regions identified in layer chicken populations, possibly due to their shared precommercialization history or to shared selection pressures between broilers and layers.     

基于全基因组SNP标记分析中国地方鸡品种的遗传多样性和种群结构

旨在对中国地方鸡品种的遗传多样性与种群结构进行分析。本研究使用Affymetrix Axiom 600K高密度鸡基因分型芯片对来自8个品种的157只地方鸡及233只商品鸡进行基因分型,以品种作为分组来计算各分组的观测杂合度、期望杂合度、次等位基因频率、近交系数及核苷酸多样性分析地方鸡群体的遗传多样性,利用进化树、主成分分析、群体结构、MDS等方法分析鸡群体的群体结构,基于状态同源(IBS)和群体分化系数(Fst)分析种群内部与种群之间的亲缘关系,利用长纯合片段(runs of homozygosity, ROH)估算得到基于ROH的近交系数。结果表明,各群体的观测杂合度均高于期望杂合度,次等位基因频率在0.175～0.236之间,近交系数在0.018～0.205之间,核苷酸多样性在0～6×10^-4之间,进化树与主成分分析表明品种间出现了明显的群体分化,地方鸡群体与商品鸡群的MDS分析发现我国地方鸡与商业肉鸡品种的遗传距离较近;IBS遗传距离在0.092 9～0.319 9之间;各品种成对Fst分析表明,群体间呈现中高分化程度(0.09～0.22);此次分析共得到了...     

基于SNP芯片的海南猪全基因组选择信号分析

旨在通过检测海南猪全基因组上的选择信号,以挖掘与海南猪重要经济性状相关的候选基因,并解析讨论海南猪在进化驯化历史中的受选择情况.本研究利用68头海南猪的GeneSeek Genomic Profiler Procine SNP 80K芯片数据,使用整合单倍型分数(integrated haplotype score,i...     

Genetic diversity and genomic inbreeding in Japanese Black cows in the islands of Okinawa Prefecture evaluated using single-nucleotide polymorphism array

Maintaining genetic diversity and inbreeding control are important in Japanese Black cattle production, especially in remote areas such as the islands of Okinawa Prefecture. Using a single-nucleotide polymorphism (SNP) array, we evaluated the genetic diversity and genomic inbreeding in Japanese Black cows from the islands of Okinawa Prefecture and compared them to those from other locations across Japan. Linkage disequilibrium decay was slower in cows in the islands of Okinawa Prefecture. The estimated effective population size declined over time in both populations. The genomic inbreeding coefficient (F_ROH) was estimated using long stretches of consecutive homozygous SNPs (runs of homozygosity; ROH). F_ROH was higher in the cows on the islands of Okinawa Prefecture than on other locations. In total, 818 ROH fragments, including those containing NCAPG and PLAG1, which are major quantitative trait loci for carcass weight in Japanese Black cattle, were present at significantly higher frequencies in cows in the islands of Okinawa Prefecture. This suggests that the ROH fragments are under strong selection and that cows in the islands of Okinawa Prefecture have low genetic diversity and high genomic inbreeding relative to those at other locations. SNP arrays are useful tools for evaluating genetic diversity and genomic inbreeding in cattle.