Gene ontology functional annotation datasets for the ITAG3.2 and ITAG4.0 tomato (Solanum lycopersicum) genome annotations.

Functional annotation based on Gene Ontology has provided a structured and comprehensive system to access the current knowledge about the function of genes. For model plants such as Arabidopsis thaliana, there is a constant updating and restructuring of the functional annotation that increases the reliability of the analyses that use it. For tomato (Solanum lycopersicum), a crop widely used as a model plant for the study of fleshy fruits, there is no functional annotation, at least not freely accessible, even though its genome has long been sequenced and annotated. In this work, we generated, using a simplified version of the maize GAMER pipeline, a tomato Gene Ontology functional annotation with 72.42% (ITAG3.2) and 74.2% (ITAG4.0) of protein-coding genes with at least one GO term association. With this dataset, we share a reliable and easy-to-use tool with the tomato community.

Multi-omic analyses reveal the unique properties of chia (Salvia hispanica) seed metabolism.

Chia (Salvia hispanica) is an emerging crop considered a functional food containing important substances with multiple potential applications. However, the molecular basis of some relevant chia traits, such as seed mucilage and polyphenol content, remains to be discovered. This study generates an improved chromosome-level reference of the chia genome, resolving some highly repetitive regions, describing methylation patterns, and refining genome annotation. Transcriptomic analysis shows that seeds exhibit a unique expression pattern compared to other organs and tissues. Thus, a metabolic and proteomic approach is implemented to study seed composition and seed-produced mucilage. The chia genome exhibits a significant expansion in mucilage synthesis genes (compared to Arabidopsis), and gene network analysis reveals potential regulators controlling seed mucilage production. Rosmarinic acid, a compound with enormous therapeutic potential, was classified as the most abundant polyphenol in seeds, and candidate genes for its complex pathway are described. Overall, this study provides important insights into the molecular basis for the unique characteristics of chia seeds.

Assembly and annotation of the Gossypium barbadense L. 'Pima-S6' genome raise questions about the chromosome structure and gene content of Gossypium barbadense genomes.

BACKGROUND: Gossypium barbadense L. Pima cotton is known for its resistance to Fusarium wilt and for producing fibers of superior quality highly prized in the textile market. We report a high-quality genome assembly and annotation of Pima-S6 cotton and its comparison at the chromosome and protein level to other ten Gossypium published genome assemblies. RESULTS: Synteny and orthogroup analyses revealed important differences on chromosome structure and annotated proteins content between our Pima-S6 and other publicly available G. barbadense assemblies, and across Gossypium assemblies in general. Detailed synteny analyses revealed chromosomal rearrangements between Pima-S6 and other Pima genomes on several chromosomes, with three major inversions in chromosomes A09, A13 and D05, raising questions about the true chromosome structure of Gossypium barbadense genomes. CONCLUSION: Analyses of the re-assembled and re-annotated genome of the close relative G. barbadense Pima 3-79 using our Pima-S6 assembly suggest that contig placement of some recent G. barbadense assemblies might have been unduly influenced by the use of the G. hirsutum TM-1 genome as the anchoring reference. The Pima-S6 reference genome provides a valuable genomic resource and offers new insights on genomic structure, and can serve as G. barbadense genome reference for future assemblies and further support FOV4-related studies and breeding efforts.

Enhanced phenylpropanoid metabolism underlies resistance to Fusarium oxysporum f. sp. vasinfectum race 4 infection in the cotton cultivar Pima-S6 (Gossypium barbadense L.).

