Rubber tree chloroplast SNP (Single Nucleotide Polymorphism) marker-based rubber tree B-type cultivated variety identification method

Abstract

The invention belongs to the technical field of molecular markers, and particularly relates to a rubber tree chloroplast SNP marker-based rubber tree B-type cultivated variety identification method, 19 SNP loci are involved in total, an SNP marker combination provided by the invention can be used for identifying rubber tree varieties with higher yield, and by utilizing the SNP marker combination provided by the invention, the yield of the rubber tree varieties is increased. And variety improvement and new variety breeding can be more accurately carried out, so that good varieties with higher yield are selected, and the planting benefits of the rubber trees are improved.

Description

Technical Field

[0001] The invention belongs to the technical field of molecular markers, and in particular relates to a method for identifying a B-type cultivar of rubber tree based on a rubber tree chloroplast SNP marker. Background Art

[0002] Natural rubber mainly comes from the Brazilian rubber tree. It is an important industrial raw material and is widely used in transportation, medical and health fields. Its excellent properties in elasticity, wear resistance, impact resistance, etc. have not been replaced by synthetic rubber so far. The Brazilian rubber tree is native to the Amazon Basin in Brazil. After many introductions, it is now widely cultivated in tropical Asia. It is cultivated in Taiwan Province, southern Fujian, Guangdong, Guangxi, Hainan and southern Yunnan, with Hainan and Yunnan being the most cultivated. There are two famous introductions of rubber trees in history. In 1876, the British Wickham collected rubber tree germplasm in Brazil, and then brought the seeds and seedlings to Kew Garden in the UK for breeding. Later, they were planted in Sri Lanka and then introduced to Malaysia, Indonesia, Thailand and other Southeast Asian countries for planting. The second time was in 1981 when the International Rubber Research and Development Board (IRRDB) organized the collection and introduction of rubber tree germplasm. In addition, it also includes other rubber tree germplasm currently available in major rubber-producing countries. After years of introduction, cultivation, hybrid breeding and screening, my country has bred a number of cold-resistant, wind-resistant and high-yield varieties, including Reyan 879, Reyan 73397 and other high-yield varieties, Reyan 917 and other wind-resistant high-yield varieties, Yunyan 774 and other cold-resistant high-yield varieties. Among them, Reyan 879 is a super high-yield variety, with an per-acre yield increased by more than 50% over the control, and is currently the variety with the highest yield in the world.

[0003] With the continuous advancement of sequencing technology, more and more rubber tree varieties, germplasms and rubber tree species genomes have been analyzed, laying a solid foundation for further genome selection and molecular-assisted breeding, and accelerating the rubber breeding process. Chloroplasts are the site of photosynthesis in green plants. They have an independent genetic system, their genetic structure is relatively stable, and gene recombination rarely occurs. Therefore, chloroplast genome sequences are widely used in plant phylogeny and species identification. Plant genome resequencing is mostly based on leaves, and mesophyll cells contain a large number of chloroplasts. In addition, the size of the chloroplast genome is also very small, only more than 100 kb. Therefore, compared with the nuclear genome, only a small amount of leaf resequencing data can be used for chloroplast genome assembly and SNP analysis. Rubber tree cultivars are all derived from the screening of wild germplasms in the Amazon Basin of Brazil. They generally have characteristics such as high yield and stress resistance. By comparing the chloroplast genomes of cultivated varieties and wild germplasms, screening effective SNP marker sites will accelerate the breeding process of high-yield and stress-resistant rubber trees. Summary of the invention

[0004] The purpose of the present invention is to provide a SNP marker combination that can be used for high-yield varieties of rubber trees, enrich the molecular markers available for breeding excellent varieties of rubber trees, and provide new ideas for high-yield and stress-resistant breeding of rubber trees.

