Brassica rapa CURLY LEAF is a major H3K27 methyltransferase regulating flowering time.

MAIN CONCLUSION: In Brassica rapa, the epigenetic modifier BraA.CLF orchestrates flowering by modulating H3K27me3 levels at the floral integrator genes FT, SOC1, and SEP3, thereby influencing their expression. CURLY LEAF (CLF) is the catalytic subunit of the plant Polycomb Repressive Complex 2 that mediates the trimethylation of histone H3 lysine 27 (H3K27me3), an epigenetic modification that leads to gene silencing. While the function of CURLY LEAF (CLF) has been extensively studied in Arabidopsis thaliana, its role in Brassica crops is barely known. In this study, we focused on the Brassica rapa homolog of CLF and found that the loss-of-function mutant braA.clf-1 exhibits an accelerated flowering together with pleiotropic phenotypic alterations compared to wild-type plants. In addition, we carried out transcriptomic and H3K27me3 genome-wide analyses to identify the genes regulated by BraA.CLF. Interestingly, we observed that several floral regulatory genes, including the B. rapa homologs of FT, SOC1 and SEP3, show reduced H3K27me3 levels and increased transcript levels compared to wild-type plants, suggesting that they are direct targets of BraA.CLF and key players in regulating flowering time in this crop. In addition, the results obtained will enhance our understanding of the epigenetic mechanisms regulating key developmental traits and will aid to increase crop yield by engineering new Brassica varieties with different flowering time requirements.

Functional Profiling of Soft Tissue Sarcoma Using Mechanistic Models.

Soft tissue sarcoma is an umbrella term for a group of rare cancers that are difficult to treat. In addition to surgery, neoadjuvant chemotherapy has shown the potential to downstage tumors and prevent micrometastases. However, finding effective therapeutic targets remains a research challenge. Here, a previously developed computational approach called mechanistic models of signaling pathways has been employed to unravel the impact of observed changes at the gene expression level on the ultimate functional behavior of cells. In the context of such a mechanistic model, RNA-Seq counts sourced from the Recount3 resource, from The Cancer Genome Atlas (TCGA) Sarcoma project, and non-diseased sarcomagenic tissues from the Genotype-Tissue Expression (GTEx) project were utilized to investigate signal transduction activity through signaling pathways. This approach provides a precise view of the relationship between sarcoma patient survival and the signaling landscape in tumors and their environment. Despite the distinct regulatory alterations observed in each sarcoma subtype, this study identified 13 signaling circuits, or elementary sub-pathways triggering specific cell functions, present across all subtypes, belonging to eight signaling pathways, which served as predictors for patient survival. Additionally, nine signaling circuits from five signaling pathways that highlighted the modifications tumor samples underwent in comparison to normal tissues were found. These results describe the protective role of the immune system, suggesting an anti-tumorigenic effect in the tumor microenvironment, in the process of tumor cell detachment and migration, or the dysregulation of ion homeostasis. Also, the analysis of signaling circuit intermediary proteins suggests multiple strategies for therapy.

Conserved and distinct roles of H3K27me3 demethylases regulating flowering time in Brassica rapa.

Epigenetic regulation is necessary for optimal organism development and preservation of gene expression profiles in the cell. In plants, the trimethylation of histone H3 lysine 27 (H3K27me3) is a silencing epigenetic mark relevant for developmental transitions like flowering. The floral transition is a key agronomic trait; however, the epigenetic mechanisms of flowering time regulation in crops remain poorly understood. Here we study the Jumonji H3K27me3 demethylases BraA.REF6 and BraA.ELF6 in Brassica rapa. Phenotypic characterization of novel mutant lines and genome-wide H3K27me3 chromatin immunoprecipitation and transcriptomic analyses indicated that BraA.REF6 plays a greater role than BraA.ELF6 in fine-tuning H3K27me3 levels. In addition, we found that braA.elf6 mutants were early flowering due to high H3K27me3 levels at B. rapa homologs of the floral repressor FLC. Unlike mutations in Arabidopsis thaliana, braA.ref6 mutants were late flowering without altering the expression of B. rapa FLC genes. Remarkably, we found that BraA.REF6 regulated a number of gibberellic acid (GA) biosynthetic genes, including a homolog of GA1, and that GA-treatment complemented the late flowering mutant phenotype. This study increases our understanding of the epigenetic regulation of flowering time in B. rapa, highlighting conserved and distinct regulatory mechanisms between model and crop species.

Regulation of gene expression in roots of the pH-sensitive Vaccinium corymbosum and the pH-tolerant Vaccinium arboreum in response to near neutral pH stress using RNA-Seq.

BACKGROUND: Blueberries are one of the few horticultural crops adapted to grow in acidic soils. Neutral to basic soil pH is detrimental to all commonly cultivated blueberry species, including Vaccinium corymbosum (VC). In contrast, the wild species V. arboreum (VA) is able to tolerate a wider range of soil pH. To assess the molecular mechanisms involved in near neutral pH stress response, plants from pH-sensitive VC (tetraploid) and pH-tolerant VA (diploid) were grown at near neutral pH 6.5 and at the preferred pH of 4.5. RESULTS: Transcriptome sequencing of root RNA was performed for 4 biological replications per species x pH level interaction, for a total of 16 samples. Reads were mapped to the reference genome from diploid V. corymbosum, transforming ~55% of the reads to gene counts. A quasi-likelihood F test identified differential expression due to pH stress in 337 and 4867 genes in VA and VC, respectively. Both species shared regulation of genes involved in nutrient homeostasis and cell wall metabolism. VA and VC exhibited differential regulation of signaling pathways related to abiotic/biotic stress, cellulose and lignin biosynthesis, and nutrient uptake. CONCLUSIONS: The specific responses in VA likely facilitate tolerance to higher soil pH. In contrast, response in VC, despite affecting a greater number of genes, is not effective overcoming the stress induced by pH. Further inspection of those genes with differential expression that are specific in VA may provide insight on the mechanisms towards tolerance.

Chemical profile and analysis of biosynthetic pathways and genes of volatile terpenes in Pityopsis ruthii, a rare and endangered flowering plant.

It is critical to gather biological information about rare and endangered plants to incorporate into conservation efforts. The secondary metabolism of Pityopsis ruthii, an endangered flowering plant that only occurs along limited sections of two rivers (Ocoee and Hiwassee) in Tennessee, USA was studied. Our long-term goal is to understand the mechanisms behind P. ruthii's adaptation to restricted areas in Tennessee. Here, we profiled the secondary metabolites, specifically in flowers, with a focus on terpenes, aiming to uncover the genomic and molecular basis of terpene biosynthesis in P. ruthii flowers using transcriptomic and biochemical approaches. By comparative profiling of the nonpolar portion of metabolites from various tissues, P. ruthii flowers were rich in terpenes, which included 4 monoterpenes and 10 sesquiterpenes. These terpenes were emitted from flowers as volatiles with monoterpenes and sesquiterpenes accounting for almost 68% and 32% of total emission of terpenes, respectively. These findings suggested that floral terpenes play important roles for the biology and adaptation of P. ruthii to its limited range. To investigate the biosynthesis of floral terpenes, transcriptome data for flowers were produced and analyzed. Genes involved in the terpene biosynthetic pathway were identified and their relative expressions determined. Using this approach, 67 putative terpene synthase (TPS) contigs were detected. TPSs in general are critical for terpene biosynthesis. Seven full-length TPS genes encoding putative monoterpene and sesquiterpene synthases were cloned and functionally characterized. Three catalyzed the biosynthesis of sesquiterpenes and four catalyzed the biosynthesis of monoterpenes. In conclusion, P. ruthii plants employ multiple TPS genes for the biosynthesis of a mixture of floral monoterpenes and sesquiterpenes, which probably play roles in chemical defense and attracting insect pollinators alike.

Species diversity and phylogeography of Cornus kousa (Asian dogwood) captured by genomic and genic microsatellites.

