Fine mapping and candidate gene analysis of a major QTL for panicle structure in rice

Key message

A gene not only control tiller and plant height, but also regulate panicle structure by QTL dissection in rice.

Abstract

An ideal panicle structure is important for improvement of plant architecture and rice yield. In this study, using recombinant inbred lines (RILs) of PA64s and 93-11, we identified a quantitative trait locus (QTL), designated qPPB3 for primary panicle branch number. With a BC₃F₂ population derived from a backcross between a resequenced RIL carrying PA64s allele and 93-11, qPPB3 was fine mapped to a 34.6-kb genomic region. Gene prediction analysis identified four putative genes, among which Os03g0203200, a previously reported gene for plant height and tiller number, Dwarf 88 (D88)/Dwarf 14 (D14), had three nucleotide substitutions in 93-11 compared with PA64s. The T to G substitution resulted in one amino acid change from valine in 93-11 to glycine in PA64s. Real-time PCR analysis showed expression level of D88 was higher in 93-11 than PA64s. The expression of APO1 and IPA1 increased, while GN1a and DST decreased in 93-11 compared with PA64s. Therefore, D88/D14 is not only a key regulator for branching, but also affects panicle structure.     

OsGA20ox1, a candidate gene for a major QTL controlling seedling vigor in rice

Seedling vigor is among the major determinants of stable stand establishment in direct-seeded rice (Oryza sativa L.) in temperate regions. Quantitative trait loci (QTL) for seedling vigor were identified using 250 recombinant inbred lines (RILs) derived from a cross between two japonica rice cultivars Kakehashi and Dunghan Shali. Seedling heights measured at 14 days after sowing were 20.3 and 29.4 cm for Kakehashi and Dunghan Shali, respectively. For the RILs, the height ranged from 14.1 to 31.7 cm. Four putative QTLs associated with seedling height were detected. qPHS3-2, the major QTL that was located on the long arm of chromosome 3, accounted for 26.2 % of the phenotypic variance. Using progeny of the near isogenic lines (NILs) produced by the backcross introduction of a chromosome segment carrying this major QTL into an elite cultivar Iwatekko, we fine-mapped qPHS3-2 to a 81-kb interval between two markers, ID_CAPS_01 and RM16227. Within this mapped region, we identified the gene OsGA20ox1, which is related to gibberellin (GA) biosynthesis. The relative expression levels of GA20ox1 in seedlings of Dunghan Shali and NILs were higher than that of Iwatekko. Concomitantly, the amount of endogenous active GA was higher in Dunghan Shali and the NILs compared to the level detected in Iwatekko. These results indicate that OsGA20ox1 is a strong candidate gene for major QTL controlling seedling vigor in rice.     

Fine mapping and identification of candidate rice genes associated with qSTV11 SG , a major QTL for rice stripe disease resistance

Rice stripe disease, caused by rice stripe virus (RSV) is a serious constraint to rice production in subtropical regions of East Asia. We performed fine mapping of a RSV resistance QTL on chromosome 11, qSTV11 ( SG ), using near-isogenic lines (NILs, BC(6)F(4)) derived from a cross between the highly resistant variety, Shingwang, and the highly susceptible variety, Ilpum, using 11 insertion and deletion (InDel) markers. qSTV11 ( SG ) was localized to a 150-kb region between InDel 11 (17.86 Mbp) and InDel 5 (18.01?Mbp). Among the two markers in this region, InDel 7 is diagnostic of RSV resistance in 55 Korean japonica and indica rice varieties. InDel 7 could also distinguish the allele type of Nagdong, Shingwang, Mudgo, and Pe-bi-hun from Zenith harboring the Stv-b ( i ) allele. As a result, qSTV11 ( SG ) is likely to be the Stv-b ( i ) allele. There were 21 genes in the 150-kb region harboring the qSTV11 ( SG ) locus. Three of these genes, LOC_Os11g31430, LOC_Os11g31450, and LOC_Os11g31470, were exclusively expressed in the susceptible variety. These expression profiles were consistent with the quantitative nature along with incomplete dominance of RSV resistance. Sequencing of these genes showed that there were several amino acid substitutions between susceptible and resistant varieties. Putative functions of these candidate genes for qSTV11 ( SG ) are discussed.     

Saturation mapping of QTL regions and identification of putative candidate genes for drought tolerance in rice