[0005] The invention provides a SNP marker combination for identifying high-yield varieties of rubber trees. The SNP marker combination includes s8165-T / G, s9952-G / A, s23507-G / A, s26511-A / G, s28467-C / A, s29302-T / G, s29303-C / A, s29894-G / T, s43759-T / C, s46966-A / C, s69276-G / T, s69443-C / A, s74109-G / A, s79616-G / T, s86288-C / T, s119902-C / A, s122256-C / T, s129310-G / A and s133745-T / C; the s8165-T / G is located at SEQ ID The sequence shown in No.1 has a T / G natural variation at the 65th position from the 5' end, specifically: 5'-AGTCAAAATACAATTATTTTCAAAAAAAAAATGCTTGTTATGCTTAATATTTTTAGTTTAATTTTTATCTGTTTTAATTCTGCCATTTTTTCAAGCAATTTTTTCTTTACAAAATTGCCCGAAGCCTACGCCTTTTTGAA-3'; the corresponding wild-type sequence is shown in SEQ ID No.2, specifically: 5'-AGTCAAAATACAATTATTTTCAAAAAAAAAATGCTTGTTATGCTTAATATTTTTAGTTTAATTTGTATCTGTTTTAATTCTGCCATTTTTTCAAGCAATTTTTTCTTTACAAAATTGCCCGAAGCCTACGCCTTTTTGAA-3'; the s9952-G / A is located at SEQ ID The sequence shown in No.3 has a G / A natural variation at position 72 from the 5' end, specifically: 5'-TAATTATAATGGATAATGACAATTAGAATGAATCTTTCATGAATATAAGAATCAAAAAATTCTATGGAATCGTGAAAGACAGAAAGATTCTTCGAATTTAGTCATAACATTTCGTTATATTGACAATTTCAAAAACTGTTCATACTATGAGCATAGTATG-3'; the corresponding wild-type sequence is shown in SEQ ID No.4, specifically: 5'-TAATTATAATGGATAATGACAATTAGAATGAATCTTTCATGAATATAAGAATCAAAAAATTCTATGGAATCATGAAAGACAGAAAGATTCTTCGAATTTAGTCATAACATTTCGTTATATTGACAATTTCAAAAACTGTTCATACTATGAGCATAGTATG-3'; the s23507-G / A is located at the 67th position from the 5' end of the sequence shown in SEQ ID No.5, and there is a natural G / A variation at this site, specifically: 5'-GTAGATACACAAATCGCTTTAAATACAAGAAGTCGGGTGGGCGGATTGGTCCGAGTGGAGAGAAAAGAAAAAAAAATGGAACTTAAAATCTTTTCTGGAGATATCCATTTTCCGGGAGAGACAGATAAAATATCCCGACA-3'; the corresponding wild-type sequence is shown in SEQ ID No.6, specifically: 5'-GTAGATACACAAATCGCTTTAAATACAAGAAGTCGGGTGGGCGGATTGGTCCGAGTGGAGAGAAAAAAAAAAAAAATGGAACTTAAAATCTTTTCTGGAGATATCCATTTTCCGGGAGAGACAGATAAAATATCCCGACA-3'; the s26511-A / G is located at the 81st position from the 5' end of the sequence shown in SEQ ID No.7, and there is a natural A / G variation at this site, specifically: 5'-TTCGATTCCAACGAATGATGACGCTATAGCTTCAATCCGATTAATTCTTAATAAATTAGTATTTGCAATTTGTGAGGGTCATTCTAGCTATATACGAAATCCCTGATTAGGAATAAGTAAGATAGATTAAATAAGATAACTAAATTATTTTGTGAACTCC-3'; the corresponding wild-type sequence is shown in SEQ ID No.8, specifically: 5'-TTCGATTCCAACGAATGATGACGCTATAGCTTCAATCCGATTAATTCTTAATAAATTAGTATTTGCAATTTGTGAGGGTCGTTCTAGCTATATACGAAATCCCTGATTAGGAATAAGTAAGATAGATTAAATAAGATAACTAAATTATTTTGTGAACTCC-3'; the s28467-C / A is located in SEQ ID No.The sequence shown in 9 has a C / A natural variation at position 87 from the 5' end, specifically: 5'-TCCTTTCTTTTCGTTTGTAACTGTGAATTGAATATTGAATAAACAGATTAAATCAAGAAAAAGAACGAAGAGTTCGAATTCAAAAACTGATTCGTAACAAAGAAAGAAATGGAAAAACAAAAAAGTTGTATGAATTCACAAAAACTTTGTGAATAGGAAC-3'; the corresponding wild-type sequence is shown in SEQ ID No.10, specifically: 5'-TCCTTTCTTTTCGTTTGTAACTGTGAATTGAATATTGAATAAACAGATTAAATCAAGAAAAAGAACGAAGAGTTCGAATTCAAAAAATGATTCGTAACAAAGAAAGAAATGGAAAAACAAAAAAGTTGTATGAATTCACAAAAACTTTGTGAATAGGAAC-3'; the s29302-T / G is located at SEQ ID The sequence shown in No.11 is at the 142nd position from the 5' end, and there is a T / G natural variation at this site. The s29303-C / A is located at the 143rd position from the 5' end of the sequence shown in SEQ ID No.11, and there is a C / A natural variation at this site, specifically: 5'-CCAATCGAAACCTTTTTTTTAAGGAAATGTTGCAACAGGCATTTTTGTTTTGCAATTAATACGAGACCTGTTACCGAATTCCCGAATTCTAAACTGAATTTTAATAAGAATAATATTTATTAGTATTAGAGATTTCAAATTTCAAATTGAAAAAAAAAATGCATTTATAATTCGAATATAAATCGAATATTTTTTTACTCTATTTTTAAATAGAATAAAAAAAAAAAACAAGAAAAAAACAAGAACAATAATAAACAATAACATAT-3'; the corresponding wild-type sequence is shown in SEQ ID No.12, specifically: 5'-CCAATCGAAACCTTTTTTTTAAGGAAATGTTGCAACAGGCATTTTTGTTTTGCAATTAATACGAGACCTGTTACCGAATTCCCGAATTCTAAACTGAATTTTAATAAGAATAATATTTATTAGTATTAGAGATTTCAAA TTGAAAATTGAAAAAAAAAATGCATTTATAATTCGAATATAAATCGAATATTTTTTTTACTCTATTTTTAAATAGAATAAAAAAAAAAAAAAACAAGAAAAAAACAAGAACAATAATAAACAATAACATAT-3'; the s29894-G / T is located in SEQ. ID The sequence shown in No. 13 has a G / T natural variation at position 134 from the 5' end, specifically: 5'-ATTCGAAATTCAGAAAAACTACGCGAGGGGGCTATTGAACAGCTGGAAAAAGCCCGGGCCCGCTTACGGAAAGTGGAAATAGAAGCAGATCAGTTTCGAACGAATGGATATTCTGAGATAGAACGAGAAAAATGGAATTTGATTAATTCAACTTATAAGACTTTGGAACAATTAGAAAATTACAAAAATGAAACCATTCATTTTGAACAACAACGAACGATTAATCAAGTCCGACAACGGGTTTTCCAAC-3'; the corresponding wild-type sequence is shown in SEQ ID As shown in No. 14, specifically: 5'-ATTCGAAATTCAGAAAAACTACGCGAGGGGCTATTGAACAGCTGGAAAAAGCCCGGGCCCGCTTACGGAAAGTGGAAATAGAAGCAGATCAGTTTCGAACGAATGGATATTCTGAGATAGAACGAGA AAAATTGAATTTGATTAATTCAACTTATAAGACTTTGGAACAATTAGAAAATTACAAAAATGAAACCATTCATTTTGAACAACAACGAACGATTAATCAAGTCCGACAACGGGTTTTTCCAAC-3'; the s43759-T / C is located in SEQ. ID No.The sequence shown in 15 has a T / C natural variation at position 139 from the 5' end, specifically: 5'-AAGATTTGCTTTATCCGGTATCAAACGAGAGCTGCGGGCAAATAGAACTCCTTTCAGAAGTATCAATACCGTCACATGAATCGTGAATGCATGAATGTGATGTACCAAAAAATCCGCGGTTCCTAATGGAATCGGTAATAAAGCAACCTTGCCACCAACTGCCACTAAATCACCACCCCCCCAAGTTAAACTGGTGCTTGCTGTTGCACCAGGAGCCGTTGCACCAGGTGCTAAAGCATGGGTGTTTTGTATCCATTGAGCAAAGACAGG-3'; the corresponding wild-type sequence is shown in SEQ ID As shown in No. 16, specifically: 5'-AAGATTTGCTTTATCCGGTATCAAACGAGAGCTGCGGGCAAATAGAACTCCTTTCAGAAGTATCAATACCGTCACATGAATCGTGAATGCATGAATGTGATGTACCAAAAAATCCGCGGTTCCTAATGGAATCGGTAA CAAAGCAACCTTGCCACCAACTGCCACTAAATCACCACCCCCCCAAGTTAAACTGGTGCTTGCTGTTGCACCAGGAGCCGTTGCACCAGGTGCTAAAGCATGGGTGTTTTGTATCCATTGAGCAAAGACAGG-3'; the s46966-A / C is located in SEQ ID The sequence shown in No.17 has a natural A / C variation at position 136 from the 5' end, specifically: 5'-AAGCCTATTTCTAGTAAATTTTTTTCTTTTTTTTTTTTTTTTTCTTTTTTCTTTCTATAGTGGAGATAGTCGCACGTAATGACAGATCACGGCCATATTATTTAAAGCTTGTGGTAAGAAGGGGTTTCGTTCTAGTGACCGAAAATAATATTCCAAAGCTTTTGTGTGTTCTCCATTACTTGTGTGAATAAGGCCTATATTATAGAGTATATAACTTCGATCATAGGGATCAATTTCTAGCC-3'; the corresponding wild-type sequence is shown in SEQ ID No.18, specifically: 5'-AAGCCTATTTCTAGTAAATTTTTTTCTTTTTTTTTTTTTTTTTTTTCTTTCTATAGTGGAGATAGTCGCACGTAATGACAGATCACGGCCATATTATTTAAAGCTTGTGGTAAGAAGGGGTT TCGTTCTAGTGCCCGAAAATAATATTCCAAAGCTTTTGTGTGTTCCCATTACTTGTGTGAATAAGGCCTATATTATAGAGTATATAACTTCGATCATAGGGATCAATTTCTAGCC-3'; the s69276-G / T is located in SEQ. ID The sequence shown in No. 19 has a G / T natural variation at position 196 from the 5' end, specifically: 5'-CAAAAAAATAATTGGAATAATCAATATTCGCGATGCACGTATTCGTGTATTAGATCCAAAAGGTCTTTCTTGACACTTAACTACAAAGAAAAGAAGTTTTTTTTAATATGGAAAAAATGTAAAATCTAGAAAAGAAAAAGATAATCGGAATTACTAGTCTTAGAATCTATTTGTACAAAAACAAAAGAAAGATTTGATTTTTTCGAGTGATCTGATCATGCGATACTTTTTTTTCTTATTCGAGATATTGTGTGAATGAATCTACTAATAAATCTAATAA-3'; the corresponding wild-type sequence is shown in SEQ ID As shown in No. 20, specifically: 5'-CAAAAAAATAATTGGAATAATCAATATTCGCGATGCACGTATTCGTGTATTAGATCCAAAAGGTCTTTCTTGACACTTAACTACAAAGAAAAGAAGTTTTTTTTTAATATGGAAAAAATGTAAAATCTAGAAAAGAAAAAGATA ATCGGAATTACTAGTCTTAGAATCTATTTGTACAAAAACAAAAGAAAGATTTTATTTTTTTCGAGTGATCTGATCATGCGATACTTTTTTTTCTTATTCGAGATATTGTGTGAATGAATCTACTAATAATCTAATAA-3'; the s69443-C / A is located in SEQ ID No.The sequence shown in 21 has a C / A natural variation at position 143 from the 5' end, specifically: 5'-CGATACTTTTTTTTCTTATTCGAGATATTGTGTGAATGAATCTACTAATAAATCTAATAAGAGTTAAATTCAATTTTAAAATGAAAAATGACACAAATTCGATACAAATATACAAATAAATAGTAAATAAATATTAAAAATCCATTTTAAAATAAACTAAAAGAAACTATCAAAATAGAAATGATTTTCCAATTTTATTTTGTAGTGATTTCCGTAGTAATTCAAAGAAAGTTTTTTGAATAAAGTTATACAACACAGTTTTTATTACTT-3'; the corresponding wild-type sequence is shown in SEQ ID As shown in No. 22, specifically: 5'-CGATACTTTTTTTTCTTATTCGAGATTTGTGTGAATGAATCTACTAATAATCTAATAAGAGTTAAATTCAATTTTAAAATGAAAAATGACACAAATTCGATACAAATATACAAATAAATAGTAAATAAATATTAAA AATCAATTTTAAAATAAACTAAAAGAAACTATCAAAATAGAAATGATTTTCCAATTTTATTTTGTAGTGATTTCCGTAGTAATTCAAAGAAAGTTTTTTGAATAAAAGTTATACAACACAGTTTTTATTACTT-3'; the s74109-G / A is located in SEQ. ID The sequence shown in No.23 has a G / A natural variation at position 129 from the 5' end, specifically: 5'-ATTTCCCTGCTTAATGGATAATCATTTCCTACCAATGGGAAATTCTTTCTCGTTTTAAATTGAAAATTGAGGTGATTGAATTTGCACCAATGAAACCATAAATTTGATACACAATAGACGAATATGATGAATCTTTTTTTTTAGATAATGAAGGGGATTCTTCCATTCTATTCTATTCATCCGTACTGATCACAGATAGTGGAAAATTTT-3'; the corresponding wild-type sequence is shown in SEQ ID No.24, specifically: 5'-ATTTCCCTGCTTAATGGATAATCATTTCCTACCAATGGGAAATTCTTTCTCGTTTTAAATTGAAAATTGAGGTGATTGAATTTGCACCAATGAAACCATAAATTTGATACACAATAGACGAATATGATAAATCTTTTTTTTTAGATAATGAAGGGGATTCTTCCATTCTATTCTATTCATCCGTACTGATCACAGATAGTGGAAAATTTT-3'; the s79616-G / T is located at SEQ ID The sequence shown in No. 25 has a G / T natural variation at position 116 from the 5' end, specifically: 5'-CGTGCAATTTCTTTTTTTTTTTCTGTATTTCTGGAATATGAGTGTGCGACTTGTTATAATTTATCCTATTGATAGTACAGAGAGTGGGTCTGTCATCTTGATAGAGATGATTTTTGACCTCGTCAGGTATTTATTCTAGTATCTTTTATTCTAGTATCTGGAGCACGGAATATATGAAATAGATTACGAAGTATTAGAAC-3'; the corresponding wild-type sequence is shown in SEQ ID No.26, specifically: 5'-CGTGCAATTTCTTTTTTTTTTTCTGTATTTCTGGAATATGAGTGTGCGACTTGTTATAATTTATCCTATTGATAGTACAGAGAGTGGGTCTGTCATCTTGATAGAGATGATTTTTTACCTCGTCAGGTATTTATTCTAGTATCTTTTATTCTAGTATCTGGAGCACGGAATATATGAAATAGATTACGAAGTATTAGAAC-3'; the s86288-C / T is located in SEQ ID No.The sequence shown in 27 has a C / T natural variation at position 128 from the 5' end, specifically: 5'-TTCTCGCTATATTTTCTGCTACTCCGCCCATTTCATAAAGTATTCTACCTGGTTTAACGACAGCTACCCAATATTCGGGAGATCCTTTCCCCGAACCCATACGTGTTTCCGTAGGTCTTAAAGTAACCGGTTTGTCGGGAAATATGCGTACCCATATTTTTCCACCGCGGCGTGCATTTCGTGTCATTGCTCGTCGCCCCGCTTCTATTT-3'; the corresponding wild-type sequence is shown in SEQ ID As shown in No. 28, specifically: 5'-TTCTCGCTATATTTTCTGCTACTCCGCCCATTTCATAAAGTATTCTACCTGGTTTAACGACAGCTACCCAATATTCGGGAGATCCTTTCCCCGAACCCATACGTGTTT CCGTAGGTCTTAAAGTAACTGGTTTGTCGGGAAATATGCGTACCCATATTTTTCCACCGCGGCGTGCATTTCGTGTCATTGCTCGTCGCCCCGCTTCTATTT-3'; the s119902-C / A is located in SEQ ID The sequence shown in No.29 has a C / A natural variation at position 132 from the 5' end, specifically: 5'-CAAAAAAAGAAAAAAGAGAAGACACAAAGTTTTATCTTTCTTTTTATTATTAGGATATTTTATCATTTTCAGGATAGGGATTATTATTTTCCCCATCGATCCATTTGTCACAATCACAATAATAAATATTCCTAAGTTTTGTCTTTTTTATATTTATACTATTTGAATATAAAAGAGCAAATCTTTGAATGATCTACTTTAATCATATGCTTGGGCTAGAAAATGGATCA-3'; the corresponding wild-type sequence is shown in SEQ ID No.30 is shown, specifically: 5'-CAAAAAAAGAAAAAAGAGAAGACACAAAGTTTTATCTTTCTTTTTATTAGGATATTTTATCATTTTCAGGATAGGGATTATTTTCCCCATCGATCCATTTGTCACAATCACAATA ATAAATATTCATAAGTTTTGTCTTTTTTATATTTATACTATTTGAATATAAAAGAGCAAATCTTTGAATGATCTACTTTAATCATATGCTTGGGCTAGAAAATGGATCA-3'; the s122256-C / T is located at SEQID The sequence shown in No.31 has a C / T natural variation at position 146 from the 5' end, specifically: 5'-GGGAAGCTAGTGATAAGATACTGAATGTCGTGAATATTTTTGGCATTGAGGTAGCCATTCCACCCATTTCATCAAGATAAACACGACGTATTCTATCATAACCCGTTCCTGCCAAGAAAAAAAGTGCGGCACCAATAAATCCATGCGATATTATTTGTAAAATGGCTCCATTGAGTCCCATATCACTTATAGAGCAAATTCCTATAATTATGAAACCCATATGAGATACAGAAGAATAGGCTATTCTTTTTTTTAAATTTCGTTGACCAGGAGATGTTGA-3'; the corresponding wild-type sequence is shown in SEQ Shown in IDNo.32, specifically: 5'-GGGAAGCTAGTGATAAGATACTGAATGTCGTGAATATTTTTGGCATTGAGGTAGCCATTCCACCCATTTCATCAAGATAAACACGACGTATTCTATCATAACCCGTTCCTGCCAAGAAAAAAAGTGCGGCACCAATAAATCC ATGTGATATTATTTGTAAAATGGCTCATTGAGTCCCATATCACTTATAGAGCAAATTCCTATAATTATGAAACCCATATGAGATACAGAAGAATAGGCTATTCTTTTTTTAAATTTCGTTGACCAGGAGATGTTGA-3'; the s129310-G / A is located in SEQ. ID No.The sequence shown in 33 has a G / A natural variation at position 170 from the 5' end, specifically: 5'-AATAATCCCAACGTGTTACATAGGGCAAATATTGTATAATTGTTCGATTTTCCGCAATTTTTTCCATTCCTCTGTGTAAATAACCTAATATTGGTTCACAGTCAATAACATCTTCCCCGTCTAGAGTAACGATGAGGCGAAGAACACCATGCATTGATGGGTGGTGGGGGCCCATATTAACTATCATAAGGTCTTTTCGTGTAGCTGGTACATTCATAGGGTTTCCTTGATCCATTCGTCCATGAATTACTGAAAACAAGTTCATCTCAAGATCTAATAAATCAAATAATTCAA-3'; the corresponding wild-type sequence is shown in SEQ ID As shown in No. 34, specifically: 5'-AATAATCCCAACGTGTTACATAGGGCAAATATTGTATAATTGTTCGATTTTCCGCAATTTTTTCCATTCCTCTGTGTAAATAACCTAATATTGGTTCACAGTCAATAACATCTTCCCGTCTAGAGTAACGATGAGGCGAAGAACACCATGCA TTGATGGGTGGTGGGGACCCATATTAACTATCATAAGGTCTTTTCGTGTAGCTGGTACATTCATAGGGTTTCCTTGATCCATTCGTCCATGAATTACTGAAAACAAAAAGAAGTTCATCTCAAGATCTAATAAATCAAATAATTCAA-3'; the s133745-T / C is located in SEQ. ID No.The sequence shown in 35 has a T / C natural variation at position 115 from the 5' end, specifically: 5'-AAAAGAGGAGAATGCACATTTGCTTGAAACAGTTCCCAAATAGCTATTTTGCGTCTTTGTGCGCGCACGGATCCTTTTATTATGTCTCGGCGAAAATCCGATTGTTGTGAATAATGTATCAAAGTCACTTCTTCTATTTGATCAGAATTCGTTGTATCTTTGGTATTATTATAAGTATCTGTATTTTTTTTATTATCAGTAAAAATCACTACCCGTTTGGCTTTTCGTGAACGAATTTCATAATCTTCCG-3'; the corresponding wild-type sequence is shown in SEQ ID No.36, specifically: 5'-AAAAGAGGAGAATGCACATTTGCTTGAAACAGTTCCCAAATAGCTATTTTGCGTCTTTGTGCGCGCACGGATCCTTTTATTATGTCTCGGCGAAAATCCGATTGTTGTGAATAACGTATCAAAGTCACTTCTTCTATTTGATCAGAATTCGTTGTATCTTTGGTATTATTATAAGTATCTGTATTTTTTTTATTATCAGTAAAAATCACTACCCGTTTGGCTTTTCGTGAACGAATTTCATAATCTTCCG-3'; wherein the bold bases represent the sites of the SNP markers; the SNP marker combination is related to the latex yield of the rubber tree; the base of the s8165-T / G site is T, the base of the s9952-G / A site is G, and the base of the s235 The base at the s07-G / A site is G, the base at the s26511-A / G site is A, the base at the s28467-C / A site is C, the base at the s29302-T / G site is T, the base at the s29303-C / A site is C, the base at the s29894-G / T site is G, the base at the s43759-T / C site is T, the base at the s46966-A / C site is A, the base at the s69276-G / T site is G, and the base at the s43759-T / C site is T. When the base of the s69443-C / A site is C, the base of the s74109-G / A site is G, the base of the s79616-G / T site is G, the base of the s86288-C / T site is C, the base of the s119902-C / A site is C, the base of the s122256-C / T site is C, the base of the s129310-G / A site is G, and the base of the s133745-T / C site is T, it indicates that the rubber tree is a high-yield variety.

