Application of Zm00001eb050010 gene in breeding of corn with coarse ear

Abstract

The invention relates to the field of corn molecular marker assisted breeding, and particularly discloses an application of a Zm00001eb050010 gene in regulation and control of corncob thickness, the invention provides a corncob thickness molecular marker, the molecular marker is a Zm00001eb050010 gene with a sequence shown as SEQ ID NO.1, the molecular marker is used for analyzing SNP255223086 which is located in a chromosome 1 and is remarkably related to the corncob thickness in a plurality of environments by utilizing GWAS analysis, and the molecular marker is a Zm00001eb050010 gene with a sequence shown as SEQ ID NO.3. The molecular marker is a Zm00001eb050010 gene with a sequence shown as SEQ ID NO.4. The molecular marker is a Zm00001eb050010 gene. The functional gene Zm00001eb050010 and SNP255223086 for regulating and controlling the ear thickness can explain 5.23% of the phenotypic variation of the ear thickness, and by means of molecular marker breeding, the generation interval can be shortened, the breeding process can be greatly accelerated, the selection precision can be improved, and the cost can be saved.

Description

Technical Field

[0001] The present invention relates to the field of molecular marker-assisted breeding of corn, and in particular to an application of a Zm00001eb050010 gene in breeding of coarse corn. Background Art

[0002] Corn, as one of the world's most important food crops, plays a crucial role in global agricultural production. With the continued growth of the global population and the increasing demand for food, corn production faces significant challenges. According to a report by the Food and Agriculture Organization of the United Nations (FAO), while some major corn-producing countries (such as the United States, China, and Brazil) have seen increases in yields per unit area, the overall growth rate has slowed, with some regions even experiencing a downward trend. Therefore, increasing corn yields per unit area has become a pressing task for global food security.

[0003] Corn yield per plant is closely related to inflorescence development. Optimizing ear traits, in particular, plays a decisive role in increasing yield per plant. Research has shown that ear diameter, a key component of ear traits, significantly impacts corn yield per plant. Thicker ears can accommodate more kernels per row or larger kernel sizes, directly increasing yield per plant. Therefore, optimizing ear diameter is considered a key breakthrough in increasing corn yield per plant and a key focus of future corn breeding.

[0004] The expression of crop traits depends not only on the effects of individual genes, but also on the effects of interactions between genes (epistasis). Under different environmental conditions, the contribution of genes to traits may be different. The present invention aims to identify new genes and SNPs associated with ear diameter through genome-wide association analysis (GWAS) of tropical maize multi-parent populations using whole genome sequencing (WGS) technology, screen out tag SNPs through haplotype analysis, and reduce redundant testing costs during GS or MAS breeding while retaining the integrity of genetic information, in order to provide innovative genetic resources and precise improvement strategies for significantly improving the yield per plant of maize, thereby promoting the optimization and efficiency of maize production. Summary of the Invention

[0005] To address the shortcomings of existing technologies, the present invention provides an application of the Zm00001eb050010 gene in regulating corn ear diameter, explores the functional gene Zm00001eb050010 closely associated with corn ear diameter, and provides a theoretical basis for molecular marker-assisted selection of corn with thicker ears. To achieve the above objectives, the technical solution of the present invention is implemented through the following technical solutions:

[0006] The present invention provides a maize ear diameter molecular marker, wherein the molecular marker is a Zm00001eb050010 gene with a sequence shown in SEQ ID NO: 1.

[0007] Furthermore, the expression level of the Zm00001eb050010 gene was positively correlated with corn ear length.

[0008] Furthermore, the gene is 255223481, 255223574, 255223613、255223628、255223724、255223793、255223873、 255223939、255224065、255224228、255224255、255224256、 255224338、255224343、255224361、255224427、255224704、 255224706、255224751、255224754、255224910、255224958、 255224996、255225071、255225192、255225569、255225614、 255225621、255225791、255225815、255226193、255226196、 255226341、255226431、255226541、255226746、255226805、 255226932、255226974、255227019、255227065、255227078、 255227143、255227176、255227395、255227427、255227431、

[0018] 255227507, 255227549, 255227672, 255227872, 255227972, the bases are presented in order

[0019] When the expression is GCCAACAGTCGCGGTAGGCGGGCGGGCTTCGCACGCATCAATG AGCGGGAGA, corn shows the dominant trait of thick ear.

