Go Suzuki

Abstract Weedy rice is a major problem in paddy fields around the world. It is well known that weedy rice appears to grow faster and mature earlier than cultivated rice. It is possible that differences in the root microbial genetics are correlated with this characteristic. This study incorporated 16S rRNA amplicon sequencing to study the microbial composition in the rhizosphere and endosphere of rice root. No significant difference was found between the microbiota associated with weedy and cultivated rice lines grown in the same field. It was found that the endosphere had less microbial diversity compared to the rhizosphere. The major groups of bacteria found in the endosphere are from the phylum Proteobacteria, Myxococcota, Chloroflexota, and Actinobacteria. In addition, by analyzing the microbiome of japonica rice grown in the field in a temperate climate, we found that despite differences in genotype and location, some bacterial taxa were found to be common and these members of the putative rice core microbiome can also be detected by in situ hybridization. The delineation of a core microbiome in the endosphere of rice suggests that these bacterial taxa might be important in the life cycle of a wide range of rice types.

Abstract Key message S²⁹ haplotype does not require the MLPK function for self-incompatibility in Brassica rapa. Abstract Self-incompatibility (SI) in Brassicaceae is regulated by the self-recognition mechanism, which is based on the S-haplotype-specific direct interaction of the pollen-derived ligand, SP11/SCR, and the stigma-side receptor, SRK. M locus protein kinase (MLPK) is known to be one of the positive effectors of the SI response. MLPK directly interacts with SRK, and is phosphorylated by SRK in Brassica rapa. In Brassicaceae, MLPK was demonstrated to be essential for SI in B. rapa and Brassica napus, whereas it is not essential for SI in Arabidopsis thaliana (with introduced SRK and SP11/SCR from related SI species). Little is known about what determines the need for MLPK in SI of Brassicaceae. In this study, we investigated the relationship between S-haplotype diversity and MLPK function by analyzing the SI phenotypes of different S haplotypes in a mlpk/mlpk mutant background. The results have clarified that in B. rapa, all the S haplotypes except the S²⁹ we tested need the MLPK function, but the S²⁹ haplotype does not require MLPK for the SI. Comparative analysis of MLPK-dependent and MLPK-independent S haplotype might provide new insight into the evolution of S-haplotype diversity and the molecular mechanism of SI in Brassicaceae.

<title>Abstract</title> Self-incompatibility (SI) is a breeding system that promotes cross-fertilization. In <italic>Brassica</italic>, pollen rejection is induced by a haplotype-specific interaction between pistil determinant SRK (<italic>S</italic> receptor kinase) and pollen determinant SP11 (<italic>S</italic>-locus Protein 11, also named SCR) from the <italic>S</italic>-locus. Although the structure of the <italic>B. rapa S</italic>₉-SRK ectodomain (eSRK) and <italic>S</italic>₉-SP11 complex has been determined, it remains unclear how SRK discriminates self- and nonself-SP11. Here, we uncover the detailed mechanism of self/nonself-discrimination in <italic>Brassica</italic> SI by determining the <italic>S</italic>₈-eSRK–<italic>S</italic>₈-SP11 crystal structure and performing molecular dynamics (MD) simulations. Comprehensive binding analysis of eSRK and SP11 structures reveals that the binding free energies are most stable for cognate eSRK–SP11 combinations. Residue-based contribution analysis suggests that the modes of eSRK–SP11 interactions differ between intra- and inter-subgroup (a group of phylogenetically neighboring haplotypes) combinations. Our data establish a model of self/nonself-discrimination in <italic>Brassica</italic> SI.

Self-compatibility in Arabidopsis thaliana represents the relatively recent disruption of ancestral obligate cross pollination, recognized as one of the prevalent evolutionary pathways in flowering plants, as noted by Darwin. Our previous study found that inversion of the male specificity gene (SP11/SCR) disrupted self-incompatibility, which was restored by overexpressing the SCR with the reversed inversion. However, SCR in A. thaliana has other mutations aside from the pivotal inversion, in both promoter and coding regions, with probable effects on transcriptional regulation. To examine the functional consequences of these mutations, we conducted reciprocal introductions of native promoters and downstream sequences from orthologous loci of self-compatible A. thaliana and self-incompatible A. halleri. Use of this inter-species pair enabled us to expand the scope of the analysis to transcriptional regulation and deletion in the intron, in addition to inversion in the native genomic background. Initial analysis revealed that A. thaliana has a significantly lower basal expression level of SCR transcripts in the critical reproductive stage compared to that of A. halleri, suggesting that the promoter was attenuated in inducing transcription in A. thaliana. However, in reciprocal transgenic experiments, this A. thaliana promoter was able to restore partial function if coupled with the functional A. halleri coding sequence, despite extensive alterations due to the self-compatible mode of reproduction in A. thaliana. This represents a synergistic effect of the promoter and the inversion resulting in fixation of self-compatibility, primarily enforced by disruption of SCR. Our findings elucidate the functional and evolutionary context of the historical transition in A. thaliana thus contributing to the understanding of the molecular events leading to development of self-compatibility.

© 2019, Genetics Society of Japan. All rights reserved. Self-incompatibility (SI) is a sophisticated system for pollen selectivity to prevent pollination by genetically identical pollen. In Brassica, it is genetically controlled by a single, highly polymorphic S-locus, and the male and female S-determinant factors have been identified as S-locus protein 11 (SP11)/S-locus cysteine-rich protein (SCR) and S-locus receptor kinase (SRK), respectively. However, the overall molecular system and identity of factors in the downstream cascade of the SI reaction remain unclear. Previously, we identified a self-compatible B. rapa mutant line, TSC28, which has a disruption in an unidentified novel factor of the SI signaling cascade. Here, in a genetic analysis of TSC28, using an F2 population from a cross with the reference B. rapa SI line Chiifu-401, the causal gene was mapped to a genetic region of DNA containing markers BrSA64 and ACMP297 in B. rapa chromosome A1. By fine mapping using an F2 population of 1,034 plants, it was narrowed down to a genetic region between DNA markers ACMP297 and BrgMS4028, with physical length approximately 1.01 Mbp. In this genomic region, 113 genes are known to be located and, among these, we identified 55 genes that were expressed in the papilla cells. These are candidates for the gene responsible for the disruption of SI in TSC28. This list of candidate genes will contribute to the discovery of a novel downstream factor in the SP11–SRK signaling cascade in the Brassica SI system.

