基于温带和热带玉米群体全基因组F_(ST)和XP-EHH的选择信号检测

【目的】玉米起源于热带地区,经过自然和人工选择,广泛的种植于温带地区。开花是玉米生长发育的中心环节,也是热带玉米向温带环境种植的主要适应性性状。鉴定玉米在驯化过程中出现的受选择基因区段,并进一步挖掘开花候选基因,为玉米的群体改良、开花遗传机理解析提供数据支撑。【方法】首先单独分析30份温带玉米自交系和21份热带玉米自交系的单倍型数据,通过过滤高缺失和等位基因频率较低的变异位点,得到高质量的SNP(single nucleotide polymorphism)标记,利用SnpEff软件对温带和热带玉米群体的基因组多态性位点进行了功能预测。其次过滤得到同时存在于温带和热带玉米的高质量SNP标记,对温带和热带玉米的基因型数据进行主成分分析(principle component analysis,PCA)以确定其群体结构,之后利用群体分化指数(fixationindex,F_(ST))和群体间扩展单倍型纯合度(cross population extended haplotype homozygosity,XP-EHH)法分析温带和热带玉米群体间的选择信号分布情况,选择F_(ST)和XP-EHH值的top 1%为阈值,筛选得到受选择位点。通过对SNP进行功能注释得到温热带玉米群体受到选择的基因。利用agriGO工具对候选驯化基因进行功能富集分析。利用相关的生物信息学数据库对候选基因进行功能注释,进一步鉴定玉米驯化过程中的开花候选基因。【结果】通过对温热带玉米群体的高测序深度的SNP进行分析,发现热带玉米群体的SNP数目为14 123 408个,温带玉米群体的SNP数目为8 791 673个,鉴定到的SNP主要分布于基因间区。2个群体中均存在的SNP标记数目是204 752个。主成分分析表明温带和热带玉米可以显著的分为两个类群。F_(ST)选择信号的top 1%是0.3593,共鉴定到557个候选驯化基因,XP-EHH选择信号法的top 1%是3.2681,共鉴定到1 913个候选基因。鉴定到多个候选基因与玉米的开花调控密切相关,包括ZmCCT9、COL1、GRMZM2G387528。如ZmCCT9抑制开花基因ZCN8的表达,导致玉米在长日照环境下出现晚花表型,是一个重要的开花调控基因;COL1与开花促进因子FT蛋白互作,加速玉米开花以适应长日照环境;GRMZM2G387528的功能注释揭示该基因是一个光敏色素互作因子,与光周期基因ZmphyB1互作。【结论】热带玉米群体具有更高的遗传多态性,筛选到一系列参与了热带玉米和温带玉米的分化候选基因,并且重点挖掘了参与其中的玉米开花调控相关基因。

Estimation of inbreeding and identification of regions under heavy selection based on runs of homozygosity in a Large White pig population

Background: Runs of homozygosity(ROHs) are homozygous segments of the genome where the two haplotypes inherited from the parents are identical. The current availability of genotypes for a very large number of single nucleotide polymorphisms(SNPs) is leading to more accurate characterization of ROHs in the whole genome. Here,we investigated the occurrence and distribution of ROHs in 3,692 Large White pigs and compared estimates of inbreeding coefficients calculated based on ROHs(FROH), homozygosity(FHOM), genomic relationship matrix(FGRM)and pedigree(FPED). Furthermore, we identified genomic regions with high ROH frequencies and annotated their candidate genes.Results: In total, 176,182 ROHs were identified from 3,569 animals, and all individuals displayed at least one ROH longer than 1 Mb. The ROHs identified were unevenly distributed on the autosomes. The highest and lowest coverages of Sus scrofa chromosomes(SSC) by ROH were on SSC14 and SSC13, respectively. The highest pairwise correlation among the different inbreeding coefficient estimates was 0.95 between FROH_totaland FHOM, while the lowest was-0.083 between FGRMand FPED. The correlations between FPEDand FROHusing four classes of ROH lengths ranged from 0.18 to 0.37 and increased with increasing ROH length, except for ROH > 10 Mb. Twelve ROH islands were located on four chromosomes(SSC1, 4, 6 and 14). These ROH islands harboured genes associated with reproduction, muscular development, fat deposition and adaptation, such as SIRT1, MYPN, SETDB1 and PSMD4.Conclusion: FROHcan be used to accurately assess individual inbreeding levels compared to other inbreeding coefficient estimators. In the absence of pedigree records, FROHcan provide an alternative to inbreeding estimates.Our findings can be used not only to effectively increase the response to selection by appropriately managing the rate of inbreeding and minimizing the negative effects of inbreeding depression but also to help detect genomic regions with an effect on traits under selection.

Genome-wide analysis of runs of homozygosity identifies new susceptibility regions of lung cancer in Han Chinese

Runs of homozygosity （ROHs） are a class of important but poorly studied genomic variations and may be in- volved in individual susceptibility to diseases. To better understand ROH and its relationship with lung cancer, we performed a genome-wide ROH analysis of a subset of a previous genome-wide case-control study （1,473 cases and 1,962 controls） in a Han Chinese population. ROHs were classified into two classes, based on lengths, intermedi- ate and long ROils, to evaluate their association with lung cancer risk using existing genome-wide single nucleofide polymorphism （SNP） data. We found that the overall level of intermediate ROHs was significantly associated with a decreased risk of lung cancer （odds ratio = 0.63; 95% confidence interval： 0.51-0.77; P = 4.78 × 10-6 ）, while the long ROHs seemed to be a risk factor of lung cancer. We also identified one ROH region at 14q23A that was con- sistently associated with lung cancer risk in the study. These results indicated that ROHs may be a new class of variation which may be associated with lung cancer risk, and genetic variants at 14q23.1 may be involved in the development of lung cancer.

Genome-wide detection for runs of homozygosity analysis in three pig breeds from Chinese Taihu Basin and Landrace pigs by SLAF-seq data

Erhualian(E),Meishan(MS)and Mi(MI)pigs are excellent indigenous pig breeds in Chinese Taihu Basin,which have made great contributions to the genetic improvement of commercial pigs.Investigation of the genetic structure and inbreeding level of the 3 pig breeds is of great significance for the sustainable breeding of commercial pigs.The length and number of runs of homozygosity(ROH)as well as the frequency of genomes covered by ROH can be used as indicators to evaluate the level of inbreeding and the origin of the population.In this study,the ROH characteristics of E,MS,MI and Landrace(L)pigs were analyzed by SLAF-seq data,and the inbreeding coefficient based on ROH(F_(ROH))was calculated.In addition,we have identified candidate genes in the genomic regions associated with ROH.A total of10568 ROH were detected in 116 individuals of 4 pig breeds.The analysis showed that there were significant differences in genetic structure between 3 Taihu Basin pig breeds and L,and the genetic structure of E and MI was similar.The results of F_(ROH)showed that the inbreeding level of MS was the highest(0.25±0.07),while E and MI were lower than L.Compared with the other 3 pig populations,MS showed a higher frequency of long ROH(>5 Mb),indicating higher inbreeding in MS in recent generations.A large number of candidate genes related to reproductive traits are located in the genomic regions with a high frequency of ROH,and these genes are expected to be used as candidate genes in marker-assisted selection(MAS)breeding programs.Our findings can provide theoretical support for genetic conservation and genetic improvement of 3 pig breeds in Chinese Taihu Basin.

Homozygosity mapping of a consanguineous Pakistani family affected with oculocutaneous albinism to Tyrosinase gene

Dear Sir,Iam Haiba Kaul,from the Department of Biochemistry,University of Health Sciences,Lahore,Pakistan.I write to present a study of oculocutaneous albinism（OCA）in consanguineous Pakistani families.OCA is a genetic defect of melanin biosynthesis that mainly affects eyes,skin and hair.

Comment on homozygosity mapping of a consanguineous Pakistani family affected with oculocutaneous albinism to Tyrosinase gene

Dear Editor,I have carefully read the article entitled ＂Homozygosity mapping of a consanguineous Pakistani family affected with oculocutaneous albinism to Tyrosinase gene＂,published by Shakil et al in 2016 and found it very interesting for the scientific community.

基于GBS测序的湘东黑山羊的亲缘关系及近交系数

