Next Article in Journal
Consensify: A Method for Generating Pseudohaploid Genome Sequences from Palaeogenomic Datasets with Reduced Error Rates
Next Article in Special Issue
Mapping-by-Sequencing via MutMap Identifies a Mutation in ZmCLE7 Underlying Fasciation in a Newly Developed EMS Mutant Population in an Elite Tropical Maize Inbred
Previous Article in Journal
A Model of Hormonal Regulation of Stamen Abortion during Pre-Meiosis of Litsea cubeba
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 1
10.3390/genes11010049
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Association Study (GWAS) for Mesocotyl Elongation in Rice (Oryza sativa L.) under Multiple Culture Conditions

by
Hongyan Liu
1,†,
Junhui Zhan
1,†,
Jiaolong Li
1,
Xiang Lu
1,
Jindong Liu
1,
Yamei Wang
1,*,
Quanzhi Zhao
2,* and
Guoyou Ye
1,3
1
CAAS-IRRI Joint Laboratory for Genomics-Assisted Germplasm Enhancement, Agricultural Genomics Institute in Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
2
Collaborative Innovation Center of Henan Grain Crops and Key Laboratory of Rice Biology in Henan Province, College of Agronomy, Henan Agricultural University, Zhengzhou 450002, China
3
Strategic Innovation Platform, International Rice Research Institute, DAPO Box 7777, Metro Manila, Philippines
*
Authors to whom correspondence should be addressed.
†
Both authors contributed equally to this work.
Genes 2020, 11(1), 49; https://doi.org/10.3390/genes11010049
Submission received: 10 December 2019 / Revised: 25 December 2019 / Accepted: 27 December 2019 / Published: 31 December 2019
(This article belongs to the Special Issue Recent Advances in Genetics and Breeding of Major Staple Food Crops)
Browse Figures
Review Reports Versions Notes

Abstract

:
Mesocotyl is a crucial organ for pushing buds out of soil, which plays a vital role in seedling emergence and establishment in dry direct-seeded rice. However, the genetic mechanisms of mesocotyl elongation remains unclear. In our study, 208 rice accessions were used to identify the SNPs significantly associated with mesocotyl length under various culture conditions, including sand, water and soil. The mesocotyl length ranges from 0 to 4.88 cm, 0 to 3.99 cm and 0 to 4.51 cm in sand, water and soil covering, respectively. A total of 2,338,336 SNPs were discovered by re-sequencing of 208 rice accessions. Genome-wide association study (GWAS) based on mixed linear model (MLM) was conducted and 16 unique loci were identified on chromosomes 1, 2 (2), 3, 4, 5 (2), 6 (2), 7, 8, 9 (2) and 12 (3), respectively, explaining phenotypic variations ranging from 6.3 to 15.9%. Among these loci, 12 were stable across two or more environments. Ten out of the sixteen loci coincided with known genes or quantitative trait locus (QTL), whereas the other six were potentially novel loci. Furthermore, five high-confidence candidate genes related to mesocotyl elongation were identified on chromosomes 1, 3, 5, 9 and 12. Moreover, qRT-PCR analysis showed that all the five genes showed significant expression difference between short-mesocotyl accessions and long-mesocotyl accessions. This study provides new insights into the genetic architecture of rice mesocotyl, the associated SNPs and germplasms with long mesocotyl could be useful in the breeding of mechanized dry direct-seeded rice.

1. Introduction

Rice (Oryza sativa) is one of the most important cereal crops grown worldwide. Dry direct-seeded rice refers to the process of establishing the crop from seeds sown on non-puddled and unsaturated soil; in contrast, seedlings from nursery are transplanted into puddled or submerged soil in transplanted rice [1]. Compared with traditional transplanted rice, dry direct-seeded rice has been proposed as a water-efficient and labor-saving approach, which can reduce the cost of water and labor at about 50% [2,3]. Also, dry direct-seeded rice could efficiently utilize early-season rainfall in drought-prone environments and complete its growth cycle within the wet season in rainfed lowlands [4]. To date, direct seeding has been adopted by more than 25% of the worldwide rice cultivation area [1]. Dry direct-seeded rice is becoming a popular option in Philippines, India, Thailand, Cambodia, Laos and Indonesia in tropical Asia as well as the US, Australia and Latin America [5,6].
However, a few constraints in dry direct-seeded rice, e.g., poor seedling establishment, weed infestation and lodging susceptibility could lead to large reduction in grain yield and quality [7,8,9,10,11]. Seedling establishment is very important for high yielding and controlling weed infestation in direct seeding rice [10]. Meanwhile, deep sowing could be good for root and basal internodes elongation in deeper soil, which could gain the stability of rice plant and increase lodging resistance of direct seeding rice [12]. However, deep sowing is always in contradiction to seedling uniformity and leads to grain yield reduction [13]. Rapid and good seedling establishment is important for weed competitiveness and good harvesting. Mesocotyl, an organ between the coleoptilar node and the basal part of seminal root in young monocot seedlings, plays a key role in pushing buds out of the deep water or soil during germination for successful seedling establishment. Longer mesocotyl can facilitate seedling establishment under deep sown conditions in dry direct-seeded rice [14,15,16,17,18]. Besides, Mgonja et al. (1988) also reported the strong association between mesocotyl length and seedling vigor [19]. However, the molecular mechanisms of rice mesocotyl length variation are poorly understood.
Genome-wide association study (GWAS) based on linkage disequilibrium (LD) has been widely adopted to identify loci significantly associated with important and complex agronomic traits in rice [20,21,22]. Extremely high resolution can be achieved by dense SNPs identified in diverse germplasm panels based on the next-generation sequencing (NGS) or SNP chip approaches [20,21,23,24]. Thousands of rice landraces or cultivars have been screened for mesocotyl length by GWAS and dozens of QTL have been reported in previous studies [17,18,25,26]. Wu et al. (2015) screened 270 rice accessions and 16 loci were identified associated with mesocotyl elongation [17]. Furthermore, 469 Indica accessions were used to measure mesocotyl length, and 23 loci were significantly associated with mesocotyl length [18]. Zhao et al. (2018) evaluated the mesocotyl length of 621 rice accessions and detected 13 QTLs [26]. Also, Sun et al. (2018) have identified three QTL for mesocotyl length from 510 rice accessions [25].
Identifying the QTLs in multi-environments could provide accurate information for gene cloning and molecular breeding. In this study, GWAS based on resequencing was conducted in a set of 208 rice accessions for mesocotyl elongation under sand, water and soil (2.0 cm, 4.0 cm and 6.0 cm). The objectives of this study were: (1) To dissect the genetic architecture of mesocotyl elongation, (2) to identify SNPs significantly associated with mesocotyl length, (3) to search for candidate mesocotyl elongation genes for further study, and (4) to select accessions with longer mesocotyl for the breeding of mechanized dry direct-seeded rice.

