Yingxin Zhang

Transcriptome analysis of rice root heterosis by RNA-Seq

BackgroundHeterosis is a phenomenon in which hybrids exhibit superior performance relative to parental phenotypes. In addition to the heterosis of above-ground agronomic traits on which most existing studies have focused, root heterosis is also an indispensable component of heterosis in the entire plant and of major importance to plant breeding. Consequently, systematic investigations of root heterosis, particularly in reproductive-stage rice, are needed. The recent advent of RNA sequencing technology (RNA-Seq) provides an opportunity to conduct in-depth transcript profiling for heterosis studies.ResultsUsing the Illumina HiSeq 2000 platform, the root transcriptomes of the super-hybrid rice variety Xieyou 9308 and its parents were analyzed at tillering and heading stages. Approximately 391 million high-quality paired-end reads (100-bp in size) were generated and aligned against the Nipponbare reference genome. We found that 38,872 of 42,081 (92.4%) annotated transcripts were represented by at least one sequence read. A total of 829 and 4186 transcripts that were differentially expressed between the hybrid and its parents (DGHP) were identified at tillering and heading stages, respectively. Out of the DGHP, 66.59% were down-regulated at the tillering stage and 64.41% were up-regulated at the heading stage. At the heading stage, the DGHP were significantly enriched in pathways related to processes such as carbohydrate metabolism and plant hormone signal transduction, with most of the key genes that are involved in the two pathways being up-regulated in the hybrid. Several significant DGHP that could be mapped to quantitative trait loci (QTLs) for yield and root traits are also involved in carbohydrate metabolism and plant hormone signal transduction pathways.ConclusionsAn extensive transcriptome dataset was obtained by RNA-Seq, giving a comprehensive overview of the root transcriptomes at tillering and heading stages in a heterotic rice cross and providing a useful resource for the rice research community. Using comparative transcriptome analysis, we detected DGHP and identified a group of potential candidate transcripts. The changes in the expression of the candidate transcripts may lay a foundation for future studies on molecular mechanisms underlying root heterosis.

Fine mapping of the lesion mimic and early senescence 1 (lmes1) in rice (Oryza sativa).

A novel rice mutant, lesion mimic and early senescence 1 (lmes1), was induced from the rice 93-11 cultivar in a γ-ray field. This mutant exhibited spontaneous disease-like lesions in the absence of pathogen attack at the beginning of the tillering stage. Moreover, at the booting stage, lmes1 mutants exhibited a significantly increased MDA but decreased chlorophyll content, soluble protein content and photosynthetic rate in the leaves, which are indicative of an early senescence phenotype. The lmes1 mutant was significantly more resistant than 93-11 against rice bacterial blight infection, which was consistent with a marked increase in the expression of three resistance-related genes. Here, we employed a map-based cloning approach to finely map the lmes1 gene. In an initial mapping with 94 F2 individuals derived from a cross between the lmes1 mutant and Nipponbare, the lmes1 gene was located in a 10.6-cM region on the telomere of the long arm of chromosome 7 using simple sequence repeat (SSR) markers. To finely map lmes1, we derived two F2 populations with 940 individuals from two crosses between the lmes1 mutant and two japonica rice varieties, Nipponbare and 02428. Finally, the lmes1 gene was mapped to an 88-kb region between two newly developed inDel markers, Zl-3 and Zl-22, which harbored 15 ORFs.

Identification of Transcriptome SNPs for Assessing Allele-Specific Gene Expression in a Super-Hybrid Rice Xieyou9308

