Tissue-Specific Transcriptomic Profiling of Sorghum propinquum using a Rice Genome Array

Sorghum (Sorghum bicolor) is one of the world''s most important cereal crops. S. propinquum is a perennial wild relative of S. bicolor with well-developed rhizomes. Functional genomics analysis of S. propinquum, especially with respect to molecular mechanisms related to rhizome growth and development, can contribute to the development of more sustainable grain, forage, and bioenergy cropping systems. In this study, we used a whole rice genome oligonucleotide microarray to obtain tissue-specific gene expression profiles of S. propinquum with special emphasis on rhizome development. A total of 548 tissue-enriched genes were detected, including 31 and 114 unique genes that were expressed predominantly in the rhizome tips (RT) and internodes (RI), respectively. Further GO analysis indicated that the functions of these tissue-enriched genes corresponded to their characteristic biological processes. A few distinct cis-elements, including ABA-responsive RY repeat CATGCA, sugar-repressive TTATCC, and GA-responsive TAACAA, were found to be prevalent in RT-enriched genes, implying an important role in rhizome growth and development. Comprehensive comparative analysis of these rhizome-enriched genes and rhizome-specific genes previously identified in Oryza longistaminata and S. propinquum indicated that phytohormones, including ABA, GA, and SA, are key regulators of gene expression during rhizome development. Co-localization of rhizome-enriched genes with rhizome-related QTLs in rice and sorghum generated functional candidates for future cloning of genes associated with rhizome growth and development.     

Genetic control of rhizomes and genomic localization of a major-effect growth habit QTL in perennial wildrye

Rhizomes are prostrate subterranean stems that provide primitive mechanisms of vegetative dispersal, survival, and regrowth of perennial grasses and other monocots. The extent of rhizome proliferation varies greatly among grasses, being absent in cereals and other annuals, strictly confined in caespitose perennials, or highly invasive in some perennial weeds. However, genetic studies of rhizome proliferation are limited and genes controlling rhizomatous growth habit have not been elucidated. Quantitative trait loci (QTLs) controlling rhizome spreading were compared in reciprocal backcross populations derived from hybrids of rhizomatous creeping wildrye (Leymus triticoides) and caespitose basin wildrye (L. cinereus), which are perennial relatives of wheat. Two recessive QTLs were unique to the creeping wildrye backcross, one dominant QTL was unique to the basin wildrye backcross, and one additive QTL was detectable in reciprocal backcrosses with high log odds (LOD = 31.6) in the basin wildrye background. The dominant QTL located on linkage group (LG)-2a was aligned to a dominant rhizome orthogene (Rhz3) of perennial rice (Oryza longistamina) and perennial sorghum (Sorghum propinquum). Nonparametric 99 % confidence bounds of the 31.6-LOD QTL were localized to a distal 3.8-centiMorgan region of LG-6a, which corresponds to a 0.7-Mb region of Brachypodium Chromosome 3 containing 106 genes. An Aux/IAA auxin signal factor gene was located at the 31.6-LOD peak, which could explain the gravitropic and aphototropic behavior of rhizomes. Findings elucidate genetic mechanisms controlling rhizome development and architectural growth habit differences among plant species. Results have possible applications to improve perennial forage and turf grasses, extend the vegetative life cycle of annual cereals, such as wheat, or control the invasiveness of highly rhizomatous weeds such as quackgrass (Elymus repens).     

Analysis of ESTs from a Normalized cDNA Library of the Rhizome Tip of <Emphasis Type="Italic">Oryza longistaminata</Emphasis>

Oryza longistaminata, a perennial wild rice species with an AA genome, is characterized by the presence of rhizomatous stems. The rhizomatous trait in rice was previously shown to be quantitatively controlled by many genes, but the molecular mechanism related to rhizome development is still unknown. In the present study, expressed sequence tags (ESTs) generated from rhizome tips of O. longistaminata were collected and analyzed. A total of 10,283 complimentary deoxyribunucleic acid clones were randomly sequenced, which generated 10,136 raw sequences, and finally, 4,419 unisequences with diverse functional categories were generated. These unisequences were mapped onto the Oryza sativa genome, which revealed that 4,285 (96.97%) and 4,151 (93.94%) of the unisequences were alignable to the japonica and indica genomic sequences, respectively, with >80% sequence identity. Additionally, 41 unisequences showed four typical types of alternative splicing patterns. More than 600 simple sequence repeats were identified in these unisequences. A subset of unisequences were physically colocalized onto rhizome-related quantitative trait locus intervals in rice and sorghum, and one gene, OLRR1, was further confirmed to be enriched in the rhizome tip and young leaf by real-time polymerase chain reaction and in situ hybridization. Unisequences reported in this study provide valuable data for molecular dissection of the rhizomatous growth habit in O. longistaminata.     

Insight into the molecular evolution of non-specific lipid transfer proteins via comparative analysis between rice and sorghum

Phylogenetic analysis was conducted on 9 kDa non-specific lipid transfer protein (nsLTP) genes from nine plant species. Each of the five classified types in angiosperms exhibited eight conserved cysteine patterns. The most abundant nsLTP genes fell into the type I category, which was particularly enriched in a grass-specific lineage of clade I.1. Six pairs of tandem copies of nsLTP genes on the distal region of rice chromosomes 11 and 12 were well-preserved under concerted evolution, which was not observed in sorghum. The transgenic promoter-reporter assay revealed that both rice and sorghum nsLTP genes of type I displayed a relatively conserved expression feature in the epidermis of growing tissue, supporting its functional roles in cutin synthesis or defence against phytopathogens. For type I, the frequent expression in the stigma and seed are indicative of functional involvement in pistil-pollen interactions and seed development. By way of contrast, several type V genes were observed, mainly in the vascular bundle of the rosette as well as the young shoots, which might be related with vascular tissue differentiation or defence signalling. Compared with sorghum, the highly redundant tissue-specific expression pattern among members of rice nsLTP genes in clade I.1 suggests that concerted evolution via gene conversion favours the preservation of crucial expression motifs via the homogenization of proximal promoter sequences under high selection constraints. However, extensive regulatory subfunctionalization might also have occurred under relative low selection constraints, resulting in functional divergence at the expression level.     

Comparative analysis of QTLs affecting plant height and flowering among closely-related diploid and polyploid genomes.

Quantitative trait loci (QTLs) affecting plant height and flowering were studied in the two Saccharum species from which modern sugarcane cultivars are derived. Two segregating populations derived from interspecific crosses between Saccharum officinarum and Saccharum spontaneum were genotyped with 735 DNA markers. Among the 65 significant associations found between these two traits and DNA markers, 35 of the loci were linked to sugarcane genetic maps and 30 were unlinked DNA markers. Twenty-one of the 35 mapped QTLs were clustered in eight genomic regions of six sugarcane homologous groups. Some of these could be divergent alleles at homologous loci, making the actual number of genes implicated in these traits much less than 35. Four QTL clusters controlling plant height in sugarcane corresponded closely to four of the six plant-height QTLs previously mapped in sorghum. One QTL controlling flowering in sugarcane corresponded to one of three flowering QTLs mapped in sorghum. The correspondence in locations of QTLs affecting plant height and flowering in sugarcane and sorghum reinforce the notion that the simple sorghum genome is a valuable "template" for molecular dissection of the much more complex sugarcane genome.     

Genome-Wide Survey and Comparative Analysis of LTR Retrotransposons and Their Captured Genes in Rice and Sorghum

Long terminal repeat (LTR) retrotransposons are the major class I mobile elements in plants. They play crucial roles in gene expansion, diversification and evolution. However, their captured genes are yet to be genome-widely identified and characterized in most of plants although many genomes have been completely sequenced. In this study, we have identified 7,043 and 23,915 full-length LTR retrotransposons in the rice and sorghum genomes, respectively. High percentages of rice full-length LTR retrotransposons were distributed near centromeric region in each of the chromosomes. In contrast, sorghum full-length LTR retrotransposons were not enriched in centromere regions. This dissimilarity could be due to the discrepant retrotransposition during and after divergence from their common ancestor thus might be contributing to species divergence. A total of 672 and 1,343 genes have been captured by these elements in rice and sorghum, respectively. Gene Ontology (GO) and gene set enrichment analysis (GSEA) showed that no over-represented GO term was identified in LTR captured rice genes. For LTR captured sorghum genes, GO terms with functions in DNA/RNA metabolism and chromatin organization were over-represented. Only 36% of LTR captured rice genes were expressed and expression divergence was estimated as 11.9%. Higher percentage of LTR captured rice genes have evolved into pseudogenes under neutral selection. On the contrary, higher percentage of LTR captured sorghum genes were under purifying selection and 72.4% of them were expressed. Thus, higher percentage of LTR captured sorghum genes was functional. Small RNA analysis suggested that some of LTR captured genes in rice and sorghum might have been involved in negative regulation. On the other hand, positive selection has been observed in both rice and sorghum LTR captured genes and some of them were still expressed and functional. The data suggest that some of these LTR captured genes might have evolved into new gene functions.     

F-box proteins in rice. Genome-wide analysis, classification, temporal and spatial gene expression during panicle and seed development, and regulation by light and abiotic stress

F-box proteins constitute a large family in eukaryotes and are characterized by a conserved F-box motif (approximately 40 amino acids). As components of the Skp1p-cullin-F-box complex, F-box proteins are critical for the controlled degradation of cellular proteins. We have identified 687 potential F-box proteins in rice (Oryza sativa), the model monocotyledonous plant, by a reiterative database search. Computational analysis revealed the presence of several other functional domains, including leucine-rich repeats, kelch repeats, F-box associated domain, domain of unknown function, and tubby domain in F-box proteins. Based upon their domain composition, they have been classified into 10 subfamilies. Several putative novel conserved motifs have been identified in F-box proteins, which do not contain any other known functional domain. An analysis of a complete set of F-box proteins in rice is presented, including classification, chromosomal location, conserved motifs, and phylogenetic relationship. It appears that the expansion of F-box family in rice, in large part, might have occurred due to localized gene duplications. Furthermore, comprehensive digital expression analysis of F-box protein-encoding genes has been complemented with microarray analysis. The results reveal specific and/or overlapping expression of rice F-box protein-encoding genes during floral transition as well as panicle and seed development. At least 43 F-box protein-encoding genes have been found to be differentially expressed in rice seedlings subjected to different abiotic stress conditions. The expression of several F-box protein-encoding genes is also influenced by light. The structure and function of F-box proteins in plants is discussed in light of these results and the published information. These data will be useful for prioritization of F-box proteins for functional validation in rice.     

Efficacy of microarray profiling data combined with QTL mapping for the identification of a QTL gene controlling the initial growth rate in rice

Seedling vigor, which is controlled by many quantitative trait loci (QTLs), is one of several important agronomic traits for direct-seedling rice systems. However, isolating these QTL genes is laborious and expensive. Here, we combined QTL mapping and microarray profiling to identify QTL genes for seedling vigor. By performing QTL mapping using 82 backcross inbred lines (BILs) of the Koshihikari (japonica) and Habataki (indica) cultivars for the rice initial growth, we identified two QTLs, early-stage plant development1/2 (qEPD1 and qEPD2), whose Koshihikari alleles promote plant height and/or leaf sheath length. Phenotypic analysis of the two substituted lines carrying the Habataki qEPD1 or qEPD2 allele revealed that qEPD2 functioned more dominantly for the initial growth of rice. From the microarray experiment, 55 and 45 candidate genes were found in the qEPD1 and qEPD2 genomic regions, which are expressed differentially between each substitution line (SL) and Koshihikari. Gibberellin 20 oxidase-2 (OsGA20ox2), which is identical to Semi Dwarf1 (SD1), was included among the 55 candidate genes of qEPD1, whereas its paralog, OsGA20ox1, was included among the 45 candidate genes of qEPD2. Consistently, introduction of the Koshihikari OsGA20ox1 allele into SL(qEPD2) increaseed its plant height and leaf sheath length significantly relative to the introduction of the Habataki OsGA20ox1 allele. Therefore, microarray profiling could be useful for rapidly screening QTL candidate genes. We concluded that OsGA20ox1 and OsGA20ox2 (SD1) function during the initial growth of rice, but OsGA20ox1 plays a dominant role in increasing plant height and leaf sheath length at the initial growth stage.     

水稻蒸煮品质相关QTL定位及候选基因分析