Introduction: Fusarium oxysporum f. sp. vasinfectum (FOV) race 4 (FOV4) is a highly pathogenic soil-borne fungus responsible for Fusarium wilt in cotton (Gossypium spp.) and represents a continuing threat to cotton production in the southwest states of the United States, including California, New Mexico, and Texas. Pima (G. barbadense L.) cotton, which is highly valued for its fiber quality, has been shown to be more susceptible to this pathogen than Upland (G. hirsutum L.) cotton. Still, some Pima cultivars present resistance to FOV4 infection. Methods: To gain insights into the FOV4-resistance mechanism, we performed comparative transcriptional and metabolomic analyses between FOV4-susceptible and FOV4-resistant Pima cotton entries. FOV4-resistant Pima-S6 and FOV4-susceptible Pima S-7 and Pima 3-79 cotton plants were infected with FOV4 in the greenhouse, and the roots harvested 11 days post-infection for further analysis. Results: We found that an enhanced root phenylpropanoid metabolism in the resistant Pima-S6 cultivar determines FOV4-resistance. Gene-ontology enrichment of phenylpropanoid biosynthesis and metabolism categories correlated with the accumulation of secondary metabolites in Pima-S6 roots. Specifically, we found esculetin, a coumarin, an inhibitor of Fusarium's growth, accumulated in the roots of Pima-S6 even under non-infected conditions. Genes related to the phenylpropanoid biosynthesis and metabolism, including phenylalanine ammonia-lyase 2 (PAL2) and pleiotropic drug resistance 12 (PDR12) transporter, were found to be upregulated in Pima-S6 roots. Discussion: Our results highlight an essential role for the phenylpropanoid synthesis pathway in FOV4 resistance in Pima-S6 cotton. These genes represent attractive research prospects for FOV4-disease resistance and breeding approaches of other cotton cultivars of economic relevance.

A comparative genomics examination of desiccation tolerance and sensitivity in two sister grass species.

Desiccation tolerance is an ancient and complex trait that spans all major lineages of life on earth. Although important in the evolution of land plants, the mechanisms that underlay this complex trait are poorly understood, especially for vegetative desiccation tolerance (VDT). The lack of suitable closely related plant models that offer a direct contrast between desiccation tolerance and sensitivity has hampered progress. We have assembled high-quality genomes for two closely related grasses, the desiccation-tolerant Sporobolus stapfianus and the desiccation-sensitive Sporobolus pyramidalis Both species are complex polyploids; S. stapfianus is primarily tetraploid, and S. pyramidalis is primarily hexaploid. S. pyramidalis undergoes a major transcriptome remodeling event during initial exposure to dehydration, while S. stapfianus has a muted early response, with peak remodeling during the transition between 1.5 and 1.0 grams of water (gH2O) g-1 dry weight (dw). Functionally, the dehydration transcriptome of S. stapfianus is unrelated to that for S. pyramidalis A comparative analysis of the transcriptomes of the hydrated controls for each species indicated that S. stapfianus is transcriptionally primed for desiccation. Cross-species comparative analyses indicated that VDT likely evolved from reprogramming of desiccation tolerance mechanisms that evolved in seeds and that the tolerance mechanism of S. stapfianus represents a recent evolution for VDT within the Chloridoideae. Orthogroup analyses of the significantly differentially abundant transcripts reconfirmed our present understanding of the response to dehydration, including the lack of an induction of senescence in resurrection angiosperms. The data also suggest that failure to maintain protein structure during dehydration is likely critical in rendering a plant desiccation sensitive.

Defective cytokinin signaling reprograms lipid and flavonoid gene-to-metabolite networks to mitigate high salinity in Arabidopsis.

Cytokinin (CK) in plants regulates both developmental processes and adaptation to environmental stresses. Arabidopsis histidine phosphotransfer ahp2,3,5 and type-B Arabidopsis response regulator arr1,10,12 triple mutants are almost completely defective in CK signaling, and the ahp2,3,5 mutant was reported to be salt tolerant. Here, we demonstrate that the arr1,10,12 mutant is also more tolerant to salt stress than wild-type (WT) plants. A comprehensive metabolite profiling coupled with transcriptome analysis of the ahp2,3,5 and arr1,10,12 mutants was conducted to elucidate the salt tolerance mechanisms mediated by CK signaling. Numerous primary (e.g., sugars, amino acids, and lipids) and secondary (e.g., flavonoids and sterols) metabolites accumulated in these mutants under nonsaline and saline conditions, suggesting that both prestress and poststress accumulations of stress-related metabolites contribute to improved salt tolerance in CK-signaling mutants. Specifically, the levels of sugars (e.g., trehalose and galactinol), amino acids (e.g., branched-chain amino acids and γ-aminobutyric acid), anthocyanins, sterols, and unsaturated triacylglycerols were higher in the mutant plants than in WT plants. Notably, the reprograming of flavonoid and lipid pools was highly coordinated and concomitant with the changes in transcriptional levels, indicating that these metabolic pathways are transcriptionally regulated by CK signaling. The discovery of the regulatory role of CK signaling on membrane lipid reprogramming provides a greater understanding of CK-mediated salt tolerance in plants. This knowledge will contribute to the development of salt-tolerant crops with the ability to withstand salinity as a key driver to ensure global food security in the era of climate crisis.

