Microarray expression analysis of meiosis and microsporogenesis in hexaploid bread wheat

Pericarp starch metabolism is associated with caryopsis development and endosperm starch accumulation in common wheat

The wheat pericarp is the main component of the caryopsis at the early development stage and ultimately converts into a tissue that covers the mature caryopsis. A large number of starch granules are accumulated in the pericarp, but the production of and the role of starch granules in caryopsis development remain- elusive. In the present study, the relationship between accumulated starch granules and starch metabolism-related genes in wheat pericarp was investigated using paraffin section observations, expression analysis, and mutant analysis. Starch synthesis is initiated before anthesis and is dependent on a sucrose uptake and conversion system similar to that in the endosperm. TaPTST2 is required to initiate the production of pericarp starch granules. Pericarp starch granules gradually disappeared at the filling stage with high expression levels of genes encoding β-amylase, sucrose-phosphate synthase, and sucrose-phosphate phosphatase. As a maternal tissue adjacent to the endosperm and embryo, the pericarp plays a temporary reservoir for excess nutrients delivered into the caryopsis during the early development stage and exported at the filling stage. The pericarp contributes to the development of the endosperm and embryo as well as the accumulation of endosperm starch. The metabolism of pericarp starch may affect the weight of the wheat caryopsis.

The coordinated regulation of early meiotic stages is dominated by non-coding RNAs and stage-specific transcription in wheat

Reproductive success hinges on precisely coordinated meiosis, yet our understanding of how structural rearrangements of chromatin and phase transitions during meiosis are transcriptionally regulated is limited. In crop plants, detailed analysis of the meiotic transcriptome could identify regulatory genes and epigenetic regulators that can be targeted to increase recombination rates and broaden genetic variation, as well as provide a resource for comparison among eukaryotes of different taxa to answer outstanding questions about meiosis. We conducted a meiotic stage-specific analysis of messenger RNA (mRNA), small non-coding RNA (sncRNA), and long intervening/intergenic non-coding RNA (lincRNA) in wheat (Triticum aestivum L.) and revealed novel mechanisms of meiotic transcriptional regulation and meiosis-specific transcripts. Amidst general repression of mRNA expression, significant enrichment of ncRNAs was identified during prophase I relative to vegetative cells. The core meiotic transcriptome was comprised of 9309 meiosis-specific transcripts, 48 134 previously unannotated meiotic transcripts, and many known and novel ncRNAs differentially expressed at specific stages. The abundant meiotic sncRNAs controlled the reprogramming of central metabolic pathways by targeting genes involved in photosynthesis, glycolysis, hormone biosynthesis, and cellular homeostasis, and lincRNAs enhanced the expression of nearby genes. Alternative splicing was not evident in this polyploid species, but isoforms were switched at phase transitions. The novel, stage-specific regulatory controls uncovered here challenge the conventional understanding of this crucial biological process and provide a new resource of requisite knowledge for those aiming to directly modulate meiosis to improve crop plants. The wheat meiosis transcriptome dataset can be queried for genes of interest using an eFP browser located at https://bar.utoronto.ca/efp_wheat/cgi-bin/efpWeb.cgi?dataSource=Wheat_Meiosis.

New insight into the sucrose biosynthesis pathway from genome-wide identification, gene expression analysis, and subcellular localization in hexaploid wheat (Triticum aestivum L.)

Sucrose, the main synthesized product and transported form of photoassimilates, moves from leaves to support plant growth and storage substance biosynthesis occurring in the heterotrophic sink organs. Enhancing sucrose biosynthesis efficiency is a top priority for crop yield breeding programs. However, the molecular mechanism of sucrose biosynthesis is still elusive, especially in wheat. We performed transcriptome sequencing, subcellular localization, and bioinformatics analysis to identify key isoforms and metabolic branches associated with sucrose biosynthesis in wheat. Our analysis identified 45 genes from 13 families that exhibited high expression in young leaves with an evident diurnal change. The carbon flux from photoassimilates to sucrose was divided into two pathways. In the cytoplasm, assimilates initiating at phosphotrioses (TPs) exported by TaTPT1 from chloroplasts flowed along the TP-Sucrose branch formed by TaALD6, TaFBP5, TacPGI, TacPGM, TaUGP1, TaSPS5, and TaSPP1. Intermediates either from the Calvin cycle or TP-Sucrose branch were converted into ADPGlc to synthesize the simple starch, which was transiently degraded by a series of enzymes, including TaBAM4 and TaSEX4 in the chloroplast. Similar to the starch-biosynthesis branch in endosperms, the TP-Sucrose branch is the most prominent in leaves because each reaction can be catalyzed by at least one highly expressed isoform with expected cytosolic localization. The key isoforms and major branches highlighted in the wheat sucrose biosynthesis pathway expand our molecular understanding of crop sucrose biosynthesis and provide clues to increase wheat yield by enhancing the sucrose synthesis efficiency of leaves.

Allelic Variation Analysis at the Vernalization Response and Photoperiod Genes in Russian Wheat Varieties Identified Two Novel Alleles of Vrn-B3

Heading time is an important agronomic trait affecting the adaptability and productivity of common wheat. In this study, 95 common wheat varieties from Russia and the late-maturing breeding line ‘Velut’ were tested for allelic diversity of genes having the strongest effect on heading. In this research, allelic variation at the Ppd-D1, Vrn-A1, Vrn-B1, Vrn-D1, and Vrn-B3 loci was tested. The Vrn-B1 and Vrn-B3 loci provided the largest contribution to genetic diversity. We found two novel allelic variants of the Vrn-B3 gene in the studied varieties. Ten varieties carried a 160 bp insertion in the promoter region, and the breeding line ‘Velut’ carried a 1617 bp insertion. These alleles were designated Vrn-B3e and Vrn-B3d, respectively. The analysis of the sequences showed the recent insertion of a retrotransposon homologous to the LTR retrotransposon (RLX_Hvul_Dacia_ RND-1) in the Vrn-B3d allele. Plants with the Vrn-B3e and the ‘Velut’ line with the Vrn-B3d allele headed later than the plants with the wild-type allele; among these plants, ‘Velut’ is the latest maturing wheat variety. Analysis of the gene expression of two groups of lines differing by the Vrn-B3 alleles (Vrn-B3d or vrn-B3) from the F₂ population with ‘Velut’ as a parental line did not reveal a significant difference in the expression level between the groups. Additional research is required to study the reasons for the late maturation of the ‘Velut’ line. However, the studied wheat varieties could be used as a potential source of natural variation in genes controlling heading times.

Occurrence of albinism during wheat androgenesis is correlated with repression of the key genes required for proper chloroplast biogenesis

The phenomenon of albinism in wheat androgenesis is linked to the transcriptional repression of specific genes involved in chloroplast biogenesis during the first weeks of in vitro culture. Isolated microspore culture is widely used to accelerate breeding programs and produce new cultivars. However, in cereals and particularly in wheat, the use of this technique is limited due to the high proportion of regenerated albino plantlets. The causes and mechanisms leading to the formation of albino plantlets in wheat remain largely unknown and, to date, no concrete solution has been found to make it possible to overcome this barrier. We performed a molecular study of proplastid-to-chloroplast differentiation within wheat microspore cultures by analyzing the expression of 20 genes specifically involved in chloroplast biogenesis. Their expression levels were compared between two wheat genotypes that exhibit differential capacities to regenerate green plantlets, i.e., Pavon and Paledor, which produce high and low rates of green plants, respectively. We observed that chloroplast biogenesis within wheat microspores was affected as of the very early stages of the androgenesis process. A successful transition from a NEP- to a PEP-dependent transcription during early plastid development was found to be strongly correlated with the formation of green plantlets, while failure of this transition was strongly correlated with the regeneration of albino plantlets. The very low expression of plastid-encoded 16S and 23S rRNAs within plastids of the recalcitrant genotype Paledor suggests a low translation activity in albino plastids. Furthermore, a delay in the activation of the transcription of nuclear encoded key genes like GLK1 related to chloroplast biogenesis was observed in multicellular structures and pro-embryos of the genotype Paledor. These data help to understand the phenomenon of albinism in wheat androgenesis, which appears to be linked to the transcriptional activation of specific genes involved in the initial steps of chloroplast biogenesis that occurs between days 7 and 21 of in vitro culture.

Identification and validation of reference genes for RT-qPCR normalization in wheat meiosis

