Transcriptomic Analysis of Differentially Expressed Genes during Flower Organ Development in Genetic Male Sterile and Male Fertile Tagetes erecta by Digital Gene-Expression Profiling

Identification and characterization of CYC2-like genes related to floral symmetric development in Tagetes erecta (Asteraceae)

Marigold (Tagetes erecta) is an annual herbaceous flower belonging to Asteraceae, whose capitulum is composed of bilateral symmetry ray florets on the outer periphery and radial symmetry disk florets on the inside. The flower symmetry evolution from radial symmetry to bilateral symmetry has changed the morphology, inflorescence architecture and function of florets among several lineages in Asteraceae. Several studies have identified that CYC2 genes in TCP transcription factor family are the key genes regulating the flower morphogenesis, such as corolla symmetry and stamen development. Here, seven TeCYC2 genes were cloned and phylogenetically grouped into the CYC2 branch of TCP transcription family. TeCYC2c and TeCYC2d were found to be expressed specifically in ray florets, TeCYC2b was strongly expressed in both ray and disk florets, TeCYC2g was significantly higher expressed in ray florets than in disk florets, while TeCYC2a, TeCYC2e1 and TeCYC2e2 were significantly expressed in disk florets, according to an examination of the expression profile. Among the ectopic expression lines of seven TeCYC2 genes in Arabidopsis thaliana, the flower symmetry of all transgenic lines was changed from radial symmetry to bilateral symmetry, and only the reproductive growth of TeCYC2c lines was affected. In TeCYC2c transgenic Arabidopsis, the pollen sac was difficult to crack, and the filaments were shorter than the pistils, resulting in a significant decrease in the seed setting rate. All TeCYC2 proteins were localized in the nucleus. Eight pairs of interactions between TeCYC2 proteins were validated by Y2H and BiFC assays, indicating the possibility of TeCYC2 proteins forming homodimers or heterodimers to improve functional specificity. Our findings verified the main regulatory role of TeCYC2c on the development of corollas and stamen in marigold, and analyzed the interaction network of the formation mechanism of floral symmetry in two florets, which provided more insights into the expansion of CYC2 genes in the evolution of Asteraceae inflorescence and contributed to elucidate the complex regulatory network, as well as the molecular breeding concerning flower form diversity in marigold.

Integrated Bulk Segregant Analysis, Fine Mapping, and Transcriptome Revealed QTLs and Candidate Genes Associated with Drought Adaptation in Wild Watermelon

Drought stress has detrimental effects on crop productivity worldwide. A strong root system is crucial for maintaining water and nutrients uptake under drought stress. Wild watermelons possess resilient roots with excellent drought adaptability. However, the genetic factors controlling this trait remain uninvestigated. In this study, we conducted a bulk segregant analysis (BSA) on an F2 population consisting of two watermelon genotypes, wild and domesticated, which differ in their lateral root development under drought conditions. We identified two quantitative trait loci (qNLR_Dr. Chr01 and qNLR_Dr. Chr02) associated with the lateral root response to drought. Furthermore, we determined that a small region (0.93 Mb in qNLR_Dr. Chr01) is closely linked to drought adaptation through quantitative trait loci (QTL) validation and fine mapping. Transcriptome analysis of the parent roots under drought stress revealed unique effects on numerous genes in the sensitive genotype but not in the tolerant genotype. By integrating BSA, fine mapping, and the transcriptome, we identified six genes, namely L-Ascorbate Oxidase (AO), Cellulose Synthase-Interactive Protein 1 (CSI1), Late Embryogenesis Abundant Protein (LEA), Zinc-Finger Homeodomain Protein 2 (ZHD2), Pericycle Factor Type-A 5 (PFA5), and bZIP transcription factor 53-like (bZIP53-like), that might be involved in the drought adaptation. Our findings provide valuable QTLs and genes for marker-assisted selection in improving water-use efficiency and drought tolerance in watermelon. They also lay the groundwork for the genetic manipulation of drought-adapting genes in watermelon and other Cucurbitacea species.

Upcoming progress of transcriptomics studies on plants: An overview