Whole-genome characterization and comparative genomics of a novel freshwater cyanobacteria species: Pseudanabaena punensis.

Cyanobacteria are emerging as a potential source of novel, beneficial bioactive compounds. However, some cyanobacteria species can harm water quality and public health through the production of toxins. Therefore, surveying the occurrence and generating genomic resources of cyanobacteria producing harmful compounds could help develop the control methods necessary to manage their growth and limit the release contaminants into the water bodies. Here, we describe a novel strain, Pseudanabaena punensis isolated from the open ends of pipelines supplying freshwater. This isolate was characterized morphologically, biochemically and by whole-genome sequence analysis. We also provide genomic information for P. punensis to help understand and highlight the features unique to this isolate. Morphological and genetic (analysis using 16S rRNA and rbcL genes) data were used to assign this novel strain to phylogenetic and taxonomic groups. The isolate was identified as a filamentous and non-heterocystous cyanobacteria. Based on morphological and 16S rRNA phylogeny, this isolate shares characteristics with the Pseudanabaenaceae family, but remains distinct from well-characterized species suggesting its polyphyletic assemblage. The whole-genome sequence analysis suggests greater genomic and phenotypic plasticity. Genome-wide sequence and comparative genomic analyses, comparing against several closely related species, revealed diverse and important genes associated with synthesizing bioactive compounds, multi-drug resistance pathway, heavy metal resistance, and virulence factors. This isolate also produces several important fatty acids with potential industrial applications. The observations described in this study emphasize both industrial applications and risks associated with the freshwater contamination, and therefore genomic resources provided in this study offer an opportunity for further investigations.

Low nitrogen availability inhibits the phosphorus starvation response in maize (Zea mays ssp. mays L.).

BACKGROUND: Nitrogen (N) and phosphorus (P) are macronutrients essential for crop growth and productivity. In cultivated fields, N and P levels are rarely sufficient, contributing to the gap between realized and potential production. Fertilizer application increases nutrient availability, but is not available to all farmers, nor are current rates of application sustainable or environmentally desirable. Transcriptomic studies of cereal crops have revealed dramatic responses to either low N or low P single stress treatments. In the field, however, levels of both N and P may be suboptimal. The interaction between N and P starvation responses remains to be fully characterized. RESULTS: We characterized growth and root and leaf transcriptomes of young maize plants under nutrient replete, low N, low P or combined low NP conditions. We identified 1555 genes to respond to our nutrient treatments, in one or both tissues. A large group of genes, including many classical P starvation response genes, were regulated antagonistically between low N and P conditions. An additional experiment over a range of N availability indicated that a mild reduction in N levels was sufficient to repress the low P induction of P starvation genes. Although expression of P transporter genes was repressed under low N or low NP, we confirmed earlier reports of P hyper accumulation under N limitation. CONCLUSIONS: Transcriptional responses to low N or P were distinct, with few genes responding in a similar way to the two single stress treatments. In combined NP stress, the low N response dominated, and the P starvation response was largely suppressed. A mild reduction in N availability was sufficient to repress the induction of P starvation associated genes. We conclude that activation of the transcriptional response to P starvation in maize is contingent on N availability.

Identification of genuine and novel miRNAs in Amaranthus hypochondriacus from high-throughput sequencing data.

Amaranth has been proposed as an exceptional alternative for food security and climate change mitigation. Information about the distribution, abundance, or specificity of miRNAs in amaranth species is scare. Here, small RNAs from seedlings under control, drought, heat, and cold stress conditions of the Amaranthus hypocondriacus variety "Gabriela" were sequenced and miRNA loci identified in the amaranth genome using the ShortStack software. Fifty-three genuine miRNA clustersthirty-nine belonging to conserved families, and fourteen novel, were identified. Identification of their target genes suggests that conserved amaranth miRNAs are involved in growth and developmental processes, as well as stress responses. MiR0005, an amaranth-specific miRNA, exhibited an unusual high level of expression, akin to that of conserved miRNAs. Overall, our results broaden our knowledge regarding the distribution, abundance and expression of miRNAs in amaranth, providing the basis for future research on miRNAs and their functions in this important species.

Assessment of the ptxD gene as a growth and selective marker in Trichoderma atroviride using Pccg6, a novel constitutive promoter.

BACKGROUND: Trichoderma species are among the most effective cell factories to produce recombinant proteins, whose productivity relies on the molecular toolkit and promoters available for the expression of the target protein. Although inducible promoter systems have been developed for producing recombinant proteins in Trichoderma, constitutive promoters are often a desirable alternative. Constitutive promoters are simple to use, do not require external stimuli or chemical inducers to be activated, and lead to purer enzyme preparations. Moreover, most of the promoters for homologous and heterologous expression reported in Trichoderma have been commonly evaluated by directly assessing production of industrial enzymes, requiring optimization of laborious protocols. RESULTS: Here we report the identification of Pccg6, a novel Trichoderma atroviride constitutive promoter, that has similar transcriptional strength as that of the commonly used pki1 promoter. Pccg6 displayed conserved arrangements of transcription factor binding sites between promoter sequences of Trichoderma ccg6 orthologues genes, potentially involved in their regulatory properties. The predicted ccg6-encoded protein potentially belongs to the SPE1/SPI1 protein family and shares high identity with CCG6 orthologue sequences from other fungal species including Trichoderma reesei, Trichoderma virens, Trichoderma asperellum, and to a lesser extent to that of Neurospora crassa. We also report the use of the Pccg6 promoter to drive the expression of PTXD, a phosphite oxidoreductase of bacterial origin, which allowed T. atroviride to utilize phosphite as a sole source of phosphorus. We propose ptxD as a growth reporter gene that allows real-time comparison of the functionality of different promoters by monitoring growth of Trichoderma transgenic lines and enzymatic activity of PTXD. Finally, we show that constitutive expression of ptxD provided T. atroviride a competitive advantage to outgrow bacterial contaminants when supplied with phosphite as a sole source of phosphorus. CONCLUSIONS: A new constitutive promoter, ccg6, for expression of homologous and heterologous proteins has been identified and tested in T. atroviride to express PTXD, which resulted in an effective and visible phenotype to evaluate transcriptional activity of sequence promoters. Use of PTXD as a growth marker holds great potential for assessing activity of other promoters and for biotechnological applications as a contamination control system.

The plant MBF1 protein family: a bridge between stress and transcription.

The Multiprotein Bridging Factor 1 (MBF1) proteins are transcription co-factors whose molecular function is to form a bridge between transcription factors and the basal machinery of transcription. MBF1s are present in most archaea and all eukaryotes, and numerous reports show that they are involved in developmental processes and in stress responses. In this review we summarize almost three decades of research on the plant MBF1 family, which has mainly focused on their role in abiotic stress responses, in particular the heat stress response. However, despite the amount of information available, there are still many questions that remain about how plant MBF1 genes, transcripts, and proteins respond to stress, and how they in turn modulate stress response transcriptional pathways.

Bioinformatic Analysis of Small RNA Sequencing Libraries.

Bioinformatic analysis of small RNA sequencing libraries consists of transforming a series of small RNA sequencing experiment fastq files into a table containing small RNA sequences and their abundance. This is achieved by cleaning the reads, aligning the cleaned reads to a reference, and parsing the alignment results. In this protocol we present the most common option, and the rationale, for each of these steps.

Entering the Next Dimension: Plant Genomes in 3D.

After linear sequences of genomes and epigenomic landscape data, the 3D organization of chromatin in the nucleus is the next level to be explored. Different organisms present a general hierarchical organization, with chromosome territories at the top. Chromatin interaction maps, obtained by chromosome conformation capture (3C)-based methodologies, for eight plant species reveal commonalities, but also differences, among them and with animals. The smallest structures, found in high-resolution maps of the Arabidopsis genome, are single genes. Epigenetic marks (histone modification and DNA methylation), transcriptional activity, and chromatin interaction appear to be correlated, and whether structure is the cause or consequence of the function of interacting regions is being actively investigated.

Allele specific expression analysis identifies regulatory variation associated with stress-related genes in the Mexican highland maize landrace Palomero Toluqueño.

BACKGROUND: Gene regulatory variation has been proposed to play an important role in the adaptation of plants to environmental stress. In the central highlands of Mexico, farmer selection has generated a unique group of maize landraces adapted to the challenges of the highland niche. In this study, gene expression in Mexican highland maize and a reference maize breeding line were compared to identify evidence of regulatory variation in stress-related genes. It was hypothesised that local adaptation in Mexican highland maize would be associated with a transcriptional signature observable even under benign conditions. METHODS: Allele specific expression analysis was performed using the seedling-leaf transcriptome of an F1 individual generated from the cross between the highland adapted Mexican landrace Palomero Toluqueño and the reference line B73, grown under benign conditions. Results were compared with a published dataset describing the transcriptional response of B73 seedlings to cold, heat, salt and UV treatments. RESULTS: A total of 2,386 genes were identified to show allele specific expression. Of these, 277 showed an expression difference between Palomero Toluqueño and B73 alleles under benign conditions that anticipated the response of B73 cold, heat, salt and/or UV treatments, and, as such, were considered to display a prior stress response. Prior stress response candidates included genes associated with plant hormone signaling and a number of transcription factors. Construction of a gene co-expression network revealed further signaling and stress-related genes to be among the potential targets of the transcription factors candidates. DISCUSSION: Prior activation of responses may represent the best strategy when stresses are severe but predictable. Expression differences observed here between Palomero Toluqueño and B73 alleles indicate the presence of cis-acting regulatory variation linked to stress-related genes in Palomero Toluqueño. Considered alongside gene annotation and population data, allele specific expression analysis of plants grown under benign conditions provides an attractive strategy to identify functional variation potentially linked to local adaptation.

The bHLH transcription factor SPATULA enables cytokinin signaling, and both activate auxin biosynthesis and transport genes at the medial domain of the gynoecium.

Fruits and seeds are the major food source on earth. Both derive from the gynoecium and, therefore, it is crucial to understand the mechanisms that guide the development of this organ of angiosperm species. In Arabidopsis, the gynoecium is composed of two congenitally fused carpels, where two domains: medial and lateral, can be distinguished. The medial domain includes the carpel margin meristem (CMM) that is key for the production of the internal tissues involved in fertilization, such as septum, ovules, and transmitting tract. Interestingly, the medial domain shows a high cytokinin signaling output, in contrast to the lateral domain, where it is hardly detected. While it is known that cytokinin provides meristematic properties, understanding on the mechanisms that underlie the cytokinin signaling pattern in the young gynoecium is lacking. Moreover, in other tissues, the cytokinin pathway is often connected to the auxin pathway, but we also lack knowledge about these connections in the young gynoecium. Our results reveal that cytokinin signaling, that can provide meristematic properties required for CMM activity and growth, is enabled by the transcription factor SPATULA (SPT) in the medial domain. Meanwhile, cytokinin signaling is confined to the medial domain by the cytokinin response repressor ARABIDOPSIS HISTIDINE PHOSPHOTRANSFERASE 6 (AHP6), and perhaps by ARR16 (a type-A ARR) as well, both present in the lateral domains (presumptive valves) of the developing gynoecia. Moreover, SPT and cytokinin, probably together, promote the expression of the auxin biosynthetic gene TRYPTOPHAN AMINOTRANSFERASE OF ARABIDOPSIS 1 (TAA1) and the gene encoding the auxin efflux transporter PIN-FORMED 3 (PIN3), likely creating auxin drainage important for gynoecium growth. This study provides novel insights in the spatiotemporal determination of the cytokinin signaling pattern and its connection to the auxin pathway in the young gynoecium.

Phosphorus acquisition efficiency in arbuscular mycorrhizal maize is correlated with the abundance of root-external hyphae and the accumulation of transcripts encoding PHT1 phosphate transporters.

