CN112652362B - Method for screening bovine plateau hypoxia adaptive gene ALDOC and functional molecular marker and application thereof - Google Patents
Method for screening bovine plateau hypoxia adaptive gene ALDOC and functional molecular marker and application thereof Download PDFInfo
- Publication number
- CN112652362B CN112652362B CN202011363669.9A CN202011363669A CN112652362B CN 112652362 B CN112652362 B CN 112652362B CN 202011363669 A CN202011363669 A CN 202011363669A CN 112652362 B CN112652362 B CN 112652362B
- Authority
- CN
- China
- Prior art keywords
- cattle
- gene
- altitude
- screening
- flk
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
- G16B30/20—Sequence assembly
-
- 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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B10/00—ICT specially adapted for evolutionary bioinformatics, e.g. phylogenetic tree construction or analysis
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
- G16B20/20—Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
- G16B20/30—Detection of binding sites or motifs
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
- G16B40/10—Signal processing, e.g. from mass spectrometry [MS] or from PCR
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Biotechnology (AREA)
- Medical Informatics (AREA)
- Evolutionary Biology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Analytical Chemistry (AREA)
- Bioinformatics & Computational Biology (AREA)
- Theoretical Computer Science (AREA)
- Molecular Biology (AREA)
- Organic Chemistry (AREA)
- Genetics & Genomics (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Signal Processing (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Crystallography & Structural Chemistry (AREA)
- Animal Behavior & Ethology (AREA)
- Bioethics (AREA)
- Immunology (AREA)
- Microbiology (AREA)
- Evolutionary Computation (AREA)
- Epidemiology (AREA)
- Databases & Information Systems (AREA)
- Biochemistry (AREA)
- General Engineering & Computer Science (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
The invention provides a method for screening a bovine plateau hypoxia adaptive gene ALDOC and a functional molecular marker and application thereof, belonging to the technical field of animal molecular breeding. The method selects and distributes the cattle varieties at high altitude and low altitude from the aspects of cattle genome selection, genetic adaptation and the like, integrates three genome selection signal analysis methods of FLK, hapFLK and XPEHH by adopting an SNP chip, compares the genome selection signals between the cattle at low altitude and in the cattle varieties, identifies candidate genes adapting to the extreme environment of high-altitude hypoxia through a combination screening strategy and bioinformatics analysis, further excavates potential functional molecular markers through selection signal verification and re-sequencing data analysis, and establishes a corresponding detection method. The invention provides scientific basis and simple and easy detection technology for molecular breeding such as breeding of plateau hypoxia special cattle variety; meanwhile, the method has important significance and value for the protection, evaluation and utilization of local cattle variety genetic resources.
Description
Technical Field
The invention belongs to the technical field of animal molecular breeding, and particularly relates to a method for screening a bovine plateau hypoxia adaptive gene ALDOC and a functional molecular marker and application thereof.
Background
The information in this background section is only for enhancement of understanding of the general background of the invention and is not necessarily to be construed as an admission or any form of suggestion that this information forms the prior art that is already known to a person of ordinary skill in the art.
Cattle are ruminants that live in different environments with the ability to efficiently convert low quality feed into high energy fats, meat and milk (Elsik et al, 2009). Under the action of rumen microorganisms, cattle digest roughage into fermentation products such as acetic acid, butyric acid and propionic acid (Russell and van soest, 1984). The high altitude of the Tibet plateau impairs rumen fermentation and increases basal metabolic rate (Qiao et al, 2013). Oxygen is a key substance for maintaining the metabolism and life of animal bodies and is necessary for normal life activities. Hypoxic environments produce a series of physiological and genomic changes that allow humans and local animals to adapt to local extreme high altitude climates (O' brien et al, 2020). Since the domestication of cattle, cattle become indispensable animal resources in plateau areas and play an important role in providing meat, milk and services. Therefore, the environmental adaptation mechanism of cattle attracts people's attention.
Yaks (Bos grunniens) and common cattle (Bos taurus) differentiate approximately 500 ten thousand years ago, but retain a high degree of homology at the genomic level (Qiu et al, 2012). Yak mainly lives in high altitude areas, especially the Xingdushushu-Himalayan area and the Qinghai-Tibet plateau area with 3000 + 5500 meters. From a physiological perspective, yaks have greater metabolic capacity than low altitude cattle (Wang et al, 2011), lower cardiopulmonary capacity, greater cardiopulmonary volume (Wiener et al, 2003). Yak and tibetan varieties have adapted to the plateau environment by expanding the number of pulmonary, bronchial arterial smooth muscles, erythrocytes, hemoglobin and hematocrit (Li et al, 2006; Ma et al, 2011). These studies lay the foundation for understanding the molecular mechanisms of human and animal adaptation to high altitude environments. Historically, ape za cows were bred 80-100 years ago by introduction of the tumor bovine (Bos indicus) ancestor from india, dan and nipol, and by cross breeding with native cattle in china (Guan et al, 2017). Tibet cattle and Japanese Ka cattle also live in plateau with altitude of more than 3500 m, and these cattle have been well adapted to high altitude anoxic environment through long-term natural selection. They also have strong foraging capacity, rough feeding resistance, high energy metabolism and strong cardiopulmonary function to ensure normal circulation and transport of blood in vivo (Dolt et al, 2007; Wang et al, 2011), which provides an ideal model for understanding the plateau hypoxia adaptive mechanism.
High-altitude areas are characterized by oxygen deficiency (the oxygen content in the air is only 50% -60% of that in plain areas), low annual average temperature (from-1 ℃ to-5 ℃), short growing season (from 6 months to 9 months) and seasonal changes in feed supply (Shao et al, 2010), resulting in lack of high-quality pasture. Therefore, under the condition of insufficient plateau nutrition supply, how the cattle provide metabolic function to ensure the growth, development, survival and production performance is important for the hypoxia metabolism regulation and adaptation of plateau cattle organisms. In plateau regions, the energy supply is insufficient due to hypoxia, lack of feed and nutrition, and the animal body often obtains the required energy by improving glycolysis or improving the utilization rate of substrates.
The major energy source for most cells is glucose, from which Adenosine Triphosphate (ATP) is produced by glycolysis and/or oxidative metabolism. Since a deficiency in glucose causes an increase in the intracellular adenosine phosphate/adenosine triphosphate ratio, adenylate activated protein kinase (AMPK) is activated, restoring energy balance by inhibiting synthetic processes (e.g., protein or lipid synthesis) and promoting catabolic processes (e.g., glycolysis) (Xiao et al, 2017). Glycolysis is a major catabolic process in response to energy stress, the metabolic pathway by which most of the pyruvate metabolism is converted from glucose to lactate. Although normal non-proliferating cells can only undergo glycolysis under anaerobic conditions, some rely primarily on glycolysis to produce ATP and the basic unit of biosynthesis even under aerobic conditions, the so-called "aerobic glycolysis" or "wobbe effect" (Vander Heiden et al, 2009). The fructose 1.6-bisphosphate Aldolase (ALDOC) gene, aldolase, can act as a monitoring system, sensing a decrease in glucose supply before the cell energy status decreases, and activating AMPK in an AMP/ADP independent manner. Aldolase (ALDOC) is an ideal sensor for glucose utilization (Zhang et al, 2018). Under hypoxic conditions, ALDOC upregulation expression promotes glycolysis to produce ATP to meet the needs of animal energy metabolism (Leiherer et al, 2014). Under hypoxic conditions, ALDOC upregulation expression promotes glycolysis to produce ATP to meet the needs of animal energy metabolism (Leiherer et al, 2014).
Disclosure of Invention
Aiming at the prior art, the invention provides a method for screening a bovine plateau hypoxia adaptive gene ALDOC and a functional molecular marker and application thereof. The method selects the cattle varieties distributed at high altitude (the altitude is more than 1800 m) and low altitude (the altitude is less than 1500 m) from the aspects of cattle genome selection, genetic adaptation and the like, integrates three genome selection signal analysis methods of FLK, hapFLK and XPEHH by adopting a cattle 777K high-density SNP chip, compares the genome selection signals between the cattle at high altitude and in the cattle varieties, identifies candidate genes adapting to the extreme environment of high-altitude hypoxia through a combined screening strategy and bioinformatics analysis, further excavates potential functional molecular markers through selection signal verification and sequence re-testing data analysis, and establishes a corresponding detection method. The invention provides scientific basis and simple and easy detection technology for molecular breeding such as breeding of plateau hypoxia special cattle varieties; meanwhile, the method has important significance and value for the protection, evaluation and utilization of local cattle variety genetic resources.
The invention aims to provide a method for screening a bovine plateau hypoxia adaptive gene and a functional molecular marker.
The second object of the present invention is to provide the application of the above method.
In order to realize the purpose, the invention relates to the following technical scheme:
in a first aspect of the present invention, a method for screening a bovine plateau hypoxia adaptive gene and a functional molecular marker is provided, the method at least comprising:
DNA samples of different altitude cattle varieties are detected and analyzed based on the SNP chip;
the method is characterized in that analysis is carried out based on FLK, hapFLK and XPEHH genome selection signal analysis methods, single nucleotide polymorphic SNPs and candidate genes subjected to significant positive selection are screened, and plateau hypoxia adaptive candidate genes and SNP molecular markers are screened by integrating gene function annotation.
Further, based on the method, the ALDOC gene is screened and identified as a key gene of the bovine plateau hypoxia adaptability; furthermore, the bovine plateau hypoxia adaptive molecular marker also comprises SNP located in the gene, including but not limited to g.20588331T > G (rs133198943A > C); g.20590099G > A (c.288C > T, aa96F > F); g, 20590664T > C; g.20591903C > A (g.20591903G > T); the single nucleotide polymorphic site g.20591903C > A has the strongest selection signal and is a functional molecular marker verified by experiments, wherein the bovine plateau hypoxia adaptive genotype is TT type or complementary sequence AA type.
In a second aspect of the present invention, there is provided the use of the above method for screening a bovine individual (population, line or breed) suitable for high altitude hypoxia survival.
The beneficial technical effects of one or more technical schemes are as follows:
the technical scheme aims at the genetic characteristic of plateau hypoxia adaptation of the cattle variety, selects local cattle varieties distributed at high altitude and low altitude from the aspects of genome selection and genetic adaptation, integrates methods and strategies of various genome selection signals, bioinformatics, promoter activity analysis, gene expression analysis and the like, efficiently and accurately screens the key gene ALDOC and the functional molecular marker which adapt to the plateau hypoxia, and establishes a CRS-PCR-RFLP method of a creative enzyme cutting site by introducing mutation to detect the genotype, has reasonable design, and has the characteristics of high accuracy, simple and convenient application operation and low cost according to the detection method designed by the key gene and the functional marker;
by adopting the method provided by the technical scheme, individuals with high plateau hypoxia adaptability can be effectively screened out, and the method has important significance for breeding work and genetic improvement of cattle varieties in high-altitude regions. The technical scheme is that the molecular breeding technology is well applied to production practice once, can provide a technology for the conservation and utilization of bovine germplasm resources and the research of breed provenance, greatly saves the breeding cost and the special germplasm cultivation time, and obtains good economic and social benefits, thereby having good practical application value.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, are included to provide a further understanding of the invention, and are included to illustrate an exemplary embodiment of the invention and not to limit the invention.
FIG. 1 shows the consensus genes identified for FLK, hapFLK and XPEHH in the present examples (Wayne diagram).
FIG. 2 shows the selection signal distribution of chromosome 19 hapFLK and XPEHH between high and low altitude bovine species in an example of the present invention.
FIG. 3 shows the SNP and haplotype evolution tree of the Chinese local cattle variety (25) differential selection region (Chr19:20.537-20.664Mb, including ALDOC gene) in the present example.
FIG. 4 is an analysis of allele frequencies of ALODC gene structure and potential functional variations thereof in different bovine populations according to an embodiment of the present invention.
FIG. 5 shows the effect of SNP (g.20591903C > A) on the promoter activity of ALDOC gene in examples of the present invention;
note: denotes P < 0.05.
FIG. 6 shows the difference in expression of liver tissue ALDOC of individual with different genotypes of SNP (g.20591903C > A) in the example of the present invention.
FIG. 7 shows the sequencing results of SNP (g.20591903G > T) individuals with different genotypes in the examples of the present invention.
Detailed Description
It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the disclosure. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the present application. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof, unless the context clearly indicates otherwise.
At present, few molecular biology researches on plateau hypoxia adaptability of local cattle are carried out, and no method for screening plateau hypoxia adaptability local cattle exists.
In view of the above, in one embodiment of the present invention, a method for screening a bovine plateau hypoxia adaptive gene and a functional molecular marker is provided, the method at least comprising:
DNA samples of different altitude cattle varieties are detected and analyzed based on the SNP chip;
the method is characterized in that the method is based on FLK, hapFLK and XPEHH genome selection signal analysis methods for analysis, single nucleotide polymorphisms SNPs and candidate genes subjected to remarkable positive selection are screened, and gene function annotation is integrated to screen plateau hypoxia adaptive candidate genes and SNP molecular markers.
In another embodiment of the present invention, the method for screening the hypoxia adaptive gene and the functional molecular marker of the bovine plateau comprises:
s1, collecting and DNA extracting samples of cattle varieties at different altitudes;
s2, detecting and analyzing the SNP chip;
s3, analyzing FLK and hapFLLK genome selection signals;
s4, XPEHH genome selection signal analysis;
s5, identifying a positively selected genetic variation and a differential selection area;
s6, screening the candidate genes based on a screening strategy of the candidate genes;
s7, verifying a positive selection signal of the difference selection area;
and S8, identifying potential functional variation in the positive selection gene by using the remeasurement data.
Wherein, the steps S3 and S4 have no precedence, so the step sequence can be S3-S4 or S4-S3;
in another embodiment of the present invention, the step S1 includes:
s1.1, selecting common cattle and tumor cattle varieties at different altitudes, and mixed varieties of the common cattle and the tumor cattle;
s1.2, blood of a cow is collected, and DNA in blood tissues is extracted.
In another embodiment of the present invention, the step S2 includes:
s2.1, analyzing the sample by using the SNP chip, and genotyping;
s2.2, filtering the SNP data, and further analyzing the remaining SNPs meeting the requirements;
s2.3 haplotype was constructed for each chromosome.
In another embodiment of the present invention, the step S3 includes:
s3.1, respectively comparing selection signal values corresponding to SNPs of higher altitude and lower altitude groups based on an FLK and hapFLK whole genome selection signal analysis method, and acquiring difference selection signals of the SNPs between the two groups;
s3.2, respectively operating FLK and hapFLK programs by using genotype data to obtain selection signals of SNPs in the variety, and taking Nelore as a distant population; in FLK analysis Nelore is defined as the distant population, and let K10 and nFit 20;
s3.3, constructing an evolutionary tree of the whole genome and a local genome region for the selected region by using hapFLK results and Python and R scripts;
s3.4 the p-value of hapFLK values were calculated by fitting a standard normal distribution of whole genomes in R.
In another embodiment of the present invention, the step S4 includes:
s4.1, estimating XPEHH values between high-altitude and low-altitude cattle varieties, and defining the relation between physical distance and genetic distance by using 1Mb ≈ 1cM in cattle genome;
s4.2 construction of haplotypes based on ReqDel-FAST.
In another embodiment of the present invention, the step S5 includes:
based on the results of FLK, hapFLK and XPEHH selection signal analysis, SNPs with positive selection signal values at the top 0.1% were screened separately.
In another embodiment of the present invention, the step S6 includes:
s6.1, defining the upstream and downstream 50Kbp of the SNPs obtained by identification as selected areas, and connecting the overlapped areas in series to obtain differential selection areas;
s6.2, annotating a reference gene in the difference selection area or overlapped with the area by using the UCSC online website;
s6.3, screening out the FLK, hapFLK and XPEHH selection signals which are compared among high and low altitude groups, ranking 10-bit differential selection regions, respectively identifying candidate genes (such as ALDOC genes), and carrying out functional annotation.
In another embodiment of the present invention, the step S7 includes:
s7.1, respectively constructing a whole genome and a local selection region evolutionary tree based on FLK and hapFLK analysis results, and analyzing the selection condition of the candidate gene in each cattle variety;
s7.2 the results of the comparisons between groups were further validated based on the results of the in-breed FLK and hapFLLK selection signal analysis.
In another embodiment of the present invention, the method of step S8 includes:
s8.1, acquiring double-end Illumina re-sequencing data, and respectively comparing the acquired data with a reference genome (UMD 3.1);
s8.2, filtering to obtain a high-quality mapping sequence, and detecting genetic variation from the high-quality mapping sequence;
s8.3, extracting genetic variation from a specific gene region, wherein the extracted regions are an upstream region and a gene region of a candidate gene;
s8.5, respectively calculating the genotype and allele frequency of each genetic variation;
s8.6 through functional annotation, the SNP sites of potential functional variation are determined.
In another embodiment of the invention, a molecular marker for bovine plateau hypoxia adaptability, which is obtained by identification based on the method, is provided, and specifically, the ALDOC gene is a key gene of bovine plateau hypoxia adaptability and can be used as a molecular marker for bovine plateau hypoxia adaptability; meanwhile, the bovine plateau hypoxia adaptive molecular marker also comprises four mononucleotide polymorphic sites positioned on the gene, such as: g.20588331T > G (rs133198943A > C); g.20590099G > A (c.288C > T, aa96F > F); g, 20590664T > C; g.20591903C > A (g.20591903G > T). The SNP (g.20591903C > A) is a functional molecular marker verified by experiments, wherein the bovine plateau hypoxia adaptive genotype is TT type or complementary sequence AA type.
In another embodiment of the present invention, a kit for detecting the above molecular marker is provided, which can be used to screen bovine individuals (population, strain or breed) suitable for high altitude hypoxia survival; more specifically, the kit comprises a primer for detecting the single nucleotide polymorphic site; the single nucleotide polymorphic site comprises g.20588331T > G (rs133198943A > C); g.20590099G > A (c.288C > T, aa96F > F); g, 20590664T > C; any one or more of g.20591903C > A (g.20591903G > T);
in another embodiment of the invention, the invention provides a kit for detecting a single nucleotide polymorphic site g.20591903G > T, wherein the kit at least comprises the following primers:
F:5‘-GTAATTGTTTACGGTGACGC-3’(SEQ ID NO.7)；
R:5‘-GGCCTTGTTCTGATTCCTGC-3’(SEQ ID NO.8)。
in yet another embodiment of the present invention, there is provided the use of the above method for screening a bovine individual (population, line or breed) suitable for high altitude hypoxia survival.
In another embodiment of the present invention, the application method specifically comprises:
extracting blood DNA of different cattle individuals;
the kit for detecting the molecular marker is utilized to identify the genotype of the marker by a CRS-PCR-RFLP method, and identify individuals with the plateau hypoxia adaptive specific molecular marker.
In yet another embodiment of the invention, the molecular marker is a single nucleotide polymorphic site g.20591903G > T; in this case, the kit comprises at least the following primers and BstUI restriction enzyme;
F:5‘-GTAATTGTTTACGGTGACGC-3’(SEQ ID NO.7)；
R:5‘-GGCCTTGTTCTGATTCCTGC-3’(SEQ ID NO.8)。
the invention actually establishes a method for screening plateau hypoxia adaptability local cattle by using a genome direct sequencing technology. According to the invention, high-density SNP chip analysis is carried out on different varieties of local cattle, a selection signal analysis method is combined, specific SNP sites related to plateau hypoxia are found, and plateau hypoxia adaptability of the cattle can be judged by verifying single SNP sites or combination of SNP sites. The method has important significance for the molecular breeding development of cattle in plateau areas.
The present invention will be further described with reference to the following examples, but the present invention is not limited thereto. It should be understood that these examples are for illustrative purposes only and are not intended to limit the scope of the present invention. The test methods in the following examples, which are not specified under specific conditions, are generally carried out under conventional conditions.
Examples
1. Acquisition and DNA extraction of high-low altitude cattle variety samples
25 domestic local cattle breeds distributed at different altitudes and 1 foreign breed are selected as distant reference populations, and 352 cattle are selected, including common cattle, rumen cattle and mixed breeds of the common cattle and the rumen cattle are respectively extracted for blood DNA. Wherein, the cattle breeds distributed in the area with the altitude lower than 1500 m form a low altitude group (LA), and the cattle breeds distributed in the area with the altitude higher than 1800 m form a high altitude group (HA). The cattle breed and grouping information are shown in table 1:
TABLE 1 cattle breed information and grouping
2. Genotyping Using Illumina high Density Bovine HD 777KSNP chips
Genotyping is carried out by using the chip, and the genotype result in AB format is extracted from the genotype detection result file Finalreport.txt by using the python program PEDDA _ ROW. The chip had a total of 777,962 Single Nucleotide Polymorphisms (SNPs), and SNPs data were filtered using Plunk1.9 (http:// zzz. bwh. harvard. edu/plink 2.shtml) software to remove first 40,497 SNPs on X, Y and mitochondrial chromosomes and SNPs that did not map uniquely to UMD 3.1. Wherein the SNP detection rate is lower than 90 percent or the SNP with the minimum allele frequency of less than 0.05 is removed, and the residual 702,622 autosomal SNPs are subjected to subsequent analysis after filtration.
3. High-low altitude group and variety internal genome FLK and hapFLK selection signal analysis
Using hapflk software (https://forge-dga.jouy.inra.fr/projects/hapflk/) A method for analyzing FLK and hapFLlk whole genome selection signals respectively compares selection signal values corresponding to SNPs of a high-altitude group and a low-altitude group to obtain difference selection signals of the SNPs between the two groups. In addition, the FLK and hapFLK programs were run using genotype data from 25 cultivars, respectively, to obtain selection signals for SNPs within the cultivars, and Nelore was used as the distant population. In the FLK analysis, Nelore is defined as the distant population, and K is 10 and ntit is 20. hapFLK results were used along with modified Python and R scripts (local _ reanalds.py and local _ trees.r) to construct evolutionary trees of whole and local genomic regions for selected regions. And the p-value of the hapFLK value was calculated by fitting a standard normal distribution of the whole genome in R.
4.XPEHH genomic selection signal analysis
XPEHH was applied in the Selscan software to estimate XPEHH values between high and low altitude varieties. The XPEHH values were normalized in each set of comparisons to have mean and unit squared error. We used 1Mb ≈ 1cM genetic relationships in the bovine genome. And haplotype construction was performed using ReqDel-FAST in plink1.9 and fastPHASE1.4.
5. Identification of positively selected genetic variants and differentially selected regions
The FLK, hapFLK and XPEHH selection signal analysis method is used for screening the top 0.1% of SNPs subjected to significant positive selection respectively compared with high-altitude and low-altitude cattle groups, and 701, 702 and 700 SNPs are identified respectively.
6. Differentially selected regions and identification of candidate genes
The upstream and downstream 50Kbp of SNPs identified above were defined as selected regions, overlapping regions were concatenated to obtain differential selection regions, and 351, 32 and 99 differential selection regions were identified for FLK, hapFLK and XPEHH, respectively. Finally, using online UCSC (https://genome.ucsc.edu/cgi-bin/hgTables) Identification of reference genes within or overlapping the differentially selected region resulted in the identification of 341, 50 and142 candidate positive selection genes, likely associated with high altitude adaptation. Among these, three methods were applied to identify 2 common (overlapping) genes; FLK and hapFLK identified 9 consensus genes; FLK and XPEHH identified 23 common genes; hapFLK and XPEHH identified 5 common genes in total, including ALDOC, SPAG5, KIAA0100, PIGS, SDF2 (fig. 1). Further, screening out FLK, hapFLK and XPEHH selection signals which are compared among high and low altitude groups and selecting the difference selection areas which are 10 bits before the ranking, respectively identifying candidate genes and carrying out function annotation. Among them, hapFLK values of the differential selection region (Chr19:20.537-20.664Mb) where ADLOC gene was found were ranked first; moreover, the gene was focused on finding that the gene is associated with anaerobic glycolysis. We performed subsequent analyses using the selected region and the included ALDOC gene as candidate genes.
ADLOC gene is located on chromosome 19, and the results of genome selection signal analysis of hapFLK and XPEHH are shown in FIG. 2. In hapFLK analysis, the differential selection region has 56 SNPs in total, 11 of which are located in the first 0.1% (hapFLK ≧ 17.2), and in particular SNP (rs133198943) located in the ALODC, which has the strongest selection signal value (hapFLK ═ 20.51, P ═ 3.97E-04). In the XPEHH analysis, only 1 SNP (rs136098191) ranked the XPEHH signal value (XPEHH 3.51, P7.77E-04) at the top 0.1% (XPEHH ≧ 3.386).
7. Positive select signal verification of differentially selected regions
Respectively constructing a whole genome and a local selection region evolutionary tree by utilizing FLK (marker kinase) and hapFLK (marker kinase) analysis results in 25 varieties, analyzing the selection condition of the gene in each cattle variety, and analyzing the result as shown in figure 3, wherein in the FLK analysis in the varieties, the ALODC gene is selected most strongly in Tibet cattle (Tibet) of the plateau variety (P is 1.40E-41); the other plateau cattle variety, Aripinza broken cattle (Apeijiaza), is also strongly positively selected (P ═ 4.4E-09). In hapFLK analysis, three plateau breeds, tibetan (P7.7E-83), camu (Shigatse Humped; P4.9E-07), and apenema (P9.8E-08), were strongly selected, with the strongest signal. The results of the intravarietal selection signal analysis further confirmed the results of the comparisons between groups.
8. Identification of potential functional variations within positively selected genes using remeasured data
Double-ended Illumina re-sequencing data (data Nos.: PRJNA285834, PRJNA422979 and PRJNA396672) were downloaded from the NCBISRA database. The data contains 20 of 29 Tibet cattle and low-altitude cattle varieties (Kazakh cattle and Mongolia cattle), 10 of yaks, 14 of ruminants, 18 of wild cattle and 22 of longnius cattle. The downloaded data were aligned to the reference genome (UMD3.1) using BWA-MEM software, respectively. High quality mapping sequences were filtered using Samtools and design parameters "-view-f 4-q 20" and duplicate reads were removed with the "rmdup" parameter. Genetic variation was detected from high quality mapped sequences using GATK software. The bam file was sorted and the repeated sequence was masked using Picard (http:// broadlisting. github. io/Picard /). And extracting genetic variation from a specific gene region by using Vcftools software, wherein the extracted region is 2000Kbp upstream of the ATG of the ALDOC gene and a gene region. As a result, we identified 275 genetic variations in the population and calculated the frequency of each genotype and allele separately. By functional annotation, a total of 4 potential functional variants were identified (g.20588331T)>G，rs133198943A>C；g.20590099G>A，c.288C>T，aa96F>F；g.20590664T>C；g.20591903C>A) They have significant differences in the high and low altitude groups (fig. 4), and are potential functional variations associated with hypoxia adaptation in cattle plateau. Among them, SNP (g.20591903C)>A) In the promoter region, using AnimalTFDB (http://bioinfo.life.hust.edu.cn/AnimalTFDB/#！/) The transcription factor binding site prediction was carried out on the wild type (acggtgactcGccttcatc) and mutant sequence (acggtgactcTccttcatc), and the SNP g.20591903C was found>A(g.-1428G>T; a of translation initiation ATG +1) mutation caused a great change in transcription factor binding site, and as shown in table 2, 5 transcription factor (FOS, HDAC2, YBX1, YBX2, AP1) binding sites were deleted while 1 transcription factor (FOXO3) binding site was added. Many of the transcription factors described above are associated with hypoxia adaptation. For example, FOXO3 has been reported to be associated with hypoxia-induced endothelial apoptosis (Hu et al, Cell Signal,2018,51: 233-242). Research has shown that hypoxia inhibits hydroxylation and degradation of the FOXO3 protein, resulting in FOXO3 in the renal tubulesAccumulation and activation in epithelial cells. Hypoxia-activated HIF-1 is involved in the activation of FOXO3 and protects the kidney. In hypoxic kidneys, the stress-responsive transcription factor FOXO3 can be activated to adapt to hypoxic conditions, slowing the progression of chronic renal disease (Li et al, J Clin Invest,2019,129(6): 2374-. Hifs alters metabolic pathways by promoting anaerobic glycolysis and inhibiting oxidative phosphorylation, hypoxia activates Hifs to confer an anaerobic metabolic program on small oocytes, hypoxia induces oocyte dormancy by expressing FOXO3 and Hifs (Shimamoto et al PNAS, 2019, 116 (25): 12321-12326). FOXO 3-interfered mice soon lost fertility due to over-activation of immature oocytes (Castrillon et al, Science, 2003, 301: 215-. FOXO3 functions in reducing mitochondrial mass and oxygen consumption by inhibiting a series of nuclear-encoded mitochondrial genes (Jensen et al, EMBO J.2011, 30: 4554-. The above studies suggest that the transcription factor FOXO3 is associated with anaerobic glycolysis and hypoxia adaptation, SNP (g.20591903C)>A) It is possible to alter the expression of ALDOC by affecting the binding of transcription factors to ALDOC.
TABLE 2 ALDOC Gene SNP (g.20591903C > A or G > T) changes in wild-type and mutant sequence transcription factor binding sites
9. Effect of SNP (g.20591903C > A) on the Activity of the promoter of the ALDOC Gene
(1) Construction of recombinant plasmid: a pair of primers ALDOC-SNP20591903 is designed by taking DNA of SNP (g.20591903C > A) wild type and mutant type individuals as templates.
5'-GGGGTACCtaccttctccatccccctct-3', SEQ ID NO.1, capital letters are restriction enzyme cutting sites and protective basic groups of restriction enzyme KpnI;
R5'-CCGCTCGAGcgatgaggtggagtgactga-3', SEQ ID NO.2, capital letters are restriction sites and protective bases of restriction enzyme XhoI
Respectively amplifying fragments containing SNP (g.20591903C > A) by PCR (polymerase chain reaction), wherein the length of the amplified fragments is 327bp, constructing a wild type plasmid (pGL3-aldoc-GG) and a mutant type plasmid (pGL3-aldoc-TT) containing SNP (g.20591903C > A), and verifying the accuracy of the sequences by using product sequencing.
(2) Cell culture and transient cell transfection: the culture medium used for culturing the human liver cancer cells (HepG2) and fetal bovine fibroblasts (BFF) is a DMEM complete culture medium containing 7% FBS and 1% double antibody, the frozen cells are revived at 37 ℃, placed in a 37 ℃ and 5% incubator for culturing for 48-72 hours and then subcultured. The cultured healthy cells are subcultured to a 24-well plate and cultured for 24-48h, and cell transfection is carried out when the cell density reaches 70% -90%, wherein the method comprises the following steps:
preparing a transfection solution:
solution (50. mu.L of Opti-MEM + 1.8. mu.L of lipofectamine 2000;
solution two 50 uL Opti-MEM +800ng pGL3-basic recombinant plasmid/no-load plasmid +20ng internal reference plasmid pRL-TK and standing at room temperature for 5 min;
the solution (III) and the solution (II) are mixed in a centrifugal tube of 1.5ml and incubated for 20 min;
the original medium in the 24-well plate was discarded, washed gently twice with Opti-MEM, and 400. mu.L of Opti-MEM was added within 20min of incubation of solution I and solution II. And placing the solution (c) in a 24-well plate, culturing for 5h in a 5% incubator at 37 ℃, discarding the transfection solution, and replacing with a complete culture medium containing serum to continue culturing for 48 h. Luciferase activity was measured.
Experimental group (transfection of recombinant plasmids with different target sequences) and control group (transfection of PGL3-basic empty plasmid) were transfected separately by the above method, and BFF and HepG2 cells were transfected separately for each plasmid, and transfection was repeated 3 times.
(3) Detecting the luciferase activity: after culturing the cells for 48 hours, the 24-well plate was taken out and the medium was discarded, and the cells were washed twice with PBS. Then, 100 μ L of 1 × cell lysate PLB was added to a 24-well plate, the plate was gently shaken to accelerate cell lysis, and the cell lysis was observed in a microscope, after lysis for 15-30min, the lysed cell suspension was removed and stored in 1.5ml centrifuge tubes, respectively. And (3) putting 20 mu L of cell lysate into a 1.5ml centrifuge tube with good light transmittance, adding 100 mu L of LARII reagent, mixing gently, quickly placing in a fluorescence detector, and analyzing the activity of the firefly luciferase. And taking out the centrifuge tube from the fluorescence detector, quickly adding 100 mu L of Stop & GloTM reagent, quickly quenching the activity of firefly luciferase, activating the renilla luciferase reaction of the internal reference plasmid pRL-TK, and placing the renilla luciferase reaction product in the fluorescence detector for analyzing the activity of the renilla luciferase. The ratio of firefly Luciferase Activity to Renilla Luciferase Activity is the Relative Luciferase Activity (RLA), with the RLA results being the average of three independent replicates. The results are shown in FIG. 5, demonstrating that TT mutations enhance ALDOC promoter activity.
10. Effect of SNP (g.20591903C > A) on hepatic tissue ALDOC expression
(1) Liver tissues of 3 cattle with GG and TT genotypes are selected respectively, and RNA is extracted.
(2) A pair of primers ALDOC-mRNA is designed, and the size of the product is 249 bp.
F:5’-agtacgttacagagaaggtcct-3’SEQ ID NO.3
R:5’-cattgaggttgagagatgcct-3’SEQ ID NO.4
And analyzing the expression condition of the ALDOC gene in the liver tissue by using fluorescence quantitative QPCR. The results are shown in FIG. 6 below.
11. SNP sequencing identification method of plateau hypoxia adaptive gene ALDOC
(1) 10 high-altitude cattle and 10 low-altitude cattle are selected to extract blood DNA.
(2) A pair of PCR amplification primers containing the SNP (g.20591903G > T) site was designed.
F:5’-taccttctccatccccctct-3’(SEQ ID NO.5)
R:5’-cgatgaggtggagtgactga-3’(SEQ ID NO.6)
(3) PCR amplification was performed in a 25. mu.L PCR reaction system including 0.5. mu.L (10. mu. mol/L) of the forward primer, 0.5. mu.L (10. mu. mol/L) of the reverse primer, 1. mu.L (50. mu. mol/L) of the DNA template, ddH 2 O10.5. mu.L, 2 XTaq PCR Master Mix 12.5. mu.L, reaction conditions were: pre-denaturation at 94 ℃ for 4min, denaturation at 94 ℃ for 30s, annealing at 60 ℃ for 30s, and extension at 72 ℃ for 30s, wherein the steps are performed for 35 cycles, and finally extension at 72 ℃ for 10min to obtain the target fragment length327 bp. The PCR product was detected by electrophoresis on a 1% agarose gel.
(4) And (3) directly sequencing the amplification product, and analyzing and comparing according to a bovine ALDOC gene sequence published by NCBI (national center for information technology) to obtain individuals with G > T mutation, namely individuals with the genotype of TT, belonging to plateau hypoxia-adaptive cattle (see figure 7).
12. SNP (g.20591903G > T) gene detection method of plateau hypoxia adaptive gene ALDOC
(1) Different cattle varieties are selected, and blood DNA is extracted.
(2) Primers are designed aiming at specific SNP (g.20591903G > T) sites, and restriction enzyme cutting sites BstUI are generated by introducing 1 mutation into an upstream primer, so that BstUI restriction enzymes can cut PCR products with different mutations into fragments with different lengths.
F:5’-GTAATTGTTTACGGTGACGC-3’(SEQ ID NO.7)
R:5’-GGCCTTGTTCTGATTCCTGC-3’(SEQ ID NO.8)
And (3) PCR amplification: PCR amplification was carried out in a 25. mu.L PCR reaction system including 0.5. mu.L (10. mu. mol/L) of the forward primer, 0.5. mu.L (10. mu. mol/L) of the reverse primer, 1. mu.L (50. mu. mol/L) of the DNA template, ddH2O 10.5.5. mu.L, 2 XTaq PCR Master Mix 12.5. mu.L under the following conditions: pre-denaturation at 94 ℃ for 4min, denaturation at 94 ℃ for 30s, annealing at 60 ℃ for 30s, and extension at 72 ℃ for 30s, wherein the steps are performed for 35 cycles, and finally extension at 72 ℃ for 10min, and the target fragment is 104bp in length.
CRS-PCR-RFLP genotyping: carrying out enzyme digestion on the PCR product by using restriction enzyme BstUI, detecting by 2.5% agarose gel electrophoresis after enzyme digestion, and separating two bands of 84bp and 20bp from a wild homozygous individual, wherein only one band, namely 84bp, is displayed in gel because the 20bp fragment is smaller); the heterozygote individual can separate four bands of 104bp, 84bp and 20bp, wherein, the 20bp segment is smaller, so two bands of 104bp and 84bp are displayed in the gel; a band of 104bp can be separated from individuals with homozygous mutant.
Secondly, the embodiment also discloses a kit containing the primer and the enzyme.
The kit also comprises a PCR amplification reaction reagent and an enzyme digestion reaction reagent.
Specifically, the PCR amplification reaction reagent comprises dNTP (25mM each), MgCl 2 (25mM)、PCR Bμffer、ddH 2 O, etc.;
the enzyme digestion reaction comprises ddH 2 O, BstUI enzyme Buffer, BstUI enzyme (10U/. mu.l).
It should be noted that the above examples are only used to illustrate the technical solutions of the present invention and not to limit them. Although the present invention has been described in detail with reference to the examples given, those skilled in the art can modify the technical solution of the present invention as needed or equivalent substitutions without departing from the spirit and scope of the technical solution of the present invention.
SEQUENCE LISTING
<110> Shandong OersX animal husbandry Co., Ltd, Dairy research center of Shandong province academy of agricultural sciences
<120> method for screening bovine plateau hypoxia adaptive gene ALDOC and functional molecular marker and application thereof
<130>
<160> 8
<170> PatentIn version 3.3
<210> 1
<211> 28
<212> DNA
<213> Artificial sequence
<400> 1
ggggtaccta ccttctccat ccccctct 28
<210> 2
<211> 29
<212> DNA
<213> Artificial sequence
<400> 2
ccgctcgagc gatgaggtgg agtgactga 29
<210> 3
<211> 22
<212> DNA
<213> Artificial sequence
<400> 3
agtacgttac agagaaggtc ct 22
<210> 4
<211> 21
<212> DNA
<213> Artificial sequence
<400> 4
cattgaggtt gagagatgcc t 21
<210> 5
<211> 20
<212> DNA
<213> Artificial sequence
<400> 5
<210> 6
<211> 20
<212> DNA
<213> Artificial sequence
<400> 6
<210> 7
<211> 20
<212> DNA
<213> Artificial sequence
<400> 7
<210> 8
<211> 20
<212> DNA
<213> Artificial sequence
<400> 8
ggccttgttc tgattcctgc 20
Claims (8)
1. A method for screening a bovine plateau hypoxia adaptive gene and a functional molecular marker, the method at least comprising:
DNA samples of different altitude cattle varieties are detected and analyzed based on the SNP chip;
performing analysis based on FLK, hapFLK and XPEHH genome selection signal analysis methods, screening single nucleotide polymorphisms SNPs and candidate genes subjected to remarkable positive selection, integrating gene function annotation, and screening plateau hypoxia adaptive candidate genes and SNP molecular markers;
the method for screening the bovine plateau hypoxia adaptive gene and the functional molecular marker comprises the following steps:
s1, collecting and DNA extracting samples of cattle varieties at different altitudes;
s2, SNP chip detection and analysis;
s3, analyzing FLK and hapFLLK genome selection signals;
s4 XPEHH genome selection signal analysis;
s5, identifying a positively selected genetic variation and a differentially selected region;
s6, screening candidate genes based on a screening strategy of the candidate genes;
s7, verifying a positive selection signal of the difference selection area;
s8, identifying potential functional variation in the positive selection gene by using the remeasured data;
wherein, the steps S3 and S4 are not divided in sequence;
the specific method of step S8 includes:
s8.1, acquiring double-end Illumina re-sequencing data, and comparing the acquired data with reference genomes;
s8.2, filtering to obtain a high-quality mapping sequence, and detecting genetic variation from the high-quality mapping sequence;
s8.3, extracting genetic variation from a specific gene region, wherein the extracted regions are an upstream region and a gene region of a candidate gene;
s8.4, respectively calculating the genotype and the allele frequency of each genetic variation;
s8.5, determining potential SNP sites of functional variation through functional annotation;
the bovine plateau hypoxia adaptive gene is ALDOC gene;
the bovine plateau hypoxia adaptive molecular marker also comprises a single nucleotide polymorphic site located in an ALDOC gene, and the single nucleotide polymorphic site comprises g.20588331T > G; g.20590099G > A; g, 20590664T > C; g.20591903G > T.
2. The method as claimed in claim 1, wherein the step S1 includes the following steps:
s1.1, selecting common cattle and rumen cattle varieties at different altitudes, and a mixed variety of the common cattle and the rumen cattle;
s1.2, blood of a cow is collected, and DNA in blood tissues is extracted.
3. The method as claimed in claim 1, wherein the step S2 includes the following steps:
s2.1, analyzing the sample by using the SNP chip, and genotyping;
s2.2, filtering the SNP data, and further analyzing the remaining SNPs meeting the requirements;
s2.3 haplotype was constructed for each chromosome.
4. The method as claimed in claim 1, wherein the step S3 includes the following steps:
s3.1, respectively comparing selection signal values corresponding to SNPs of higher altitude and lower altitude groups based on an FLK and hapFLK whole genome selection signal analysis method, and acquiring difference selection signals of the SNPs between the two groups;
s3.2, respectively operating FLK and hapFLK programs by using genotype data to obtain selection signals of SNPs in the variety, and taking Nelore as a distant population; in the FLK analysis, Nelore is defined as the distant population, and K10 and nti 20;
s3.3, constructing an evolutionary tree of the whole genome and a local genome region for the selected region by using hapFLK results and Python and R scripts;
s3.4 p-value of hapFLK values were calculated by fitting a standard normal distribution of the whole genome in R.
5. The method as claimed in claim 1, wherein the step S4 includes the following steps:
s4.1, estimating XPEHH values between high-altitude and low-altitude cattle varieties, and defining the relation between physical distance and genetic distance by using 1Mb ≈ 1cM in cattle genome;
s4.2 construction of haplotypes based on ReqDel-FAST.
6. The method as claimed in claim 1, wherein step S5 is embodied by the method comprising:
respectively screening SNPs with positive selection signal values of 0.1% in the front based on FLK, hapFLK and XPEHH selection signal analysis results; or the like, or a combination thereof,
the specific method of step S6 includes:
s6.1, defining 50Kbp of the upstream and downstream SNPs obtained by identification as selected areas, and connecting overlapping areas in series to obtain a differential selection area;
s6.2, annotating a reference gene in the difference selection area or overlapped with the area by using the UCSC online website;
s6.3, screening out the FLK, hapFLK and XPEHH selection signals which are compared among high and low altitude groups and selecting the difference selection areas which are 10 th from the top in rank, respectively identifying candidate genes and carrying out functional annotation.
7. The method as claimed in claim 1, wherein the step S7 includes the following steps:
s7.1, respectively constructing a whole genome and a local selection region evolutionary tree based on FLK and hapFLK analysis results, and analyzing the selection condition of the candidate gene in each cattle variety;
s7.2 the results of the comparisons between groups were further validated based on the results of the in-breed FLK and hapFLLK selection signal analysis.
8. Use of the method of any one of claims 1-7 for screening bovine individuals for suitability for high altitude hypoxia survival;
the application mode specifically comprises the following steps:
extracting blood DNA of different cattle individuals;
identifying the genotype of the marker by a CAS-PCR-RFLP method by using the kit for detecting the molecular marker, and identifying individuals with the plateau hypoxia adaptation specific molecular marker;
the molecular marker is a single nucleotide polymorphic site g.20591903G > T;
in this case, the kit comprises at least the following primers and BstUI restriction enzyme;
F: 5‘- GTAATTGTTTACGGTGACGC-3’ （SEQ ID NO.7）；
R: 5‘- GGCCTTGTTCTGATTCCTGC-3’ （SEQ ID NO.8）。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011363669.9A CN112652362B (en) | 2020-11-27 | 2020-11-27 | Method for screening bovine plateau hypoxia adaptive gene ALDOC and functional molecular marker and application thereof |
| AU2021100816A AU2021100816A4 (en) | 2020-11-27 | 2021-02-10 | Method for screening bovine high-altitude hypoxic adaptation gene ALDOC and functional molecular marker and application thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011363669.9A CN112652362B (en) | 2020-11-27 | 2020-11-27 | Method for screening bovine plateau hypoxia adaptive gene ALDOC and functional molecular marker and application thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112652362A CN112652362A (en) | 2021-04-13 |
| CN112652362B true CN112652362B (en) | 2022-08-30 |
Family
ID=75349711
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011363669.9A Active CN112652362B (en) | 2020-11-27 | 2020-11-27 | Method for screening bovine plateau hypoxia adaptive gene ALDOC and functional molecular marker and application thereof |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112652362B (en) |
| AU (1) | AU2021100816A4 (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113774154B (en) * | 2021-10-26 | 2023-09-22 | 山东省农业科学院畜牧兽医研究所 | Method for screening bovine body high mutation related molecular marker and application thereof |
| CN114606325B (en) * | 2022-01-11 | 2024-03-15 | 佛山科学技术学院 | SNP (Single nucleotide polymorphism) marker locus related to bovine thrombopoiesis and application thereof |
| CN117457068B (en) * | 2023-06-30 | 2024-05-24 | 上海睿璟生物科技有限公司 | Multi-genetics-based functional biomarker screening method, system, terminal and medium |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015077717A1 (en) * | 2013-11-25 | 2015-05-28 | The Broad Institute Inc. | Compositions and methods for diagnosing, evaluating and treating cancer by means of the dna methylation status |
| CN105930904A (en) * | 2016-04-12 | 2016-09-07 | 同济大学 | Population evolution-based complex product system adaptive change method |
| CN106164291A (en) * | 2014-03-28 | 2016-11-23 | 中国科学院昆明动物研究所 | For detecting method and the test kit of the genetic marker of litter size in pig |
| CN109295052A (en) * | 2018-11-09 | 2019-02-01 | 云南农业大学 | Wooden dipper chicken is without coda gene and detects the method for wooden dipper chicken anury character, primer, kit |
| CN109360606A (en) * | 2018-11-19 | 2019-02-19 | 广西壮族自治区农业科学院水稻研究所 | A kind of method of low-density SNP genome area Accurate Prediction BSA-seq candidate gene |
| CN109994153A (en) * | 2019-04-09 | 2019-07-09 | 山东省农业科学院奶牛研究中心 | A kind of method and its application for screening ox high altitude hypoxia adaptation molecular labeling |
-
2020
- 2020-11-27 CN CN202011363669.9A patent/CN112652362B/en active Active
-
2021
- 2021-02-10 AU AU2021100816A patent/AU2021100816A4/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015077717A1 (en) * | 2013-11-25 | 2015-05-28 | The Broad Institute Inc. | Compositions and methods for diagnosing, evaluating and treating cancer by means of the dna methylation status |
| CN106164291A (en) * | 2014-03-28 | 2016-11-23 | 中国科学院昆明动物研究所 | For detecting method and the test kit of the genetic marker of litter size in pig |
| CN105930904A (en) * | 2016-04-12 | 2016-09-07 | 同济大学 | Population evolution-based complex product system adaptive change method |
| CN109295052A (en) * | 2018-11-09 | 2019-02-01 | 云南农业大学 | Wooden dipper chicken is without coda gene and detects the method for wooden dipper chicken anury character, primer, kit |
| CN109360606A (en) * | 2018-11-19 | 2019-02-19 | 广西壮族自治区农业科学院水稻研究所 | A kind of method of low-density SNP genome area Accurate Prediction BSA-seq candidate gene |
| CN109994153A (en) * | 2019-04-09 | 2019-07-09 | 山东省农业科学院奶牛研究中心 | A kind of method and its application for screening ox high altitude hypoxia adaptation molecular labeling |
Non-Patent Citations (3)
| Title |
|---|
| Admixture facilitates genetic adaptations to high altitude in Tibet;Choongwon Jeong;《nature communications》;20140210;第1-7页 * |
| Genome‑wide scans identify known and novel regions associated with prolifcacy and reproduction traits in a sub‑Saharan African indigenous sheep (Ovis aries);Asrat Tera Dolebo;《Mammalian Genome》;20191122;第1-14页 * |
| 绵羊TYRP基因快速进化以适应海拔变化;徐方正;《中国优秀博硕士学位论文全文数据库(硕士) 农业科技辑》;20170215;第1-28页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| AU2021100816A4 (en) | 2021-04-22 |
| CN112652362A (en) | 2021-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112652362B (en) | Method for screening bovine plateau hypoxia adaptive gene ALDOC and functional molecular marker and application thereof | |
| WO2020206896A1 (en) | Method for screening molecular marker of cattle adapting to high altitude hypoxia and application thereof | |
| Ghasemi et al. | Genome-wide association study of birth weight in sheep | |
| Tait Jr et al. | µ-Calpain, calpastatin, and growth hormone receptor genetic effects on preweaning performance, carcass quality traits, and residual variance of tenderness in Angus cattle selected to increase minor haplotype and allele frequencies | |
| US20090269741A1 (en) | Method for assessing traits selected from longissimus dorsi peak force, intramuscular fat, retail beef yield and net feed intake in bovine animals | |
| Zhang et al. | Identification of signatures of selection by whole-genome resequencing of a Chinese native pig | |
| CN109694916B (en) | Molecular marker related to sheep feed conversion rate and application thereof | |
| KR101770075B1 (en) | Pig production and traceability systems selected as a single nucleotide polymorphism markers for the introduction of the method | |
| Allan et al. | Genomic mapping of direct and correlated responses to long-term selection for rapid growth rate in mice | |
| Lam et al. | Identification of functional candidate variants and genes for feed efficiency in Holstein and Jersey cattle breeds using RNA-sequencing | |
| Vohra et al. | Genome-wide association studies in Indian Buffalo revealed genomic regions for lactation and fertility | |
| CN107988379A (en) | Genetic marker relevant with tibetan sheep high altitude hypoxia adaptation and its application | |
| Ding et al. | Identify known and novel candidate genes associated with backfat thickness in Duroc pigs by large-scale genome-wide association analysis | |
| Li et al. | The first sheep graph-based pan-genome reveals the spectrum of structural variations and their effects on tail phenotypes | |
| Wang et al. | Genome-wide association study reveals a quantitative trait locus and two candidate genes on Sus scrofa chromosome 5 affecting intramuscular fat content in Suhuai pigs | |
| Zhang et al. | Effects of genetic variations within goat PITX2 gene on growth traits and mRNA expression | |
| Geldermann et al. | Genome-wide mapping of quantitative trait loci for fatness, fat cell characteristics and fat metabolism in three porcine F2 crosses | |
| Amaral et al. | Combining genome-wide association analyses and gene interaction networks to reveal new genes associated with carcass traits, meat quality and fatty acid profiles in pigs | |
| Duchemin et al. | Identification of QTL on chromosome 18 associated with non-coagulating milk in Swedish Red cows | |
| CN112941198A (en) | SNP marker for detecting pig eye muscle area and application thereof | |
| KR100774849B1 (en) | Development of molecular marker of fatty acid-binding proteinFABP gene associated with longissimus muscle area and backfat thickness in Korean cattle | |
| CN103320429A (en) | Method for detecting Qinchuan cattle Wnt7a gene single nucleotide polymorphism, and application thereof | |
| Liu et al. | Identification of LTBP2 gene polymorphisms and their association with thoracolumbar vertebrae number, body size, and carcass traits in Dezhou donkeys | |
| CN111269989B (en) | Pig MID1 gene as mortality-related molecular marker and application thereof | |
| Zegeye et al. | Genome-wide diversity and admixture of five indigenous cattle populations from the Tigray region of northern Ethiopia |
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 |