[0020] Furthermore, the present invention provides a product for detecting the molecular marker, the product comprising a kit, a reagent or a gene chip, and the product detecting the genotype or expression level of the molecular marker.

[0021] Furthermore, the present invention provides products for detecting the genotype of the molecular marker, including products prepared by PCR, qPCR, Sanger sequencing, high-throughput sequencing, fluorescence in situ hybridization, TaqMan probe method, ARMS-PCR method or KASP method.

[0022] Furthermore, the present invention provides the use of the product in identifying or assisting in identifying the coarseness trait of corn ears.

[0023] Furthermore, the present invention provides the application of the molecular markers in any of the following: a) genetic diversity analysis of coarse-ear corn; b) construction of a molecular genetic map of coarse-ear corn; c) genome-wide association analysis of coarse-ear corn; d) variety identification or assisted identification of coarse-ear corn; e) molecular marker-assisted selection breeding of coarse-ear or high-yield corn; f) genome-wide selection breeding of coarse-ear corn.

[0024] Furthermore, the present invention provides a method for screening coarse-ear corn germplasm, which comprises taking a corn sample to be tested, detecting the molecular markers, and screening germplasm that meets the genotype of the molecular markers, namely, coarse-ear corn germplasm.

[0025] Furthermore, the present invention provides a method for screening high-yield corn germplasm, taking a corn sample to be tested, detecting the molecular markers, and screening germplasm that meets the genotype of the molecular markers, namely, the high-yield corn germplasm.

[0026] The term "Molecular Marker-assisted selection (MAS)" refers to a breeding technique that uses molecular markers of target traits to select offspring lines and thereby obtain superior individual plants containing the target gene.

[0027] The term "Genomic Selection (GS)" refers to a modern breeding technology based on genetic evaluation and selection based on whole-genome marker information. It aims to predict individual breeding values ​​or phenotypic performance through high-density molecular markers, thereby accelerating the breeding process and improving selection efficiency.

[0028] The term "GS-MAS" refers to a breeding strategy that combines genome-wide selection (GS) with marker-assisted selection (MAS). GS-MAS uses high-density marker information covering the entire genome and phenotypes to estimate individual breeding values, while simultaneously linking major and minor genes. This allows for early prediction and selection of complex traits (such as those with low heritability and difficult to measure) based on breeding values, thereby shortening generation intervals, accelerating the breeding process, improving selection accuracy, and reducing costs.

[0029] The technical effects achieved by the present invention are:

[0030] This study constructed a multi-parent maize population with significant differences in ear diameter by using the temperate maize inbred line Ye107, which exhibits a relatively thin ear diameter, as a common parent and crossing it with six tropical and subtropical maize inbred lines with relatively thick ear diameters. GWAS analysis pinpointed SNP_255223086 on chromosome 1, which is significantly associated with ear diameter. Furthermore, the functional gene Zm00001eb050010 regulating ear diameter was identified. SNP_255223086 can explain 5.23% of the phenotypic variation in ear diameter. Haplotype analysis showed that among 1,107 RILs (recombinant inbred lines), Zm00001eb050010 had a total of 13 haplotypes (Hap1, Hap2, ..., Hap13), among which Hap10 had significantly higher ear diameter than the other 12 haplotypes. Therefore, Hap10 of the Zm00001eb050010 gene is a haplotype type that significantly increases ear diameter. The results of the present invention will help further study the regulatory mechanism of corn ear diameter and also provide new major effect genes for the development of corn varieties with larger ear diameters, thereby shortening the generation interval for high-yield corn breeding, accelerating the breeding process, improving selection accuracy, and saving costs. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] Figure 1 : The population structure diagram of 1107 RILs in the embodiment of the present invention; wherein: (a) principal component analysis; (b) unrooted tree; (c) Bayesian clustering diagram of 1107 RILs when K=6; (d) LD decay diagram;

[0032] Figure 2 GWAS results of ear thickness in the embodiment of the present invention; wherein: (a) GWAS results based on the mean ear thickness phenotype in Yanshan in 2022; (b) GWAS results based on the mean ear thickness phenotype in Yanshan in 2023; (c) GWAS results based on the mean ear thickness phenotype in Jinghong in 2024; (d) GWAS results based on the BLUP value of ear thickness;

