CN114350845A - Rice storage-resistant gene excavation and molecular marker development - Google Patents
Rice storage-resistant gene excavation and molecular marker development Download PDFInfo
- Publication number
- CN114350845A CN114350845A CN202210085264.6A CN202210085264A CN114350845A CN 114350845 A CN114350845 A CN 114350845A CN 202210085264 A CN202210085264 A CN 202210085264A CN 114350845 A CN114350845 A CN 114350845A
- Authority
- CN
- China
- Prior art keywords
- rice
- storage
- loc
- tolerant
- snp
- 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.)
- Withdrawn
Links
- 235000007164 Oryza sativa Nutrition 0.000 title claims abstract description 58
- 235000009566 rice Nutrition 0.000 title claims abstract description 58
- 108090000623 proteins and genes Proteins 0.000 title claims abstract description 43
- 239000003147 molecular marker Substances 0.000 title claims abstract description 13
- 238000011161 development Methods 0.000 title claims abstract description 11
- 238000009412 basement excavation Methods 0.000 title abstract description 4
- 240000007594 Oryza sativa Species 0.000 title 1
- 241000209094 Oryza Species 0.000 claims abstract description 57
- 238000001514 detection method Methods 0.000 claims abstract description 16
- 238000012216 screening Methods 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 13
- 239000000463 material Substances 0.000 claims abstract description 11
- 230000032683 aging Effects 0.000 claims abstract description 10
- 230000018109 developmental process Effects 0.000 claims abstract description 10
- 240000008467 Oryza sativa Japonica Group Species 0.000 claims abstract description 7
- 230000035784 germination Effects 0.000 claims abstract description 7
- 240000002582 Oryza sativa Indica Group Species 0.000 claims abstract description 4
- 238000011282 treatment Methods 0.000 claims abstract description 4
- 238000006243 chemical reaction Methods 0.000 claims description 21
- 210000000349 chromosome Anatomy 0.000 claims description 21
- 239000000523 sample Substances 0.000 claims description 18
- 238000012163 sequencing technique Methods 0.000 claims description 18
- 230000002596 correlated effect Effects 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 9
- 238000012098 association analyses Methods 0.000 claims description 9
- 238000004925 denaturation Methods 0.000 claims description 9
- 230000036425 denaturation Effects 0.000 claims description 9
- 239000003550 marker Substances 0.000 claims description 9
- 108700028369 Alleles Proteins 0.000 claims description 6
- WSMYVTOQOOLQHP-UHFFFAOYSA-N Malondialdehyde Chemical compound O=CCC=O WSMYVTOQOOLQHP-UHFFFAOYSA-N 0.000 claims description 6
- 238000012408 PCR amplification Methods 0.000 claims description 6
- 238000000137 annealing Methods 0.000 claims description 6
- 238000010219 correlation analysis Methods 0.000 claims description 6
- 235000014113 dietary fatty acids Nutrition 0.000 claims description 6
- 229930195729 fatty acid Natural products 0.000 claims description 6
- 239000000194 fatty acid Substances 0.000 claims description 6
- 150000004665 fatty acids Chemical class 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 5
- 238000005065 mining Methods 0.000 claims description 5
- 238000007400 DNA extraction Methods 0.000 claims description 3
- 241000196324 Embryophyta Species 0.000 claims description 3
- 102000011426 Enoyl-CoA hydratase Human genes 0.000 claims description 3
- 108010023922 Enoyl-CoA hydratase Proteins 0.000 claims description 3
- 102000004195 Isomerases Human genes 0.000 claims description 3
- 108090000769 Isomerases Proteins 0.000 claims description 3
- 101100478969 Oryza sativa subsp. japonica SUS2 gene Proteins 0.000 claims description 3
- 108010059820 Polygalacturonase Proteins 0.000 claims description 3
- 239000013614 RNA sample Substances 0.000 claims description 3
- 108010043934 Sucrose synthase Proteins 0.000 claims description 3
- 238000000246 agarose gel electrophoresis Methods 0.000 claims description 3
- 238000007844 allele-specific PCR Methods 0.000 claims description 3
- 238000009395 breeding Methods 0.000 claims description 3
- 230000001488 breeding effect Effects 0.000 claims description 3
- 230000002860 competitive effect Effects 0.000 claims description 3
- 230000000875 corresponding effect Effects 0.000 claims description 3
- 238000012217 deletion Methods 0.000 claims description 3
- 230000037430 deletion Effects 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 3
- 102000054766 genetic haplotypes Human genes 0.000 claims description 3
- 229940118019 malondialdehyde Drugs 0.000 claims description 3
- 239000000843 powder Substances 0.000 claims description 3
- 238000012257 pre-denaturation Methods 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims description 3
- 102000004169 proteins and genes Human genes 0.000 claims description 3
- 238000003908 quality control method Methods 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 230000035882 stress Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000005382 thermal cycling Methods 0.000 claims description 3
- 238000011222 transcriptome analysis Methods 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- 108700026220 vif Genes Proteins 0.000 claims description 3
- 235000013339 cereals Nutrition 0.000 abstract description 11
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000008303 genetic mechanism Effects 0.000 description 1
- 150000002632 lipids Chemical class 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002773 nucleotide Substances 0.000 description 1
- 125000003729 nucleotide group Chemical group 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
- C12Q1/6895—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for plants, fungi or algae
-
- 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
- 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/6869—Methods for sequencing
-
- 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
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
- G16B25/20—Polymerase chain reaction [PCR]; Primer or probe design; Probe optimisation
-
- 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/10—Sequence alignment; Homology search
-
- 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
-
- 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
- G16B50/00—ICT programming tools or database systems specially adapted for bioinformatics
- G16B50/10—Ontologies; Annotations
-
- 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
- G16B50/00—ICT programming tools or database systems specially adapted for bioinformatics
- G16B50/30—Data warehousing; Computing architectures
-
- 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/13—Plant traits
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Biotechnology (AREA)
- Organic Chemistry (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Genetics & Genomics (AREA)
- General Engineering & Computer Science (AREA)
- Bioethics (AREA)
- Biochemistry (AREA)
- Databases & Information Systems (AREA)
- Microbiology (AREA)
- Immunology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Artificial Intelligence (AREA)
- Mycology (AREA)
- Botany (AREA)
- Software Systems (AREA)
- Public Health (AREA)
- Evolutionary Computation (AREA)
- Epidemiology (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
The invention discloses a rice storage-resistant gene excavation and molecular marker development method, which comprises the following steps: s1, rice germplasm resource collection: selecting 56 parts of main rice variety materials at the middle and lower reaches of Yangtze river, wherein the main rice variety materials comprise 25 indica rice and 31 japonica rice; s2, aging: the sample which is not aged is taken as a reference and is marked as CK; samples aged at a constant temperature and humidity (37 ℃ C., 85% RH) for 5 weeks and 7 weeks were designated as CD5 and CD7, respectively, and the germination percentage and other quantitative properties of the samples were measured for different treatments. By adopting the detection method designed by the invention, KASP molecular markers can be designed according to the identified SNP, so that a detection kit for identifying the storable varieties can be developed conveniently in the future, the storable rice varieties with better quality can be planted by screening the genotype of the storable rice and cultivating the new genotype of the storable quality, and therefore, the accurate regulation and control of grain storage are realized, and the detection method has important significance for reducing the loss of stored grain and ensuring the quality of the stored grain.
Description
Technical Field
The invention relates to the technical field of agriculture, in particular to paddy rice storage-resistant gene excavation and molecular marker development.
Background
At present, the storage technology for slowing down the deterioration of the paddy by controlling the storage environment conditions in the market is mature, the application of the genotype of the paddy in the grain storage safety is still in the initial exploration stage, the molecular marker related to the storage-resistant lipid of the paddy is developed, the genotype of the paddy with storage resistance can be screened, the new genotype with storage-resistant quality can be cultivated, the paddy variety which is suitable for storage and has better quality can be planted, so that the accurate regulation and control of the grain storage can be realized, the accurate regulation and control of the grain storage loss can be realized, the important significance is realized for reducing the grain storage loss and ensuring the grain storage quality, and the influence mechanism research on the storage resistance of the paddy is not clarified because the genetic mechanism of the storage resistance among the paddy varieties is complex, more genes determining the storage resistance are mainly caused by the regulation and control of micro-effect polygene.
Disclosure of Invention
The invention aims to provide the rice storage-resistant gene mining and molecular marker development, has the advantage of being convenient for developing and identifying a detection kit of a storage-resistant variety in the future, and solves the problems provided by the background technology.
In order to achieve the purpose, the invention provides the following technical scheme: the method for mining the storage-resistant gene of the rice and developing the molecular marker comprises the following steps:
s1, rice germplasm resource collection:
selecting 56 parts of main rice variety materials at the middle and lower reaches of Yangtze river, wherein the main rice variety materials comprise 25 indica rice and 31 japonica rice;
s2, aging:
the sample which is not aged is taken as a reference and is marked as CK; accelerating and aging in a constant temperature and humidity box (37 ℃ and 85% RH) for 5 weeks and 7 weeks, respectively marking as CD5 and CD7, and simultaneously measuring the germination percentage, fatty acid value, malonaldehyde content, conductivity and other quantitative properties of the samples under different treatments;
s3, storage-tolerant phenotype identification:
extracting DNA of a rice sample, constructing a library after the DNA sample is detected to be qualified, carrying out machine sequencing after the library is detected to be qualified, adopting Illumina HiSeq/MiSeqe sequencing, carrying out sequencing data quality control analysis on original sequence data obtained by sequencing through FASTQ, comparing the data after sequencing filtration with a reference genome japonica rice variety Nipponbin through Hisat2 software, completing SNP detection development on samples in a group by utilizing GATK software, adopting Plink software, eliminating gene marker SNP with the minimum allele frequency of less than 5 percent and the deletion rate of more than 20 percent, and screening 353196 high-quality SNP markers for next phenotype association analysis;
s4, genome wide association analysis (GWAS) and transcriptome analysis:
(1) selecting a univariate mixed linear model by adopting Gemma software, and screening the obviously associated SNP sites by combining the genetic relationship analysis results among samples;
(2) performing correlation analysis on the phenotype data under the three environments respectively to detect storage-resistant sites which stably appear at different points;
(3) the threshold P-value of the significant associated SNP sites is 10 < -4 >, the total number of SNPs used for association analysis is 254815 (353196 SNPs are developed in the experiment), the functions of SNPs are annotated according to GFF3 annotation files of a Japanese sunny genome, and the Halpoview software is adopted to perform haplotype analysis on target candidate genes in combination with phenotypic data;
(4) significantly associated with the SNP site: through correlation analysis, the rice genotype germination rate is commonly correlated to 2071 SNP sites in the whole genome range, wherein strong signals are arranged on the correlation regions of No. 6 and No. 8 chromosomes, a plurality of SNP sites in the 1.2MB-1.5MB region on the No. 8 chromosome exceed a threshold value, the rice genotype fatty acid character is commonly correlated to 510 SNP sites in the whole genome range, wherein the correlation regions have strong signals on the No. 1, No. 7 and No. 9 chromosomes, the No. 1 chromosome and the No. 9 chromosome signals display strong correlation signals, the rice genotype malondialdehyde character is commonly correlated to 4034 SNP sites in the whole genome range, wherein the strong correlation signals are arranged on the No. 1, 3, 4, 5 and No. 8 chromosomes, the rice genotype conductivity character is commonly correlated to 1510 SNP sites in the whole genome range, and the strong correlation signals appear on the No. 1 chromosome 2;
(5) selecting 4 rice genotypes with extreme storage tolerance according to phenotype data to perform transcriptome sequencing;
(6) storage-tolerant genotype: xiushui 134 and Zhejiang 96; the non-storage-tolerant gene types of Nanjing 9108 and Suxiu 867 adopt a Trizol method to extract total RNA in the rice powder before and after aging, and strictly control the quality of an RNA sample: analyzing the RNA integrity of the sample and the existence of DNA pollution by agarose gel electrophoresis; detecting the purity of RNA by using a NanoPhotometer; then, the RNA concentration and the integrity are accurately quantified by using the Qubit2.0 and Agilent 2100, the sequence is carried out on the machine after the RNA quality is qualified, raw data 300G is obtained in total, a Nipponbare genome is used as a reference genome, Hisat2 software is selected for comparison, R package Deseq2 is used for counting differential expression genes, a rice genome annotation project database is compared, and 11039 differential genes are annotated in total;
(7) combining with transcriptome results, the test analyzes the expression patterns of 14 storage-tolerant candidate genes obtained by GWAS screening, and finds that three of the genes LOC _ Os06g09450.8, LOC _ Os01g45060.1 and LOC _ Os01g47350.1 respectively encode sucrose synthase, polygalacturonase and enoyl-CoA hydratase isomerase family proteins, the expression quantity difference is obvious in storage-tolerant and non-storage-tolerant genotypes, and the storage-tolerant genotypes can be up-regulated to express LOC _ Os06g09450.8 and LOC _ Os01g45060.1 to provide energy for stress reaction;
(8) significantly different genes: the gene is normally expressed in the non-storable genotype, 7 sites (chr1_25581589, chr1_270544634, chr6_4796306, chr6_4800231, chr6_4802754, chr6_4802754 and chr6_4802768) which are related to the storable shape are found in the LOC _ Os01g45060, LOC _ Os01g47350 and LOC _ Os06g09450 genes, and one related site (chr1_25579324) is dug in the upstream 1000bp range of the LOC _ Os01g 45060;
s5, screening of storage-resistant candidate genes:
the method comprises the steps of finding 8 strongly-associated SNP sites near 3 genes associated with the storage-tolerant traits, developing 8 corresponding pairs of fluorescent molecular markers based on the 8 sites, and further verifying and typing SNPs associated with the DNA samples of rice genomes of different varieties by combining a competitive Allele Specific PCR (KASP) technology;
s6, designing a universal fluorescent primer:
KASP requires preparation:
(1) the Primer Mix consists of two allele forward primers with different terminal bases and a reverse Primer, and the 5' ends of the two forward primers are respectively connected with detection Primer sequences with different sequences;
(2) two detection primers Master Mix with different fluorescence; 3) a DNA template;
the KASP reaction was carried out on a LightCycler480(LC480), and the PCR amplification reaction system and reaction program were determined by the method of Weiyu et al [25 ]:
adding 50-100 ng. mu.L-1 of the diluted DNA template into a 96-well plate respectively, and preparing the following reaction system for PCR amplification reaction: 5 μ L of template DNA, 2 XKASP Master mix5 μ L of KASP marker primer, 0.14 μ L of KASP marker primer, and the total reaction volume of each well is 10.14 μ L;
the PCR reaction program on LightCycler480(LC480) was:
(1) KASP thermal cycling: pre-denaturation at 94 ℃ for 15 min; denaturation at 94 ℃ for 20 s; annealing at 61-55 deg.C for 60s, 10 cycles (0.6 deg.C reduction per cycle); denaturation at 94 ℃ for 20 s; annealing and extending at 55 ℃ for 60s for 26 cycles; reading at 37 ℃ for 1 min;
(2) and (3) KASP recovery: denaturation at 94 ℃ for 20 s; extending for 60s at 57 ℃, and repeating the steps (1) - (2) twice (totally 3 times); reading at 37 ℃ for 1 min;
(3) after the reaction is finished, judging the sample typing condition according to the detected 2 fluorescent signals, wherein different genotypes are obtained by different fluorescent signals;
s7, KASP molecular marker typing:
the experimental result proves that the 8 SNP markers exist in the storage-resistant material and are available molecular markers of the rice storage-resistant trait QTL, and the molecular markers can provide basis for obtaining individual plants of storage-resistant varieties later and lay a foundation for the subsequent breeding of the storage-resistant varieties.
Preferably, in step S3, the DNA of the rice sample is extracted by using a DNA extraction kit.
Compared with the prior art, the invention has the following beneficial effects:
by adopting the detection method designed by the invention, KASP molecular markers can be designed according to the identified SNP, so that a detection kit for identifying the storable varieties can be developed conveniently in the future, the storable rice varieties with better quality can be planted by screening the genotype of the storable rice and cultivating the new genotype of the storable quality, and therefore, the accurate regulation and control of grain storage are realized, and the detection method has important significance for reducing the loss of stored grain and ensuring the quality of the stored grain.
Drawings
FIG. 1 is a schematic flow chart of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The components of the present invention are all standard components or components known to those skilled in the art, and the structure and principle thereof can be known to those skilled in the art through technical manuals or through routine experimentation.
Referring to fig. 1, the mining of rice storage-resistant genes and the development of molecular markers includes the following steps:
s1, rice germplasm resource collection:
selecting 56 parts of main rice variety materials at the middle and lower reaches of Yangtze river, wherein the main rice variety materials comprise 25 indica rice and 31 japonica rice;
s2, aging:
the sample which is not aged is taken as a reference and is marked as CK; accelerating and aging in a constant temperature and humidity box (37 ℃ and 85% RH) for 5 weeks and 7 weeks, respectively marking as CD5 and CD7, and simultaneously measuring the germination percentage, fatty acid value, malonaldehyde content, conductivity and other quantitative properties of the samples under different treatments;
s3, storage-tolerant phenotype identification:
extracting DNA of a rice sample, constructing a library after the DNA sample is detected to be qualified, performing machine sequencing after the library is detected to be qualified, performing Illumina HiSeq/MiSeqe sequencing, performing sequencing data quality control analysis on original sequence data obtained by sequencing through FASTQ, comparing the data after sequencing filtration with a reference genome japonica rice variety Nipponbin through Hisat2 software, completing SNP detection development on samples in a group by utilizing GATK software, removing gene marker SNP (single nucleotide polymorphism) with the minimum allele frequency of less than 5% and the deletion rate of more than 20% by utilizing Plink software, and screening 353196 high-quality SNP markers for next phenotype association analysis in total, wherein a DNA extraction kit is adopted to extract the DNA of the rice sample in the step S3;
s4, genome wide association analysis (GWAS) and transcriptome analysis:
(1) selecting a univariate mixed linear model by adopting Gemma software, and screening the obviously associated SNP sites by combining the genetic relationship analysis results among samples;
(2) performing correlation analysis on the phenotype data under the three environments respectively to detect storage-resistant sites which stably appear at different points;
(3) the threshold P-value of the significant associated SNP sites is 10 < -4 >, the total number of SNPs used for association analysis is 254815 (353196 SNPs are developed in the experiment), the functions of SNPs are annotated according to GFF3 annotation files of a Japanese sunny genome, and the Halpoview software is adopted to perform haplotype analysis on target candidate genes in combination with phenotypic data;
(4) significantly associated with the SNP site: through correlation analysis, the rice genotype germination rate is commonly correlated to 2071 SNP sites in the whole genome range, wherein strong signals are arranged on the correlation regions of No. 6 and No. 8 chromosomes, a plurality of SNP sites in the 1.2MB-1.5MB region on the No. 8 chromosome exceed a threshold value, the rice genotype fatty acid character is commonly correlated to 510 SNP sites in the whole genome range, wherein the correlation regions have strong signals on the No. 1, No. 7 and No. 9 chromosomes, the No. 1 chromosome and the No. 9 chromosome signals display strong correlation signals, the rice genotype malondialdehyde character is commonly correlated to 4034 SNP sites in the whole genome range, wherein the strong correlation signals are arranged on the No. 1, 3, 4, 5 and No. 8 chromosomes, the rice genotype conductivity character is commonly correlated to 1510 SNP sites in the whole genome range, and the strong correlation signals appear on the No. 1 chromosome 2;
(5) selecting 4 rice genotypes with extreme storage tolerance according to phenotype data to perform transcriptome sequencing;
(6) storage-tolerant genotype: xiushui 134 and Zhejiang 96; the non-storage-tolerant gene types of Nanjing 9108 and Suxiu 867 adopt a Trizol method to extract total RNA in the rice powder before and after aging, and strictly control the quality of an RNA sample: analyzing the RNA integrity of the sample and the existence of DNA pollution by agarose gel electrophoresis; detecting the purity of RNA by using a NanoPhotometer; then, the RNA concentration and the integrity are accurately quantified by using the Qubit2.0 and Agilent 2100, the sequence is carried out on the machine after the RNA quality is qualified, raw data 300G is obtained in total, a Nipponbare genome is used as a reference genome, Hisat2 software is selected for comparison, R package Deseq2 is used for counting differential expression genes, a rice genome annotation project database is compared, and 11039 differential genes are annotated in total;
(7) combining with transcriptome results, the test analyzes the expression patterns of 14 storage-tolerant candidate genes obtained by GWAS screening, and finds that three of the genes LOC _ Os06g09450.8, LOC _ Os01g45060.1 and LOC _ Os01g47350.1 respectively encode sucrose synthase, polygalacturonase and enoyl-CoA hydratase isomerase family proteins, the expression quantity difference is obvious in storage-tolerant and non-storage-tolerant genotypes, and the storage-tolerant genotypes can be up-regulated to express LOC _ Os06g09450.8 and LOC _ Os01g45060.1 to provide energy for stress reaction;
(8) significantly different genes: the gene is normally expressed in the non-storable genotype, 7 sites (chr1_25581589, chr1_270544634, chr6_4796306, chr6_4800231, chr6_4802754, chr6_4802754 and chr6_4802768) which are related to the storable shape are found in the LOC _ Os01g45060, LOC _ Os01g47350 and LOC _ Os06g09450 genes, and one related site (chr1_25579324) is dug in the upstream 1000bp range of the LOC _ Os01g 45060;
s5, screening of storage-resistant candidate genes:
the method comprises the steps of finding 8 strongly-associated SNP sites near 3 genes associated with the storage-tolerant traits, developing 8 corresponding pairs of fluorescent molecular markers based on the 8 sites, and further verifying and typing SNPs associated with the DNA samples of rice genomes of different varieties by combining a competitive Allele Specific PCR (KASP) technology;
s6, designing a universal fluorescent primer:
KASP requires preparation:
(1) the Primer Mix consists of two allele forward primers with different terminal bases and a reverse Primer, and the 5' ends of the two forward primers are respectively connected with detection Primer sequences with different sequences;
(2) two detection primers Master Mix with different fluorescence; 3) a DNA template;
the KASP reaction was carried out on a LightCycler480(LC480), and the PCR amplification reaction system and reaction program were determined by the method of Weiyu et al [25 ]:
adding 50-100 ng. mu.L-1 of the diluted DNA template into a 96-well plate respectively, and preparing the following reaction system for PCR amplification reaction: 5 μ L of template DNA, 2 XKASP Master mix5 μ L of KASP marker primer, 0.14 μ L of KASP marker primer, and the total reaction volume of each well is 10.14 μ L;
the PCR reaction program on LightCycler480(LC480) was:
(1) KASP thermal cycling: pre-denaturation at 94 ℃ for 15 min; denaturation at 94 ℃ for 20 s; annealing at 61-55 deg.C for 60s, 10 cycles (0.6 deg.C reduction per cycle); denaturation at 94 ℃ for 20 s; annealing and extending at 55 ℃ for 60s for 26 cycles; reading at 37 ℃ for 1 min;
(2) and (3) KASP recovery: denaturation at 94 ℃ for 20 s; extending for 60s at 57 ℃, and repeating the steps (1) - (2) twice (totally 3 times); reading at 37 ℃ for 1 min;
(3) after the reaction is finished, judging the sample typing condition according to the detected 2 fluorescent signals, wherein different genotypes are obtained by different fluorescent signals;
s7, KASP molecular marker typing:
the experimental result proves that the 8 SNP markers exist in the storage-resistant material and are available molecular markers of the rice storage-resistant property QTL, the molecular markers can be used for obtaining single plants of storage-resistant varieties later, a basis is provided for the subsequent breeding of the storage-resistant varieties, and a foundation is laid.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.
Claims (2)
1. The method for mining the rice storage-resistant gene and developing the molecular marker is characterized in that: the method comprises the following steps:
s1, rice germplasm resource collection:
selecting 56 parts of main rice variety materials at the middle and lower reaches of Yangtze river, wherein the main rice variety materials comprise 25 indica rice and 31 japonica rice;
s2, aging:
the sample which is not aged is taken as a reference and is marked as CK; accelerating and aging in a constant temperature and humidity box (37 ℃ and 85% RH) for 5 weeks and 7 weeks, respectively marking as CD5 and CD7, and simultaneously measuring the germination percentage, fatty acid value, malonaldehyde content, conductivity and other quantitative properties of the samples under different treatments;
s3, storage-tolerant phenotype identification:
extracting DNA of a rice sample, constructing a library after the DNA sample is detected to be qualified, carrying out machine sequencing after the library is detected to be qualified, adopting Illumina HiSeq/MiSeqe sequencing, carrying out sequencing data quality control analysis on original sequence data obtained by sequencing through FASTQ, comparing the data after sequencing filtration with a reference genome japonica rice variety Nipponbin through Hisat2 software, completing SNP detection development on samples in a group by utilizing GATK software, adopting Plink software, eliminating gene marker SNP with the minimum allele frequency of less than 5 percent and the deletion rate of more than 20 percent, and screening 353196 high-quality SNP markers for next phenotype association analysis;
s4, genome wide association analysis (GWAS) and transcriptome analysis:
(1) selecting a univariate mixed linear model by adopting Gemma software, and screening the obviously associated SNP sites by combining the genetic relationship analysis results among samples;
(2) performing correlation analysis on the phenotype data under the three environments respectively to detect storage-resistant sites which stably appear at different points;
(3) the threshold P-value of the significant associated SNP sites is 10 < -4 >, the total number of SNPs used for association analysis is 254815 (353196 SNPs are developed in the experiment), the functions of SNPs are annotated according to GFF3 annotation files of a Japanese sunny genome, and the Halpoview software is adopted to perform haplotype analysis on target candidate genes in combination with phenotypic data;
(4) significantly associated with the SNP site: through correlation analysis, the rice genotype germination rate is commonly correlated to 2071 SNP sites in the whole genome range, wherein strong signals are arranged on the correlation regions of No. 6 and No. 8 chromosomes, a plurality of SNP sites in the 1.2MB-1.5MB region on the No. 8 chromosome exceed a threshold value, the rice genotype fatty acid character is commonly correlated to 510 SNP sites in the whole genome range, wherein the correlation regions have strong signals on the No. 1, No. 7 and No. 9 chromosomes, the No. 1 chromosome and the No. 9 chromosome signals display strong correlation signals, the rice genotype malondialdehyde character is commonly correlated to 4034 SNP sites in the whole genome range, wherein the strong correlation signals are arranged on the No. 1, 3, 4, 5 and No. 8 chromosomes, the rice genotype conductivity character is commonly correlated to 1510 SNP sites in the whole genome range, and the strong correlation signals appear on the No. 1 chromosome 2;
(5) selecting 4 rice genotypes with extreme storage tolerance according to phenotype data to perform transcriptome sequencing;
(6) storage-tolerant genotype: xiushui 134 and Zhejiang 96; the non-storage-tolerant gene types of Nanjing 9108 and Suxiu 867 adopt a Trizol method to extract total RNA in the rice powder before and after aging, and strictly control the quality of an RNA sample: analyzing the RNA integrity of the sample and the existence of DNA pollution by agarose gel electrophoresis; detecting the purity of RNA by using a NanoPhotometer; then, the RNA concentration and the integrity are accurately quantified by using the Qubit2.0 and Agilent 2100, the sequence is carried out on the machine after the RNA quality is qualified, raw data 300G is obtained in total, a Nipponbare genome is used as a reference genome, Hisat2 software is selected for comparison, R package Deseq2 is used for counting differential expression genes, a rice genome annotation project database is compared, and 11039 differential genes are annotated in total;
(7) combining with transcriptome results, the test analyzes the expression patterns of 14 storage-tolerant candidate genes obtained by GWAS screening, and finds that three of the genes LOC _ Os06g09450.8, LOC _ Os01g45060.1 and LOC _ Os01g47350.1 respectively encode sucrose synthase, polygalacturonase and enoyl-CoA hydratase isomerase family proteins, the expression quantity difference is obvious in storage-tolerant and non-storage-tolerant genotypes, and the storage-tolerant genotypes can be up-regulated to express LOC _ Os06g09450.8 and LOC _ Os01g45060.1 to provide energy for stress reaction;
(8) significantly different genes: the gene is normally expressed in the non-storable genotype, 7 sites (chr1_25581596, chr1_27054464, chr6_4796306, chr6_4800231, chr6_4802745, chr6_4802754 and chr6_4802768) which are related to the storable shape are found in the LOC _ Os01g45060, LOC _ Os01g47350 and LOC _ Os06g09450 genes, and one related site (chr1_25579324) is dug in the upstream 1000bp range of the LOC _ Os01g 45060;
s5, screening of storage-resistant candidate genes:
the method comprises the steps of finding 8 strongly-associated SNP sites near 3 genes associated with the storage-tolerant traits, developing 8 corresponding pairs of fluorescent molecular markers based on the 8 sites, and further verifying and typing SNPs associated with the DNA samples of rice genomes of different varieties by combining a competitive Allele Specific PCR (KASP) technology;
s6, designing a universal fluorescent primer:
KASP requires preparation:
(1) the Primer Mix consists of two allele forward primers with different terminal bases and a reverse Primer, and the 5' ends of the two forward primers are respectively connected with detection Primer sequences with different sequences;
(2) two detection primers Master Mix with different fluorescence; 3) a DNA template;
the KASP reaction was carried out on a LightCycler480(LC480), and the PCR amplification reaction system and reaction program were determined by the method of Weiyu et al [25 ]:
adding 50-100 ng. mu.L-1 of the diluted DNA template into a 96-well plate respectively, and preparing the following reaction system for PCR amplification reaction: 5 μ L of template DNA, 2 XKASP Master mix5 μ L of KASP marker primer, 0.14 μ L of KASP marker primer, and the total reaction volume of each well is 10.14 μ L;
the PCR reaction program on LightCycler480(LC480) was:
(1) KASP thermal cycling: pre-denaturation at 94 ℃ for 15 min; denaturation at 94 ℃ for 20 s; annealing at 61-55 deg.C for 60s, 10 cycles (0.6 deg.C reduction per cycle); denaturation at 94 ℃ for 20 s; annealing and extending at 55 ℃ for 60s for 26 cycles; reading at 37 ℃ for 1 min;
(2) and (3) KASP recovery: denaturation at 94 ℃ for 20 s; extending for 60s at 57 ℃, and repeating the steps (1) - (2) twice (totally 3 times); reading at 37 ℃ for 1 min;
(3) after the reaction is finished, judging the sample typing condition according to the detected 2 fluorescent signals, wherein different genotypes are obtained by different fluorescent signals;
(4) development of molecular marker sequences:
s7, KASP molecular marker typing:
the experimental result proves that the 8 SNP markers exist in the storage-resistant material and are available molecular markers of the rice storage-resistant trait QTL, and the molecular markers can provide basis for obtaining individual plants of storage-resistant varieties later and lay a foundation for the subsequent breeding of the storage-resistant varieties.
2. The rice storage-tolerant gene mining and molecular marker development according to claim 1, wherein: in the step S3, DNA of the rice sample is extracted by using a DNA extraction kit.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210085264.6A CN114350845A (en) | 2022-01-25 | 2022-01-25 | Rice storage-resistant gene excavation and molecular marker development |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210085264.6A CN114350845A (en) | 2022-01-25 | 2022-01-25 | Rice storage-resistant gene excavation and molecular marker development |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114350845A true CN114350845A (en) | 2022-04-15 |
Family
ID=81093693
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210085264.6A Withdrawn CN114350845A (en) | 2022-01-25 | 2022-01-25 | Rice storage-resistant gene excavation and molecular marker development |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114350845A (en) |
-
2022
- 2022-01-25 CN CN202210085264.6A patent/CN114350845A/en not_active Withdrawn
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115029451B (en) | Sheep liquid phase chip and application thereof | |
CN108103235B (en) | SNP molecular marker and primer for identifying cold resistance of apple rootstock and application of SNP molecular marker and primer | |
CN111088382B (en) | Corn whole genome SNP chip and application thereof | |
CN114657238B (en) | Medlar 40K liquid phase chip and application | |
CN109554493A (en) | The SNP marker and its detection method of a kind of plum blossom weeping branch character close linkage and application | |
CN110846429A (en) | Corn whole genome InDel chip and application thereof | |
CN114350841A (en) | Polymorphic molecular marker based on whole genome sequencing, preparation method and application | |
CN112289384A (en) | Construction method and application of whole citrus genome KASP marker library | |
CN113564266B (en) | SNP typing genetic marker combination, detection kit and application | |
CN108642201B (en) | SNP (Single nucleotide polymorphism) marker related to millet plant height character as well as detection primer and application thereof | |
CN107815502B (en) | dCAPS marker development and application for identifying waterlogging tolerance of chrysanthemum | |
CN117683927A (en) | Functional KASP molecular marker of rice blast resistance gene and application thereof | |
CN109762925B (en) | KASP molecular marker of rice blast-resistant broad-spectrum gene pi9, detection method and application thereof | |
CN107447022B (en) | SNP molecular marker for predicting corn heterosis and application thereof | |
CN114350845A (en) | Rice storage-resistant gene excavation and molecular marker development | |
CN115141893A (en) | Molecular marker group containing 7 molecular markers and used for predicting dry matter content of kiwi fruit, application and kit thereof | |
CN108715901B (en) | SNP marker related to millet plant height character and detection primer and application thereof | |
CN108642199B (en) | SNP (Single nucleotide polymorphism) marker related to growth of millet flag leaves as well as detection primer and application thereof | |
CN117821636B (en) | SNP molecular marker of rice blast resistance gene Pib, primer set and application | |
CN111705157A (en) | PCR marker and primer for screening I-type S haplotypes of cabbage self-incompatible line | |
CN118064638B (en) | SNP molecular marker locus related to drought tolerance of corn and application thereof | |
CN116814841B (en) | Primer group for identifying rice black brown glume gene HK4, and method and application thereof | |
CN110616256B (en) | Multi-site cynoglossus semilaevis true and false male fish screening system based on SNaPshot technology and application | |
CN108728566B (en) | SNP (Single nucleotide polymorphism) marker related to thousand grain weight traits of millet as well as detection primer and application thereof | |
CN117385093A (en) | SNP molecular marker closely linked with brassica napus anti-cracking gene BnMAN7A07 and application thereof |
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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20220415 |
|
WW01 | Invention patent application withdrawn after publication |