Cornus kousa (Asian dogwood), an East Asia native tree, is the most economically important species of the dogwood genus, owing to its desirable horticultural traits and ability to hybridize with North America-native dogwoods. To assess the species genetic diversity and to better inform the ongoing and future breeding efforts, we assembled an herbarium and arboretum collection of 131 noncultivated C. kousa specimens. Genotyping and capillary electrophoresis analyses of our C. kousa collection with the newly developed genic and published nuclear genomic microsatellites permitted assessment of genetic diversity and evolutionary history of the species. Regardless of the microsatellite type used, the study yielded generally similar insights into the C. kousa diversity with subtle differences deriving from and underlining the marker used. The accrued evidence pointed to the species distinct genetic pools related to the plant country of origin. This can be helpful in the development of the commercial cultivars for this important ornamental crop with increased pyramided utility traits. Analyses of the C. kousa evolutionary history using the accrued genotyping datasets pointed to an unsampled ancestor population, possibly now extinct, as per the phylogeography of the region. To our knowledge, there are few studies utilizing the same gDNA collection to compare performance of genomic and genic microsatellites. This is the first detailed report on C. kousa species diversity and evolutionary history inference.

Genome-wide analysis of the H3K27me3 epigenome and transcriptome in Brassica rapa.

BACKGROUND: Genome-wide maps of histone modifications have been obtained for several plant species. However, most studies focus on model systems and do not enforce FAIR data management principles. Here we study the H3K27me3 epigenome and associated transcriptome of Brassica rapa, an important vegetable cultivated worldwide. FINDINGS: We performed H3K27me3 chromatin immunoprecipitation followed by high-throughput sequencing and transcriptomic analysis by 3'-end RNA sequencing from B. rapa leaves and inflorescences. To analyze these data we developed a Reproducible Epigenomic Analysis pipeline using Galaxy and Jupyter, packaged into Docker images to facilitate transparency and reuse. We found that H3K27me3 covers roughly one-third of all B. rapa protein-coding genes and its presence correlates with low transcript levels. The comparative analysis between leaves and inflorescences suggested that the expression of various floral regulatory genes during development depends on H3K27me3. To demonstrate the importance of H3K27me3 for B. rapa development, we characterized a mutant line deficient in the H3K27 methyltransferase activity. We found that braA.clf mutant plants presented pleiotropic alterations, e.g., curly leaves due to increased expression and reduced H3K27me3 levels at AGAMOUS-like loci. CONCLUSIONS: We characterized the epigenetic mark H3K27me3 at genome-wide levels and provide genetic evidence for its relevance in B. rapa development. Our work reveals the epigenomic landscape of H3K27me3 in B. rapa and provides novel genomics datasets and bioinformatics analytical resources. We anticipate that this work will lead the way to further epigenomic studies in the complex genome of Brassica crops.

Taraxacum kok-saghyz (rubber dandelion) genomic microsatellite loci reveal modest genetic diversity and cross-amplify broadly to related species.

Taraxacum kok-saghyz (TKS) carries great potential as alternative natural rubber source. To better inform future breeding efforts with TKS and gain a deeper understanding of its genetic diversity, we utilized de novo sequencing to generate novel genomic simple sequence repeats markers (gSSRs). We utilized 25 gSSRs on a collection of genomic DNA (gDNA) samples from germplasm bank, and two gDNA samples from historical herbarium specimens. PCR coupled with capillary electrophoresis and an array of population genetics tools were employed to analyze the dataset of our study as well as a dataset of the recently published genic SSRs (eSSRs) generated on the same germplasm. Our results using both gSSRs and eSSRs revealed that TKS has low- to- moderate genetic diversity with most of it partitioned to the individuals and individuals within populations, whereas the species lacked population structure. Nineteen of the 25 gSSR markers cross-amplified to other Taraxacum spp. collected from Southeastern United States and identified as T. officinale by ITS sequencing. We used a subset of 14 gSSRs to estimate the genetic diversity of the T. officinale gDNA collection. In contrast to the obligatory outcrossing TKS, T. officinale presented evidence for population structure and clonal reproduction, which agreed with the species biology. We mapped the molecular markers sequences from this study and several others to the well-annotated sunflower genome. Our gSSRs present a functional tool for the biodiversity analyses in Taraxacum, but also in the related genera, as well as in the closely related tribes of the Asteraceae.

Comprehensive evaluation of RNA-seq analysis pipelines in diploid and polyploid species.