Plant interactions with arbuscular mycorrhizal fungi have long attracted interest for their potential to promote more efficient use of mineral resources in agriculture. Their use, however, remains limited by a lack of understanding of the processes that determine the outcome of the symbiosis. In this study, the impact of host genotype on growth response to mycorrhizal inoculation was investigated in a panel of diverse maize lines. A panel of 30 maize lines was evaluated with and without inoculation with arbuscular mycorrhizal fungi. The line Oh43 was identified to show superior response and, along with five other reference lines, was characterized in greater detail in a split-compartment system, using 33 P to quantify mycorrhizal phosphorus uptake. Changes in relative growth indicated variation in host capacity to profit from the symbiosis. Shoot phosphate content, abundance of root-internal and -external fungal structures, mycorrhizal phosphorus uptake, and accumulation of transcripts encoding plant PHT1 family phosphate transporters varied among lines. Superior response in Oh43 is correlated with extensive development of root-external hyphae, accumulation of specific Pht1 transcripts and high phosphorus uptake by mycorrhizal plants. The data indicate that host genetic factors influence fungal growth strategy with an impact on plant performance.

Selection of Reference Genes for Quantitative Real-Time RT-PCR Studies in Tomato Fruit of the Genotype MT-Rg1.

Quantitative real-time RT-PCR (qRT-PCR) has become one of the most widely used methods for accurate quantification of gene expression. Since there are no universal reference genes for normalization, the optimal strategy to normalize raw qRT-PCR data is to perform an initial comparison of a set of independent reference genes to assess the most stable ones in each biological model. Normalization of a qRT-PCR experiment helps to ensure that the results are both statistically significant and biologically meaningful. Tomato is the model of choice to study fleshy fruit development. The miniature tomato (Solanum lycopersicum L.) cultivar Micro-Tom (MT) is considered a model system for tomato genetics and functional genomics. A new genotype, containing the Rg1 allele, improves tomato in vitro regeneration. In this work, we evaluated the expression stability of four tomato reference genes, namely CAC, SAND, Expressed, and ACTIN2. We showed that the genes CAC and Exp are the best reference genes of the four we tested during fruit development in the MT-Rg1 genotype. Furthermore, we validated the reference genes by showing that the expression profiles of the transcription factors FRUITFULL1 and APETALA2c during fruit development are comparable to previous reports using other tomato cultivars.

Regulatory network analysis reveals novel regulators of seed desiccation tolerance in Arabidopsis thaliana.

Desiccation tolerance (DT) is a remarkable process that allows seeds in the dry state to remain viable for long periods of time that in some instances exceed 1,000 y. It has been postulated that seed DT evolved by rewiring the regulatory and signaling networks that controlled vegetative DT, which itself emerged as a crucial adaptive trait of early land plants. Understanding the networks that regulate seed desiccation tolerance in model plant systems would provide the tools to understand an evolutionary process that played a crucial role in the diversification of flowering plants. In this work, we used an integrated approach that included genomics, bioinformatics, metabolomics, and molecular genetics to identify and validate molecular networks that control the acquisition of DT in Arabidopsis seeds. Two DT-specific transcriptional subnetworks were identified related to storage of reserve compounds and cellular protection mechanisms that act downstream of the embryo development master regulators LEAFY COTYLEDON 1 and 2, FUSCA 3, and ABSCICIC ACID INSENSITIVE 3. Among the transcription factors identified as major nodes in the DT regulatory subnetworks, PLATZ1, PLATZ2, and AGL67 were confirmed by knockout mutants and overexpression in a desiccation-intolerant mutant background to play an important role in seed DT. Additionally, we found that constitutive expression of PLATZ1 in WT plants confers partial DT in vegetative tissues.

Altered expression of the bZIP transcription factor DRINK ME affects growth and reproductive development in Arabidopsis thaliana.

Here we describe an uncharacterized gene that negatively influences Arabidopsis growth and reproductive development. DRINK ME (DKM; bZIP30) is a member of the bZIP transcription factor family, and is expressed in meristematic tissues such as the inflorescence meristem (IM), floral meristem (FM), and carpel margin meristem (CMM). Altered DKM expression affects meristematic tissues and reproductive organ development, including the gynoecium, which is the female reproductive structure and is determinant for fertility and sexual reproduction. A microarray analysis indicates that DKM overexpression affects the expression of cell cycle, cell wall, organ initiation, cell elongation, hormone homeostasis, and meristem activity genes. Furthermore, DKM can interact in yeast and in planta with proteins involved in shoot apical meristem maintenance such as WUSCHEL, KNAT1/BP, KNAT2 and JAIBA, and with proteins involved in medial tissue development in the gynoecium such as HECATE, BELL1 and NGATHA1. Taken together, our results highlight the relevance of DKM as a negative modulator of Arabidopsis growth and reproductive development.