Software for computing and annotating genomic ranges

Significant improvement of cardiac outflow tract septation defects in a DiGeorge syndrome model after minoxidil treatment

The T-BOX transcription factor TBX1 is essential for the development of the pharyngeal apparatus and it is haploinsufficient in DiGeorge syndrome (DGS), a developmental anomaly associated with congenital heart disease and other abnormalities. The murine model recapitulates the heart phenotype and showed collagen accumulation. We first used a cellular model to study gene expression during cardiogenic differentiation of WT and Tbx1-/- mouse embryonic stem cells. Then we used a mouse model of DGS to test whether interfering with collagen accumulation using an inhibitor of lysyl hydroxylase would modify the cardiac phenotype of the mutant. We found that loss of Tbx1 in a precardiac differentiation model was associated with up regulation of a subset of ECM-related genes, including several collagen genes. In the in vivo model, early prenatal treatment with Minoxidil, a lysyl hydroxylase inhibitor, ameliorated the cardiac outflow tract septation phenotype in Tbx1 mutant fetuses, but it had no effect on septation in WT fetuses. We conclude that TBX1 suppresses a defined subset of ECM-related genes. This function is critical for OFT septation because the inhibition of collagen cross-linking in the mutant reduces significantly the penetrance of septation defects.

Reduced protein-coding transcript diversity in severe dengue emphasises the role of alternative splicing

Dengue fever, a neglected tropical arboviral disease, has emerged as a global health concern in the past decade. Necessitating a nuanced comprehension of the intricate dynamics of host-virus interactions influencing disease severity, we analysed transcriptomic patterns using bulk RNA-seq from 112 age- and gender-matched NS1 antigen-confirmed hospital-admitted dengue patients with varying severity. Severe cases exhibited reduced platelet count, increased lymphocytosis, and neutropenia, indicating a dysregulated immune response. Using bulk RNA-seq, our analysis revealed a minimal overlap between the differentially expressed gene and transcript isoform, with a distinct expression pattern across the disease severity. Severe patients showed enrichment in retained intron and nonsense-mediated decay transcript biotypes, suggesting altered splicing efficiency. Furthermore, an up-regulated programmed cell death, a haemolytic response, and an impaired interferon and antiviral response at the transcript level were observed. We also identified the potential involvement of the RBM39 gene among others in the innate immune response during dengue viral pathogenesis, warranting further investigation. These findings provide valuable insights into potential therapeutic targets, underscoring the importance of exploring transcriptomic landscapes between different disease sub-phenotypes in infectious diseases.

The importance of DNA sequence for nucleosome positioning in transcriptional regulation

Nucleosome positioning is a key factor for transcriptional regulation. Nucleosomes regulate the dynamic accessibility of chromatin and interact with the transcription machinery at every stage. Influences to steer nucleosome positioning are diverse, and the according importance of the DNA sequence in contrast to active chromatin remodeling has been the subject of long discussion. In this study, we evaluate the functional role of DNA sequence for all major elements along the process of transcription. We developed a random forest classifier based on local DNA structure that assesses the sequence-intrinsic support for nucleosome positioning. On this basis, we created a simple data resource that we applied genome-wide to the human genome. In our comprehensive analysis, we found a special role of DNA in mediating the competition of nucleosomes with cis-regulatory elements, in enabling steady transcription, for positioning of stable nucleosomes in exons, and for repelling nucleosomes during transcription termination. In contrast, we relate these findings to concurrent processes that generate strongly positioned nucleosomes in vivo that are not mediated by sequence, such as energy-dependent remodeling of chromatin.

Act1 drives chemoresistance via regulation of antioxidant RNA metabolism and redox homeostasis

The IL-17 receptor adaptor molecule Act1, an RNA-binding protein, plays a critical role in IL-17-mediated cancer progression. Here, we report a novel mechanism of how IL-17/Act1 induces chemoresistance by modulating redox homeostasis through epitranscriptomic regulation of antioxidant RNA metabolism. Transcriptome-wide mapping of direct Act1-RNA interactions revealed that Act1 binds to the 5'UTR of antioxidant mRNAs and Wilms' tumor 1-associating protein (WTAP), a key regulator in m6A methyltransferase complex. Strikingly, Act1's binding sites are located in proximity to m6A modification sites, which allows Act1 to promote the recruitment of elF3G for cap-independent translation. Loss of Act1's RNA binding activity or Wtap knockdown abolished IL-17-induced m6A modification and translation of Wtap and antioxidant mRNAs, indicating a feedforward mechanism of the Act1-WTAP loop. We then developed antisense oligonucleotides (Wtap ASO) that specifically disrupt Act1's binding to Wtap mRNA, abolishing IL-17/Act1-WTAP-mediated antioxidant protein production during chemotherapy. Wtap ASO substantially increased the antitumor efficacy of cisplatin, demonstrating a potential therapeutic strategy for chemoresistance.

Cell type- and transcription-independent spatial proximity between enhancers and promoters

Cell type-specific enhancers are critically important for lineage specification. The mechanisms that determine cell-type specificity of enhancer activity, however, are not fully understood. Most current models for how enhancers function invoke physical proximity between enhancer elements and their target genes. Here, we use an imaging-based approach to examine the spatial relationship of cell type-specific enhancers and their target genes with single-cell resolution. Using high-throughput microscopy, we measure the spatial distance from target promoters to their cell type-specific active and inactive enhancers in individual pancreatic cells derived from distinct lineages. We find increased proximity of all promoter-enhancer pairs relative to non-enhancer pairs separated by similar genomic distances. Strikingly, spatial proximity between enhancers and target genes was unrelated to tissue-specific enhancer activity. Furthermore, promoter-enhancer proximity did not correlate with the expression status of target genes. Our results suggest that promoter-enhancer pairs exist in a distinctive chromatin environment but that genome folding is not a universal driver of cell-type specificity in enhancer function.

DNA mismatch and damage patterns revealed by single-molecule sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other diseases1,2. Most mutations begin as nucleotide mismatches or damage in one of the two strands of the DNA before becoming double-strand mutations if unrepaired or misrepaired3,4. However, current DNA-sequencing technologies cannot accurately resolve these initial single-strand events. Here we develop a single-molecule, long-read sequencing method (Hairpin Duplex Enhanced Fidelity sequencing (HiDEF-seq)) that achieves single-molecule fidelity for base substitutions when present in either one or both DNA strands. HiDEF-seq also detects cytosine deamination-a common type of DNA damage-with single-molecule fidelity. We profiled 134 samples from diverse tissues, including from individuals with cancer predisposition syndromes, and derive from them single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumours deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples that are deficient in only polymerase proofreading. We also define a single-strand damage signature for APOBEC3A. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. As double-strand DNA mutations are only the end point of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable studies of how mutations arise in a variety of contexts, especially in cancer and ageing.

Deletion of the mRNA endonuclease Regnase-1 promotes NK cell anti-tumor activity via OCT2-dependent transcription of Ifng

Limited infiltration and activity of natural killer (NK) and T cells within the tumor microenvironment (TME) correlate with poor immunotherapy responses. Here, we examined the role of the endonuclease Regnase-1 on NK cell anti-tumor activity. NK cell-specific deletion of Regnase-1 (Reg1ΔNK) augmented cytolytic activity and interferon-gamma (IFN-γ) production in vitro and increased intra-tumoral accumulation of Reg1ΔNK-NK cells in vivo, reducing tumor growth dependent on IFN-γ. Transcriptional changes in Reg1ΔNK-NK cells included elevated IFN-γ expression, cytolytic effectors, and the chemokine receptor CXCR6. IFN-γ induced expression of the CXCR6 ligand CXCL16 on myeloid cells, promoting further recruitment of Reg1ΔNK-NK cells. Mechanistically, Regnase-1 deletion increased its targets, the transcriptional regulators OCT2 and IκBζ, following interleukin (IL)-12 and IL-18 stimulation, and the resulting OCT2-IκBζ-NF-κB complex induced Ifng transcription. Silencing Regnase-1 in human NK cells increased the expression of IFNG and POU2F2. Our findings highlight NK cell dysfunction in the TME and propose that targeting Regnase-1 could augment active NK cell persistence for cancer immunotherapy.