[0006] Preferably, the high yield means latex yield>50mL / knife.

[0007] The present invention also provides the application of the SNP marker combination described in the above technical solution in the breeding of high-yield rubber tree varieties.

[0008] The present invention also provides a method for identifying high-yield varieties of rubber trees based on SNP markers of rubber tree chloroplasts, comprising the following steps: amplifying DNA of a sample to be tested, performing SNP site polymorphism detection using the SNP marker combination described in the above technical solution, and performing genotyping on the sample to be tested according to the detection result, wherein the base of the s8165-T / G site is T, the base of the s9952-G / A site is G, the base of the s23507-G / A site is G, the base of the s26511-A / G site is A, the base of the s28467-C / A site is C, the base of the s29302-T / G site is T, the base of the s29303-C / A site is C, and the base of the s29894-G / A site is T. -G / T site is G, the base at the s43759-T / C site is T, the base at the s46966-A / C site is A, the base at the s69276-G / T site is G, the base at the s69443-C / A site is C, the base at the s74109-G / A site is G, the base at the s79616-G / T site is G, the base at the s86288-C / T site is C, the base at the s119902-C / A site is C, the base at the s122256-C / T site is C, the base at the s129310-G / A site is G and the base at the s133745-T / C site is T, it indicates that the sample is a high-yield variety of rubber tree.