Transcriptome sequencing or RNA-Sequencing is a high-resolution, sensitive and high-throughput next-generation sequencing (NGS) approach used to study non-model plants and other organisms. In other words, it is an assembly of RNA transcripts from individual or whole samples of functional and developmental stages. RNA-Seq is a significant technique for identifying gene predictions and mining functional analysis that improves gene ontology understanding mechanisms of biological processes, molecular functions, and cellular components, but there is limited information available on this topic. Transcriptomics research on different types of plants can assist researchers to understand functional genes in better ways and regulatory processes to improve breeding selection and cultivation practices. In recent years, several advancements in RNA-Seq technology have been made for the characterization of the transcriptomes of distinct cell types in biological tissues in an efficient manner. RNA-Seq technologies are briefly introduced and examined in terms of their scientific applications. In a nutshell, it introduces all transcriptome sequencing and analysis techniques, as well as their applications in plant biology research. This review will focus on numerous existing and forthcoming strategies for improving transcriptome sequencing technologies for functional gene mining in various plants using RNA- Seq technology, based on the principles, development, and applications.

Multiomics reveals an essential role of long-distance translocation in regulating plant cadmium resistance and grain accumulation in allohexaploid wheat (Triticum aestivum)

Cadmium (Cd) is a highly toxic heavy metal that readily enters cereals, such as wheat, via the roots and is translocated to the shoots and grains, thereby posing high risks to human health. However, the vast and complex genome of allohexaploid wheat makes it challenging to understand Cd resistance and accumulation. In this study, a Cd-resistant cultivar of wheat, 'ZM1860', and a Cd-sensitive cultivar, 'ZM32', selected from a panel of 442 accessions, exhibited significantly different plant resistance and grain accumulation. We performed an integrated comparative analysis of the morpho-physiological traits, ionomic and phytohormone profiles, genomic variations, transcriptomic landscapes, and gene functionality in order to identify the mechanisms underlying these differences. Under Cd toxicity, 'ZM1860' outperformed 'ZM32', which showed more severe leaf chlorosis, poorer root architecture, higher accumulation of reactive oxygen species, and disordered phytohormone homeostasis. Ionomics showed that 'ZM32' had a higher root-to-shoot translocation coefficient of Cd and accumulated more Cd in the grains than 'ZM1860'. Whole-genome re-sequencing (WGS) and transcriptome sequencing identified numerous DNA variants and differentially expressed genes involved in abiotic stress responses and ion transport between the two genotypes. Combined ionomics, transcriptomics, and functional gene analysis identified the plasma membrane-localized heavy metal ATPase TaHMA2b-7A as a crucial Cd exporter regulating long-distance Cd translocation in wheat. WGS- and PCR-based analysis of sequence polymorphisms revealed a 25-bp InDel site in the promoter region of TaHMA2b-7A, and this was probably responsible for the differential expression. Our multiomics approach thus enabled the identification of a core transporter involved in long-distance Cd translocation in wheat, and it may provide an elite genetic resource for improving plant Cd resistance and reducing grain Cd accumulation in wheat and other cereal crops.

Comparative analysis of long noncoding RNA and mRNA expression provides insights into adaptation to hypoxia in Tibetan sheep

Tibetan sheep have lived on the Qinghai-Tibetan Plateau for thousands of years and have good adaptability to the hypoxic environment and strong disease resistance. However, the molecular mechanism by which Tibetan sheep adapt to this extreme environment, especially the role of genetic regulation, is still unknown. Emerging evidence suggests that long noncoding RNAs (lncRNAs) participate in the regulation of a diverse range of biological processes. To explore the potential lncRNAs involved in the adaptation to high-altitude hypoxia of Tibetan sheep, we analysed the expression profile of lncRNAs and mRNAs in the liver and lung tissues of sheep using comparative transcriptome analysis between four Tibetan sheep populations (high altitude) and one Hu sheep population (low altitude). The results showed a total of 7848 differentially expressed (DE) lncRNA transcripts, and 22,971 DE mRNA transcripts were detected by pairwise comparison. The expression patterns of selected mRNAs and lncRNAs were validated by qRT-PCR, and the results correlated well with the transcriptome data. Moreover, the functional annotation analysis based on the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases showed that DE mRNAs and the target genes of the lncRNAs were significantly enriched in organ morphogenesis, response to stimulus, haem binding, the immune system, arginine and proline metabolism, and fatty acid biosynthesis. The prediction of mRNA–mRNA and lncRNA–mRNA interaction networks further revealed transcripts potentially involved in adaptation to high-altitude hypoxia, and the hub genes DDX24, PDCD11, EIF4A3, NDUFA11, SART1, PRPF8 and TCONS_00306477, TCONS_00306029, TCONS_00139593, TCONS_00293272, and TCONS_00313398 were selected. Additionally, a set of target genes, PIK3R1, IGF1R, FZD6, IFNB2, ATF3, MB, CYP2B4, PSMD13, and TGFB1, were also identified as candidate genes associated with high-altitude hypoxia adaptation. In conclusion, a collection of novel expressed lncRNAs, a set of target genes and biological pathways known to be relevant for altitude adaptation were identified by comparative transcriptome analysis between Tibetan sheep and Hu sheep. Our results are the first to identify the characterization and expression profile of lncRNAs between Tibetan sheep and Hu sheep and provide insights into the genetic regulation mechanisms by which Tibetan sheep adapt to high-altitude hypoxic environments.

Functional conservation and divergence of SEPALLATA-like genes in the development of two-type florets in marigold

Marigold (Tagetes erecta), as one member of Asteraceae family, bears a typical capitulum with two morphologically distinct florets. The SEPALLATA genes are involved in regulating the floral meristem determinacy, organ identity, fruit maturation, seed formation, and plant architecture. Here, five SEP-like genes were cloned and identified from marigold. Sequence alignment and phylogenetic analysis demonstrated that TeSEP3-1, TeSEP3-2, and TeSEP3-3 proteins were grouped into SEP3 clade, and TeSEP1 and TeSEP4 proteins were clustered into SEP1/2/4 clade. Quantitative real-time PCR analysis revealed that TeSEP1 and TeSEP3-3 were broadly expressed in floral organs, and that TeSEP3-2 and TeSEP4 were mainly expressed in pappus and corollas, while TeSEP3-1 was mainly expressed in two inner whorls. Ectopic expression of TeSEP1, TeSEP3-2, TeSEP3-3, and TeSEP4 in arabidopsis and tobacco resulted in early flowering. However, overexpression of TeSEP3-1 in arabidopsis and tobacco caused no visible phenotypic changes. Notably, overexpression of TeSEP4 in tobacco decreased the number of petals and stamens. Overexpression of TeSEP1 in tobacco led to longer sepals and simpler inflorescence architecture. The comprehensive pairwise interaction analysis suggested that TeSEP proteins had a broad interaction with class A, C, D, E proteins to form dimers. The yeast three-hybrid analysis suggested that in ternary complexes, class B proteins interacted with TeSEP3 by forming heterodimer TePI-TeAP3-2. The regulatory network analysis of MADS-box genes in marigold further indicated that TeSEP proteins played a "glue" role in regulating floral organ development, implying functional conservation and divergence of MADS box genes in regulating two-type floret developments. This study provides an insight into the formation mechanism of floral organs of two-type florets, thus broadening our knowledge of the genetic basis of flower evolution.

Comparative transcriptome between male fertile and male sterile alfalfa (Medicago varia)

Male sterility is an important factor in improving crop quality and yield through heterosis breeding. In this study, we analyzed the transcriptomes of male fertile (MF) and male sterile (MS) alfalfa flower buds using the Illumina HiSeq™ 4000 platform. A total of 54.05 million clean reads were generated and assembled into 65,777 unigenes with an average length of 874 bp. The differentially expressed genes (DEGs) between the MF and MS flowers at three stages of pollen development were identified, and there were 3832, 5678 and 5925 DEGs respectively in stages 1, 2 and 3. GO and KEGG functional enrichment analysis revealed 12, 12, 6 and 12 key branch-point genes involved in circadian rhythm, transcription factors, pollen development and flavonoid biosynthesis. Our findings provide novel insights into the mechanism of male sterility in alfalfa.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01026-x.

Integrative analysis of transcriptomic and proteomic changes related to male sterility in Tagetes erecta

Male sterile and male fertile two-type lines are important in heterosis utilization and breeding in Tagetes erecta, but the genes and pathways involved in male sterility are poorly understood. To explore these topics, transcriptome data (by RNA-seq) and proteome data (by iTRAQ) were gathered from flower buds of the male sterile line 'MS2-2' and male fertile line 'MF2-2' and integrated for a better understanding of the underlying molecular mechanisms of male sterility in T. erecta. The RNA-seq procedure generated 285,139,740 clean reads and 63359 unigenes and 6640 differentially expressed genes (DEGs) were identified, of which 4136 were downregulated and 2504 were upregulated in 'MS2-2'. DEGs related to flower development, pollen development, pollen wall assembly, endogenous hormones and transcription factors were identified. The iTRAQ analysis identified 3950 proteins in total; 789 were differentially expressed proteins (381 upregulated, 408 downregulated), which were mainly annotated to the Ribosome, Carbon metabolism and Biosynthesis of amino acids pathways. An association analysis revealed strong correlation (r Pearson = 0.6019) between the transcriptomic and proteomic data, and 256 and 34 proteins showed the same and opposite expression patterns with regard to their transcripts, respectively. Pathways such as photosynthesis, fatty acid biosynthesis and phenylpropanoid biosynthesis which influence tapetum and pollen development in male sterile plants, were significantly enriched at the transcript and protein levels. Most genes involved in these pathways were downregulated in 'MS2-2'. The low expression of these genes or functional loss of proteins could be associated with flower development, pollen development and related to changes in fertility in T. erecta. This study provided transcriptomic and proteomic information for T. erecta that could illuminate the mechanism of male sterility.

High level of intraspecific divergence and low frequency of RNA editing in the chloroplast genome sequence of Tagetes erecta

Tagetes erecta L. is an important commercial and medicinal plant. In this study, we reported the complete chloroplast genome sequence of T. erecta. The genome has a circular structure of 152,076 bp containing a large single-copy region (LSC) of 83,914 bp, a small copy region (SSC) of 18,064 bp, and two inverted repeats (IR) of 25,049 bp by each. It harbors 111 unique genes, including 79 protein-coding genes, 4 ribosomal RNA genes, and 28 transfer RNA genes. A total of 41 microsatellite, 20 tandem, and 37 interspersed repeats were detected in the genome. The phylogenomic analysis shows that T. erecta is a single phylogenetic cluster. The complete chloroplast genome of T. erecta lays the foundation for the phylogenetic, evolutionary, and conservation studies of the genus Tagetes. Furthermore, the intergenic region of atpB-rbcL was variable among the species T. erecta. This suggests that this region might be a mutation hotspot and will be useful for phylogenetic study and the development of molecular markers. At last, we systematically identified the RNA editing sites in the chloroplast genome of T. erecta based on the transcriptome downloaded from the SRA database. This study identified the characteristics of the T. erecta chloroplast genome, SNPs, and RNA editing sites, which will facilitate species identification and phylogenetic analysis within T. erecta.

Transcriptome and digital gene expression analysis unravels the novel mechanism of early flowering in Angelica sinensis

Angelica sinensis (Oliv.) Diels is a widely used medicinal plant mainly originated in Gansu, China. Angelica sinensis is greatly demanded in the clinical practice of Chinese medicine due to its broad pharmacological activities of hematopoietic and anti-inflammatory properties. But, the percentage of early flowering in Angelica sinensis arrives to 20%~30%, which severely affects its quality and quantity. Here, transcriptome profiling and digital gene expression analysis were applied to study the mechanism of early flowering in Angelica sinensis. A total of 49,183,534 clean reads were obtained and assembled into 68,262 unigenes, and 49,477 unigenes (72.5%) could be annotated to a minimum of one database in the Nr, Nt, Swiss-Pro, GO, COG and KEGG. Taking the above transcriptome data as a reference, digital gene expression result showed that 5,094 genes expression level were significant changed during early flowering. These annotated genes offered much information promoting that the biosynthesis of secondary metabolites pathway, the hormone signal transduction pathway, and the transcription regulation system may be closely related to the early flowering phenomenon of Angelica sinensis. Further expression patterns of key genes contribute to early flowering were analyzed using quantitative real-time PCR. The transcriptome result offered important gene expression information about early flowering in Angelica sinensis.

Propensity of Tagetes erecta L., a Medicinal Plant Commonly Used in Diabetes Management, to Accumulate Perfluoroalkyl Substances

It has been extensively demonstrated that plants accumulate organic substances emanating from various sources, including soil and water. This fact suggests the potentiality of contamination of certain vital bioresources, such as medicinal plants, by persistent contaminants, such as perfluorooctanoic acid (PFOA), perfluorooctane sulfonate (PFOS), and perfluorobutane sulfonate (PFBS). Hence, in this study, the propensity of Tagetes erecta L. (a commonly used medicinal plant) to accumulate PFOA, PFOS, and PFBS was determined using liquid chromatography/tandem mass spectrometry (LC⁻MS/MS-8030). From the results, PFOA, PFOS, and PFBS were detected in all the plant samples and concentration levels were found to be 94.83 ng/g, 5.03 ng/g, and 1.44 ng/g, respectively, with bioconcentration factor (BCF) ranges of 1.30 to 2.57, 13.67 to 72.33, and 0.16 to 0.31, respectively. Little evidence exists on the bioaccumulative susceptibility of medicinal plants to these persistent organic pollutants (POPs). These results suggest that these medicinal plants (in particular, Tagetes erecta L., used for the management of diabetes) are also potential conduits of PFOA, PFOS, and PFBS into humans.

The transcriptome analyses of Tagetes erecta provides novel insights into secondary metabolite biosynthesis during flower development

Genomic and transcriptomic sequences of Tagetes erecta are very limited in public databases, despite its nutritional and economical value. In this study, nine cDNA libraries were constructed from leaves, immature and mature flowers and sequenced using Illumina Hiseq 2000. Over 45.66 Gb nucleotides were generated and a total of 72,499 unigenes were assembled, 34,892 (48.13%) of which were annotated in NR, Swiss-prot, COG, GO, KOG, KEGG and Pfam. 11,721 differentially expressed genes were identified in leaves, immature flowers and mature flowers. The differentially expressed genes between immature flowers and mature flowers were mainly involved in photosynthesis and secondary metabolite biosynthesis. Additionally, the catalysis-related unigenes and their expression involved in carotenoids biosynthesis were determined. Using our newly identified reference genes as internal control, the expression profiles of carotenoids biosynthetic genes were verified by real-time qRT-PCR, and four of the unigenes were cloned for full length sequences confirmation. Furthermore, gene expansions occurred among the catalytic gene families in carotenoids biosynthesis pathway, which might explain the high pigment content in T. erecta.

Cloning and Functional Analysis of MADS-box Genes, TaAG-A and TaAG-B, from a Wheat K-type Cytoplasmic Male Sterile Line

Wheat (Triticum aestivum L.) is a major crop worldwide. The utilization of heterosis is a promising approach to improve the yield and quality of wheat. Although there have been many studies on wheat cytoplasmic male sterility, its mechanism remains unclear. In this study, we identified two MADS-box genes from a wheat K-type cytoplasmic male sterile (CMS) line using homology-based cloning. These genes were localized on wheat chromosomes 3A and 3B and named TaAG-A and TaAG-B, respectively. Analysis of TaAG-A and TaAG-B expression patterns in leaves, spikes, roots, and stems of Chinese Spring wheat determined using quantitative RT-PCR revealed different expression levels in different tissues. TaAG-A had relatively high expression levels in leaves and spikes, but low levels in roots, while TaAG-B had relatively high expression levels in spikes and lower expression in roots, stems, and leaves. Both genes showed downregulation during the mononucleate to trinucleate stages of pollen development in the maintainer line. In contrast, upregulation of TaAG-B was observed in the CMS line. The transcript levels of the two genes were clearly higher in the CMS line compared to the maintainer line at the trinucleate stage. Overexpression of TaAG-A and TaAG-B in Arabidopsis resulted in phenotypes with earlier reproductive development, premature mortality, and abnormal buds, stamens, and stigmas. Overexpression of TaAG-A and TaAG-B gives rise to mutants with many deformities. Silencing TaAG-A and TaAG-B in a fertile wheat line using the virus-induced gene silencing (VIGS) method resulted in plants with green and yellow striped leaves, emaciated spikes, and decreased selfing seed set rates. These results demonstrate that TaAG-A and TaAG-B may play a role in male sterility in the wheat CMS line.

Characterization and Functional Analysis of Five MADS-Box B Class Genes Related to Floral Organ Identification in Tagetes erecta