Differential leaf flooding resilience in Arabidopsis thaliana is controlled by ethylene signaling-activated and age-dependent phosphorylation of ORESARA1

The phytohormone ethylene is a major regulator of plant adaptive responses to flooding. In flooded plant tissues, ethylene quickly increases to high concentrations owing to its low solubility and diffusion rates in water. Ethylene accumulation in submerged plant tissues makes it a reliable cue for triggering flood acclimation responses, including metabolic adjustments to cope with flood-induced hypoxia. However, persistent ethylene accumulation also accelerates leaf senescence. Stress-induced senescence hampers photosynthetic capacity and stress recovery. In submerged Arabidopsis, senescence follows a strict age-dependent pattern starting with the older leaves. Although mechanisms underlying ethylene-mediated senescence have been uncovered, it is unclear how submerged plants avoid indiscriminate breakdown of leaves despite high systemic ethylene accumulation. We demonstrate that although submergence triggers leaf-age-independent activation of ethylene signaling via EIN3 in Arabidopsis, senescence is initiated only in old leaves. EIN3 stabilization also leads to overall transcript and protein accumulation of the senescence-promoting transcription factor ORESARA1 (ORE1) in both old and young leaves during submergence. However, leaf-age-dependent senescence can be explained by ORE1 protein activation via phosphorylation specifically in old leaves, independent of the previously identified age-dependent control of ORE1 via miR164. A systematic analysis of the roles of the major flooding stress cues and signaling pathways shows that only the combination of ethylene and darkness is sufficient to mimic submergence-induced senescence involving ORE1 accumulation and phosphorylation. Hypoxia, most often associated with flooding stress in plants, appears to have no role in these processes. Our results reveal a mechanism by which plants regulate the speed and pattern of senescence during environmental stresses such as flooding. Age-dependent ORE1 activity ensures that older, expendable leaves are dismantled first, thus prolonging the life of younger leaves and meristematic tissues that are vital to whole-plant survival.

Evidence for gene flow and trait reversal during radiation of Mexican Goodeid fish

Understanding the phylogeographic history of a group and identifying the factors contributing to speciation is an important challenge in evolutionary biology. The Goodeinae are a group of live-bearing fishes endemic to Mexico. Here, we develop genomic resources for species within the Goodeinae and use phylogenomic approaches to characterise their evolutionary history. We sequenced, assembled and annotated the genomes of four Goodeinae species, including Ataeniobius toweri, the only matrotrophic live-bearing fish without a trophotaenia in the group. We estimated timings of species divergence and examined the extent and timing of introgression between the species to assess if this may have occurred during an early radiation, or in more recent episodes of secondary contact. We used branch-site models to detect genome-wide positive selection across Goodeinae, and we specifically asked whether this differs in A. toweri, where loss of placental viviparity has recently occurred. We found evidence of gene flow between geographically isolated species, suggesting vicariant speciation was supplemented by limited post-speciation gene flow, and gene flow may explain previous uncertainties about Goodeid phylogeny. Genes under positive selection in the group are likely to be associated with the switch to live-bearing. Overall, our studies suggest that both volcanism-driven vicariance and changes in reproductive mode influenced radiation in the Goodeinae.

Uncovering genes involved in pollinator-driven mating system shifts and selfing syndrome evolution in Brassica rapa

Shifts in pollinator occurrence and their pollen transport effectiveness drive the evolution of mating systems in flowering plants. Understanding the genomic basis of these changes is essential for predicting the persistence of a species under environmental changes. We investigated the genomic changes in Brassica rapa over nine generations of pollination by hoverflies associated with rapid morphological evolution toward the selfing syndrome. We combined a genotyping-by-sequencing (GBS) approach with a genome-wide association study (GWAS) to identify candidate genes, and assessed their functional role in the observed morphological changes by studying mutations of orthologous genes in the model plant Arabidopsis thaliana. We found 31 candidate genes involved in a wide range of functions from DNA/RNA binding to transport. Our functional assessment of orthologous genes in A. thaliana revealed that two of the identified genes in B. rapa are involved in regulating the size of floral organs. We found a protein kinase superfamily protein involved in petal width, an important trait in plant attractiveness to pollinators. Moreover, we found a histone lysine methyltransferase (HKMT) associated with stamen length. Altogether, our study shows that hoverfly pollination leads to rapid evolution toward the selfing syndrome mediated by polygenic changes.

DBD-α4 helix of EWSR1::FLI1 is required for GGAA microsatellite binding that underlies genome regulation in Ewing sarcoma

Ewing sarcoma is the second most common bone cancer in children and young adults. In 85% of patients, a translocation between chromosomes 11 and 22 results in a potent fusion oncoprotein, EWSR1::FLI1. EWSR1::FLI1 is the only genetic alteration in an otherwise unaltered genome of Ewing sarcoma tumors. The EWSR1 portion of the protein is an intrinsically disordered domain involved in transcriptional regulation by EWSR1::FLI1. The FLI portion of the fusion contains a DNA binding domain shown to bind core GGAA motifs and GGAA repeats. A small alpha-helix in the DNA binding domain of FLI1, DBD-𝛼4 helix, is critical for the transcription function of EWSR1::FLI1. In this study, we aimed to understand the mechanism by which the DBD-𝛼4 helix promotes transcription, and therefore oncogenic transformation. We utilized a multi-omics approach to assess chromatin organization, active chromatinmarks, genome binding, and gene expression in cells expressing EWSR1::FLI1 constructs with and without the DBD-𝛼4 helix. Our studies revealed DBD-𝛼4 helix is crucial for cooperative binding of EWSR1::FLI1 at GGAA microsatellites. This binding underlies many aspects of genome regulation by EWSR1::FLI1 such as formation of TADs, chromatin loops, enhancers and productive transcription hubs.

Genome-wide DNA methylation changes in human spermatogenesis

Sperm production and function require the correct establishment of DNA methylation patterns in the germline. Here, we examined the genome-wide DNA methylation changes during human spermatogenesis and its alterations in disturbed spermatogenesis. We found that spermatogenesis is associated with remodeling of the methylome, comprising a global decline in DNA methylation in primary spermatocytes followed by selective remethylation, resulting in a spermatids/sperm-specific methylome. Hypomethylated regions in spermatids/sperm were enriched in specific transcription factor binding sites for DMRT and SOX family members and spermatid-specific genes. Intriguingly, while SINEs displayed differential methylation throughout spermatogenesis, LINEs appeared to be protected from changes in DNA methylation. In disturbed spermatogenesis, germ cells exhibited considerable DNA methylation changes, which were significantly enriched at transposable elements and genes involved in spermatogenesis. We detected hypomethylation in SVA and L1HS in disturbed spermatogenesis, suggesting an association between the abnormal programming of these regions and failure of germ cells progressing beyond meiosis.

Pre-hypertrophic chondrogenic enhancer landscape of limb and axial skeleton development

