RiceAntherNet: a gene co-expression network for identifying anther and pollen development genes

Recent advances in exploring transcriptional regulatory landscape of crops

Crop breeding entails developing and selecting plant varieties with improved agronomic traits. Modern molecular techniques, such as genome editing, enable more efficient manipulation of plant phenotype by altering the expression of particular regulatory or functional genes. Hence, it is essential to thoroughly comprehend the transcriptional regulatory mechanisms that underpin these traits. In the multi-omics era, a large amount of omics data has been generated for diverse crop species, including genomics, epigenomics, transcriptomics, proteomics, and single-cell omics. The abundant data resources and the emergence of advanced computational tools offer unprecedented opportunities for obtaining a holistic view and profound understanding of the regulatory processes linked to desirable traits. This review focuses on integrated network approaches that utilize multi-omics data to investigate gene expression regulation. Various types of regulatory networks and their inference methods are discussed, focusing on recent advancements in crop plants. The integration of multi-omics data has been proven to be crucial for the construction of high-confidence regulatory networks. With the refinement of these methodologies, they will significantly enhance crop breeding efforts and contribute to global food security.

StCoExpNet: a global co-expression network analysis facilitates identifying genes underlying agronomic traits in potatoes

KEY MESSAGE

We constructed a gene expression atlas and co-expression network for potatoes and identified several novel genes associated with various agronomic traits. This resource will accelerate potato genetics and genomics research. Potato (Solanum tuberosum L.) is the world's most crucial non-cereal food crop and ranks third in food production after wheat and rice. Despite the availability of several potato transcriptome datasets at public databases like NCBI SRA, an effort has yet to be put into developing a global transcriptome atlas and a co-expression network for potatoes. The objectives of our study were to construct a global expression atlas for potatoes using publicly available transcriptome datasets, identify housekeeping and tissue-specific genes, construct a global co-expression network and identify co-expression clusters, investigate the transcriptional complexity of genes involved in various essential biological processes related to agronomic traits, and provide a web server (StCoExpNet) to easily access the newly constructed expression atlas and co-expression network to investigate the expression and co-expression of genes of interest. In this study, we used data from 2299 publicly available potato transcriptome samples obtained from 15 different tissues to construct a global transcriptome atlas. We found that roughly 87% of the annotated genes exhibited detectable expression in at least one sample. Among these, we identified 281 genes with consistent and stable expression levels, indicating their role as housekeeping genes. Conversely, 308 genes exhibited marked tissue-specific expression patterns. We exemplarily linked some co-expression clusters to important agronomic traits of potatoes, such as self-incompatibility, anthocyanin biosynthesis, tuberization, and defense responses against multiple pathogens. The dataset compiled here constitutes a new resource (StCoExpNet), which can be accessed at https://stcoexpnet.julius-kuehn.de . This transcriptome atlas and the co-expression network will accelerate potato genetics and genomics research.

Dynamic gene regulatory networks improving spike fertility through regulation of floret primordia fate in wheat

The developmental process of spike is critical for spike fertility through affecting floret primordia fate in wheat; however, the genetic regulation of this dynamic and complex developmental process remains unclear. Here, we conducted a high temporal-resolution analysis of spike transcriptomes and monitored the number and morphology of floret primordia within spike. The development of all floret primordia in a spike was clearly separated into three distinct phases: differentiation, pre-dimorphism and dimorphism. Notably, we identified that floret primordia with meiosis ability at the pre-dimorphism phase usually develop into fertile floret primordia in the next dimorphism phase. Compared to control, increasing plant space treatment achieved the maximum increasement range (i.e., 50%) in number of fertile florets by accelerating spike development. The process of spike fertility improvement was directed by a continuous and dynamic regulatory network involved in transcription factor and genes interaction. This was based on the coordination of genes related to heat shock protein and jasmonic acid biosynthesis during differentiation phase, and genes related to lignin, anthocyanin and chlorophyll biosynthesis during dimorphism phase. The multi-dimensional association with high temporal-resolution approach reported here allows rapid identification of genetic resource for future breeding studies to realise the maximum spike fertility potential in more cereal crops.