We have developed 85 new markers (50 RFLPs, 5 SSRs, 12 DD cDNAs, 9 ESTs, 8 HSP-encoding cDNAs and one BSA-derived AFLP marker) for saturation mapping of QTL regions for drought tolerance in rice, in our efforts to identify putative candidate genes. Thirteen of the markers were localized in the close vicinity of the targeted QTL regions. Fifteen of the additional markers mapped, respectively, inside one QTL region controlling osmotic adjustment on chromosome 3 ( oa3.1) and 14 regions that affect root traits on chromosomes 1, 2, 4, 5, 6, 7, 8, 9, 10 and 12. Differential display was used to identify more putative candidate genes and to saturate the QTL regions of the genetic map. Eleven of the isolated cDNA clones were found to be derived from drought-inducible genes. Two of them were unique and did not match any genes in the GenBank, while nine were highly similar to cDNAs encoding known proteins, including a DnaJ-related protein, a zinc-finger protein, a protease inhibitor, a glutathione-S-transferase, a DNA recombinase, and a protease. Twelve new cDNA fragments were mapped onto the genetic linkage map; seven of these mapped inside, or in close proximity to, the targeted QTL regions determining root thickness and osmotic adjustment capacity. The gene I12A1, which codes for a UDP-glucose 4-epimerase homolog, was identified as a putative target gene within the prt7.1/brt7.1 QTL region, as it is involved in the cell wall biogenesis pathway and hence may be implicated in modulating the ability of rice roots to penetrate further into the substratum when exposed to drought conditions. RNAs encoding elongation factor 1, a DnaJ-related protein, and a homolog of wheat zinc-finger protein were more prominently induced in the leaves of IR62266 (the lowland rice parent of the mapping materials used) than in those of CT9993 (the upland rice parent) under drought conditions. Homologs of 18S ribosomal RNA, and mRNAs for a multiple-stress induced zinc-finger protein, a protease inhibitor, and a glutathione-S-transferase were expressed at significantly higher levels in CT9993 than in IR62266. Thus several genes involved in the regulation of DNA structure and mRNA translation were found to be drought-regulated, and may be implicated in drought resistance.Communicated by R. Hagemann     

QTL analysis and mapping of pi21, a recessive gene for field resistance to rice blast in Japanese upland rice

Field resistance is defined as the resistance that allows effective control of a parasite under natural field conditions and is durable when exposed to new races of that parasite. To identify the genes for field resistance to rice blast, quantitative trait loci (QTLs) conferring field resistance to rice blast in Japanese upland rice were detected and mapped using RFLP and SSR markers. QTL analysis was carried out in F₄ progeny lines from the cross between Nipponbare (moderately susceptible, lowland) and Owarihatamochi (resistant, upland). Two QTLs were detected on chromosome 4 and one QTL was detected on each of chromosomes 9 and 12. The phenotypic variation explained by each QTL ranged from 7.9 to 45.7% and the four QTLs explained 66.3% of the total phenotypic variation. Backcrossed progeny lines were developed to transfer the QTL with largest effect using the susceptible cultivar Aichiasahi as a recurrent parent. Among 82 F₃ lines derived from the backcross, resistance segregated in the expected ratio of resistant 1 : heterozygous 2 : susceptible 1. The average score for blast resistance measured in the field was 4.2 ± 0.67, 7.5 ± 0.51and 8.2 ± 0.66, for resistant, heterozygous and susceptible groups, respectively. The resistance gene, designated pi21, was mapped on chromosome 4 as a single recessive gene between RFLP marker loci G271 and G317 at a distance of 5.0 cM and 8.5 cM, respectively. The relationship to previously reported major genes and QTLs conferring resistance to blasts, and the significance of marker-assisted selection to improve field resistance, are discussed. Received: 8 June 2000 / Accepted: 24 November 2000     

DFP: a Bioconductor package for fuzzy profile identification and gene reduction of microarray data

Background  

Expression profiling assays done by using DNA microarray technology generate enormous data sets that are not amenable to simple analysis. The greatest challenge in maximizing the use of this huge amount of data is to develop algorithms to interpret and interconnect results from different genes under different conditions. In this context, fuzzy logic can provide a systematic and unbiased way to both (i) find biologically significant insights relating to meaningful genes, thereby removing the need for expert knowledge in preliminary steps of microarray data analyses and (ii) reduce the cost and complexity of later applied machine learning techniques being able to achieve interpretable models.     

Fine mapping a major QTL for flag leaf size and yield-related traits in rice

