Regulation of regeneration in Arabidopsis thaliana

We employed several algorithms with high efficacy to analyze the public transcriptomic data, aiming to identify key transcription factors (TFs) that regulate regeneration in Arabidopsis thaliana. Initially, we utilized CollaborativeNet, also known as TF-Cluster, to construct a collaborative network of all TFs, which was subsequently decomposed into many subnetworks using the Triple-Link and Compound Spring Embedder (CoSE) algorithms. Functional analysis of these subnetworks led to the identification of nine subnetworks closely associated with regeneration. We further applied principal component analysis and gene ontology (GO) enrichment analysis to reduce the subnetworks from nine to three, namely subnetworks 1, 12, and 17. Searching for TF-binding sites in the promoters of the co-expressed and co-regulated (CCGs) genes of all TFs in these three subnetworks and Triple-Gene Mutual Interaction analysis of TFs in these three subnetworks with the CCGs involved in regeneration enabled us to rank the TFs in each subnetwork. Finally, six potential candidate TFs-WOX9A, LEC2, PGA37, WIP5, PEI1, and AIL1 from subnetwork 1-were identified, and their roles in somatic embryogenesis (GO:0010262) and regeneration (GO:0031099) were discussed, so were the TFs in Subnetwork 12 and 17 associated with regeneration. The TFs identified were also assessed using the CIS-BP database and Expression Atlas. Our analyses suggest some novel TFs that may have regulatory roles in regeneration and embryogenesis and provide valuable data and insights into the regulatory mechanisms related to regeneration. The tools and the procedures used here are instrumental for analyzing high-throughput transcriptomic data and advancing our understanding of the regulation of various biological processes of interest.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-023-00121-9.

Epitranscriptomic subtyping, visualization, and denoising by global motif visualization

Advances in sequencing technologies have empowered epitranscriptomic profiling at the single-base resolution. Putative RNA modification sites identified from a single high-throughput experiment may contain one type of modification deposited by different writers or different types of modifications, along with false positive results because of the challenge of distinguishing signals from noise. However, current tools are insufficient for subtyping, visualization, and denoising these signals. Here, we present iMVP, which is an interactive framework for epitranscriptomic analysis with a nonlinear dimension reduction technique and density-based partition. As exemplified by the analysis of mRNA m5C and ModTect variant data, we show that iMVP allows the identification of previously unknown RNA modification motifs and writers and the discovery of false positives that are undetectable by traditional methods. Using putative m6A/m6Am sites called from 8 profiling approaches, we illustrate that iMVP enables comprehensive comparison of different approaches and advances our understanding of the difference and pattern of true positives and artifacts in these methods. Finally, we demonstrate the ability of iMVP to analyze an extremely large human A-to-I editing dataset that was previously unmanageable. Our work provides a general framework for the visualization and interpretation of epitranscriptomic data.

Evolutionary conservation of maternal RNA localization in fishes and amphibians revealed by TOMO-Seq

Asymmetrical localization of biomolecules inside the egg, results in uneven cell division and establishment of many biological processes, cell types and the body plan. However, our knowledge about evolutionary conservation of localized transcripts is still limited to a few models. Our goal was to compare localization profiles along the animal-vegetal axis of mature eggs from four vertebrate models, two amphibians (Xenopus laevis, Ambystoma mexicanum) and two fishes (Acipenser ruthenus, Danio rerio) using the spatial expression method called TOMO-Seq. We revealed that RNAs of many known important transcripts such as germ layer determinants, germ plasm factors and members of key signalling pathways, are localized in completely different profiles among the models. It was also observed that there was a poor correlation between the vegetally localized transcripts but a relatively good correlation between the animally localized transcripts. These findings indicate that the regulation of embryonic development within the animal kingdom is highly diverse and cannot be deduced based on a single model.

The two clock proteins CCA1 and LHY activate VIN3 transcription during vernalization through the vernalization-responsive cis-element