Chondrocyte differentiation controls skeleton development and stature. Here we provide a comprehensive map of chondrocyte-specific enhancers and show that they provide a mechanistic framework through which non-coding genetic variants can influence skeletal development and human stature. Working with fetal chondrocytes isolated from mice bearing a Col2a1 fluorescent regulatory sensor, we identify 780 genes and 2'704 putative enhancers specifically active in chondrocytes using a combination of RNA-seq, ATAC-seq and H3K27ac ChIP-seq. Most of these enhancers (74%) show pan-chondrogenic activity, with smaller populations being restricted to limb (18%) or trunk (8%) chondrocytes only. Notably, genetic variations overlapping these enhancers better explain height differences than those overlapping non-chondrogenic enhancers. Finally, targeted deletions of identified enhancers at the Fgfr3, Col2a1, Hhip and, Nkx3-2 loci confirm their role in regulating cognate genes. This enhancer map provides a framework for understanding how genes and non-coding variations influence bone development and diseases.

A single-cell transcriptome atlas of human euploid and aneuploid blastocysts

Aneuploidy is frequently detected in early human embryos as a major cause of early pregnancy failure. However, how aneuploidy affects cellular function remains elusive. Here, we profiled the transcriptomes of 14,908 single cells from 203 human euploid and aneuploid blastocysts involving autosomal and sex chromosomes. Nearly all of the blastocysts contained four lineages. In aneuploid chromosomes, 19.5% ± 1.2% of the expressed genes showed a dosage effect, and 90 dosage-sensitive domains were identified. Aneuploidy leads to prevalent genome-wide transcriptome alterations. Common effects, including apoptosis, were identified, especially in monosomies, partially explaining the lower cell numbers in autosomal monosomies. We further identified lineage-specific effects causing unstable epiblast development in aneuploidies, which was accompanied by the downregulation of TGF-β and FGF signaling, which resulted in insufficient trophectoderm maturation. Our work provides crucial insights into the molecular basis of human aneuploid blastocysts and may shed light on the cellular interaction during blastocyst development.

Interferon-γ inducible factor 16 (IFI16) restricts adeno-associated virus type 2 (AAV2) transduction in an immune-modulatory independent way

We determined the transcription profile of adeno-associated virus type 2 (AAV2)-infected primary human fibroblasts. Subsequent analysis revealed that cells respond to AAV infection through changes in several significantly affected pathways, including cell cycle regulation, chromatin modulation, and innate immune responses. Various assays were performed to validate selected differentially expressed genes and to confirm not only the quality but also the robustness of the raw data. One of the genes upregulated in AAV2-infected cells was interferon-γ inducible factor 16 (IFI16). IFI16 is known as a multifunctional cytosolic and nuclear innate immune sensor for double-stranded as well as single-stranded DNA, exerting its effects through various mechanisms, such as interferon response, epigenetic modifications, or transcriptional regulation. IFI16 thereby constitutes a restriction factor for many different viruses among them, as shown here, AAV2 and thereof derived vectors. Indeed, the post-transcriptional silencing of IFI16 significantly increased AAV2 transduction efficiency, independent of the structure of the virus/vector genome. We also show that IFI16 exerts its inhibitory effect on AAV2 transduction in an immune-modulatory independent way by interfering with Sp1-dependent transactivation of wild-type AAV2 and AAV2 vector promoters.

IMPORTANCE

Adeno-associated virus (AAV) vectors are among the most frequently used viral vectors for gene therapy. The lack of pathogenicity of the parental virus, the long-term persistence as episomes in non-proliferating cells, and the availability of a variety of AAV serotypes differing in their cellular tropism are advantageous features of this biological nanoparticle. To deepen our understanding of virus-host interactions, especially in terms of antiviral responses, we present here the first transcriptome analysis of AAV serotype 2 (AAV2)-infected human primary fibroblasts. Our findings indicate that interferon-γ inducible factor 16 acts as an antiviral factor in AAV2 infection and AAV2 vector-mediated cell transduction in an immune-modulatory independent way by interrupting the Sp1-dependent gene expression from viral or vector genomes.

Identification and analysis of short-term and long-term salt-associated lncRNAs in the leaf of Avicennia marina

As a highly salt-resistant mangrove, Avicennia marina can thrive in the hypersaline water. The leaves of Avicennia marina play a crucial role in salinity stress adaptability by secreting salt. Although the functions of long non-coding RNAs (lncRNAs) in leaves remain unknown, they have emerged as regulators in leaf development, aging and salt response. In this study, we employed transcriptomic data of both short-term and long-term salt treated leaves to identify salt-associated lncRNAs of leaf tissue. As a result, 687 short-term and 797 long-term salt-associated lncRNAs were identified. Notably, both short-term and long-term salt-associated lncRNAs exhibited slightly longer lengths and larger exons, but smaller introns compared with salt-non-associated lncRNAs. Furthermore, salt-associated lncRNAs also displayed higher tissue-specificity than salt-non-associated lncRNAs. Most of the salt-associated lncRNAs were common to short- and long-term salt treatments. And about one fifth of the downregulated salt-associated lncRNAs identified both in two terms were leaf tissue-specific lncRNAs. Besides, these leaf-specific lncRNAs were found to be involved in the oxidation-reduction and photosynthesis processes, as well as several metabolic processes, suggesting the noticeable functions of salt-associated lncRNAs in regulating salt responses of Avicennia marina leaves.

Predicting gene expression state and prioritizing putative enhancers using 5hmC signal

BACKGROUND

Like its parent base 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC) is a direct epigenetic modification of cytosines in the context of CpG dinucleotides. 5hmC is the most abundant oxidized form of 5mC, generated through the action of TET dioxygenases at gene bodies of actively-transcribed genes and at active or lineage-specific enhancers. Although such enrichments are reported for 5hmC, to date, predictive models of gene expression state or putative regulatory regions for genes using 5hmC have not been developed.

RESULTS

Here, by using only 5hmC enrichment in genic regions and their vicinity, we develop neural network models that predict gene expression state across 49 cell types. We show that our deep neural network models distinguish high vs low expression state utilizing only 5hmC levels and these predictive models generalize to unseen cell types. Further, in order to leverage 5hmC signal in distal enhancers for expression prediction, we employ an Activity-by-Contact model and also develop a graph convolutional neural network model with both utilizing Hi-C data and 5hmC enrichment to prioritize enhancer-promoter links. These approaches identify known and novel putative enhancers for key genes in multiple immune cell subsets.

CONCLUSIONS

Our work highlights the importance of 5hmC in gene regulation through proximal and distal mechanisms and provides a framework to link it to genome function. With the recent advances in 6-letter DNA sequencing by short and long-read techniques, profiling of 5mC and 5hmC may be done routinely in the near future, hence, providing a broad range of applications for the methods developed here.

A repertoire of intronic lariat RNAs reveals tissue-specific regulation and target mimicry potential in plants

Lariat RNA is concomitantly produced by excised intron during RNA splicing, which is usually debranched by DBR1, an RNA debranching enzyme. However, increasing evidence showed that some lariat RNA could escape debranching. Little is known about how and why these lariat RNAs could be retained. By comparing the atlas of lariat RNAs between the non-dividing cell (mature pollen) and three actively dividing tissues (young shoot apex, young seeds, and young roots), we identified hundreds to thousands of lariat RNA naturally retained in each tissue, and the incidence of lariat RNA retention is much less in shoot apex while much more in pollen. Many lariat RNAs derived from the same intron or different lariat RNAs from the same pre-mRNA could be retained in one tissue while degraded in the other tissues. By deciphering lariat RNA sequences, we identified an AG-rich (RAAAAVAAAR) motif and a UC-rich (UCUCUYUCUC) motif for pollen-specific and the other three tissues-retained lariat RNAs, respectively. Reconstitution of the pollen-specific AG-rich motif indeed enhanced lariat RNA retention in plants. Biologically, hundreds of lariat RNAs harbored miRNA binding sites, and dual-luciferase reporter assay showed that these natural lariat RNAs had the potential to protect expression of miRNA target genes. Collectively, our results uncover that selective retention of lariat RNA is an actively regulatory process, and provide new insights into understanding how lariat RNA metabolism may impact miRNA activity.