[0009] The present invention also provides the use of the SNP marker combination described in the above technical solution in the preparation of a product for identifying high-yield varieties of rubber trees.

[0010] Preferably, the product comprises a kit and / or a gene chip.

[0011] The present invention also provides a kit for identifying high-yield varieties of rubber trees, wherein the kit comprises the SNP marker combination described in the above technical solution.

[0012] The present invention also provides the application of the SNP marker combination described in the above technical solution in constructing a rubber tree gene library.

[0013] Beneficial effects of the present invention: The present invention retrieves NCBI and my country's National Genome Science Data Center, and finds that there are hundreds of different varieties and wild germplasm genome resequencing data, among which the project number PRJCA004986 released in September 2023 has the largest data volume, including more than 100 cultivars and more than 200 wild germplasm materials, and by mining and analyzing these data, it is conducive to screening effective SNP marker sites, which lays the foundation for further assisting rubber tree high-yield resistance breeding. Using the SNP marker combination provided by the present invention, variety improvement and new variety breeding can be carried out more accurately, and then the higher-yield improved varieties are selected to improve the planting efficiency of rubber trees. BRIEF DESCRIPTION OF THE DRAWINGS

[0014] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required to be used in the embodiments are briefly introduced below.

[0015] Figure 1 This is the evolutionary tree of 670 chloroplast genome sequences of 335 samples in Example 1; Figure 2 This is a cluster analysis heat map of 299 core SNP sites in 335 samples in Example 2; Figure 3 This is a heat map of the cluster analysis of 25 SNP sites in Example 2. DETAILED DESCRIPTION