利用基因分型测序(GBS)对60只湘东黑山羊(9只公羊、51只母羊)进行基因组测序,获得湘东黑山羊染色体上的SNP数据和连续纯合子片段(ROH)数据。通过SNP数据进行基于G矩阵的基因组亲缘关系分析和NJ聚类分析,构建湘东黑山羊的群体家系结构,并通过ROH数据得到群体平均近交系数。结果表明,60只湘东黑山羊的平均测序深度为9.15X,与参考基因组比对率平均为99.59%;在SNP检测过程中,3只母羊不符合基因型数据质控标准,将其过滤后获得其他57只黑山羊29条染色体上153046个SNPs位点和1937个ROH片段;构建的湘东黑山羊8个家系,其中9只公羊和6只母羊被分为7个家系,另42只母羊被单独分类;基于ROH片段分析的湘东黑山羊群体中每个个体的近交系数均值为0.047,说明湘东黑山羊公羊在繁育过程中有较大的利用空间。

黄鳝减数分裂雌核发育及子代遗传纯合度分析

为了培育性状优良、遗传稳定的黄鳝(Monopterus albus)群体,研究探索了人工诱导黄鳝减数分裂雌核发育方法。针对养殖性状优良的深黄大斑鳝进行种质纯合,创制出3个深黄大斑鳝雌核发育群体。流式细胞术和染色体计数分析表明雌核发育黄鳝的细胞DNA含量、染色体数目均与野生型相同。性腺组织学切片观察表明雌核发育黄鳝卵巢发育正常。微卫星多样性分析表明雌核发育黄鳝群体的有效等位基因(Ne)、平均香农指数(I)、平均多态信息指数(PIC)、平均观测杂合度(Ho)及平均期望杂合度(He)等参数均极显著低于野生群体,雌核发育群体的遗传纯合度显著提高。雌核发育群体之间的遗传距离增大,遗传相似系数变小。深黄大斑鳝雌核发育群体为培育性状优良的黄鳝养殖品种提供了育种材料。

利用“黄海芯1号”55K SNP芯片评估凡纳滨对虾选育群体的遗传多样性与基因组近交水平

凡纳滨对虾的主要选育目标分为两个方面:一是培育具有较强抗病、抗逆性的“高抗系”(GK),二是培育具有快速生长特性的“快大系”(KD)。然而,国内缺少针对这两个选育群体的遗传多样性特别是基因组近交水平的调查分析研究。基于液相芯片“黄海芯1号”(55 K SNP)的基因分型数据,首次分析了GK(1064尾个体)和KD(564尾个体)选育群体的遗传结构和遗传多样性,调查了连续性纯合片段(ROH)的基因组分布特征,并重点评估了两个群体的基因组近交水平。PCA及进化树分析表明GK及KD群体可明确分层,亲缘关系热图表明KD群体内个体间的亲缘关系比GK群体更近。GK群体包括的家系数量更多,导致其遗传多样性高于KD群体;两群体间的F_(st)为0.09,存在中等遗传分化。GK和KD群体每个ROH的平均长度分别为(1.70±0.34)Mb和(1.65±0.38)Mb,每个样本ROH的平均数量分别为1.98±1.30和2.07±1.37。GK和KD群体0.8~1.25 Mb长度的ROH占比分别为11.41%和19.17%,表明KD群体的选育历史比GK群体更长。两个群体>2.25 Mb长度的ROH片段占比分别为10.26和9.74%,表明两个群体短期内未发生近亲交配。七种基因组近交系数评估结果表明,KD群体的近交水平高于GK群体。不依赖基础群体等位基因频率的F_(ROH)和F_(HOM)方法可准确地评价育种群体的真实近交水平,而F_(VR1)、F_(YA1)和F_(LH1)等依赖基础群体等位基因频率的方法可以用来比较群体及个体间的相对近交水平。上述结果为准确地评估育种群体的近交水平和优化育种方案提供了重要参考依据。

基于简化基因组测序评估陇南山羊群体遗传多样性和群体结构

该试验旨在利用简化基因组测序(specific-locus amplified fragment sequencing, SLAF-seq)分析陇南山羊群体遗传多样和群体结构,为陇南山羊后续保种计划的制订提供依据。该研究对50只陇南山羊(公、母各25只)进行全基因组SNP检测;计算观测杂合度(Ho)、期望杂合度(He)、多态信息含量(PIC)、香浓维纳指数(SHI)、基因多样性指数(Nei)及次要等位基因频率(MAF)等6个指标进行遗传多样性评估。结果表明,50只陇南山羊共检测到655514个SNPs位点。陇南山羊的Ho、 He、 PIC、 SHI、 Nei及MAF指标值分别为0.193、0.286、 0.236、 0.449、 0.290、 0.200,说明该群体遗传多样性较低。陇南山羊群体LD值较低,衰退速度较快,说明该群体并未受到过强烈的人工选择压力。此外,PCA与系统发育树结果均表明陇南山羊群体内部出现分化,可分为2个大的亚群。陇南山羊群体平均IBS遗传距离为0.7391,结合亲缘关系G矩阵热图,表明陇南山羊群体大部分个体亲缘关系较远。50只陇南山羊个体共检测到47206个ROH,平均ROH长度37.86 Mb,基于ROH的近交系数FROH范围为0.0535~0.2574,平均FROH值为0.1472,说明陇南山羊群体存在近交风险。可见,陇南山羊内部存在分化,可分为2个大的亚群,陇南山羊群体遗传多样性较低,部分个体间存在较大的近交风险,可能需要采取保种措施保护陇南山羊遗传资源。

