Chengjun Zhang

Genetic composition of yield heterosis in an elite rice hybrid

Heterosis refers to the superior performance of hybrids relative to the parents. Utilization of heterosis has contributed tremendously to the increased productivity in many crops for decades. Although there have been a range of studies on various aspects of heterosis, the key to understanding the biological mechanisms of heterotic performance in crop hybrids is the genetic basis, much of which is still uncharacterized. In this study, we dissected the genetic composition of yield and yield component traits using data of replicated field trials of an “immortalized F2” population derived from an elite rice hybrid. On the basis of an ultrahigh-density SNP bin map constructed with population sequencing, we calculated single-locus and epistatic genetic effects in the whole genome and identified components pertaining to heterosis of the hybrid. The results showed that the relative contributions of the genetic components varied with traits. Overdominance/pseudo-overdominance is the most important contributor to heterosis of yield, number of grains per panicle, and grain weight. Dominance × dominance interaction is important for heterosis of tillers per plant and grain weight and has roles in yield and grain number. Single-locus dominance has relatively small contributions in all of the traits. The results suggest that cumulative effects of these components may adequately explain the genetic basis of heterosis in the hybrid.

Natural variation in Ghd7.1 plays an important role in grain yield and adaptation in rice

Natural variation in Ghd7.1 plays an important role in grain yield and adaptation in rice

Sequence and expression analysis of the thioredoxin protein gene family in rice.

Thioredoxin (Trx) proteins play important biological functions in cells by changing redox via thioldisulfied exchange. This system is especially widespread in plants. Through database search, we identified 30 potential Trx protein-encoding genes (OsTrx) in rice (Oryza sativa L.). An analysis of the complete set of OsTrx proteins is presented here, including chromosomal location, conserved motifs, domain duplication, and phylogenetic relationships. Our findings suggest that the expansion of the Trx gene family in rice, in large part, occurred due to gene duplication. A comprehensive expression profile of Trx genes family was investigated by analyzing the signal data of this family extracted from the whole genome microarray analysis of Minghui 63 and Zhenshan 97, two indica parents, and their hybrid Shanyou 63, using 27 different tissues representing the entire life cycle of rice. Results revealed specific expression of some members at germination transition as well as the 3-leaf stage during the vegetative growth phase of rice. OsTrx genes were also found to be differentially up- or down-regulated in rice seedlings subjected to treatments of phytohormones and light/dark conditions. The expression levels of the OsTrx genes in the different tissues and under different treatments were also checked by RT-PCR analysis. The identification of OsTrx genes showing differential expression in specific tissues among different genotypes or in response to different environmental cues could provide a new avenue for functional analyses in rice.

CHD3 protein recognizes and regulates methylated histone H3 lysines 4 and 27 over a subset of targets in the rice genome

Histone lysine methylation is an important component of the epigenetic system demarcating transcriptionally active and inactive chromatin domains. It is of primary importance in understanding how different histone lysine methylation marks and a specific combination of them are read and interpreted by chromatin proteins to regulate gene expression. In this paper, we report that the rice CHD3 protein CHR729 that was required for many aspects of plant development can interact with dimethylated histone H3 lysine 4 (H3K4me2, a mark associated with moderately expressed or repressed genes) and with trimethylated histone H3 lysine 27 (H3K27me3, a mark associated with repressed genes), respectively, through the chromodomains and the plant homeodomain (PHD) finger of the protein. A mutation or down-regulation of the gene provoked a decrease of H3K27me3 and H3K4me3 (a mark associated with active genes). Genome-wide analysis revealed that H3K27me3 and H3K4me3, respectively, were lost from about 56 and 23% of marked loci, which correspond mostly to under-expressed or repressed genes. In the mutant, a higher-than-expected proportion of down-regulated genes lost H3K4me3, among which many encode DNA-binding transcription factors. These results suggest that the rice CHD3 protein is a bifunctional chromatin regulator able to recognize and modulate H3K4 and H3K27 methylation over repressed or tissue-specific genes, which may be associated with regulation of a gene transcription program of plant development.

gKaKs: the pipeline for genome-level Ka/Ks calculation

SUMMARY gKaKs is a codon-based genome-level Ka/Ks computation pipeline developed and based on programs from four widely used packages: BLAT, BLASTALL (including bl2seq, formatdb and fastacmd), PAML (including codeml and yn00) and KaKs_Calculator (including 10 substitution rate estimation methods). gKaKs can automatically detect and eliminate frameshift mutations and premature stop codons to compute the substitution rates (Ka, Ks and Ka/Ks) between a well-annotated genome and a non-annotated genome or even a poorly assembled scaffold dataset. It is especially useful for newly sequenced genomes that have not been well annotated. We applied gKaKs to estimate the genome-wide substitution rates in five pairs of closely related species. The average Ka and Ks computed by gKaKs were consistent with previous studies. We also compared the Ka, Ks and Ka/Ks of mouse and rat orthologous protein-coding genes estimated by gKaKs and based on the alignments generated by PAL2NAL. Results from two methods are compatible. AVAILABILITY AND IMPLEMENTATION gKaKs is implemented in Perl and is freely available on http://longlab.uchicago.edu/?q=gKaKs. The detailed user manual is available on the website.

Comprehensive sequence and expression profile analysis of PEX11 gene family in rice.