Rice co-expression network analysis identifies gene modules associated with agronomic traits

Identifying trait-associated genes is critical for rice (Oryza sativa) improvement, which usually relies on map-based cloning, quantitative trait locus analysis, or genome-wide association studies. Here we show that trait-associated genes tend to form modules within rice gene co-expression networks, a feature that can be exploited to discover additional trait-associated genes using reverse genetics. We constructed a rice gene co-expression network based on the graphical Gaussian model using 8,456 RNA-seq transcriptomes, which assembled into 1,286 gene co-expression modules functioning in diverse pathways. A number of the modules were enriched with genes associated with agronomic traits, such as grain size, grain number, tiller number, grain quality, leaf angle, stem strength, and anthocyanin content, and these modules are considered to be trait-associated gene modules. These trait-associated gene modules can be used to dissect the genetic basis of rice agronomic traits and to facilitate the identification of trait genes. As an example, we identified a candidate gene, OCTOPUS-LIKE 1 (OsOPL1), a homolog of the Arabidopsis (Arabidopsis thaliana) OCTOPUS gene, from a grain size module and verified it as a regulator of grain size via functional studies. Thus, our network represents a valuable resource for studying trait-associated genes in rice.

Natural variation MeMYB108 associated with tolerance to stress-induced leaf abscission linked to enhanced protection against reactive oxygen species in cassava

Natural variation of the MeMYB108 exon was associated with reactive oxygen scavengers led to alleviate leaf abscission under drought in cassava. The reactive oxygen scavengers play important roles in regulating the cassava (Manihot esculenta Crantz) leaf abscission induced by stresses. To date, the relationship between natural variations of MYB genes and reactive oxygen scavengers under drought in cassava genotypes remains unclear. Here, we reported the transcription factor MeMYB108 played an important role in regulating leaf abscission exposed to drought in cassava. The expression levels of MeMYB108 in abscission zones of cassava leaf pulvinus were higher in cassava genotype SC124, which were less easy to shed leaves under stress than cassava genotype SC8 when the leaf abscission induced by the same drought condition. Compared with wild type and interference expression plants, overexpression of MeMYB108 significantly reduced the drought-induced leaf abscission rate under drought. The consecutively 2-year analysis of reactive oxygen scavengers showed significant differences among different cassava genotypes under drought-induced leaf abscission, indicating the relevance between reactive oxygen scavengers and leaf abscission. Correlation analysis revealed the natural variation of the MeMYB108 exon was associated with reactive oxygen scavengers during drought-induced leaf abscission. Association analysis between pairwise LD of DNA polymorphism indicated the MeMYB108 allele enhanced the tolerance of cassava to drought-induced leaf abscission. Complementation transgenic lines containing the elite allele of MeMYB108 ^SC124 decreased the leaf abscission rate induced by drought conditions, demonstrating natural variation in MeMYB108 contributed to leaf abscission tolerance induced by drought in cassava. Further studies showed MeMYB108 played an active role in the tolerance of cassava to drought-induced leaf abscission by inducing scavenging of reactive oxygen species.

Comparative transcriptome analysis of pollen and anther wall reveals novel insights into the regulatory mechanisms underlying anther wall development and its dehiscence in rice

To further understand the regulatory mechanism for anther dehiscence in rice, we carried out transcriptome analysis for the following two tissues: the anther wall and pollen at the anthesis stage. With the anatomical meta-expression data, in addition to these tissues, the differentially expressed genes (DEGs) between the two tissues were further refined to identify 1,717 pollen-preferred genes and 534 anther wall-preferred genes. A GUS transgenic line and RT-qPCR analysis for anther wall-preferred genes supported the fidelity of our gene candidates for further analysis. The refined DEGs were functionally classified through Gene Ontology (GO) enrichment and MapMan analyses. Through the analysis of cis-acting elements and alternative splicing variants, we also suggest the feature of regulatory sequences in promoter regions for anther wall-preferred expression and provide information of the unique splicing variants in anther wall. Subsequently, it was found that hormone signaling and the resulting transcriptional regulation pathways may play an important role in anther dehiscence and anther wall development. Our results could provide useful insights into future research to broaden the molecular mechanism of anther dehiscence or anther wall development in rice.