According to the floral organ development ABC model, B class genes specify petal and stamen identification. In order to study the function of B class genes in flower development of Tagetes erecta, five MADS-box B class genes were identified and their expression and putative functions were studied. Sequence comparisons and phylogenetic analyses indicated that there were one PI-like gene-TePI, two euAP3-like genes-TeAP3-1 and TeAP3-2, and two TM6-like genes-TeTM6-1 and TeTM6-2 in T. erecta. Strong expression levels of these genes were detected in stamens of the disk florets, but little or no expression was detected in bracts, receptacles or vegetative organs. Yeast hybrid experiments of the B class proteins showed that TePI protein could form a homodimer and heterodimers with all the other four B class proteins TeAP3-1, TeAP3-2, TeTM6-1 and TeTM6-2. No homodimer or interaction was observed between the euAP3 and TM6 clade members. Over-expression of five B class genes of T. erecta in Nicotiana rotundifolia showed that only the transgenic plants of 35S::TePI showed altered floral morphology compared with the non-transgenic line. This study could contribute to the understanding of the function of B class genes in flower development of T. erecta, and provide a theoretical basis for further research to change floral organ structures and create new materials for plant breeding.

Digital Gene Expression Analysis Provides Insight into the Transcript Profile of the Genes Involved in Aporphine Alkaloid Biosynthesis in Lotus (Nelumbo nucifera)

The predominant alkaloids in lotus leaves are aporphine alkaloids. These are the most important active components and have many pharmacological properties, but little is known about their biosynthesis. We used digital gene expression (DGE) technology to identify differentially-expressed genes (DEGs) between two lotus cultivars with different alkaloid contents at four leaf development stages. We also predicted potential genes involved in aporphine alkaloid biosynthesis by weighted gene co-expression network analysis (WGCNA). Approximately 335 billion nucleotides were generated; and 94% of which were aligned against the reference genome. Of 22 thousand expressed genes, 19,000 were differentially expressed between the two cultivars at the four stages. Gene Ontology (GO) enrichment analysis revealed that catalytic activity and oxidoreductase activity were enriched significantly in most pairwise comparisons. In Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, dozens of DEGs were assigned to the categories of biosynthesis of secondary metabolites, isoquinoline alkaloid biosynthesis, and flavonoid biosynthesis. The genes encoding norcoclaurine synthase (NCS), norcoclaurine 6-O-methyltransferase (6OMT), coclaurine N-methyltransferase (CNMT), N-methylcoclaurine 3'-hydroxylase (NMCH), and 3'-hydroxy-N-methylcoclaurine 4'-O-methyltransferase (4'OMT) in the common pathways of benzylisoquinoline alkaloid biosynthesis and the ones encoding corytuberine synthase (CTS) in aporphine alkaloid biosynthetic pathway, which have been characterized in other plants, were identified in lotus. These genes had positive effects on alkaloid content, albeit with phenotypic lag. The WGCNA of DEGs revealed that one network module was associated with the dynamic change of alkaloid content. Eleven genes encoding proteins with methyltransferase, oxidoreductase and CYP450 activities were identified. These were surmised to be genes involved in aporphine alkaloid biosynthesis. This transcriptomic database provides new directions for future studies on clarifying the aporphine alkaloid pathway.

Doubled Haploid 'CUDH2107' as a Reference for Bulb Onion (Allium cepa L.) Research: Development of a Transcriptome Catalogue and Identification of Transcripts Associated with Male Fertility

Researchers working on model plants have derived great benefit from developing genomic and genetic resources using ‘reference’ genotypes. Onion has a large and highly heterozygous genome making the sharing of germplasm and analysis of sequencing data complicated. To simplify the discovery and analysis of genes underlying important onion traits, we are promoting the use of the homozygous double haploid line ‘CUDH2107’ by the onion research community. In the present investigation, we performed transcriptome sequencing on vegetative and reproductive tissues of CUDH2107 to develop a multi-organ reference transcriptome catalogue. A total of 396 million 100 base pair paired reads was assembled using the Trinity pipeline, resulting in 271,665 transcript contigs. This dataset was analysed for gene ontology and transcripts were classified on the basis of putative biological processes, molecular function and cellular localization. Significant differences were observed in transcript expression profiles between different tissues. To demonstrate the utility of our CUDH2107 transcriptome catalogue for understanding the genetic and molecular basis of various traits, we identified orthologues of rice genes involved in male fertility and flower development. These genes provide an excellent starting point for studying the molecular regulation, and the engineering of reproductive traits.