Vernalization, a long-term cold-mediated acquisition of flowering competence, is critically regulated by VERNALIZATION INSENSITIVE 3 (VIN3), a gene induced by vernalization in Arabidopsis. Although the function of VIN3 has been extensively studied, how VIN3 expression itself is upregulated by long-term cold is not well understood. In this study, we identified a vernalization-responsive cis-element in the VIN3 promoter, VREVIN3, composed of a G-box and an evening element (EE). Mutations in either the G-box or the EE prevented VIN3 expression from being fully induced upon vernalization, leading to defects in the vernalization response. We determined that the core clock proteins CIRCADIAN CLOCK-ASSOCIATED 1 (CCA1) and LATE-ELONGATED HYPOCOTYL (LHY) associate with the EE of VREVIN3, both in vitro and in vivo. In a cca1 lhy double mutant background harboring a functional FRIGIDA allele, long-term cold-mediated VIN3 induction and acceleration of flowering were impaired, especially under mild cold conditions such as at 12°C. During prolonged cold exposure, oscillations of CCA1/LHY transcripts were altered, while CCA1 abundance increased at dusk, coinciding with the diurnal peak of VIN3 transcripts. We propose that modulation of the clock proteins CCA1 and LHY participates in the systems involved in sensing long-term cold for the activation of VIN3 transcription.

TSSFinder-fast and accurate ab initio prediction of the core promoter in eukaryotic genomes

Promoter annotation is an important task in the analysis of a genome. One of the main challenges for this task is locating the border between the promoter region and the transcribing region of the gene, the transcription start site (TSS). The TSS is the reference point to delimit the DNA sequence responsible for the assembly of the transcribing complex. As the same gene can have more than one TSS, so to delimit the promoter region, it is important to locate the closest TSS to the site of the beginning of the translation. This paper presents TSSFinder, a new software for the prediction of the TSS signal of eukaryotic genes that is significantly more accurate than other available software. We currently are the only application to offer pre-trained models for six different eukaryotic organisms: Arabidopsis thaliana, Drosophila melanogaster, Gallus gallus, Homo sapiens, Oryza sativa and Saccharomyces cerevisiae. Additionally, our software can be easily customized for specific organisms using only 125 DNA sequences with a validated TSS signal and corresponding genomic locations as a training set. TSSFinder is a valuable new tool for the annotation of genomes. TSSFinder source code and docker container can be downloaded from http://tssfinder.github.io. Alternatively, TSSFinder is also available as a web service at http://sucest-fun.org/wsapp/tssfinder/.

NoPeak: k-mer-based motif discovery in ChIP-Seq data without peak calling

MOTIVATION

The discovery of sequence motifs mediating DNA-protein binding usually implies the determination of binding sites using high-throughput sequencing and peak calling. The determination of peaks, however, depends strongly on data quality and is susceptible to noise.

RESULTS

Here, we present a novel approach to reliably identify transcription factor-binding motifs from ChIP-Seq data without peak detection. By evaluating the distributions of sequencing reads around the different k-mers in the genome, we are able to identify binding motifs in ChIP-Seq data that yield no results in traditional pipelines.

AVAILABILITY AND IMPLEMENTATION

NoPeak is published under the GNU General Public License and available as a standalone console-based Java application at https://github.com/menzel/nopeak.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Different Cis-Regulatory Elements Control the Tissue-Specific Contribution of Plastid ω-3 Desaturases to Wounding and Hormone Responses