Hybridization, a common process in nature, can give rise to a vast reservoir of allelic variants. Combination of these allelic variants may result in novel patterns of gene action and is thought to contribute to heterosis. In this study, we analyzed genome-wide allele-specific gene expression (ASGE) in the super-hybrid rice variety Xieyou9308 using RNA sequencing technology (RNA-Seq). We identified 9325 reliable single nucleotide polymorphisms (SNPs) distributed throughout the genome. Nearly 68% of the identified polymorphisms were CT and GA SNPs between R9308 and Xieqingzao B, suggesting the existence of DNA methylation, a heritable epigenetic mark, in the parents and their F1 hybrid. Of 2793 identified transcripts with consistent allelic biases, only 480 (17%) showed significant allelic biases during tillering and/or heading stages, implying that trans effects may mediate most transcriptional differences in hybrid offspring. Approximately 67% and 62% of the 480 transcripts showed R9308 allelic expression biases at tillering and heading stages, respectively. Transcripts with higher levels of gene expression in R9308 also exhibited R9308 allelic biases in the hybrid. In addition, 125 transcripts were identified with significant allelic expression biases at both stages, of which 74% showed R9308 allelic expression biases. R9308 alleles may tend to preserve their characteristic states of activity in the hybrid and may play important roles in hybrid vigor at both stages. The allelic expression of 355 transcripts was highly stage-specific, with divergent allelic expression patterns observed at different developmental stages. Many transcripts associated with stress resistance were differently regulated in the F1 hybrid. The results of this study may provide valuable insights into molecular mechanisms of heterosis.

OsCUL3a Negatively Regulates Cell Death and Immunity by Degrading OsNPR1 in Rice

OsCUL3a interacts with and degrades OsNPR1 to regulate innate immunity in rice. Cullin3-based RING E3 ubiquitin ligases (CRL3), composed of Cullin3 (CUL3), RBX1, and BTB proteins, are involved in plant immunity, but the function of CUL3 in the process is largely unknown. Here, we show that rice (Oryza sativa) OsCUL3a is important for the regulation of cell death and immunity. The rice lesion mimic mutant oscul3a displays a significant increase in the accumulation of flg22- and chitin-induced reactive oxygen species, and in pathogenesis-related gene expression as well as resistance to Magnaporthe oryzae and Xanthomonas oryzae pv oryzae. We cloned the OsCUL3a gene via a map-based strategy and found that the lesion mimic phenotype of oscul3a is associated with the early termination of OsCUL3a protein. Interaction assays showed that OsCUL3a interacts with both OsRBX1a and OsRBX1b to form a multisubunit CRL in rice. Strikingly, OsCUL3a interacts with and degrades OsNPR1, which acts as a positive regulator of cell death in rice. Accumulation of OsNPR1 protein is greater in the oscul3a mutant than in the wild type. Furthermore, the oscul3a osnpr1 double mutant does not exhibit the lesion mimic phenotype of the oscul3a mutant. Our data demonstrate that OsCUL3a negatively regulates cell death and immunity by degrading OsNPR1 in rice.

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.

Mapping and Comparative Analysis of QTL for Rice Plant Height Based on Different Sample Sizes within a Single Line in RIL Population

To clarify the most appropriate sample size for obtaining phenotypic data for a single line, we investigated the main-effect QTL (M-QTL) of a quantitative trait plant height (ph) in a recombinant inbred line (RIL) population of rice (derived from the cross between Xieqingzao B and Zhonghui 9308) using five individual plants in 2006 and 2009. Twenty-six ph phenotypic datasets from the completely random combinations of 2, 3, 4, and 5 plants in a single line, and five ph phenotypic datasets from five individual plants were used to detect the QTLs. Fifteen M-QTLs were detected by 1 to 31 datasets. Of these, qph7a was detected repeatedly by all the 31 ph datasets in 2006 and explained 11.67% to 23.93% of phenotypic variation; qph3 was detected repeatedly by all the 31 datasets and explained 5.21% to 7.93% and 11.51% to 24.46% of phenotypic variance in 2006 and 2009, respectively. The results indicate that the M-QTL for a quantitative trait could be detected repeatedly by the phenotypic values from 5 individual plants and 26 sets of completely random combinations of phenotypic data within a single line in an RIL population under different environments. The sample size for a single line of the RIL population did not affect the efficiency for identification of stably expressed M-QTLs.

The Rice AAA-ATPase OsFIGNL1 Is Essential for Male Meiosis

