Liyong Cao

Dissection of genetic architecture of rice plant height and heading date by multiple-strategy-based association studies.

Xieyou9308 is a certified super hybrid rice cultivar with a high grain yield. To investigate its underlying genetic basis of high yield potential, a recombinant inbred line (RIL) population derived from the cross between the maintainer line XieqingzaoB (XQZB) and the restorer line Zhonghui9308 (ZH9308) was constructed for identification of quantitative trait SNPs (QTSs) associated with two important agronomic traits, plant height (PH) and heading date (HD). By re-sequencing of 138 recombinant inbred lines (RILs), a total of ~0.7 million SNPs were identified for the association studies on the PH and HD. Three association mapping strategies (including hypothesis-free genome-wide association and its two complementary hypothesis-engaged ones, QTL-based association and gene-based association) were adopted for data analysis. Using a saturated mixed linear model including epistasis and environmental interaction, we identified a total of 31 QTSs associated with either the PH or the HD. The total estimated heritability across three analyses ranged from 37.22% to 45.63% and from 37.53% to 55.96% for the PH and HD, respectively. In this study we examined the feasibility of association studies in an experimental population (RIL) and identified several common loci through multiple strategies which could be preferred candidates for further research.

OsCOL16, encoding a CONSTANS-like protein, represses flowering by up-regulating Ghd7 expression in rice

Flowering time is an important agronomic trait that coordinates the plant life cycle with regional adaptability and thereby impacts yield potentials for cereal crops. The CONSTANS (CO)-like gene family plays vital roles in the regulation of flowering time. CO-like proteins are typically divided into four phylogenetic groups in rice. Several genes from groups I, III, and IV have been functionally characterized, though little is known about the genes of group II in rice. We report the functional characterization in rice of a constitutive floral inhibitor, OsCOL16, encoding a group-II CO-like protein that delays flowering time and increases plant height and grain yield. Overexpression of OsCOL16 resulted in late heading under both long-day and short-day conditions. OsCOL16 expression exhibits a diurnal oscillation and serves as a transcription factor with transcriptional activation activity. We determined that OsCOL16 up-regulates the expression of the floral repressor Ghd7, leading to down-regulation of the expression of Ehd1, Hd3a, and RFT1. Moreover, genetic diversity and evolutionary analyses suggest that remarkable differences in flowering times correlate with two major alleles of OsCOL16. Our combined molecular biology and phylogeographic analyses revealed that OsCOL16 plays an important role in regulating rice photoperiodic flowering, allowing for environmental adaptation of rice.

Dissecting Genetic Architecture Underlying Seed Traits in Multiple Environments

The seeds of flowering plants develop from double fertilization and play a vital role in reproduction and supplying human and animal food. The genetic variation of seed traits is influenced by multiple genetic systems, e.g., maternal, embryo, and/or endosperm genomes. Understanding the genetic architecture of seed traits is a major challenge because of this complex mechanism of multiple genetic systems, especially the epistasis within or between different genomes and their interactions with the environment. In this study, a statistical model was proposed for mapping QTL with epistasis and QTL-by-environment (QE) interactions underlying endosperm and embryo traits. Our model integrates the maternal and the offspring genomes into one mapping framework and can accurately analyze maternal additive and dominant effects, endosperm/embryo additive and dominant effects, and epistatic effects of two loci in the same or two different genomes, as well as interaction effects of each genetic component of QTL with environment. Intensive simulations under different sampling strategies, heritabilities, and model parameters were performed to investigate the statistical properties of the model. A set of real cottonseed data was analyzed to demonstrate our methods. A software package, QTLNetwork-Seed-1.0.exe, was developed for QTL analysis of seed traits.

Genetic dissection of QTL against phosphate deficiency in the hybrid rice ‘Xieyou9308’

Phosphate (Pi) plays important roles in plant development and architecture. With the goal of identifying genomic regions that influence tolerance to Pi deficiency (TPDE) in hybrid rice (Oryza sativa L.), quantitative trait loci (QTL) were mapped using recombinant inbred lines (RILs) that were derived from a cross between tolerant ‘XieqingzaoB’ (XB) and susceptible ‘Zhonghui9308’. Six TPDE-related traits, including the root length, root dry weight, tillers number, shoot dry weight, total plant dry weight and root-to-shoot ratio, were evaluated for QTL analysis during both the tillering and heading stages. A correlation analysis showed that most of the traits were correlated with each other. Twenty-one additive QTL were detected and jointly explained between 10–49% of the trait variance, tending to cluster on chromosomes 4, 6, 10 and 11. Three QTL, qTPDE4XB, qTPDE10XB and qTPDE11.3XB, were validated by the phenotypic evaluation using near isogenic lines (NILs, BC4F3) during the seedling stage. qTPDE4XB showed the most stable tolerance against Pi deficiency. These QTL will enrich the genetic resources and accelerate hybrid rice breeding against Pi deficiency.