Trienoic fatty acids are essential constituents of biomembranes and precursors of jasmonates involved in plant defense responses. Two ω-3 desaturases, AtFAD7 and AtFAD8, synthetize trienoic fatty acids in the plastid. Promoter:GUS and mutagenesis analysis was used to identify cis-elements controlling AtFAD7 and AtFAD8 basal expression and their response to hormones or wounding. AtFAD7 promoter GUS activity was much higher than that of AtFAD8 in leaves, with specific AtFAD7 expression in the flower stamen and pistil and root meristem and vasculature. This specific tissue and organ expression of AtFAD7 was controlled by different cis-elements. Thus, promoter deletion and mutagenesis analysis indicated that WRKY proteins might be essential for basal expression of AtFAD7 in leaves. Two MYB target sequences present in the AtFAD7 promoter might be responsible for its expression in the flower stamen and stigma of the pistil and in the root meristem, and for the AtFAD7 wound-specific response. Two MYB target sequences detected in the distal region of the AtFAD8 gene promoter seemed to negatively control AtFAD8 expression, particularly in true leaves and flowers, suggesting that MYB transcription factors act as repressors of AtFAD8 gene basal expression, modulating the different relative abundance of both plastid ω-3 desaturases at the transcriptional level. Our data showed that the two ABA repression sequences detected in the AtFAD7 promoter were functional, suggesting an ABA-dependent mechanism involved in the different regulation of both ω-3 plastid desaturases. These results reveal the implication of different signaling pathways for the concerted regulation of trienoic fatty acid content in Arabidopsis.

Amino Acid and Carbohydrate Metabolism Are Coordinated to Maintain Energetic Balance during Drought in Sugarcane

The ability to expand crop plantations without irrigation is a major goal to increase agriculture sustainability. To achieve this end, we need to understand the mechanisms that govern plant growth responses under drought conditions. In this study, we combined physiological, transcriptomic, and genomic data to provide a comprehensive picture of drought and recovery responses in the leaves and roots of sugarcane. Transcriptomic profiling using oligoarrays and RNA-seq identified 2898 (out of 21,902) and 46,062 (out of 373,869) transcripts as differentially expressed, respectively. Co-expression analysis revealed modules enriched in photosynthesis, small molecule metabolism, alpha-amino acid metabolism, trehalose biosynthesis, serine family amino acid metabolism, and carbohydrate transport. Together, our findings reveal that carbohydrate metabolism is coordinated with the degradation of amino acids to provide carbon skeletons to the tricarboxylic acid cycle. This coordination may help to maintain energetic balance during drought stress adaptation, facilitating recovery after the stress is alleviated. Our results shed light on candidate regulatory elements and pave the way to biotechnology strategies towards the development of drought-tolerant sugarcane plants.

Assembly of the 373k gene space of the polyploid sugarcane genome reveals reservoirs of functional diversity in the world's leading biomass crop

BACKGROUND

Sugarcane cultivars are polyploid interspecific hybrids of giant genomes, typically with 10-13 sets of chromosomes from 2 Saccharum species. The ploidy, hybridity, and size of the genome, estimated to have >10 Gb, pose a challenge for sequencing.

RESULTS

Here we present a gene space assembly of SP80-3280, including 373,869 putative genes and their potential regulatory regions. The alignment of single-copy genes in diploid grasses to the putative genes indicates that we could resolve 2-6 (up to 15) putative homo(eo)logs that are 99.1% identical within their coding sequences. Dissimilarities increase in their regulatory regions, and gene promoter analysis shows differences in regulatory elements within gene families that are expressed in a species-specific manner. We exemplify these differences for sucrose synthase (SuSy) and phenylalanine ammonia-lyase (PAL), 2 gene families central to carbon partitioning. SP80-3280 has particular regulatory elements involved in sucrose synthesis not found in the ancestor Saccharum spontaneum. PAL regulatory elements are found in co-expressed genes related to fiber synthesis within gene networks defined during plant growth and maturation. Comparison with sorghum reveals predominantly bi-allelic variations in sugarcane, consistent with the formation of 2 "subgenomes" after their divergence ∼3.8-4.6 million years ago and reveals single-nucleotide variants that may underlie their differences.

CONCLUSIONS