基于Illumina Bovine SNP 50K芯片对比利时蓝牛的遗传规律解析

为了研究大型肉牛比利时蓝牛生长发育的遗传规律,筛选优异基因,试验基于Illumina Bovine SNP 50K芯片数据,采用PLINK软件对270头比利时蓝牛常染色体数据进行基因组长纯合片段(ROH)检测并基于选择信号分析,通过核苷酸多态性检测取前5%的单核苷酸多态性(single nucleotide polymorphism,SNP)位点,基于牛参考基因组(ARS-UCD1.2)对结果SNPs进行基因注释,对候选基因进行GO功能注释与KEGG信号通路富集分析,并计算染色体上ROH长度占基因组总长度的比例(FROH)。结果表明:在全部270个个体数据中共检测出1893个ROH片段,平均长度13.2311 Mb,平均FROH为0.0392;得到与生长发育相性状相关的基因有NEB、TET2、NEK11、NCKAP1、MYH15、EIF4A2、bta-miR-1248-1、DCAF8、PRORP、DOCK3、SYT15、MYEF2、ZDHHC13,与公牛生育能力相关的基因有CFAP61、DNAL1、BAG1。说明通过对比利时蓝牛生长发育性状相关分子标记的解析可以为比利时蓝牛遗传改良提供理论指导。

Conservation genomics provides insights into genetic resilience and adaptation of the endangered Chinese hazelnut, Corylus chinensis

Global climate change has increased concerns regarding biodiversity loss.However,many key conservation issues still required further research,including demographic history,deleterious mutation load,adaptive evolution,and putative introgression.Here we generated the first chromosome-level genome of the endangered Chinese hazelnut,Corylus chinensis,and compared the genomic signatures with its sympatric widespread C.kwechowensis-C yunnanensis complex.We found large genome rearrangements across all Corylus species and identified species-specific expanded gene families that may be involved in adaptation.Population genomics revealed that both C.chinensis and the C.kwechowensis-C.yunnanensis complex had diverged into two genetic lineages,forming a consistent pattern of southwestern-northern differentiation.Population size of the narrow southwestern lineages of both species have decreased continuously since the late Miocene,whereas the widespread northern lineages have remained stable(C.chinensis) or have even recovered from population bottlenecks(C.kwechowensis-C.yunnanensis complex) during the Quaternary.Compared with C.kwechowensis-C. yunnanensis complex,C.chinensis showed significantly lower genomic diversity and higher inbreeding level.However,C.chinensis carried significantly fewer deleterious mutations than C.kwechowensis-C. yunnanensis complex,as more effective purging selection reduced the accumulation of homozygous variants.We also detected signals of positive selection and adaptive introgression in different lineages,which facilitated the accumulation of favorable variants and formation of local adaptation.Hence,both types of selection and exogenous introgression could have mitigated inbreeding and facilitated survival and persistence of C.chinensis.Overall,our study provides critical insights into lineage differentiation,local adaptation,and the potential for future recovery of endangered trees.

基于全基因组重测序检测中国地方猪的体型选择信号