The rice phosphate transporter OsPHT1;7 plays a dual role in phosphorus redistribution and anther development

Inorganic phosphate (Pi) is the predominant form of phosphorus (P) readily accessible to plants, and Pi Transporter 1 (PHT1) genes are the major contributors to root Pi uptake. However, the mechanisms underlying the transport and recycling of Pi within plants, which are vital for optimizing P use efficiency, remain elusive. Here, we characterized a functionally unknown rice (Oryza sativa) PHT1 member barely expressed in roots, OsPHT1;7. Yeast complementation and Xenopus laevis oocyte assay demonstrated that OsPHT1;7 could mediate Pi transport. Reverse-transcription quantitative polymerase chain reaction and histochemical analyses showed that OsPHT1;7 was preferentially expressed in source leaves and nodes. A further fine-localization analysis by immunostaining showed that OsPHT1;7 expression was restricted in the vascular bundle (VB) sheath and phloem of source leaves as well as in the phloem of regular/diffuse- and enlarged-VBs of nodes. In accordance with this expression pattern, mutation of OsPHT1;7 led to increased and decreased P distribution in source (old leaves) and sink organs (new leaves/panicles), respectively, indicating that OsPHT1;7 is involved in P redistribution. Furthermore, OsPHT1;7 showed an overwhelmingly higher transcript abundance in anthers than other PHT1 members, and ospht1;7 mutants were impaired in P accumulation in anthers but not in pistils or husks. Moreover, the germination of pollen grains was significantly inhibited upon OsPHT1;7 mutation, leading to a >80% decrease in seed-setting rate and grain yield. Taken together, our results provide evidence that OsPHT1;7 is a crucial Pi transporter for Pi transport and recycling within rice plants, stimulating both vegetative and reproductive growth.

The proteome of developing barley anthers during meiotic prophase I

Flowering plants reproduce sexually by combining a haploid male and female gametophyte during fertilization. Male gametophytes are localized in the anthers, each containing reproductive (meiocyte) and non-reproductive tissue necessary for anther development and maturation. Meiosis, where chromosomes pair and exchange their genetic material during a process called recombination, is one of the most important and sensitive stages in breeding, ensuring genetic diversity. Most anther development studies have focused on transcript variation, but very few have been correlated with protein abundance. Taking advantage of a recently published barley anther transcriptomic (BAnTr) dataset and a newly developed sensitive mass spectrometry-based approach to analyse the barley anther proteome, we conducted high-resolution mass spectrometry analysis of barley anthers, collected at six time points and representing their development from pre-meiosis to metaphase. Each time point was carefully staged using immunocytology, providing a robust and accurate staging mirroring our previous BAnTr dataset. We identified >6100 non-redundant proteins including 82 known and putative meiotic proteins. Although the protein abundance was relatively stable throughout prophase I, we were able to quantify the dynamic variation of 336 proteins. We present the first quantitative comparative proteomics study of barley anther development during meiotic prophase I when the important process of homologous recombination is taking place.

A putative SUBTILISIN-LIKE SERINE PROTEASE 1 (SUBSrP1) regulates anther cuticle biosynthesis and panicle development in rice

INTRODUCTION

Panicle abortion is a severe physiological defect and causes a reduction in grain yield.

OBJECTIVES

In this study, we aim to provide the characterization and functional analysis of a mutant apa1331 (apical panicle abortion1331).

METHODS

The isolated mutant from an EMS-mutagenized population was subjected to SSR analysis and Mutmap assay for candidate gene mapping. We performed phenotypic analysis, anthers cross-sections morphology, wax and cutin profiling, biochemical assays and phylogenetic analysis for characterization and evaluation of apa1331. We used CRISPR/Cas9 disruption for functional validation of its candidate gene. Furthermore, comparative RNA-seq and relative expression analysis were performed to get further insights into mechanistic role of the candidate gene.