Leaf size is a major determinant of plant architecture and yield potential in crops. A previous study showed that the genomic region of chromosome 1 contains a major quantitative trait locus (QTL) for flag leaf size in a set of backcross recombinant inbred lines derived from two elite parental lines (Zhenshan 97 and 93-11). In the present study, the QTL (qFL1) was shown to explain a large proportion of the variation in flag leaf size (leaf length, width and area) in derived populations (BC₂F₃ and BC₃F₂) in multiple environments. Using a large segregating population, we narrowed the location of qFL1 to a 31 kb region containing four predicted genes. Expression of one of these genes, OsFTL1, differed between leaves in near-isogenic lines carrying alleles of Zhenshan 97 and 93-11. qFL1 had a pleiotropic effect on flag leaf size and yield-related traits. Conditional QTL analysis of the derived population (BC₃F₂) supports the assertion that qFL1 is the QTL for flag leaf length and exhibits pleiotropy. Pyramiding of qFL1 with two known genes (GS3 and Wx) from 93-11 into Zhenshan 97 enlarged flag leaves, improved grain size and amylose content, and increased yield per plant, but slightly delayed heading date. These results provide a foundation for the functional characterization of the gene underlying the pleiotropic effects of qFL1 and for genetic improvement of the plant architecture and yield potential of rice.     

QTL analysis for rice grain length and fine mapping of an identified QTL with stable and major effects

Grain length in rice plays an important role in determining rice appearance, milling, cooking and eating quality. In this study, the genetic basis of grain length was dissected into six main-effect quantitative trait loci (QTLs) and twelve pairs of epistatic QTLs. The stability of these QTLs was evaluated in four environments using an F₇ recombinant inbred line (RIL) population derived from the cross between a Japonica variety, Asominori, and an Indica variety, IR24. Moreover, chromosome segment substitution lines (CSSLs) harboring each of the six main-effect QTLs were used to evaluate gene action of QTLs across eight environments. A major QTL denoted as qGL-3a, was found to express stably not only in the isogenic background of Asominori but also in the recombinant background of Asominori and IR24 under multiple environments. The IR24 allele at qGL-3a has a positive effect on grain length. Based on the test of advanced backcross progenies, qGL-3a was dissected as a single Mendelian factor, i.e., long rice grain was controlled by a recessive gene gl-3. High-resolution genetic and physical maps were further constructed for fine mapping gl-3 by using 11 simple sequence repeat (SSR) markers designed using sequence information from seven BAC/PAC clones and a BC₄F₂ population consisting of 2,068 individuals. Consequently, the gl-3 gene was narrowed down to a candidate genomic region of 87.5 kb long defined by SSR markers RMw357 and RMw353 on chromosome 3, which provides a basis for map-based cloning of this gene and for marker-aided QTL pyramiding in rice quality breeding.     

Efficacy of RNA amplification is dependent on sequence characteristics: implications for gene expression profiling using a cDNA microarray

Minute tissue samples or single cells increasingly provide the starting material for gene expression profiling, which often requires RNA amplification. Although much effort has been put into optimizing amplification protocols, the relative abundance of RNA templates in the amplified product is frequently biased. We applied a T7 polymerase-based technique to amplify RNA from two tissues of a cichlid fish and compared expression levels of unamplified and amplified RNA on a cDNA microarray. Amplification bias was generally minor and comprised features that were lost (1.3%) or gained (2.5%) through amplification and features that were scored as regulated before but unregulated after amplification (4.2%) or vice versa (19.5%). We examined 10 sequence-specific properties and found that GC content, folding energy, hairpin length and number, and lengths of poly(A) and poly(T) stretches significantly affected RNA amplification. We conclude that, if RNA amplification is used in gene expression studies, preceding experiments controlling for amplification bias should be performed.     

QTL mapping and molecular marker analysis for the resistance of rice to ozone

The resistance of rice to ozone (O3) is a quantitative trait controlled by nuclear genes. The identification of quantitative trait loci (QTL) and analysis of molecular markers of O3 resistance is important for increasing the resistance of rice to O3 stress. QTL associated with the O3 resistance of rice were mapped on chromosomes 1, 7 and 11 using 164 recombinant inbred (RI) lines from a cross between 'Milyang 23' and 'Gihobyeo'. The quantitative trait loci were tightly linked to the markers RG109, C507 and RG1094 and were detected in each of three replications. The association between these markers and O3 resistance in 26 rice cultivars and doubled haploid (DH) populations was analysed. The markers permit the screening of rice germplasm for O3 resistance and the introduction of resistance into elite lines in breeding programs.     

Fine mapping of the qCTS12 locus, a major QTL for seedling cold tolerance in rice