2. Materials and Methods

2.1. Plant Materials

The panel for GWAS consisted of two parts of rice germplasms, one part (114 accessions) includes advanced lines from International Rice Research Institute (IRRI) and mega varieties released in southeast Asian countries, such as India, Philippines and Bangladesh. The others (94 accessions) are landraces and released varieties mainly from China. All the 208 accessions could be divided into two ecotypes, Indica (170) and Japonica (38) (Table S1).

2.2. Genotyping

Total genomic DNA for re-sequencing was extracted from young leaves according to the CTAB method. The 208 accessions were genotyped using the Illumina HiSeq 2000 (PE150) (50X) by Berry Genomics Corporation (https://www.berrygenomics.com/) (Beijing, China). The average sequencing depth of each accession genome was 50×. Reads were aligned to the Nipponbare RefSeq (IRGSP-1.0) using BWA-MEM (release 0.7.10). Then, the mapped reads were sorted and duplicated were removed by Picard tools (http://broadinstitute.github.io/picard/). The variants for each accession were called by the GATK best practices (release 3.2-2).

2.3. Measurement of Mesocotyl Length

To evaluate the variation of mesocotyl length (the distance from the basal part of seminal root to the coleoptilar node) at varied sowing depths, two replications were set up, and in each replication 15 good-quality seeds from each accession were sown at a depth of 2, 4 and 6 cm in soil contained by plastic trays. The size of plastic tray contained 50 holes with the size of 9.5 cm depth, 4.5 cm top diameter and 2.1 cm bottom diameter. After sowing, the plastic tray was kept in a plastic pallet with 3-cm-deep soil, and the whole system were maintained in a 30 °C dark incubator and the soil in each pallet was kept water saturated for seed germination and seedling growth. Emergence rate was recorded every day, until the emergence rate of any accession reached 100%. Three days after that, seedlings from each hole were carefully excavated and washed for measuring mesocotyl length using Image J (https://imagej.en.softonic.com/).
To quantify the mesocotyl length under different culture media, sand and water culture were used, in addition to soil culture. Fifteen good-quality seeds were sown at a depth of 6 cm in sand, the culture method and sampling timing in sand was the same as that in soil. As for water culture, 15 good-quality seeds were put on a gauze in each hole of a plastic tray (50 holes with a depth of 4.8 cm, top diameter at 4.1 cm and bottom diameter at 2.1 cm). The plastic tray was put on a plastic pallet containing water in level of 2 cm from the bottom. Each plastic tray was wrapped with tin foil paper to avoid illumination and kept in phytotron with constant temperature of 30 °C for 10 days. After that, ten uniform seedlings out of 15 seeds of each accession were used for mesocotyl length measurement with Image J. Two replications were set up in both culture methods and mean of the two replications was used for GWAS analysis.

2.4. Population Structure and LD Decay Analysis

The population structure of the 208 accessions was analyzed using Admixture 1.3.0 [27] based on 2000 SNP markers with default cross-validation (K value ranged from 2 to 5). An adhoc quantity statistic ΔK was used to predict the real number of subpopulations and k = 2 was chosen. Besides, the principal component analysis (PCA) and neighbor-jointing trees were also used to validate population stratification with the software Tassel v5.0 (https://www.maizegenetics.net/tassel).
LD decay analysis was done for the whole population. LD decay was measured by correlation coefficients (r2) for all pairs of SNPs within 500 Kb that were calculated using PopLDdecay v3.27 (https://github.com/BGI-shenzhen/PopLDdecay) with the following parameters: -MaxDist 500-MAF 0.05-Het 0.88-Miss 0.999.

2.5. Genome-Wide Association Analysis and Candidate Genes Identification

Associations between genotypic and phenotypic data were analyzed using the kinship matrix in an MLM (Mixed linear model) by GAPIT (http://www.zzlab.net/GAPIT/) based on R 3.6.1 to control background variation and eliminate the spurious MTAs. Due to Bonferroni-Holm correction for multiple testing (α = 0.05) was too conserved and only a few significant MTAs were detected with this criterion, markers with an adjusted –log10 (p-value) ≥ 6.0 were regarded as the significant ones.
Candidate genes for the loci consistently identified in two or more environments were identified. The following steps were conducted to identify candidate genes for important QTL. Firstly, found all the genes located in LD block region around the peak SNP (± 150 Kb based on LD decay analysis) of each important QTL from the MSU Rice Genome Annotation Project (http://rice.plantbiology.msu.edu/cgi-bin/gbrowse/rice/). Then, all available SNPs located inside of these genes were searched. The genes (except for expressed protein, hypothetical protein, transposon protein and retrotransposon protein) with SNP in coding region and which could further lead to sense mutation were considered as candidate genes. Besides, mesocotyl elongation is highly regulated by various phytohormones, including strigolactones (SLs), cytokinin (CTK), brassinosteroid (BR), abscisic acid (ABA), jasmonates (JAs), gibberellins (GA) and auxin (IAA). Thus, these genes involved in the metabolism of the phytohormones talked above were selected as high-confidence candidate genes for mesocotyl elongation.

2.6. Gene Expression Analysis

Quantitative real-time PCR (qRT-PCR) was conducted to test the expression differences of candidate genes between accessions with extreme mesocotyl length (Table S4). The mesocotyl part was sampled at the 52 h after germination before the coleoptile excavated. Total RNA was extracted according to the Trizol method, cDNA was synthesized with the HiScript II 1st Strand cDNA Synthesis Kit (Vazyme, China). Primers (Table S5) were designed with the primer premier 5.0 software (http://www.premierbiosoft.com/; last accessed Jan 2019). The PCR procedure was conducted in a volume of 20 μL, containing 2μL cDNA, 0.4 μL of each primer (μM), 10 μL ChamQ Universal SYBR qPCR Master Mix. OsActin1 was used as the internal control to normalize the expression level of different samples. All assays were performed in two independent experiments, each with three repetitions.

3. Results

3.1. Phenotypic Variation of Mesocotyl Elongation

Mesocotyl length of 208 rice accessions were evaluated in water, sand and soil (2, 4 and 6 cm). Continuous variation was observed across all environments, presenting a wide range of mesocotyl length and indicating that this diversity panel was ideal for conducting GWAS (Figure 1). In sand culture, the mesocotyl lengths ranged from 0 to 5.03 cm, with an average of 1.67 cm; while in water, the mesocotyl lengths ranged from 0 to 3.57 cm, with an average of 1.21 cm. The seeds under soil were sown at a depth of 2, 4 and 6 cm, with the average mesocotyl lengths of 0.33 cm (0–0.71 cm), 0.51 cm (0–1.66 cm) and 2.47 cm (0–4.51 cm), respectively. Among different sowing depth, most mesocotyl lengths of 6 cm sowing depth were significantly longer than those of 2 cm or 4 cm sowing depth (Figure 1, Figure 2a). The maximum of mesocotyl lengths were 0.70 cm, 1.66 cm and 4.51 cm at 2 cm, 4 cm and 6 cm sowing depths in soil, respectively, indicating that mesocotyl elongation of these accessions had larger variation at 6 cm sowing depth than those at 2 cm or 4 cm. Furthermore, there was only a small proportion, about 32.4%, 7.2% and 70.0% accessions with longer mesocotyls (≥ 2.0 cm) in sand, water and soil (6 cm), respectively.
The correlation analysis showed that moderate correlations were observed among different culture media and high correlations were found among different sowing depth under soil. Mesocotyl length evaluated in sand showed active correlation with that under water (r = 0.70), soil (2 cm) (r = 0.74), soil (4 cm) (r = 0.75) and soil (6 cm) (r = 0.78). Also, mesocotyl length evaluated in water showed active correlation with that under soil (2 cm) (r = 0.76), soil (4 cm) (r = 0.78) and soil (6 cm) (r = 0.75). High correlations were observed among different sowing depth under soil. Mesocotyl length evaluated in soil (2 cm) showed significant active correlation with mesocotyl length at soil (4 cm) (r = 0.93) and soil (6 cm) (r = 0.85), whereas the correlation coefficient (r) between soil (4 cm) and soil (6 cm) was 0.86.
Analysis of variance (ANOVA) for mesocotyl length in 208 rice accessions revealed significant differences (P ≤ 0.001) among genotypes (G), environments (E), and genotype × environment interactions (GEI) (Table 1). The VE is much larger than VG and VGEI, indicating that mesocotyl length were greatly affected by environments. Besides, the broad sense heritability (hb2) estimate for mesocotyl length across all environments was 0.49. Thus, identifying the loci significantly associated with mesocotyl length under multiple environments is crucial for MAS breeding.

3.2. Marker Coverage

The sequencing data were mapped to the Nipponbare reference genome (IRGSP 1.0). After removing the SNPs with minor allele frequency (MAF) < 5% and missing data > 10%, 2,338,386 SNPs were left and employed for GWAS. The chromosome size varied from 22.8 Mb for chromosome 9 to 43.2 Mb for chromosome 1. These markers spanned a physical distance of 373 Mb, with an average density of 0.16 Kb per marker.

3.3. Population Structure and Linkage Disequilibrium

The diversity panel could be divided into two subgroups, the Japonica subpopulation and the other derived from the Indica accessions, whose characterization were largely consistent with geographic origins. Besides, admixture between Indica and Japonica were also observed in the present study (Figure 3c). Numerous studies have shown that the lack of appropriate correction for population structure can lead to spurious maker-trait associations (MTAs). Consequently, to eliminate spurious MTAs resulting from population structure, an MLM implemented in GAPIT were adopted for association analysis in the current study. PCA analysis indicated that the top three PCs could explain 28.5%, 8.2% and 3.2% of the total variation of population structure, respectively, and this panel consists of two subgroups (Figure 3a). The neighbor-joining (NJ) tree showing two clades also represented the two subpopulations (Figure 3b).
The decay of LD along physical distances was computed for the 208 rice accessions. A scatter r2 against physical distance showed a clean pattern of LD decay in the 208 rice accessions. A critical value of the determination coefficients r2 > 0.2 was determined to be the threshold for LD decay and the LD decay distance was about 150 Kb (Figure 3d).

3.4. GWAS of Mesocotyl Lengths

As Bonferroni correction was extremely conservative, a compromised threshold of –log10(P) ≥ 6.0 was used as the threshold for significantly associated SNPs. The Manhattan plots for the markers significantly associated with mesocotyl length under different conditions were shown in Figure 4. A total of 16 unique loci were detected, which explained the phenotypic variations ranged from 6.3% to 15.9%. Totally, 12 loci for mesocotyl length were identified on chromosomes 1, 2, 3, 5, 6, 8, 9 and 12 in sand, and each explained phenotypic variation ranged from 6.6 to 15.7%. Besides, seven loci were detected on chromosomes 2, 3, 5, 6 and 9 under water culture, which explain phenotypic variations ranging from 7.8 to 15.5%. There were remarkable differences among three sowing depths under soil culture. As far as the 2 cm sowing depth, a total of seven notable loci were detected on chromosomes 1, 4, 5, 7, 9 and 12, each explaining 7.6 to 15.7% of the phenotypic variation. When coming to the sowing depth of 4 and 6 cm, only two (chromosomes 1 and 12) and three (chromosomes 2, 6 and 12) loci were detected respectively, with an explanation rate of 14.4–14.8% and 6.3–15.9%.
Among all the above QTLs, 12 on chromosomes 1, 2, 3, 5 and 9 were identified under two or more environments. The loci located on chromosomes 1 (14.1–17.3 Mb) and 12 (13.3–15.0 Mb) were identified in sand, soil (2 cm) and soil (6 cm); whereas the locus on chromosome 12 (4.4-4.5 Mb) was detected in sand, soil (4 cm) and soil (6 cm). Besides, five loci on chromosomes 2 (5.6–8.7 Mb), 3 (15.2–15.3 Mb), 5 (5.8–6.3 Mb), 6 (7.3–9.7 Mb) and 9 (1.3–2.8 Mb) were identified in water and sand; two loci identified on chromosomes 5 (10.6–12.0 Mb) and 9 (7.1–9.1 Mb) were detected in sand and soil (6 cm); the locus on chromosome 6 (15.3–16.6 Mb) was identified on sand and soil (4 cm); the locus on chromosomes 2 (11.7–15.4 Mb) was identified on water and soil (4 cm) (Table 2, Table S2, Figure 4).

3.5. Candidate Genes Related with Mesocotyl Elongation

According to the results of GWAS, 107 candidate genes were found to be related with mesocotyl elongation (Table S3). Also, considering that mesocotyl elongation is highly regulated by various phytohormones, such as SLs, CTK, BR, ABA, JAs, GA and IAA, five genes involved in the biological metabolism of the phytohormones, cell elongation and division were selected as the high-confidence candidate genes for mesocotyl elongation (Table 3). That is the zinc finger CCCH type family protein (LOC_Os05g10670) and transcription factor jumonji (LOC_Os05g10770) on chromosome 5 (5.6–6.1 Mb), a gene encoding cytokinin-O-glucosyl transferase 2 on chromosome 9 (1.2–1.5 Mb) and two genes (4.3–4.5 Mb and 14.1–14.5 Mb) encoding flavin monooxygenase and 9-cis-epoxycarotenoid dioxygenase 1, respectively.
The expression of five candidate genes were detected using qRT-PCR (Figure 5). The results showed that there was significant difference between accessions with extreme mesocotyl length, with the longer one exhibited higher expression abundance. In detail, five genes showed more than 4.4-fold (LOC_Os05g10670), 2.8-fold (LOC_Os05g10770), 3.0-fold (LOC_Os09g03140), 8.0-fold (LOC_Os12g08780) and 32.6-fold (LOC_Os12g24800) higher expression in long-mesocotyl accessions compared to short-mesocotyl accessions, respectively.

4. Discussion

The GWAS panel, including released cultivars, advanced lines and landraces from different ecological regions, has a high genetic diversity with a wide range of mesocotyl length. In this study, Indica accessions had 32.4% and 70.0% accessions with longer mesocotyls (≥2.0 cm) in sand and 6 cm soil culture, respectively, but only 30.0% and 45.0% accessions in the Japonica, indicating that mesocotyl length had a great potential for improvement. Although mesocotyl is a crucial developmental trait of rice and is imminently required to improve crop adaptability to modern cultivation mode, the molecular mechanisms of its elongation and domestication remains unclear. Here, a GWAS based on re-sequencing were employed to identify associated SNP markers and understand the genetic basis of mesocotyl length.

4.1. Mesocotyl Elongation Varied among Culture Conditions

The elongation of mesocotyl was influenced by various environmental factors, such as light, temperature, moisture [28,29]. Thus, rice mesocotyl elongation were varied among different culture media, such as soil [26], sand [17], water [30]. In our study, mesocotyl length ranged from 0 to 3.99 cm under water culture, which was close to 0-5.05 cm in a previous study [17]. While, the mesocotyl length under 6 cm sand culture ranged from 0 to 4.88 cm, which differed from the previous result ranging from 0 to 2.05 cm under 5 cm sand culture [17]. Moreover, more than half of the 208 accessions had mesocotyl length longer than 1 cm under water culture. However, Wu et al. (2015) reported that only 10.5% accessions had mesocotyl length longer than 1 cm under water culture [17]. It is probably that the diverse results were mainly caused by different genotypes. Previous studies suggested that under varied sowing depths of soil cover, mesocotyl length was promoted by the depth in soil culture [11,18,26]. Lee et al. (1999) also reported that mesocotyl length under 3-cm sowing depth was shorter than that under 5-cm sowing depth [31]. Similarly, our results showed that mesocotyl length had a large variation among sowing depth, with an order of 6 cm > 4 cm >2 cm in soil. These results indicated that the variation of mesocotyl length existed in varied culture media and sowing depths.

4.2. LD Pattern

LD decay affects the precision of GWAS and was influenced by population structure, allele frequency, recombination rate and selection. Previous studies reported that LD decay in rice ranged between 120 Kb and 270 Kb [18,23,31,32]. In this panel, the LD decay was about 150 Kb for the whole genome (Figure 3e), consisting with previous reports [18,23,32,33].

4.3. GWAS Was Effective to Identify QTLs Controlled Complex Quantitative Traits

So far, several previous studies have reported the associated QTLs and candidate genes for rice mesocotyl elongation using linkage mapping or GWAS [20,25,34,35,36]. In this study, the loci identified on chromosome 7 (10.03–13.61 Mb) at the 6 cm sowing depth of soil was overlapped with QTLs qFML7-1, qFML7-2, qFML7-3, qFML7-4 and qFML7-5 [26]. The QTL on chromosome 1 (14.1–17.3 Mb) were overlapped with qml1-1 [37,38] and qML1 [23]. Furthermore, two QTLs on chromosome 2 (5.6–8.7 Mb) and 6 (7.3–9.7 Mb) identified under water and sand were overlapped with the QTL qML2 and qML6 detected from Shennong 265/Lijiangxintuanheigu RIL population [20]. Ouyang et al. (2015) have reported qml3-1, qml4-1 and qml5-1 by a RIL population [39], which nearly coincide with the QTL identified on chromosome 3 (15.2–15.3 Mb), 4 (8.3 Mb) and 5 (5.8–6.3 Mb) in our study. Zhang et al. (2006) detected four QTLs for submergence tolerance in seedling stage from a DH population, which explained phenotypic variations from 10.7% to 29.3% [35]. Among these, qLOE-8 were overlapped with the QTL regions located on chromosome 8 (208–5.2 Mb) identified in sand. Lee et al. (2012) reported two QTLs associated with mesocotyl elongation, qMel-1 (R2 = 37.3%) and qMel-3 (R2 = 6.5%) in weedy rice on chromosomes 1 and 3, respectively [40]. Moreover, Lu et al. (2016) reported 17 loci and 289 candidate genes for rice mesocotyl by GWAS on chromosomes 1 (4), 4 (4), 5 (1), 6 (3), 7 (1), 9 (3) and 11 (1), respectively [18]. QTLs identified on chromosomes 9 (1.3–2.8 Mb) and 12 (13.3–15.0 Mb) were overlapped with those QTLs reported by previous studies [18,40].
Wu et al. (2015) found that when under water culture and at a significant level of –log10(P) ≥8.0, a total of 13 loci were detected to be significantly associated with mesocotyl length, while only three significant SNPs were declared when the medium was merely changed to sand, and two of the three SNPs were co-localized under both conditions [17]. Lu et al. (2016) reported that six and seven loci were detected under two environments, respectively [18]. Only two of them on chromosome 6 were identified under both environments. In different studies, nearly 40 loci significantly associated with mesocotyl length have been detected on all twelve chromosomes by GWAS, and most of these loci in each study cannot be validated by others [17,18,25,26]. The interaction of G (genotype) × E (environment) can influence QTL and association mapping results, indicating that the effect of identified genomic regions was required to be estimated for each environment [41,42]. Only the major genomic regions could be co-localized under multi-environments.

4.4. Potential Candidate Genes for Mesocotyl Elongation were Identified

According to the result of GWAS, 107 candidate genes were found to be related to mesocotyl elongation. Among them, five involved in the biological metabolism of phytohormones, cell elongation and division were selected as the high-confidence candidate genes. Bioinformatics analysis indicated that loci on chromosome 5 (5.6–6.1 Mb) corresponding to zinc finger CCCH type family protein (LOC_Os05g10670) and transcription factor jumonji (LOC_Os05g10770). CCCH-type zinc finger proteins comprise a large family that is induced by drought, high temperature stress and hydrogen peroxide, and is also induced by abscisic acid, methyl jasmonate and salicylic acid [43,44]. Noh et al., (2004) reported that jumonji class transcription factor controls stem elongation, transposon activity and panicle development [45]. Besides, Feng et al. (2017) have identified several candidate genes for mesocotyl in rice including zinc-finger protein genes, which involved in the JA biosynthesis and signaling pathways [29]. Besides, a candidate gene on chromosome 9 (1.2–1.5 Mb) encoding cytokinin-O-glucosyltransferase 2 were identified. CK is a class of plant hormones that were first identified as cell division-promoting factors and were subsequently identified as factors that control various processes in plant growth and development, including mesocotyl elongation [46]. Besides, CK played an important role in the biosynthesis of BRs, a group of steroid plant hormones essential in plant growth and development [47]. Two candidate genes in loci 4.3–4.5 Mb and 14.1–14.5 Mb on chromosomes 12 were identified, which encodes flavin monooxygenase and 9-cis-epoxycarotenoid dioxygenase 1, respectively. Flavin monooxygenase catalyzes hydroxylation of the amino group of tryptamines, a rate-limiting step in tryptophan-dependent auxin biosynthesis, which regulate many processes in plant development [48]. The 9-cis-epoxycarotenoid dioxygenase is essential for the biosynthesis of ABA [49]. The 9-cis-epoxycarotenoid dioxygenase gene could lead to the over-production of abscisic acid, which plays an important role in mesocotyl elongation [50]. Expression of the five candidate genes in different accessions indicated that five genes were all functional in regulating mesocotyl elongation. The most significant difference in expression level was observed in gene, LOC_Os12g24800, which indicated that this gene could be the most likely functional gene for mesocotyl length. These results also proved that the loci identified by GWAS were reliable. Mesocotyl elongation is a consequence of complicated biological processes and its mechanism remains unclear; more detailed experimental analyses are needed to confirm the function of candidate genes in mesocotyl elongation.

5. Conclusions

In the present study, a GWAS for mesocotyl length in 208 rice accessions was conducted. Mesocotyl length varied in different culture conditions and sowing depths. GWAS have identified 16 loci significantly associated with mesocotyl elongation, which explained phenotypic variations ranged from 7.1% to 10.0%. Furthermore, five high-confidence candidate genes were identified on chromosomes 1, 3, 5, 9 and 12. The accessions with longer mesocotyl and the markers significantly associated with mesocotyl length can be used to promote the process of rice breeding. This study will facilitate our uanderstanding on the genetic architecture of rice mesocotyl length.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/1/49/s1, Table S1: Detailed information of 208 rice accessions, Table S2: SNPs significantly associated with mesocotyl length under different environments, Table S3: Candidate genes for mesocotyl length in 208 rice landraces or varieties, Table S4: Accessions used in qRT-PCR analysis, Table S5: Primers used in this study. The raw Illumina sequencing data will be available upon the request.

Author Contributions

Conceived and designed the experiments: H.L., Y.W., G.Y. Performed the experiments: H.L., J.Z., X.L. Analyzed the data: H.L., J.L. (Jiaolong Li), G.Y. Wrote and revised the paper: H.L., J.Z., J.L. (Jindong Liu), Y.W., G.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Agricultural Science and Technology Innovation Program Cooperation and Innovation Mission (Grant No. CAAS-XTCX2016001) and Shenzhen Science and Technology Projects (No. JSGG20160608160725473).

Acknowledgments

We thank Guosheng Xiong for providing some of the germplasms used in this study. We are grateful to the field research assistant Bingmei Lin in phenotyping and other field work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kumar, V.; Ladha, J.K. Direct seeding of rice: Recent developments and future research needs. Adv. Agron. 2011, 111, 297–413. [Google Scholar]
  2. Kato, Y.; Katsura, K. Rice adaptation to aerobic soils: Physiological considerations and implications for agronomy. Plant Prod. Sci. 2014, 17, 1–12. [Google Scholar] [CrossRef]
  3. Liu, H.; Hussain, S.; Zheng, M.; Peng, S.; Huang, J.; Cui, K.; Nie, L. Dry direct-seeded rice as an alternative to transplanted-flooded rice in Central China. Agron. Sustain. Dev. 2015, 35, 285–294. [Google Scholar] [CrossRef] [Green Version]
  4. Ohno, H.; Banayo, N.P.; Bueno, C.; Kashiwagi, J.I.; Nakashima, T.; Iwama, K.; Corales, A.M.; Garcia, R.; Kato, Y. On-farm assessment of a new early-maturing drought-tolerant rice cultivar for dry direct seeding in rainfed lowlands. Field Crop. Res. 2018, 219, 222–228. [Google Scholar] [CrossRef]
  5. GRiSP. Rice Almanac, 4th ed.; International Rice Research Institute: Manila, Philippines, 2013; p. 283. [Google Scholar]
  6. Haefele, S.M.; Kato, Y.; Singh, S. Climate ready rice: Augmenting drought tolerance with best management practices. Field Crop. Res. 2013, 190, 60–69. [Google Scholar] [CrossRef]
  7. Kono, M. Physiological aspects of lodging. In Science of the Rice Plant, Physiology, 2; Matsuo, Y., Kumazawa, K., Ishii, R., Ishihara, K., Hirata, H., Eds.; Food and Agricultural Policy Research Center: Tokyo, Japan, 1995; pp. 971–982. [Google Scholar]
  8. Ramzan, M. Evaluation of various planting methods in rice-wheat cropping system, Punjab, Pakistan. Rice Crop Rep. 2003, 4, 4–5. [Google Scholar]
  9. Kashiwagi, T.; Sasaki, H.; Ishimaru, K. Factors responsible for decreasing sturdiness of the lower part in lodging of rice (Oryza sativa L.). Plant Prod. Sci. 2005, 8, 166–172. [Google Scholar] [CrossRef]
  10. Mahender, A.; Anandan, A.; Pradhan, S.K. Early seedling vigour, an imperative trait for direct-seeded rice: An overview on physio-morphological parameters and molecular markers. Planta 2015, 241, 1027–1050. [Google Scholar] [CrossRef]
  11. Lee, H.S.; Sasaki, K.; Kang, J.W.; Sato, T.; Song, W.Y.; Ahn, S.N. Mesocotyl elongation is essential for seedling emergence under deep-seeding condition in rice. Rice 2017, 10, 32. [Google Scholar] [CrossRef] [Green Version]
  12. Ekanayake, I.J.; Garrity, D.P.; O’Toole, J.C. Influence of Deep Root Density on Root Pulling Resistance in Rice. Crop Sci. 1986, 26, 1181–1186. [Google Scholar] [CrossRef]
  13. Kirby, E. Effect of sowing depth on seedling emergence, growth and development in barley and wheat. Field Crop. Res. 1993, 35, 101–111. [Google Scholar] [CrossRef]
  14. Turner, F.T.; Chen, C.C.; Bollich, C.N. Coleoptile and mesocotyl lengths in semi-dwarf rice seedlings. Crop Sci. 1982, 22, 43–46. [Google Scholar] [CrossRef]
  15. Chung, N.J. Elongation habit of mesocotyls and coleoptiles in weedy rice with high emergence ability in direct-seeding on dry paddy fields. Crop Pasture Sci. 2010, 61, 911–917. [Google Scholar] [CrossRef]
  16. Trachsel, S.R.; Messmer, P.; Stamp, N.R.; Hund, A. QTLs for early vigor of tropical maize. Mol. Breed. 2010, 25, 91–103. [Google Scholar] [CrossRef] [Green Version]
  17. Wu, J.; Feng, F.; Lian, X.; Teng, X.; Wei, H.; Yu, H.; Xie, W.; Yan, M.; Fan, P.; Li, Y.; et al. Genome-wide Association Study (GWAS) of mesocotyl elongation based on re-sequencing approach in rice. BMC Plant Biol. 2015, 15, 218. [Google Scholar] [CrossRef] [Green Version]
  18. Lu, Q.; Zhang, M.; Niu, X.; Wang, C.; Xu, Q.; Feng, Y.; Wang, S.; Yuan, X.; Yu, H.; Wang, Y. Uncovering novel loci for mesocotyl elongation and shoot length in indica rice through genome-wide association mapping. Planta 2016, 243, 645–657. [Google Scholar] [CrossRef] [Green Version]
  19. Mgonja, M.; Dilday, R.; Skinner, S.; Collins, F. Association of mesocotyl elongation with seedling vigor in rice. J. Ark. Acad. Sci. 1988, 42, 52–55. [Google Scholar]
  20. Huang, C.; Jiang, S.K.; Feng, L.; Xu, Z.J.; Chen, W.F. QTL analysis for mesocotyl length in rice (Oryza sativa L.). Acta Agron. Sin. 2010, 36, 1108–1113. [Google Scholar]
  21. Huang, X.; Zhao, Y.; Li, C.; Wang, A.; Zhao, Q.; Li, W.; Guo, Y.; Deng, L.; Zhu, C.; Fan, D.; et al. Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat. Genet. 2012, 44, 32–39. [Google Scholar] [CrossRef]
  22. Chen, W.; Gao, Y.; Xie, W.; Gong, L.; Lu, K.; Wang, W.; Li, Y.; Liu, X.; Zhang, H.; Dong, H.; et al. Genome-wide association analyses provide genetic and biochemical insights into natural variation in rice metabolism. Nat. Genet. 2014, 46, 714–721. [Google Scholar] [CrossRef]
  23. Huang, X.; Sang, T.; Zhao, Q.; Feng, Q.; Zhao, Y.; Li, C.; Zhu, C.; Lu, T.; Zhang, Z.; Li, M.; et al. Genome-wide association studies of 14 agronomic traits in rice landraces. Nat. Genet. 2010, 42, 961–967. [Google Scholar] [CrossRef] [PubMed]
  24. Zhao, K.; Tung, C.W.; Eizenga, G.C.; Wright, M.H.; Ali, M.L.; Price, A.H.; Norton, G.J.; Islam, M.R.; Reynolds, A.; Mezey, J.; et al. Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nat. Commun. 2011, 2, 467. [Google Scholar] [CrossRef] [PubMed]
  25. Sun, S.; Wang, T.; Wang, L.; Li, X.; Jia, Y.; Liu, C.; Huang, X.; Xie, W.; Wang, X. Natural selection of a GSK3 determines rice mesocotyl domestication by coordinating strigolactone and brassinosteroid signaling. Nat. Commun. 2018, 9, 2523. [Google Scholar] [CrossRef] [PubMed]
  26. Zhao, Y.; Zhao, W.; Jiang, C.; Wang, X.; Xiong, H.; Todorovska, E.G.; Yin, Z.; Chen, Y.; Wang, X.; Xie, J.; et al. Genetic architecture and candidate genes for deep-sowing tolerance in rice revealed by non-syn GWAS. Front Plant Sci. 2018, 9, 332. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Alexander, D.H.; Novembre, J.; Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res 2009, 19, 1655–1664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Radford, B.J.; Henzell, R.G. Temperature affects the mesocotyl and coleoptile length of grain sorghum genotypes. Aust. J. Agric. Res. 1990, 41, 79–87. [Google Scholar] [CrossRef]
  29. Feng, F.; Mei, H.; Fan, P.; Li, Y.; Xu, X.; Wei, H.; Yan, M.; Luo, L. Dynamic transcriptome and phytohormone profiling along the time of light exposure in the mesocotyl of rice seedling. Sci. Rep. 2017, 7, 11961. [Google Scholar] [CrossRef] [Green Version]
  30. Cai, H.; Morishima, H. QTL clusters reflect character associations in wild and cultivated rice. Theor. Appl. Genet. 2002, 104, 1217–1228. [Google Scholar] [CrossRef]
  31. Lee, S.; Kim, J.; Hong, S. Effects of priming and growth regulator treatment of seed on emergence and seedling growth of rice. Korean J. Crop Sci. 1999, 44, 134–137. [Google Scholar]
  32. Qiu, X.; Pang, Y.; Yuan, Z.; Xing, D.; Xu, J.; Dingkuhn, M.; Li, Z.; Ye, G. Genome-wide association study of grain appearance and milling quality in a worldwide collection of indica rice germplasm. PLoS ONE 2015, 10, e0145577. [Google Scholar] [CrossRef] [Green Version]
  33. McCouch, S.R.; Wright, M.H.; Tung, C.W.; Maron, L.G.; McNally, K.L.; Fitzgerald, M.; Singh, N.; DeClerck, G.; Agosto-Perez, F.; Korniliev, P.; et al. Open access resources for genome-wide association mapping in rice. Nat. Commun. 2016, 7, 10532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Redona, E.D.; Mackill, D.J. Mapping quantitative trait loci for seedling vigor in rice using RFLPs. Theor. Appl. Genet. 1996, 92, 395–402. [Google Scholar] [CrossRef] [PubMed]
  35. Zhang, G.; Zeng, D.; Hu, S.; Su, Y.; Ajia, L.; Guo, L.; Qian, Q. QTL analysis of traits concerned submergence tolerance at seedling stage in rice (Oryza sativa L.). Acta Agron. Sin. 2006, 32, 1280–1286. [Google Scholar]
  36. Lim, J.H.; Yang, H.J.; Jung, K.H.; Yoo, S.C.; Paek, N.C. Identification of chromosomal regions influencing mesocotyl elongation by bulked segregation analysis based on genome re-sequencing in rice. Acta Agric. Shanghai 2017, 33, 10–15. [Google Scholar]
  37. Cao, L.; Zhu, J.; Ren, L.; Zhao, S.; Yan, Q. Mapping QTLs and epistasis for seedling vigor in rice (Oryza sativa L.). Acta Agron. Sin. 2002, 28, 809–815. [Google Scholar]
  38. Cao, L.; Zhu, J.; Yan, Q.; He, L.; Wei, X.; Cheng, S. Mapping QTLs with epistasis for mesocotyl length in a DH population from indica-japonica cross of rice (Oryza sativa). Chin. J. Rice Sci. 2002, 16, 221–224. [Google Scholar]
  39. Ouyang, Y.N.; Zhang, Q.Y.; Zhang, K.Q.; Yu, S.M.; Zhuang, J.Y.; Jin, Q.Y.; Cheng, S.H. QTL mapping and interaction analysis of genotype × environment (Fe2+-concentrations) for mesocotyl length in rice (Oryza sativa L.). Acta Genet. Sin. 2005, 32, 712–718. [Google Scholar]
  40. Lee, H.S.; Sasaki, K.; Higashitani, A.; Ahn, S.N.; Sato, T. Mapping and characterization of quantitative trait loci for mesocotyl elongation in rice (Oryza sativa L.). Rice 2012, 5, 1–10. [Google Scholar] [CrossRef] [Green Version]
  41. Zhang, K.; Tian, J.; Zhao, L.; Wang, S. Mapping QTLs with epistatic effects and QTL x environment interactions for plant height using a doubled haploid population in cultivated wheat. J Genet Genom. 2008, 35, 119–127. [Google Scholar] [CrossRef]
  42. Tetard-Jones, C.; Kertesz, M.A.; Preziosi, R.F. Identification of plant quantitative trait loci modulating a rhizobacteria-aphid indirect effect. PLoS ONE 2012, 7, e41524. [Google Scholar] [CrossRef]
  43. Jan, A.; Maruyama, K.; Todaka, D.; Kidokoro, S.; Abo, M.; Yoshimura, E.; Shinozaki, K.; Nakashima, K.; Yamaguchi-Shinozaki, K. OsTZF1, a CCCH-tandem zinc finger protein, confers delayed senescence and stress tolerance in rice by regulating stress-related genes. Plant Physiol. 2013, 161, 1202–1216. [Google Scholar] [CrossRef] [Green Version]
  44. Xiong, Q.; Ma, B.; Lu, X.; Huang, Y.H.; He, S.J.; Yang, C.; Yin, C.C.; Zhao, H.; Zhou, Y.; Zhang, W.K.; et al. Ethylene-Inhibited Jasmonic Acid Biosynthesis Promotes Mesocotyl/Coleoptile Elongation of Etiolated Rice Seedlings. Plant Cell 2017, 29, 1053–1072. [Google Scholar] [CrossRef] [Green Version]
  45. Noh, B.; Lee, S.H.; Kim, H.J.; Yi, G.; Shin, E.A.; Lee, M.; Jung, K.J.; Doyle, M.R.; Amasino, R.M.; Noh, Y.S. Divergent roles of a pair of homologous Jumonji/Zinc-Finger–class transcription factor proteins in the regulation of Arabidopsis flowering time. Plant Cell 2004, 16, 2601–2613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Hu, Z.; Yamauchi, T.; Yang, J.; Jikumaru, Y.; Tsuchida-Mayama, T.; Ichikawa, H.; Takamure, I.; Nagamura, Y.; Tsutsumi, N.; Yamaguchi, S.; et al. Strigolactone and cytokinin act antagonistically in regulating rice mesocotyl elongation in darkness. Plant Cell Physiol. 2014, 55, 30–41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Yuldashev, R.; Avalbaev, A.; Bezrukova, M.; Vysotskaya, L.; Khripach, V.; Shakirova, F. Cytokinin oxidase is involved in the regulation of cytokinin content by 24-epibrassinolide in wheat seedlings. Plant Physiol. Biochem. 2012, 55, 1–6. [Google Scholar] [CrossRef] [PubMed]
  48. Zhao, Y.; Christensen, S.K.; Fankhauser, C.; Cashman, J.R.; Cohen, J.D.; Weigel, D.; Chory, J. A role for flavin monooxygenase-like enzymes in auxin biosynthesis. Science 2001, 291, 306–309. [Google Scholar] [CrossRef] [PubMed]
  49. Wan, X.R.; Li, L. Regulation of ABA level and water-stress tolerance of Arabidopsis by ectopic expression of a peanut 9-cis-epoxycarotenoid dioxygenase gene. Biochem. Biophys. Res. Commun. 2006, 347, 1030–1038. [Google Scholar] [CrossRef]
  50. Thompson, A.J.; Jackson, A.C.; Symonds, R.C.; Mulholland, B.J.; Dadswell, A.R.; Blake, P.S.; Burbidge, A.; Taylor, I.B. Ectopic expression of a tomato 9-cis-epoxycarotenoid dioxygenase gene causes over-production of abscisic acid. Plant J. 2000, 23, 363–374. [Google Scholar] [CrossRef]
Genes 11 00049 g001 550
Figure 1. Summary of the phenotypic variations of mesocotyl length among different culture conditions, including correlations between mesocotyl lengths under different culture conditions and mesocotyl length distribution under varied culture conditions. Soil_2: soil culture with 2 cm sowing depth; Soil_4: soil culture with 4 cm sowing depth; Soil_6: soil culture with 6 cm sowing depth, numbers in blue boxes are correlation coefficients (r), * indicates the significance at p ≤ 0.05.
Figure 1. Summary of the phenotypic variations of mesocotyl length among different culture conditions, including correlations between mesocotyl lengths under different culture conditions and mesocotyl length distribution under varied culture conditions. Soil_2: soil culture with 2 cm sowing depth; Soil_4: soil culture with 4 cm sowing depth; Soil_6: soil culture with 6 cm sowing depth, numbers in blue boxes are correlation coefficients (r), * indicates the significance at p ≤ 0.05.
Genes 11 00049 g001
Genes 11 00049 g002 550
Figure 2. Phenotypic variations of mesocotyl elongation among three culture media. (a) Mesocotyl length of 208 rice accessions measured in soil culture with three sowing depths. (b) Mesocotyl length variation in sand, soil culture with 6 cm sowing depth and water.
Figure 2. Phenotypic variations of mesocotyl elongation among three culture media. (a) Mesocotyl length of 208 rice accessions measured in soil culture with three sowing depths. (b) Mesocotyl length variation in sand, soil culture with 6 cm sowing depth and water.
Genes 11 00049 g002
Genes 11 00049 g003 550
Figure 3. Population structure and LD decay of 208 rice accessions. (a) Principal component plot, different colors represent different ingredients; (b) neighbor-Jointing tree; (c) structure; (d) LD decay analysis.
Figure 3. Population structure and LD decay of 208 rice accessions. (a) Principal component plot, different colors represent different ingredients; (b) neighbor-Jointing tree; (c) structure; (d) LD decay analysis.
Genes 11 00049 g003
Genes 11 00049 g004 550
Figure 4. Genome-wide association studies for mesocotyl length among three culture media. Manhattan plots and Quantile-quantile plots of mesocotyl lengths in sand culture (a), water (b) and soil culture with a sowing depths of 2 cm (c), 4 cm (d) and 6 cm (e), respectively. Negative log10-transformed p values from a genome-wide scan are plotted against position on each of 12 chromosomes. The red horizontal dashed lines indicate the genome-wide significance threshold. Blue stripes show important signals in different GWAS.
Figure 4. Genome-wide association studies for mesocotyl length among three culture media. Manhattan plots and Quantile-quantile plots of mesocotyl lengths in sand culture (a), water (b) and soil culture with a sowing depths of 2 cm (c), 4 cm (d) and 6 cm (e), respectively. Negative log10-transformed p values from a genome-wide scan are plotted against position on each of 12 chromosomes. The red horizontal dashed lines indicate the genome-wide significance threshold. Blue stripes show important signals in different GWAS.
Genes 11 00049 g004
Genes 11 00049 g005 550
Figure 5. Expression analysis of five candidate genes. Relative expressions of the five genes in long- and short -mesocotyl length rice accessions were tested by qRT-PCR. Transcription levels were normalized to that of Osactin1 and compared to the short control. 1 to 5 represent five candidate genes, that is 1, LOC_Os05g10670; 2, LOC_Os05g10770; 3, LOC_Os09g03140; 4, LOC_Os12g08780; 5, LOC_Os12g24800. Values of fold change are showed in mean ± SEM. Relative transcription levels were calculated by the log2-△△CT method. * and ** indicate the significance at p ≤ 0.05 and p ≤ 0.01, respectively.
Figure 5. Expression analysis of five candidate genes. Relative expressions of the five genes in long- and short -mesocotyl length rice accessions were tested by qRT-PCR. Transcription levels were normalized to that of Osactin1 and compared to the short control. 1 to 5 represent five candidate genes, that is 1, LOC_Os05g10670; 2, LOC_Os05g10770; 3, LOC_Os09g03140; 4, LOC_Os12g08780; 5, LOC_Os12g24800. Values of fold change are showed in mean ± SEM. Relative transcription levels were calculated by the log2-△△CT method. * and ** indicate the significance at p ≤ 0.05 and p ≤ 0.01, respectively.
Genes 11 00049 g005
Table
Table 1. ANOVA for mesocotyl length under different environment.
Table 1. ANOVA for mesocotyl length under different environment.
SourcedfSum of SquareMean SquareF-Value
Block (Environment)50.90630.18132.5698 ***
Genotype207583.15712.817239.9407 ***
Environment41270.5662317.64154503.3657 ***
Genotyped × Environment813533.98770.65689.3120 ***
*** indicates the significance at p ≤ 0.001.
Table
Table 2. Loci significantly associated with mesocotyl length under different environment.
Table 2. Loci significantly associated with mesocotyl length under different environment.
Peak MarkerChr.Genomic Regions (Mb)P ValueR2 (%)EnvironmentReported LocusReferences
rs_1_16752644114.1~17.31.9 E-08~7.7 E-077.2–15.2Sand, Soil (2 cm), Soil (6 cm)qml1-1 and qML1Cao et al. 2002a,
Huang et al., 2010
rs_2_564711625.6-8.74.8~9.4 E-077.6–8.1Water, SandqML2Huang et al., 2010
rs_2_11740641211.7-15.42.8-7.3 E-078.1–15.5Water, Soil (4 cm)
rs_3_15245138315.2-15.32.5-5.7 E-0714.8–14.9Water, Sandqml3-1Ouyang et al., 2005
rs_4_723428447.28.2 E-0814.9Soil (6 cm)Qml4-1Ouyang et al., 2005
rs_5_620613655.8-6.39.5E-09~9.7E-076.6–15.3Water, SandQml5-5Ouyang et al. 2005
rs_5_11950679510.6-12.06.6~7.6E-076.6–8.1Sand, Soil (6 cm)
rs_6_731300567.3-9.77.7E-08~8.4E-077.1–8.1Water, SandqML6Huang et al., 2010
rs_6_16643928615.3-16.61.2~3.7 E-076.7–14.4Sand, Soil (4 cm)
rs_7_10769188710.0-13.62.8E-09~9.9E-076.4–15.9Soil (6 cm)qFML7-1Zhao et al., 2018
rs_8_521944982.8-5.26.8E-08~3.5E-0714.6–15.3SandqLOE-8Zhang et al., 2006
rs_9_144682491.3-2.84.8E-08~4.8E-0714.9–15.2Water, Sandseq-rs4080Lu et al., 2016
rs_9_705924697.1-9.18.5 E-08~1.2 E-0714.7–14.8Sand, Soil (6 cm)
rs_12_4444693124.4-4.54.3E-09~5.9E-0714.8–15.7Sand, Soil (4 cm), Soil (6 cm)
rs_12_6008455126.0-7.81.2~4.6 E-0714.7–15.4Sand
rs_12_141554861213.3-15.02.3E-08~9.9E-077.5–15.3Sand, Soil (2 cm), Soil (6 cm)qMel-12Lee et al. 2012b
Table
Table 3. High confidence candidate genes for mesocotyl lengths in 208 rice landraces or varieties.
Table 3. High confidence candidate genes for mesocotyl lengths in 208 rice landraces or varieties.
ChromosomeGenomic RegionsCandidate GeneFunctional Annotation
55688534-6067598LOC_Os05g10670Zinc finger CCCH type family protein, putative, expressed
55688534-6067598LOC_Os05g10770Transcription factor jumonji, putative, expressed
91245230-1521824LOC_Os09g03140Cytokinin-O-glucosyltransferase 2, putative, expressed
124369693-4551492LOC_Os12g08780Flavin monooxygenase, putative, expressed
1214080486-14520450LOC_Os12g248009-cis-epoxycarotenoid dioxygenase 1, chloroplast precursor, putative, expressed

Share and Cite

MDPI and ACS Style

Liu, H.; Zhan, J.; Li, J.; Lu, X.; Liu, J.; Wang, Y.; Zhao, Q.; Ye, G. Genome-Wide Association Study (GWAS) for Mesocotyl Elongation in Rice (Oryza sativa L.) under Multiple Culture Conditions. Genes 2020, 11, 49. https://doi.org/10.3390/genes11010049

AMA Style

Liu H, Zhan J, Li J, Lu X, Liu J, Wang Y, Zhao Q, Ye G. Genome-Wide Association Study (GWAS) for Mesocotyl Elongation in Rice (Oryza sativa L.) under Multiple Culture Conditions. Genes. 2020; 11(1):49. https://doi.org/10.3390/genes11010049

Chicago/Turabian Style

Liu, Hongyan, Junhui Zhan, Jiaolong Li, Xiang Lu, Jindong Liu, Yamei Wang, Quanzhi Zhao, and Guoyou Ye. 2020. "Genome-Wide Association Study (GWAS) for Mesocotyl Elongation in Rice (Oryza sativa L.) under Multiple Culture Conditions" Genes 11, no. 1: 49. https://doi.org/10.3390/genes11010049

APA Style

Liu, H., Zhan, J., Li, J., Lu, X., Liu, J., Wang, Y., Zhao, Q., & Ye, G. (2020). Genome-Wide Association Study (GWAS) for Mesocotyl Elongation in Rice (Oryza sativa L.) under Multiple Culture Conditions. Genes, 11(1), 49. https://doi.org/10.3390/genes11010049

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop