CN118109605A - SNP molecular marker combination related to growth traits of Nile-Lafei buffalo and application - Google Patents
SNP molecular marker combination related to growth traits of Nile-Lafei buffalo and application Download PDFInfo
- Publication number
- CN118109605A CN118109605A CN202410150893.1A CN202410150893A CN118109605A CN 118109605 A CN118109605 A CN 118109605A CN 202410150893 A CN202410150893 A CN 202410150893A CN 118109605 A CN118109605 A CN 118109605A
- Authority
- CN
- China
- Prior art keywords
- buffalo
- snp
- chr2
- raffia
- marker
- 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
Links
- 239000003147 molecular marker Substances 0.000 title claims abstract description 36
- 238000012098 association analyses Methods 0.000 claims abstract description 23
- 238000009395 breeding Methods 0.000 claims description 36
- 230000001488 breeding effect Effects 0.000 claims description 36
- 239000002773 nucleotide Substances 0.000 claims description 25
- 125000003729 nucleotide group Chemical group 0.000 claims description 25
- 239000003550 marker Substances 0.000 claims description 24
- 238000000034 method Methods 0.000 claims description 17
- 238000012163 sequencing technique Methods 0.000 claims description 15
- 108020004414 DNA Proteins 0.000 claims description 14
- 230000003321 amplification Effects 0.000 claims description 8
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 8
- 238000003205 genotyping method Methods 0.000 claims description 6
- 108091028043 Nucleic acid sequence Proteins 0.000 claims description 4
- 230000035772 mutation Effects 0.000 claims description 4
- 238000012216 screening Methods 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 13
- 238000005516 engineering process Methods 0.000 description 12
- 230000002068 genetic effect Effects 0.000 description 12
- 230000000694 effects Effects 0.000 description 11
- 108090000623 proteins and genes Proteins 0.000 description 11
- 238000011161 development Methods 0.000 description 9
- 238000011160 research Methods 0.000 description 9
- 239000008280 blood Substances 0.000 description 8
- 210000004369 blood Anatomy 0.000 description 8
- 241001465754 Metazoa Species 0.000 description 6
- 238000001514 detection method Methods 0.000 description 6
- 230000006872 improvement Effects 0.000 description 6
- 238000003908 quality control method Methods 0.000 description 6
- 239000000047 product Substances 0.000 description 5
- 238000007400 DNA extraction Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000012795 verification Methods 0.000 description 4
- 108700028369 Alleles Proteins 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 244000144972 livestock Species 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 241000283690 Bos taurus Species 0.000 description 2
- 238000012408 PCR amplification Methods 0.000 description 2
- 238000000246 agarose gel electrophoresis Methods 0.000 description 2
- 210000000349 chromosome Anatomy 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000008303 genetic mechanism Effects 0.000 description 2
- 210000004731 jugular vein Anatomy 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 235000013372 meat Nutrition 0.000 description 2
- 235000013336 milk Nutrition 0.000 description 2
- 239000008267 milk Substances 0.000 description 2
- 210000004080 milk Anatomy 0.000 description 2
- 108020004707 nucleic acids Proteins 0.000 description 2
- 150000007523 nucleic acids Chemical class 0.000 description 2
- 102000039446 nucleic acids Human genes 0.000 description 2
- 244000144977 poultry Species 0.000 description 2
- 238000013517 stratification Methods 0.000 description 2
- 238000001712 DNA sequencing Methods 0.000 description 1
- 208000035240 Disease Resistance Diseases 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000000137 annealing Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000013365 dairy product Nutrition 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004925 denaturation Methods 0.000 description 1
- 230000036425 denaturation Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000000499 gel Substances 0.000 description 1
- 238000012252 genetic analysis Methods 0.000 description 1
- 230000007614 genetic variation Effects 0.000 description 1
- 229920000669 heparin Polymers 0.000 description 1
- ZFGMDIBRIDKWMY-PASTXAENSA-N heparin Chemical compound CC(O)=N[C@@H]1[C@@H](O)[C@H](O)[C@@H](COS(O)(=O)=O)O[C@@H]1O[C@@H]1[C@@H](C(O)=O)O[C@@H](O[C@H]2[C@@H]([C@@H](OS(O)(=O)=O)[C@@H](O[C@@H]3[C@@H](OC(O)[C@H](OS(O)(=O)=O)[C@H]3O)C(O)=O)O[C@@H]2O)CS(O)(=O)=O)[C@H](O)[C@H]1O ZFGMDIBRIDKWMY-PASTXAENSA-N 0.000 description 1
- 229960001008 heparin sodium Drugs 0.000 description 1
- 238000012165 high-throughput sequencing Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 235000013622 meat product Nutrition 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000012257 pre-denaturation Methods 0.000 description 1
- 239000012264 purified product Substances 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000009394 selective breeding Methods 0.000 description 1
- 238000011451 sequencing strategy Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
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
- 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/6844—Nucleic acid amplification reactions
- C12Q1/6858—Allele-specific amplification
-
- 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
-
- 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)
- Organic Chemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biophysics (AREA)
- Immunology (AREA)
- Microbiology (AREA)
- Molecular Biology (AREA)
- Biotechnology (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Genetics & Genomics (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
The invention belongs to the technical field of molecular biology, and discloses SNP molecular marker combinations related to growth traits of Nile-Laffy buffalo and application thereof based on whole genome association analysis screening, wherein the SNP molecular marker combinations comprise at least one or more of 8 SNP markers which are ：chr1_54231：A/G、chr9_306973：G/A、chr9_999337：G/A、chr9_1963317：T/C、chr2_173473：C/T、chr2_1811972：C/T、chr2_789901：C/A、chr2_280007：C/T. respectively.
Description
Technical Field
The invention belongs to the technical field of molecular biology, and particularly relates to SNP molecular marker combination related to growth traits of Nile-Lafei buffalo based on whole genome association analysis screening and application thereof.
Background
The Nile-Lafei buffalo has the characteristics of heat resistance, coarse feeding resistance, strong disease resistance, normal growth, strong adaptability and the like in the subtropical climate environment in the south of China, and has excellent milk and meat dual-purpose performance compared with the local swamp buffalo. The hybrid improvement is carried out by using the Nile-Lafei buffalo and the local buffalo, and the meat and milk production performance of the hybrid offspring can be greatly improved.
The growth characteristics of buffalo are generally measured by the height, the body oblique length, the chest circumference and the like, and are one of the most important economic characteristics of buffalo, and the growth and development speed of buffalo directly influences the development and the growth of buffalo industry. Since the growth traits belong to quantitative traits and are controlled by a plurality of genes, the final phenotype data can be collected only in the adult period of buffalo by conventional phenotype-based breeding, and early evaluation and breeding of the growth traits are difficult to realize. With the development of high-throughput sequencing technology, the whole genome resequencing technology is increasingly widely applied to genetic analysis and genetic improvement of important economic traits of livestock and poultry. Through molecular marker assisted selective breeding, the breeding period of buffalo can be greatly shortened, and the genetic progress speed of buffalo is greatly improved. Because the molecular marker assisted selection is not easily affected by the environment and has no sex and age limitation, the molecular marker assisted selection is used for early seed selection, the generation interval can be shortened, the selection intensity is improved, and the seed selection efficiency and accuracy are improved. Therefore, screening molecular markers related to development of buffalo growth traits and developing a molecular marker assisted breeding technology become effective means for improving the production performance of buffalo in China, and have important significance for accelerating breeding, seed selection and seed production of improved breed buffalo meat and dairy products.
In early studies on complex trait phenotypes of livestock and poultry, QTL localization (QTL MAPPING) was mainly used to find the association between phenotypic and genomic variations. Although the method is a classical gene positioning method, the method is more suitable for genetic research of single-gene diseases or single-gene control traits, has limited detection effect on complex traits and traits with low genetic power, obtains a large QTL confidence interval, possibly contains hundreds of genes, is not beneficial to accurate positioning of subsequent functional genes, and has poor repeatability in different groups.
Compared with the QTL linkage analysis method, the whole Genome linkage analysis (Genome-wide association study, GWAS) is a method for carrying out overall linkage analysis on common genetic variation (single nucleotide polymorphism and copy number) in the whole Genome range, the method takes natural population as a research object, combines the diversity of a target character phenotype and the polymorphism of genes (or marker loci) based on linkage disequilibrium (linkage disequilibrium, LD) among genes (loci) remained after long-term recombination, and can directly identify the gene loci or marker loci which are closely related to the phenotype variation and have specific functions. The GWAS technology is adopted to conduct research in the whole genome range, multiple characters can be located at one time, the method is suitable for research in the aspects of locating character association intervals, functional gene research, development character breeding, functional marking and the like, and the obtained results are more reliable. The GWAS technology is taken as a new method and can be widely applied to the field of livestock breeding.
Through the above analysis, the problems and defects existing in the prior art are as follows: by adopting the technology, a genetic mechanism affecting the growth traits of the Nile-Lafei buffalo is revealed, and the application of a molecular marker with a forward effect in genome selection of the Nile-Lafei buffalo is a problem to be solved in the current urgent need for establishing an efficient and accurate gene breeding technology.
Disclosure of Invention
Aiming at the problems existing in the prior art, the invention provides SNP molecular marker combination related to the growth traits of Nile-Lafei buffalo and application thereof based on whole genome association analysis screening; can simply and quickly identify the Nili-Lafei buffalo individuals with good growth characteristics.
The invention is realized in such a way that SNP molecular marker combinations related to the growth traits of the Nile-Lafei buffalo screened based on whole genome association analysis comprise at least one or a plurality of combinations of the following 8 SNP markers, and the positions and variation information of the SNP molecular marker loci adopt chromosome_physical positions: the information of the reference genotype/variant genotype, which is ：chr1_54231：A/G、chr9_306973：G/A、chr9_999337：G/A、chr9_1963317：T/C、chr2_173473：C/T、chr2_1811972：C/T、chr2_789901：C/A、chr2_280007：C/T,8 SNP molecular markers respectively, is shown in Table 4, the base sequence of 101bp before and after each SNP molecular marker is provided in Table 4, and the base variation exists at 51 bp.
Further, the SNP marker is numbered chr1_54231: the nucleotide sequence of A/G is shown as SEQ ID NO. 1;
the SNP marker is numbered chr9_306973: the nucleotide sequence of the G/A is shown as SEQ ID NO. 2;
The SNP marker is numbered chr9_999337: the nucleotide sequence of the G/A is shown as SEQ ID NO. 3;
the SNP marker is numbered chr9_1963317: the nucleotide sequence of the T/C is shown as SEQ ID NO. 4;
SNP markers numbered chr2_173473: the nucleotide sequence of the C/T is shown as SEQ ID NO. 5;
the SNP marker is numbered chr2_1811972: the nucleotide sequence of the C/T is shown as SEQ ID NO. 6;
The SNP marker is numbered chr2_789901: the nucleotide sequence of the C/A is shown as SEQ ID NO. 7;
The SNP marker is numbered chr2_280007: the nucleotide sequence of the C/T is shown as SEQ ID NO. 8.
The invention also aims to provide the application of the SNP molecular marker combination related to the growth traits of the Nile-Laffy buffalo based on whole genome association analysis screening in the auxiliary breeding of the growth traits of the Nile-Laffy buffalo.
Another object of the present invention is to provide a method for detecting the genotype of niri-raffia buffalo using molecular biological technology based on a combination of SNP molecular markers related to growth traits of niri-raffia buffalo screened by whole genome association analysis, comprising the steps of: and respectively designing amplification primers according to nucleotide sequences of two flanks of the 8 SNP molecular marker loci, amplifying by taking the DNA of the Nile-Lafei buffalo individuals as a template, sequencing the first generation of amplification products, and genotyping the Nile-Lafei buffalo individuals according to a sequencing result, thereby identifying the genotypes of the Nile-Lafei buffalo individuals to be detected.
The invention also aims to provide application of the method for detecting the genotype of the Nile-Lafei buffalo by utilizing a molecular biological technology in the molecular marker assisted breeding of the growth trait of the Nile-Lafei buffalo.
In combination with the technical scheme and the technical problems to be solved, the technical scheme to be protected has the following advantages and positive effects:
first, the invention provides a SNP molecular marker combination related to the growth trait phenotype of Nile-Laffy buffalo, and uses whole genome association analysis to screen SNP molecular markers obviously related to the growth trait of Nile-Laffy buffalo, which can be used for auxiliary breeding of Nile-Laffy buffalo molecular markers and quicken breeding of improved variety of Nile-Laffy buffalo with good growth trait. The SNP molecular marker is used for identifying the growth traits of the Nile-Lafei buffalo, so that the breeding selection efficiency can be remarkably improved, the breeding period can be shortened, the breeding cost can be reduced, the breeding efficiency can be improved, and the breeding process can be accelerated.
Secondly, the invention based on SNP molecular marker combination related to the growth trait of Nile-Lafei buffalo screened by whole genome association analysis brings the following remarkable technical effects:
1. Improving accuracy of growth trait prediction
By precisely screening 8 SNP molecular markers remarkably related to the growth traits of the Nile-Lafei buffalo, the method can more accurately predict the growth traits of buffalo. The high-precision prediction provides scientific basis for breeding selection, and is helpful for improving breeding efficiency and breeding precision.
2. Accelerating the progress of genetic improvement
The SNP markers are used as tools for molecular auxiliary selection, so that the genetic improvement process of the Nile-Lafei buffalo can be accelerated. By screening and applying the key SNP loci, individuals with excellent growth traits can be rapidly identified, so that aggregation of excellent genes in a population is accelerated.
3. Promote the development of accurate animal husbandry
The application of the invention promotes the development of accurate animal husbandry, and realizes the accurate management and utilization of animal genetic resources through accurate identification on molecular level. The technology content of the animal husbandry is improved, and a powerful technical support is provided for sustainable development of the animal husbandry.
4. Optimizing breeding strategies and management
Through intensive research on SNP markers related to growth traits, the invention is helpful for optimizing breeding strategies and management. Breeders can select a breeding scheme more pertinently according to the information of the SNP markers, and more efficient genetic management and breeding decision can be implemented.
5. Providing new study objects for bioinformatics and genomics study
The SNP molecular marker identified in the invention not only has application value in practical breeding work, but also provides new research objects and data resources for the research in the fields of bioinformatics, genomics and the like. These data are of great significance for a thorough understanding of the genetic background and genetic mechanisms of growth traits of nii-raffia buffalo.
In conclusion, the invention provides an effective tool for buffalo breeding by identifying SNP molecular marker combinations related to the growth traits of Nile-Lafei buffalo, and has the remarkable technical effects of improving prediction accuracy, accelerating genetic improvement, developing accurate animal husbandry, optimizing breeding strategies, providing new resources for scientific research and the like.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are needed in the embodiments of the present invention will be briefly described below, and it is obvious that the drawings described below are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a high Manhattan diagram provided by an embodiment of the present invention;
FIG. 2 is a diagram of a body height Quantile-quantile provided by an embodiment of the present invention;
FIG. 3 is a Manhattan diagram of body diagonal length provided by an embodiment of the present invention;
FIG. 4 is a diagram of a body diagonal Quantile-quantile provided by an embodiment of the present invention;
FIG. 5 is a Manhattan diagram of a chest circumference provided by an embodiment of the present invention;
Fig. 6 is a diagram Quantile-quantile of a chest circumference provided by an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the following examples in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
The SNP molecular marker combination related to the growth trait of Nile-Laffy buffalo, which is screened based on whole genome association analysis, comprises at least one or a combination of the following 8 SNP markers, wherein the position and variation information of the SNP molecular marker loci adopt chromosome_physical positions: the information of the reference genotype/variant genotype, which is ：chr1_54231：A/G、chr9_306973：G/A、chr9_999337：G/A、chr9_1963317：T/C、chr2_173473：C/T、chr2_1811972：C/T、chr2_789901：C/A、chr2_280007：C/T,8 SNP molecular markers respectively, is shown in Table 4, the base sequence of 101bp before and after each SNP molecular marker is provided in Table 4, and the base variation exists at 51 bp.
Further, the SNP marker is numbered chr1_54231: the nucleotide sequence of A/G is shown as SEQ ID NO. 1;
the SNP marker is numbered chr9_306973: the nucleotide sequence of the G/A is shown as SEQ ID NO. 2;
The SNP marker is numbered chr9_999337: the nucleotide sequence of the G/A is shown as SEQ ID NO. 3;
the SNP marker is numbered chr9_1963317: the nucleotide sequence of the T/C is shown as SEQ ID NO. 4;
SNP markers numbered chr2_173473: the nucleotide sequence of the C/T is shown as SEQ ID NO. 5;
the SNP marker is numbered chr2_1811972: the nucleotide sequence of the C/T is shown as SEQ ID NO. 6;
The SNP marker is numbered chr2_789901: the nucleotide sequence of the C/A is shown as SEQ ID NO. 7;
The SNP marker is numbered chr2_280007: the nucleotide sequence of the C/T is shown as SEQ ID NO. 8.
The method for detecting the Nile-Lafei buffalo genotype by utilizing a molecular biology technology based on SNP molecular marker combination related to Nile-Lafei buffalo growth traits screened by whole genome association analysis provided by the embodiment of the invention comprises the following steps: and respectively designing amplification primers according to nucleotide sequences of two flanks of the 8 SNP molecular marker loci, amplifying by taking the DNA of the Nile-Lafei buffalo individuals as a template, sequencing the first generation of amplification products, and genotyping the Nile-Lafei buffalo individuals according to a sequencing result, thereby identifying the genotypes of the Nile-Lafei buffalo individuals to be detected.
Nile-Lafei buffalo growth trait whole genome association analysis
1. Population selection and growth trait phenotype data collection
And selecting 99 nieis-Lafei buffalo with healthy body and good development and 24 months of age, and collecting corresponding data related to the growth characters of the height, the body oblique length and the chest circumference.
2. Growth trait phenotype data collection and sorting and statistical analysis
Statistical analysis of the collected phenotypic data, including minimum, maximum, mean, standard deviation, coefficient of variation, was performed using SPSS20 software and the results are shown in table 1.
TABLE 1 statistical analysis of growth trait phenotypes of Nile-Lafei buffalo
3. Blood genomic DNA extraction and detection
5ML of jugular blood of each cattle is collected, heparin sodium is anticoagulated, the mixture is gently shaken to be uniformly mixed, the genomic DNA of the cattle is extracted by adopting a blood genomic DNA extraction kit (Tiangen), the DNA concentration is detected by a nucleic acid concentration detector, and the mass of the DNA sample is determined by agarose gel electrophoresis of 1 percent. The DNA sample that was detected was sent to Beijing-Nodezafirland Bioinformation technologies Inc. for genomic sequencing. And after quality inspection is qualified, the obtained product is used for resequencing and library building. The sequencing strategy was Illumina PE150.
4. Sequencing data quality control
The depth of resequencing is more than 5X, quality control is carried out on raw data obtained through sequencing, and filtering standards comprise: (1) removing reads containing the linker sequence; (2) removing PAIRED READS having an N content exceeding 10%; (3) And removing reads with the low-quality base number (quality less than or equal to 5) accounting for more than 50 percent. Relevant data statistics are carried out on CLEAN DATA obtained, wherein the sequencing data yield, the Q20 content, the Q30 content and the GC content are shown in table 2, and the high-quality CLEAN DATA data size is 101.75Gb.
TABLE 2 statistics of Nile-Lafei buffalo genomic DNA sequencing results
5. Genomic data alignment
Efficient high quality sequencing data was aligned to the reference genome by BWA bioinformatics software (parameters: mem-t 4-k 32-M). The average comparison rate of the population samples with reference to the genome download address ：ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/471/725/GCF_000471725.1_UMD_CASPUR_WB_2.0/GCF_000471725.1_UMD_CASPUR_WB_2.0_genomic.fna.gz is 99.49919192%, the average sequencing depth of the genome is 10.98767677X, and the average 1X coverage (coverage of at least one base) is 3.877676768%, and the information is shown in Table 3.
TABLE 3 Nile-Lafei buffalo genome data alignment statistics
6. Mutation detection, quality control and filtration
Detection of population SNPs was performed using SAMTOOLS software. Polymorphic sites in the population are detected using a bayesian model. The quality control method comprises the following steps: (1) Filtering SNPs with the quality value of more than 20 (error rate more than 1 percent); (2) If the distance between the two SNPs is detected to be within 5bp, removing the SNPs; (3) The depth of coverage of SNPs is between 1/3 and 5 times the average depth. The SNP obtained through the quality control preliminary detection is further subjected to strict quality control, and is used for guaranteeing the accuracy of mutation information detection. The SNPs obtained were filtered to obtain high quality SNPs under conditions of dp2, miss0.2, maf0.01, mapping quality < 40.0, MQRankSum < -12.5, readPosRankSum < -8.0, SOR > 3.0 and FISHER STRAND > 60.0, and finally 143575 SNP sites were retained for subsequent association analysis.
7. Whole genome association analysis
And carrying out association analysis by adopting GEMMA software and combining phenotype information and genome SNP information, and screening potential candidate SNPs when the-log 10 (P) >5 time difference is different and obvious through the associated significance (P-value). FIG. 1 is a high Manhattan diagram; FIG. 2 is a view of the body height Quantile-quantile; FIG. 3 is a Manhattan diagram of body diagonal length; FIG. 4 is a diagram of a body diagonal Quantile-quantile; FIG. 5 is a Manhattan diagram of a chest circumference; FIG. 6 is a diagram Quantile-quantile of the chest circumference; wherein the Manhattan diagram is a diagram of sorting the genetic marker effect value, namely the P value of the whole genome subjected to F test according to the physical position on the chromosome, the abscissa is the physical position on the chromosome of the genome, the ordinate is-log 10P, and the smaller the P value, the stronger the correlation is, and the larger the ordinate is. The horizontal dashed line in the Manhattan plot represents the level of significance, and when-log 10 (P) >5, the SNP is considered significantly associated with the trait. qq plot (Quantile-quantileplot) representing a distribution of actual P values and unassociated zero hypothesis expected P values for detecting the impact of population stratification and individual affinity on association analysis, if observed P values and expected P values occur only at the far right end of the distribution, indicating that the trait is not caused by population stratification. QQ plot is mainly used to estimate the difference between quantitative trait observations and predictions. The quantitative trait data obtained are normally distributed data. The X and Y axes of QQ plot in the GWAS study are primarily-lgPvalues representing the individual SNPs. The predicted line is a dashed line at a 45 angle from the origin. Carrying out character association analysis by adopting a mixed linear model, wherein the population genetic structure is used as a fixed effect, and the individual relationship is used as a random effect so as to correct the influence of the population structure and the individual relationship:
y＝Xα+Zβ+Wμ+e
Wherein y is a phenotypic trait, X is an indication matrix of a fixed effect, and alpha is an estimated parameter of the fixed effect; z is an indication matrix of SNP, and beta is the effect of SNP; w is an indicator matrix of random effects, μ is the predicted random individuals, e is the random residual, subject to e to (0, δe2).
8. Screening and extracting SNP molecular markers related to Nile-Lafei buffalo body height, body oblique length and chest circumference growth characteristics
And carrying out whole genome association analysis and screening 8 obviously associated SNP loci which are associated with the high, oblique body length and chest circumference growth traits of the Nile-Latifolian buffalo, and can be used for breeding the growth traits of the Nile-Latifolian buffalo. Details are shown in Table 4.
1 Locus was selected from 8 significantly associated SNPs for verifying the accuracy of the association results. The site is chr2_1811972: C/T.
TABLE 4 Nile-Lafei buffalo 8 molecular marker information with significant association to growth traits
9. Verification of Nile-Lafei buffalo SNP molecular markers
The Nile-Lafei buffalo 62 heads were selected. The jugular vein collects blood and extracts genomic DNA by using a blood genomic DNA extraction kit, and detects DNA concentration by using a nucleic acid concentration measuring instrument, and the mass of the DNA sample is measured by agarose gel electrophoresis of 1%. Amplification primers were designed for the 1 candidate SNP site selected. The primer sequences are as follows:
PCR amplification was performed using the above primers and the Nile-Lafei buffalo blood genomic DNA as a template. A50. Mu.L PCR reaction system was used: ddH 2O 19μL,Premix TaqTM. Mu.L, DNA template 2.0. Mu.L, primers (10. Mu. Mol/L for both upstream and downstream primers) 2.0. Mu.L each.
PCR reaction conditions:
Pre-denaturation at 95 ℃ for 4min; denaturation at 94℃for 10s, annealing for 30s (60 ℃) and extension at 72℃for 1min for 35 cycles; extending at 72℃for 5min.
The PCR amplified product was purified by Gel Extraction Kit kit from Shanghai Biotechnology Co., ltd, and specific steps are shown in the kit instruction. The PCR purified product obtained above is recovered and directly sent to Huada gene (Shenzhen) Biotech company for first generation sequencing. And (5) genotyping the individual according to the sequencing result. Genotyping results are shown in Table 5.
Table 5 Chr2_1811972 locus was at Nile-Lafei buffalo genotype frequencies and allele frequencies
It can be seen from Table 5 that the C allele frequency at mutation site chr2-181197 is significantly greater than the T allele frequency.
TABLE 6 verification of significant associations SNP of Nile-Lafei buffalo chest circumference
The chest circumference growth characteristic phenotype value is expressed by 'least square mean value +/-standard deviation', and different letters of the same-column data shoulder marks show that the difference is obvious (P < 0.05); the same letters of the shoulder marks or no letter labels indicate that the difference is not significant (P > 0.05); the SNP locus genotypes are sequentially arranged according to mutant type, heterozygous type and reference type.
The result shows that the candidate chr2_1811972 has obvious difference with the niry-Lafei buffalo chest circumference through verification, the CC type chest circumference is obviously larger than the CT and TT types, the CT and TT types are not obvious, and the growth is good when the chest circumference is large. The verification result is consistent with the whole genome association analysis result, which shows that the SNP molecular marker screened by the whole genome association analysis is accurate and reliable, and can be directly used as molecular marker auxiliary breeding for identifying growth characters such as the chest circumference of the Nile-Lafei buffalo and the like, and can be applied to production.
During breeding of the Nile-Lafei buffalo, the 8 SNP molecular markers can be utilized, amplification primers are respectively designed on nucleotide sequences of two lateral wings of the Nile-Lafei buffalo, jugular vein blood of the Nile-Lafei buffalo is collected when the Nile-Lafei buffalo is 3 months old, a blood genome DNA extraction kit is adopted to extract genome DNA, PCR amplification is carried out by taking individual DNA of the Nile-Lafei buffalo as a template, first generation sequencing of amplified products is carried out, genotyping is carried out according to a sequencing result, individuals with target character genotypes are reserved, the breeding period can be shortened, the breeding cost is reduced, the breeding efficiency is improved, and the breeding process is accelerated.
The foregoing is merely illustrative of specific embodiments of the present invention, and the scope of the invention is not limited thereto, but any modifications, equivalents, improvements and alternatives falling within the spirit and principles of the present invention will be apparent to those skilled in the art within the scope of the present invention.
Claims (5)
1. The SNP molecular marker combination related to the growth trait of Nile-Laffy buffalo screened based on whole genome association analysis is characterized by comprising at least one or a combination of the following 8 SNP markers, wherein the positions of the SNP molecular marker loci and mutation information are chromosome_physical positions: formal representations of reference genotype/variant genotype, respectively ：chr1_54231：A/G、chr9_306973：G/A、chr9_999337：G/A、chr9_1963317：T/C、chr2_173473：C/T、chr2_1811972：C/T、chr2_789901：C/A、chr2_280007：C/T.
2. The combination of SNP molecular markers associated with the growth trait of niri-raffia buffalo screened based on whole genome association analysis as set forth in claim 1, wherein the SNP marker is numbered chr1_54231: the nucleotide sequence of A/G is shown as SEQ ID NO. 1;
the SNP marker is numbered chr9_306973: the nucleotide sequence of the G/A is shown as SEQ ID NO. 2;
The SNP marker is numbered chr9_999337: the nucleotide sequence of the G/A is shown as SEQ ID NO. 3;
the SNP marker is numbered chr9_1963317: the nucleotide sequence of the T/C is shown as SEQ ID NO. 4;
SNP markers numbered chr2_173473: the nucleotide sequence of the C/T is shown as SEQ ID NO. 5;
the SNP marker is numbered chr2_1811972: the nucleotide sequence of the C/T is shown as SEQ ID NO. 6;
The SNP marker is numbered chr2_789901: the nucleotide sequence of the C/A is shown as SEQ ID NO. 7;
The SNP marker is numbered chr2_280007: the nucleotide sequence of the C/T is shown as SEQ ID NO. 8.
3. The method for detecting the genotype of the nii-raffia buffalo based on the SNP molecular marker combination related to the growth traits of the nii-raffia buffalo screened by the whole genome association analysis according to any one of claims 1 to 2, wherein the amplification primers are respectively designed according to the nucleotide sequences of two flanks of the 8 SNP molecular marker loci, the DNA of the nii-raffia buffalo individual is used as a template for amplification, the first generation of the amplified product is sequenced, and the nii-raffia buffalo individual is subjected to genotyping according to the sequencing result, so that the genotype of the nii-raffia buffalo individual to be detected is identified.
4. Use of a combination of SNP molecular markers related to the growth trait of niri-raffia buffalo screened based on whole genome association analysis as set forth in claim 1 in molecular marker-assisted breeding of niri-raffia buffalo growth traits.
5. Use of the method for detecting the genotype of niry-raffia buffalo by using molecular biological technique as defined in claim 3 in molecular marker assisted breeding of niry-raffia buffalo growth traits.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410150893.1A CN118109605B (en) | 2024-02-02 | 2024-02-02 | SNP molecular marker combination related to growth traits of Nile-Lafei buffalo and application |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410150893.1A CN118109605B (en) | 2024-02-02 | 2024-02-02 | SNP molecular marker combination related to growth traits of Nile-Lafei buffalo and application |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118109605A true CN118109605A (en) | 2024-05-31 |
| CN118109605B CN118109605B (en) | 2024-09-10 |
Family
ID=91217501
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410150893.1A Active CN118109605B (en) | 2024-02-02 | 2024-02-02 | SNP molecular marker combination related to growth traits of Nile-Lafei buffalo and application |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118109605B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113308554A (en) * | 2021-06-28 | 2021-08-27 | 广西壮族自治区水牛研究所 | SNP molecular marker related to bovine growth traits and application thereof |
| CN116397032A (en) * | 2023-01-30 | 2023-07-07 | 广西壮族自治区水牛研究所 | SNP molecular marker related to buffalo weight growth traits and application thereof |
| CN116837112A (en) * | 2023-07-12 | 2023-10-03 | 中国农业科学院兰州畜牧与兽药研究所 | SNP molecular marker related to yak growth traits and application thereof |
-
2024
- 2024-02-02 CN CN202410150893.1A patent/CN118109605B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113308554A (en) * | 2021-06-28 | 2021-08-27 | 广西壮族自治区水牛研究所 | SNP molecular marker related to bovine growth traits and application thereof |
| CN116397032A (en) * | 2023-01-30 | 2023-07-07 | 广西壮族自治区水牛研究所 | SNP molecular marker related to buffalo weight growth traits and application thereof |
| CN116837112A (en) * | 2023-07-12 | 2023-10-03 | 中国农业科学院兰州畜牧与兽药研究所 | SNP molecular marker related to yak growth traits and application thereof |
Non-Patent Citations (2)
| Title |
|---|
| SAHER ISLAM等: "Population demographic history and population structure for Pakistani Nili-Ravi breeding bulls based on SNP genotyping to identify genomic regions associated with male effects for milk yield and body weight", PLOS ONE, vol. 15, no. 11, 24 November 2020 (2020-11-24), pages 0242500 * |
| 黄萌;郑海英;杨春艳;黄加祥;鄢胜飞;李舒露;于农淇;李孟琪;尚江华;: "摩拉水牛和尼里-拉菲水牛PRKAA2基因SNPs检测及遗传多样性分析", 中国畜牧兽医, no. 05, 20 May 2018 (2018-05-20), pages 1274 - 1282 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118109605B (en) | 2024-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115029451B (en) | Sheep liquid phase chip and application thereof | |
| CN117385045B (en) | Application of goat AMHR2 gene as lambing number-associated molecular marker | |
| CN112126690B (en) | SNP molecular marker influencing thoracic vertebra number character of sheep and application | |
| US11542562B2 (en) | Single nucleotide polymorphism marker related to Chinese horse short stature trait and use thereof | |
| CN117757952A (en) | Sheep whole genome 45K SNP liquid phase chip and application thereof | |
| CN118109605B (en) | SNP molecular marker combination related to growth traits of Nile-Lafei buffalo and application | |
| CN115927667A (en) | Molecular marker and primer related to pig intramuscular fat character and application of molecular marker and primer | |
| CN118086517B (en) | SNP molecular marker combination for detecting propagation traits of Mola buffalo and application | |
| CN117904317B (en) | SNP molecular marker combination for detecting propagation traits of Nile-Lafei buffalo and application | |
| CN109750106B (en) | Long-chain non-coding RNA combination for evaluating bull sperm motility and detection method and application thereof | |
| CN117701727B (en) | SNP molecular marker combination related to size and birth weight of Mora buffalo based on whole genome sequencing screening and application | |
| CN118957090A (en) | SNP molecular marker for detecting correlation of calving interval reproduction traits of Mola buffalo and application of SNP molecular marker | |
| CN106755370B (en) | Method for detecting sheep FTH-1 gene single nucleotide polymorphism by using PCR-RFLP and application thereof | |
| CN115807100B (en) | SNP molecular marker related to abdominal fat rate of broiler chickens and application thereof | |
| CN113930521A (en) | C7H15orf39 gene SNP molecular marker associated with boar semen quality traits and application | |
| CN118421803A (en) | SNP (Single nucleotide polymorphism) marker primer pair on pig chromosome 1 and related to pig vulva width trait and application thereof | |
| CN117867133A (en) | Application of PDGFD gene upstream SNP marker in sheep variety tail type selection | |
| CN117867132A (en) | Application of downstream SNP marker of BMP2 gene in tail type selection of sheep variety | |
| CN118562966A (en) | SNP molecular marker related to silurus meridionalis antibacterial septicemia trait and application thereof | |
| CN116640860A (en) | Molecular marker related to sheep residual feed intake traits and application thereof | |
| CN117904307A (en) | SNP molecular marker related to lean meat percentage of pig on chromosome 6 of pig and application thereof | |
| CN116694779A (en) | SNP molecular marker related to Duroc pig eye muscle area and application thereof | |
| CN115820875A (en) | Molecular marker for evaluating pig eye muscle area and screening method and application thereof | |
| CN118581233A (en) | Commercial pig liquid phase chip and preparation method and application thereof | |
| CN117004739A (en) | Method for judging black spot characters of red tilapia in wintering period |
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 |