挖掘与稻米蒸煮品质相关的数量性状基因座(quantitative trait locus, QTL)，分析候选基因，并通过遗传育种手段改良稻米蒸煮品质相关性状，可有效提升稻米的口感。以籼稻华占(Huazhan, HZ)、粳稻热研2号(Nekken2)及由其构建的120个重组自交系(recombinant inbred lines, RILs)群体为实验材料，测定成熟期稻米的糊化温度(gelatinization temperature, GT)、胶稠度(gel consistency, GC)和直链淀粉含量(amylose content, AC)。结合高密度分子遗传图谱进行QTL定位，共检测到26个与稻米蒸煮品质相关的QTLs (糊化温度相关位点1个、胶稠度相关位点13个、直链淀粉含量相关位点12个)，其中最高奇数的可能性(likelihood of odd, LOD)值达30.24。通过实时荧光定量PCR (quantitative real-time polymerase chain reaction, qRT-PCR)分析定位区间内候选基因的表达量，发现6个基因在双亲间的表达量差异显著，推测LOC_Os04g20270和LOC_Os11g40100的高表达可能会极大地提高稻米的胶稠度，而LOC_Os01g04920和LOC_Os02g17500的高表达以及LOC_Os03g02650和LOC_Os05g25840的低表达有助于降低直链淀粉含量。这些结果为培育优质水稻新品种奠定了分子基础，并为揭示稻米蒸煮品质的分子调控机制提供了重要的遗传资源。     

Mapping QTLs for Plant Phenology and Production Traits Using Indica Rice (Oryza sativa L.) Lines Adapted to Rainfed Environment

Drought is a major abiotic stress limiting rice production and yield stability in rainfed ecosystems. Identifying quantitative trait loci (QTL) for rice yield and yield components under water limited environments will help to develop drought resilient cultivars using marker assisted breeding (MAB) strategy. A total of 232 recombinant inbred lines of IR62266/Norungan were used to map QTLs for plant phenology and production traits under rainfed condition in target population of environments. A total of 79 QTLs for plant phenology and production traits with phenotypic variation ranging from 4.4 to 72.8% were detected under non-stress and drought stress conditions across two locations. Consistent QTLs for phenology and production traits were detected across experiments and water regimes. The QTL region, RM204-RM197-RM217 on chromosome 6 was linked to days to 50% flowering and grain yield per plant under both rainfed and irrigated conditions. The same genomic region, RM585-RM204-RM197 was also linked to harvest index under rainfed condition with positive alleles from Norungan, a local landrace. QTLs for plant production and drought resistance traits co-located near RM585-RM204-RM197-RM217 region on chromosome 6 in several rice genotypes. Thus with further fine mapping, this region may be useful as a candidate QTL for MAB, map-based cloning of genes and functional genomics studies for rainfed rice improvement.     

Gene expression polymorphisms and ESTs associated with gravitropic response of subterranean branch meristems and growth habit in Leymus wildryes

Negatively orthogeotropic (NOGT) tiller and diageotropic (DGT) rhizome meristems develop from the same type of lateral axillary meristems and phytomer structure. Although subterranean NOGT and DGT buds appear similar, they display different responses to gravity and perhaps other cues governing branch angle and overall growth habit (GH). Leymus wildryes show remarkable variation in GH and include some of the largest native grasses in western North America. Previous studies detected GH QTLs on homoeologous regions of LG3a and LG3b controlling differences between caespitose Leymus cinereus and rhizomatous Leymus triticoides allotetraploids. Heterologous barley and wheat microarrays in conjunction with bulk segregate analysis were used to find gene expression polymorphisms associated with GH QTLs. Approximately 34% and 25% of the probe sets showed detectable signals on the barley and wheat arrays, respectively. Overall gene expression patterns of NOGT and DGT meristems were remarkably similar, consistent with the assertion that Leymus NOGT and DGT buds develop from homologous meristems. Only 28 and 27 genes on barley and wheat gene chips, respectively, showed more than twofold differential expressions between NOGT and DGT tissues. One expression polymorphism genetically mapped in the Leymus LG3 rhizome QTL region.     

Identification of genes contributing to quantitative disease resistance in rice