Meiosis is crucial in reproduction of plants and ensuring genetic diversity. Although several genes involved in homologous recombination and DNA repair have been reported, their functions in rice (Oryza sativa) male meiosis remain poorly understood. Here, we isolated and characterized the rice OsFIGNL1 (OsFidgetin-like 1) gene, encoding a conserved AAA-ATPase, and explored its function and importance in male meiosis and pollen formation. The rice Osfignl1 mutant exhibited normal vegetative growth, but failed to produce seeds and displayed pollen abortion phenotype. Phenotypic comparisons between the wild-type and Osfignl1 mutant demonstrated that OsFIGNL1 is required for anther development, and that the recessive mutation of this gene causes male sterility in rice. Complementation and CRISPR/Cas9 experiments demonstrated that wild-type OsFIGNL1 is responsible for the male sterility phenotype. Subcellular localization showed that OsFIGNL1-green fluorescent protein was exclusively localized in the nucleus of rice protoplasts. Male meiosis in the Osfignl1 mutant exhibited abnormal chromosome behavior, including chromosome bridges and multivalent chromosomes at diakinesis, lagging chromosomes, and chromosome fragments during meiosis. Yeast two-hybrid assays demonstrated OsFIGNL1 could interact with RAD51A1, RAD51A2, DMC1A, DMC1B, and these physical interactions were further confirmed by BiFC assay. Taken together, our results suggest that OsFIGNL1 plays an important role in regulation of male meiosis and anther development.

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.

Genetic mapping of a QTL controlling source–sink size and heading date in rice

Source size, sink size and heading date (HD) are three important classes of traits that determine the productivity of rice. In this study, a set of recombinant inbred lines (RILs) derived from the cross between an elite indica line Big Grain1 (BG1) and a japonica line Xiaolijing (XLJ) were used to map quantitative trait loci (QTLs) for source-sink size and heading date. Totally, thirty-one QTLs for source size, twenty-two for sink size, four for heading date and seven QTL clusters which included QTLs for multiple traits were identified in three environmental trials. Thirty QTLs could be consistently detected in at least two trials and generally located in the clusters. Using a set of BC4F2 lines, the QTL cluster in C5-1-C5-2 on chromosome 5 was validated to be a major QTL pleiotropically affecting heading date, source size (flag leaf area) and panicle type (neck length of panicle, primary branching number and the ratio of secondary branching number to primary branching number), and was narrowed down to a 309.52Kb region. QTL clusters described above have a large effect on source-sink size and/or heading date, therefore they should be good resources to improve the adaptability and high yield potential of cultivars genetically.

ES7, encoding a ferredoxin-dependent glutamate synthase, functions in nitrogen metabolism and impacts leaf senescence in rice

Glutamate synthase (GOGAT) is a key enzyme for nitrogen metabolism and ammonium assimilation in plants. In this study, an early senescence 7 (es7) mutant was identified and characterized. The leaves of the es7 mutant begin to senesce at the tillering stage about 60day after sowing, and become increasingly senescent as the plants develop at the heading stage. When es7 plants are grown under photorespiration-suppressed conditions (high CO2), the senescence phenotype and chlorophyll content are rescued. qRT-PCR analysis showed that senescence- associated genes were up-regulated significantly in es7. A map-based cloning strategy was used to identify ES7, which encodes a ferredoxin-dependent glutamate synthase (Fd-GOGAT). ES7 was expressed constitutively, and the ES7 protein was localized in chloroplast. qRT-PCR analysis indicated that several genes related to nitrogen metabolism were differentially expressed in es7. Further, we also demonstrated that chlorophyll synthesis-associated genes were significantly down-regulated in es7. In addition, when seedlings are grown under increasing nitrogen concentrations (NH4NO3) for 15days, the contents of chlorophyll a, chlorophyll b and total chlorophyll were significantly lower in es7. Our results demonstrated that ES7 is involved in nitrogen metabolism, effects chlorophyll synthesis, and may also associated with photorespiration, impacting leaf senescence in rice.