The rice CONSTANS-like protein OsCOL15 suppresses flowering by promoting Ghd7 and repressing RID1

The photoperiodic flowering pathway is one of the most important regulatory networks controlling flowering time in rice (Oryza sativa L.). Rice is a facultative short-day (SD) plant; flowering is promoted under inductive SD conditions and delayed under non-inductive long-day (LD) conditions. In rice, flowering inhibitor genes play an important role in maintaining the trade-off between reproduction and yield. In this study, we identified a novel floral inhibitor, OsCOL15, which encodes a CONSTANS-like transcription factor. Consistent with a function in transcriptional regulation, OsCOL15 localized to the nucleus. Moreover, OsCOL15 had transcriptional activation activity, and the central region of the protein between the B-box and CCT domains was required for this activity. We determined that OsCOL15 is most highly expressed in young organs and exhibits a diurnal expression pattern typical of other floral regulators. Overexpression of OsCOL15 resulted in a delayed flowering phenotype under both SD and LD conditions. Real-time quantitative RT-PCR analysis of flowering regulator gene expression suggested that OsCOL15 suppresses flowering by up-regulating the flowering repressor Grain number, plant height and heading date 7 (Ghd7) and down-regulating the flowering activator Rice Indeterminate 1 (RID1), thus leading to the down-regulation of the flowering activators Early heading date 1, Heading date 3a, and RICE FLOWERING LOCUS T1. These results demonstrate that OsCOL15 is an important floral regulator acting upstream of Ghd7 and RID1 in the rice photoperiodic flowering-time regulatory network.

Detection of Epistatic and Gene-Environment Interactions Underlying Three Quality Traits in Rice Using High-Throughput Genome-Wide Data

With development of sequencing technology, dense single nucleotide polymorphisms (SNPs) have been available, enabling uncovering genetic architecture of complex traits by genome-wide association study (GWAS). However, the current GWAS strategy usually ignores epistatic and gene-environment interactions due to absence of appropriate methodology and heavy computational burden. This study proposed a new GWAS strategy by combining the graphics processing unit- (GPU-) based generalized multifactor dimensionality reduction (GMDR) algorithm with mixed linear model approach. The reliability and efficiency of the analytical methods were verified through Monte Carlo simulations, suggesting that a population size of nearly 150 recombinant inbred lines (RILs) had a reasonable resolution for the scenarios considered. Further, a GWAS was conducted with the above two-step strategy to investigate the additive, epistatic, and gene-environment associations between 701,867 SNPs and three important quality traits, gelatinization temperature, amylose content, and gel consistency, in a RIL population with 138 individuals derived from super-hybrid rice Xieyou9308 in two environments. Four significant SNPs were identified with additive, epistatic, and gene-environment interaction effects. Our study showed that the mixed linear model approach combining with the GPU-based GMDR algorithm is a feasible strategy for implementing GWAS to uncover genetic architecture of crop complex traits.

qSE7 is a major quantitative trait locus (QTL) influencing stigma exsertion rate in rice ( Oryza sativa L.)

Stigma exsertion is a key determinant to increase the efficiency of commercial hybrid rice seed production. The major quantitative trait locus (QTL) qSE7 for stigma exsertion rate was previously detected on the chromosome 7 using 75 Chromosome Segment Substitution Lines (CSSLs) derived from a cross between the high stigma exsertion indica maintainer XieqingzaoB (XQZB) and low stigma exsertion indica restorer Zhonghui9308 (ZH9308). The C51 line, a CSSL population with an introgression from XQZB, was backcrossed with ZH9308 to produce the secondary F2 (BC5F2) and F2:3 (BC5F2:3) populations. As a result, the Near Isogenic Line (NIL qSE7XB) was developed. Analysis indicated qSE7 acted as a single Mendelian factor and decreased the stigma exsertion. We hypothesized qSE7 regulates single, dual, and total stigma exsertion rate, provided experimental support. qSE7 was mapped and localized between RM5436 and RM5499 markers, within a physical distance of 1000-kb. With use of new insertion-deletion (InDel) markers and analysis of the heterozygous and phenotypic data, it was ultimately dissected to a 322.9-kb region between InDel SER4-1 and RM5436. The results are useful for additional identification and isolation of this candidate gene controlling stigma exsertion rate, and provide a basis for further fine mapping, gene cloning, and Marker Assisted Selection (MAS) breeding later.

Genomic Regions Analysis of Seedling Root Traits and Their Regulation in Responses to Phosphorus Deficiency Tolerance in CSSL Population of Elite Super Hybrid Rice