[0016] To further illustrate the present invention, a method for identifying a type B cultivar of rubber tree based on rubber tree chloroplast SNP markers provided by the present invention is described in detail below in conjunction with the accompanying drawings and examples, but they should not be construed as limiting the scope of protection of the present invention.

[0017] The present invention downloaded 335 genome resequencing data of rubber tree cultivars and wild germplasms from the National Genome Science Data Center, used GetOrganelle software to assemble chloroplasts of 127 cultivars and 208 wild germplasms, a total of 335 materials, and then used the chloroplast genome of the main rubber tree cultivar Reyan 73397 as a reference, used snippy software to screen and analyze core SNPs, and finally used the Heatmap analysis tool of TBtools software to cluster analyze the core SNPs, and screened out the SNP sites between the cultivars and the wild germplasms. The methods described in the following embodiments are all conventional methods unless otherwise specified.

[0018] Example 1 335 materials chloroplast assembly The genome resequencing data of 335 rubber tree cultivars and wild germplasms were downloaded from the National Genome Science Data Center (https: / / ngdc.cncb.ac.cn / ) to the local server (Project No.: PRJCA004986). The chloroplast genomes of each sample were assembled using the GetOrganelle software. The specific command was: get_organelle_from_reads.py -1 CRR286406_f1.fastq.gz -2 CRR286406_r2.fastq.gz -F embplant_pt -oCRR286406 -R 15 -t 20 --reverse-lsc. The chloroplast genomes of 335 samples were obtained. The chloroplast genome of each sample included two configurations, namely, embplant_pt.K85.complete.graph1.1.path_sequence.fasta and embplant_pt.K85.complete.graph1.2. path_sequence.fasta, and named them CRR286406_1 and CRR286406_2 respectively (other samples were named in the same way). The specific commands were: ll|grepCRR|awk '{print"cat "$NF" / *complete.graph1.1.path_sequence.fasta>Result / "$NF"_1"}'|sh and ll|grep CRR|awk '{print"cat "$NF" / *complete.graph1.2.path_sequence. fasta>Result / "$NF"_2"}'|sh. Then, the cat command was used to merge the 670 sequences of 335 samples into a fasta file named CRR286-cp.fasta. The mafft software was used for multiple sequence alignment. The command was: mafft CRR286-cp.fasta>CRR286-cp.fasta.mafft. Then, the fasttree software was used to construct the evolutionary tree. The command was: fasttree CRR286-cp.fasta.mafft>CRR286-cp.fasta.mafft.tree. Transfer the obtained evolutionary tree file CRR286-cp.fasta.mafft.tree to a personal computer and use Figtree software to view and display the evolutionary tree, such as Figure 1 As shown, a total of 335 chloroplast genome sequences from the branch where CRR286406_2 was located were selected for subsequent SNP analysis.