[0033] Figure 3 Schematic diagram of the relationship between significant SNPs and candidate genes in the embodiment of the present invention;

[0034] Figure 4 Candidate gene LD block analysis diagram;

[0035] Figure 5 Schematic diagram of the mutation sites that form different haplotypes in the Zm00001eb050010 gene;

[0036] Figure 6 The differences between different haplotypes in subpopulations are shown, where A is the difference percentage and B is the difference box plot.

[0037] Figure 7 The expression level of Zm00001eb032370 gene in different tissues. DETAILED DESCRIPTION

[0038] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention are clearly and completely described below in conjunction with the embodiments of the present invention. For purchased commodities in the test method, if no specific conditions are specified, they are carried out according to conventional conditions or the conditions recommended by the manufacturer. If the manufacturer of the reagents or instruments used is not specified, they can all be conventional products purchased from the market.

[0039] The gene sequence of Zm00001eb050010 marked by the present invention is shown in SEQ ID NO: 1, which corresponds to positions 255223470-255228057 from the 5′ end of chromosome 1 of the reference genome version Zm-B73-REFERENCE-NAM-5.0.

[0040] Example 1

[0041] 1. Experiment

[0042] 1.1 Plant materials and experimental design

[0043] The present invention used the temperate, high-quality, backbone maize inbred line Ye107 as the common male parent and hybridized it with six tropical and subtropical maize inbred lines with relatively thick ears (CML312, CML395, YML46, YML32, NK40-1, and YML1218) to form six F1 hybrid combinations. These combinations were then selfed continuously to F9 lines using single-kernel descent, resulting in a multi-parent population containing 1,107 recombinant inbred lines. The number of recombinant inbred lines contained in each subpopulation (pop1-pop6) was 142, 148, 303, 178, 200, and 136, respectively (Table 1).

[0044] This experiment was conducted in Yanshan County, Yunnan Province, China, in 2022 and 2023, and in Jinghong City, Yunnan Province, China, in 2024. A completely randomized block design was used, with three replications per site. The experimental material was planted in rows with a length of 3 m, a spacing of 0.70 m, 14 plants per row, and two rows per plot. Planting was managed according to local standard agronomic practices. After ripening, nine ears were randomly selected from each row to record their ear diameter, and the average value was calculated.

[0045] Table 1 Parental information

[0046]

[0047]

[0048] 1.2 Phenotypic data analysis

[0049] First, the phenotypic data were preprocessed using Excel 2019. SPSS 26.0 was used to assess whether the phenotypic data conformed to a normal distribution, and basic statistics were calculated to help intuitively understand the central tendency and dispersion of the data. R4.3.2 was used to draw heat maps and assess the correlation between phenotypic data under different subpopulations and environmental conditions. A two-way analysis of variance was performed using the aov() function to assess whether the effects of subpopulations and environments on the phenotypic data were statistically significant. The broad-sense heritability was calculated using the lme4 software package using the formula H = Vg / (Vg + (Ve / L)), where Vg represents the genetic variance, Ve is the residual variance, and L is the number of environments.

[0050] 1.3 Whole-genome sequencing

[0051] 1.3.1 Sample collection and DNA extraction

[0052] Healthy maize leaves were sampled between the seedling and jointing stages and immediately frozen in liquid nitrogen to prevent tissue and DNA degradation. Leaf DNA was extracted using the Tiangen Plant Genomic DNA Extraction Kit (Tiangen Biotech Co., Ltd., Beijing, China). The extracted DNA samples were subjected to quality control using agarose gels to confirm concentration (above 100 ng / μl) and purity (A260 / A280 ratio, typically between 1.8 and 2.0), ensuring that the extracted DNA was suitable for subsequent whole-genome sequencing.

[0053] 1.3.2 Library construction and sequencing

[0054] DNA fragmentation was performed using a Vibra-Cell series ultrasonic breaker (Sonics & Materials, Inc., Connecticut, America), and the DNA samples after the break were subjected to gel electrophoresis to assess quality (ideal DNA fragment size should be between 200-500 bp). End repair was performed using T4 DNA polymerase (Thermo Fisher Scientific Inc., Waltham, Massachusetts, America), and short DNA adapters were connected to the two ends of the DNA fragments. In order to ensure the quality of the library, DNA fragments without adapters were removed using Exonuclease I enzyme (Thermo Fisher Scientific Inc., Waltham, Massachusetts, America). DNA fragments of 200-500 bp in length (including the portion of the adapter sequence) were selected by the Agencourt AMPure XP magnetic bead method. The library was PCR amplified to increase the number of target DNA fragments. Finally, the concentration of the library was quantitatively detected using Qubit, and the target concentration was usually between 10-20 nM. The size distribution of the library was detected using a Bioanalyzer (Agilent Technologies, Santa Clara, California, America), and the fragment size was generally 200-500 bp.