PEX11 gene family has been shown to be involved in peroxisome biogenesis but very little is known about this gene family in rice. Here we show that five putative PEX11 genes (OsPEX11-1-5) present in rice genome and each contain three conserved motifs. The PEX11 sequences from rice and other species can be classified into three major groups. Among the five rice PEX11 genes, OsPEX11-2 and -3 are most likely duplicated. Expression profile and RT-PCR analysis suggested that the members of PEX11 family in rice had differential expression patterns: OsPEX11-1 and OsPEX11-4 had higher expression levels in leaf tissues than in the other tissues, OsPEX11-2 was detected only in germinated seeds, OsPEX11-3 was expressed predominantly in endosperm and germinated seeds, and OsPEX11-5 was expressed in all the tissues investigated. We also observed that the rice PEX11 genes had differential expression patterns under different abiotic stresses. OsPEX11-1 and OsPEX11-4 were induced by abscisic acid (ABA), hydrogen peroxide (H2O2), salt and low nitrogen stress conditions. OsPEX11-3 was responsive to ABA and H2O2 treatments, and OsPEX11-5 was responsive to ABA, H2O2, and salt treatments. However, OsPEX11-2 had no response to any of the stresses. Our results suggest that the rice PEX11 genes have diversification not only in sequences but also in expression patterns under normal and various stress conditions.

Genomes of 13 domesticated and wild rice relatives highlight genetic conservation, turnover and innovation across the genus Oryza

The genus Oryza is a model system for the study of molecular evolution over time scales ranging from a few thousand to 15 million years. Using 13 reference genomes spanning the Oryza species tree, we show that despite few large-scale chromosomal rearrangements rapid species diversification is mirrored by lineage-specific emergence and turnover of many novel elements, including transposons, and potential new coding and noncoding genes. Our study resolves controversial areas of the Oryza phylogeny, showing a complex history of introgression among different chromosomes in the young ‘AA’ subclade containing the two domesticated species. This study highlights the prevalence of functionally coupled disease resistance genes and identifies many new haplotypes of potential use for future crop protection. Finally, this study marks a milestone in modern rice research with the release of a complete long-read assembly of IR 8 ‘Miracle Rice’, which relieved famine and drove the Green Revolution in Asia 50 years ago.Genome assemblies of 13 domesticated and wild rice relatives reveal salient features of genome evolution across the genus Oryza, especially rapid species diversification and turnover of transposons. This study also releases a complete long-read assembly of IR 8 ‘Miracle Rice’.

Single feature polymorphisms between two rice cultivars detected using a median polish method

Expression levels measured in microarrays of oligonucleotide probes have now been adapted as a high throughput approach for identifying DNA sequence variation between genotypes, referred to as single feature polymorphisms (SFPs). Although there have been increasing interests in this approach, there is still need for improving the algorithm in order to achieve high sensitivity and specificity especially with complex genome and large datasets, while maintaining optimal computational performance. We obtained microarray datasets for expression profiles of two rice cultivars and adapted a median polish method to detect SFPs. The analysis identified 6,655 SFPs between two the rice varieties representing 3,131 rice unique genes. We showed that the median polish method has the advantage of avoiding fitting complex linear models thus can be used to analyze complex transcriptome datasets like the ones in this study. The method is also superior in sensitivity, accuracy and computing time requirement compared with two previously used methods. A comparison with data from a resequencing project indicated that 75.6% of the SFPs had SNP supports in the probe regions. Further comparison revealed that SNPs in sequences immediately flanking the probes also had contributions to the detection of SFPs in cases where the probes and the targets had perfectly matched sequences. It was shown that differences in minimum free energies caused by flanking SNPs, which may change the stability of RNA secondary structure, may partly explain the SFPs as detected. These SFPs may facilitate gene discovery in future studies.

Evolutionary Character of Alternative Splicing in Plants

Alternative splicing (AS) is one of the most important ways to enhance the functional diversity of genes. Huge amounts of data have been produced by microarray, expressed sequence tag, and RNA-seq, and plenty of methods have been developed specifically for this task. The most frequently asked questions in previous research were as follows. What is the content rate of AS genes among the whole gene set? How many AS types are presented in the genome, and which type is dominant? How about the conservation ability of AS among different species? Which kinds of isoforms from some genes have the environmental response to help individual adaptation? Based on this background, we collected analysis results from 17 species to try to map out the landscape of AS studies in plants. We have noted the shortages of previous results, and we appeal to all scientists working in the AS field to make a standard protocol so that analyses between different projects are comparable.

Evolution of gene structural complexity: An alternative-splicing based model accounts for intron-containing retrogenes

An alternative-splicing-based model explains the formation of retrogenes that retained the parental intron structure and indicates that plants have a much higher percentage of this kind of complex retrogene than animals. The structure of eukaryotic genes evolves extensively by intron loss or gain. Previous studies have revealed two models for gene structure evolution through the loss of introns: RNA-based gene conversion, dubbed the Fink model and retroposition model. However, retrogenes that experienced both intron loss and intron-retaining events have been ignored; evolutionary processes responsible for the variation in complex exon-intron structure were unknown. We detected hundreds of retroduplication-derived genes in human (Homo sapiens), fly (Drosophila melanogaster), rice (Oryza sativa), and Arabidopsis (Arabidopsis thaliana) and categorized them either as duplicated genes that have all introns lost or as duplicated genes that have at least lost one and retained one intron compared with the parental copy (intron-retaining [IR] type). Our new model attributes intron retention alternative splicing to the generation of these IR-type gene pairs. We presented 25 parental genes that have an intron retention isoform and have retained introns in the same locations in the IR-type duplicate genes, which directly support our hypothesis. Our alternative-splicing-based model in conjunction with the retroposition and Fink models can explain the IR-type gene observed. We discovered a greater percentage of IR-type genes in plants than in animals, which may be due to the abundance of intron retention cases in plants. Given the prevalence of intron retention in plants, this new model gives a support that plant genomes have very complex gene structures.