RESULTS

The anthers from the apical spikelets of apa1331 were degenerated, pollen-less and showed defects in cuticle formation. Transverse sections of apa1331 anthers showed defects in post-meiotic microspore development at stage 8-9. Gas Chromatography showed a significant reduction of wax and cutin in anthers of apa1331 compared to Wildtype (WT). Quantification of H₂O₂ and MDA has indicated the excessive ROS (reactive oxygen species) in apa1331. Trypan blue staining and TUNEL assay revealed cell death and excessive DNA fragmentation in apa1331. Map-based cloning and Mutmap analysis revealed that LOC_Os04g40720, encoding a putative SUBTILISIN-LIKE SERINE PROTEASE (OsSUBSrP1), harbored an SNP (A > G) in apa1331. Phenotypic defects were only seen in apical spikelets due to highest expression of OsSUBSrP1 in upper panicle portion. CRISPR-mediated knock-out lines of OsSUBSrP1 displayed spikelet abortion comparable to apa1331. Global gene expression analysis revealed a significant downregulation of wax and cutin biosynthesis genes.

CONCLUSIONS

Our study reports the novel role of SUBSrP1 in anther cuticle biosynthesis by ROS-mediated programmed cell death in rice.

Gene Co-expression Network Analysis and Linking Modules to Phenotyping Response in Plants

Environmental factors, including different stresses, can have an impact on the expression of genes and subsequently the phenotype and development of plants. Since a large number of genes are involved in response to the perturbation of the environment, identifying groups of co-expressed genes is meaningful. The gene co-expression network models can be used for the exploration, interpretation, and identification of genes responding to environmental changes. Once a gene co-expression network is constructed, one can determine gene modules and the association of gene modules to the phenotypic response. To link modules to phenotype, one approach is to find the correlated eigengenes of given modules or to integrate all eigengenes in regularized linear model. This manuscript describes the method from construction of co-expression network, module discovery, association between modules and phenotypic data, and finally to annotation/visualization.

Conserved co-functional network between maize and Arabidopsis aid in the identification of seed defective genes in maize

BACKGROUND

The biological pathways related to Arabidopsis seed development have been well studied and functional genes involved in it have been discovered. However, functional studies about maize seed development were more limited compared to Arabidopsis.

OBJECTIVE

Therefore, transferring knowledge from Arabidopsis into maize would facilitate functional studies about maize seed development.

METHOD

In this study, public transcriptome data of the two species related to seed development were obtained. Co-expression network in each species was compared by integrating orthology information.

RESULTS

This conserved co-functional network contained 4510 maize and 4808 Arabidopsis genes, respectively. Most of these genes were expressed in throughout embryo, early or later endosperm/seed. These conserved co-functional genes were significantly enriched for members of PPR protein family, which was consistent with that PPR proteins play an important role in maize seed development. Spatial-temporally co-functional genes were discovered in the seed coat and embryo. Furthermore, 66 well-studied genes involved in Arabidopsis seed development were co-functional with 319 maize genes and one maize gene (GRMZM2G036050) was further confirmed using an EMS-induced seed defective mutant by bulked segregating RNA sequencing (BSR) analysis.

CONCLUSIONS

Altogether, these results showed the potential of this approach to support functional studies in maize seed development by transferring knowledge from Arabidopsis.

Transcriptome-based approaches for clarification of nutritional responses and improvement of crop production

Genome-wide transcriptome profiling is a powerful tool for identifying key genes and pathways involved in plant development and physiological processes. This review summarizes studies that have used transcriptome profiling mainly in rice to focus on responses to macronutrients such as nitrogen, phosphorus and potassium, and spatio-temporal root profiling in relation to the regulation of root system architecture as well as nutrient uptake and transport. We also discuss strategies based on meta- and co-expression analyses with different attributed transcriptome data, which can be used for investigating the regulatory mechanisms and dynamics of nutritional responses and adaptation, and speculate on further advances in transcriptome profiling that could have potential application to crop breeding and cultivation.