This assembly represents a large step towards a whole-genome assembly of a commercial sugarcane cultivar. It includes a rich diversity of genes and homo(eo)logous resolution for a representative fraction of the gene space, relevant to improve biomass and food production.

Spontaneous mutation in the AgrRS two-component regulatory system ofCupriavidus metalliduransresults in enhanced silver resistance

Cupriavidus metalliduransis able to adapt to toxic silver concentrations and previously uncharacterized periplasmic proteins played a crucial role in this adaptation process.

Identification of Single Nucleotide Non-coding Driver Mutations in Cancer

Recent whole-genome sequencing studies have identified millions of somatic variants present in tumor samples. Most of these variants reside in non-coding regions of the genome potentially affecting transcriptional and post-transcriptional gene regulation. Although a few hallmark examples of driver mutations in non-coding regions have been reported, the functional role of the vast majority of somatic non-coding variants remains to be determined. This is because the few driver variants in each sample must be distinguished from the thousands of passenger variants and because the logic of regulatory element function has not yet been fully elucidated. Thus, variants prioritized based on mutational burden and location within regulatory elements need to be validated experimentally. This is generally achieved by combining assays that measure physical binding, such as chromatin immunoprecipitation, with those that determine regulatory activity, such as luciferase reporter assays. Here, we present an overview of in silico approaches used to prioritize somatic non-coding variants and the experimental methods used for functional validation and characterization.

Allelic Imbalance in Regulation of ANRIL through Chromatin Interaction at 9p21 Endometriosis Risk Locus

Genome-wide association studies (GWASs) have discovered numerous single nucleotide polymorphisms (SNPs) associated with human complex disorders. However, functional characterization of the disease-associated SNPs remains a formidable challenge. Here we explored regulatory mechanism of a SNP on chromosome 9p21 associated with endometriosis by leveraging “allele-specific” functional genomic approaches. By re-sequencing 1.29 Mb of 9p21 region and scrutinizing DNase-seq data from the ENCODE project, we prioritized rs17761446 as a candidate functional variant that was in perfect linkage disequilibrium with the original GWAS SNP (rs10965235) and located on DNase I hypersensitive site. Chromosome conformation capture followed by high-throughput sequencing revealed that the protective G allele of rs17761446 exerted stronger chromatin interaction with ANRIL promoter. We demonstrated that the protective allele exhibited preferential binding affinities to TCF7L2 and EP300 by bioinformatics and chromatin immunoprecipitation (ChIP) analyses. ChIP assays for histone H3 lysine 27 acetylation and RNA polymerase II reinforced the enhancer activity of the SNP site. The allele specific expression analysis for eutopic endometrial tissues and endometrial carcinoma cell lines showed that rs17761446 was a cis-regulatory variant where G allele was associated with increased ANRIL expression. Our work illuminates the allelic imbalances in a series of transcriptional regulation from factor binding to gene expression mediated by chromatin interaction underlie the molecular mechanism of 9p21 endometriosis risk locus. Functional genomics on common disease will unlock functional aspect of genotype-phenotype correlations in the post-GWAS stage.

A large number of variants associated with human complex diseases have been discovered by genome-wide association studies (GWASs). These discoveries have been anticipated to be translated into the definitive understanding of disease pathogeneses; however, functional characterization of the disease-associated SNPs remains a formidable challenge. Here we explored regulatory mechanism of a variant on chromosome 9p21 associated with endometriosis, a common gynecological disorder. By scrutinizing linkage disequilibrium structure and DNase I hypersensitive sites across the risk locus, we prioritized rs17761446 as a candidate causal variant. The results of our “allele-specific” functional genomic approaches sheds light on regulatory mechanisms underlying 9p21 endometriosis risk locus, in which preferential bindings of TCF7L2 and its coactivator EP300 to the protective G allele of rs17761446 lead to stronger chromatin interaction with the promoter of ANRIL, which in turn activate transcription of the non-coding RNA. Motivated by the fact that TCF7L2 was a key transcription factor of Wnt signaling pathway, we postulated that the induction of Wnt signaling activated expression levels of ANRIL and cell cycle inhibitors, CDKN2A/2B. Functional genomics on common disease will unlock functional aspect of genotype-phenotype correlations in the post-GWAS stage.