Phosphorus (P) is the essential macro-element supporting rice productivity. Quantitative trait loci (QTL) underlying related traits at the seedling stage under two different phosphorus levels was investigated in rice using a population of 76 Chromosomal Sequence Substitution Lines (CSSLs) derived from a cross between the maintainer variety XieqingzaoB (P stress tolerant) and the restorer variety Zhonghui9308 (P stress sensitive); the parents of super hybrid rice Xieyou9308. A genetic linkage map with 120 DNA marker loci was constructed. At logarithmic odd (LOD) value of 2.0, a total of seven QTLs were detected for studied traits under two P levels and their relative ratio. The LOD values ranged from 2.00 to 3.32 and explaining 10.82% to 18.46% of phenotypic variation. Three QTLs were detected under low phosphorus (P−), one under normal (P+) and three under their relative ratio (P−/P+) on the rice chromosomes 3, 5, 6, 8 and 10. No significant QTLs were found for shoot dry weight (SDW) and total dry weight (TDW). The pleiotropic QTLs influencing root number (qRN5) and root dry weight (qRDW5) as novel QTLs under P− level were detected near marker RM3638 on chromosome 5, which considered to directly contributing to phosphorus deficiency tolerance in rice. These QTLs need further analysis, including the fine mapping and cloning, which may use in molecular marker assisted breeding.

WB1, a Regulator of Endosperm Development in Rice, Is Identified by a Modified MutMap Method

Abnormally developed endosperm strongly affects rice (Oryza sativa) appearance quality and grain weight. Endosperm formation is a complex process, and although many enzymes and related regulators have been identified, many other related factors remain largely unknown. Here, we report the isolation and characterization of a recessive mutation of White Belly 1 (WB1), which regulates rice endosperm development, using a modified MutMap method in the rice mutant wb1. The wb1 mutant develops a white-belly endosperm and abnormal starch granules in the inner portion of white grains. Representative of the white-belly phenotype, grains of wb1 showed a higher grain chalkiness rate and degree and a lower 1000-grain weight (decreased by ~34%), in comparison with that of Wild Type (WT). The contents of amylose and amylopectin in wb1 significantly decreased, and its physical properties were also altered. We adopted the modified MutMap method to identify 2.52 Mb candidate regions with a high specificity, where we detected 275 SNPs in chromosome 4. Finally, we identified 19 SNPs at 12 candidate genes. Transcript levels analysis of all candidate genes showed that WB1 (Os04t0413500), encoding a cell-wall invertase, was the most probable cause of white-belly endosperm phenotype. Switching off WB1 with the CRISPR/cas9 system in Japonica cv. Nipponbare demonstrates that WB1 regulates endosperm development and that different mutations of WB1 disrupt its biological function. All of these results taken together suggest that the wb1 mutant is controlled by the mutation of WB1, and that the modified MutMap method is feasible to identify mutant genes, and could promote genetic improvement in rice.

Identification of cyp703a3-3 and analysis of regulatory role of CYP703A3 in rice anther cuticle and pollen exine development

Anther cuticle and pollen exine are two elaborated lipid-soluble barriers protecting pollen grains from environmental and biological stresses. However, less is known about the mechanisms underlying the synthesis of these lipidic polymers. Here, we identified a no-pollen male-sterility mutant cyp703a3-3 from the indica restorer line Zhonghui 8015 (Zh8015) mutant library treated with 60Coγ-ray radiation. Histological analysis indicated that cyp703a3-3 underwent abnormal tapetal cells development, produced few orbicules and secreted less sporopollenin precursors to anther locule, as well as cutin monomers on anther. Genetic analysis revealed that cyp703a3-3 was controlled by a single recessive gene. Map-based cloning was performed to narrow down the mutant gene to a 47.78-kb interval on the chromosome 8 between two markers S15-29 and S15-30. Sequence analysis detected three bases (GAA) deletion in the first exon of LOC_Os08g03682, annotated as CYP703A3 with homologous sequences related to male sterility in Arabidopsis, causing the Asparagine deletion in the mutant site. Moreover, we transformed genomic fragment of CYP703A3 into cyp703a3-3, which male-sterility phenotype was recovered. Both the wild-type and cyp703a3-3 mutant 3D structure of CYP703A3 protein were modeled. Results of qPCR suggested CYP703A3 mainly expressed in anthers with greatest abundance at microspore stage, and genes involved in sporopollenin precursors formation and transportation, such as GAMYB, TDR, CYP704B2, DPW2, OsABCG26 and OsABCG15, were significantly reduced in cyp703a3-3. Collectively, our results further elaborated CYP703A3 plays vital role in anther cuticle and pollen exine development in rice (Oryza sativa L.).