Meiosis is a specialized type of cell division occurring in sexually reproducing organisms to generate haploid cells known as gametes. In flowering plants, male gametes are produced in anthers, being encased in pollen grains. Understanding the genetic regulation of meiosis key events such as chromosome recognition and pairing, synapsis and recombination, is needed to manipulate chromosome associations for breeding purposes, particularly in important cereal crops like wheat. Reverse transcription-quantitative PCR (RT-qPCR) is widely used to analyse gene expression and to validate the results obtained by other transcriptomic analyses, like RNA-seq. Selection and validation of appropriate reference genes for RT-qPCR normalization is essential to obtain reproducible and accurate expression data. In this work, twelve candidate reference genes were evaluated using the mainstream algorithms geNorm, Normfinder, BestKeeper and ΔCt, then ranked from most to least suitable for normalization with RefFinder. Different sets of reference genes were recommended to normalize gene expression data in anther meiosis of bread and durum wheat, their corresponding genotypes in the absence of the Ph1 locus and for comparative studies among wheat genotypes. Comparisons between meiotic (anthers) and somatic (leaves and roots) wheat tissues were also carried out. To the best of our knowledge, our study provides the first comprehensive list of reference genes for robust RT-qPCR normalization to study differentially expressed genes during male meiosis in wheat in a breeding framework.

Interaction between serine carboxypeptidase-like protein TtGS5 and Annexin D1 in developing seeds of Triticum timopheevi

GS5 encoding a serine carboxypeptidase-like protein positively regulates grain size and weight through the regulation of grain width and filling and is helpful in improving cereal yields. Grain width variation determined by GS5 is associated with cell number and size, but the actual underlying mechanism is still unclear. Two orthologs of GS5, TtGS5-3A-G and TtGS5-3G-G, were cloned from the Triticum timopheevi accession no. CWI17006. To identify the proteins that interacted with TtGS5-3A-G and TtGS5-3G-G in premature grains, we performed pull-down assays followed by liquid chromatography-mass spectrometry/mass spectrometry analysis. The analyses revealed 18 proteins were present in both the TtGS5-3A-G and TtGS5-3G-G interactomes. Among five candidates selected, only Annexin D1 interacted with both TtGS5-3A-G and TtGS5-3G-G in yeast. Annexin D1, TtGS5-3A-G, and TtGS5-3G-G were located on the cytoplasmic membranes of Arabidopsis protoplasts and onion epidermal cells, and interactions between Annexin D1 and TtGS5-3A-G, as well as TtGS5-3G-G, were shown by bimolecular fluorescence complementation assays. Annexin D1 was expressed widely in different tissues, and it co-expressed with TtGS5-3A-G/TtGS5-3G-G at the grain enlargement phase. These results indicated that Annexin D1 interacted with TtGS5-3A-G and TtGS5-3G-G in premature grains. Together with the structural similarities of Annexin D1 to known fiber elongation factors, we proposed that TtGS5 might regulate the cell size by interacting with Annexin D1. The results provide significant new information for understanding the roles that GS5 plays in regulating grain size, which may be useful in improving crop yields.

Tapetum-Dependent Male Meiosis Progression in Plants: Increasing Evidence Emerges

In higher plants, male meiosis is a key process during microsporogenesis and is crucial for male fertility and seed set. Meiosis involves a highly dynamic organization of chromosomes and cytoskeleton and specifically takes place within sexual cells. However, studies in multiple plant species have suggested that the normal development of tapetum, the somatic cell layer surrounding the developing male meiocytes, is indispensable for the completion of the male meiotic cell cycle. Disrupted tapetum development causes alterations in the expression of a large range of genes involved in male reproduction. Moreover, recent experiments suggest that small RNAs (sRNAs) present in the anthers, including microRNAs (miRNAs) and phased, secondary, small interfering RNAs (phasiRNAs), play a potential but important role in controlling male meiosis, either by influencing the expression of meiotic genes in the meiocytes or through other unclear mechanisms, supporting the hypothesis that male meiosis is non-cell autonomously regulated. In this mini review, we summarize the recorded meiotic defects that occur in plants with defective tapetum development in both Arabidopsis and crops. Thereafter, we outline the latest understanding on the molecular mechanisms that potentially underpin the tapetum-dependent regulation of male meiosis, and we especially discuss the regulatory role of sRNAs. At the end, we propose several outstanding questions that should be addressed in future studies.

The novel function of the Ph1 gene to differentiate homologs from homoeologs evolved in Triticum turgidum ssp. dicoccoides via a dramatic meiosis-specific increase in the expression of the 5B copy of the C-Ph1 gene

The Ph1 gene is the principal regulator of homoeologous chromosome pairing control (HECP) that ensures the diploid-like meiotic chromosome pairing behavior of polyploid wheat. The HECP control was speculated to have evolved after the first event of polyploidization. With the objective to accurately understand the evolution of the HECP control, wild emmer wheat accessions previously known to differ for HECP control were characterized for the structure and expression of the candidate Ph1 gene, C-Ph1. The C-TdPh1-5A and 5B gene copies of emmer wheat showed 98 and 99% DNA sequence similarity respectively with the corresponding hexaploid wheat copies. Further, the C-TdPh1-5B carried the C-Ph1-5B specific structural changes and transcribed three splice variants as observed in the hexaploid wheat. Further, single nucleotide changes differentiating accessions varying for HECP control were identified. Analyzed by quantitative expression analysis, the wild emmer accessions with HECP control showed ~ 10,000-fold higher transcript abundance of the C-TdPh1-5B copy during prophase-I compared to accessions lacking the control. Differential transcriptional regulation of C-TdPh1-5B splice variants further revealed that C-Ph1-5B^alt1 variant is mainly responsible for differential accumulation of C-Ph1-5B copy in accessions with HECP control. Taken together, these results showed that the HECP control evolved via transcriptional regulation of splice variants during meiosis.

Identification of conserved genes involved in nitrogen metabolic activities in wheat

Nitrogen (N) plays a very important role in crop growth and development. Many N-metabolism-related genes responsive to N application have been identified in many plants such as Arabidopsis, rice and maize; however, few genes have been reported in wheat, which is one of the most widely grown crops in the world. In this study, a wheat wild type with N dependent lesion mimic (LM) and its mutants without LM were used to identify conserved N-metabolism-related genes. TaPAP, TaUPS and TaNMR were differentially expressed among N levels both in the wild type and two of its mutants, and the expression patterns of these genes were further studied under application of three chemotypes of N (NH_4⁺, NO_3^- and NH₄NO₃). The results showed that these genes are conserved N-metabolism-related genes and TaNMR is a novel player in N-metabolism.

Profiling of indole metabolic pathway in thermo-sensitive Bainong male sterile line in wheat (Triticum aestivum L.)

Bainong male sterile (BNS) wheat (Triticum aestivum L.) is a thermo-sensitive genic male sterile line with excellent sterility and self-restoration. We focused on transcriptional profiles of differentially expressed probes between BNS sterile and fertile anthers. Anthers, rachis and spikes from sterile line and fertile line were collected. Extracted RNA was assayed using wheat expression microarray and Gene Ontology was analyzed using Cytoscape with ClueGO. An indole (indole-3-acetic acid: IAA) metabolism pathway sub-network was almost formed in all differentially expressed profiles between sterile and fertile samples. IAA sub-network contained four nodes of indole and alkaloid metabolism connecting main network via indole compounds. This sub-network was absent in rachis and intact in transformed fertile anthers, which was the main differently expressed metabolism pathway in F₁ anthers with restorer genes. Alkaloid metabolism was absent in sterile anthers. Abnormal metabolism of IAA may be involved in BNS sterility. BNS transformation may be regulated by the production of IAA and alkaloid metabolism pathway, which favor the safe utilization of the sterile line in hybrid wheat production.

Validation of reference genes for accurate normalization of gene expression with quantitative real-time PCR in Haloxylon ammodendron under different abiotic stresses

Haloxylon ammodendron plays an important role in maintaining the structure and function of the entire ecosystem where it grows. No suitable reference genes have been reported in H. ammodendron plants to date. In this study, a total of 8 reference genes (18S, ACT1, ACT7, UBC18, TUA5, GAPDH, EF-1α and UBQ10) were selected from the available trancriptome database, and the expression stability of these 8 candidate genes was validated under different abiotic stress with three different statistical algorithms (geNorm, NormFinder and BestKeeper). The results produced from different models were in agreement with each other essentially: 18S and TUA5 were the most stable genes under drought stress, 18S, the most stable gene under heat stress and mechanical damage, ACT7 and UBC18, stable under salt stress while TUA5 and GAPDH expressed constantly under mechanical damage, and ACT1 expressed steadily under cold conditions. Expression profiles of several stress response genes, including FT-5, FT-9, DREB2A and DREB2C, were further confirmed with various candidate reference genes. None of the candidate genes showed a constant expression among all tested samples. Hence, it's essential to use more than one reference gene in order to guarantee the accuracy of quantitative real-time PCR. The results of this study will contribute to the accuracy and reliability in transcripts quantification, which is of significance to transcription-based studies and applications in this important shrub H. ammodendron.