[0019] Example 2 SNP analysis of chloroplast genomes in 335 samples According to the phylogenetic tree of the chloroplast genome sequence constructed in Example 1, a total of 335 chloroplast genome sequences in the branch where CRR286406_2 is located are selected, and core SNP analysis is performed using the software snippy. The specific command is: snippy--outdirCRR286364_1a --ref . / CRR286406_2 --ctgs CRR286364_1 --cpus 10. The command takes CRR286364_1 as an example, and the chloroplast genome sequence CRR286406_2 of the cultivated variety Reyan 73397 is used as the reference genome. The above command analysis is performed on the 335 chloroplast genome sequences respectively; then the core SNP analysis is performed, and the specific command is: snippy-core CRR286348_1a CRR286348_2a CRR286349_1a CRR286349_2a CRR286350_1a...(335 folder names in total), the analysis results are the core.tab file, which is the core SNP matrix file, including 299 core SNP sites. Use the sed command to replace A, T, C, and G with numbers 1, 2, 3, and 4. The specific command is: morecore.tab|sed -e 's / A / 1 / g' -e 's / T / 2 / g' -e 's / C / 3 / g' -e 's / G / 4 / g'>core.tab.xlsx, transfer core.tab.xlsx to a personal desktop, open it with Excel, select all and selectively paste in a new table, check transpose, delete the first row to get the SNP site matrix. Use the Heatmap analysis tool of the TBtools software to draw a heat map for the transposed core.tab.xlsx file and perform cluster analysis on the columns, such as Figure 2 As shown in the figure, Cul is a high-yield resistant material selected and cultivated from wild germplasm, and cultivated on a large scale; WAC, WRO and WMG represent wild germplasm materials collected from three regions of Brazil. Through comparative analysis, it was found that 25 SNP loci can divide high-yield varieties into two complementary types, such as Figure 3As shown, the two complementary types are named A class and B class, and a SNP site covering both A and B types is named AB class SNP site, of which A class includes 5 SNP sites, B class includes 19 SNP sites, and AB class includes 1 SNP site. Among them, the 19 SNP sites of type B are represented by the structure of s-Reyan 73397 chloroplast genome physical position-cultivated variety genotype / wild germplasm genotype, and the variation information is as follows: s8165-T / G, s9952-G / A, s23507-G / A, s26511-A / G, s28467-C / A, s29302-T / G, s29303-G / A, s29304-G / A, s29305-G / A, s29306-G / A, s29307-G / A, s29308-G / A, s29309-G / A, s29303 ... -C / A, s29894-G / T, s43759-T / C, s46966-A / C, s69276-G / T, s69443-C / A, s74109-G / A , s79616-G / T, s86288-C / T, s119902-C / A, s122256-C / T, s129310-G / A and s133745-T / C.