Total Binding Affinity Profiles of Regulatory Regions Predict Transcription Factor Binding and Gene Expression in Human Cells

Transcription factors regulate gene expression by binding regulatory DNA. Understanding the rules governing such binding is an essential step in describing the network of regulatory interactions, and its pathological alterations. We show that describing regulatory regions in terms of their profile of total binding affinities for transcription factors leads to increased predictive power compared to methods based on the identification of discrete binding sites. This applies both to the prediction of transcription factor binding as revealed by ChIP-seq experiments and to the prediction of gene expression through RNA-seq. Further significant improvements in predictive power are obtained when regulatory regions are defined based on chromatin states inferred from histone modification data.

Evolutionary Conservation and Diversification of Puf RNA Binding Proteins and Their mRNA Targets

Reprogramming of a gene’s expression pattern by acquisition and loss of sequences recognized by specific regulatory RNA binding proteins may be a major mechanism in the evolution of biological regulatory programs. We identified that RNA targets of Puf3 orthologs have been conserved over 100–500 million years of evolution in five eukaryotic lineages. Focusing on Puf proteins and their targets across 80 fungi, we constructed a parsimonious model for their evolutionary history. This model entails extensive and coordinated changes in the Puf targets as well as changes in the number of Puf genes and alterations of RNA binding specificity including that: 1) Binding of Puf3 to more than 200 RNAs whose protein products are predominantly involved in the production and organization of mitochondrial complexes predates the origin of budding yeasts and filamentous fungi and was maintained for 500 million years, throughout the evolution of budding yeast. 2) In filamentous fungi, remarkably, more than 150 of the ancestral Puf3 targets were gained by Puf4, with one lineage maintaining both Puf3 and Puf4 as regulators and a sister lineage losing Puf3 as a regulator of these RNAs. The decrease in gene expression of these mRNAs upon deletion of Puf4 in filamentous fungi (N. crassa) in contrast to the increase upon Puf3 deletion in budding yeast (S. cerevisiae) suggests that the output of the RNA regulatory network is different with Puf4 in filamentous fungi than with Puf3 in budding yeast. 3) The coregulated Puf4 target set in filamentous fungi expanded to include mitochondrial genes involved in the tricarboxylic acid (TCA) cycle and other nuclear-encoded RNAs with mitochondrial function not bound by Puf3 in budding yeast, observations that provide additional evidence for substantial rewiring of post-transcriptional regulation. 4) Puf3 also expanded and diversified its targets in filamentous fungi, gaining interactions with the mRNAs encoding the mitochondrial electron transport chain (ETC) complex I as well as hundreds of other mRNAs with nonmitochondrial functions. The many concerted and conserved changes in the RNA targets of Puf proteins strongly support an extensive role of RNA binding proteins in coordinating gene expression, as originally proposed by Keene. Rewiring of Puf-coordinated mRNA targets and transcriptional control of the same genes occurred at different points in evolution, suggesting that there have been distinct adaptations via RNA binding proteins and transcription factors. The changes in Puf targets and in the Puf proteins indicate an integral involvement of RNA binding proteins and their RNA targets in the adaptation, reprogramming, and function of gene expression.

A map of the evolutionary history of Puf proteins and their RNA targets shows that reprogramming of global gene expression programs via adaptive mutations that affect protein-RNA interactions is an important source of biological diversity.