Gene Expression Dynamics in Rice Peduncles at the Heading Stage

Improving grain yield in the staple food crop rice has been long sought goal of plant biotechnology. One of the traits with significant impact on rice breeding programs is peduncle elongation at the time of heading failing which leads to significant reduction in grain yield due to incomplete panicle exsertion. To decipher transcriptional dynamics and molecular players underlying peduncle elongation, we performed RNA sequencing analysis of elongating and non-elongating peduncles in two Indian cultivars, Swarna and Pokkali, at the time of heading. Along with genes associated with cell division and cell wall biosynthesis, we observed significant enrichment of genes associated with auxins, gibberellins, and brassinosteroid biosynthesis/signaling in the elongating peduncles before heading in both the genotypes. Similarly, genes associated with carbohydrate metabolism and mobilization, abiotic stress response along with cytokinin, abscisic acid, jasmonic acid, and ethylene biosynthesis/signaling were enriched in non-elongating peduncles post heading. Significant enrichment of genes belonging to key transcription factor families highlights their specialized roles in peduncle elongation and grain filling before and after heading, respectively. A comparison with anther/pollen development-related genes provided 76 candidates with overlapping roles in anther/pollen development and peduncle elongation. Some of these are important for carbohydrate remobilization to the developing grains. These can be engineered to combat with incomplete panicle exsertion in male sterile lines and manipulate carbohydrate dynamics in grasses. Overall, this study provides baseline information about potential target genes for engineering peduncle elongation with implications on plant height, biomass composition and grain yields in rice.

Cytological and molecular characterizations of a novel 2A nullisomic line derived from a widely-grown wheat cultivar Zhoumai 18 conferring male sterility

A dwarf, multi-pistil and male sterile dms mutant was previously reported by us. However, the genetic changes in this dms are unclear. To examine the genetic changes, single nucleotide polymorphism (SNP) association, chromosome counting, and high-resolution chromosome fluorescence in situ hybridization (FISH) techniques were employed. By comparing tall plants (T) with dwarf plants (D) in the offspring of dms mutant plants, SNP association analysis indicated that most SNPs were on chromosome 2A. There were three types in offspring of dms plants, with 42, 41 and 40 chromosomes respectively. High-resolution chromosome painting analysis demonstrated that T plants had all 42 wheat chromosomes; the medium plants (M) had 41 chromosomes, lacking one chromosome 2A; while D plants had 40 wheat chromosomes, and lacked both 2A chromosomes. These data demonstrated that dms resulted from a loss of chromosome 2A. We identified 23 genes on chromosome 2A which might be involved in the development of stamens or pollen grains. These results lay a solid foundation for further analysis of the molecular mechanisms of wheat male sterility. Because D plants can be used as a female parent to cross with other wheat genotypes, dms is a unique germplasm for any functional study of chromosome 2A and wheat breeding specifically targeting genes on 2A.

Gene Regulatory Network Inference: Connecting Plant Biology and Mathematical Modeling

Plant responses to environmental and intrinsic signals are tightly controlled by multiple transcription factors (TFs). These TFs and their regulatory connections form gene regulatory networks (GRNs), which provide a blueprint of the transcriptional regulations underlying plant development and environmental responses. This review provides examples of experimental methodologies commonly used to identify regulatory interactions and generate GRNs. Additionally, this review describes network inference techniques that leverage gene expression data to predict regulatory interactions. These computational and experimental methodologies yield complex networks that can identify new regulatory interactions, driving novel hypotheses. Biological properties that contribute to the complexity of GRNs are also described in this review. These include network topology, network size, transient binding of TFs to DNA, and competition between multiple upstream regulators. Finally, this review highlights the potential of machine learning approaches to leverage gene expression data to predict phenotypic outputs.

Transcriptional trajectories of anther development provide candidates for engineering male fertility in sorghum

Sorghum is a self-pollinated crop with multiple economic uses as cereal, forage, and biofuel feedstock. Hybrid breeding is a cornerstone for sorghum improvement strategies that currently relies on cytoplasmic male sterile lines. To engineer genic male sterility, it is imperative to examine the genetic components regulating anther/pollen development in sorghum. To this end, we have performed transcriptomic analysis from three temporal stages of developing anthers that correspond to meiotic, microspore and mature pollen stages. A total of 5286 genes were differentially regulated among the three anther stages with 890 of them exhibiting anther-preferential expression. Differentially expressed genes could be clubbed into seven distinct developmental trajectories using K-means clustering. Pathway mapping revealed that genes involved in cell cycle, DNA repair, regulation of transcription, brassinosteroid and auxin biosynthesis/signalling exhibit peak expression in meiotic anthers, while those regulating abiotic stress, carbohydrate metabolism, and transport were enriched in microspore stage. Conversely, genes associated with protein degradation, post-translational modifications, cell wall biosynthesis/modifications, abscisic acid, ethylene, cytokinin and jasmonic acid biosynthesis/signalling were highly expressed in mature pollen stage. High concurrence in transcriptional dynamics and cis-regulatory elements of differentially expressed genes in rice and sorghum confirmed conserved developmental pathways regulating anther development across species. Comprehensive literature survey in conjunction with orthology analysis and anther-preferential accumulation enabled shortlisting of 21 prospective candidates for in-depth characterization and engineering male fertility in sorghum.

Cytological and transcriptome analyses reveal OsPUB73 defect affects the gene expression associated with tapetum or pollen exine abnormality in rice

BACKGROUND

As one of the main crops in the world, sterility of rice (Oryza sativa L.) significantly affects the production and leads to yield decrease. Our previous research showed that OsPUB73, which encodes U-box domain-containing protein 73, may be associated with male sterility. However, little information is available on this gene that is required for anther development. In the present study, we knocked out OsPUB73 by using the CRISPR/Cas9 system and studied the cytological and transcriptome of the gene-defect associated with pollen development and sterility in the rice variety (Taichung 65).

RESULTS

The sequence analysis indicated that OsPUB73 was comprised of 3 exons and 2 introns, of which CDS encoded 586 amino acids including a U-box domain. The expression pattern of OsPUB73 showed that it was highly expressed in the anther during meiosis stage. The ospub73 displayed low pollen fertility (19.45%), which was significantly lower than wild type (WT) (85.37%). Cytological observation showed tapetum vacuolated at the meiosis stage and pollen exine was abnormal at the bi-cellular pollen stage of ospub73. RNA-seq analysis detected 2240 down and 571 up-regulated genes in anther of ospub73 compared with WT during meiosis stage. Among of 2240 down-regulated genes, seven known genes were associated with tapetal cell death or pollen exine development, including CYP703A3 (Cytochrome P450 Hydroxylase703A3), CYP704B2 (Cytochrome P450 Hydroxylase704B2), DPW (Defective Pollen Wall), PTC1 (Persistant Tapetal Cell1), UDT1 (Undeveloped Tapetum1), OsAP37 (Aspartic protease37) and OsABCG15 (ATP binding cassette G15), which were validated by quantitative real-time polymerase chain reaction (qRT-PCR). These results suggested OsPUB73 may play an important role in tapetal or pollen exine development and resulted in pollen partial sterility.

CONCLUSION

Our results revealed that OsPUB73 plays an important role in rice male reproductive development, which provides valuable information about the molecular mechanisms of the U-box in rice male reproductive development.

Cytological and transcriptome analysis reveal that interaction at Sb pollen sterility locus cause down-regulation of important meiosis-related genes associated with high pollen sterility in autotetraploid rice hybrids

Polyploidy could increase the interactions of pollen sterility loci and Sb locus interaction cause higher pollen abortion than other loci. Therefore, we focused on the interaction at Sb pollen sterility locus in autotetraploid rice compared to diploid rice hybrid using the near-isogenic lines in the present study. Cytological observations indicated that interaction at Sb locus cause high pollen sterility (69.9%) and abnormal chromosome behavior (37.02%) at Metaphase II in autotetraploid rice hybrid. A total of 139 meiosis-related or meiosis stage-specific genes were detected in the autotetraploid rice hybrid harboring interaction at Sb locus and 27 of these meiosis-related or specific genes displayed significant down-regulation, including four pollen fertility related genes (Rad51, XRI1, PSS1 and MIL1). These results revealed a stronger interaction at Sb pollen sterility locus than other loci, which cause down-regulation of many important meiosis-related genes that were associated with higher pollen sterility in autotetraploid rice hybrids.

Proteomic analysis reveals that sugar and fatty acid metabolisms play a central role in sterility of the male-sterile line 1355A of cotton

Cotton (Gossypium spp.) is one of the most important economic crops and exhibits yield-improving heterosis in specific hybrid combinations. The genic male-sterility system is the main strategy used for producing heterosis in cotton. To better understand the mechanisms of male sterility in cotton, we carried out two-dimensional electrophoresis (2-DE) and label-free quantitative proteomics analysis in the anthers of two near-isogenic lines, the male-sterile line 1355A and the male-fertile line 1355B. We identified 39 and 124 proteins that were significantly differentially expressed between these two lines in the anthers at the tetrad stage (stage 7) and uninucleate pollen stage (stage 8), respectively. Gene ontology-based analysis revealed that these differentially expressed proteins were mainly associated with pyruvate, carbohydrate, and fatty acid metabolism. Biochemical analysis revealed that in the anthers of line 1355A, glycolysis was activated, which was caused by a reduction in fructose, glucose, and other soluble sugars, and that accumulation of acetyl-CoA was increased along with a significant increase in C14:0 and C18:1 free fatty acids. However, the activities of pyruvate dehydrogenase and fatty acid biosynthesis were inhibited and fatty acid β-oxidation was activated at the translational level in 1355A. We speculate that in the 1355A anther, high rates of glucose metabolism may promote fatty acid synthesis to enable anther growth. These results provide new insights into the molecular mechanism of genic male sterility in upland cotton.

Recent advances in gene function prediction using context-specific coexpression networks in plants

Predicting gene functions from genome sequence alone has been difficult, and the functions of a large fraction of plant genes remain unknown. However, leveraging the vast amount of currently available gene expression data has the potential to facilitate our understanding of plant gene functions, especially in determining complex traits. Gene coexpression networks-created by integrating multiple expression datasets-connect genes with similar patterns of expression across multiple conditions. Dense gene communities in such networks, commonly referred to as modules, often indicate that the member genes are functionally related. As such, these modules serve as tools for generating new testable hypotheses, including the prediction of gene function and importance. Recently, we have seen a paradigm shift from the traditional "global" to more defined, context-specific coexpression networks. Such coexpression networks imply genetic correlations in specific biological contexts such as during development or in response to a stress. In this short review, we highlight a few recent studies that attempt to fill the large gaps in our knowledge about cellular functions of plant genes using context-specific coexpression networks.

Wheat TaMs1 is a glycosylphosphatidylinositol-anchored lipid transfer protein necessary for pollen development

BACKGROUND

In flowering plants, lipid biosynthesis and transport within anthers is essential for male reproductive success. TaMs1, a dominant wheat fertility gene located on chromosome 4BS, has been previously fine mapped and identified to encode a glycosylphosphatidylinositol (GPI)-anchored non-specific lipid transfer protein (nsLTP). Although this gene is critical for pollen exine development, details of its function remains poorly understood.

RESULTS

In this study, we report that TaMs1 is only expressed from the B sub-genome, with highest transcript abundance detected in anthers containing microspores undergoing pre-meiosis through to meiosis. β-glucuronidase transcriptional fusions further revealed that TaMs1 is expressed throughout all anther cell-types. TaMs1 was identified to be expressed at an earlier stage of anther development relative to genes reported to be necessary for sporopollenin precursor biosynthesis. In anthers missing a functional TaMs1 (ms1c deletion mutant), these same genes were not observed to be mis-regulated, indicating an independent function for TaMs1 in pollen development. Exogenous hormone treatments on GUS reporter lines suggest that TaMs1 expression is increased by both indole-3-acetic acid (IAA) and abscisic acid (ABA). Translational fusion constructs showed that TaMs1 is targeted to the plasma membrane.