Identification of Wheat Inflorescence Development-Related Genes Using a Comparative Transcriptomics Approach

Inflorescence represents the highly specialized plant tissue producing the grains. Although key genes regulating flower initiation and development are conserved, the mechanism regulating fertility is still not well explained. To identify genes and gene network underlying inflorescence morphology and fertility of bread wheat, expressed sequence tags (ESTs) from different tissues were analyzed using a comparative transcriptomics approach. Based on statistical comparison of EST frequencies of individual genes in EST pools representing different tissues and verification with RT-PCR and RNA-seq data, 170 genes of 59 gene sets predominantly expressed in the inflorescence were obtained. Nearly one-third of the gene sets displayed differentiated expression profiles in terms of their subgenome orthologs. The identified genes, most of which were predominantly expressed in anthers, encode proteins involved in wheat floral identity determination, anther and pollen development, pollen-pistil interaction, and others. Particularly, 25 annotated gene sets are associated with pollen wall formation, of which 18 encode enzymes or proteins participating in lipid metabolic pathway, including fatty acid ω-hydroxylation, alkane and fatty alcohol biosynthesis, and glycerophospholipid metabolism. We showed that the comparative transcriptomics approach was effective in identifying genes for reproductive development and found that lipid metabolism was particularly active in wheat anthers.

Reference genes validation in Phenacoccus solenopsis under various biotic and abiotic stress conditions

Real-time PCR (RT-qPCR) expression analysis is a powerful analytical technique, but for normalization of data requires the use of stable reference genes. However, suitable reference genes are still not known in the case of Phenacoccus solenopsis under variable experimental treatments. The present study focused on the identification of stable housekeeping genes as a reference for analysis under different abiotic and biotic factors in P. solenopsis. We analyzed the relative expression of six commonly used candidate reference genes in different developmental stages, host-feeding assay, temperature treatments and field distribution conditions. Expression stabilities were analyzed by geNorm, NormFinder, and RefFinder. Under developmental and field distribution conditions, β-Tubulin was found to be most stable reference genes followed by rpl32 and α-Tubulin. In the case host feeding treatment conditions, β-Tubulin and α-tubulin identified to be the most stable reference genes, while in temperature stress, a combination of α-Tubulin and rpl32 found to be suitable for normalizing the RT-qPCR data. Further, the above-identified genes were validated using RT-qPCR based gene expression analysis of four objective genes namely, Myoinhibitory peptides (MIPs), Zinc_metalloprotease (Zn_Mp), fatty acid synthase (fas) and alpha-glucosidase. Identified reference genes will facilitate gene expression studies in future under different stress treatments in P. solenopsis.

Analysis of the meiotic transcriptome reveals the genes related to the regulation of pollen abortion in cytoplasmic male-sterile pepper (Capsicum annuum L.)

CMS, which refers to the inability to generate functional pollen grains while still producing a normal gynoecium, has been widely used for pepper hybrid seed production. Pepper line 8214A is an excellent CMS line exhibiting 100% male sterility and superior economic characteristics. A TUNEL assay revealed the nuclear DNA is damaged in 8214A PMCs during meiosis. TEM images indicated that the 8214A PMCs exhibited asynchronous meiosis after prophase I, and some PMCs degraded prematurely with morphological features typical of PCD. Additionally, at the end of meiosis, the 8214A PMCs formed abnormal non-tetrahedral tetrads that degraded in situ. To identify the genes involved in the pollen abortion of line 8214A, the transcriptional profiles of the 8214A and the 8214B anthers (i.e., from the fertile maintainer line) during meiosis were analyzed using an RNA-seq approach. A total of 1355 genes were determined to be differentially expressed, including 424 and 931 up- and down- regulated genes, respectively, in the 8214A anthers during meiosis relative to the expression levels in the 8214B. The expression levels of ubiquitin ligase and cell cycle-related genes were apparently down-regulated, while the expression of methyltransferase genes was up-regulated in the 8214A anthers during meiosis, which likely contributed to the PCD of these PMCs during meiosis. Thus, our results may be useful for revealing the molecular mechanism regulating the pollen abortion of CMS pepper.

Plant transcriptomics and responses to environmental stress: an overview

Different stresses include nutrient deficiency, pathogen attack, exposure to toxic chemicals etc. Transcriptomic studies have been mainly applied to only a few plant species including the model plant, Arabidopsis thaliana. These studies have provided valuable insights into the genetic networks of plant stress responses. Transcriptomics applied to cash crops including barley, rice, sugarcane, wheat and maize have further helped in understanding physiological and molecular responses in terms of genome sequence, gene regulation, gene differentiation, posttranscriptional modifications and gene splicing. On the other hand, comparative transcriptomics has provided more information about plant's response to diverse stresses. Thus, transcriptomics, together with other biotechnological approaches helps in development of stress tolerance in crops against the climate change.

Novel Meiotic miRNAs and Indications for a Role of PhasiRNAs in Meiosis

Small RNAs (sRNA) add additional layers to the regulation of gene expression, with siRNAs directing gene silencing at the DNA level by RdDM (RNA-directed DNA methylation), and micro RNAs (miRNAs) directing post-transcriptional regulation of specific target genes, mostly by mRNA cleavage. We used manually isolated male meiocytes from maize (Zea mays) to investigate sRNA and DNA methylation landscapes during zygotene, an early stage of meiosis during which steps of meiotic recombination and synapsis of paired homologous chromosomes take place. We discovered two novel miRNAs from meiocytes, zma-MIR11969 and zma-MIR11970, and identified putative target genes. Furthermore, we detected abundant phasiRNAs of 21 and 24 nt length. PhasiRNAs are phased small RNAs which occur in 21 or 24 nt intervals, at a few hundred loci, specifically in male reproductive tissues in grasses. So far, the function of phasiRNAs remained elusive. Data from isolated meiocytes now revealed elevated DNA methylation at phasiRNA loci, especially in the CHH context, suggesting a role for phasiRNAs in cis DNA methylation. In addition, we consider a role of these phasiRNAs in chromatin remodeling/dynamics during meiosis. However, this is not well supported yet and will need more additional data. Here, we only lay out the idea due to other relevant literature and our additional observation of a peculiar GC content pattern at phasiRNA loci. Chromatin remodeling is also indicated by the discovery that histone genes were enriched for sRNA of 22 nt length. Taken together, we gained clues that lead us to hypothesize sRNA-driven DNA methylation and possibly chromatin remodeling during male meiosis in the monocot maize which is in line with and extends previous knowledge.

Polyploidy Enhances F1 Pollen Sterility Loci Interactions That Increase Meiosis Abnormalities and Pollen Sterility in Autotetraploid Rice

Intersubspecific autotetraploid rice (Oryza sativa ssp. indica × japonica) hybrids have greater biological and yield potentials than diploid rice. However, the low fertility of intersubspecific autotetraploid hybrids, which is largely caused by high pollen abortion rates, limits their commercial utility. To decipher the cytological and molecular mechanisms underlying allelic interactions in autotetraploid rice, we developed an autotetraploid rice hybrid that was heterozygous (S(i)S(j)) at F1 pollen sterility loci (Sa, Sb, and Sc) using near-isogenic lines. Cytological studies showed that the autotetraploid had higher percentages (>30%) of abnormal chromosome behavior and aberrant meiocytes (>50%) during meiosis than did the diploid rice hybrid control. Analysis of gene expression profiles revealed 1,888 genes that were differentially expressed between the autotetraploid and diploid hybrid lines at the meiotic stage, among which 889 and 999 were up- and down-regulated, respectively. Of the 999 down-regulated genes, 940 were associated with the combined effect of polyploidy and pollen sterility loci interactions (IPE). Gene Ontology enrichment analysis identified a prominent functional gene class consisting of seven genes related to photosystem I (Gene Ontology 0009522). Moreover, 55 meiosis-related or meiosis stage-specific genes were associated with IPE in autotetraploid rice, including Os02g0497500, which encodes a DNA repair-recombination protein, and Os02g0490000, which encodes a component of the ubiquitin-proteasome pathway. These results suggest that polyploidy enhances epistatic interactions between alleles of pollen sterility loci, thereby altering the expression profiles of important meiosis-related or meiosis stage-specific genes and resulting in high pollen sterility.