We set out to trace the evolutionary history of an RNA binding protein and how its interactions with targets change over evolution. Identifying this natural history is a step toward understanding the critical differences between organisms and how gene expression programs are rewired during evolution. Using bioinformatics and experimental approaches, we broadly surveyed the evolution of binding targets of a particular family of RNA binding proteins—the Puf proteins, whose protein sequences and target RNA sequences are relatively well-characterized—across 99 eukaryotic species. We found five groups of species in which targets have been conserved for at least 100 million years and then took advantage of genome sequences from a large number of fungal species to deeply investigate the conservation and changes in Puf proteins and their RNA targets. Our analyses identified multiple and extensive reconfigurations during the natural history of fungi and suggest that RNA binding proteins and their RNA targets are profoundly involved in evolutionary reprogramming of gene expression and help define distinct programs unique to each organism. Continuing to uncover the natural history of RNA binding proteins and their interactions will provide a unique window into the gene expression programs of present day species and point to new ways to engineer gene expression programs.

Directed targeting of chromatin to the nuclear lamina is mediated by chromatin state and A-type lamins

Nuclear organization has been implicated in regulating gene activity. Recently, large developmentally regulated regions of the genome dynamically associated with the nuclear lamina have been identified. However, little is known about how these lamina-associated domains (LADs) are directed to the nuclear lamina. We use our tagged chromosomal insertion site system to identify small sequences from borders of fibroblast-specific variable LADs that are sufficient to target these ectopic sites to the nuclear periphery. We identify YY1 (Ying-Yang1) binding sites as enriched in relocating sequences. Knockdown of YY1 or lamin A/C, but not lamin A, led to a loss of lamina association. In addition, targeted recruitment of YY1 proteins facilitated ectopic LAD formation dependent on histone H3 lysine 27 trimethylation and histone H3 lysine di- and trimethylation. Our results also reveal that endogenous loci appear to be dependent on lamin A/C, YY1, H3K27me3, and H3K9me2/3 for maintenance of lamina-proximal positioning.

Reversal of senescence by N resupply to N-starved Arabidopsis thaliana: transcriptomic and metabolomic consequences

Leaf senescence is a developmentally controlled process, which is additionally modulated by a number of adverse environmental conditions. Nitrogen shortage is a well-known trigger of precocious senescence in many plant species including crops, generally limiting biomass and seed yield. However, leaf senescence induced by nitrogen starvation may be reversed when nitrogen is resupplied at the onset of senescence. Here, the transcriptomic, hormonal, and global metabolic rearrangements occurring during nitrogen resupply-induced reversal of senescence in Arabidopsis thaliana were analysed. The changes induced by senescence were essentially in keeping with those previously described; however, these could, by and large, be reversed. The data thus indicate that plants undergoing senescence retain the capacity to sense and respond to the availability of nitrogen nutrition. The combined data are discussed in the context of the reversibility of the senescence programme and the evolutionary benefit afforded thereby. Future prospects for understanding and manipulating this process in both Arabidopsis and crop plants are postulated.

Scan for Motifs: a webserver for the analysis of post-transcriptional regulatory elements in the 3' untranslated regions (3' UTRs) of mRNAs

Background

Gene expression in vertebrate cells may be controlled post-transcriptionally through regulatory elements in mRNAs. These are usually located in the untranslated regions (UTRs) of mRNA sequences, particularly the 3′UTRs.

Results

Scan for Motifs (SFM) simplifies the process of identifying a wide range of regulatory elements on alignments of vertebrate 3′UTRs. SFM includes identification of both RNA Binding Protein (RBP) sites and targets of miRNAs. In addition to searching pre-computed alignments, the tool provides users the flexibility to search their own sequences or alignments. The regulatory elements may be filtered by expected value cutoffs and are cross-referenced back to their respective sources and literature. The output is an interactive graphical representation, highlighting potential regulatory elements and overlaps between them. The output also provides simple statistics and links to related resources for complementary analyses. The overall process is intuitive and fast. As SFM is a free web-application, the user does not need to install any software or databases.

Conclusions

Visualisation of the binding sites of different classes of effectors that bind to 3′UTRs will facilitate the study of regulatory elements in 3′ UTRs.