[目的]利用全基因组选择信号检测方法探索中国地方猪基因组中与体型相关的候选基因,分析中国地方猪在进化和驯化过程中的受选择情况。[方法]基于310头地方猪的全基因组重测序数据,利用跨群体扩展单倍型纯合(cross population extended haplotype homozygosity,XP-EHH)和固定分化指数F ST统计量(fixation index)两种方法进行全基因组选择信号检测。取两种方法各前0.1%位点的交集作为候选位点,分别向上、下游延伸200 kb作为潜在受选择区域。通过富集分析进一步探索选择信号的生物学功能。[结果]使用XP-EHH和F ST两种方法共检测到633个潜在受选择位点,获得633个潜在区域,共注释到133个基因。数量性状基因座(quantitative trait locus,QTL)富集分析发现,不同体型中国地方猪之间的选择信号与肉质和胴体性状相关;染色质状态富集结果发现,不同体型中国地方猪群体之间的差异组织为内脏、消化系统和小脑。GO功能和KEGG通路富集分析发现16个基因在6个功能条目(P<0.05,count>3)上显著富集,主要集中在生长代谢相关通路,其中ACSM 5、ACSM 4、FXR 1和FOXA 2作为候选基因主要富集于脂肪酸合成与代谢、肌肉生长发育等信号通路。[结论]本研究采用两种选择信号检测方法对不同体型的中国地方猪进行了分析,筛选到一系列候选位点和基因,重点挖掘了参与猪生长发育的候选基因,如ACSM 5、ACSM 4、FXR 1和FOXA 2。

微卫星评价牙鲆雌核发育二倍体纯合性

采用8个微卫星座位分别对牙鲆减数雌核发育二倍体家系和卵裂雌核发育二倍体家系的纯合性进行检验。卵裂雌核发育二倍体在所有检测座位全部纯合。减数雌核发育二倍体在部分座位发生纯合,但未发现在所有座位全部纯合的个体,在Poli9TUF、Poli9-8TUF、Poli11TUF、Poli13TUF、Poli23TUF、Poli30TUF、Poli123TUF和Poli130TUF座位,杂合子比例分别为1·0000、1·0000、0·1944、0·9459、0·8611、1·0000、0·7778和0·8000,平均杂合子比例为0·8224。由此表明,牙鲆除了Poli11TUF外,在其余7个座位均具有很高的重组率。研究结果显示,牙鲆卵裂雌核发育二倍体一代即可形成纯合子;而减数雌核发育二倍体由于具有较高的重组率,使其与母本的遗传同质性较高。

有丝分裂雌核发育牙鲆的微卫星鉴定

通过静水压抑制卵裂方法获得多个牙鲆(Paralichthys olivaceus)有丝分裂雌核发育家系,在仔鱼阶段经历大批死亡后,将各个家系残存的个体混合饲养,24个月后获得1批存活个体。利用11对微卫星分子标记,对173尾有丝分裂雌核发育牙鲆进行分析,结果表明,其中111尾在11位点全部纯合,62尾在部分位点杂合,平均杂合子比例为0.2338。11个位点的等位基因数为4～8,平均等位基因数为5.0336,有效等位基因数为2.0421～5.1268,平均有效等位基因数为3.2815。采用NJ法对111尾纯合牙鲆的亲缘关系进行聚类分析,可分为4群。微卫星标记鉴定结果表明,通过静水压抑制卵裂法可以获得完全纯合的牙鲆,但有部分杂合个体出现。分析其原因,可能是由于某些卵子发育不同步,静水压抑制其第二极体排出而非全部抑制卵裂,导致减数分裂雌核发育所致。

相对Shannon信息量与基因变异的测量

在ｓｈａｎｎｏｎ信息量的基础上，建立了相对纯合度信息量Ｓ′Ｊ（Ａ２）、相对杂合度信息量Ｓ′Ｈ（Ａ２）和相对信息量Ｓ′（Ａ２）的概念，并赋予它们以遗传学意义。与纯合度Ｊ和杂合度Ｈ进行了理论比较，结果表明，二者在数量规律性上有很好的一致性，而各相对信息量有更好的性质。１－Ｓ′Ｈ主要反映基因的确定度，Ｓ′Ｈ（Ａ２）主要反映基因的变异度；Ｓ′（Ａ２）与Ｋ无关，它能反映群体的遗传变异程度，亦能比较不同位点间的遗传变异程度。

牙鲆连续两代减数分裂雌核发育家系的遗传特征