Meiosis, unreduced gametes, and parthenogenesis: implications for engineering clonal seed formation in crops

Meiosis and unreduced gametes. Sexual flowering plants produce meiotically derived cells that give rise to the male and female haploid gametophytic phase. In the ovule, usually a single precursor (the megaspore mother cell) undergoes meiosis to form four haploid megaspores; however, numerous mutants result in the formation of unreduced gametes, sometimes showing female specificity, a phenomenon reminiscent of the initiation of gametophytic apomixis. Here, we review the developmental events that occur during female meiosis and megasporogenesis at the light of current possibilities to engineer unreduced gamete formation. We also provide an overview of the current understanding of mechanisms leading to parthenogenesis and discuss some of the conceptual implications for attempting the induction of clonal seed production in cultivated plants.

The meiotic transcriptome architecture of plants

Although a number of genes that play key roles during the meiotic process have been characterized in great detail, the whole process of meiosis is still not completely unraveled. To gain insight into the bigger picture, large-scale approaches like RNA-seq and microarray can help to elucidate the transcriptome landscape during plant meiosis, discover co-regulated genes, enriched processes, and highly expressed known and unknown genes which might be important for meiosis. These high-throughput studies are gaining more and more popularity, but their beginnings in plant systems reach back as far as the 1960's. Frequently, whole anthers or post-meiotic pollen were investigated, while less data is available on isolated cells during meiosis, and only few studies addressed the transcriptome of female meiosis. For this review, we compiled meiotic transcriptome studies covering different plant species, and summarized and compared their key findings. Besides pointing to consistent as well as unique discoveries, we finally draw conclusions what can be learned from these studies so far and what should be addressed next.

The transcriptome landscape of early maize meiosis

BACKGROUND

A major step in the higher plant life cycle is the decision to leave the mitotic cell cycle and begin the progression through the meiotic cell cycle that leads to the formation of gametes. The molecular mechanisms that regulate this transition and early meiosis remain largely unknown. To gain insight into gene expression features during the initiation of meiotic recombination, we profiled early prophase I meiocytes from maize (Zea mays) using capillary collection to isolate meiocytes, followed by RNA-seq.

RESULTS

We detected ~2,000 genes as preferentially expressed during early meiotic prophase, most of them uncharacterized. Functional analysis uncovered the importance of several cellular processes in early meiosis. Processes significantly enriched in isolated meiocytes included proteolysis, protein targeting, chromatin modification and the regulation of redox homeostasis. The most significantly up-regulated processes in meiocytes were processes involved in carbohydrate metabolism. Consistent with this, many mitochondrial genes were up-regulated in meiocytes, including nuclear- and mitochondrial-encoded genes. The data were validated with real-time PCR and in situ hybridization and also used to generate a candidate maize homologue list of known meiotic genes from Arabidopsis.

CONCLUSIONS

Taken together, we present a high-resolution analysis of the transcriptome landscape in early meiosis of an important crop plant, providing support for choosing genes for detailed characterization of recombination initiation and regulation of early meiosis. Our data also reveal an important connection between meiotic processes and altered/increased energy production.

Transcriptome analysis during seed germination of elite Chinese bread wheat cultivar Jimai 20

BACKGROUND

Wheat seed germination directly affects wheat yield and quality. Although transcriptome and proteome analyses during seed germination have been reported in some crop plant species, dynamic transcriptome characterization during wheat seed germination has not been conducted. We performed the first comprehensive dynamic transcriptome analysis during different seed germination stages of elite Chinese bread wheat cultivar Jimai 20 using the Affymetrix Wheat Genome Array.

RESULTS

A total of 61,703 probe sets representing 51,411 transcripts were identified during the five seed germination stages of Jimai 20, of which 2,825 differential expression probe sets corresponding to 2,646 transcripts with different functions were declared by ANOVA and a randomized variance model. The seed germination process included a rapid initial uptake phase (0-12 hours after imbibition [HAI]), a plateau phase (12-24 HAI), and a further water uptake phase (24-48 HAI), corresponding to switches from the degradation of small-molecule sucrose to the metabolism of three major nutrients and to photosynthesis. Hierarchical cluster and MapMan analyses revealed changes in several significant metabolism pathways during seed germination as well as related functional groups. The signal pathway networks constructed with KEGG showed three important genes encoding the phosphofructokinase family protein, with fructose-1, 6-bisphosphatase, and UTP-glucose-1-phosphate uridylyltransferase located at the center, indicating their pivotal roles in the glycolytic pathway, gluconeogenesis, and glycogenesis, respectively. Several significant pathways were selected to establish a metabolic pathway network according to their degree value, which allowed us to find the pathways vital to seed germination. Furthermore, 51 genes involved in transport, signaling pathway, development, lipid metabolism, defense response, nitrogen metabolism, and transcription regulation were analyzed by gene co-expression network with a k-core algorithm to determine which play pivotal roles in germination. Twenty-three meaningful genes were found, and quantitative RT-PCR analysis validated the expression patterns of 12 significant genes.

CONCLUSIONS

Wheat seed germination comprises three distinct phases and includes complicated regulation networks involving a large number of genes. These genes belong to many functional groups, and their co-regulations guarantee regular germination. Our results provide new insight into metabolic changes during seed germination and interactions between some significant genes.

Regulation of meiotic gene expression in plants

With the recent advances in genomics and sequencing technologies, databases of transcriptomes representing many cellular processes have been assembled. Meiotic transcriptomes in plants have been studied in Arabidopsis thaliana, rice (Oryza sativa), wheat (Triticum aestivum), petunia (Petunia hybrida), sunflower (Helianthus annuus), and maize (Zea mays). Studies in all organisms, but particularly in plants, indicate that a very large number of genes are expressed during meiosis, though relatively few of them seem to be required for the completion of meiosis. In this review, we focus on gene expression at the RNA level and analyze the meiotic transcriptome datasets and explore expression patterns of known meiotic genes to elucidate how gene expression could be regulated during meiosis. We also discuss mechanisms, such as chromatin organization and non-coding RNAs that might be involved in the regulation of meiotic transcription patterns.

Transcriptomic landscape of prophase I sunflower male meiocytes

Meiosis is a form of specialized cell division that generates gametes, allowing recombination of alleles and halving the chromosome number. Arabidopsis and maize are the plant models that have been most extensively studied to determine the genes involved in meiosis. Here we present an RNA-seq study in which gene expression in male meiocytes isolated during prophase I was compared to that in somatic tissues of the sunflower HA89 line. We sampled more than 490 million gene tags from these libraries, assembled them de novo into a sunflower transcriptome. We obtained expression data for 36,304 sunflower genes, of which 19,574 (54%) were differentially expressed (DE) between meiocytes and somatic tissue. We also determined the functional categories and metabolic pathways that are DE in these libraries. As expected, we found large differences between the meiotic and somatic transcriptomes, which is in accordance with previous studies in Arabidopsis and maize. Furthermore, most of the previously implicated meiotic genes were abundantly and DE in meiocytes and a large repertoire of transcription factors (TF) and genes related to silencing are expressed in the sunflower meiocytes. We detected TFs which appear to be exclusively expressed in meiocytes. Our results allow for a better understanding of the conservation and differences in the meiotic transcriptome of plants.

Developmental transcriptional profiling reveals key insights into Triticeae reproductive development

Despite their importance, there remains a paucity of large-scale gene expression-based studies of reproductive development in species belonging to the Triticeae. As a first step to address this deficiency, a gene expression atlas of triticale reproductive development was generated using the 55K Affymetrix GeneChip(®) wheat genome array. The global transcriptional profiles of the anther/pollen, ovary and stigma were analyzed at concurrent developmental stages, and co-expressed as well as preferentially expressed genes were identified. Data analysis revealed both novel and conserved regulatory factors underlying Triticeae floral development and function. This comprehensive resource rests upon detailed gene annotations, and the expression profiles are readily accessible via a web browser.

Analyzing the meiotic transcriptome using isolated meiocytes of Arabidopsis thaliana

Improved transcriptome sequencing technologies (RNA-seq) have advanced our understanding of the tissue-specific transcriptome landscapes, including those of messenger RNAs, noncoding RNAs and small RNAs. However, transcriptome profiles of plant meiocytes remain challenging due to the lack of efficient methods to enrich meiocytes for the analysis of temporal and spatial gene expression patterns during meiosis. In this chapter, we describe a method to analyze the Arabidopsis meiotic transcriptome using isolated male meiocytes.

