MatK impacts differential chloroplast translation by limiting spliced tRNA-K(UUU) abundance.

The protein levels of chloroplast photosynthetic genes and genes related to the chloroplast genetic apparatus vary to adapt to different conditions. However, the underlying mechanisms governing these variations remain unclear. The chloroplast intron Maturase K is encoded within the trnK intron and has been suggested to be required for splicing several group IIA introns, including the trnK intron. In this study, we used RNA immunoprecipitation followed by high-throughput sequencing (RIP-Seq) to identify MatK's preference for binding to group IIA intron domains I and VI within target transcripts. Importantly, these domains are crucial for splice site selection, and we discovered alternative 5'-splice sites in three MatK target introns. The resulting alternative trnK lariat structure showed increased accumulation during heat acclimation. The cognate codon of tRNA-K(UUU) is highly enriched in mRNAs encoding ribosomal proteins and a trnK-matK over-expressor exhibited elevated levels of the spliced tRNA-K(UUU). Ribosome profiling analysis of the overexpressor revealed a significant up-shift in the translation of ribosomal proteins compared to photosynthetic genes. Our findings suggest the existence of a novel regulatory mechanism linked to the abundance of tRNA-K(UUU), enabling the differential expression of functional chloroplast gene groups.

Liam tackles complex multimodal single-cell data integration challenges.

Multi-omics characterization of single cells holds outstanding potential for profiling the dynamics and relations of gene regulatory states of thousands of cells. How to integrate multimodal data is an open problem, especially when aiming to combine data from multiple sources or conditions containing both biological and technical variation. We introduce liam, a flexible model for the simultaneous horizontal and vertical integration of paired single-cell multimodal data and mosaic integration of paired with unimodal data. Liam learns a joint low-dimensional representation of the measured modalities, which proves beneficial when the information content or quality of the modalities differ. Its integration accounts for complex batch effects using a tunable combination of conditional and adversarial training, which can be optimized using replicate information while retaining selected biological variation. We demonstrate liam's superior performance on multiple paired multimodal data types, including Multiome and CITE-seq data, and in mosaic integration scenarios. Our detailed benchmarking experiments illustrate the complexities and challenges remaining for integration and the meaningful assessment of its success.

Genome-wide screening identifies Trim33 as an essential regulator of dendritic cell differentiation.

The development of dendritic cells (DCs), including antigen-presenting conventional DCs (cDCs) and cytokine-producing plasmacytoid DCs (pDCs), is controlled by the growth factor Flt3 ligand (Flt3L) and its receptor Flt3. We genetically dissected Flt3L-driven DC differentiation using CRISPR-Cas9-based screening. Genome-wide screening identified multiple regulators of DC differentiation including subunits of TSC and GATOR1 complexes, which restricted progenitor growth but enabled DC differentiation by inhibiting mTOR signaling. An orthogonal screen identified the transcriptional repressor Trim33 (TIF-1γ) as a regulator of DC differentiation. Conditional targeting in vivo revealed an essential role of Trim33 in the development of all DCs, but not of monocytes or granulocytes. In particular, deletion of Trim33 caused rapid loss of DC progenitors, pDCs, and the cross-presenting cDC1 subset. Trim33-deficient Flt3+ progenitors up-regulated pro-inflammatory and macrophage-specific genes but failed to induce the DC differentiation program. Collectively, these data elucidate mechanisms that control Flt3L-driven differentiation of the entire DC lineage and identify Trim33 as its essential regulator.

Stress-induced nuclear speckle reorganization is linked to activation of immediate early gene splicing.

Current models posit that nuclear speckles (NSs) serve as reservoirs of splicing factors and facilitate posttranscriptional mRNA processing. Here, we discovered that ribotoxic stress induces a profound reorganization of NSs with enhanced recruitment of factors required for splice-site recognition, including the RNA-binding protein TIAR, U1 snRNP proteins and U2-associated factor 65, as well as serine 2 phosphorylated RNA polymerase II. NS reorganization relies on the stress-activated p38 mitogen-activated protein kinase (MAPK) pathway and coincides with splicing activation of both pre-existing and newly synthesized pre-mRNAs. In particular, ribotoxic stress causes targeted excision of retained introns from pre-mRNAs of immediate early genes (IEGs), whose transcription is induced during the stress response. Importantly, enhanced splicing of the IEGs ZFP36 and FOS is accompanied by relocalization of the corresponding nuclear mRNA foci to NSs. Our study reveals NSs as a dynamic compartment that is remodeled under stress conditions, whereby NSs appear to become sites of IEG transcription and efficient cotranscriptional splicing.

Toward Identification of Functional Sequences and Variants in Noncoding DNA.

Understanding the noncoding part of the genome, which encodes gene regulation, is necessary to identify genetic mechanisms of disease and translate findings from genome-wide association studies into actionable results for treatments and personalized care. Here we provide an overview of the computational analysis of noncoding regions, starting from gene-regulatory mechanisms and their representation in data. Deep learning methods, when applied to these data, highlight important regulatory sequence elements and predict the functional effects of genetic variants. These and other algorithms are used to predict damaging sequence variants. Finally, we introduce rare-variant association tests that incorporate functional annotations and predictions in order to increase interpretability and statistical power.

Brassinosteroid gene regulatory networks at cellular resolution in the Arabidopsis root.

Brassinosteroids are plant steroid hormones that regulate diverse processes, such as cell division and cell elongation, through gene regulatory networks that vary in space and time. By using time series single-cell RNA sequencing to profile brassinosteroid-responsive gene expression specific to different cell types and developmental stages of the Arabidopsis root, we identified the elongating cortex as a site where brassinosteroids trigger a shift from proliferation to elongation associated with increased expression of cell wall-related genes. Our analysis revealed HOMEOBOX FROM ARABIDOPSIS THALIANA 7 (HAT7) and GT-2-LIKE 1 (GTL1) as brassinosteroid-responsive transcription factors that regulate cortex cell elongation. These results establish the cortex as a site of brassinosteroid-mediated growth and unveil a brassinosteroid signaling network regulating the transition from proliferation to elongation, which illuminates aspects of spatiotemporal hormone responses.

A computational map of the human-SARS-CoV-2 protein-RNA interactome predicted at single-nucleotide resolution.

RNA-binding proteins (RBPs) are critical host factors for viral infection, however, large scale experimental investigation of the binding landscape of human RBPs to viral RNAs is costly and further complicated due to sequence variation between viral strains. To fill this gap, we investigated the role of RBPs in the context of SARS-CoV-2 by constructing the first in silico map of human RBP-viral RNA interactions at nucleotide-resolution using two deep learning methods (pysster and DeepRiPe) trained on data from CLIP-seq experiments on more than 100 human RBPs. We evaluated conservation of RBP binding between six other human pathogenic coronaviruses and identified sites of conserved and differential binding in the UTRs of SARS-CoV-1, SARS-CoV-2 and MERS. We scored the impact of mutations from 11 variants of concern on protein-RNA interaction, identifying a set of gain- and loss-of-binding events, as well as predicted the regulatory impact of putative future mutations. Lastly, we linked RBPs to functional, OMICs and COVID-19 patient data from other studies, and identified MBNL1, FTO and FXR2 RBPs as potential clinical biomarkers. Our results contribute towards a deeper understanding of how viruses hijack host cellular pathways and open new avenues for therapeutic intervention.

Altered and allele-specific open chromatin landscape reveals epigenetic and genetic regulators of innate immunity in COVID-19.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection causes severe COVID-19 in some patients and mild COVID-19 in others. Dysfunctional innate immune responses have been identified to contribute to COVID-19 severity, but the key regulators are still unknown. Here, we present an integrative single-cell multi-omics analysis of peripheral blood mononuclear cells from hospitalized and convalescent COVID-19 patients. In classical monocytes, we identified genes that were potentially regulated by differential chromatin accessibility. Then, sub-clustering and motif-enrichment analyses revealed disease condition-specific regulation by transcription factors and their targets, including an interaction between C/EBPs and a long-noncoding RNA LUCAT1, which we validated through loss-of-function experiments. Finally, we investigated genetic risk variants that exhibit allele-specific open chromatin (ASoC) in COVID-19 patients and identified a SNP rs6800484-C, which is associated with lower expression of CCR2 and may contribute to higher viral loads and higher risk of COVID-19 hospitalization. Altogether, our study highlights the diverse genetic and epigenetic regulators that contribute to COVID-19.

A critical period of translational control during brain development at codon resolution.

Translation modulates the timing and amplification of gene expression after transcription. Brain development requires uniquely complex gene expression patterns, but large-scale measurements of translation directly in the prenatal brain are lacking. We measure the reactants, synthesis and products of mRNA translation spanning mouse neocortex neurogenesis, and discover a transient window of dynamic regulation at mid-gestation. Timed translation upregulation of chromatin-binding proteins like Satb2, which is essential for neuronal subtype differentiation, restricts protein expression in neuronal lineages despite broad transcriptional priming in progenitors. In contrast, translation downregulation of ribosomal proteins sharply decreases ribosome biogenesis, coinciding with a major shift in protein synthesis dynamics at mid-gestation. Changing activity of eIF4EBP1, a direct inhibitor of ribosome biogenesis, is concurrent with ribosome downregulation and affects neurogenesis of the Satb2 lineage. Thus, the molecular logic of brain development includes the refinement of transcriptional programs by translation. Modeling of the developmental neocortex translatome is provided as an open-source searchable resource at https://shiny.mdc-berlin.de/cortexomics .

Protocol for fast scRNA-seq raw data processing using scKB and non-arbitrary quality control with COPILOT.

We describe a protocol to perform fast and non-arbitrary quality control of single-cell RNA sequencing (scRNA-seq) raw data using scKB and COPILOT. scKB is a wrapper script of kallisto and bustools for accelerated alignment and transcript count matrix generation, which runs significantly faster than the popular tool Cell Ranger. COPILOT then offers non-arbitrary background noise removal by comparing distributions of low-quality and high-quality cells. Together, this protocol streamlines the processing workflow and provides an easy entry for new scRNA-seq users. For complete details on the use and execution of this protocol, please refer to Shahan et al. (2022).

Identifying interpretable gene-biomarker associations with functionally informed kernel-based tests in 190,000 exomes.

Here we present an exome-wide rare genetic variant association study for 30 blood biomarkers in 191,971 individuals in the UK Biobank. We compare gene-based association tests for separate functional variant categories to increase interpretability and identify 193 significant gene-biomarker associations. Genes associated with biomarkers were ~ 4.5-fold enriched for conferring Mendelian disorders. In addition to performing weighted gene-based variant collapsing tests, we design and apply variant-category-specific kernel-based tests that integrate quantitative functional variant effect predictions for missense variants, splicing and the binding of RNA-binding proteins. For these tests, we present a computationally efficient combination of the likelihood-ratio and score tests that found 36% more associations than the score test alone while also controlling the type-1 error. Kernel-based tests identified 13% more associations than their gene-based collapsing counterparts and had advantages in the presence of gain of function missense variants. We introduce local collapsing by amino acid position for missense variants and use it to interpret associations and identify potential novel gain of function variants in PIEZO1. Our results show the benefits of investigating different functional mechanisms when performing rare-variant association tests, and demonstrate pervasive rare-variant contribution to biomarker variability.

Control of immediate early gene expression by CPEB4-repressor complex-mediated mRNA degradation.

BACKGROUND: Cytoplasmic polyadenylation element-binding protein 4 (CPEB4) is known to associate with cytoplasmic polyadenylation elements (CPEs) located in the 3' untranslated region (UTR) of specific mRNAs and assemble an activator complex promoting the translation of target mRNAs through cytoplasmic polyadenylation. RESULTS: Here, we find that CPEB4 is part of an alternative repressor complex that mediates mRNA degradation by associating with the evolutionarily conserved CCR4-NOT deadenylase complex. We identify human CPEB4 as an RNA-binding protein (RBP) with enhanced association to poly(A) RNA upon inhibition of class I histone deacetylases (HDACs), a condition known to cause widespread degradation of poly(A)-containing mRNA. Photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) analysis using endogenously tagged CPEB4 in HeLa cells reveals that CPEB4 preferentially binds to the 3'UTR of immediate early gene mRNAs, at G-containing variants of the canonical U- and A-rich CPE located in close proximity to poly(A) sites. By transcriptome-wide mRNA decay measurements, we find that the strength of CPEB4 binding correlates with short mRNA half-lives and that loss of CPEB4 expression leads to the stabilization of immediate early gene mRNAs. Akin to CPEB4, we demonstrate that CPEB1 and CPEB2 also confer mRNA instability by recruitment of the CCR4-NOT complex. CONCLUSIONS: While CPEB4 was previously known for its ability to stimulate cytoplasmic polyadenylation, our findings establish an additional function for CPEB4 as the RNA adaptor of a repressor complex that enhances the degradation of short-lived immediate early gene mRNAs.

Cluster-independent marker feature identification from single-cell omics data using SEMITONES.

Identification of cell identity markers is an essential step in single-cell omics data analysis. Current marker identification strategies typically rely on cluster assignments of cells. However, cluster assignment, particularly for developmental data, is nontrivial, potentially arbitrary, and commonly relies on prior knowledge. In response, we present SEMITONES, a principled method for cluster-free marker identification. We showcase and evaluate its application for marker gene and regulatory region identification from single-cell data of the human haematopoietic system. Additionally, we illustrate its application to spatial transcriptomics data and show how SEMITONES can be used for the annotation of cells given known marker genes. Using several simulated and curated data sets, we demonstrate that SEMITONES qualitatively and quantitatively outperforms existing methods for the retrieval of cell identity markers from single-cell omics data.

Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis-regulatory elements.

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

The BTB transcription factors ZBTB11 and ZFP131 maintain pluripotency by repressing pro-differentiation genes.