An integrative epigenome-based strategy for unbiased functional profiling of clinical kinase inhibitors

More than 500 kinases are implicated in the control of most cellular process in mammals, and deregulation of their activity is linked to cancer and inflammatory disorders. 80 clinical kinase inhibitors (CKIs) have been approved for clinical use and hundreds are in various stages of development. However, CKIs inhibit other kinases in addition to the intended target(s), causing both enhanced clinical effects and undesired side effects that are only partially predictable based on in vitro selectivity profiling. Here, we report an integrative approach grounded on the use of chromatin modifications as unbiased, information-rich readouts of the functional effects of CKIs on macrophage activation. This approach exceeded the performance of transcriptome-based approaches and allowed us to identify similarities and differences among CKIs with identical intended targets, to recognize novel CKI specificities and to pinpoint CKIs that may be repurposed to control inflammation, thus supporting the utility of this strategy to improve selection and use of CKIs in clinical settings.

A disease-associated gene desert directs macrophage inflammation through ETS2

Increasing rates of autoimmune and inflammatory disease present a burgeoning threat to human health1. This is compounded by the limited efficacy of available treatments1 and high failure rates during drug development2, highlighting an urgent need to better understand disease mechanisms. Here we show how functional genomics could address this challenge. By investigating an intergenic haplotype on chr21q22-which has been independently linked to inflammatory bowel disease, ankylosing spondylitis, primary sclerosing cholangitis and Takayasu's arteritis3-6-we identify that the causal gene, ETS2, is a central regulator of human inflammatory macrophages and delineate the shared disease mechanism that amplifies ETS2 expression. Genes regulated by ETS2 were prominently expressed in diseased tissues and more enriched for inflammatory bowel disease GWAS hits than most previously described pathways. Overexpressing ETS2 in resting macrophages reproduced the inflammatory state observed in chr21q22-associated diseases, with upregulation of multiple drug targets, including TNF and IL-23. Using a database of cellular signatures7, we identified drugs that might modulate this pathway and validated the potent anti-inflammatory activity of one class of small molecules in vitro and ex vivo. Together, this illustrates the power of functional genomics, applied directly in primary human cells, to identify immune-mediated disease mechanisms and potential therapeutic opportunities.

Aberrant DNA methylation distorts developmental trajectories in atypical teratoid/rhabdoid tumors

Atypical teratoid/rhabdoid tumors (AT/RTs) are pediatric brain tumors known for their aggressiveness and aberrant but still unresolved epigenetic regulation. To better understand their malignancy, we investigated how AT/RT-specific DNA hypermethylation was associated with gene expression and altered transcription factor binding and how it is linked to upstream regulation. Medulloblastomas, choroid plexus tumors, pluripotent stem cells, and fetal brain were used as references. A part of the genomic regions, which were hypermethylated in AT/RTs similarly as in pluripotent stem cells and demethylated in the fetal brain, were targeted by neural transcriptional regulators. AT/RT-unique DNA hypermethylation was associated with polycomb repressive complex 2 and linked to suppressed genes with a role in neural development and tumorigenesis. Activity of the several NEUROG/NEUROD pioneer factors, which are unable to bind to methylated DNA, was compromised via the suppressed expression or DNA hypermethylation of their target sites, which was also experimentally validated for NEUROD1 in medulloblastomas and AT/RT samples. These results highlight and characterize the role of DNA hypermethylation in AT/RT malignancy and halted neural cell differentiation.

Conserved role of hnRNPL in alternative splicing of epigenetic modifiers enables B cell activation

The multifunctional RNA-binding protein hnRNPL is implicated in antibody class switching but its broader function in B cells is unknown. Here, we show that hnRNPL is essential for B cell activation, germinal center formation, and antibody responses. Upon activation, hnRNPL-deficient B cells show proliferation defects and increased apoptosis. Comparative analysis of RNA-seq data from activated B cells and another eight hnRNPL-depleted cell types reveals common effects on MYC and E2F transcriptional programs required for proliferation. Notably, while individual gene expression changes are cell type specific, several alternative splicing events affecting histone modifiers like KDM6A and SIRT1, are conserved across cell types. Moreover, hnRNPL-deficient B cells show global changes in H3K27me3 and H3K9ac. Epigenetic dysregulation after hnRNPL loss could underlie differential gene expression and upregulation of lncRNAs, and explain common and cell type-specific phenotypes, such as dysfunctional mitochondria and ROS overproduction in mouse B cells. Thus, hnRNPL is essential for the resting-to-activated B cell transition by regulating transcriptional programs and metabolism, at least in part through the alternative splicing of several histone modifiers.

BCAT1, IKZF1 and SEPT9: methylated DNA biomarkers for detection of pan-gastrointestinal adenocarcinomas

INTRODUCTION

Methylated circulating tumour DNA (ctDNA) blood tests for BCAT1/IKZF1 (COLVERA) and SEPT9 (Epi proColon) are used to detect colorectal cancer (CRC). However, there are no ctDNA assays approved for other gastrointestinal adenocarcinomas. We aimed to characterize BCAT1, IKZF1 and SEPT9 methylation in different gastrointestinal adenocarcinoma and non-gastrointestinal tumours to determine if these validated CRC biomarkers might be useful for pan-gastrointestinal adenocarcinoma detection.

METHODS

Tissue DNA methylation data from colorectal (COAD, READ), gastroesophageal (ESCA, STAD), pancreatic (PAAD) and cholangiocarcinoma (CHOL) adenocarcinoma cohorts within The Cancer Genome Atlas were used for differential methylation analyses. Clinicodemographic predictors of BCAT1, IKZF1 and SEPT9 methylation, and the selectivity of hypermethylated BCAT1, IKZF1 and SEPT9 for colorectal adenocarcinomas in comparison to other cancers were each explored with beta regression.

RESULTS

Hypermethylated BCAT1, IKZF1 and SEPT9 were each differentially methylated in colorectal and gastroesophageal adenocarcinomas. IKZF1 was differentially methylated in pancreatic adenocarcinoma. Hypermethylated DNA biomarkers BCAT1, IKZF1 and SEPT9 were largely stable across different stages of disease and were highly selective for gastrointestinal adenocarcinomas relative to other cancer types.

DISCUSSION

Existing CRC methylated ctDNA blood tests for BCAT1/IKZF1 and SEPT9 might be usefully repurposed for use in other gastrointestinal adenocarcinomas and warrant further prospective ctDNA studies.

The chromatin landscape of pathogenic transcriptional cell states in rheumatoid arthritis

Synovial tissue inflammation is a hallmark of rheumatoid arthritis (RA). Recent work has identified prominent pathogenic cell states in inflamed RA synovial tissue, such as T peripheral helper cells; however, the epigenetic regulation of these states has yet to be defined. Here, we examine genome-wide open chromatin at single-cell resolution in 30 synovial tissue samples, including 12 samples with transcriptional data in multimodal experiments. We identify 24 chromatin classes and predict their associated transcription factors, including a CD8 + GZMK+ class associated with EOMES and a lining fibroblast class associated with AP-1. By integrating with an RA tissue transcriptional atlas, we propose that these chromatin classes represent 'superstates' corresponding to multiple transcriptional cell states. Finally, we demonstrate the utility of this RA tissue chromatin atlas through the associations between disease phenotypes and chromatin class abundance, as well as the nomination of classes mediating the effects of putatively causal RA genetic variants.

Environment-induced heritable variations are common in Arabidopsis thaliana