Stigmatic papillae develop at the apex of the gynoecium and play an important role as a site of pollination. The papillae in Brassicaceae are of the dry and unicellular type, and more than 15,000 genes are expressed in the papillae; however, the molecular and physiological mechanisms of their development remain unknown. We found that the papillae in Arabidopsis thaliana change their length in response to altered ambient humidity: papillae of flowers incubated under high humidity elongated more than those under normal humidity conditions. Genetic analysis and transcriptome data suggest that an abscisic acid-mediated abiotic stress response mechanism regulates papilla length. Our data suggest a flexible regulation of papilla elongation at the post-anthesis stage, in response to abiotic stress, as an adaptation to environmental conditions.

In plants, cell-cell recognition is a crucial step in the selection of optimal pairs of gametes to achieve successful propagation of progeny. Flowering plants have evolved various genetic mechanisms, mediated by cell-cell recognition, to enable their pistils to reject self-pollen, thus preventing inbreeding and the consequent reduced fitness of progeny (self-incompatibility, SI), and to reject foreign pollen from other species, thus maintaining species identity (interspecific incompatibility) 1. In the genus Brassica, the SI system is regulated by an S-haplotype-specific interaction between a stigma-expressed female receptor (S receptor kinase, SRK) and a tapetum cell-expressed male ligand (S locus protein 11, SP11), encoded by their respective polymorphic genes at the S locus2-6. However, the molecular mechanism for recognition of foreign pollen, leading to reproductive incompatibility, has not yet been identified. Here, we show that recognition between a novel pair of proteins, a pistil receptor SUI1 (STIGMATIC UNILATERAL INCOMPATIBILITY 1) and a pollen ligand PUI1 (POLLEN UNILATERAL INCOMPATIBILITY 1), triggers unilateral reproductive incompatibility between plants of two geographically distant self-incompatible Brassica rapa lines, even though crosses would be predicted to be compatible based on the S haplotypes of pollen and stigma. Interestingly, SUI1 and PUI1 are similar to the SI genes, SRK and SP11, respectively, and are maintained as cryptic incompatibility genes in these two populations. The duplication of the SRK and SP11 followed by reciprocal loss in different populations would provide a molecular mechanism of the emergence of a reproductive barrier in allopatry.

In diploid organisms, phenotypic traits are often biased by effects known as Mendelian dominant-recessive interactions between inherited alleles. Phenotypic expression of SP11 alleles, which encodes the male determinants of self-incompatibility in Brassica rapa, is governed by a complex dominance hierarchy1-3. Here, we show that a single polymorphic 24 nucleotide small RNA, named SP11 methylation inducer 2 (Smi2), controls the linear dominance hierarchy of the four SP11 alleles (S44 &gt S60 &gt S40 &gt S29). In all dominant-recessive interactions, small RNA variants derived from the linked region of dominant SP11 alleles exhibited high sequence similarity to the promoter regions of recessive SP11 alleles and acted in trans to epigenetically silence their expression. Together with our previous study4, we propose a new model: sequence similarity between polymorphic small RNAs and their target regulates mono-allelic gene expression, which explains the entire five-phased linear dominance hierarchy of the SP11 phenotypic expression in Brassica.

Background: Application of apomixis, or asexual seed formation, in crop breeding would allow rapid fixation of complex traits, economizing improved crop delivery. Identification of apomixis genes is confounded by the polyploid nature, high genome complexity and lack of genomic sequence integration with reproductive tissue transcriptomes in most apomicts. Results: A genomic and transcriptomic resource was developed for Hieracium subgenus Pilosella (Asteraceae) which incorporates characterized sexual, apomictic and mutant apomict plants exhibiting reversion to sexual reproduction. Apomicts develop additional female gametogenic cells that suppress the sexual pathway in ovules. Disrupting small RNA pathways in sexual Arabidopsis also induces extra female gametogenic cells; therefore, the resource was used to examine if changes in small RNA pathways correlate with apomixis initiation. An initial characterization of small RNA pathway genes within Hieracium was undertaken, and ovary-expressed ARGONAUTE genes were identified and cloned. Comparisons of whole ovary transcriptomes from mutant apomicts, relative to the parental apomict, revealed that differentially expressed genes were enriched for processes involved in small RNA biogenesis and chromatin silencing. Small RNA profiles within mutant ovaries did not reveal large-scale alterations in composition or length distributions; however, a small number of differentially expressed, putative small RNA targets were identified. Conclusions: The established Hieracium resource represents a substantial contribution towards the investigation of early sexual and apomictic female gamete development, and the generation of new candidate genes and markers. Observed changes in small RNA targets and biogenesis pathways within sexual and apomictic ovaries will underlie future functional research into apomixis initiation in Hieracium.

In quantitative gene expression analysis, normalization using a reference gene as an internal control is frequently performed for appropriate interpretation of the results. Efforts have been devoted to exploring superior novel reference genes using microarray transcriptomic data and to evaluating commonly used reference genes by targeting analysis. However, because the number of specifically detectable genes is totally dependent on probe design in the microarray analysis, exploration using microarray data may miss some of the best choices for the reference genes. Recently emerging RNA sequencing (RNA-seq) provides an ideal resource for comprehensive exploration of reference genes since this method is capable of detecting all expressed genes, in principle including even unknown genes. We report the results of a comprehensive exploration of reference genes using public RNA-seq data from plants such as Arabidopsis thaliana (Arabidopsis), Glycine max (soybean), Solanum lycopersicum (tomato) and Oryza sativa (rice). To select reference genes suitable for the broadest experimental conditions possible, candidates were surveyed by the following four steps: (1) evaluation of the basal expression level of each gene in each experiment; (2) evaluation of the expression stability of each gene in each experiment; (3) evaluation of the expression stability of each gene across the experiments; and (4) selection of top-ranked genes, after ranking according to the number of experiments in which the gene was expressed stably. Employing this procedure, 13, 10, 12 and 21 top candidates for reference genes were proposed in Arabidopsis, soybean, tomato and rice, respectively. Microarray expression data confirmed that the expression of the proposed reference genes under broad experimental conditions was more stable than that of commonly used reference genes. These novel reference genes will be useful for analyzing gene expression profiles across experiments carried out under various experimental conditions.

Plants subjected to abiotic stress can regulate gene expression post-transcriptionally by means of small RNAs such as microRNAs. Cool-temperature stress causes abnormal tapetum hypertrophy in rice anthers, leading to pollen sterility. As a first step toward understanding the molecular mechanisms of cool tolerance in developing anthers of rice, we report here a comprehensive comparative analysis of microRNAs between cool-sensitive Sasanishiki and cool-tolerant Hitomebore cultivars. High-throughput Illumina sequencing revealed 241 known and 46 novel microRNAs. Interestingly, 15 of these microRNAs accumulated differentially in the two cultivars at the uninucleate microspore stage under cool conditions. Inverse correlations between expression patterns of microRNAs and their target genes were confirmed by quantitative RT-PCR analysis, and cleavage sites of some of the target genes were determined by 5' RNA ligase-mediated RACE experiments. Thus, our data are useful resources to elucidate microRNA-mediated mechanism(s) of cool tolerance in rice anthers at the booting stage.

Pollination is an important early step in sexual plant reproduction. In Arabidopsis thaliana, sequential pollination events, from pollen adhesion onto the stigma surface to pollen tube germination and elongation, occur on the stigmatic papilla cells. Following successful completion of these events, the pollen tube penetrates the stigma and finally fertilizes a female gametophyte. The pollination events are thought to be initiated and regulated by interactions between papilla cells and pollen. Here, we report the characterization of gene expression profiles of unpollinated (UP), compatible pollinated (CP) and incompatible pollinated (IP) papilla cells in A. thaliana. Based on cell type-specific transcriptome analysis from a combination of laser microdissection and RNA sequencing, 15,475, 17,360 and 16,918 genes were identified as expressed in UP, CP and IP papilla cells, respectively, and, of these, 14,392 genes were present in all three data sets. Differentially expressed gene (DEG) analyses identified 147 and 71 genes up-regulated in CP and IP papilla cells, respectively, and 115 and 46 genes down-regulated. Gene Ontology and metabolic pathway analyses revealed that papilla cells play an active role as the female reproductive component in pollination, particularly in information exchange, signal transduction, internal physiological changes and external morphological modification. This study provides fundamental information on the molecular mechanisms involved in pollination in papilla cells, furthering our understanding of the reproductive role of papilla cells.