[0055] All samples were sequenced using the Illumina NovaSeq 6000 platform (Illumina, Inc., San Diego, California, America) in paired-end sequencing mode at a sequencing depth of 5×. Sequencing results were stored as FASTQ files. The quality of sequencing data was assessed using FastQC (Babraham Institute, Bioinformatics Group, Cambridge, United Kingdom) based on the following criteria: (1) Q-score ≥ 20, (2) genomic GC content of 40% to 60%, (3) repetitive sequence proportion less than 30%, (4) N-containing proportion less than 1–2%, and (5) read length consistent with expected length.

[0056] 1.4 Group stratification analysis

[0057] 1.4.1 Principal Component Analysis

[0058] To assess population genetic structure and control for the effects of population stratification on GWAS results, we performed principal component analysis (PCA) using GCTA software. First, the genotype data were standardized, and a genetic similarity matrix between samples was calculated. The first few principal components were then extracted from this matrix, which was used to characterize population structure.

[0059] 1.4.2 Unrooted Tree Construction

[0060] To further analyze the genetic relationships between subpopulations, we used MEGA (Molecular Evolutionary Genetics Analysis) software to construct unrooted trees to assess the genetic distances and relatedness between subpopulations. During the construction process, we used the neighbor-joining (NJ) method, which efficiently and accurately reflects the genetic differentiation of populations by minimizing the sum of genetic distances between samples.

[0061] 1.4.3 Analysis of population genetic structure

[0062] We used ADMIXTURE software for population stratification analysis. This method assumes that each individual is composed of a mixture of components from multiple subpopulations and infers the proportion of contributions from different subpopulations to each individual's genome. We first selected and set the expected number of subpopulations (K value) based on the study objectives. We then tried different K values ​​and selected the optimal K value through cross-validation to avoid overfitting or underfitting. We analyzed the population structure under different K values ​​and determined the genetic background of the samples based on the population stratification pattern.

[0063] 1.4.4 Linkage disequilibrium analysis

[0064] We used PopLDdecay software to assess the genetic linkage disequilibrium (LD) of the population. The purpose of LD analysis is to measure the degree of linkage disequilibrium between SNPs, that is, whether the genotypes between a pair of markers are independently inherited. In particular, r 2 The value is often used to indicate the degree of LD between SNPs, r 2 The higher the value, the stronger the linkage between the two SNPs and the closer the genetic association. PopLDdecay software was used to calculate the r between different SNP markers. 2 We further analyzed LD decay at different genetic distances using the Plot_OnePop.pl script included with PopLDdecay to generate visualizations of LD decay. These plots clearly demonstrate the LD structure at different positions within a population, identify regions of high LD within the genome, and further aid in marker selection and mapping of genetic regions for GWAS.

[0065] 1.5 Genome-wide association analysis

[0066] A mixed linear model (MLM) was used in GEMMA software to analyze the ear diameter trait in maize. In this model, population genetic structure was treated as a fixed effect, and individual kinship was treated as a random effect to correct for the effects of population structure and kinship on trait expression. To reduce false positives and false negatives, a Bonferroni correction was used for multiple comparisons. The P values ​​of all individual SNPs were compared to the Bonferroni-corrected significance level. If the P value was less than p, the SNP was considered statistically significant. Furthermore, to ensure the reliability and accuracy of the results, an FDR correction was used to reduce the error caused by multiple comparisons. Generally, a p-value less than 5×10^-8 was considered significant. The identified significant SNPs were functionally annotated using the ANNOVAR tool to determine the genomic location, region, and mutation type of the variant loci. Furthermore, PLINK software was used to generate Manhattan plots and QQ plots to visualize the association analysis results. Manhattan plots assess the strength of the association between SNPs and the trait, while QQ plots assess the degree to which the p-value distribution deviates from the expected value.