Parental or ancestral environments can induce heritable phenotypic changes, but whether such environment-induced heritable changes are a common phenomenon remains unexplored. Here, we subject 14 genotypes of Arabidopsis thaliana to 10 different environmental treatments and observe phenotypic and genome-wide gene expression changes over four successive generations. We find that all treatments caused heritable phenotypic and gene expression changes, with a substantial proportion stably transmitted over all observed generations. Intriguingly, the susceptibility of a genotype to environmental inductions could be predicted based on the transposon abundance in the genome. Our study thus challenges the classic view that the environment only participates in the selection of heritable variation and suggests that the environment can play a significant role in generating of heritable variations.

Chromatin profiling identifies putative dual roles for H3K27me3 in regulating transposons and cell type-specific genes in choanoflagellates

Gene expression is tightly controlled during animal development to allow the formation of specialized cell types. Our understanding of how animals evolved this exquisite regulatory control remains elusive, but evidence suggests that changes in chromatin-based mechanisms may have contributed. To investigate this possibility, here we examine chromatin-based gene regulatory features in the closest relatives of animals, choanoflagellates. Using Salpingoeca rosetta as a model system, we examined chromatin accessibility and histone modifications at the genome scale and compared these features to gene expression. We first observed that accessible regions of chromatin are primarily associated with gene promoters and found no evidence of distal gene regulatory elements resembling the enhancers that animals deploy to regulate developmental gene expression. Remarkably, a histone modification deposited by polycomb repressive complex 2, histone H3 lysine 27 trimethylation (H3K27me3), appeared to function similarly in S. rosetta to its role in animals, because this modification decorated genes with cell type-specific expression. Additionally, H3K27me3 marked transposons, retaining what appears to be an ancestral role in regulating these elements. We further uncovered a putative new bivalent chromatin state at cell type-specific genes that consists of H3K27me3 and histone H3 lysine 4 mono-methylation (H3K4me1). Together, our discoveries support the scenario that gene-associated histone modification states that underpin development emerged before the evolution of animal multicellularity.

An amphioxus neurula stage cell atlas supports a complex scenario for the emergence of vertebrate head mesoderm

The emergence of new structures can often be linked to the evolution of novel cell types that follows the rewiring of developmental gene regulatory subnetworks. Vertebrates are characterized by a complex body plan compared to the other chordate clades and the question remains of whether and how the emergence of vertebrate morphological innovations can be related to the appearance of new embryonic cell populations. We previously proposed, by studying mesoderm development in the cephalochordate amphioxus, a scenario for the evolution of the vertebrate head mesoderm. To further test this scenario at the cell population level, we used scRNA-seq to construct a cell atlas of the amphioxus neurula, stage at which the main mesodermal compartments are specified. Our data allowed us to validate the presence of a prechordal-plate like territory in amphioxus. Additionally, the transcriptomic profile of somite cell populations supports the homology between specific territories of amphioxus somites and vertebrate cranial/pharyngeal and lateral plate mesoderm. Finally, our work provides evidence that the appearance of the specific mesodermal structures of the vertebrate head was associated to both segregation of pre-existing cell populations, and co-option of new genes for the control of myogenesis.

Prediction of tumor-specific splicing from somatic mutations as a source of neoantigen candidates

Motivation

Neoantigens are promising targets for cancer immunotherapies and might arise from alternative splicing. However, detecting tumor-specific splicing is challenging because many non-canonical splice junctions identified in tumors also appear in healthy tissues. To increase tumor-specificity, we focused on splicing caused by somatic mutations as a source for neoantigen candidates in individual patients.

Results

We developed the tool splice2neo with multiple functionalities to integrate predicted splice effects from somatic mutations with splice junctions detected in tumor RNA-seq and to annotate the resulting transcript and peptide sequences. Additionally, we provide the tool EasyQuant for targeted RNA-seq read mapping to candidate splice junctions. Using a stringent detection rule, we predicted 1.7 splice junctions per patient as splice targets with a false discovery rate below 5% in a melanoma cohort. We confirmed tumor-specificity using independent, healthy tissue samples. Furthermore, using tumor-derived RNA, we confirmed individual exon-skipping events experimentally. Most target splice junctions encoded neoepitope candidates with predicted major histocompatibility complex (MHC)-I or MHC-II binding. Compared to neoepitope candidates from non-synonymous point mutations, the splicing-derived MHC-I neoepitope candidates had lower self-similarity to corresponding wild-type peptides. In conclusion, we demonstrate that identifying mutation-derived, tumor-specific splice junctions can lead to additional neoantigen candidates to expand the target repertoire for cancer immunotherapies.

Availability and implementation

The R package splice2neo and the python package EasyQuant are available at https://github.com/TRON-Bioinformatics/splice2neo and https://github.com/TRON-Bioinformatics/easyquant, respectively.

The landscape of transcription factor promoter activity during vegetative development in Marchantia

Transcription factors (TFs) are essential for the regulation of gene expression and cell fate determination. Characterizing the transcriptional activity of TF genes in space and time is a critical step toward understanding complex biological systems. The vegetative gametophyte meristems of bryophytes share some characteristics with the shoot apical meristems of flowering plants. However, the identity and expression profiles of TFs associated with gametophyte organization are largely unknown. With only ∼450 putative TF genes, Marchantia (Marchantia polymorpha) is an outstanding model system for plant systems biology. We have generated a near-complete collection of promoter elements derived from Marchantia TF genes. We experimentally tested reporter fusions for all the TF promoters in the collection and systematically analyzed expression patterns in Marchantia gemmae. This allowed us to build a map of expression domains in early vegetative development and identify a set of TF-derived promoters that are active in the stem-cell zone. The cell markers provide additional tools and insight into the dynamic regulation of the gametophytic meristem and its evolution. In addition, we provide an online database of expression patterns for all promoters in the collection. We expect that these promoter elements will be useful for cell-type-specific expression, synthetic biology applications, and functional genomics.

Methylation of the chromatin modifier KMT2D by SMYD2 contributes to therapeutic response in hormone-dependent breast cancer

Activating mutations in PIK3CA are frequently found in estrogen-receptor-positive (ER+) breast cancer, and the combination of the phosphatidylinositol 3-kinase (PI3K) inhibitor alpelisib with anti-ER inhibitors is approved for therapy. We have previously demonstrated that the PI3K pathway regulates ER activity through phosphorylation of the chromatin modifier KMT2D. Here, we discovered a methylation site on KMT2D, at K1330 directly adjacent to S1331, catalyzed by the lysine methyltransferase SMYD2. SMYD2 loss attenuates alpelisib-induced KMT2D chromatin binding and alpelisib-mediated changes in gene expression, including ER-dependent transcription. Knockdown or pharmacological inhibition of SMYD2 sensitizes breast cancer cells, patient-derived organoids, and tumors to PI3K/AKT inhibition and endocrine therapy in part through KMT2D K1330 methylation. Together, our findings uncover a regulatory crosstalk between post-translational modifications that fine-tunes KMT2D function at the chromatin. This provides a rationale for the use of SMYD2 inhibitors in combination with PI3Kα/AKT inhibitors in the treatment of ER+/PIK3CA mutant breast cancer.

Repression of pervasive antisense transcription is the primary role of fission yeast RNA polymerase II CTD serine 2 phosphorylation

The RNA polymerase II carboxy-terminal domain (CTD) consists of conserved heptapeptide repeats that can be phosphorylated to influence distinct stages of the transcription cycle, including RNA processing. Although CTD-associated proteins have been identified, phospho-dependent CTD interactions have remained elusive. Proximity-dependent biotinylation (PDB) has recently emerged as an alternative approach to identify protein-protein associations in the native cellular environment. In this study, we present a PDB-based map of the fission yeast RNAPII CTD interactome in living cells and identify phospho-dependent CTD interactions by using a mutant in which Ser2 was replaced by alanine in every repeat of the fission yeast CTD. This approach revealed that CTD Ser2 phosphorylation is critical for the association between RNAPII and the histone methyltransferase Set2 during transcription elongation, but is not required for 3' end processing and transcription termination. Accordingly, loss of CTD Ser2 phosphorylation causes a global increase in antisense transcription, correlating with elevated histone acetylation in gene bodies. Our findings reveal that the fundamental role of CTD Ser2 phosphorylation is to establish a chromatin-based repressive state that prevents cryptic intragenic transcription initiation.

Organization of corticocortical and thalamocortical top-down inputs in the primary visual cortex

Unified visual perception requires integration of bottom-up and top-down inputs in the primary visual cortex (V1), yet the organization of top-down inputs in V1 remains unclear. Here, we used optogenetics-assisted circuit mapping to identify how multiple top-down inputs from higher-order cortical and thalamic areas engage V1 excitatory and inhibitory neurons. Top-down inputs overlap in superficial layers yet segregate in deep layers. Inputs from the medial secondary visual cortex (V2M) and anterior cingulate cortex (ACA) converge on L6 Pyrs, whereas ventrolateral orbitofrontal cortex (ORBvl) and lateral posterior thalamic nucleus (LP) inputs are processed in parallel in Pyr-type-specific subnetworks (Pyr←ORBvl and Pyr←LP) and drive mutual inhibition between them via local interneurons. Our study deepens understanding of the top-down modulation mechanisms of visual processing and establishes that V2M and ACA inputs in L6 employ integrated processing distinct from the parallel processing of LP and ORBvl inputs in L5.

Effect of genomic and cellular environments on gene expression noise

BACKGROUND

Individual cells from isogenic populations often display large cell-to-cell differences in gene expression. This "noise" in expression derives from several sources, including the genomic and cellular environment in which a gene resides. Large-scale maps of genomic environments have revealed the effects of epigenetic modifications and transcription factor occupancy on mean expression levels, but leveraging such maps to explain expression noise will require new methods to assay how expression noise changes at locations across the genome.

RESULTS

To address this gap, we present Single-cell Analysis of Reporter Gene Expression Noise and Transcriptome (SARGENT), a method that simultaneously measures the noisiness of reporter genes integrated throughout the genome and the global mRNA profiles of individual reporter-gene-containing cells. Using SARGENT, we perform the first comprehensive genome-wide survey of how genomic locations impact gene expression noise. We find that the mean and noise of expression correlate with different histone modifications. We quantify the intrinsic and extrinsic components of reporter gene noise and, using the associated mRNA profiles, assign the extrinsic component to differences between the CD24+ "stem-like" substate and the more "differentiated" substate. SARGENT also reveals the effects of transgene integrations on endogenous gene expression, which will help guide the search for "safe-harbor" loci.

CONCLUSIONS

Taken together, we show that SARGENT is a powerful tool to measure both the mean and noise of gene expression at locations across the genome and that the data generatd by SARGENT reveals important insights into the regulation of gene expression noise genome-wide.

A nuclear genome assembly of an extinct flightless bird, the little bush moa

We present a draft genome of the little bush moa (Anomalopteryx didiformis)-one of approximately nine species of extinct flightless birds from Aotearoa, New Zealand-using ancient DNA recovered from a fossil bone from the South Island. We recover a complete mitochondrial genome at 249.9× depth of coverage and almost 900 megabases of a male moa nuclear genome at ~4 to 5× coverage, with sequence contiguity sufficient to identify more than 85% of avian universal single-copy orthologs. We describe a diverse landscape of transposable elements and satellite repeats, estimate a long-term effective population size of ~240,000, identify a diverse suite of olfactory receptor genes and an opsin repertoire with sensitivity in the ultraviolet range, show that the wingless moa phenotype is likely not attributable to gene loss or pseudogenization, and identify potential function-altering coding sequence variants in moa that could be synthesized for future functional assays. This genomic resource should support further studies of avian evolution and morphological divergence.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

A high-resolution map of functional miR-181 response elements in the thymus reveals the role of coding sequence targeting and an alternative seed match

MicroRNAs (miRNAs) are critical post-transcriptional regulators in many biological processes. They act by guiding RNA-induced silencing complexes to miRNA response elements (MREs) in target mRNAs, inducing translational inhibition and/or mRNA degradation. Functional MREs are expected to predominantly occur in the 3' untranslated region and involve perfect base-pairing of the miRNA seed. Here, we generate a high-resolution map of miR-181a/b-1 (miR-181) MREs to define the targeting rules of miR-181 in developing murine T cells. By combining a multi-omics approach with computational high-resolution analyses, we uncover novel miR-181 targets and demonstrate that miR-181 acts predominantly through RNA destabilization. Importantly, we discover an alternative seed match and identify a distinct set of targets with repeat elements in the coding sequence which are targeted by miR-181 and mediate translational inhibition. In conclusion, deep profiling of MREs in primary cells is critical to expand physiologically relevant targetomes and establish context-dependent miRNA targeting rules.

p300/CBP degradation is required to disable the active AR enhanceosome in prostate cancer

Prostate cancer is an exemplar of an enhancer-binding transcription factor-driven disease. The androgen receptor (AR) enhanceosome complex comprised of chromatin and epigenetic coregulators assembles at enhancer elements to drive disease progression. The paralog lysine acetyltransferases p300 and CBP deposit histone marks that are associated with enhancer activation. Here, we demonstrate that p300/CBP are determinant cofactors of the active AR enhanceosome in prostate cancer. Histone H2B N-terminus multisite lysine acetylation (H2BNTac), which is exclusively reliant on p300/CBP catalytic function, marked active enhancers and was notably elevated in prostate cancer lesions relative to the adjacent benign epithelia. Degradation of p300/CBP rapidly depleted acetylation marks associated with the active AR enhanceosome, which was only partially phenocopied by inhibition of their reader bromodomains. Notably, H2BNTac was effectively abrogated only upon p300/CBP degradation, which led to a stronger suppression of p300/CBP-dependent oncogenic gene programs relative to bromodomain inhibition or the inhibition of its catalytic domain. In vivo experiments using an orally active p300/CBP proteolysis targeting chimera (PROTAC) degrader (CBPD-409) showed that p300/CBP degradation potently inhibited tumor growth in preclinical models of castration-resistant prostate cancer and synergized with AR antagonists. While mouse p300/CBP orthologs were effectively degraded in host tissues, prolonged treatment with the PROTAC degrader was well tolerated with no significant signs of toxicity. Taken together, our study highlights the pivotal role of p300/CBP in maintaining the active AR enhanceosome and demonstrates how target degradation may have functionally distinct effects relative to target inhibition, thus supporting the development of p300/CBP degraders for the treatment of advanced prostate cancer.

Exploring the gene expression network involved in the heat stress response of a thermotolerant tomato genotype

BACKGROUND

The increase in temperatures due to the current climate change dramatically affects crop cultivation, resulting in yield losses and altered fruit quality. Tomato is one of the most extensively grown and consumed horticultural products, and although it can withstand a wide range of climatic conditions, heat stress can affect plant growth and development specially on the reproductive stage, severely influencing the final yield. In the present work, the heat stress response mechanisms of one thermotolerant genotype (E42) were investigated by exploring its regulatory gene network. This was achieved through a promoter analysis based on the identification of the heat stress elements (HSEs) mapping in the promoters, combined with a gene co-expression network analysis aimed at identifying interactions among heat-related genes.

RESULTS

Results highlighted 82 genes presenting HSEs in the promoter and belonging to one of the 52 gene networks obtained by the GCN analysis; 61 of these also interact with heat shock factors (Hsfs). Finally, a list of 13 candidate genes including two Hsfs, nine heat shock proteins (Hsps) and two GDSL esterase/lipase (GELPs) were retrieved by focusing on those E42 genes exhibiting HSEs in the promoters, interacting with Hsfs and showing variants, compared to Heinz reference genome, with HIGH and/or MODERATE impact on the translated protein. Among these, the Gene Ontology annotation analysis evidenced that only LeHsp100 (Solyc02g088610) belongs to a network specifically involved in the response to heat stress.