Background: The usual analysis of RNA sequencing (RNA-seq) reads is based on an existing reference genome and annotated gene models. However, when a reference for the sequenced species is not available, alternatives include using a reference genome from a related species or reconstructing transcript sequences with de novo assembly. In addition, researchers are faced with many options for RNA-seq data processing and limited information on how their decisions will impact the final outcome. Using both a diploid and polyploid species with a distant reference genome, we have tested the influence of different tools at various steps of a typical RNA-seq analysis workflow on the recovery of useful processed data available for downstream analysis. Findings: At the preprocessing step, we found error correction has a strong influence on de novo assembly but not on mapping results. After trimming, a greater percentage of reads could be used in downstream analysis by selecting gentle quality trimming performed with Skewer instead of strict quality trimming with Trimmomatic. This availability of reads correlated with size, quality, and completeness of de novo assemblies and with number of mapped reads. When selecting a reference genome from a related species to map reads, outcome was significantly improved when using mapping software tolerant of greater sequence divergence, such as Stampy or GSNAP. Conclusions: The selection of bioinformatic software tools for RNA-seq data analysis can maximize quality parameters on de novo assemblies and availability of reads in downstream analysis.

Sunflower HaGPAT9-1 is the predominant GPAT during seed development.

In oil crops, triacylglycerol biosynthesis is an important metabolic pathway in which glycerol-3-phosphate acyltransferase (GPAT) performs the first acylation step. Mass spectrometry analysis of developing sunflower (Helianthus annuus) seed membrane fractions identified an abundant GPAT, HaGPAT9 isoform 1, with a N-terminal peptide that possessed two phosphorylated residues with possible regulatory function. HaGPAT9-1 belongs to a broad eukaryotic GPAT family, similar to mammalian GPAT3, and it represents one of the two sunflower GPAT9 isoforms, sharing 90% identity with HaGPAT9-2. Both sunflower genes are expressed during seed development and in vegetative tissues, with HaGPAT9-1 transcripts accumulating at relatively higher levels than those for HaGPAT9-2. Green fluorescent protein tagging of HaGPAT9-1 confirmed its subcellular accumulation in the endoplasmic reticulum. Despite their overall sequence similarities, the two sunflower isoforms displayed significant differences in their enzymatic activities. For instance, HaGPAT9-1 possesses in vivo GPAT activity that rescues the lethal phenotype of the cmy228 yeast strain, while in vitro assays revealed a preference of HaGPAT9-1 for palmitoyl-, oleoyl- and linoleoyl-CoAs of one order of magnitude, with the highest increase in yield for oleoyl- and linoleoyl-CoAs. By contrast, no enzymatic activity could be detected for HaGPAT9-2, even though its over-expression modified the TAG profile of yeast.

Cloning, heterologous expression and biochemical characterization of plastidial sn-glycerol-3-phosphate acyltransferase from Helianthus annuus.

The acyl-[acyl carrier protein]:sn-1-glycerol-3-phosphate acyltransferase (GPAT; E.C. 2.3.1.15) catalyzes the first step of glycerolipid assembly within the stroma of the chloroplast. In the present study, the sunflower (Helianthus annuus, L.) stromal GPAT was cloned, sequenced and characterized. We identified a single ORF of 1344base pairs that encoded a GPAT sharing strong sequence homology with the plastidial GPAT from Arabidopsis thaliana (ATS1, At1g32200). Gene expression studies showed that the highest transcript levels occurred in green tissues in which chloroplasts are abundant. The corresponding mature protein was heterologously overexpressed in Escherichia coli for purification and biochemical characterization. In vitro assays using radiolabelled acyl-ACPs and glycerol-3-phosphate as substrates revealed a strong preference for oleic versus palmitic acid, and weak activity towards stearic acid. The positional fatty acid composition of relevant chloroplast phospholipids from sunflower leaves did not reflect the in vitro GPAT specificity, suggesting a more complex scenario with mixed substrates at different concentrations, competition with other acyl-ACP consuming enzymatic reactions, etc. In summary, this study has confirmed the affinity of this enzyme which would partly explain the resistance to cold temperatures observed in sunflower plants.