对牙鲆(Paralichthys olivaceus)减数分裂雌核发育二倍体(Meio-G1)再次诱导减数分裂雌核发育,获得连续两代雌核发育二倍体家系(Meio-G2),以Meio-G1亲本与1尾普通牙鲆雄鱼人工授精获得的家系作为对照组(control)。利用15个微卫星标记对3个家系进行遗传特征分析。结果显示,15个微卫星位点在Meio-G1、Meio-G2和对照组3个家系中,分别扩增到30、28、50个等位基因,平均等位基因数为2.0、1.9、3.3,平均观测杂合度(Ho)分别为0.875 3、0.774 2、0.908 3,平均纯合度分别为0.124 7、0.221 5、0.091 7。3个家系个体间的平均遗传相似系数分别为0.891 7、0.923 8、0.520 2,亲代与子代之间的平均相似系数分别为0.916 6、0.930 4、0.560 3。高重组率的Poli9-8tuf、Poli18tuf、Poli107tuf 3个位点在Meio-G1和Meio-G2中观测杂合度均为1.0,低重组率的Poli33tuf、Poli24MHFS两个位点在Meio-G1和Meio-G2中均全部纯合,观测杂合度为0。结果表明,Meio-G2的纯合度、个体间平均相似系数以及亲子之间的平均相似度均略高于Meio-G1,显著高于对照组家系,证明连续两代诱导减数分裂雌核发育,能提高鱼类纯合度和遗传相似度,具有固定母本遗传性状的作用。

大黄鱼连续两代雌核发育群体的微卫星标记分析

通过对大黄鱼(Pseudosciaena crocea)异质雌核发育一代群体(meio-G1)与二代群体(meio-G2)微卫星位点的纯合度进行分析,研究异质雌核发育对大黄鱼基因纯化的效率。结果显示:meio-G1和meio-G215个微卫星座位的平均纯合度分别为0.661和0.803,纯合位点比例最高个体分别为0.867(13/15)和0.933(14/15),两个群体内个体间的平均相似系数分别为0.5903和0.8672,最高分别达0.9286和1.0(遗传距离为0.0741和0),远高于两性交配繁殖群体(平均纯合度0.376,平均相似系数0.4687,个体间最小遗传距离0.2288);其中meio-G2群体有7个位点(46.7%)已经完全纯合固定,并与普通养殖群体产生较明显的遗传分化;表明人工诱导异质雌核发育可大大加速大黄鱼大多数基因位点的纯合,是快速建立高纯品系的有效手段。但不同位点的纯合度差异很大,部分位点在异质雌核发育后代中迅速纯合,在meio-G1中就达到很高的纯合度,而有些位点则在meio-G1和meio-G2中仍保持很高的杂合度;meio-G1和meio-G2群体中不同个体纯合位点比例差异也很大。研究培育的雌核发育群体为大黄鱼进一步选育提供了良好的遗传材料。

国家大豆区域试验品种(系)SSR位点纯合度分析

利用30对SSR引物,检测2005—2009年参加国家区域试验的1068份大豆品种（系）的位点纯合度。每年区试品种（系）的平均纯合度为94.9%~97.6%,纯合度高于85%和90%的品种（系）所占比例分别为95.0%和91.4%,纯合度低于85%的品种（系）有42份（占3.93%）,主要为北方春大豆和黄淮夏大豆。位点纯合度低于85%的参试品种（系）的产量分析结果表明,只有11份品种（系）比对照增产5%以上,20份比对照减产0.04%~13.08%;与纯合度为100%的品种（系）比较发现,位点纯合度较低的品种（系）在区试中产量也较低。建议国家区试品种（系）纯合度标准不低于90%,以保证审定品种的特征特性,为大豆新品种的持续推广利用提供科技支撑。

枫泾猪、梅山猪和沙乌头猪全基因组ROH检测及候选基因鉴定

通过检测枫泾猪、梅山猪和沙乌头猪的长纯合片段(Runs of homozygosity,ROH),鉴定其与3个猪种经济性状相关的候选基因。利用全基因组单核苷酸多态性(Single Nucleotide Polymorphisms,SNP)芯片数据对这3个猪群体进行全基因组ROH检测,统计ROH数量、长度及分布,计算基于ROH的近交系数(FROH),并在高频ROH区域鉴定候选功能基因。结果显示,从枫泾猪、梅山猪和沙乌头猪群体中分别检测到1588个、3314个和4419个ROH。枫泾猪的ROH平均长度为12.84 Mb,平均FROH为0.084;梅山猪的ROH平均长度为10.55 Mb,平均FROH为0.164;沙乌头猪的ROH平均长度为10.52 Mb,平均FROH为0.141。在高频ROH区域共鉴定到9个枫泾猪候选功能基因和8个梅山猪候选功能基因。研究结果为枫泾猪、梅山猪和沙乌头猪的资源保护和分子标记辅助选择奠定了基础。