Preliminary characterisation of two early meiotic wheat proteins after identification through 2D gel electrophoresis proteomics

Various genetic-based approaches including mutant population screens, microarray analyses, cloning and transgenesis have broadened our knowledge of gene function during meiosis in plants. Nonetheless, these genetic tools are not without inherent limitations. One alternative approach to studying plant meiosis, especially in polyploids such as Triticum aestivum L. (bread wheat), is proteomics. However, protein-based approaches using proteomics have seldom been described, with only two attempts at studying early plant meiosis reported. Here, we report the investigation of early bread wheat meiosis using proteomics. Five differentially expressed protein spots were identified using 2D gel electrophoresis (2DGE) on protein extracts from four pooled stages of meiosis and three genotypes (Chinese Spring wild-type, ph1b and ph2a wheat mutant lines). Tandem mass spectrometry (MS/MS) identification of peptides from these protein spots led to the isolation and characterisation of the full-length clones of a wheat Speckle-type POZ protein, an SF21-like protein and HSP70, and a partial coding sequence of a hexose transporter. Significantly, the putative functions of the Speckle-type POZ protein and HSP70 were confirmed using in vitro DNA binding assays. Through the use of a 2DGE proteomics approach, we show that proteomics is a viable alternative to genetic-based approaches when studying meiosis in wheat. More significantly, we report a potential role for a Speckle-type POZ protein and a HSP70 in chromosome pairing during the early stages of meiosis in bread wheat.

Transcriptional profile analysis of E3 ligase and hormone-related genes expressed during wheat grain development

BACKGROUND

Wheat grains are an important source of food, stock feed and raw materials for industry, but current production levels cannot meet world needs. Elucidation of the molecular mechanisms underlying wheat grain development will contribute valuable information to improving wheat cultivation. One of the most important mechanisms implicated in plant developmental processes is the ubiquitin-proteasome system (UPS). Among the different roles of the UPS, it is clear that it is essential to hormone signaling. In particular, E3 ubiquitin ligases of the UPS have been shown to play critical roles in hormone perception and signal transduction.

RESULTS

A NimbleGen microarray containing 39,179 UniGenes was used to study the kinetics of gene expression during wheat grain development from the early stages of cell division to the mid-grain filling stage. By comparing 11 consecutive time-points, 9284 differentially expressed genes were identified and annotated during this study. A comparison of the temporal profiles of these genes revealed dynamic transcript accumulation profiles with major reprogramming events that occurred during the time intervals of 80-120 and 220-240°Cdays. The list of the genes expressed differentially during these transitions were identified and annotated. Emphasis was placed on E3 ligase and hormone-related genes. In total, 173 E3 ligase coding genes and 126 hormone-related genes were differentially expressed during the cell division and grain filling stages, with each family displaying a different expression profile.

CONCLUSIONS

The differential expression of genes involved in the UPS and plant hormone pathways suggests that phytohormones and UPS crosstalk might play a critical role in the wheat grain developmental process. Some E3 ligase and hormone-related genes seem to be up- or down-regulated during the early and late stages of the grain development.

The isolation and characterisation of the wheat molecular ZIPper I homologue, TaZYP1

BACKGROUND

The synaptonemal complex (SC) is a proteinaceous tripartite structure used to hold homologous chromosomes together during the early stages of meiosis. The yeast ZIP1 and its homologues in other species have previously been characterised as the transverse filament protein of the synaptonemal complex. Proper installation of ZYP1 along chromosomes has been shown to be dependent on the axial element-associated protein, ASY1 in Arabidopsis.

RESULTS

Here we report the isolation of the wheat (Triticum aestivum) ZYP1 (TaZYP1) and its expression profile (during and post-meiosis) in wild-type, the ph1b deletion mutant as well as in Taasy1 RNAi knock-down mutants. TaZYP1 has a putative DNA-binding S/TPXX motif in its C-terminal region and we provide evidence that TaZYP1 interacts non-preferentially with both single- and double-stranded DNA in vitro. 3-dimensional dual immunofluorescence localisation assays conducted with an antibody raised against TaZYP1 show that TaZYP1 interacts with chromatin during meiosis but does not co-localise to regions of chromatin where TaASY1 is present. The TaZYP1 signal lengthens into regions of chromatin where TaASY1 has been removed in wild-type but this appears delayed in the ph1b mutant. The localisation profile of TaZYP1 in four Taasy1 knock-down mutants is similar to wild-type but TaZYP1 signal intensity appears weaker and more diffused.

CONCLUSIONS

In contrast to previous studies performed on plant species where ZYP1 signal is sandwiched by ASY1 signal located on both axial elements of the SC, data from the 3-dimensional dual immunofluorescence localisation assays conducted in this study show that TaZYP1 signal only lengthens into regions of chromatin after TaASY1 signal is being unloaded. However, the observation that TaZYP1 loading appears delayed in both the ph1b and Taasy1 mutants suggests that TaASY1 may still be essential for TaZYP1 to play a role in SC formation during meiosis. These data further suggest that the temporal installation of ZYP1 onto pairing homologous chromosomes in wheat is different to that of other plant species and highlights the need to study this synaptonemal complex protein on a species to species basis.

Poor Homologous Synapsis 1 Interacts with Chromatin but Does Not Colocalise with ASYnapsis 1 during Early Meiosis in Bread Wheat

Chromosome pairing, synapsis, and DNA recombination are three key processes that occur during early meiosis. A previous study of Poor Homologous Synapsis 1 (PHS1) in maize suggested that PHS1 has a role in coordinating these three processes. Here we report the isolation of wheat (Triticum aestivum) PHS1 (TaPHS1), and its expression profile during and after meiosis. While the TaPHS1 protein has sequence similarity to other plant PHS1/PHS1-like proteins, it also possesses a unique region of oligopeptide repeat units. We show that TaPHS1 interacts with both single- and double-stranded DNA in vitro and provide evidence of the protein region that imparts the DNA-binding ability. Immunolocalisation data from assays conducted using antisera raised against TaPHS1 show that TaPHS1 associates with chromatin during early meiosis, with the signal persisting beyond chromosome synapsis. Furthermore, TaPHS1 does not appear to colocalise with the asynapsis protein (TaASY1) suggesting that these proteins are probably independently coordinated. Significantly, the data from the DNA-binding assays and 3-dimensional immunolocalisation of TaPHS1 during early meiosis indicates that TaPHS1 interacts with DNA, a function not previously observed in either the Arabidopsis or maize PHS1 homologues. As such, these results provide new insight into the function of PHS1 during early meiosis in bread wheat.

Comparative Transcriptomics Reveals 129 Transcripts That Are Temporally Regulated during Anther Development and Meiotic Progression in Both Bread Wheat (Triticum aestivum) and Rice (Oryza sativa)

Meiosis is a specialised type of cell division in sexually reproducing organisms that generates genetic diversity and prevents chromosome doubling in successive generations. The last decade has seen forward and reverse genetic approaches identifying many genes in the plant kingdom which highlight similarities and differences in the mechanics of meiosis between taxonomic kingdoms. We present here a high throughput in silico analysis, using bread wheat and rice, which has generated a list of 129 transcripts containing genes with meiotic roles and some which are currently unknown.

Gene expression profiling of Arabidopsis meiocytes

Meiosis is a special type of cell division present in all organisms that reproduce by sexual reproduction. It ensures the transition between the sporophytic and gametophytic state and allows gamete production through meiotic recombination and chromosome number reduction. In this paper, we describe a technique for the isolation of Arabidopsis thaliana male meiocytes. From this cellular material, it was then possible to develop large-scale transcriptome studies using CATMA microarrays and thus to obtain an overview of genes expressed during Arabidopsis meiosis. The expression profiles were studied with either stringent statistical criteria or by performing clustering. Both methods resulted in gene clusters enriched in meiosis-specific genes (from 14- to 55-fold). Analysis of these data provided a unique set of genes that will be pivotal to further analysis aimed at understanding the meiotic process.

Overexpression of TaNAC69 leads to enhanced transcript levels of stress up-regulated genes and dehydration tolerance in bread wheat

