CN110295236A - The SNP molecular genetic marker of pannage conversion ratio - Google Patents

The SNP molecular genetic marker of pannage conversion ratio Download PDF

Info

Publication number
CN110295236A
CN110295236A CN201910489840.1A CN201910489840A CN110295236A CN 110295236 A CN110295236 A CN 110295236A CN 201910489840 A CN201910489840 A CN 201910489840A CN 110295236 A CN110295236 A CN 110295236A
Authority
CN
China
Prior art keywords
genetic marker
snp
weight gain
feed
molecular genetic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910489840.1A
Other languages
Chinese (zh)
Other versions
CN110295236B (en
Inventor
赵云翔
邝伟键
李智丽
喻维维
朱晓萍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Foshan University
Original Assignee
Foshan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Foshan University filed Critical Foshan University
Priority to CN201910489840.1A priority Critical patent/CN110295236B/en
Publication of CN110295236A publication Critical patent/CN110295236A/en
Application granted granted Critical
Publication of CN110295236B publication Critical patent/CN110295236B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6888—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/124—Animal traits, i.e. production traits, including athletic performance or the like
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/156—Polymorphic or mutational markers
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/172—Haplotypes
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P60/00—Technologies relating to agriculture, livestock or agroalimentary industries
    • Y02P60/80—Food processing, e.g. use of renewable energies or variable speed drives in handling, conveying or stacking
    • Y02P60/87—Re-use of by-products of food processing for fodder production

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Organic Chemistry (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Microbiology (AREA)
  • Immunology (AREA)
  • Molecular Biology (AREA)
  • Biotechnology (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Genetics & Genomics (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

The SNP molecular genetic marker of disclosure offer pannage conversion ratio, the molecular labeling WU_10.2_3_116703757 of boar feed weight gain ratio is influenced by identification one, there were significant differences for the feed weight gain ratio of the label different genotype boar, and increase weight by boar feed than the association analysis with full-length genome molecular genetic marker, the big effector molecule genetic marker for influencing pannage weight gain ratio is screened in success, and there were significant differences for the feed weight gain ratio of the label different genotype boar；By detecting the molecular labeling, it can be applied to the breeding of boar, select and remain the homozygous pig of low feed weight gain ratio, it can effectively reduce feed consumption in production process, reduce pig production cost, Business Economic Benefit and competitiveness are improved, accelerates the breeding progress of high food utilization efficiency strain, effectively reduces feed consumption and aquaculture cost.

Description

The SNP molecular genetic marker of pannage conversion ratio
Technical field
This disclosure relates to pig gene technical field, in particular to the SNP molecular genetic marker of pannage conversion ratio.
Background technique
The feed efficiency economic characters important as one, are constantly subjected to pig raising enterprise both at home and abroad and boar improves company It pays close attention to.20th century mid-term, continue to use both at home and abroad feed weight gain than (Feed/Gain ratio, F/G) study feed efficiency, should Character is medium heritability quantitative character.(Mignon G S, rideau N, the Gabriel I, et such as Mignon in 2015 al.Detection of QTL controlling feed efficiency and excretion in chickens fed A wheat-based diet.Genetics Selection Evolution, 47 (1): 74 (2015)) it is detected in chicken 13 QTLs related with feed conversion rate.It is reported that (the Suppressor of of cell signalling inhibiting factor 2 of pig Cytokine Signalling 2, CRADD) gene, Melanocortin receptor 4 (Melanocortin 4 receptor, MC4R) Gene exist to the significant relevant single nucleotide polymorphism of feed conversion rate (single nucleotide polymorphism, SNP) site.(Ma Y, Yu C, Mohamed E M, the et al.A causal link from ALK to such as Ma hexokinase II overexpression and hyperactive glycolysis in EML4-ALK-positive Lung cancer.Oncogene, 35 (47), 6132-6142 (2016) .doi:10.1038/onc.2016.150) research table Bright, ALK gene (ALK Receptor Tyrosine Kinase) participates in reconciling body glucose metabolic process.Therefore, feed Medium heritability character of the conversion ratio as evaluation food utilization efficiency, has good selection to react.
Have using SNP marker assist-breeding food utilization efficiency correlated traits to pig production management and Business Economic Benefit It has a major impact.(Vigors S, Sweeney T, Oshea C J, et higher than the food utilization efficiency of pig 1. low feed increases weight al.Pigs that are divergent in feed efficiency,differ in intestinal enzyme and nutrient transporter gene expression,nutrient digestibility and microbial Activity.Animal, 10 (11): 1848-1855 (2016)), the feed usage amount and production cost in production can be reduced, into And feed resource has been saved, also the blowdown flow rate of pig can be reduced to a certain degree, to slow down pig and human competition grain resource and support The pressure of pig industry environmental issue.2. developing effective molecular labeling to work for the breeding of feed efficiency correlated traits, greatly shorten Cultivation period reduces and cultivates cost, improves seed selection accuracy, accelerates genetic progress, is avoided that and introduces a fine variety that-degeneration-introduces a fine variety again Phenomenon.
Therefore, it excavates and is had great significance using new increasing weight with feed than genetic breeding of the related gene for pig. The trait phenotypes record of high density SNP data and big group based on covering full-length genome, can pass through whole-genome association Technology (GWAS) (Hirschhorn, J.N.&Daly, M.J.Genome-wide association studies for Common diseases and complex traits.Nat.Rev.Genet.6,95-108 (2005)) control is accurately positioned The candidate gene of character.Although the technology still have some defects (De, R., Bush, W.S.&Moore, J.H.Bioinformatics challenges in genome-wide association studies(GWAS) .Methods Mol.Biol.1168,63-81 (2014)), it has been widely used in the excavation of mankind's complex disease candidate gene With the positioning of livestock and poultry important economical trait key gene.Classical GWAS is generally basede on Plink (Purcell, S.et al.PLINK:A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses.Am.J.Hum.Genet.813,559-575 (2007)) etc. softwares single label carried out to all labels one by one return Analysis sets a remarkable threshold then to screen significant site.Such methods, which often face, calculates that intensity is big, excessively high estimation is marked Remember that effect, conspicuousness threshold value set the problems such as unreasonable.In order to further increase the efficiency of GWAS, new method and software constantly quilt It proposes.Wherein, one-step method whole-genome association (wssGWAS, (WANG, H., MISZTAL, I., AGUILAR, I., LEGARRA,A.&MUIR,W.M.Genome-wide association mapping including phenotypes from relatives without genotypes.Genet Res 94,73–83(2012))；(Wang,H.et al.Genome- wide association mapping including phenotypes from relatives without genotypes in a single-step(ssGWAS)for 6-week body weight in broiler Chickens.Front.Genet.5,1-10 (2014)) simultaneously using pedigree, history individual phenotypic record and genotype data into Row association analysis, possess phenotypic record suitable for a large amount of individuals and the case where only a small amount of individual possesses genotype data, especially Whole-genome association suitable for livestock and poultry important economical trait.Based on GBLUPf90 software (Misztal, I.et al.BLUPF90 and related programs(BGF90).in Proc.7th World Congr.Genet.Appl.Livest.Prod.21-22 (2002) .doi:9782738010520), it can realize easily wssGWAS.It is filtered out to pannage weight gain using wssGWAS than relevant SNP marker, is domestic pig feed efficiency character Hereditary choosing amount provides feasible way, is of great significance to pig breeding industry.
Summary of the invention
In view of the above technical problems, the disclosure provides the SNP molecular genetic marker of pannage conversion ratio, passes through identification one The molecular labeling WU_10.2_3_116703757 of boar feed weight gain ratio is influenced, the feed of the label different genotype boar increases Again than there were significant differences, and increased weight by boar feed than the association analysis with full-length genome molecular genetic marker, success Screening influences the big effector molecule genetic marker of pannage weight gain ratio, the WU_10.2_3_116703757 heredity mark being directed to The mutational site that note i.e. No. SNP is WU_10.2_3_116703757, sees pig genome database (Sscrofa11.1) in NCBI.
The SNP molecular genetic marker (WU_10.2_3_116703757) of the disclosure, refering to Ensembl database (http://asia.ensembl.org/Sus_scrofa/Search/New? db=core), obtaining accession number is WU_10.2_ The genetic fragment (or No. RS is rs343880869) of 3_116703757, SNP molecular genetic marker is located at No. 3 chromosomes of pig The position 109882181bp, WU_10.2_3_116703757 belong to the intron sequences of ALK gene, which is C > T prominent Become (mutational site), C > T, that is, C be big frequency allele, T be small frequency allele, symbol > be gene frequency Size, the SNP molecular genetic marker are as follows in mutational site upstream and downstream 100bp sequence:
5’-ATACAATTCAACCCAAAACAGAACCTGAGGAAGAGCAATCTGTGCCCTATTTGTCCTGTAAGCGG TGACATTGTTTCAGTAAACAGCACTTCCATCACAAR(T/C)GAATTGTCGCTTAATTTGAATTGCATGGTTTTGAA CATTTGTATTTATATGCAAATAATTATTGAATTATTTTGTTTGGTTTTGATAGTCATGCAAGAGC-3'；R is mutation Site when the R at 101 nucleotide of above-mentioned sequence is C or T, i.e. when R (T/C), leads to above-mentioned sequence polymorphism；When above-mentioned core When 101st nucleotide of nucleotide sequence is T, pig has lower feed weight gain ratio, and 5 '-and -3 ' respectively indicate nucleotides sequence 5 ' the ends and 3 ' ends of column.
Feed between above-mentioned WU_10.2_3_116703757 marker genetype TT genotype individuals and CC genotype boar individual Weight gain reduces 8.41% than CC body feedstuff weight gain ratio of difference 0.18, TT genotype individuals ratio, so, T is to be conducive to significantly It is raw to can effectively reduce pig by selecting and remain the TT genotype homozygosis pig of low feed weight gain ratio for the allele for reducing feed weight gain ratio (since the DNA of pig is reverse acting spiral duplex structure, the mutational site nucleotide of two chains is T to feed consumption during production When be TT genotype homozygosis pig, wherein each chain has a nucleotide sequence, and T indicates that a mutational site is T, TT gene Type is that the mutational site of double-strand is all the homozygous pig of T, and similarly, CC genotype is that the mutational site of double-strand is all the homozygous pig of C；CT It is the pig of C that genotype, which is T another mutational site for the mutational site of a chain).
Screen pannage conversion ratio SNP molecular genetic marker method specifically includes the following steps:
1, the process step of the acquisition of molecular labeling
1.1, the ear tissue sample or blood sample of boar are acquired, extracts total DNA, and quality testing is carried out to DNA.Using GGP 50k SNP (GeneSeek, US) chip carries out Genotyping, obtains the SNP marker genotype of covering full-length genome.
1.2, genome (Sscrofa11.1) is referred to according to the pig of latest edition, using NCBI genome alignment program (https: //www.ncbi.nlm.nih.gov/) is updated the physical location of all SNP markers.Genomic locations are unknown SNP be not used in association analysis.
1.3, for the SNP marker on all autosomes, quality control, standard are carried out using Plink software are as follows: individual Recall rate >=90%；SNP recall rate >=90%；Small gene frequency >=0.01；Hardy-Weinberg equilibrium p value >=10-6.It is right In deletion Genotype, it is filled using Beagle software (version 4.1).
2, the process step of the verifying of molecular labeling
2.1, boar pedigree is arranged, mainly includes the information such as boar individual number, father, mother and nascent date.UsingThe growth data that formula records full-automatic boar performance test system (FIRE, the U.S.) carries out analysis acquisition Feed weight gain is used for phenotype-genotype association analysis than phenotypic data.Wherein, FCR is feed weight gain ratio；WaFor living body weight gain Amount；WfFor feed consumption.
2.2, statistical model, using one-step method whole-genome association method (the weighted single step of weighting Genome-wide association study, wssGWAS) carry out whole-genome association.This method is primarily based on mixing Model equation group estimates individual breeding value, and the equivalence relation then based on breeding value model and marker effect model is by breeding value Be converted to marker effect.The whole-genome association model that the present invention uses is as follows:
Y=Xb+Za+Wp+e,
Wherein, y is that observation vector is compared in feed weight gain；X, Z and W are design matrix；B is fixed effect vector (environment, day Age)；For breeding value vector；For the permanent environmental effect of individual；For residual error.H To integrate the affiliation matrix of pedigree and SNP marker simultaneously, inverse matrix calculation formula is as follows:
Wherein, A is the affiliation matrix based on pedigree；A22To there is the corresponding matrix in block form of genotype individuals in A；Gω= 0.9G+0.1A22,For the affiliation square based on full-length genome SNP marker, Z is small gene frequency Genotype matrix after (minor allele frequency, MAF) correction, wherein 0-2p, 1-2p and 2-2p respectively represent AA, Tri- kinds of genotype of Aa and aa, p are small gene frequency；D is diagonal matrix, indicates the weight of SNP；piIt is marked for i-th Small gene frequency；M is marker number.
For above-mentioned mixed model, using AI-REML (average information restricted maximum Likelihood) method estimate variance component, and breeding value is obtained by solving Mixed model mixed.It is obtained by way of iteration Weight must be marked, key step is as follows:
Step 1: initialization (t=1), D(t)=I, G(t)=λ ZD(t)Z ',
Step 2: individual breeding value is calculated by ssGBLUP；
Step 3: pass through formulaIndividual breeding value is converted into SNP effect, whereinTo there is gene The breeding value of type individual；
Step 4: formula is utilizedIt calculates SNP weight and is used for next round iteration；
Step 5: formula is utilizedSNP weight is standardized, to guarantee that variance is consistent；
Step 6: formula G is utilized(t+1)=λ ZD(t+1)Z ' calculating affiliation matrix is used for next round iteration；
Step 7: t=t+1, and the next round iteration since step 2 are enabled.
Above-mentioned steps iteration is three times, final to obtain SNP marker effect.The marker effect that third round iteration is exported is as most Whole result.Calculating process mainly calls BLUPF90 software to realize by statisticalling analyze platform programming in R, wherein AIREMLF90 program is used for variance component estimate, and BLUPF90 program is for calculating breeding value, and postGSf90 is for calculating label Effect.
2.3, label screening, effect value markd for institute, takes its absolute value to draw Manhattan figure, shows and screen big effect The SNP marker answered.And using variance analysis and Multiple range test (R statisticallys analyze platform), analysis WU_10.2_3_116703757 mark Difference condition is compared in note different genotype group boar feed weight gain.
The disclosure has the beneficial effect that the disclosure provides the SNP molecular genetic marker of pannage conversion ratio, and the label is different There were significant differences for the feed weight gain ratio of genotype boar；By detecting the molecular labeling, it can be applied to the breeding of boar, select and remain The homozygous pig of low feed weight gain ratio, can effectively reduce feed consumption in production process, reduces pig production cost, improve enterprise Economic benefit and competitiveness accelerate the breeding progress of high food utilization efficiency strain, effectively reduce feed consumption and cultivation at This.
Detailed description of the invention
By the way that the embodiment in conjunction with shown by attached drawing is described in detail, above-mentioned and other features of the disclosure will More obvious, identical reference label indicates the same or similar element in disclosure attached drawing, it should be apparent that, it is described below Attached drawing be only some embodiments of the present disclosure, for those of ordinary skill in the art, do not making the creative labor Under the premise of, it is also possible to obtain other drawings based on these drawings, in the accompanying drawings:
Fig. 1 show the method work flow diagram of the SNP molecular genetic marker of the screening pannage conversion ratio of the disclosure；
Fig. 2 show the marker gene group position WU_10.2_3_116703757 of the disclosure and full genome is compared in feed weight gain Group SNP effect distribution.
Specific embodiment
It is carried out below with reference to technical effect of the embodiment and attached drawing to the design of the disclosure, specific structure and generation clear Chu, complete description, to be completely understood by the purpose, scheme and effect of the disclosure.It should be noted that the case where not conflicting Under, the features in the embodiments and the embodiments of the present application can be combined with each other.
As shown in Figure 1 for according to the method workflow of the SNP molecular genetic marker of the screening pannage conversion ratio of the disclosure Cheng Tu illustrates the method for the SNP molecular genetic marker of the screening pannage conversion ratio according to the disclosure below with reference to Fig. 1.
The method that the disclosure screens the SNP molecular genetic marker of pannage conversion ratio, specifically includes the following steps:
(1) phenotype-pedigree data acquisition
The basic research group of the disclosure is Duroc boars, all is from Guangxi Yang Xiang Farming Ltd. kind pig farm. It include 4 generations, 735 boars in complete pedigree, the feed that 370 Duroc boars are wherein had recorded between 2015-2018 increases Again than trait phenotypes data.Feed weight gain is than usingFormula is to full-automatic boar performance test system The growth data of (FIRE, the U.S.) record carries out analysis acquisition, is used for phenotype-genotype association analysis.Wherein, FCR is feed Weight gain ratio；WaFor living body gain in weight；WfFor feed consumption.
(2) Genotyping and quality control
The ear tissue sample or blood sample of 1733 boars are acquired, extracts total DNA, and use GGP 50k SNP (GeneSeek, US) chip carries out Genotyping, obtains 50705 SNP markers of covering full-length genome.According to the pig of latest edition With reference to genome (Sscrofa11.1), using NCBI genome alignment program (https: //www.ncbi.nlm.nih.gov/) The physical location of all SNP markers is updated.The unknown SNP of genomic locations is not used in association analysis.For all normal SNP marker on chromosome carries out quality control, standard are as follows: individual recall rate >=90% using Plink software；SNP recall rate >=90%；Small gene frequency >=0.01；Hardy-Weinberg equilibrium p value >=10-6.For deletion Genotype, using Beagle Software (version 4.1) is filled.Based on the above quality control standard, remaining 1623 boars and 28289 SNP markers For association analysis.
(3) statistical model
In order to make full use of all phenotypic datas and genotype data, the present invention discloses the one-step method full genome using weighting Group correlation fractal dimension (weighted single step genome-wide association study, wssGWAS) carries out Whole-genome association.This method is primarily based on Mixed model mixed to estimate individual breeding value, is then based on breeding value Breeding value is converted to marker effect by the equivalence relation of model and marker effect model.The full-length genome association point that the present invention uses It is as follows to analyse model:
Y=Xb+Za+Wp+e,
Wherein, y is that observation vector is compared in feed weight gain；X, Z and W are design matrix；B is fixed effect vector (environment, day Age)；For breeding value vector；For the permanent environmental effect of individual；For residual error.H To integrate the affiliation matrix of pedigree and SNP marker simultaneously, inverse matrix calculation formula is as follows:
Wherein, A is the affiliation matrix based on pedigree；A22To there is the corresponding matrix in block form of genotype individuals in A；Gω= 0.9G+0.1A22,For the affiliation square based on full-length genome SNP marker, Z is small gene frequency Genotype matrix after (minor allele frequency, MAF) correction, wherein 0-2p, 1-2p and 2-2p respectively represent AA, Tri- kinds of genotype of Aa and aa, p are small gene frequency；D is diagonal matrix, indicates the weight of SNP；piIt is marked for i-th Small gene frequency；M is marker number.
Corresponding above-mentioned mixed model, using AI-REML (average information restricted maximum Likelihood) method estimate variance component, and breeding value is obtained by solving Mixed model mixed.It is obtained by way of iteration Weight must be marked, key step is as follows:
Step 1: initialization (t=1), D(t)=I, G(t)=λ ZD(t)Z ',
Step 2: individual breeding value is calculated by ssGBLUP；
Step 3: pass through formulaIndividual breeding value is converted into SNP effect, whereinTo there is gene The breeding value of type individual；
Step 4: formula is utilizedIt calculates SNP weight and is used for next round iteration；
Step 5: formula is utilizedSNP weight is standardized, to guarantee that variance is consistent；
Step 6: formula G is utilized(t+1)=λ ZD(t+1)Z ' calculating affiliation matrix is used for next round iteration；
Step 7: t=t+1, and the next round iteration since step 2 are enabled.
Above-mentioned steps iteration is three times, final to obtain SNP marker effect to get SNP marker effect is arrived.Third round iteration is defeated Marker effect out is as final result.Calculating process mainly calls BLUPF90 software by statisticalling analyze platform programming in R It realizes, wherein AIREMLF90 program is used for variance component estimate, and BLUPF90 program is for calculating breeding value, postGSf90 For calculating marker effect.
(4) label screening
Effect value markd for institute, takes its absolute value to draw Manhattan figure, shows and screen the SNP marker of big effect. And using variance analysis and Multiple range test (R statisticallys analyze platform), analysis WU_10.2_3_116703757 marks different genotype Difference condition is compared in the weight gain of group's boar feed.
Analyze different genotype boar feed weight gain ratio
Effect value markd for institute, takes its absolute value to draw Manhattan figure, shows and screen the SNP marker of big effect (as shown in Fig. 2, Fig. 2 show the marker gene group position WU_10.2_3_116703757 of the disclosure and feed increases weight than full base Because of a group SNP effect distribution).And using variance analysis and Multiple range test (R statisticallys analyze platform), analysis different genotype group is public Pannage increases weight than difference condition (table 1).
The SNP molecular genetic marker (WU_10.2_3_116703757 marker gene sequence) of pannage conversion ratio is raised in pig Expect the application in conversion ratio assisted Selection:
The disclosure identifies the molecular labeling WU_10.2_3_116703757 for influencing boar feed weight gain ratio, by table 1 It is found that the feed weight gain of the label different genotype boar is than there were significant differences, (table 1 is WU_10.2_3_116703757 label Different genotype boar feed weight gain ratio)；
1 WU_10.2_3_116703757 of table marks different genotype boar feed weight gain ratio
Boar breeding is assisted by detection WU_10.2_3_116703757 marker genetype, it can be by selecting and remain homozygous kind of TT Pig enters core group, reduces feed weight gain ratio, effectively reduces feed consumption and aquaculture cost；
The SNP molecular genetic marker (WU_10.2_3_116703757) of the disclosure is located at No. 3 chromosomes of pig The position 109882181bp belongs to the intron sequences of ALK gene, which is that a C > T is mutated (Sscrofa11.1), should SNP molecular genetic marker is as follows in mutational site upstream and downstream 100bp sequence:
5’-ATACAATTCAACCCAAAACAGAACCTGAGGAAGAGCAATCTGTGCCCTATTTGTCCTGTAAGCGG TGACATTGTTTCAGTAAACAGCACTTCCATCACAAR(T/C)GAATTGTCGCTTAATTTGAATTGCATGGTTTTGAA CATTTGTATTTATATGCAAATAATTATTGAATTATTTTGTTTGGTTTTGATAGTCATGCAAGAGC-3'；R is mutation Site when the R at 101 nucleotide of above-mentioned sequence is C or T, i.e. when R (T/C), leads to above-mentioned sequence polymorphism；When above-mentioned core When 101st nucleotide of nucleotide sequence is T, pig has lower feed weight gain ratio, and 5 '-and -3 ' respectively indicate nucleotides sequence 5 ' the ends and 3 ' ends of column.
(sequence such as sequence table SEQ IDNo.1 shown in nucleotide sequence of the above-mentioned sequence when being mutated point and being T), sequence List SEQ IDNo.1 be the present invention screen genetic marker (i.e. No. SNP be WU_10.2_3_116703757, No. RS is Rs343880869 the nucleotide sequence of mutational site upstream and downstream 100bp).
Feed weight gain ratio differs 0.18 between above-mentioned WU_10.2_3_116703757 marker genetype TT and CC boar individual, CC body feedstuff weight gain ratio of TT individual ratio reduces 8.41%, so, T is the equipotential base for being conducive to significantly reduce feed weight gain ratio Cause.
Leading reference:
1.Mignon G S,rideau N,Gabriel I,et al.Detection of QTL controlling feed efficiency and excretion in chickens fed a wheat-based diet.Genetics Selection Evolution,47(1):74(2015)。
2.Ma Y,Yu C,Mohamed E M,et al.A causal link from ALK to hexokinase II overexpression and hyperactive glycolysis in EML4-ALK-positive lung cancer.Oncogene,35(47),6132–6142(2016).doi:10.1038/onc.2016.150。
3.Vigors S,Sweeney T,Oshea C J,et al.Pigs that are divergent in feed efficiency,differ in intestinal enzyme and nutrient transporter gene expression,nutrientdigestibility and microbial activity.Animal,10(11):1848- 1855(2016)。
4.Hirschhorn,J.N.&Daly,M.J.Genome-wide association studies for common diseases and complex traits.Nat.Rev.Genet.6,95–108(2005)。
5.De,R.,Bush,W.S.&Moore,J.H.Bioinformatics challenges in genome-wide association studies(GWAS).Methods Mol.Biol.1168,63–81(2014)。
6.Purcell,S.et al.PLINK:A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses.Am.J.Hum.Genet.813,559–575(2007)。
7.WANG,H.,MISZTAL,I.,AGUILAR,I.,LEGARRA,A.&MUIR,W.M.Genome-wide association mapping including phenotypes from relatives without genotypes.Genet Res 94,73–83(2012)。
8.Wang,H.et al.Genome-wide association mapping including phenotypes from relatives without genotypes in a single-step(ssGWAS)for 6-week body weight in broiler chickens.Front.Genet.5,1–10(2014)。
9.Misztal,I.et al.BLUPF90 and related programs(BGF90).in Proc.7th World Congr.Genet.Appl.Livest.Prod.21–22(2002).doi:9782738010520。
Sequence table
<110>Foshan Science &. Technology College
<120>the SNP molecular genetic marker of pannage conversion ratio
<141> 2019-06-05
<160> 1
<170> SIPOSequenceListing 1.0
<210> 1
<211> 201
<212> DNA
<213> Sscrofa11.1
<220>
<221> gene
<222> (1)..(201)
<220>
<221> mutation
<222> (100)..(101)
<400> 1
atacaattca acccaaaaca gaacctgagg aagagcaatc tgtgccctat ttgtcctgta 60
agcggtgaca ttgtttcagt aaacagcact tccatcacaa tgaattgtcg cttaatttga 120
attgcatggt tttgaacatt tgtatttata tgcaaataat tattgaatta ttttgtttgg 180
ttttgatagt catgcaagag c 201

Claims (6)

1. the SNP molecular genetic marker of pannage conversion ratio, which is characterized in that the molecular genetic marker is SNP molecular genetic Label, the SNP molecular genetic marker are located at No. 3 positions chromosome 109882181bp of pig, belong to the introne sequence of ALK gene Column, the position are C > T mutation, and pig is Sscrofa11.1 with reference to genome.
2. the SNP molecular genetic marker of pannage conversion ratio according to claim 1, which is characterized in that the SNP molecule The sequence of genetic marker is the upstream and downstream 100bp sequence in mutational site.
3. the SNP molecular genetic marker of pannage conversion ratio according to claim 2, which is characterized in that the SNP molecule The sequence of genetic marker is as shown below:
5’-ATACAATTCAACCCAAAACAGAACCTGAGGAAGAGCAATCTGTGCCCTATTTGTCCTGTAAGCGGTGAC ATTGTTTCAGTAAACAGCACTTCCATCACAARGAATTGTCGCTTAATTTGAATTGCATGGTTTTGAACATTTGTAT TTATATGCAAATAATTATTGAATTATTTTGTTTGGTTTTGATAGTCATGCAAGAGC-3'；R is mutational site, when R is When T, pig significantly reduces feed weight gain ratio.
4. the SNP molecular genetic marker of pannage conversion ratio according to claim 1, which is characterized in that the SNP molecule No. SNP of genetic marker is WU_10.2_3_116703757.
5. the SNP molecular genetic marker of pannage conversion ratio as claimed in claim 3 is in pannage conversion ratio assisted Selection Using.
6. the SNP molecular genetic marker of pannage conversion ratio according to claim 5 is in pannage conversion ratio assisted Selection In application, which is characterized in that raised between the TT genotype individuals of SNP molecular genetic marker genotype and CC genotype boar individual Expect that weight gain reduces 8.41% than CC body feedstuff weight gain ratio of difference 0.18, TT individual ratio, by selecting and remain low feed weight gain ratio TT genotype homozygosis pig, can effectively reduce feed consumption in pig production process.
CN201910489840.1A 2019-06-06 2019-06-06 SNP molecular genetic marker for pig feed conversion rate Active CN110295236B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910489840.1A CN110295236B (en) 2019-06-06 2019-06-06 SNP molecular genetic marker for pig feed conversion rate

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910489840.1A CN110295236B (en) 2019-06-06 2019-06-06 SNP molecular genetic marker for pig feed conversion rate

Publications (2)

Publication Number Publication Date
CN110295236A true CN110295236A (en) 2019-10-01
CN110295236B CN110295236B (en) 2023-05-30

Family

ID=68027594

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910489840.1A Active CN110295236B (en) 2019-06-06 2019-06-06 SNP molecular genetic marker for pig feed conversion rate

Country Status (1)

Country Link
CN (1) CN110295236B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111500746A (en) * 2020-05-22 2020-08-07 华中农业大学 SNP molecular marker related to feed conversion efficiency of pigs
CN113699246A (en) * 2021-07-26 2021-11-26 华南农业大学 SNP molecular marker influencing pig feed conversion efficiency traits and application thereof
CN115341045A (en) * 2022-10-19 2022-11-15 佛山科学技术学院 Method for predicting pig feed conversion rate by using microorganisms and related SNP sites thereof

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005112544A2 (en) * 2004-02-19 2005-12-01 The Governors Of The University Of Alberta Leptin promoter polymorphisms and uses thereof
WO2007073160A2 (en) * 2005-12-19 2007-06-28 Nutreco Nederland B.V. Methods for improving turkey meat production
CN104250646A (en) * 2013-06-27 2014-12-31 华中农业大学 Molecular marker correlated with pig feed conversion efficiency characters and detection method and application
CN104774836A (en) * 2015-04-15 2015-07-15 兰州大学 Polygene pyramiding early-breeding method for raising litter size of sheep
CN105316412A (en) * 2015-11-16 2016-02-10 中国农业科学院北京畜牧兽医研究所 Method for identification or auxiliary identification of day age of pig at 100kg weight and special kit of method
CN105624155A (en) * 2016-02-29 2016-06-01 华南农业大学 Molecular marker influencing feed conversion ratio character of pig and application
CN107937556A (en) * 2017-11-14 2018-04-20 中国农业大学 One and the relevant SNP site of pannage conversion ratio and its application
CN108559781A (en) * 2018-03-28 2018-09-21 中国农业科学院北京畜牧兽医研究所 A method of cultivating high food utilization efficiency pig
CN109402270A (en) * 2018-12-07 2019-03-01 佛山科学技术学院 One kind SNP marker relevant to Large White growth traits and its application

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005112544A2 (en) * 2004-02-19 2005-12-01 The Governors Of The University Of Alberta Leptin promoter polymorphisms and uses thereof
WO2007073160A2 (en) * 2005-12-19 2007-06-28 Nutreco Nederland B.V. Methods for improving turkey meat production
CN104250646A (en) * 2013-06-27 2014-12-31 华中农业大学 Molecular marker correlated with pig feed conversion efficiency characters and detection method and application
CN104774836A (en) * 2015-04-15 2015-07-15 兰州大学 Polygene pyramiding early-breeding method for raising litter size of sheep
CN105316412A (en) * 2015-11-16 2016-02-10 中国农业科学院北京畜牧兽医研究所 Method for identification or auxiliary identification of day age of pig at 100kg weight and special kit of method
CN105624155A (en) * 2016-02-29 2016-06-01 华南农业大学 Molecular marker influencing feed conversion ratio character of pig and application
CN107937556A (en) * 2017-11-14 2018-04-20 中国农业大学 One and the relevant SNP site of pannage conversion ratio and its application
CN108559781A (en) * 2018-03-28 2018-09-21 中国农业科学院北京畜牧兽医研究所 A method of cultivating high food utilization efficiency pig
CN109402270A (en) * 2018-12-07 2019-03-01 佛山科学技术学院 One kind SNP marker relevant to Large White growth traits and its application

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
BARBARA U METZLER-ZEBELI等: "Finishing pigs that are divergent in feed efficiency show small differences in intestinal functionality and structure", 《PLOS ONE》 *
CHUNYAN等: "Genomic evaluation of feed efficiency component traits in Duroc pigs using 80K, 650K and whole-genome sequence variants", 《GENET SEL EVOL》 *
ENSEMBL RELEASE 96: ""rs343880869"", 《EMBL-EBI》 *
SAHANA等: "A genome-wide association scan in pig identifies novel regions associated with feed efficiency trait", 《JOURNAL OF ANIMAL SCIENCE》 *
蒲蕾等: "杜洛克猪HMGA1基因多态位点与生长、饲料利用性状的关联分析", 《中国畜牧兽医》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111500746A (en) * 2020-05-22 2020-08-07 华中农业大学 SNP molecular marker related to feed conversion efficiency of pigs
CN113699246A (en) * 2021-07-26 2021-11-26 华南农业大学 SNP molecular marker influencing pig feed conversion efficiency traits and application thereof
CN113699246B (en) * 2021-07-26 2023-07-11 华南农业大学 SNP molecular marker affecting pig feed conversion efficiency character and application thereof
CN115341045A (en) * 2022-10-19 2022-11-15 佛山科学技术学院 Method for predicting pig feed conversion rate by using microorganisms and related SNP sites thereof

Also Published As

Publication number Publication date
CN110295236B (en) 2023-05-30

Similar Documents

Publication Publication Date Title
CN110191965B (en) Pig whole genome 50K SNP chip and application
CN110218799A (en) The molecular genetic marker of pig residue feed intake character and application
WO2023165128A1 (en) Chicken low-density snp liquid-phase chip based on target capture sequencing and use thereof
CN112002371B (en) Genome selection method for residual feed intake of white-feather broilers
CN110358840A (en) The SNP molecular genetic marker of TPP2 gene relevant to remaining feed intake
CN110295236A (en) The SNP molecular genetic marker of pannage conversion ratio
CN110358839A (en) The SNP molecular genetic marker of GCKR gene relevant to pannage conversion ratio
CN110358838A (en) SNP genetic marker relevant to pannage conversion in FA2H genetic fragment
CN111926085B (en) Molecular marker influencing chicken muscle brightness and application thereof
CN111926086B (en) Molecular marker influencing oblique growth of chicken body and application thereof
CN114941033A (en) Method for breeding local high-quality white-feather chicken high-egg-yield strain based on SNP locus assistance
CN114921561B (en) Duroc whole genome low-density SNP chip and preparation method and application thereof
CN117500943A (en) SNP molecular marker combination and chip for Beijing black pig genotyping and preparation method and application thereof
CN110273006A (en) The relevant molecular genetic marker of the effective sperm count of one herd boar
CN109735633A (en) The detection method and its application of FSHR gene specific SNP marker, the black sheep litter size character in Turfan
Lv et al. Genomic differentiation and selection signatures of two elite varieties of Yesso scallop Mizuhopecten yessoensis
US20240043912A1 (en) Genomic selection (gs) breeding chip of huaxi cattle and use thereof
CN111199773B (en) Evaluation method for fine positioning character associated genome homozygous fragments
CN111370058B (en) Method for tracing buffalo blood line source and carrying out genome matching based on whole genome SNP information
CN114752678B (en) SNP molecular marker related to backfat thickness of pig reaching 115kg body weight and application thereof
CN111910009B (en) Molecular marker influencing chicken bursal disease index and application thereof
CN112176073B (en) PROS1 gene molecular marker related to chicken carcass traits and application
CN113549699A (en) Genome selection method for egg number of white feather broilers
Kang et al. Genome-wide association study for sow lifetime productivity related traits in a Landrace purebred population
Hang et al. Genetic Diversity and Breeding Signatures for Regional Indica Rice Improvement in Guangdong of Southern China

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant