RGeasy: a reference gene analysis tool for gene expression studies via RT-qPCR.

Gene expression through RT-qPCR can be performed by the relative quantification method, which requires the expression normalization through reference genes. Therefore, it is essential to validate, experimentally, the candidate reference genes. Thus, although there are several studies that are performed to identify the most stable reference genes, most them validate genes for very specific conditions, not exploring the whole potential of the research since not all possible combinations of treatments and/or conditions of the study are explored. For this reason, new experiments must be conducted by researchers that have interest in analyzing gene expression of treatments and/or conditions present, but not explored, in these studies. Here, we present the RGeasy tool, which aims to facilitate the selection of reference genes, allowing the user to choose genes for a greater number of combinations of treatments/conditions, compared to the ones present in the original articles, through just a few clicks. RGeasy was validated with RT-qPCR data from gene expression studies performed in two coffee species, Coffea arabica and Coffea canephora, and it can be used for any animal, plant or microorganism species. In addition to displaying a rank of the most stable reference genes for each condition or treatment, the user also has access to the primer pairs for the selected reference genes.

Studying the influence of TaGW8 and TaGS5-3A genes on the yield of soft spring wheat in arid climate conditions of the Republic of Kazakhstan.

Drought is a primary ecological stress limiting wheat yield in water-deficient regions. Conducting targeted genetic selection of wheat cultivars can expedite the adaptation process of wheat to the climatic conditions of the region, allowing for the identification of high-yielding varieties with stable genetic traits. This study investigated the impact of the TaGW8 and TaGS3A genes, known for their contribution to wheat productivity. The effective productivity genes TaGW8-B1b/B1a and the TaGS5-3A-T genome exert a 32.8% influence on the variability of the 1000 grain weight (TGW) trait. This influence stems from both individual genes and their interactions, with at least 17.5% of TGW variability explained by the gene combinations examined in the study. Notably, the TaGS5-3A-T gene exhibits a significant positive correlation with total yield, exceeding 63%. The integration of these productivity genes, based on field phenotypic data, has resulted in an overall yield increase of selected samples by 0.8 tons/ha compared to the country's average multi-year indicator.

The plastomes of Lepismium cruciforme (Vell.) Miq and Schlumbergera truncata (Haw.) Moran reveal tribe-specific rearrangements and the first loss of the trnT-GGU gene in Cactaceae.

BACKGROUND: Recent studies have revealed atypical features in the plastomes of the family Cactaceae, the largest lineage of succulent species adapted to arid and semi-arid regions. Most plastomes sequenced to date are from short-globose and cylindrical cacti, while little is known about plastomes of epiphytic cacti. Published cactus plastomes reveal reduction and complete loss of IRs, loss of genes, pseudogenization, and even degeneration of tRNA structures. Aiming to contribute with new insights into the plastid evolution of Cactaceae, particularly within the tribe Rhipsalideae, we de novo assembled and analyzed the plastomes of Lepismium cruciforme and Schlumbergera truncata, two South American epiphytic cacti. METHODS AND RESULTS: Our data reveal many gene losses in both plastomes and the first loss of functionality of the trnT-GGU gene in Cactaceae. The trnT-GGU is a pseudogene in L. cruciforme plastome and appears to be degenerating in the tribe Rhipsalideae. Although the plastome structure is conserved among the species of the tribe Rhipsalideae, with tribe-specific rearrangements, we mapped around 200 simple sequence repeats and identified nine nucleotide polymorphism hotspots, useful to improve the phylogenetic resolutions of the Rhipsalideae. Furthermore, our analysis indicated high gene divergence and rapid evolution of RNA editing sites in plastid protein-coding genes in Cactaceae. CONCLUSIONS: Our findings show that some characteristics of the Rhipsalideae tribe are conserved, such as plastome structure with IRs containing only the ycf2 and two tRNA genes, structural degeneration of the trnT-GGU gene and ndh complex, and lastly, pseudogenization of rpl33 and rpl23 genes, both plastid translation-related genes.

Reference genes for RT-qPCR analysis in Musa acuminata genotypes contrasting in resistance to Fusarium oxysporum f. sp. cubense subtropical race 4.

Banana (Musa spp.) is the most widely consumed fruit globally. Fusarium wilt, caused by Fusarium oxysporum f. sp. cubense (Foc), is a highly threatening disease to banana production. Resistance genes to Foc exist in wild Musa genotypes such as Musa acuminata subsp. burmannicoides var. Calcutta 4. Whilst real-time PCR (RT-qPCR) is appropriate for accurate analysis of gene expression in pathways involved in host defence responses, reference genes with stable expression under specific biotic stress conditions and host tissue types are necessary for normalization of sample variation. In this context, the stability in potential host reference genes ACT1, APT, EF1α, GAPDH, αTUB, RAN, UBIQ1, UBIQ2, ßTUB1, ßTUB3, L2 and ACTA1 was evaluated in total RNA samples from root tissues in Calcutta 4 (resistant) and Musa sp. cultivar Prata-anã (susceptible) extracted during interaction with Foc subtropical race 4 (STR4). Expression stability was calculated using the algorithms geNorm, NormFinder and BestKeeper. ßTUB3 and L2 were identified as the most stable in Calcutta 4, with ACTA1 and GAPDH the most stable in Prata-anã. These reference genes for analysis of gene expression modulation in the Musa-Foc STR4 pathosystem are fundamental for advancing understanding of host defence responses to this important pathogen.

Carica papaya L. sex chromosome review and physical mapping of the serk 2, svp-like and mdar 4 sequences.

Physical mapping evidences the chromosome organization and structure. Despite the data about plant cytogenomics, physical mapping has been conducted from single-copy and/or low-copy genes for few species. Carica papaya cytogenomics has been accomplished from BAC-FISH and repeatome sequences. We aimed to map the serk 2, svp-like and mdar 4 sequences in C. papaya. The sequences were amplified and the amplicons sequenced, showing similarity in relation to serk 2, svp-like and mdar 4 genes. Carica papaya diploidy was confirmed and the mitotic chromosomes characterized. The chromosome 1 exhibited the secondary constriction pericentromeric to the centromere of the long arm. So, we concluded that it is the sex chromosomes. serk 2 was mapped in the long arm interstitial portion of the sex chromosomes, and the interphase nuclei showed two fluorescence signals. Considering these results and the sequencing data from the C. papaya sex chromosomes, svp-like and mdar 4 genes were mapped in the interstitial region of the sex chromosome long arm. Both sequences showed only one fluorescence signal in the interphase nuclei. The procedure adopted here can be reproduced for other single-copy and/or low-copy genes, allowing the construction of cytogenetic maps. In addition, we revisited the cytogenomics data about C. papaya sex chromosomes, presenting a revised point of view about the structure and evolution to these chromosomes.

Genome-wide association analyses using multilocus models on bananas (Musa spp.) reveal candidate genes related to morphology, fruit quality, and yield.

Bananas (Musa spp.) are an essential fruit worldwide and rank as the fourth most significant food crop for addressing malnutrition due to their rich nutrients and starch content. The potential of their genetic diversity remains untapped due to limited molecular breeding tools. Our study examined a phenotypically diverse group of 124 accessions from the Colombian Musaceae Collection conserved in AGROSAVIA. We assessed 12 traits categorized into morphology, fruit quality, and yield, alongside sequence data. Our sequencing efforts provided valuable insights, with an average depth of about 7× per accession, resulting in 187,133 single-nucleotide polymorphisms (SNPs) against Musa acuminata (A genome) and 220,451 against Musa balbisiana (B genome). Population structure analysis grouped samples into four and five clusters based on the reference genome. By using different association models, we identified marker-trait associations (MTAs). The mixed linear model revealed four MTAs, while the Bayesian-information and linkage-disequilibrium iteratively nested keyway and fixed and random model for circulating probability unification models identified 82 and 70 MTAs, respectively. We identified 38 and 40 candidate genes in linkage proximity to significant MTAs for the A genome and B genome, respectively. Our findings provide insights into the genetic underpinnings of morphology, fruit quality, and yield. Once validated, the SNP markers and candidate genes can potentially drive advancements in genomic-guided breeding strategies to enhance banana crop improvement.

Designer TALEs enable discovery of cell death-inducer genes.

Transcription activator-like effectors (TALEs) in plant-pathogenic Xanthomonas bacteria activate expression of plant genes and support infection or cause a resistance response. PthA4AT is a TALE with a particularly short DNA-binding domain harboring only 7.5 repeats which triggers cell death in Nicotiana benthamiana; however, the genetic basis for this remains unknown. To identify possible target genes of PthA4AT that mediate cell death in N. benthamiana, we exploited the modularity of TALEs to stepwise enhance their specificity and reduce potential target sites. Substitutions of individual repeats suggested that PthA4AT-dependent cell death is sequence specific. Stepwise addition of repeats to the C-terminal or N-terminal end of the repeat region narrowed the sequence requirements in promoters of target genes. Transcriptome profiling and in silico target prediction allowed the isolation of two cell death inducer genes, which encode a patatin-like protein and a bifunctional monodehydroascorbate reductase/carbonic anhydrase protein. These two proteins are not linked to known TALE-dependent resistance genes. Our results show that the aberrant expression of different endogenous plant genes can cause a cell death reaction, which supports the hypothesis that TALE-dependent executor resistance genes can originate from various plant processes. Our strategy further demonstrates the use of TALEs to scan genomes for genes triggering cell death and other relevant phenotypes.

Comparative analysis of infected cassava root transcriptomics reveals candidate genes for root rot disease resistance.

Cassava root-rot incited by soil-borne pathogens is one of the major diseases that reduces root yield. Although the use of resistant cultivars is the most effective method of management, the genetic basis for root-rot resistance remains poorly understood. Therefore, our work analyzed the transcriptome of two contrasting genotypes (BRS Kiriris/resistant and BGM-1345/susceptible) using RNA-Seq to understand the molecular response and identify candidate genes for resistance. Cassava seedlings (resistant and susceptible to root-rot) were both planted in infested and sterilized soil and samples from Initial-time and Final-time periods, pooled. Two controls were used: (i) seedlings collected before planting in infested soil (absolute control) and, (ii) plants grown in sterilized soil (mock treatments). For the differentially expressed genes (DEGs) analysis 23.912 were expressed in the resistant genotype, where 10.307 were differentially expressed in the control treatment, 15 DEGs in the Initial Time-period and 366 DEGs in the Final Time-period. Eighteen candidate genes from the resistant genotype were related to plant defense, such as the MLP-like protein 31 and the peroxidase A2-like gene. This is the first model of resistance at the transcriptional level proposed for the cassava × root-rot pathosystem. Gene validation will contribute to screening for resistance of germplasm, segregating populations and/or use in gene editing in the pursuit to develop most promising cassava clones with resistance to root-rot.

Micropropagation of Chayote (Sechium edule L.) var. virens levis in RITA®.

Chayote (Sechium edule) belongs to the Cucurbitaceae family, an important family at the nutritional and medicinal levels, that has been covering international markets. Having vigorous and healthy plants is important for producers, who are very interested in cultivating chayote plants obtained from in vitro tissue culture in their orchards. Bioreactors have become an alternative with high potential for plant propagation, showing significant advantages over micropropagation in semisolid medium, by generating more plant material, larger, and more vigorous. In this chapter, a micropropagation protocol of S. edule in RITA® bioreactors is reported.

Differential expression of genes potentially related to the callogenesis and in situ hybridization of SERK gene in macaw palm (Acrocomia aculeata Jacq.) Lodd. ex Mart.

For the purpose of understanding the molecular processes triggered during callus formation in macaw palm, the expression of seven genes potentially involved in this process, identified in previous studies and from the literature, was investigated by RT-qPCR. In addition, in situ hybridization of the SERK gene was performed. Leaf tissues from adult plants from two macaw palm accession were inoculated in a medium combined with Picloram at a concentration of 450 µM to induce callus. The expression analysis was performed from leaf samples from two accessions of different origins (Municipalities of Tiros, MG, and Buriti Vermelho, DF, Brazil), which are characterized as non-responsive (NR) and responsive (R), respectively. The material was collected before callus induction (0 DAI, initial day) and 120 days after callus induction (120 DAI). Genes related to development (SERK, OASA, EF1, ANN1) and stress (LEA, CAT2, and MDAR5) were evaluated. The results obtained showed that all the genes involved with the development had their expressions downregulated at 0 DAI when the accession R was compared with the accession NR. On the other hand, it was possible to observe that these genes were upregulated at 120 DAI. The LEA stress gene showed a tendency to increase expression in the NR accession, while the R accession showed decreased expression and the CAT2 and MDAR5 genes showed upregulation in both accessions. In situ hybridization showed SERK transcripts in the vascular bundles, indicating the expression of SERK in this region, in addition to its expression in calluses. The results obtained in this study support our hypothesis that the regulation of genes involved in the control of oxidative stress and development is crucial for the formation of calluses in macaw palm.