NAC proteins are plant-specific transcription factors and enriched with members involved in plant response to drought stress. In this study, we analyzed the expression profiles of TaNAC69 in bread wheat using Affymetrix Wheat Genome Array datasets and quantitative RT-PCR. TaNAC69 expression was positively associated with wheat responses to both abiotic and biotic stresses and was closely correlated with a number of stress up-regulated genes. The functional analyses of TaNAC69 in transgenic wheat showed that TaNAC69 driven by a barley drought-inducible HvDhn4s promoter led to marked drought-inducible overexpression of TaNAC69 in the leaves and roots of transgenic lines. The HvDhn4s:TaNAC69 transgenic lines produced more shoot biomass under combined mild salt stress and water-limitation conditions, had longer root and more root biomass under polyethylene glycol-induced dehydration. Analysis of transgenic lines with constitutive overexpression of TaNAC69 showed the enhanced expression levels of several stress up-regulated genes. DNA-binding assays revealed that TaNAC69 and its rice homolog (ONAC131) were capable of binding to the promoter elements of three rice genes (chitinase, ZIM, and glyoxalase I) and an Arabidopsis glyoxalase I family gene, which are homologs of TaNAC69 up-regulated stress genes. These data suggest that TaNAC69 is involved in regulating stress up-regulated genes and wheat adaptation to drought stress.

Analysis of anther transcriptomes to identify genes contributing to meiosis and male gametophyte development in rice

BACKGROUND

In flowering plants, the anther is the site of male gametophyte development. Two major events in the development of the male germline are meiosis and the asymmetric division in the male gametophyte that gives rise to the vegetative and generative cells, and the following mitotic division in the generative cell that produces two sperm cells. Anther transcriptomes have been analyzed in many plant species at progressive stages of development by using microarray and sequence-by synthesis-technologies to identify genes that regulate anther development. Here we report a comprehensive analysis of rice anther transcriptomes at four distinct stages, focusing on identifying regulatory components that contribute to male meiosis and germline development. Further, these transcriptomes have been compared with the transcriptomes of 10 stages of rice vegetative and seed development to identify genes that express specifically during anther development.

RESULTS

Transcriptome profiling of four stages of anther development in rice including pre-meiotic (PMA), meiotic (MA), anthers at single-celled (SCP) and tri-nucleate pollen (TPA) revealed about 22,000 genes expressing in at least one of the anther developmental stages, with the highest number in MA (18,090) and the lowest (15,465) in TPA. Comparison of these transcriptome profiles to an in-house generated microarray-based transcriptomics database comprising of 10 stages/tissues of vegetative as well as reproductive development in rice resulted in the identification of 1,000 genes specifically expressed in anther stages. From this sub-set, 453 genes were specific to TPA, while 78 and 184 genes were expressed specifically in MA and SCP, respectively. The expression pattern of selected genes has been validated using real time PCR and in situ hybridizations. Gene ontology and pathway analysis of stage-specific genes revealed that those encoding transcription factors and components of protein folding, sorting and degradation pathway genes dominated in MA, whereas in TPA, those coding for cell structure and signal transduction components were in abundance. Interestingly, about 50% of the genes with anther-specific expression have not been annotated so far.

CONCLUSIONS

Not only have we provided the transcriptome constituents of four landmark stages of anther development in rice but we have also identified genes that express exclusively in these stages. It is likely that many of these candidates may therefore contribute to specific aspects of anther and/or male gametophyte development in rice. In addition, the gene sets that have been produced will assist the plant reproductive community in building a deeper understanding of underlying regulatory networks and in selecting gene candidates for functional validation.

TaNF-YB3 is involved in the regulation of photosynthesis genes in Triticum aestivum

Nuclear factor Y (NF-Y) transcription factor is a heterotrimer comprised of three subunits: NF-YA, NF-YB and NF-YC. Each of the three subunits in plants is encoded by multiple genes with differential expression profiles, implying the functional specialisation of NF-Y subunit members in plants. In this study, we investigated the roles of NF-YB members in the light-mediated regulation of photosynthesis genes. We identified two NF-YB members from Triticum aestivum (TaNF-YB3 & 7) which were markedly upregulated by light in the leaves and seedling shoots using quantitative RT-PCR. A genome-wide coexpression analysis of multiple Affymetrix Wheat Genome Array datasets revealed that TaNF-YB3-coexpressed transcripts were highly enriched with the Gene Ontology term photosynthesis. Transgenic wheat lines constitutively overexpressing TaNF-YB3 had a significant increase in the leaf chlorophyll content, photosynthesis rate and early growth rate. Quantitative RT-PCR analysis showed that the expression levels of a number of TaNF-YB3-coexpressed transcripts were elevated in the transgenic wheat lines. The mRNA level of TaGluTR encoding glutamyl-tRNA reductase, which catalyses the rate-limiting step of the chlorophyll biosynthesis pathway, was significantly increased in the leaves of the transgenic wheat. Significant increases in the expression level in the transgenic plant leaves were also observed for four photosynthetic apparatus genes encoding chlorophyll a/b-binding proteins (Lhca4 and Lhcb4) and photosystem I reaction centre subunits (subunit K and subunit N), as well as for a gene coding for chloroplast ATP synthase γ subunit. These results indicate that TaNF-YB3 is involved in the positive regulation of a number of photosynthesis genes in wheat.

Meiosis-specific gene discovery in plants: RNA-Seq applied to isolated Arabidopsis male meiocytes

BACKGROUND

Meiosis is a critical process in the reproduction and life cycle of flowering plants in which homologous chromosomes pair, synapse, recombine and segregate. Understanding meiosis will not only advance our knowledge of the mechanisms of genetic recombination, but also has substantial applications in crop improvement. Despite the tremendous progress in the past decade in other model organisms (e.g., Saccharomyces cerevisiae and Drosophila melanogaster), the global identification of meiotic genes in flowering plants has remained a challenge due to the lack of efficient methods to collect pure meiocytes for analyzing the temporal and spatial gene expression patterns during meiosis, and for the sensitive identification and quantitation of novel genes.

RESULTS

A high-throughput approach to identify meiosis-specific genes by combining isolated meiocytes, RNA-Seq, bioinformatic and statistical analysis pipelines was developed. By analyzing the studied genes that have a meiosis function, a pipeline for identifying meiosis-specific genes has been defined. More than 1,000 genes that are specifically or preferentially expressed in meiocytes have been identified as candidate meiosis-specific genes. A group of 55 genes that have mitochondrial genome origins and a significant number of transposable element (TE) genes (1,036) were also found to have up-regulated expression levels in meiocytes.

CONCLUSION

These findings advance our understanding of meiotic genes, gene expression and regulation, especially the transcript profiles of MGI genes and TE genes, and provide a framework for functional analysis of genes in meiosis.

Histological and microarray analysis of the direct effect of water shortage alone or combined with heat on early grain development in wheat (Triticum aestivum)

Based on the in silico analysis of the representation of expressed sequence tags (ESTs) in wheat grain-related cDNA libraries, a specific 15k oligonucleotide microarray has been developed in order to monitor environmental stress-dependent gene expression changes in the wheat caryopses. Using this array, the effect of water withdrawal, with and without additional heat stress, has been investigated during the first five days of kernel development on two wheat cultivars differing in their drought sensitivity. Water shortage affected (more than twofold change) the expression of only 0.5% of the investigated genes. A parallel heat treatment increased the ratio of responding genes to 5-7% because of the temperature stress and/or the increased water deficit because of enhanced evaporation. It could be established that the two cultivars, differing in their long-term adaptation capabilities to drought, responded to the short and direct stress treatments on the same way. In response to the combined drought and heat treatment, the coordinately altered expression of genes coding for storage proteins, enzymes involved in sugar/starch metabolism, histone proteins, heat shock proteins, proteases, tonoplast aquaporins as well as several transcription factors has been observed. These gene expression changes were in agreement with histological data that demonstrated the accelerated development of the embryo as well as the endosperm.

The RAD51 and DMC1 homoeologous genes of bread wheat: cloning, molecular characterization and expression analysis

BACKGROUND

Meiotic recombination in eukaryotes requires two homologues of the E. coli RecA proteins: Rad51 and Dmc1. Both proteins play important roles in the binding of single stranded DNA, homology search, strand invasion and strand exchange. Meiotic recombination has been well studied in Arabidopsis, rice, maize and the orthologues of RAD51 and DMC1 have been characterized. However genetic analysis of the RAD51 and DMC1 genes in bread wheat has been hampered due to the absence of complete sequence information and because of the existence of multiple copies of each gene in the hexaploid wheat genome.

FINDINGS