The temperate japonica rice cultivar M202 is the predominant variety grown in California due to its tolerance to low temperature stress, good grain quality and high yield. Earlier analysis of a recombinant inbred line mapping population derived from a cross between M202 and IR50, an indica cultivar that is highly sensitive to cold stress, resulted in the identification of a number of QTL conferring tolerance to cold-induced wilting and necrosis. A major QTL, qCTS12, located on the short arm of chromosome 12, contributes over 40% of the phenotypic variance. To identify the gene(s) underlying qCTS12, we have undertaken the fine mapping of this locus. Saturating the short arm of chromosome 12 with microsatellite markers revealed that qCTS12 is closest to RM7003. Using RM5746 and RM3103, which are immediately outside of RM7003, we screened 1,954 F₅-F₁₀ lines to find recombinants in the qCTS12 region. Additional microsatellite markers were identified from publicly available genomic sequence and used to fine map qCTS12 to a region of approximately 87 kb located on the BAC clone OSJNBb0071I17. This region contains ten open reading frames (ORFs) consisting of five hypothetical and expressed proteins of unknown function, a transposon protein, a putative NBS-LRR disease resistance protein, two zeta class glutathione S-transferases (OsGSTZ1 and OsGSTZ2), and a DAHP synthetase. Further fine mapping with markers developed from the ORFs delimited the QTL to a region of about 55 kb. The most likely candidates for the gene(s) underlying qCTS12 are OsGSTZ1 and OsGSTZ2.The mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture.     

QTL mapping for the paste viscosity characteristics in rice (Oryza sativa L.)

In order to understand the genetic basis of the paste viscosity characteristics (RVA profile, which is tested on the Rapid Visco Analyser) of the rice grain, we mapped QTLs for RVA profile parameters using a DH population derived from a cross between an indica variety, Zai-Ye-Qing 8 (ZYQ8), and a japonica variety, Jing-Xi 17 (JX17). Evidence of genotype-by-environment interaction was found by comparing the mapped QTLs between two locations, Hainan (HN) and Hangzhou (HZ). A total of 20 QTLs for six parameters of the RVA profiles were identified at least one location. Only the waxy locus (wx) located on chromosome 6 was detected significantly at both environments for five traits, i.e. hot paste viscosity (HPV), cool paste viscosity (CPV), breakdown viscosity (BDV), consistency viscosity (CSV) and setback viscosity (SBV). This locus explained 19.5%–63.7% of the total variations at both environments, suggesting that the RVA profiles were mainly controlled by the wx gene. HPV, CPV, BDV, CSV and SBV were also controlled by other QTLs whose effects on the respective parameter were detected only in one environment, while for the peak viscosity (PKV), only 2 QTLs, 1 at HN,the other at HZ, were identified. These results indicate that RVA profiles are obviously affected by environment. Received: 18 July 1999 / Accepted: 27 August 1999     

Fine mapping SPP1, a QTL controlling the number of spikelets per panicle, to a BAC clone in rice (Oryza sativa)

Near isogenic lines (NILs) can be used to efficiently handle a target quantitative trait locus (QTL) by blocking genetic background noise. One QTL, SPP1, which controls the number of spikelets per panicle (SPP), was located on chromosome 1, near Gn1a, a cloned gene for rice production in a recombinant inbred line population. NILs of the SPP1 regions were quickly obtained by self-crossing recombinant inbred line 30 which is heterozygous around SPP1. Using a random NIL-F₂ population of 210 individuals, we mapped SPP1 to a 2.2-cM interval between RM1195 and RM490, which explained 51.1% of SPP variation. The difference in SPP between the two homozygotes was 44. F₂-1456, one NIL-F₂ plant, was heterozygous in the SPP1 region but was fixed in the region of Gn1a gene. This plant F₃ family showed a very wide variation in SPP, which suggested that it was SPP1 but Gn1a affected the variation of SPP in this population. In a word, SPP1 is a novel gene distinct from Gn1a. Four newly developed InDel markers were used for high-resolution mapping of SPP1 with a large NIL-F₂ population. Finally, it was narrowed down to a bacterial artificial chromosome clone spanning 107 kb; 17 open reading frames have been identified in the region. Of them, LOC_Os01g12160, which encodes an IAA synthetase, is the most interesting candidate gene.     

Validation and fine mapping of a QTL for ovulation rate on swine chromosome 3