CONCLUSIONS

As a whole, the combination of bioinformatic analyses carried out on genomic and trascriptomic data available for tomato, together with polymorphisms detected in HS-related genes of the thermotolerant E42 allowed to determine a subset of candidate genes involved in the HS response in tomato. This study provides a novel approach in the investigation of abiotic stress response mechanisms and further studies will be conducted to validate the role of the highlighted genes.

Molecular regulation of PPARγ/RXRα signaling by the novel cofactor ZFP407

Cofactors interacting with PPARγ can regulate adipogenesis and adipocyte metabolism by modulating the transcriptional activity and selectivity of PPARγ signaling. ZFP407 was previously demonstrated to regulate PPARγ target genes such as GLUT4, and its overexpression improved glucose homeostasis in mice. Here, using a series of molecular assays, including protein-interaction studies, mutagenesis, and ChIP-seq, ZFP407 was found to interact with the PPARγ/RXRα protein complex in the nucleus of adipocytes. Consistent with this observation, ZFP407 ChIP-seq peaks significantly overlapped with PPARγ ChIP-seq peaks, with more than half of ZFP407 peaks overlapping with PPARγ peaks. Transcription factor binding motifs enriched in these overlapping sites included CTCF, RARα/RXRγ, TP73, and ELK1, which regulate cellular development and function within adipocytes. Site-directed mutagenesis of frequent PPARγ phosphorylation or SUMOylation sites did not prevent its regulation by ZFP407, while mutagenesis of ZFP407 domains potentially necessary for RXR and PPARγ binding abrogated any impact of ZFP407 on PPARγ activity. These data suggest that ZFP407 controls the activity of PPARγ, but does so independently of post-translational modifications, likely by direct binding, establishing ZFP407 as a newly identified PPARγ cofactor. In addition, ZFP407 ChIP-seq analyses identified regions that did not overlap with PPARγ peaks. These non-overlapping peaks were significantly enriched for the transcription factor binding motifs of TBX19, PAX8, HSF4, and ZKSCAN3, which may contribute to the PPARγ-independent functions of ZFP407 in adipocytes and other cell types.

Regional-specific calibration enables application of computational evidence for clinical classification of 5' cis-regulatory variants in Mendelian disease

To date, clinical genetic testing for Mendelian disease variants has focused heavily on exonic coding and intronic gene regions. This multi-step study was undertaken to provide an evidence base for selecting and applying computational approaches for use in clinical classification of 5' cis-regulatory region variants. Curated datasets of clinically reported disease-causing 5' cis-regulatory region variants and variants from matched genomic regions in population controls were used to calibrate six bioinformatic tools as predictors of variant pathogenicity. Likelihood ratio estimates were aligned to code weights following ClinGen recommendations for application of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) classification scheme. Considering code assignment across all reference dataset variants, performance was best for CADD (81.2%) and REMM (81.5%). Optimized thresholds provided moderate evidence toward pathogenicity (CADD, REMM) and moderate (CADD) or supporting (REMM) evidence against pathogenicity. Both sensitivity and specificity of prediction were improved when further categorizing variants based on location in an EPDnew-defined promoter region. Combining predictions (CADD, REMM, and location in a promoter region) increased specificity at the expense of sensitivity. Importantly, the optimal CADD thresholds for assigning ACMG/AMP codes PP3 (≥10) and BP4 (≤8) were vastly different from recommendations for protein-coding variants (PP3 ≥25.3; BP4 ≤22.7); CADD <22.7 would incorrectly assign BP4 for >90% of reported disease-causing cis-regulatory region variants. Our results demonstrate the need to consider a tiered approach and tailored score thresholds to optimize bioinformatic impact prediction for clinical classification of 5' cis-regulatory region variants.

ATRX guards against aberrant differentiation in mesenchymal progenitor cells

Alterations in the tumor suppressor ATRX are recurrently observed in mesenchymal neoplasms. ATRX has multiple epigenetic functions including heterochromatin formation and maintenance and regulation of transcription through modulation of chromatin accessibility. Here, we show in murine mesenchymal progenitor cells (MPCs) that Atrx deficiency aberrantly activated mesenchymal differentiation programs. This includes adipogenic pathways where ATRX loss induced expression of adipogenic transcription factors and enhanced adipogenic differentiation in response to differentiation stimuli. These changes are linked to loss of heterochromatin near mesenchymal lineage genes together with increased chromatin accessibility and gains of active chromatin marks. We additionally observed depletion of H3K9me3 at transposable elements, which are derepressed including near mesenchymal genes where they could serve as regulatory elements. Finally, we demonstrated that loss of ATRX in a mesenchymal malignancy, undifferentiated pleomorphic sarcoma, results in similar epigenetic disruption and de-repression of transposable elements. Together, our results reveal a role for ATRX in maintaining epigenetic states and transcriptional repression in mesenchymal progenitors and tumor cells and in preventing aberrant differentiation in the progenitor context.

RNA-RNA interactions between respiratory syncytial virus and miR-26 and miR-27 are associated with regulation of cell cycle and antiviral immunity

microRNAs (miRNAs) regulate nearly all physiological processes but our understanding of exactly how they function remains incomplete, particularly in the context of viral infections. Here, we adapt a biochemical method (CLEAR-CLIP) and analysis pipeline to identify targets of miRNAs in lung cells infected with Respiratory syncytial virus (RSV). We show that RSV binds directly to miR-26 and miR-27 through seed pairing and demonstrate that these miRNAs target distinct gene networks associated with cell cycle and metabolism (miR-27) and antiviral immunity (miR-26). Many of the targets are de-repressed upon infection and we show that the miR-27 targets most sensitive to miRNA inhibition are those associated with cell cycle. Finally, we demonstrate that high confidence chimeras map to long noncoding RNAs (lncRNAs) and pseudogenes in transcriptional regulatory regions. We validate that a proportion of miR-27 and Argonaute 2 (AGO2) is nuclear and identify a long non-coding RNA (lncRNA) as a miR-27 target that is linked to transcriptional regulation of nearby genes. This work expands the target networks of miR-26 and miR-27 to include direct interactions with RSV and lncRNAs and implicate these miRNAs in regulation of key genes that impact the viral life cycle associated with cell cycle, metabolism, and antiviral immunity.

The tidyomics ecosystem: Enhancing omic data analyses

The growth of omic data presents evolving challenges in data manipulation, analysis, and integration. Addressing these challenges, Bioconductor1 provides an extensive community-driven biological data analysis platform. Meanwhile, tidy R programming2 offers a revolutionary standard for data organisation and manipulation. Here, we present the tidyomics software ecosystem, bridging Bioconductor to the tidy R paradigm. This ecosystem aims to streamline omic analysis, ease learning, and encourage cross-disciplinary collaborations. We demonstrate the effectiveness of tidyomics by analysing 7.5 million peripheral blood mononuclear cells from the Human Cell Atlas3, spanning six data frameworks and ten analysis tools.

m6ACali: machine learning-powered calibration for accurate m6A detection in MeRIP-Seq