In this study we have identified that TaRAD51 and TaDMC1 homoeologues are located on group 7 and group 5 chromosomes of hexaploid wheat, respectively. Comparative sequence analysis of cDNA derived from the TaRAD51 and TaDMC1 homoeologues revealed limited sequence divergence at both the nucleotide and the amino acid level. Indeed, comparisons between the predicted amino acid sequences of TaRAD51 and TaDMC1 and those of other eukaryotes reveal a high degree of evolutionary conservation. Despite the high degree of sequence conservation at the nucleotide level, genome-specific primers for cDNAs of TaRAD51 and TaDMC1 were developed to evaluate expression patterns of individual homoeologues during meiosis. QRT-PCR analysis showed that expression of the TaRAD51 and TaDMC1 cDNA homoeologues was largely restricted to meiotic tissue, with elevated levels observed during the stages of prophase I when meiotic recombination occurs. All three homoeologues of both strand-exchange proteins (TaRAD51 and TaDMC1) are expressed in wheat.

CONCLUSIONS

Bread wheat contains three expressed copies of each of the TaRAD51 and TaDMC1 homoeologues. While differences were detected between the three cDNA homoeologues of TaRAD51 as well as the three homoeologues of TaDMC1, it is unlikely that the predicted amino acid substitutions would have an effect on the protein structure, based on our three-dimensional structure prediction analyses. There are differences in the levels of expression of the three homoeologues of TaRAD51 and TaDMC1 as determined by QRT-PCR and if these differences are reflected at the protein level, bread wheat may be more dependent upon a particular homoeologue to achieve full fertility than all three equally.

Analysis of the Basidiomycete Coprinopsis cinerea reveals conservation of the core meiotic expression program over half a billion years of evolution

Coprinopsis cinerea (also known as Coprinus cinereus) is a multicellular basidiomycete mushroom particularly suited to the study of meiosis due to its synchronous meiotic development and prolonged prophase. We examined the 15-hour meiotic transcriptional program of C. cinerea, encompassing time points prior to haploid nuclear fusion though tetrad formation, using a 70-mer oligonucleotide microarray. As with other organisms, a large proportion (∼20%) of genes are differentially regulated during this developmental process, with successive waves of transcription apparent in nine transcriptional clusters, including one enriched for meiotic functions. C. cinerea and the fungi Saccharomyces cerevisiae and Schizosaccharomyces pombe diverged ∼500–900 million years ago, permitting a comparison of transcriptional programs across a broad evolutionary time scale. Previous studies of S. cerevisiae and S. pombe compared genes that were induced upon entry into meiosis; inclusion of C. cinerea data indicates that meiotic genes are more conserved in their patterns of induction across species than genes not known to be meiotic. In addition, we found that meiotic genes are significantly more conserved in their transcript profiles than genes not known to be meiotic, which indicates a remarkable conservation of the meiotic process across evolutionarily distant organisms. Overall, meiotic function genes are more conserved in both induction and transcript profile than genes not known to be meiotic. However, of 50 meiotic function genes that were co-induced in all three species, 41 transcript profiles were well-correlated in at least two of the three species, but only a single gene (rad50) exhibited coordinated induction and well-correlated transcript profiles in all three species, indicating that co-induction does not necessarily predict correlated expression or vice versa. Differences may reflect differences in meiotic mechanisms or new roles for paralogs. Similarities in induction, transcript profiles, or both, should contribute to gene discovery for orthologs without currently characterized meiotic roles.