[0020] The 19 B-type SNP loci were used to identify the wild-type germplasms numbered WAC-1~WAC-33, WRO-1~WRO-53, WMG-1~WMG-29 and the cultivated variety samples numbered Cul-1~Cul-127. The base of the s8165-T / G site was T, the base of the s9952-G / A site was G, the base of the s23507-G / A site was G, the base of the s26511-A / G site was A, the base of the s28467-C / A site was C, the base of the s29302-T / G site was T, the base of the s29303-C / A site was C, the base of the s29894-G / T site was G, and the base of the s43 When the base at the 759-T / C site is T, the base at the s46966-A / C site is A, the base at the s69276-G / T site is G, the base at the s69443-C / A site is C, the base at the s74109-G / A site is G, the base at the s79616-G / T site is G, the base at the s86288-C / T site is C, the base at the s119902-C / A site is C, the base at the s122256-C / T site is C, the base at the s129310-G / A site is G and the base at the s133745-T / C site is T, it indicates that the sample is a high-yield variety of rubber tree, type B; the rest are type b, and the latex yield of type B rubber tree is higher than that of type b. The identification results are shown in Tables 1 and 2, where Cul is the cultivated variety material, and WAC, WRO and WMG represent wild germplasm materials collected in three regions of Brazil.

[0021] Table 1 Identification results of WAC-1~WAC-33, WRO-1~WRO-53, WMG-1~WMG-29

[0022] Table 2 Identification results of Cul-1~Cul-127

[0023] As shown in Table 1, among the 115 wild-type samples, 2 samples were detected as type B. As shown in Table 2, among the cultivated variety samples numbered Cul-1 to Cul-127, 87 samples were type B. The results are consistent with Figure 3 It can be seen that the SNP marker combination provided by the present invention can be used to identify the B-type rubber tree cultivars, and can also identify potential high-yield and resistant germplasm materials from wild germplasm, which can greatly improve the breeding process.

[0024] Test Example 1 Some varieties were extracted, and their latex yields were collected, which were expressed as the volume of latex that could be collected per cut, and then the screening results of the present invention were verified, as shown in Tables 3 and 4 below.

[0025] Table 3 Latex yield of samples WAC-1~WAC-33, WRO-1~WRO-53, WMG-1~WMG-29

[0026] Table 4 Latex yield of samples Cul-1~Cul-93

[0027] From Tables 1 to 4, it can be seen that among the 242 samples, 89 B-type cultivars were obtained by screening the SNP loci of the present invention, of which 57 samples had a latex yield of >50 mL / knife, belonging to high-yield varieties.

[0028] There are certain discrepancies between the screening results and the yield results. This is because there are many factors affecting the yield, such as climatic conditions, pests and diseases, planting management methods and rubber tapping techniques, which can lead to deviations in the yield results. Among the germplasms screened using the SNP markers of the present invention, 64.1% are high-yield varieties, which can indicate that the SNP marker combination provided by the present invention is related to the latex yield of rubber trees, and the SNP marker combination provided by the present invention can be used to screen high-yield germplasm of rubber trees.

[0029] As can be seen from the above-described embodiments, the present invention collects and downloads chloroplast genome assembly and SNP screening analysis by resequencing data, and obtains a chloroplast genome SNP marker combination that can effectively distinguish cultivar and wild germplasm. The present invention, by cluster analysis of the screened SNP, can divide cultivar into A type, B type and AB type three types, and the SNP marker combination provided by the present invention can be used to identify the rubber tree cultivar of B type. In conjunction with Table 1 to Table 4, it can be seen that in the rubber tree of cultivar B type, 64% are high-yield varieties, and the SNP marker combination provided by the present invention can be used for screening rubber tree high-yield germplasm, and can assist rubber tree high-yield variety breeding.

[0030] Although the above embodiment describes the present invention in detail, it is only a part of the embodiments of the present invention, not all of the embodiments. People can also obtain other embodiments based on this embodiment without creativity, and these embodiments all fall within the protection scope of the present invention.

Claims