[0067] 1.6 Candidate gene haplotype analysis

[0068] Haplotype analysis of candidate genes was performed using Haploview v4.2 software. For each gene, we analyzed the frequency distribution, phenotypic differences, and statistical significance of different haplotypes. Based on frequency distribution considerations, boxplots only display haplotypes with high frequency in the sample. Haplotypes with low frequency (≤1%) are not individually displayed; these low-frequency haplotypes have a minimal statistical impact on the final results and do not affect the interpretation of the main conclusions.

[0069] 2. Results

[0070] 2.1 Analysis of corn ear thickness phenotype data

[0071] The results showed that pop4 consistently exhibited a high and stable ear thickness trait in all environments, while pop1 exhibited the lowest average ear thickness (Table 2). The ear thickness phenotypic data of the six subpopulations were statistically analyzed, as shown in Table 2 below. The coefficient of variation (CV%) of ear thickness among the six subpopulations in each environment had a large range of data, indicating significant differences in ear thickness among the subpopulations. The absolute values ​​of the skewness and kurtosis of ear thickness in different environments were generally less than 1. The broad-sense heritability of the number of ear rows in the six subpopulations ranged from 70.15% to 97.28%, indicating a strong genetic influence on this trait. The consistency of phenotypic variation in different environments highlights the reliability of the phenotypic data used for subsequent analysis. Significant differences in ear thickness were observed among subpopulations, influenced by genetic and environmental factors. This variability provides a solid foundation for further GWAS analysis and identification of important loci associated with ear thickness.

[0072] Table 2 Analysis of ear thickness phenotype data

[0073]

[0074]

[0075] Note: 22YS refers to Yanshan in 2022, 23YS refers to Yanshan in 2023, and 24JH refers to Jinghong in 2024

[0076] 2.2 Group stratification analysis

[0077] By observing the 3D PCA graph, the cluster distribution of the samples is relatively clear. The data can be roughly divided into 6 clusters, each representing a potential subgroup. Although most subgroups form a clear cluster structure in the 3D space, there are overlapping areas between some clusters. For example, there is a significant overlap between pop3 and pop6, which may indicate that the two subgroups have a high similarity in the main feature dimensions ( Figure 1 a).

[0078] In the unrooted tree analysis, although there is a certain degree of confusion between different clusters, the overall data can still be clearly divided into 6 main clusters ( Figure 1 b) This overlap may indicate that these clusters are close in the feature space, reflecting their similarity in some key features.

[0079] In the population structure analysis, when K = 6, the data is clearly divided into 6 subgroups. Although there is still a certain degree of overlap between some samples, the 6 subgroups are relatively clear overall ( Figure 1 c) This overlap indicates that there may be overlapping areas between different subgroups, and some samples may have similar characteristics, making their attribution to certain subgroups unclear. Overall, despite this overlap, these subgroups still reasonably reflect the main structural characteristics of the population.

[0080] LD decay in MPP was evaluated using 16,223,980 valid SNPs. The LD decay plot showed that LD decreased rapidly with increasing physical distance between markers. 2 When the value drops to 0.24, the physical distance of LD decay is about 10 kb ( Figure 1 d).

[0081] 2.3 Genome-wide association analysis of ear diameter traits

