Demethylation of the Protein Phosphatase PP2A Promotes Demethylation of Histones to Enable Their Function as a Methyl Group Sink.

Dysregulation of chromatin methylation is associated with defects in cellular differentiation as well as a variety of cancers. How cells regulate the opposing activities of histone methyltransferase and demethylase enzymes to set the methylation status of the epigenome for proper control of gene expression and metabolism remains poorly understood. Here, we show that loss of methylation of the major phosphatase PP2A in response to methionine starvation activates the demethylation of histones through hyperphosphorylation of specific demethylase enzymes. In parallel, this regulatory mechanism enables cells to preserve SAM by increasing SAH to limit SAM consumption by methyltransferase enzymes. Mutants lacking the PP2A methyltransferase or the effector H3K36 demethylase Rph1 exhibit elevated SAM levels and are dependent on cysteine due to reduced capacity to sink the methyl groups of SAM. Therefore, PP2A directs the methylation status of histones by regulating the phosphorylation status of histone demethylase enzymes in response to SAM levels.

A Metabolic Function for Phospholipid and Histone Methylation.

S-adenosylmethionine (SAM) is the methyl donor for biological methylation modifications that regulate protein and nucleic acid functions. Here, we show that methylation of a phospholipid, phosphatidylethanolamine (PE), is a major consumer of SAM. The induction of phospholipid biosynthetic genes is accompanied by induction of the enzyme that hydrolyzes S-adenosylhomocysteine (SAH), a product and inhibitor of methyltransferases. Beyond its function for the synthesis of phosphatidylcholine (PC), the methylation of PE facilitates the turnover of SAM for the synthesis of cysteine and glutathione through transsulfuration. Strikingly, cells that lack PE methylation accumulate SAM, which leads to hypermethylation of histones and the major phosphatase PP2A, dependency on cysteine, and sensitivity to oxidative stress. Without PE methylation, particular sites on histones then become methyl sinks to enable the conversion of SAM to SAH. These findings reveal an unforeseen metabolic function for phospholipid and histone methylation intrinsic to the life of a cell.

PmiREN2.0: from data annotation to functional exploration of plant microRNAs.