CONCLUSIONS

In summary, TaMs1 is a wheat fertility gene, expressed early in anther development and encodes a GPI-LTP targeted to the plasma membrane. The work presented provides a new insight into the process of wheat pollen development.

An integrated omics analysis reveals molecular mechanisms that are associated with differences in seed oil content between Glycine max and Brassica napus

BACKGROUND

Rapeseed (Brassica napus L.) and soybean (Glycine max L.) seeds are rich in both protein and oil, which are major sources of biofuels and nutrition. Although the difference in seed oil content between soybean (~ 20%) and rapeseed (~ 40%) exists, little is known about its underlying molecular mechanism.

RESULTS

An integrated omics analysis was performed in soybean, rapeseed, Arabidopsis (Arabidopsis thaliana L. Heynh), and sesame (Sesamum indicum L.), based on Arabidopsis acyl-lipid metabolism- and carbon metabolism-related genes. As a result, candidate genes and their transcription factors and microRNAs, along with phylogenetic analysis and co-expression network analysis of the PEPC gene family, were found to be largely associated with the difference between the two species. First, three soybean genes (Glyma.13G148600, Glyma.13G207900 and Glyma.12G122900) co-expressed with GmPEPC1 are specifically enriched during seed storage protein accumulation stages, while the expression of BnPEPC1 is putatively inhibited by bna-miR169, and two genes BnSTKA and BnCKII are co-expressed with BnPEPC1 and are specifically associated with plant circadian rhythm, which are related to seed oil biosynthesis. Then, in de novo fatty acid synthesis there are rapeseed-specific genes encoding subunits β-CT (BnaC05g37990D) and BCCP1 (BnaA03g06000D) of heterogeneous ACCase, which could interfere with synthesis rate, and β-CT is positively regulated by four transcription factors (BnaA01g37250D, BnaA02g26190D, BnaC01g01040D and BnaC07g21470D). In triglyceride synthesis, GmLPAAT2 is putatively inhibited by three miRNAs (gma-miR171, gma-miR1516 and gma-miR5775). Finally, in rapeseed there was evidence for the expansion of gene families, CALO, OBO and STERO, related to lipid storage, and the contraction of gene families, LOX, LAH and HSI2, related to oil degradation.

CONCLUSIONS

The molecular mechanisms associated with differences in seed oil content provide the basis for future breeding efforts to improve seed oil content.

MCRiceRepGP: a framework for the identification of genes associated with sexual reproduction in rice

Rice is an important cereal crop, being a staple food for over half of the world's population, and sexual reproduction resulting in grain formation underpins global food security. However, despite considerable research efforts, many of the genes, especially long intergenic non-coding RNA (lincRNA) genes, involved in sexual reproduction in rice remain uncharacterized. With an increasing number of public resources becoming available, information from different sources can be combined to perform gene functional annotation. We report the development of MCRiceRepGP, a machine learning framework which integrates heterogeneous evidence and employs multicriteria decision analysis and machine learning to predict coding and lincRNA genes involved in sexual reproduction in rice. The rice genome was reannotated using deep-sequencing transcriptomic data from reproduction-associated tissue/cell types identifying previously unannotated putative protein-coding genes and lincRNAs. MCRiceRepGP was used for genome-wide discovery of sexual reproduction associated coding and lincRNA genes. The protein-coding and lincRNA genes identified have distinct expression profiles, with a large proportion of lincRNAs reaching maximum expression levels in the sperm cells. Some of the genes are potentially linked to male- and female-specific fertility and heat stress tolerance during the reproductive stage. MCRiceRepGP can be used in combination with other genome-wide studies, such as genome-wide association studies, giving greater confidence that the genes identified are associated with the biological process of interest. As more data, especially about mutant plant phenotypes, become available, the power of MCRiceRepGP will grow, providing researchers with a tool to identify candidate genes for future experiments. MCRiceRepGP is available as a web application (http://mcgplannotator.com/MCRiceRepGP/).