A global-temporal analysis on Phytophthora sojae resistance-gene efficacy.

Plant disease resistance genes are widely used in agriculture to reduce disease outbreaks and epidemics and ensure global food security. In soybean, Rps (Resistance to Phytophthora sojae) genes are used to manage Phytophthora sojae, a major oomycete pathogen that causes Phytophthora stem and root rot (PRR) worldwide. This study aims to identify temporal changes in P. sojae pathotype complexity, diversity, and Rps gene efficacy. Pathotype data was collected from 5121 isolates of P. sojae, derived from 29 surveys conducted between 1990 and 2019 across the United States, Argentina, Canada, and China. This systematic review shows a loss of efficacy of specific Rps genes utilized for disease management and a significant increase in the pathotype diversity of isolates over time. This study finds that the most widely deployed Rps genes used to manage PRR globally, Rps1a, Rps1c and Rps1k, are no longer effective for PRR management in the United States, Argentina, and Canada. This systematic review emphasizes the need to widely introduce new sources of resistance to P. sojae, such as Rps3a, Rps6, or Rps11, into commercial cultivars to effectively manage PRR going forward.

More than a spiny morphology: plastome variation in the prickly pear cacti (Opuntieae).

BACKGROUND: Plastid genomes (plastomes) have long been recognized as highly conserved in their overall structure, size, gene arrangement and content among land plants. However, recent studies have shown that some lineages present unusual variations in some of these features. Members of the cactus family are one of these lineages, with distinct plastome structures reported across disparate lineages, including gene losses, inversions, boundary movements or loss of the canonical inverted repeat (IR) region. However, only a small fraction of cactus diversity has been analysed so far. METHODS: Here, we investigated plastome features of the tribe Opuntieae, the remarkable prickly pear cacti, which represent one of the most diverse and important lineages of Cactaceae. We assembled de novo the plastome of 43 species, representing a comprehensive sampling of the tribe, including all seven genera, and analysed their evolution in a phylogenetic comparative framework. Phylogenomic analyses with different datasets (full plastome sequences and genes only) were performed, followed by congruence analyses to assess signals underlying contentious nodes. KEY RESULTS: Plastomes varied considerably in length, from 121 to 162 kbp, with striking differences in the content and size of the IR region (contraction and expansion events), including a lack of the canonical IR in some lineages and the pseudogenization or loss of some genes. Overall, nine different types of plastomes were reported, deviating in the presence of the IR region or the genes contained in the IR. Overall, plastome sequences resolved phylogenetic relationships within major clades of Opuntieae with high bootstrap values but presented some contentious nodes depending on the dataset analysed (e.g. whole plastome vs. genes only). Congruence analyses revealed that most plastidial regions lack phylogenetic resolution, while few markers are supporting the most likely topology. Likewise, alternative topologies are driven by a handful of plastome markers, suggesting recalcitrant nodes in the phylogeny. CONCLUSIONS: Our study reveals a dynamic nature of plastome evolution across closely related lineages, shedding light on peculiar features of plastomes. Variation of plastome types across Opuntieae is remarkable in size, structure and content and can be important for the recognition of species in some major clades. Unravelling connections between the causes of plastome variation and the consequences for species biology, physiology, ecology, diversification and adaptation is a promising and ambitious endeavour in cactus research. Although plastome data resolved major phylogenetic relationships, the generation of nuclear genomic data is necessary to confront these hypotheses and assess the recalcitrant nodes further.

Agricultural weeds: the contribution of domesticated species to the origin and evolution of feral weeds.