Nearly 200 plant genomes have been sequenced over the last two years, and new functions of plant microRNAs (miRNAs) have been revealed. Therefore, timely update of the plant miRNA databases by incorporating miRNAs from the newly sequenced species and functional information is required to provide useful resources for advancing plant miRNA research. Here we report the update of PmiREN2.0 (https://pmiren.com/) with an addition of 19 363 miRNA entries from 91 plants, doubling the amount of data in the original version. Meanwhile, abundant regulatory information centred on miRNAs was added, including predicted upstream transcription factors through binding motifs scanning and elaborate annotation of miRNA targets. As an example, a genome-wide regulatory network centred on miRNAs was constructed for Arabidopsis. Furthermore, phylogenetic trees of conserved miRNA families were built to expand the understanding of miRNA evolution across the plant lineages. These data are helpful to deduce the regulatory relationships concerning miRNA functions in diverse plants. Beside the new data, a suite of design tools was incorporated to facilitate experimental practice. Finally, a forum named 'PmiREN Community' was added for discussion and resource and new discovery sharing. With these upgrades, PmiREN2.0 should serve the community better and accelerate miRNA research in plants.

POU2AF3-rearranged sarcomas: A novel tumor defined by fusions of EWSR1 or FUS to a gene formerly designated COLCA2.

Gene fusions involving EWSR1 or FUS as the 5' partner have been reported in a diverse array of sarcomas. Here, we characterize the histopathology and genomics of six tumors harboring a gene fusion between EWSR1 or FUS and POU2AF3, an understudied, putative colorectal cancer predisposition gene. Striking morphologic features reminiscent of synovial sarcoma were observed including a biphasic appearance with variable fusiform to epithelioid cytomorphology and staghorn-type vasculature. RNA sequencing demonstrated variable breakpoints in EWSR1/FUS along with similar breakpoints in POU2AF3 that encompassed a 3' portion of this gene. For cases in which additional information was available, the behavior of these neoplasms was aggressive with local spread and/or distant metastases. Although further studies are needed to confirm the functional significance of our findings, POU2AF3 fusions to EWSR1 or FUS may define a novel type of POU2AF3-rearranged sarcomas with aggressive, malignant behavior.

Miniature Inverted-repeat Transposable Elements Drive Rapid MicroRNA Diversification in Angiosperms.

MicroRNAs (miRNAs) are fast evolving endogenous small RNAs that regulate organism function and behavior in both animals and plants. Although models for de novo miRNA biogenesis have been proposed, the genomic mechanisms driving swift diversification of the miRNA repertoires in plants remain elusive. Here, by comprehensively analyzing 21 phylogenetically representative plant species, ranging from green algae to angiosperms, we systematically identified de novo miRNA events associated with 8,649 miRNA loci. We found that 399 (4.6%), 466 (5.4%), and 1,402 (16.2%) miRNAs were derived from inverted gene duplication events, long terminal repeats of retrotransposons, and miniature inverted-repeat transposable elements (MITEs), respectively. Among the miRNAs of these origins, MITEs, especially those belonging to the Mutator, Tc1/Mariner, and PIF/Harbinger superfamilies, were the predominant genomic source for de novo miRNAs in the 15 examined angiosperms but not in the six non-angiosperms. Our data further illustrated a transposition-transcription process by which MITEs are converted into new miRNAs (termed MITE-miRNAs) whereby properly sized MITEs are transcribed and therefore become potential substrates for the miRNA processing machinery by transposing into introns of active genes. By analyzing the 58,038 putative target genes for the 8,095 miRNAs, we found that the target genes of MITE-miRNAs were preferentially associated with response to environmental stimuli such as temperature, suggesting that MITE-miRNAs are pertinent to plant adaptation. Collectively, these findings demonstrate that molecular conversion of MITEs is a genomic mechanism leading to rapid and continuous changes to the miRNA repertoires in angiosperm.

MicroRNA annotation in plants: current status and challenges.

Last two decades, the studies on microRNAs (miRNAs) and the numbers of annotated miRNAs in plants and animals have surged. Herein, we reviewed the current progress and challenges of miRNA annotation in plants. Via the comparison of plant and animal miRNAs, we pinpointed out the difficulties on plant miRNA annotation and proposed potential solutions. In terms of recalling the history of methods and criteria in plant miRNA annotation, we detailed how the major progresses made and evolved. By collecting and categorizing bioinformatics tools for plant miRNA annotation, we surveyed their advantages and disadvantages, especially for ones with the principle of mimicking the miRNA biogenesis pathway by parsing deeply sequenced small RNA (sRNA) libraries. In addition, we summarized all available databases hosting plant miRNAs, and posted the potential optimization solutions such as how to increase the signal-to-noise ratio (SNR) in these databases. Finally, we discussed the challenges and perspectives of plant miRNA annotations, and indicated the possibilities offered by an all-in-one tool and platform according to the integration of artificial intelligence.

Optimal audio beam pattern synthesis for an enhanced parametric array loudspeaker.

A parametric array loudspeaker (PAL) generates highly directional audible sound in air with a small aperture size compared to a conventional loudspeaker. But in indoor applications, the long propagation distance of a PAL causes reflections, which disturbs the reproduction of narrow audio beams. Moreover, sound distortion appears along the off-axis direction due to the frequency dependence of the beam width. This study proposed an optimal audio beam pattern synthesis for a PAL-based convex optimization, which can design the audio beam of a PAL with an optimal solution. The proposed method overcame the mentioned limitations by applying it to a length-limited PAL for audio spot control and a multichannel PAL array for a constant beam width audio beam. In a length-limited PAL, the proposed method restricts the audio spot to a smaller region and weakens the sound leakage along the off-axis direction. Whereas in a multichannel PAL array, the proposed method also achieves a constant beam width near the radiator axis. Simulations and experiments verify the effectiveness of the proposed method, which will enhance the performance of a PAL in scenarios where control of the audio beam is required.

PmiREN: a comprehensive encyclopedia of plant miRNAs.

MicroRNAs (miRNAs) are small non-coding RNA molecules that function as diverse endogenous gene regulators at the post-transcriptional level. In the past two decades, as research effort on miRNA identification, function and evolution has soared, so has the demand for miRNA databases. However, the current plant miRNA databases suffer from several typical drawbacks, including a lack of entries for many important species, uneven annotation standards across different species, abundant questionable entries, and limited annotation. To address these issues, we developed a knowledge-based database called Plant miRNA Encyclopedia (PmiREN, http://www.pmiren.com/), which was based on uniform processing of sequenced small RNA libraries using miRDeep-P2, followed by manual curation using newly updated plant miRNA identification criteria, and comprehensive annotation. PmiREN currently contains 16,422 high confidence novel miRNA loci in 88 plant species and 3,966 retrieved from miRBase. For every miRNA entry, information on precursor sequence, precursor secondary structure, expression pattern, clusters and synteny in the genome, potential targets supported by Parallel Analysis of RNA Ends (PARE) sequencing, and references is attached whenever possible. PmiREN is hierarchically accessible and has eight built-in search engines. We believe PmiREN is useful for plant miRNA cataloguing and data mining, therefore a resource for data-driven miRNA research in plants.

Epithelial retinoic acid receptor ß regulates serum amyloid A expression and vitamin A-dependent intestinal immunity.

Vitamin A is a dietary component that is essential for the development of intestinal immunity. Vitamin A is absorbed and converted to its bioactive derivatives retinol and retinoic acid by the intestinal epithelium, yet little is known about how epithelial cells regulate vitamin A-dependent intestinal immunity. Here we show that epithelial cell expression of the transcription factor retinoic acid receptor ß (RARß) is essential for vitamin A-dependent intestinal immunity. Epithelial RARß activated vitamin A-dependent expression of serum amyloid A (SAA) proteins by binding directly to Saa promoters. In accordance with the known role of SAAs in regulating Th17 cell effector function, epithelial RARß promoted IL-17 production by intestinal Th17 cells. More broadly, epithelial RARß was required for the development of key vitamin A-dependent adaptive immune responses, including CD4+ T-cell homing to the intestine and the development of IgA-producing intestinal B cells. Our findings provide insight into how the intestinal epithelium senses dietary vitamin A status to regulate adaptive immunity, and highlight the role of epithelial cells in regulating intestinal immunity in response to diet.

Identification of Species-Specific MicroRNAs Provides Insights into Dynamic Evolution of MicroRNAs in Plants.

MicroRNAs (miRNAs) are an important class of regulatory small RNAs that program gene expression, mainly at the post-transcriptional level. Although sporadic examples of species-specific miRNAs (termed SS-miRNAs) have been reported, a genome-scale study across a variety of distant species has not been assessed. Here, by comprehensively analyzing miRNAs in 81 plant species phylogenetically ranging from chlorophytes to angiosperms, we identified 8048 species-specific miRNAs from 5499 families, representing over 61.2% of the miRNA families in the examined species. An analysis of the conservation from different taxonomic levels supported the high turnover rate of SS-miRNAs, even over short evolutionary distances. A comparison of the intrinsic features between SS-miRNAs and NSS-miRNAs (non-species-specific miRNAs) indicated that the AU content of mature miRNAs was the most striking difference. Our data further illustrated a significant bias of the genomic coordinates towards SS-miRNAs lying close to or within genes. By analyzing the 125,267 putative target genes for the 7966 miRNAs, we found the preferentially regulated functions of SS-miRNAs related to diverse metabolic processes. Collectively, these findings underscore the dynamic evolution of miRNAs in the species-specific lineages.

A toolbox of immunoprecipitation-grade monoclonal antibodies to human transcription factors.

A key component of efforts to address the reproducibility crisis in biomedical research is the development of rigorously validated and renewable protein-affinity reagents. As part of the US National Institutes of Health (NIH) Protein Capture Reagents Program (PCRP), we have generated a collection of 1,406 highly validated immunoprecipitation- and/or immunoblotting-grade mouse monoclonal antibodies (mAbs) to 737 human transcription factors, using an integrated production and validation pipeline. We used HuProt human protein microarrays as a primary validation tool to identify mAbs with high specificity for their cognate targets. We further validated PCRP mAbs by means of multiple experimental applications, including immunoprecipitation, immunoblotting, chromatin immunoprecipitation followed by sequencing (ChIP-seq), and immunohistochemistry. We also conducted a meta-analysis that identified critical variables that contribute to the generation of high-quality mAbs. All validation data, protocols, and links to PCRP mAb suppliers are available at http://proteincapture.org.

The dynamic landscape of fission yeast meiosis alternative-splice isoforms.

Alternative splicing increases the diversity of transcriptomes and proteomes in metazoans. The extent to which alternative splicing is active and functional in unicellular organisms is less understood. Here, we exploit a single-molecule long-read sequencing technique and develop an open-source software program called SpliceHunter to characterize the transcriptome in the meiosis of fission yeast. We reveal 14,353 alternative splicing events in 17,669 novel isoforms at different stages of meiosis, including antisense and read-through transcripts. Intron retention is the major type of alternative splicing, followed by alternate "intron in exon." Seven hundred seventy novel transcription units are detected; 53 of the predicted proteins show homology in other species and form theoretical stable structures. We report the complexity of alternative splicing along isoforms, including 683 intra-molecularly co-associated intron pairs. We compare the dynamics of novel isoforms based on the number of supporting full-length reads with those of annotated isoforms and explore the translational capacity and quality of novel isoforms. The evaluation of these factors indicates that the majority of novel isoforms are unlikely to be both condition-specific and translatable but consistent with the possibility of biologically functional novel isoforms. Moreover, the co-option of these unusual transcripts into newly born genes seems likely. Together, the results of this study highlight the diversity and dynamics at the isoform level in the sexual development of fission yeast.

miRDeep-P2: accurate and fast analysis of the microRNA transcriptome in plants.

MOTIVATION: Two major challenges arise when employing next-generation sequencing methods to comprehensively identify microRNAs (miRNAs) in plants: (i) how to minimize the false-positive inheritable to computational predictions and (ii) how to minimize the computational time required for analyzing the miRNA transcriptome in plants with complex and large genomes. RESULTS: We updated miRDeep-P to miRDeep-P2 (miRDP2) by employing a new filtering strategy and overhauling the algorithm. miRDP2 has been tested against miRNA transcriptomes in plants with increasing genome sizes that included Arabidopsis, rice, tomato, maize and wheat. Compared with miRDeep-P and several other computational tools, miRDP2 processes next-generation sequencing data with superior speed. By incorporating newly updated plant miRNA annotation criteria and developing a new scoring system, the accuracy of miRDP2 outperformed other programs. Taken together, our results demonstrate miRDP2 as a fast and accurate tool for analyzing the miRNA transcriptome in plants. AVAILABILITY AND IMPLEMENTATION: The miRDP2 are freely available from https://sourceforge.net/projects/mirdp2/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Selection of suitable reference genes for qRT-PCR expression analysis of Codonopsis pilosula under different experimental conditions.

Codonopsis pilosula is a well-known medicinal plant. Although its transcriptome sequence has been published, suitable reference genes have not been systematically identified for conducting expression analyses via quantitative real-time polymerase chain reaction (qRT-PCR). To screen appropriate genes for use with this species, we applied four different methods-GeNorm, NormFinder, BestKeeper, and RefFinder-to evaluate the stability of 13 candidates: CpiEF1Bb, CpiCACS, CpiF-Box, Cpiß-Tubulin, CpiGAPDH, CpiActin2, CpiAPT1, CpiActin7, CpiActin8, CpiRPL6, CpiHAF1, CpiTubulin6, and CpiUBQ12. Expression was examined by qRT-PCR for various tissue types, chemical treatments, and developmental stages. For all tested samples, CpiGAPDH proved to be the most stable. Comprehensive analysis indicated that the most stable internal reference genes were CpiF-Box and CpiCACS in different tissues and at different developmental stages, respectively. Under NaCl stress, CpiAPT1 was the best internal reference gene. For methyl jasmonate and abscisic acid treatments, CpiGAPDH and CpiF-Box, respectively, presented the highest degree of expression stability. Based on these findings, we chose CpiSPL9 as the target gene for validating the suitability of these selected reference genes. All of these results provide a foundation for accurate quantification of expression levels by genes of interest in C. pilosula.

Dynamic motif occupancy (DynaMO) analysis identifies transcription factors and their binding sites driving dynamic biological processes.

Biological processes are usually associated with genome-wide remodeling of transcription driven by transcription factors (TFs). Identifying key TFs and their spatiotemporal binding patterns are indispensable to understanding how dynamic processes are programmed. However, most methods are designed to predict TF binding sites only. We present a computational method, dynamic motif occupancy analysis (DynaMO), to infer important TFs and their spatiotemporal binding activities in dynamic biological processes using chromatin profiling data from multiple biological conditions such as time-course histone modification ChIP-seq data. In the first step, DynaMO predicts TF binding sites with a random forests approach. Next and uniquely, DynaMO infers dynamic TF binding activities at predicted binding sites using their local chromatin profiles from multiple biological conditions. Another landmark of DynaMO is to identify key TFs in a dynamic process using a clustering and enrichment analysis of dynamic TF binding patterns. Application of DynaMO to the yeast ultradian cycle, mouse circadian clock and human neural differentiation exhibits its accuracy and versatility. We anticipate DynaMO will be generally useful for elucidating transcriptional programs in dynamic processes.

Stress response factors drive regrowth of quiescent cells.

Quiescent cells exploit an array of transcription factors to activate stress response machinery and maintain survival under nutrient-limited conditions. Our recent findings reveal that these transcription factors also play an important role in the exit of quiescence and regrowth. By studying Saccharomyces cerevisiae under a continuous, nutrient-limited condition, we found that Msn2 and Msn4 function as master regulators of glycolytic genes in the quiescent-like phase. They control the timing of transition from quiescence to growth by regulating the accumulation rate of acetyl-CoA, a key metabolite that is downstream of glycolysis and drives growth. These findings suggest a model that Msn2/4 not only protect the cells from starvation but also facilitate their regrowth from quiescence. Thus, understanding the functions of stress response transcription factors in metabolic regulation will provide deeper insight into how quiescent cells manage the capacity of regrowth.

Functional dissection of lysine deacetylases reveals that HDAC1 and p300 regulate AMPK.

First identified as histone-modifying proteins, lysine acetyltransferases (KATs) and deacetylases (KDACs) antagonize each other through modification of the side chains of lysine residues in histone proteins. Acetylation of many non-histone proteins involved in chromatin, metabolism or cytoskeleton regulation were further identified in eukaryotic organisms, but the corresponding enzymes and substrate-specific functions of the modifications are unclear. Moreover, mechanisms underlying functional specificity of individual KDACs remain enigmatic, and the substrate spectra of each KDAC lack comprehensive definition. Here we dissect the functional specificity of 12 critical human KDACs using a genome-wide synthetic lethality screen in cultured human cells. The genetic interaction profiles revealed enzyme-substrate relationships between individual KDACs and many important substrates governing a wide array of biological processes including metabolism, development and cell cycle progression. We further confirmed that acetylation and deacetylation of the catalytic subunit of the adenosine monophosphate-activated protein kinase (AMPK), a critical cellular energy-sensing protein kinase complex, is controlled by the opposing catalytic activities of HDAC1 and p300. Deacetylation of AMPK enhances physical interaction with the upstream kinase LKB1, leading to AMPK phosphorylation and activation, and resulting in lipid breakdown in human liver cells. These findings provide new insights into previously underappreciated metabolic regulatory roles of HDAC1 in coordinating nutrient availability and cellular responses upstream of AMPK, and demonstrate the importance of high-throughput genetic interaction profiling to elucidate functional specificity and critical substrates of individual human KDACs potentially valuable for therapeutic applications.

Overexpression of microRNA408 enhances photosynthesis, growth, and seed yield in diverse plants.

The ability of a plant to produce grain, fruit, or forage depends ultimately on photosynthesis. There have been few attempts, however, to study microRNAs, which are a class of endogenous small RNAs post-transcriptionally programming gene expression, in relation to photosynthetic traits. We focused on miR408, one of the most conserved plant miRNAs, and overexpressed it in parallel in Arabidopsis, tobacco, and rice. The transgenic plants all exhibited increased copper content in the chloroplast, elevated abundance of plastocyanin, and an induction of photosynthetic genes. By means of gas exchange and optical spectroscopy analyses, we showed that higher expression of miR408 leads to enhanced photosynthesis through improving efficiency of irradiation utilization and the capacity for carbon dioxide fixation. Consequently, miR408 hyper-accumulating plants exhibited higher rate of vegetative growth. An enlargement of seed size was also observed in all three species overproducing miR408. Moreover, we conducted a 2-year-two-location field trial and observed miR408 overexpression in rice significantly increased yield, which was primarily attributed to an elevation in grain weight. Taken together, these results demonstrate that miR408 is a positive regulator of photosynthesis and that its genetic engineering is a promising route for enhancing photosynthetic performance and yield in diverse plants.

Imaging-based amplitude laser beam shaping for material processing by 2D reflectivity tuning of a spatial light modulator.

We have demonstrated an imaging-based amplitude laser-beam-shaping technique for material processing by 2D reflectivity tuning of a spatial light modulator. Intensity masks with 256 gray levels were designed to shape the input laser beam in the outline profile and inside intensity distribution. Squared and circular flattop beam shapes were obtained at the diffractive near-field and then reconstructed at an image plane of an f-theta lens (f∼100 mm). The observed intensity distribution inside the beam-shaping geometry was much more even than using binary masks. The ablation footprint well matches the desired beam shape.

Three-Dimensional Biologically Relevant Spectrum (BRS-3D): Shape Similarity Profile Based on PDB Ligands as Molecular Descriptors.

The crystallized ligands in the Protein Data Bank (PDB) can be treated as the inverse shapes of the active sites of corresponding proteins. Therefore, the shape similarity between a molecule and PDB ligands indicated the possibility of the molecule to bind with the targets. In this paper, we proposed a shape similarity profile that can be used as a molecular descriptor for ligand-based virtual screening. First, through three-dimensional (3D) structural clustering, 300 diverse ligands were extracted from the druggable protein-ligand database, sc-PDB. Then, each of the molecules under scrutiny was flexibly superimposed onto the 300 ligands. Superimpositions were scored by shape overlap and property similarity, producing a 300 dimensional similarity array termed the "Three-Dimensional Biologically Relevant Spectrum (BRS-3D)". Finally, quantitative or discriminant models were developed with the 300 dimensional descriptor using machine learning methods (support vector machine). The effectiveness of this approach was evaluated using 42 benchmark data sets from the G protein-coupled receptor (GPCR) ligand library and the GPCR decoy database (GLL/GDD). We compared the performance of BRS-3D with other 2D and 3D state-of-the-art molecular descriptors. The results showed that models built with BRS-3D performed best for most GLL/GDD data sets. We also applied BRS-3D in histone deacetylase 1 inhibitors screening and GPCR subtype selectivity prediction. The advantages and disadvantages of this approach are discussed.