1. A SNP marker combination for identifying high-yield varieties of rubber trees, characterized in that: The SNP marker combination includes s8165-T / G, s9952-G / A, s23507-G / A, s26511-A / G, s28467-C / A, s29302-T / G, s29303-C / A, s29894-G / T, s43759-T / C, s46966-A / C, s69276-G / T, s69443-C / A, s74109-G / A, s79616-G / T, s86288-C / T, s119902-C / A, s122256-C / T, s129310-G / A and s133745-T / C; The s8165-T / G is located at the 65th position from the 5' end of the sequence shown in SEQ ID No. 1, and there is a natural T / G variation at this site; The s9952-G / A is located at the 72nd position from the 5' end of the sequence shown in SEQ ID No. 3, and there is a natural G / A variation at this position; The s23507-G / A is located at the 67th position from the 5' end of the sequence shown in SEQ ID No.5, and there is a natural G / A variation at this position; The s26511-A / G is located at the 81st position from the 5' end of the sequence shown in SEQ ID No.7, and there is a natural A / G variation at this position; The s28467-C / A is located at position 87 from the 5' end of the sequence shown in SEQ ID No. 9, and there is a natural variation of C / A at this position; The s29302-T / G is located at position 142 from the 5' end of the sequence shown in SEQ ID No. 11, and there is a natural T / G variation at this site; The s29303-C / A is located at position 143 from the 5' end of the sequence shown in SEQ ID No. 11, and there is a natural variation of C / A at this position; The s29894-G / T is located at position 134 from the 5' end of the sequence shown in SEQ ID No. 13, and there is a natural G / T variation at this position; The s43759-T / C is located at position 139 from the 5' end of the sequence shown in SEQ ID No. 15, and there is a natural T / C variation at this position; The s46966-A / C is located at position 136 from the 5' end of the sequence shown in SEQ ID No. 17, and there is a natural A / C variation at this position; The s69276-G / T is located at position 196 from the 5' end of the sequence shown in SEQ ID No. 19, and there is a natural G / T variation at this position; The s69443-C / A is located at position 143 from the 5' end of the sequence shown in SEQ ID No. 21, and there is a natural variation of C / A at this position; The s74109-G / A is located at position 129 from the 5' end of the sequence shown in SEQ ID No. 23, and there is a natural G / A variation at this position; The s79616-G / T is located at the 116th position from the 5' end of the sequence shown in SEQ ID No. 25, and there is a natural G / T variation at this position; The s86288-C / T is located at position 128 from the 5' end of the sequence shown in SEQ ID No. 27, and there is a natural C / T variation at this position; The s119902-C / A is located at position 132 from the 5' end of the sequence shown in SEQ ID No. 29, and there is a natural variation of C / A at this position; The s122256-C / T is located at position 146 from the 5' end of the sequence shown in SEQ ID No. 31, and there is a natural C / T variation at this position; The s129310-G / A is located at position 170 from the 5' end of the sequence shown in SEQ ID No. 33, and there is a natural G / A variation at this position; The s133745-T / C is located at the 115th position from the 5' end of the sequence shown in SEQ ID No.35, and there is a natural T / C variation at this site; the SNP marker combination is related to the latex yield of rubber trees; the base of the s8165-T / G site is T, the base of the s9952-G / A site is G, the base of the s23507-G / A site is G, the base of the s26511-A / G site is A, the base of the s28467-C / A site is C, the base of the s29302-T / G site is T, the base of the s29303-C / A site is C, the base of the s29894-G / T site is G, and the s43759-T / C When the base of the site is T, the base of the s46966-A / C site is A, the base of the s69276-G / T site is G, the base of the s69443-C / A site is C, the base of the s74109-G / A site is G, the base of the s79616-G / T site is G, the base of the s86288-C / T site is C, the base of the s119902-C / A site is C, the base of the s122256-C / T site is C, the base of the s129310-G / A site is G and the base of the s133745-T / C site is T, it indicates that the rubber tree is a high-yield variety.

2. The SNP marker combination according to claim 1, wherein The high yield means that the latex yield is greater than 50 mL / knife.

3. Application of the SNP marker combination according to claim 1 in the breeding of high-yield rubber tree varieties.

4. A method for identifying high-yield varieties of rubber trees based on rubber tree chloroplast SNP markers, characterized in that, The following steps are involved: The DNA of the sample to be tested is amplified, and the SNP marker combination according to claim 1 is used to detect the polymorphism of the SNP site. According to the detection result, the base of the s8165-T / G site is T, the base of the s9952-G / A site is G, the base of the s23507-G / A site is G, the base of the s26511-A / G site is A, the base of the s28467-C / A site is C, the base of the s29302-T / G site is T, the base of the s29303-C / A site is C, the base of the s29894-G / T site is G, and the base of the s43759- When the base at the T / C site is T, the base at the s46966-A / C site is A, the base at the s69276-G / T site is G, the base at the s69443-C / A site is C, the base at the s74109-G / A site is G, the base at the s79616-G / T site is G, the base at the s86288-C / T site is C, the base at the s119902-C / A site is C, the base at the s122256-C / T site is C, the base at the s129310-G / A site is G and the base at the s133745-T / C site is T, it indicates that the sample is a high-yield variety of rubber tree.

5. Use of the SNP marker combination according to claim 1 in preparing a product for identifying high-yield varieties of rubber trees.

6. The use according to claim 5, characterized in that: The product includes a kit and / or a gene chip.

7. A kit for identifying high-yield varieties of rubber trees, characterized in that: The kit comprises the SNP marker combination according to claim 1.

8. Use of the SNP marker combination according to claim 1 in constructing a rubber tree gene library.

Images