In pluripotent cells, a delicate activation-repression balance maintains pro-differentiation genes ready for rapid activation. The identity of transcription factors (TFs) that specifically repress pro-differentiation genes remains obscure. By targeting ∼1,700 TFs with CRISPR loss-of-function screen, we found that ZBTB11 and ZFP131 are required for embryonic stem cell (ESC) pluripotency. ESCs without ZBTB11 or ZFP131 lose colony morphology, reduce proliferation rate, and upregulate transcription of genes associated with three germ layers. ZBTB11 and ZFP131 bind proximally to pro-differentiation genes. ZBTB11 or ZFP131 loss leads to an increase in H3K4me3, negative elongation factor (NELF) complex release, and concomitant transcription at associated genes. Together, our results suggest that ZBTB11 and ZFP131 maintain pluripotency by preventing premature expression of pro-differentiation genes and present a generalizable framework to maintain cellular potency.

A single-cell Arabidopsis root atlas reveals developmental trajectories in wild-type and cell identity mutants.

In all multicellular organisms, transcriptional networks orchestrate organ development. The Arabidopsis root, with its simple structure and indeterminate growth, is an ideal model for investigating the spatiotemporal transcriptional signatures underlying developmental trajectories. To map gene expression dynamics across root cell types and developmental time, we built a comprehensive, organ-scale atlas at single-cell resolution. In addition to estimating developmental progressions in pseudotime, we employed the mathematical concept of optimal transport to infer developmental trajectories and identify their underlying regulators. To demonstrate the utility of the atlas to interpret new datasets, we profiled mutants for two key transcriptional regulators at single-cell resolution, shortroot and scarecrow. We report transcriptomic and in vivo evidence for tissue trans-differentiation underlying a mixed cell identity phenotype in scarecrow. Our results support the atlas as a rich community resource for unraveling the transcriptional programs that specify and maintain cell identity to regulate spatiotemporal organ development.

Deep learning for prediction of population health costs.

BACKGROUND: Accurate prediction of healthcare costs is important for optimally managing health costs. However, methods leveraging the medical richness from data such as health insurance claims or electronic health records are missing. METHODS: Here, we developed a deep neural network to predict future cost from health insurance claims records. We applied the deep network and a ridge regression model to a sample of 1.4 million German insurants to predict total one-year health care costs. Both methods were compared to existing models with various performance measures and were also used to predict patients with a change in costs and to identify relevant codes for this prediction. RESULTS: We showed that the neural network outperformed the ridge regression as well as all considered models for cost prediction. Further, the neural network was superior to ridge regression in predicting patients with cost change and identified more specific codes. CONCLUSION: In summary, we showed that our deep neural network can leverage the full complexity of the patient records and outperforms standard approaches. We suggest that the better performance is due to the ability to incorporate complex interactions in the model and that the model might also be used for predicting other health phenotypes.

Intricacies of single-cell multi-omics data integration.

A wealth of single-cell protocols makes it possible to characterize different molecular layers at unprecedented resolution. Integrating the resulting multimodal single-cell data to find cell-to-cell correspondences remains a challenge. We argue that data integration needs to happen at a meaningful biological level of abstraction and that it is necessary to consider the inherent discrepancies between modalities to strike a balance between biological discovery and noise removal. A survey of current methods reveals that a distinction between technical and biological origins of presumed unwanted variation between datasets is not yet commonly considered. The increasing availability of paired multimodal data will aid the development of improved methods by providing a ground truth on cell-to-cell matches.

Single-cell-resolved dynamics of chromatin architecture delineate cell and regulatory states in zebrafish embryos.

DNA accessibility of cis-regulatory elements (CREs) dictates transcriptional activity and drives cell differentiation during development. While many genes regulating embryonic development have been identified, the underlying CRE dynamics controlling their expression remain largely uncharacterized. To address this, we produced a multimodal resource and genomic regulatory map for the zebrafish community, which integrates single-cell combinatorial indexing assay for transposase-accessible chromatin with high-throughput sequencing (sci-ATAC-seq) with bulk histone PTMs and Hi-C data to achieve a genome-wide classification of the regulatory architecture determining transcriptional activity in the 24-h post-fertilization (hpf) embryo. We characterized the genome-wide chromatin architecture at bulk and single-cell resolution, applying sci-ATAC-seq on whole 24-hpf stage zebrafish embryos, generating accessibility profiles for ∼23,000 single nuclei. We developed a genome segmentation method, ScregSeg (single-cell regulatory landscape segmentation), for defining regulatory programs, and candidate CREs, specific to one or more cell types. We integrated the ScregSeg output with bulk measurements for histone post-translational modifications and 3D genome organization and identified new regulatory principles between chromatin modalities prevalent during zebrafish development. Sci-ATAC-seq profiling of npas4l/cloche mutant embryos identified novel cellular roles for this hematovascular transcriptional master regulator and suggests an intricate mechanism regulating its expression. Our work defines regulatory architecture and principles in the zebrafish embryo and establishes a resource of cell-type-specific genome-wide regulatory annotations and candidate CREs, providing a valuable open resource for genomics, developmental, molecular, and computational biology.