Agricultural weeds descended from domesticated ancestors, directly from crops (endoferality) and/or from crop-wild hybridization (exoferality), may have evolutionary advantages by rapidly acquiring traits pre-adapted to agricultural habitats. Understanding the role of crops on the origin and evolution of agricultural weeds is essential to develop more effective weed management programs, minimize crop losses due to weeds, and accurately assess the risks of cultivated genes escaping. In this review, we first describe relevant traits of weediness: shattering, seed dormancy, branching, early flowering and rapid growth, and their role in the feralization process. Furthermore, we discuss how the design of "super-crops" can affect weed evolution. We then searched for literature documenting cases of agricultural weeds descended from well-domesticated crops, and describe six case studies of feral weeds evolved from major crops: maize, radish, rapeseed, rice, sorghum, and sunflower. Further studies on the origin and evolution of feral weeds can improve our understanding of the physiological and genetic mechanisms underpinning the adaptation to agricultural habitats and may help to develop more effective weed-control practices and breeding better crops. © 2022 Society of Chemical Industry.

Molecular profiling of bacterial blight resistance in Malaysian rice cultivars.

Bacteria blight is one of the most serious bacterial diseases of rice worldwide. The identification of genetic potential against bacterial blight in the existing rice resources is a prerequisite to develop multigenic resistance to combat the threat of climate change. This investigation was conducted to evaluate alleles variation in 38 Malaysian cultivars using thirteen Simple Sequences Repeats markers and one Sequence Tagged Sites (STS) marker which were reported to be linked with the resistance to bacterial blight. Based on molecular data, a dendrogram was constructed which classified the rice cultivars into seven major clusters at 0.0, 0.28 and 0.3 of similarity coefficient. Cluster 5 was the largest group comprised of ten rice cultivars where multiple genes were identified. However, xa13 could not be detected in the current rice germplasm, whereas xa2 was detected in 25 cultivars. Molecular analysis revealed that Malaysian rice cultivars possess multigenic resistance.

An adaptive teosinte mexicana introgression modulates phosphatidylcholine levels and is associated with maize flowering time.

Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp. parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern United States, Canada, and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time.

Population structure, allelic variation at Rht-B1 and Ppd-A1 loci and its effects on agronomic traits in Argentinian durum wheat.

Exploring the genetic variability in yield and yield-related traits is essential to continue improving genetic gains. Fifty-nine Argentinian durum wheat cultivars were analyzed for important agronomic traits in three field experiments. The collection was genotyped with 3565 genome-wide SNPs and functional markers in order to determine the allelic variation at Rht-B1 and Ppd-A1 genes. Population structure analyses revealed the presence of three main groups, composed by old, modern and genotypes with European or CIMMYT ancestry. The photoperiod sensitivity Ppd-A1b allele showed higher frequency (75%) than the insensitivity one Ppd-A1a (GS105). The semi-dwarfism Rht-B1b and the Ppd-A1a (GS105) alleles were associated with increases in harvest index and decreases in plant height, grain protein content and earlier heading date, although only the varieties carrying the Rht-B1 variants showed differences in grain yield. Out of the two main yield components, grain number per plant was affected by allelic variants at Rht-B1 and Ppd-A1 loci, while no differences were observed in thousand kernel weight. The increases in grain number per spike associated with Rht-B1b were attributed to a higher grain number per spikelet, whereas Ppd-A1a (GS105) was associated with higher grain number per spikelet, but also with lower spikelets per spike.

Characterization of cassava ORANGE proteins and their capability to increase provitamin A carotenoids accumulation.