Ovulation rate (OR) is an important component of litter size, but mutation(s) in gene(s) underlying OR QTL have yet to be identified in pigs. Markers within an OR QTL on SSC3 were genotyped in three white composite lines selected for ten generations for increased OR or uterine capacity (UC), with one line being an unselected control. Numbers of corpora lutea (CL) and UC (number of fully formed fetuses) were collected at approximately 105 days of gestation, as well as ovary weight (OW), uterine length (UL) and uterine weight (UW) measurements at 160 d of age in generation 12 and 13 females from all three lines. Six microsatellites and ten single nucleotide polymorphisms (SNPs; 0–42 cM) were genotyped in pigs from all lines of generations 11 through 13. The allele frequencies of 24269.1, SW2429, 7907.2 and 7637.2 were different (P < 0.01) in the OR line compared to the control line. A significant (P < 0.05) association of CL with 24269.1 (additive effect 0.65 ± 0.32) was detected, and additive genotypic effects approached significance for markers at 28 through 35 cM (16963.2, 27514.1 and SWR1637). Haplotyping of 7637.2 and 16963.2 (31 through 32 cM) identified a significant additive association of haplotype 1 with CL (?0.62 ± 0.30). These markers were also associated with OW (24296.1 and SWR1637), UL (16963.2, 27514.1 and haplotypes of 7637.2/16963.2) and UW (haplotypes of 7637.2/16963.2). This study verifies an OR QTL on SSC3. However, based on the data, it was concluded that there may be two genes, at 13 through 18 cM and 28 through 35 cM, controlling OR on SSC3p.     

QTL mapping of grain weight in rice and the validation of the QTL qTGW3.2

A recombinant inbred line (RIL) population bred from a cross between a javanica type (cv. D50) and an indica type (cv. HB277) rice was used to map seven quantitative trait loci (QTLs) for thousand grain weight (TGW). The loci were distributed on chromosomes 2, 3, 5, 6, 8 and 10. The chromosome 3 QTL qTGW3.2 was stably expressed over two years, and contributed 9–10% of the phenotypic variance. A residual heterozygous line (RHL) was selected from the RIL population and its selfed progeny was used to fine map qTGW3.2. In this “F₂” population, the QTL explained about 23% of the variance, rising to nearly 33% in the subsequent “F_2:3” generation. The physical location of qTGW3.2 was confined to a ~ 556 kb region flanked by the microsatellite loci RM16162 and RM16194. The region also contains other factors influencing certain yield-related traits, although it is also possible that qTGW3.2 affects these in a pleiotropic fashion.     

A mixture model approach to the mapping of QTL controlling endosperm traits with bulked samples

Endosperm traits are of triploid inheritance and have become a focus of breeding effort for their close relations with the grain quality. Current methods for mapping quantitative trait loci (QTL) underlying endosperm traits are restricted to the use of the phenotypes of single grain samples as input data set, which are often not available in practice due to the small size of the cereal seeds. This paper proposed a statistical model for one specially tailored mapping strategy, where the marker genotypes are obtained from the maternal plants in the segregation population and the phenotypic responses are replaced by the trait means of composite endosperm samples pooled from each plant. It should therefore be more practical and have wide applicability in mapping endosperm traits. The method was implemented by fitting the phenotypic means of endosperms into a Gaussian mixture model. Both the exact and approximate Expectation-Maximization algorithms were proposed to estimate the model parameters. The presence of the QTL was determined by likelihood ratio test statistics. Statistical power and other properties of the new method were investigated and compared to the current single-seed method under a variety of scenarios through simulation studies. The simulations suggest a reasonable sample size should be used to ensure reliable results. The proposed method was also applied to a simulated genome data for further evaluation. As an illustration, a real data of maize was analyzed to find the loci responsible for the popping expansion volume.     

A novel brassinolide-enhanced gene identified by cDNA microarray is involved in the growth of rice

Brassinosteroids (BRs) are growth-promoting natural substances required for normal plant growth and development. To understand the molecular mechanism of BR action, a cDNA microarray containing 1265 rice genes was analyzed for expression differences in rice lamina joint treated with brassinolide (BL). A novel BL-enhanced gene, designated OsBLE2, was identified and cloned. The full-length cDNA is 3243 bp long, encoding a predicted polypeptide of 761 amino acid residues and nine possible transmembrane regions. OsBLE2 expression was most responsive to BL in the lamina joint and leaf sheath in rice seedlings. Besides, auxin and gibberellins also increased its expression. OsBLE2 expressed more, as revealed by in situ hybridization, in vascular bundles and root primordia, where the cells are actively undergoing division, elongation, and differentiation. Transgenic rice expressing antisense OsBLE2 exhibits various degrees of repressed growth. BL could not enhance its expression in transgenic rice expressing antisense BRI1, a BR receptor, indicating that BR signaling to the enhanced expression of OsBLE2 is through BRI1. BL effect in the d1 mutant rice was much weaker than that in its wild-type control, indicating that heterotrimeric G protein may be a component of BRs signaling. These results suggest that OsBLE2 is involved in BL-regulated growth and development processes in rice.