Comprehensive integration of large-scale omics resources such as genomes, transcriptomes and metabolomes will provide deeper insights into broader aspects of molecular biology. For better understanding of plant biology, we aim to construct a next-generation sequencing (NGS)-derived gene expression network (GEN) repository for a broad range of plant species. So far we have incorporated information about 745 high-quality mRNA sequencing (mRNA-Seq) samples from eight plant species (Arabidopsis thaliana, Oryza sativa, Solanum lycopersicum, Sorghum bicolor, Vitis vinifera, Solanum tuberosum, Medicago truncatula and Glycine max) from the public short read archive, digitally profiled the entire set of gene expression profiles, and drawn GENs by using correspondence analysis (CA) to take advantage of gene expression similarities. In order to understand the evolutionary significance of the GENs from multiple species, they were linked according to the orthology of each node (gene) among species. In addition to other gene expression information, functional annotation of the genes will facilitate biological comprehension. Currently we are improving the given gene annotations with natural language processing (NLP) techniques and manual curation. Here we introduce the current status of our analyses and the web database, PODC (Plant Omics Data Center; http://bioinf.mind.meiji.ac.jp/podc/), now open to the public, providing GENs, functional annotations and additional comprehensive omics resources.

Apomixis or asexual seed formation in Hieracium praealtum (Asteraceae) is controlled by two independent dominant loci. One of these, the LOSS OF APOMEIOSIS (LOA) locus, controls apomixis initiation, mitotic embryo sac formation (apospory) and suppression of the sexual pathway. The LOA locus is found near the end of a hemizygous chromosome surrounded by extensive repeats extending along the chromosome arm. Similar apomixis-carrying chromosome structures have been found in some apomictic grasses, suggesting that the extensive repetitive sequences may be functionally relevant to apomixis. Fluorescence in situ hybridization (FISH) was used to examine chromosomes of apomeiosis deletion mutants and rare recombinants in the critical LOA region arising from a cross between sexual Hieracium pilosella and apomictic H. praealtum. The combined analyses of aposporous and nonaposporous recombinant progeny and chromosomal karyotypes were used to determine that the functional LOA locus can be genetically separated from the very extensive repeat regions found on the LOA-carrying chromosome. The large-scale repetitive sequences associated with the LOA locus in H. praealtum are not essential for apospory or suppression of sexual megasporogenesis (female meiosis).

Puroindoline-a (Pina), puroindoline-b (Pinb) genes in the wheat hardness-locus region encode 15-kDa friabilin proteins, whose accumulation in the endosperm leads to grain softness texture. In wheat, the PINA and PINB friabilins are associated with starch granules in the endosperm cell, while there is no friabilin in rice. The rice endosperm structure consisting of compound starch granules is fundamentally different from that in wheat. We previously produced two different lines of transgenic rice plants with the large genomic fragment including Pina or Pinb of Aegilops tauschii. However, localization of exogenous friabilins in the rice endosperm cell still remains to be determined. In the present study, we stacked the two different transgenic rice lines. The F-4 seeds of the stacked line, in which the homozygosity of the Pina and Pinb transgenes was checked by FISH analysis, were used for histochemical analysis of the endosperm cell. Immunodetection of PINA and PINB proteins using the Durotest antibody showed that wheat-derived friabilins were localized between compound starch granules as well as between starch granules in the rice endosperm cell. This suggests that such localization of the friabilins might prevent tight interaction between the compound starch granules and between the starch granules in the rice endosperm, leading to its soft texture.

Pollination is an early and critical step in plant reproduction, leading to successful fertilization. It consists of many sequential processes, including adhesion of pollen grains onto the surface of stigmatic papilla cells, foot formation to strengthen pollen-stigma interaction, pollen hydration and germination, and pollen tube elongation and penetration. We have focused on an examination of the expressed genes in papilla cells, to increase understanding of the molecular systems of pollination. From three representative species of Brassicaceae (Arabidopsis thaliana, A. halleri and Brassica rapa), stigmatic papilla cells were isolated precisely by laser microdissection, and cell type-specific gene expression in papilla cells was determined by RNA sequencing. As a result, 17,240, 19,260 and 21,026 unigenes were defined in papilla cells of A. thaliana, A. halleri and B. rapa, respectively, and, among these, 12,311 genes were common to all three species. Among the17,240 genes predicted in A. thaliana, one-third were papilla specific while approximately half of the genes were detected in all tissues examined. Bioinformatics analysis revealed that genes related to a wide range of reproduction and development functions are expressed in papilla cells, particularly metabolism, transcription and membrane-mediated information exchange. These results reflect the conserved features of general cellular function and also the specific reproductive role of papilla cells, highlighting a complex cellular system regulated by a diverse range of molecules in these cells. This study provides fundamental biological knowledge to dissect the molecular mechanisms of pollination in papilla cells and will shed light on our understanding of plant reproduction mechanisms.

In plant reproduction, pollination is the initial key process in bringing together the male and female gametophytes. When a pollen grain lands on the surface of the stigma, information is exchanged between the pollen and stigmatic cell to determine whether the pollen grain will be accepted or rejected. If it is accepted, the stigmatic papilla cell supplies water and other resources to the pollen for germination and pollen tube elongation. Cellular processes involving actin are essential for pollen germination and tube growth, and actin-binding proteins regulate these processes by interacting with actin filaments to assemble cytoskeletal structures and actin networks. LIM proteins, which belong to a subfamily of cysteine-rich proteins, are a family of actin-binding proteins in plants, and are considered to be important for formation of the actin cytoskeleton and maintenance of its dynamics. Although the physiological and biochemical characteristics of LIMs have been elucidated in vitro in a variety of cell types, their exact role in pollen germination and pollen tube growth during pollination remained unclear. In this manuscript, we focus on the pollen-specific LIM proteins, AtPLIM2a and AtPLIM2c, and define their biological function during pollination in Arabidopsis thaliana. The atplim2a/atplim2c double knockdown RNAi plants showed a reduced pollen germination, approximately one-fifth of wild type, and slower pollen tube growth in the pistil, that is 80.4 mu m/hr compared to 140.8 mu m/hr in wild type. These defects led to an occasional unfertilized ovule at the bottom of the silique in RNAi plants. Our data provide direct evidence of the biological function of LIM proteins during pollination as actin-binding proteins, modulating cytoskeletal structures and actin networks, and their consequent importance in seed production.

Background and AimsPollination is an important process in the life cycle of plants and is the first step in bringing together the male and female gametophytes for plant reproduction. While pollination has been studied for many years, accurate knowledge of the morphological aspects of this process is still far from complete. This study therefore focuses on a morphological characterization of pollination, using time-series image analysis of self-and cross-pollinations in Brassica rapa.MethodsTime-lapse imaging of pollen behaviour during self-and cross-pollinations was recorded for 90 min, at 1 min intervals, using a stereoscopic microscope. Using time-series digital images of pollination, characteristic features of pollen behaviours during self-and cross-pollinations were studied.Key ResultsPollen exhibited various behaviours in both self-and cross-pollinations, and these were classified into six representative patterns: germination, expansion, contraction, sudden contraction, pulsation and no change. It is noteworthy that in 'contraction' pollen grains shrunk within a short period of 30-50 min, and in 'pulsation' repeated expansion and contraction occurred with an interval of 10 min, suggesting that a dehydration system is operating in pollination. All of the six patterns were observed on an individual stigma with both self-and cross-pollinations, and the difference between self-and cross-pollinations was in the ratios of the different behaviours. With regard to water transport to and from pollen grains, this occurred in multiple steps, before, during and after hydration. Thus, pollination is regulated by a combination of multiple components of hydration, rehydration and dehydration systems. ConclusionsRegulated hydration of pollen is a key process for both pollination and self-incompatibility, and this is achieved by a balanced complex of hydration, dehydration and nutrient supply to pollen grains from stigmatic papilla cells. © 2013 The Author 2013.

Onion can be used in experimental observation of mitotic cell division in plant science because its chromosome is large and easy to observe. However, molecular genetic studies are difficult in onion because of its large genome size, and only limited information of onion genes has been available to date. Here we cloned and characterized an onion homologue of mitotic RAD21 gene, AcRAD21-1, to develop a molecular marker of mitosis. The N-terminal, middle, and C-terminal regions of deduced AcRAD21-1 protein sequence were conserved with Arabidopsis SYN4/AtRAD213 and rice OsRAD21-1, whereas three characteristic types of repetitive motifs (Repeat-1, Repeat-2/2', and Repeat-3) were observed between the conserved regions. Such inserted repetitive amino acid sequences enlarge the AcRAD21-1 protein into almost 200 kDa, which belongs to the largest class of plant proteins. Genomic organization of the AcRAD21-1 locus was also determined, and the possibility of tandem exon duplication in Repeat-2 was revealed. Subsequently, the polyclonal antiserum was raised against the N-terminal region of AcRAD21-1, and purified by affinity chromatography. Immunohistochemical analysis with the purified antibody successfully showed localization of AcRAD21-1 in onion mitosis, suggesting that it can be used as a molecular marker visualizing dynamic movement of cohesin. (C) 2012 Elsevier B.V. All rights reserved.

Plants have evolved many systems to prevent undesirable fertilization. Among these, incompatibility is a well-organized system in which pollen germination or pollen tube growth is inhibited in pistils. We previously found that a novel one-way pollen-stigma incompatibility response [unilateral incompatibility (UI)] occurred between two self-incompatible Brassica rapa plants, a Turkish line, and a Japanese cultivated hybrid variety, "Osome." Pollen from the Turkish line is rejected on the stigma of the Osome line, but the reverse cross is compatible such a UI phenotype closely resembles self-incompatibility (SI). The pollen factor of this UI has been genetically explained by a single locus which is different from the S-locus. In this study, we performed further genetic analyses of this intraspecies UI and showed that the stigma factor was also controlled by a single locus, and we named the loci corresponding to the stigma and pollen factors of the intraspecies UI, stigmatic unilateral incompatibility (SUI), and pollen unilateral incompatibility (PUI) loci, respectively. Interestingly, segregation analyses of SUI and PUI indicated that they are closely linked to each other and behave as a single unit. To investigate the effect of an SI-related gene, MLPK in this UI, we produced segregation lines for SUI and mlpk. A distorted segregation ratio of SUI phenotype in an mlpk background indicated involvement of MLPK in SUI, suggesting the existence of an MLPK-dependent novel pollen-stigma recognition mechanism. © 2013 Takada et al.

Lachrymatory factor synthase (LFS) catalyzes the formation of lachrymatory factor, one of the most distinctive traits of bulb onion (Allium cepa L.). Therefore, we used LFS as a model for a functional gene in a huge genome, and we examined the chromosomal organization of LFS in A. cepa by multiple approaches. The first-level analysis completed the chromosomal assignment of LFS gene to chromosome 5 of A. cepa via the use of a complete set of A. fistulosum-shallot (A. cepa L. Aggregatum group) monosomic addition lines. Subsequent use of an F-2 mapping population from the interspecific cross A. cepa x A. roylei confirmed the assignment of an LFS locus to this chromosome. Sequence comparison of two BAC clones bearing LFS genes, LFS amplicons from diverse germplasm, and expressed sequences from a doubled haploid line revealed variation consistent with duplicated LFS genes. Furthermore, the BAC-FISH study using the two BAC clones as a probe showed that LFS genes are localized in the proximal region of the long arm of the chromosome. These results suggested that LFS in A. cepa is transcribed from at least two loci and that they are localized on chromosome 5.

Genome evolution is a continuous process and genomic rearrangement occurs both within and between species. With the sequencing of the Arabidopsis thaliana genome, comparative genetics and genomics offer new insights into plant biology. The genus Brassica offers excellent opportunities with which to compare genomic synteny so as to reveal genome evolution. During a previous genetic analysis of clubroot resistance in Brassica rapa, we identified a genetic region that is highly collinear with Arabidopsis chromosome 4. This region corresponds to a disease resistance gene cluster in the A. thaliana genome. Relying on synteny with Arabidopsis, we fine-mapped the region and found that the location and order of the markers showed good correspondence with those in Arabidopsis. Microsynteny on a physical map indicated an almost parallel correspondence, with a few rearrangements such as inversions and insertions. The results show that this genomic region of Brassica is conserved extensively with that of Arabidopsis and has potential as a disease resistance gene cluster, although the genera diverged 20 million years ago.

BAC FISH (fluorescence in situ hybridization using bacterial artificial chromosome probes) is a useful cytogenetic technique for physical mapping, chromosome marker screening, and comparative genomics. As a large genomic fragment with repetitive sequences is inserted in each BAC clone, random BAC FISH without adding competitive DNA can unveil complex chromosome organization of the repetitive elements in plants. Here we performed the comparative analysis of the random BAC FISH in monocot plants including species having small chromosomes (rice and asparagus) and those having large chromosomes (hexaploid wheat, onion, and spider lily) in order to understand a whole view of the repetitive element organization in Poales and Asparagales monocots. More unique and less dense dispersed signals of BAC FISH were observed in species with smaller chromosomes in both the Poales and Asparagales species. In the case of large-chromosome species, 75-85% of the BAC clones were detected as dispersed repetitive FISH signals along entire chromosomes. The BAC FISH of Lycoris did not even show localized repetitive patterns (e.g., centromeric localization) of signals.

Self-incompatibility (SI) is defined as the inability to produce zygotes after self-pollination in a fertile hermaphrodite plant, which has stamens and pistils in the same flower. This structural organization of the hermaphrodite flower increases the risk of self-pollination, leading to low genetic diversity. To avoid this problem plants have established several pollination systems, among which the most elegant system is surely SI. The SI trait can be observed in Brassica crops, including cabbage, broccoli, turnip and radish. To produce hybrid seed of these crops efficiently, the SI trait has been employed in an agricultural context. From another point of view, the recognition reaction of SI during pollen-stigma interaction is an excellent model system for cell-cell communication and signal transduction in higher plants. In this review, we describe the molecular mechanisms of SI in Brassicaceae, which have been dissected by genetic, physiological, and biological approaches, and we discuss the future prospects in relation to associated scientific fields and new technologies.<BR><BR>(Communicated by Tsuneyoshi KUROIWA, M.J.A.)

Grain hardness of wheat is determined by the hardness (Ha)-locus region, which contains three friabilin-related genes: puroindoline-a (Pina), puroindoline-b (Pinb) and GSP-1. In our previous study, we produced the transgenic rice plants harboring the large genomic fragment of the Ha-locus region of Aegilops tauschii containing Pina and GSP-1 genes by Agrobacterium-mediated transformation. To examine the effects of the transgenes in the rice endosperms, we firstly confirmed the homozygosity of the T-DNAs in four independent T(2) lines by using fluorescence in situ hybridization (FISH) and DNA gel blot analyses. The transgenes, Pina and GSP-1, were stably expressed in endosperms of the T(3) and T(4) seeds at RNA and protein levels, indicating that the promoters and other regulatory elements on the wheat Ha-locus region function in rice, and that multigene transformation using a large genomic fragment is a useful strategy. The functional contribution of the transgene-derived friabilins to the rice endosperm structure was considered as an increase of spaces between compound starch granules, resulting in a high proportion of white turbidity seeds.

The LOSS OF APOMEIOSIS (LOA) locus is one of two dominant loci known to control apomixis in the eudicot Hieracium praealtum. LOA stimulates the differentiation of somatic aposporous initial cells after the initiation of meiosis in ovules. Aposporous initial cells undergo nuclear proliferation close to sexual megaspores, forming unreduced aposporous embryo sacs, and the sexual program ceases. LOA-linked genetic markers were used to isolate 1.2 Mb of LOA-associated DNAs from H. praealtum. Physical mapping defined the genomic region essential for LOA function between two markers, flanking 400 kb of identified sequence and central unknown sequences. Cytogenetic and sequence analyses revealed that the LOA locus is located on a single chromosome near the tip of the long arm and surrounded by extensive, abundant complex repeat and transposon sequences. Chromosomal features and LOA-linked markers are conserved in aposporous Hieracium caespitosum and Hieracium piloselloides but absent in sexual Hieracium pilosella. Their absence in apomictic Hieracium aurantiacum suggests that meiotic avoidance may have evolved independently in aposporous subgenus Pilosella species. The structure of the hemizygous chromosomal region containing the LOA locus in the three Hieracium subgenus Pilosella species resembles that of the hemizygous apospory-specific genomic regions in monocot Pennisetum squamulatum and Cenchrus ciliaris. Analyses of partial DNA sequences at these loci show no obvious conservation, indicating that they are unlikely to share a common ancestral origin. This suggests convergent evolution of repeat-rich hemizygous chromosomal regions containing apospory loci in these monocot and eudicot species, which may be required for the function and maintenance of the trait.

Although a centromeric DNA fragment of tobacco (Nicotiana tabacum), Nt2-7, has been reported, the overall structure of the centromeres remains unknown. To characterize the centromeric DNA sequences, we conducted a chromatin immunoprecipitation assay using anti-NtCENH3 antibody and chromatins isolated from two ancestral diploid species (Nicotiana sylvestris and Nicotiana tomentosiformis) of N. tabacum and isolated a 178-pb fragment, Nto1 from N. tomentosiformis, as a novel centromeric DNA. Fluorescence in situ hybridization (FISH) showed that Nto1 localizes on 24 out of 48 chromosomes in some cells of a BY-2 cell line. To identify the origins of the Nt2-7 and Nto1, a tobacco bacterial artificial chromosome (BAC) library was constructed from N. tabacum, and then screened by polymerase chain reaction (PCR) with primer sets designed from the Nt2-7 and Not1 DNA sequences. Twelve BAC clones were found to localize on the centromeric regions by FISH. We selected three BAC clones for sequencing and identified two centromeric retrotransposons, NtCR and NtoCR, the DNA sequences of which are similar to that of Nt2-7 and Nto1, respectively. Quantitative PCR analysis using coprecipitated DNA with anti-NtCENH3 clearly showed coexistence of NtCENH3 with both retrotransposons. These results indicate the possibility that these two retrotransposons act as centromeric DNA sequences in tobacco. NtoCR was found to be specific to N. tomentosiformis and T genome of N. tabacum, and a NtCR-like centromeric retrotransposon (TGRIV) exists in tomato. This specificity suggests that the times of amplification of these centromeric retrotransposons were different.

P&gt;Asexual seed formation, or apomixis, in the Hieracium subgenus Pilosella is controlled by two dominant independent genetic loci, LOSS OF APOMEIOSIS (LOA) and LOSS OF PARTHENOGENESIS (LOP). We examined apomixis mutants that had lost function in one or both loci to establish their developmental roles during seed formation. In apomicts, sexual reproduction is initiated first. Somatic aposporous initial (AI) cells differentiate near meiotic cells, and the sexual pathway is terminated as AI cells undergo mitotic embryo sac formation. Seed initiation is fertilization-independent. Using a partially penetrant cytotoxic reporter to inhibit meioisis, we showed that developmental events leading to the completion of meiotic tetrad formation are required for AI cell formation. Sexual initiation may therefore stimulate activity of the LOA locus, which was found to be required for AI cell formation and subsequent suppression of the sexual pathway. AI cells undergo nuclear division to form embryo sacs, in which LOP functions gametophytically to stimulate fertilization-independent embryo and endosperm formation. Loss of function in either locus results in partial reversion to sexual reproduction, and loss of function in both loci results in total reversion to sexual reproduction. Therefore, in these apomicts, sexual reproduction is the default reproductive mode upon which apomixis is superimposed. These loci are unlikely to encode genes essential for sexual reproduction, but may function to recruit the sexual machinery at specific time points to enable apomixis.

Co-expression networks systematically constructed from large-scale transcriptome data reflect the interactions and functions of genes with similar expression patterns and are a powerful tool for the comprehensive understanding of biological events and mining of novel genes. In Arabidopsis (a model dicot plant), high-resolution co-expression networks have been constructed from very large microarray datasets and these are publicly available as online information resources. However, the available transcriptome data of rice (a model monocot plant) have been limited so far, making it difficult for rice researchers to achieve reliable co-expression analysis. In this study, we performed co-expression network analysis by using combined 44 K agilent microarray datasets of rice, which consisted of 33 laser microdissection (LM)-microarray datasets of anthers, and 143 spatiotemporal transcriptome datasets deposited in RicexPro. The entire data of the rice co-expression network, which was generated from the 176 microarray datasets by the Pearson correlation coefficient (PCC) method with the mutual rank (MR)-based cut-off, contained 24,258 genes and 60,441 genes pairs. Using these datasets, we constructed high-resolution co-expression subnetworks of two specific biological events in the anther, "meiosis" and "pollen wall synthesis". The meiosis network contained many known or putative meiotic genes, including genes related to meiosis initiation and recombination. In the pollen wall synthesis network, several candidate genes involved in the sporopollenin biosynthesis pathway were efficiently identified. Hence, these two subnetworks are important demonstrations of the efficiency of co-expression network analysis in rice. Our co-expression analysis included the separated transcriptomes of pollen and tapetum cells in the anther, which are able to provide precise information on transcriptional regulation during male gametophyte development in rice. The co-expression network data presented here is a useful resource for rice researchers to elucidate important and complex biological events.

Onion, Allium cepa, is a model plant for experimental observation of somatic cell division, whose mitotic chromosome is extremely large, and contains the characteristic terminal heterochromatin. Epigenetic status of the onion chromosome is a matter of deep interest from a molecular cytogenetic point of view, because epigenetic marks regulate chromatin structure and gene expression. Here we examined chromosomal distribution of DNA methylation and histone modification in A. cepa in order to reveal the chromatin structure in detail. Immunodetection of 5-methylcytosine (5mC) and in situ nick-translation analysis showed that onion genomic DNA was highly methylated, and the methylated CG dinucleotides were distributed in entire chromosomes. In addition, distributions of histone methylation codes, which occur in close association with DNA methylation, were similar to those of other large genome species. From these results, a highly heterochromatic and less euchromatic state of large onion chromosomes were demonstrated at an epigenetic level.

A diploid organism has two copies of each gene, one inherited from each parent. The expression of two inherited alleles is sometimes biased by the effects known as dominant/recessive relationships, which determine the final phenotype of the organism. To explore the mechanisms underlying these relationships, we have examined the monoallelic expression of S-locus protein 11 genes (SP11), which encode the male determinants of self-incompatibility in Brassica. We previously reported that SP11 expression was monoallelic in some S heterozygotes, and that the promoter regions of recessive SP11 alleles were specifically methylated in the anther tapetum. Here we show that this methylation is controlled by trans-acting small non-coding RNA (sRNA). We identified inverted genomic sequences that were similar to the recessive SP11 promoters in the flanking regions of dominant SP11 alleles. These sequences were specifically expressed in the anther tapetum and processed into 24-nucleotide sRNA, named SP11 methylation inducer (Smi). Introduction of the Smi genomic region into the recessive S homozygotes triggered the methylation of the promoter of recessive SP11 alleles and repressed their transcription. This is an example showing sRNA encoded in the flanking region of a dominant allele acts in trans to induce transcriptional silencing of the recessive allele. Our finding may provide new insights into the widespread monoallelic gene expression systems.

UDP-glucose pyrophosphorylase (UGPase) is an important enzyme in the metabolism of UDP-glucose, a precursor for the synthesis of carbohydrate cell wall components, such as cellulose and callose. The Arabidopsis thaliana genome contains two putative genes encoding UGPase, AtUGP1 and AtUGP2. These genes are expressed in all organs. In order to determine the role of UGPase in vegetative and reproductive organs, we employed a reverse genetic approach using the T-DNA insertion mutants, atugp1 and atugp2. Despite a significant decrease in UGPase activity in both the atugp7 and atugp2 single mutants, no decrease in normal growth and reproduction was observed. In contrast, the atugp1/atugp2 double mutant displayed drastic growth defects and male sterility. At the reproductive phase, in the anthers of atugp1/atugp2, pollen mother cells developed normally, but callose deposition around microspores was absent. Genes coding for enzymes at the subsequent steps in the cellulose and callose synthesis pathway were also down-regulated in the double mutant. Taken together, these results demonstrate that the AtUGP1 and AtUGP2 genes are functionally redundant and UGPase activity is essential for both vegetative and reproductive phases in Arabidopsis. Importantly, male fertility was not restored in the double knockout mutant by an application of external sucrose, whereas vegetative growth was comparable in size with that of the wild type. In contrast, an application of external UDP-glucose recovered male fertility in the double mutant, suggesting that control of UGPase in carbohydrate metabolism is different in the vegetative phase as compared with the reproductive phase in A. thaliana.

Ever since Darwin&apos;s pioneering research, the evolution of self-fertilisation (selfing) has been regarded as one of the most prevalent evolutionary transitions in flowering plants(1,2). A major mechanism to prevent selfing is the self-incompatibility (SI) recognition system, which consists of male and female specificity genes at the S-locus and SI modifier genes(2-4). Under conditions that favour selfing, mutations disabling the male recognition component are predicted to enjoy a relative advantage over those disabling the female component, because male mutations would increase through both pollen and seeds whereas female mutations would increase only through seeds(5,6). Despite many studies on the genetic basis of loss of SI in the predominantly selfing plant Arabidopsis thaliana(7-15), it remains unknown whether selfing arose through mutations in the female specificity gene (S-receptor kinase, SRK), male specificity gene (S-locus cysteine-rich protein, SCR; also known as S-locus protein 11, SP11) or modifier genes, and whether any of them rose to high frequency across large geographic regions. Here we report that a disruptive 213-base-pair (bp) inversion in the SCR gene (or its derivative haplotypes with deletions encompassing the entire SCR-A and a large portion of SRK-A) is found in 95% of European accessions, which contrasts with the genome-wide pattern of polymorphism in European A. thaliana(16,17). Importantly, interspecific crossings using Arabidopsis halleri as a pollen donor reveal that some A. thaliana accessions, including Wei-1, retain the female SI reaction, suggesting that all female components including SRK are still functional. Moreover, when the 213-bp inversion in SCR was inverted and expressed in transgenic Wei-1 plants, the functional SCR restored the SI reaction. The inversion within SCR is the first mutation disrupting SI shown to be nearly fixed in geographically wide samples, and its prevalence is consistent with theoretical predictions regarding the evolutionary advantage of mutations in male components.

Self-incompatibility (SI) in Brassicaceae is sporophytically controlled by a single S-locus with multi allelic variety. The male S determinant, SP11/SCR (S-locus protein 11/S-locus cysteine-rich protein), is a small cysteine-rich protein, and the female S determinant, SRK (S-locus receptor kinase), functions as a receptor for SP11 at the surface of stigma papilla cells. Although a few of the following downstream factors in the SP11-SRK signaling cascade have been identified, a comprehensive understanding of the SI mechanism still remains unexplained in Brassicaceae. Analysis of self-compatible (SC) mutants is significant for understanding the molecular mechanism in SI reactions, thus we screened SC lines from a variety of Japanese bulk-populations of B. rapa vegetables. Two lines, TSC4 and TSC28, seem to have disruptions in the SI signaling cascade, while the other line, TSC2, seems to have a deficiency in a female S determinant, SRK. In TSC4 and TSC28, known SI-related factors, i.e. SRK, SP11, MLPK (M-locus protein kinase), THL (thioredoxin-h-like), and ARC1 (arm repeat containing 1), were expressed normally, and their expression levels were comparable with those in SI lines. On a B. rapa genetic linkage map, potential SC genes in TSC4 and TSC28 were mapped on linkage groups A3 and A1, respectively, whereas MLPK, ARC1, and THL were mapped on A3, A4, and A6, respectively. Although potential SC genes of TSC4 and MLPK were on the same linkage group, their positions were apparently independent. These results indicate that the SC genes of TSC4 and TSC28 are independent from the S-locus or known SI-related genes. Thus, the SC lines selected here have mutations in novel factors of the SI signaling cascade, and they will contribute to fill pieces in a signal transduction pathway of the SI system in Brassicaceae.

Cool temperature conditions are known to lead to pollen sterility in rice. Pollen sterility is an agriculturally important phenomenon because it imparts a large influence directly on rice yield. However, cool temperature stress tolerance varies among rice cultivars and avoidance of cool temperature stress is difficult by practical method of agriculture. In this study using two rice cultivars, Hitomebore (high tolerance) and Sasanishiki (low tolerance), we analyzed morphological features and gene expression profiles, under cool temperature stress, in anther development of rice. Hitomebore was given cool temperature stress (19 degrees C) at flowering stage, and showed 87.3% seed fertility. Meanwhile, the seed fertility decreased to 41.7% in the case of Sasanishiki. A transverse section of Hitomebore anther revealed that the degradation of the tapetum started at the uninucleate microspore stage, and the tapetum had completely vanished at mature stage. The tapetum provides nutrients for pollen development, and its degradation occurs at a late stage in pollen development. In contrast, degradation of the tapetum did not occur at the uninucleate microspore stage in Sasanishiki, and the tapetum was clearly intact at mature stage, suggesting that tapetum degradation is critical for accurate pollen development and cool temperature tolerance correlated with the degree of tapetum degeneration. In gene expression analysis of anther, 356 genes that showed different expression levels between two cultivars at cool temperatures were found. These genes will lead to understanding the mechanism of cool temperature stress response in rice pollen development and the identification of genes involved in accurate tapetum degradation.

Endosperm texture is an important factor governing the end-product quality of cereals. The texture of wheat (Triticum aestivum L) endosperm is controlled by puroindoline a and b genes which are both absent in rice (Oryza sativa L). It has been reported that the endosperm texture of rice can be modified by puroindoline genes. The mechanism, however, by which puroindolines affect the ultrastructure of rice endosperm cells remains to be investigated. In this study, we observed the ultrastructure of endosperm cells and the morphology of isolated starch granules of the transgenic rice expressing the puroindoline b gene. SEM and TEM observations indicated that compound starch granules were embedded within the matrix material in non-transgenic rice, Nipponbare, whereas they were surrounded by spaces in the transgenic rice. The morphology and size of each starch granule were not different between nontransgenic and the transgenic rice. However, the transgenic rice flour showed smaller particle size, higher starch damage, and lower viscosity during gelatinization than that of non-transgenic rice. These results confirm that puroindoline b reduces the grain hardness in rice. Moreover, the results also suggest that puroindoline b functions at the surface of compound starch granules, and not on polygonal starch granules in rice endosperm. (C) 2009 Elsevier Ltd. All rights reserved.

Introduction of T-DNAs with foreign genes by Agrobacterium-mediated transformation is widely used in plants. Multiple-introduced complex T-DNA loci, however, are difficult to clarify by conventional DNA gel blot analysis. We performed fluorescence in situ hybridization oil extended DNA fibers (fiber FISH) of transgenic tobacco plants harboring multiple 37-kb T-DNA constructs. Five and seven types of integrations were Successfully visualized in two transgenic lines. Most of the loci suffered duplication, deletion and/or translocation, indicating the complex integration events of the medium-size T-DNA. We Concluded that fiber-FISH analysis is a powerful tool to analyze organization Of multiple T-DNA loci in detail.

Transformation with large DNA molecules enables multiple genes to be introduced into plants simultaneously to produce transgenic plants with complex phenotypes. In this study, a large DNA fragment (ca. 100 kb) containing a set of Aegilops tauschii hardness genes was introduced into rice plants using a novel transformation method, called bioactive beads-mediated transformation. Nine transgenic rice plants were obtained and the presence of transgenes in the rice genome was confirmed by PCR and FISH analyses. The results suggested that multiple transgenes were successfully integrated in all transgenic plants. The expression of one of the transgenes, puroindoline b, was confirmed at the mRNA and protein levels in the T(2) generation. Our study clearly demonstrates that the bioactive bead method is capable of producing transgenic rice plants carrying large DNA fragments. This method will facilitate the production of useful transgenic plants by introducing multiple genes simultaneously.

Sequencing of the onion (Allium cepa) genome is challenging because it has one of the largest nuclear genomes among cultivated plants. We undertook pilot sequencing of onion genomic DNA to estimate gene densities and investigate the nature and distribution of repetitive DNAs. Complete sequences from two onion BACs were AT rich (64.8%) and revealed long tracts of degenerated retroviral elements and transposons, similar to other larger plant genomes. Random BACs were end sequenced and only 3 of 460 ends showed significant (e &lt; -25) non-organellar hits to the protein databases. The BAC-end sequences were AT rich (63.4%), similar to the completely sequenced BACs. A total of 499,997 bp of onion genomic DNA yielded an estimated mean density of one gene per 168 kb, among the lowest reported to date. Methyl filtration was highly effective relative to random shotgun reads in reducing frequencies of anonymous sequences from 82 to 55% and increasing non-organellar protein hits from 4 to 42%. Our results revealed no evidence for gene-dense regions and indicated that sequencing of methyl-filtered genomic fragments should be an efficient approach to reveal genic sequences in the onion genome.

In flowering plants, the male gametophyte, the pollen, develops in the anther. Complex patterns of gene expression in both the gametophytic and sporophytic tissues of the anther regulate this process. The gene expression profiles of the microspore/pollen and the sporophytic tapetum are of particular interest. In this study, a microarray technique combined with laser microdissection (44K LM-microarray) was developed and used to characterize separately the transcriptomes of the microspore/pollen and tapetum in rice. Expression profiles of 11 known tapetum specific-genes were consistent with previous reports. Based on their spatial and temporal expression patterns, 140 genes which had been previously defined as anther specific were further classified as male gametophyte specific (71 genes, 51%), tapetum-specific (seven genes, 5%) or expressed in both male gametophyte and tapetum (62 genes, 44%). These results indicate that the 44K LM-microarray is a reliable tool to analyze the gene expression profiles of two important cell types in the anther, the microspore/pollen and tapetum.

The male gametophyte and tapetum play different roles during anther development although they are differentiated from the same cell lineage, the L2 layer. Until now, it has not been possible to delineate their transcriptomes due to technical difficulties in separating the two cell types. In the present study, we characterized the separated transcriptomes of the rice microspore/pollen and tapetum using laser microdissection (LM)-mediated microarray. Spatiotemporal expression patterns of 28,141 anther-expressed genes were classified into 20 clusters, which contained 3,468 (12.3%) anther-enriched genes. In some clusters, synchronous gene expression in the microspore and tapetum at the same developmental stage was observed as a novel characteristic of the anther transcriptome. Noteworthy expression patterns are discussed in connection with gene ontology (GO) categories and gene annotations, which are related to important biological events in anther development, such as pollen maturation, pollen germination, pollen tube elongation and pollen wall formation.

Small RNAs including microRNA (miRNA) and small interfering RNA (siRNA) are known as repressors of gene expression. There are many plant proteins involved in small RNA-mediated gene silencing, such as Dicer ribonucleases and RNA-dependent RNA polymerases. However, most of these proteins have been reported to be absent in the late developmental stage of the plant male gamete, pollen. In order to clarify the existence of the small RNAs during maturation of pollen, we cloned and sequenced small RNAs from rice anthers including tricellular pollen. From fifty six candidates of small RNAs, we identified two known miRNAs (miR166 and miR167), eight potential miRNAs, and ten putative heterochromatic siRNAs (hc-siRNAs). RNA gel blot analyses clearly showed that miR166 and miR167 were accumulated in the uninuclear pollen stage of anther development and remained until the tricellular pollen stage. Our cloning and RNA gel blot analyses of small RNAs led us to propose a possible function of small RNA-mediated gene regulation for the development of male gametes in rice.