Despite the importance of quantitative disease resistance during a plant’s life, little is known about the molecular basis of this type of host-pathogen interaction, because most of the genes underlying resistance quantitative trait loci (QTLs) are unknown. To identify genes contributing to resistance QTLs in rice, we analyzed the colocalization of a set of characterized rice defense-responsive genes and resistance QTLs against different pathogens. We also examined the expression patterns of these genes in response to pathogen infection in the parents of the mapping populations, based on the strategy of validation and functional analysis of the QTLs. The results suggest that defense-responsive genes are important resources of resistance QTLs in rice. OsWRKY45-1 is the gene contributing to a major resistance QTL.NRR,OsGH3-1,and OsGLP members on chromosome 8 contribute alone or collectively to different minor resistance QTLs. These genes function in a basal resistance pathway or in major disease resistance gene-mediated race-specific pathways.     

Genome-wide analysis of WOX gene family in rice, sorghum, maize, Arabidopsis and poplar

WUSCHEL-related homeobox(WOX)genes form a large gene family specifically expressed in plants.They are known to play important roles in regulating the development of plant tissues and organs by determining cell fate.Recent available whole genome sequences allow us to do more comprehensive phylogenetic analysis of the WOX genes in plants.In the present study,we identified 11 and 21 WOXs from sorghum(Sorghum bicolor)and maize(Zea mays),respectively.The 72 WOX genes from rice(Oryza sativa),sorghum,maize,Arabidopsis(Arabidopsis thaliana)and poplar(Populus trichocarpa)were grouped into three well supported clades with nine subgroups according to the amino acid sequences of their homodomains.Their phylogenetic relationship was also supported by the observation of the motifs outside the homodomain.We observed the variation of duplication events among the nine sub-groups between monocots and eudicots,for instance,more gene duplication events of WOXs within subgroup A for monocots,while,less for dicots in this subgroup.Furthermore,we observed the conserved intron/exon structural patterns of WOX genes in rice,sorghum and Arabidopsis.In addition,WUS(Wuschel)-box and EAR(the ERF-associated amphiphilic repression)-like motif were observed to be conserved among several WOX subgroups in these five plants.Comparative analysis of expression patterns of WOX genes in rice and Arabidopsis suggest that the WOX genes play conserved and various roles in plants.This work provides insights into the evolution of the WOX gene family and is useful for future research.     

云南元江普通野生稻株高和抽穗期QTL定位研究

以云南元江普通野生稻为供体亲本,在特青的遗传背景下构建了一套BC3高代回交群体。利用117个SSR标记分析383个BC3F2株系的基因型,采用单标记分析法对控制元江普野株高和抽穗期的QTL进行分析。在北京和合肥两个地点试验结果表明,控制株高的QTL分布在第1染色体上,在RM104附近有一个QTL,与sd-1位置相当,其对表现型变异的贡献率在两个地点分别为27％和28％,其加性效应值分别为26．24cm和26．28cm,来自野生稻的等位基因显著提高回交群体的株高;在第1、3、7、8、11染色体共检测到6个控制抽穗期QTL,其中第8染色体RM25附近控制抽穗期的QTL在两个地点的贡献率分别为13％和15％,加性效应值为4．60d和3．65d,来自野生稻的等位基因使回交群体抽穗期延迟。     

A method for identification of the expression mode and mapping of QTL underlying embryo-specific characters

Embryos of crop seeds are one of the major sources of the plant protein and lipid for human nutrition. The genetic expression for embryo-specific characters in crop seeds can be controlled exclusively by the embryo or the maternal genotypes and sometimes by both simultaneously. However, current methods for mapping quantitative trait loci (QTLs) underlying characters of maternal plants have not been effective in dealing with the QTL analysis of embryo characters. On the basis of the expression feature of embryo, a statistical method was proposed for the identification of expression mode and mapping of QTL controlling embryo traits. The maximum likelihood method implemented via the expectation maximization algorithm was used to estimate parameters of a putative embryo-specific QTL. The QTL expression mode was identified by the likelihood ratio test statistic. Statistical power and other properties of the proposed method were investigated under a variety of scenarios through simulation studies. The results showed that the mapping method neglecting the effects of embryo genotype or maternal effects could neither identify the expression mode of QTL nor estimate its genetic effects accurately, whereas the proposed method could effectively map the embryo-specific QTL of various expression modes.