We present m6ACali, a novel machine-learning framework aimed at enhancing the accuracy of N6-methyladenosine (m6A) epitranscriptome profiling by reducing the impact of non-specific antibody enrichment in MeRIP-Seq. The calibration model serves as a genomic feature-based classifier that refines the identification of m6A sites, distinguishing those genuinely present from those that can be detected in in-vitro transcribed (IVT) control experiments. We find that m6ACali effectively identifies non-specific binding peaks reported by exomePeak2 and MACS2 in novel MeRIP-Seq datasets without the need for paired IVT controls. The model interpretation revealed that off-target antibody binding sites commonly occur at short exons and short mRNAs, originating from high read coverage regions that share the motif sequence with true m6A sites. We also reveal that the ML strategy can efficiently adjust differentially methylated peaks and other antibody-dependent, base-resolution m6A detection techniques. As a result, m6ACali offers a promising method for the universal enhancement of m6A profiles generated by MeRIP-Seq experiments, elevating the benchmark for omics-level m6A data integration.

Comprehensive genomic features indicative for Notch responsiveness

Transcription factor RBPJ is the central component in Notch signal transduction and directly forms a coactivator complex together with the Notch intracellular domain (NICD). While RBPJ protein levels remain constant in most tissues, dynamic expression of Notch target genes varies depending on the given cell-type and the Notch activity state. To elucidate dynamic RBPJ binding genome-wide, we investigated RBPJ occupancy by ChIP-Seq. Surprisingly, only a small set of the total RBPJ sites show a dynamic binding behavior in response to Notch signaling. Compared to static RBPJ sites, dynamic sites differ in regard to their chromatin state, binding strength and enhancer positioning. Dynamic RBPJ sites are predominantly located distal to transcriptional start sites (TSSs), while most static sites are found in promoter-proximal regions. Importantly, gene responsiveness is preferentially associated with dynamic RBPJ binding sites and this static and dynamic binding behavior is repeatedly observed across different cell types and species. Based on the above findings we used a machine-learning algorithm to predict Notch responsiveness with high confidence in different cellular contexts. Our results strongly support the notion that the combination of binding strength and enhancer positioning are indicative of Notch responsiveness.

Reciprocal conversion between annual and polycarpic perennial flowering behavior in the Brassicaceae

The development of perennial crops holds great promise for sustainable agriculture and food security. However, the evolution of the transition between perenniality and annuality is poorly understood. Here, using two Brassicaceae species, Crucihimalaya himalaica and Erysimum nevadense, as polycarpic perennial models, we reveal that the transition from polycarpic perennial to biennial and annual flowering behavior is a continuum determined by the dosage of three closely related MADS-box genes. Diversification of the expression patterns, functional strengths, and combinations of these genes endows species with the potential to adopt various life-history strategies. Remarkably, we find that a single gene among these three is sufficient to convert winter-annual or annual Brassicaceae plants into polycarpic perennial flowering plants. Our work delineates a genetic basis for the evolution of diverse life-history strategies in plants and lays the groundwork for the generation of diverse perennial Brassicaceae crops in the future.

Determinants of gastric cancer immune escape identified from non-coding immune-landscape quantitative trait loci

The landscape of non-coding mutations in cancer progression and immune evasion is largely unexplored. Here, we identify transcrptome-wide somatic and germline 3' untranslated region (3'-UTR) variants from 375 gastric cancer patients from The Cancer Genome Atlas. By performing gene expression quantitative trait loci (eQTL) and immune landscape QTL (ilQTL) analysis, we discover 3'-UTR variants with cis effects on expression and immune landscape phenotypes, such as immune cell infiltration and T cell receptor diversity. Using a massively parallel reporter assay, we distinguish between causal and correlative effects of 3'-UTR eQTLs in immune-related genes. Our approach identifies numerous 3'-UTR eQTLs and ilQTLs, providing a unique resource for the identification of immunotherapeutic targets and biomarkers. A prioritized ilQTL variant signature predicts response to immunotherapy better than standard-of-care PD-L1 expression in independent patient cohorts, showcasing the untapped potential of non-coding mutations in cancer.

Biologically informed NeuralODEs for genome-wide regulatory dynamics

BACKGROUND

Gene regulatory network (GRN) models that are formulated as ordinary differential equations (ODEs) can accurately explain temporal gene expression patterns and promise to yield new insights into important cellular processes, disease progression, and intervention design. Learning such gene regulatory ODEs is challenging, since we want to predict the evolution of gene expression in a way that accurately encodes the underlying GRN governing the dynamics and the nonlinear functional relationships between genes. Most widely used ODE estimation methods either impose too many parametric restrictions or are not guided by meaningful biological insights, both of which impede either scalability, explainability, or both.

RESULTS

We developed PHOENIX, a modeling framework based on neural ordinary differential equations (NeuralODEs) and Hill-Langmuir kinetics, that overcomes limitations of other methods by flexibly incorporating prior domain knowledge and biological constraints to promote sparse, biologically interpretable representations of GRN ODEs. We tested the accuracy of PHOENIX in a series of in silico experiments, benchmarking it against several currently used tools. We demonstrated PHOENIX's flexibility by modeling regulation of oscillating expression profiles obtained from synchronized yeast cells. We also assessed the scalability of PHOENIX by modeling genome-scale GRNs for breast cancer samples ordered in pseudotime and for B cells treated with Rituximab.

CONCLUSIONS

PHOENIX uses a combination of user-defined prior knowledge and functional forms from systems biology to encode biological "first principles" as soft constraints on the GRN allowing us to predict subsequent gene expression patterns in a biologically explainable manner.

New Drosophila promoter-associated architectural protein Mzfp1 interacts with CP190 and is required for housekeeping gene expression and insulator activity

In Drosophila, a group of zinc finger architectural proteins recruits the CP190 protein to the chromatin, an interaction that is essential for the functional activity of promoters and insulators. In this study, we describe a new architectural C2H2 protein called Madf and Zinc-Finger Protein 1 (Mzfp1) that interacts with CP190. Mzfp1 has an unusual structure that includes six C2H2 domains organized in a C-terminal cluster and two tandem MADF domains. Mzfp1 predominantly binds to housekeeping gene promoters located in both euchromatin and heterochromatin genome regions. In vivo mutagenesis studies showed that Mzfp1 is an essential protein, and both MADF domains and the CP190 interaction region are required for its functional activity. The C2H2 cluster is sufficient for the specific binding of Mzfp1 to regulatory elements, while the second MADF domain is required for Mzfp1 recruitment to heterochromatin. Mzfp1 binds to the proximal part of the Fub boundary that separates regulatory domains of the Ubx and abd-A genes in the Bithorax complex. Mzfp1 participates in Fub functions in cooperation with the architectural proteins Pita and Su(Hw). Thus, Mzfp1 is a new architectural C2H2 protein involved in the organization of active promoters and insulators in Drosophila.

Genome-wide association study of copy number variations with shank traits in a F2 crossbred chicken population

Copy number variations (CNVs) are large-scale changes in the DNA sequence that can affect the genetic structure and phenotype of an organism. The purpose of this study was to investigate the existing CNVs and their associations with the shank diameter (ShD) and shank length (ShL) traits using data from an F2 crossbred chicken population. To carry out the study, 312 chickens were genotyped using the Illumina 60k SNP Beadchip. The shank traits of the birds were measured from day 1 to 12 weeks of age. penncnv and cnvruler tools were used to find copy numbers and regions with copy number changes (CNVR), respectively. The CNVRanger package was used to perform a genome-wide association study between shank traits and CNVs. Gene ontology research in CNVRs was carried out using the david database. In this investigation, 966 CNVs and 606 regions with copy number changes were discovered. The copy number states and variations were randomly distributed along the length of the autosomal chromosomes. Weeks 1-4, 9 and 12 of growth revealed a significant association of copy number variations with shank traits, false discovery rate (FDR-corrected p-value < 0.01), and the majority of CNVs that were statistically significant were found on chromosomes 1-3. These CNV segments are nearby genes such as KCNJ12, FGF6 and MYF5, which are fundamental to growth and development. In addition, gene set analyses revealed terms related to muscle physiology, regulation of cellular processes and potassium channels.