[0082] The present invention uses 16,223,980 effective SNPs to conduct GWAS analysis of corn ear diameter. The results show that a total of 371 significant SNPs associated with corn ear diameter were identified on 10 chromosomes in all environments ( Figure 2 In particular, a significant SNP associated with corn ear diameter was consistently identified in all environments (22YS, 23YS, 24JH, BLUP): 1_255223086. Gene screening was performed within the 10kb range upstream and downstream of the associated significant SNP: 1_255223086. We consistently identified a gene Zm00001eb050010 ( Figure 3 , Table 3), and the trait correlation was evaluated with reference to its existing functional annotation. The present invention believes that this is a candidate gene worthy of further study.

[0083] Table 3 GWAS co-localized candidate genes

[0084]

[0085] 2.4 Haplotype analysis of candidate genes

[0086] like Figure 4-6 Haplotype analysis of Zm00001eb050010 (the gene sequence of Zm00001eb050010 is shown in SEQ ID NO:1, corresponding to positions 255223470-255228057 from the 5′ end of chromosome 1 in the reference genome version Zm-B73-REFERENCE-NAM-5.0) revealed 13 distinct haplotypes (Hap1, Hap2, …, Hap13) among 1107 RILs. The boxplot only displays the three haplotypes with the highest frequency in the samples. These haplotypes exhibited varying frequency distributions across the different populations studied. Hap10 exhibited the highest phenotypic values ​​and was widely distributed in Population 2. Furthermore, Hap10 was significantly different from Hap1 and Hap6 (P < 0.05), concluding that Hap10 is the dominant haplotype regulating ear diameter in maize.

[0087] The specific haplotype sites, using Zm-B73-REFERENCE-NAM-5.0 as the reference genome, are as follows from left to right on chromosome 1 starting from the 5′ end:

[0088] 255223481(REF:T / ALT:G)255223574(REF:C / ALT:G)

[0089] 255223613(REF:T / ALT:C)255223628(REF:A / ALT:G)

[0090] 255223724(REF:G / ALT:A)255223793(REF:G / ALT:C)

[0091] 255223873(REF:G / ALT:A)255223939(REF:G / ALT:A)

[0092] 255224065(REF:T / ALT:C)255224228(REF:A / ALT:C)

[0093] 255224255(REF:G / ALT:A)255224256(REF:C / ALT:A)

[0094] 255224338(REF:G / ALT:C)255224343(REF:A / ALT:G)

[0095] 255224361(REF:G / ALT:T)255224427(REF:A / ALT:C)

[0096] 255224704(REF:G / ALT:T)255224706(REF:G / ALT:A)

[0097] 255224751(REF:C / ALT:T)255224754(REF:G / ALT:T)

[0098] 255224910(REF:C / ALT:G)255224958(REF:A / ALT:G)

[0099] 255224996(REF:C / ALT:A)255225071(REF:G / ALT:C)

[0100] 255225192(REF:G / ALT:C)255225569(REF:G / ALT:A)

[0101] 255225614(REF:C / ALT:G)255225621(REF:T / ALT:C)

[0102] 255225791(REF:T / ALT:C)255225815(REF:C / ALT:T)

[0103] 255226193(REF:G / ALT:T)255226196(REF:C / ALT:T)

[0104] 255226341(REF:A / ALT:G)255226431(REF:C / ALT:T)

[0105] 255226541(REF:G / ALT:A)255226746(REF:C / ALT:T)

[0106] 255226805(REF:C / ALT:A)255226932(REF:G / ALT:T)

[0107] 255226974(REF:T / ALT:C)255227019(REF:A / ALT:G)

[0108] 255227065(REF:C / ALT:A)255227078(REF:C / ALT:T)

[0109] 255227143(REF:G / ALT:C)255227176(REF:C / ALT:A)

[0110] 255227395(REF:G / ALT:A)255227427(REF:C / ALT:T)

[0111] 255227431(REF:G / ALT:T)255227507(REF:G / ALT:A)

[0112] 255227549(REF:G / ALT:C)255227672(REF:A / ALT:G)

[0113] 255227872(REF:T / ALT:G)255227972(REF:T / ALT:A)

[0114] Hap1:

[0115] TCTGACAGTCGCCAGCTATGGGCGGGCTCTGCGTGCCGTGCTCCATTACGTA

[0116] Hap2:

[0117] TGTAGCAGCCGCGAGCTATGCACCGGCTTCGCACGTCGTACCGCGCGGGAGA

[0118] Hap3:

[0119] TCTAGGGGCCGCCAGCGGCGCACGGGCTTCGCACGTCGTACCGCATTACGTA

[0120] Hap4:

[0121] TGTAGCAGCCGCGAGCTATGGGAGGGCTTCGCACGTCGTACCGCGTGAGATA

[0122] Hap5:

[0123] TCTAGGGGTCGCGGTAGGCGCACGCGCTTCGCACGCATCAATGAATGACATA

[0124] Hap6:

[0125] TCTAGCAGTCGCGGTAGGCGCACGCGCCTCGCACGTCGTACCGCGCGGGAGA

[0126] Hap7:

[0127] TCTGACAGTCGCCAGCTATGGGCGGAGTTCGCACGTCGTACCGCGCGGGAGA

[0128] Hap8:

[0129] TCTAGGGGTAGCGAGAGGCGCACGGGGTCTGCGTGCCGTGCTCCATTACGTT

[0130] Hap9:

[0131] TCTGGCAGCCGCCAGCGACGCACGCGCTTCGCACGTCGTACCGCGCGGGAGA

[0132] Hap10:

[0133] GCCAACAGTCGCGGTAGGCGGGCGGGCTTCGCACGCATCAATGAGCGGGAGA

[0134] Hap11:

[0135] TCTGGCGATCAACAGCTATTGGAGGAGTTCGCACGTCGTACCGCGCGGGAGA

[0136] Hap12:

[0137] TCTAGGGGTAGCGAGAGGCGCACGGGCTTCGCACGCCGTA CCCGCGGGATT

[0138] Hap13:

[0139] TCTAGGGGTCGCGGTAGGTGCACCCGCTTCGCACGTCGTAC CGCGCGGGAGA

[0140] The expression of candidate gene Zm00001eb050010 in various tissues of maize was analyzed ( Figure 7 The results showed that the gene was highly expressed in maize cob primordium, which further suggested that it was related to the development of maize cob diameter.

[0141] In summary, the identification of the Zm00001eb050010 gene as a molecular marker functional gene for maize ear diameter can be used to develop allele-specific markers. Molecular markers can directly detect crop genotypes and accurately identify the presence or absence of target genes or traits, avoiding the blindness and uncertainty in traditional selection and ensuring the accuracy of offspring selection. It can be applied to genetic resource identification, MAS, genetic map construction, linkage mapping, and GS-MAS in breeding. MAS can quickly and accurately identify crop genotypes, avoiding the tedious, time-consuming, and difficult screening of traditional selection, greatly improving selection efficiency, and accelerating the breeding process.

[0142] The above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present invention.

Claims

1. A molecular marker for corn ear diameter, characterized in that: The molecular marker is the Zm00001eb050010 gene with a sequence as shown in SEQ ID NO:

1.

2. The molecular marker according to claim 1, characterized in that The expression level of Zm00001eb050010 gene was positively correlated with corn ear diameter.

3. The molecular marker according to claim 1, characterized in that The genes are 255223481, 255223574, 255223613, 255223628, 255223724, 255223793, 255223873, 255223939, 255224065, 255224228, 255224255, 255224256, 255 224338, 255224343, 255224361, 255224427, 255224704, 255224706, 255224751, 255224754, 255224910, 255224958, 255224996, 255225071, 255225192, 25522556 9、255225614、255225621、255225791、255225815、255226193、255226196、255226341、255226431、255226541、255226746、255226805、255226932、255226974、25 5227019, 255227065, 255227078, 255227143, 255227176, 255227395, 255227427, 255227431, 255227507, 255227549, 255227672, 255227872, 255227972, the bases are presented in order When GCCAACAGTCGCGGTAGGCGGGCGGGCTTCGCACGCATCAATG AGCGGGAGA, corn shows the dominant trait of thick ear.

4. Detection of the molecular marker product according to claim 1 or 2, characterized in that: The product includes a kit, a reagent or a gene chip.

5. A product for detecting the molecular marker genotype according to claim 3, characterized in that: The products include products prepared by PCR, qPCR, Sanger sequencing, high-throughput sequencing, fluorescence in situ hybridization, TaqMan probe method, ARMS-PCR method or KASP method.

6. Use of the product according to claim 4 or 5 in identifying or assisting in identifying the coarseness trait of corn ears.

7. Use of the molecular marker according to any one of claims 1 to 3 in any of the following: a) Genetic diversity analysis of coarse-ear corn; b) Construction of molecular genetic map of coarse-ear corn; c) Whole-genome association analysis of coarse-ear corn; d) Identification or auxiliary identification of coarse-ear corn varieties; e) Molecular marker-assisted selection breeding of coarse-ear or high-yield corn; f) Whole-genome selection breeding of coarse-ear corn.

8. A method for screening coarse corn germplasm, characterized in that: Take a corn sample to be tested, detect the molecular marker as described in claim 2 or 3, and screen the germplasm that meets the genotype of the molecular marker, that is, the coarse corn germplasm.

9. A method for screening high-yield corn germplasm, characterized in that: Take a corn sample to be tested, detect the molecular marker as described in claim 2 or 3, and screen the germplasm that meets the genotype of the molecular marker, which is the high-yield corn germplasm.

Images