Cassava (Manihot esculenta Crantz) biofortification with provitamin A carotenoids is an ongoing process that aims to alleviate vitamin A deficiency. The moderate content of provitamin A carotenoids achieved so far limits the contribution to providing adequate dietary vitamin A levels. Strategies to increase carotenoid content focused on genes from the carotenoids biosynthesis pathway. In recent years, special emphasis was given to ORANGE protein (OR), which promotes the accumulation of carotenoids and their stability in several plants. The aim of this work was to identify, characterize and investigate the role of OR in the biosynthesis and stabilization of carotenoids in cassava and its relationship with phytoene synthase (PSY), the rate-limiting enzyme of the carotenoids biosynthesis pathway. Gene and protein characterization of OR, expression levels, protein amounts and carotenoids levels were evaluated in roots of one white (60444) and two yellow cassava cultivars (GM5309-57 and GM3736-37). Four OR variants were found in yellow cassava roots. Although comparable expression was found for three variants, significantly higher OR protein amounts were observed in the yellow varieties. In contrast, cassava PSY1 expression was significantly higher in the yellow cultivars, but PSY protein amount did not vary. Furthermore, we evaluated whether expression of one of the variants, MeOR_X1, affected carotenoid accumulation in cassava Friable Embryogenic Callus (FEC). Overexpression of maize PSY1 alone resulted in carotenoids accumulation and induced crystal formation. Co-expression with MeOR_X1 led to greatly increase of carotenoids although PSY1 expression was high in the co-expressed FEC. Our data suggest that posttranslational mechanisms controlling OR and PSY protein stability contribute to higher carotenoid levels in yellow cassava. Moreover, we showed that cassava FEC can be used to study the efficiency of single and combinatorial gene expression in increasing the carotenoid content prior to its application for the generation of biofortified cassava with enhanced carotenoids levels.

Integrating omics approaches to discover and prioritize candidate genes involved in oil biosynthesis in soybean.

Soybean is a major source of edible protein and oil. Oil content is a quantitative trait that is significantly determined by genetic and environmental factors. Over the past 30 years, a large volume of soybean genetic, genomic, and transcriptomic data have been accumulated. Nevertheless, integrative analyses of such data remain scarce, in spite of their importance for crop improvement. We hypothesized that the co-occurrence of genomic regions for oil-related traits in different studies may reveal more stable regions encompassing important genetic determinants of oil content and quality in soybean. We integrated publicly available data, obtained with distinct techniques, to discover and prioritize candidate genes involved in oil biosynthesis and regulation in soybean. We detected key fatty acid biosynthesis genes (e.g., BCCP2 and ACCase, FADs, KAS family proteins) and several transcription factors, which are likely regulators of oil biosynthesis. In addition, we identified new candidates for seed oil accumulation and quality, such as Glyma.03G213300 and Glyma.19G160700, which encode a translocator protein homolog and a histone acetyltransferase, respectively. Further, oil and protein genomic hotspots are strongly associated with breeding and not with domestication, suggesting that soybean domestication prioritized other traits. The genes identified here are promising targets for breeding programs and for the development of soybean lines with increased oil content and quality.

A novel strategy to uncover specific GO terms/phosphorylation pathways in phosphoproteomic data in Arabidopsis thaliana.

BACKGROUND: Proteins are the workforce of the cell and their phosphorylation status tailors specific responses efficiently. One of the main challenges of phosphoproteomic approaches is to deconvolute biological processes that specifically respond to an experimental query from a list of phosphoproteins. Comparison of the frequency distribution of GO (Gene Ontology) terms in a given phosphoproteome set with that observed in the genome reference set (GenRS) is the most widely used tool to infer biological significance. Yet, this comparison assumes that GO term distribution between the phosphoproteome and the genome are identical. However, this hypothesis has not been tested due to the lack of a comprehensive phosphoproteome database. RESULTS: In this study, we test this hypothesis by constructing three phosphoproteome databases in Arabidopsis thaliana: one based in experimental data (ExpRS), another based in in silico phosphorylation protein prediction (PredRS) and a third that is the union of both (UnRS). Our results show that the three phosphoproteome reference sets show default enrichment of several GO terms compared to GenRS, indicating that GO term distribution in the phosphoproteomes does not match that of the genome. Moreover, these differences overshadow the identification of GO terms that are specifically enriched in a particular condition. To overcome this limitation, we present an additional comparison of the sample of interest with UnRS to uncover GO terms specifically enriched in a particular phosphoproteome experiment. Using this strategy, we found that mRNA splicing and cytoplasmic microtubule compounds are important processes specifically enriched in the phosphoproteome of dark-grown Arabidopsis seedlings. CONCLUSIONS: This study provides a novel strategy to uncover GO specific terms in phosphoproteome data of Arabidopsis that could be applied to any other organism. We also highlight the importance of specific phosphorylation pathways that take place during dark-grown Arabidopsis development.