Meiosis is the part of the sexual reproduction process in which the number of chromosomes in an organism is halved. This occurs in most plants, animals, and fungi; and many of the proteins involved are the same in the different organisms that have been studied. We wanted to ask whether the genes involved in the meiotic process are turned on and off at the same stages of meiosis in organisms that separated a long time ago. To do this we looked at three fungal species, Saccharomyces cerevisiae (baker's yeast), Schizosaccharomyces pombe (a very distantly related fungus of the same phylum), and Coprinopsis cinerea (a mushroom-forming fungus of a different phylum), which had a common ancestor 500–900 million years ago (in comparison, rats and mice separated ∼23 million years ago). We lined up meiotic stages and found that gene expression during the meiotic process was more conserved for meiotic genes than for non-meiotic genes, indicating ancient conservation of the meiotic process.

TaNF-YC11, one of the light-upregulated NF-YC members in Triticum aestivum, is co-regulated with photosynthesis-related genes

Nuclear factor Y (NF-Y) is a heterotrimeric transcription factor complex. Each of the NF-Y subunits (NF-YA, NF-YB and NF-YC) in plants is encoded by multiple genes. Quantitative RT-PCR analysis revealed that five wheat NF-YC members (TaNF-YC5, 8, 9, 11 and 12) were upregulated by light in both the leaf and seedling shoot. Co-expression analysis of Affymetrix wheat genome array datasets revealed that transcript levels of a large number of genes were consistently correlated with those of the TaNF-YC11 and TaNF-YC8 genes in three to four separate Affymetrix array datasets. TaNF-YC11-correlated transcripts were significantly enriched with the Gene Ontology term photosynthesis. Sequence analysis in the promoters of TaNF-YC11-correlated genes revealed the presence of putative NF-Y complex binding sites (CCAAT motifs). Quantitative RT-PCR analysis of a subset of potential TaNF-YC11 target genes showed that ten out of the 13 genes were also light-upregulated in both the leaf and seedling shoot and had significantly correlated expression profiles with TaNF-YC11. The potential target genes for TaNF-YC11 include subunit members from all four thylakoid membrane-bound complexes required for the conversion of solar energy into chemical energy and rate-limiting enzymes in the Calvin cycle. These data indicate that TaNF-YC11 is potentially involved in regulation of photosynthesis-related genes.

Characterization and comparative sequence analysis of the DNA mismatch repair MSH2 and MSH7 genes from tomato

DNA mismatch repair proteins play an essential role in maintaining genomic integrity during replication and genetic recombination. We successfully isolated a full length MSH2 and partial MSH7 cDNAs from tomato, based on sequence similarity between MutS and plant MSH homologues. Semi-quantitative RT-PCR reveals higher levels of mRNA expression of both genes in young leaves and floral buds. Genetic mapping placed MSH2 and MSH7 on chromosomes 6 and 7, respectively, and indicates that these genes exist as single copies in the tomato genome. Analysis of protein sequences and phylogeny of the plant MSH gene family show that these proteins are evolutionarily conserved, and follow the classical model of asymmetric protein evolution. Genetic manipulation of the expression of these MSH genes in tomato will provide a potentially useful tool for modifying genetic recombination and hybrid fertility between wide crosses.

Selection of reference genes for quantitative real-time PCR expression studies in the apomictic and sexual grass Brachiaria brizantha

BACKGROUND

Brachiaria brizantha is an important forage grass. The occurrence of both apomictic and sexual reproduction within Brachiaria makes it an interesting system for understanding the molecular pathways involved in both modes of reproduction. Quantitative real time PCR (qRT-PCR) has emerged as an important technique to compare expression profile of target genes and, in order to obtain reliable results, it is important to have suitable reference genes. In this work, we evaluated eight potential reference genes for B. brizantha qRT-PCR experiments, isolated from cDNA ovary libraries. Vegetative and reproductive tissues of apomictic and sexual B. brizantha were tested to validate the reference genes, including the female gametophyte, where differences in the expression profile between sexual and apomictic plants must occur.

RESULTS

Eight genes were selected from a cDNA library of ovaries of B. brizantha considering the similarity to reference genes: EF1 (elongation factor 1 alpha), E1F4A (eukaryotic initiation factor 4A), GAPDH (glucose-6-phosphate dehydrogenase), GDP (glyceroldehyde-3-phosphate dehydrogenase), SUCOA (succinyl-CoA ligase), TUB (tubulin), UBCE (ubiquitin conjugating enzyme), UBI (ubiquitin). For the analysis, total RNA was extracted from 22 samples and raw Ct data after qRT-PCR reaction was analyzed for primer efficiency and for an overall analysis of Ct range among the different samples. Elongation factor 1 alpha showed the highest expression levels, whereas succinyl-CoA ligase showed the lowest within the chosen set of samples. GeNorm application was used for evaluation of the best reference genes, and according to that, the least stable genes, with the highest M values were tubulin and succinyl-CoA ligase and the most stable ones, with the lowest M values were elongation factor 1 alpha and ubiquitin conjugating enzyme, when both reproductive and vegetative samples were tested. For ovaries and spikelets of both sexual and apomictic B. brizantha the genes with the lowest M values were BbrizUBCE, BbrizE1F4A and BbrizEF1.

CONCLUSION

In total, eight genes belonging to different cellular processes were tested. Out of them, BbrizTUB was the less stable while BbrizEF1 followed by BbrizUBCE were the more stable genes considering male and female reproductive tissues, spikelets, roots and leaves. Regarding the best reference genes for ovary tissues, where apomictic and sexual reproduction must occur, the best reference genes were BbrizUBCE, BbrizE1F4A and BbrizEF1. Our results provide crucial information for transcriptional analysis in the Brachiaria ssp, helping to improve the quality of gene expression data in these species, which constitute an excellent plant system for the study of apomixis.

Understanding meiosis and the implications for crop improvement

Over the past 50 years, the understanding of meiosis has aged like a fine bottle of wine: the complexity is developing but the wine itself is still young. While emphasis in the plant kingdom has been placed on the model diploids Arabidopsis (Arabidopsis thaliana L.) and rice (Orzya sativa L.), our research has mainly focussed on the polyploid, bread wheat (Triticum aestivum L.). Bread wheat is an important food source for nearly two-thirds of the world's population. While creating new varieties can be achieved using existing or advanced breeding lines, we would also like to introduce beneficial traits from wild related species. However, expanding the use of non-adapted and wild germplasm in cereal breeding programs will depend on the ability to manipulate the cellular process of meiosis. Three important and tightly-regulated events that occur during early meiosis are chromosome pairing, synapsis and recombination. Which key genes control these events in meiosis (and how they do so) remains to be completely answered, particularly in crops such as wheat. Although the majority of published findings are from model organisms including yeast (Saccharomyces cerevisiae) and the nematode Caenorhabditis elegans, information from the plant kingdom has continued to grow in the past decade at a steady rate. It is with this new knowledge that we ask how meiosis will contribute to the future of cereal breeding. Indeed, how has it already shaped cereal breeding as we know it today?

Comparative transcriptomics in the Triticeae

Background

Barley and particularly wheat are two grass species of immense agricultural importance. In spite of polyploidization events within the latter, studies have shown that genotypically and phenotypically these species are very closely related and, indeed, fertile hybrids can be created by interbreeding. The advent of two genome-scale Affymetrix GeneChips now allows studies of the comparison of their transcriptomes.

Results

We have used the Wheat GeneChip to create a "gene expression atlas" for the wheat transcriptome (cv. Chinese Spring). For this, we chose mRNA from a range of tissues and developmental stages closely mirroring a comparable study carried out for barley (cv. Morex) using the Barley1 GeneChip. This, together with large-scale clustering of the probesets from the two GeneChips into "homologous groups", has allowed us to perform a genomic-scale comparative study of expression patterns in these two species. We explore the influence of the polyploidy of wheat on the results obtained with the Wheat GeneChip and quantify the correlation between conservation in gene sequence and gene expression in wheat and barley. In addition, we show how the conservation of expression patterns can be used to elucidate, probeset by probeset, the reliability of the Wheat GeneChip.

Conclusion

While there are many differences in expression on the level of individual genes and tissues, we demonstrate that the wheat and barley transcriptomes appear highly correlated. This finding is significant not only because given small evolutionary distance between the two species it is widely expected, but also because it demonstrates that it is possible to use the two GeneChips for comparative studies. This is the case even though their probeset composition reflects rather different design principles as well as, of course, the present incomplete knowledge of the gene content of the two species. We also show that, in general, the Wheat GeneChip is not able to distinguish contributions from individual homoeologs. Furthermore, the comparison between the two species leads us to conclude that the conservation of both gene sequence as well as gene expression is positively correlated with absolute expression levels, presumably reflecting increased selection pressure on genes coding for proteins present at high levels. In addition, the results indicate the presence of a correlation between sequence and expression conservation within the Triticeae.

Identification and validation of reference genes for quantitative RT-PCR normalization in wheat

BACKGROUND

Usually the reference genes used in gene expression analysis have been chosen for their known or suspected housekeeping roles, however the variation observed in most of them hinders their effective use. The assessed lack of validated reference genes emphasizes the importance of a systematic study for their identification. For selecting candidate reference genes we have developed a simple in silico method based on the data publicly available in the wheat databases Unigene and TIGR.

RESULTS

The expression stability of 32 genes was assessed by qRT-PCR using a set of cDNAs from 24 different plant samples, which included different tissues, developmental stages and temperature stresses. The selected sequences included 12 well-known HKGs representing different functional classes and 20 genes novel with reference to the normalization issue. The expression stability of the 32 candidate genes was tested by the computer programs geNorm and NormFinder using five different data-sets. Some discrepancies were detected in the ranking of the candidate reference genes, but there was substantial agreement between the groups of genes with the most and least stable expression. Three new identified reference genes appear more effective than the well-known and frequently used HKGs to normalize gene expression in wheat. Finally, the expression study of a gene encoding a PDI-like protein showed that its correct evaluation relies on the adoption of suitable normalization genes and can be negatively affected by the use of traditional HKGs with unstable expression, such as actin and alpha-tubulin.

CONCLUSION

The present research represents the first wide screening aimed to the identification of reference genes and of the corresponding primer pairs specifically designed for gene expression studies in wheat, in particular for qRT-PCR analyses. Several of the new identified reference genes outperformed the traditional HKGs in terms of expression stability under all the tested conditions. The new reference genes will enable more accurate normalization and quantification of gene expression in wheat and will be helpful for designing primer pairs targeting orthologous genes in other plant species.

TaASY1 promotes homologous chromosome interactions and is affected by deletion of Ph1

During meiosis, chromosomes are sorted into homologous pairs as a preface to their intimate association via recombination and synapsis. However, little is known about the mechanism used to distinguish homologous chromosomes from other chromosomes present in the nucleus. Studies in wheat (Triticum aestivum) have shown that the Pairing homoeologous 1 (Ph1) locus is required to suppress interactions between genetically similar homoeologous chromosomes. Here we show that absence of Ph1 causes increased transcription of Asynapsis 1 (ASY1), a gene that encodes an axial-element-associated protein that is essential for synapsis and cross-over formation in Arabidopsis and rice. Localisation of ASY1 during meiosis is also affected by deletion of Ph1. In addition, transgenic wheat mutants with decreased activity of TaASY1 display reduced synapsis during prophase I and exhibit pairing between homoeologous chromosomes at metaphase I. These results suggest that ASY1 is required to promote interactions between homologous chromosomes in bread wheat, and that Ph1 has a gene regulatory role, which is consistent with its suggested genetic identity as a Cdk-like gene. Broader implications of this research suggest that we could use the Taasy1 mutants to assess their efficacy in alien chromatin introgression studies, as seen with the ph1b mutant.

The RAD51 gene family in bread wheat is highly conserved across eukaryotes, with RAD51A upregulated during early meiosis

The RADiation sensitive protein 51 (RAD51) recombinase is a eukaryotic homologue of the bacterial Recombinase A (RecA). It is required for homologous recombination of DNA during meiosis where it plays a role in processes such as homology searching and strand invasion. RAD51 is well conserved in eukaryotes with as many as four paralogues identified in vertebrates and some higher plants. Here we report the isolation and preliminary characterisation of four RAD51 gene family members in hexaploid (bread) wheat (Triticum aestivum L.). RAD51A1, RAD51A2 and RAD51D were located on chromosome group 7, and RAD51C was on chromosome group 2. Q-PCR gene expression profiling revealed that RAD51A1 was upregulated during meiosis with lower expression levels seen in mitotic tissue, and bioinformatics analysis demonstrated the evolutionary linkages of this gene family to other eukaryotic RAD51 sequences. Western blot analysis of heterologously expressed RAD51 from bread wheat has shown that it is detectable using anti-human RAD51 antibodies and that molecular modelling of the same protein revealed structural conservation when compared with yeast, human, Arabidopsis and maize RAD51A orthologues. This report has widened the knowledge base of this important protein family in plants, and highlighted the high level of structural conservation among RAD51 proteins from various species.

Whole genome approaches to identify early meiotic gene candidates in cereals

Early events during meiotic prophase I underpin not only viability but the variation of a species from generation to generation. Understanding and manipulating processes such as chromosome pairing and recombination are integral for improving plant breeding. This study uses comparative genetics, quantitative trait locus (QTL) analysis and a transcriptomics-based approach to identify genes that might have a role in genome-wide recombination control. Comparative genetics and the analysis of the yeast and Arabidopsis sequenced genomes has allowed the identification of early meiotic candidates that are conserved in wheat, rice and barley. Secondly, scoring recombination frequency as a phenotype for QTL analysis across wheat, rice and barley mapping populations has enabled us to identify genomic regions and candidate genes that could be involved in genome-wide recombination. Transcriptome data for candidate genes indicate that they are expressed in meiotic tissues. Candidates identified included a non-annotated expressed protein, a DNA topoisomerase 2-like candidate, RecG, RuvB and RAD54 homologues.