We propose the Model-based Analysis of Genome-wide CRISPR/Cas9 Knockout (MAGeCK) method for prioritizing single-guide RNAs, genes and pathways in genome-scale CRISPR/Cas9 knockout screens. MAGeCK demonstrates better performance compared with existing methods, identifies both positively and negatively selected genes simultaneously, and reports robust results across different experimental conditions. Using public datasets, MAGeCK identified novel essential genes and pathways, including EGFR in vemurafenib-treated A375 cells harboring a BRAF mutation. MAGeCK also detected cell type-specific essential genes, including BCR and ABL1, in KBM7 cells bearing a BCR-ABL fusion, and IGF1R in HL-60 cells, which depends on the insulin signaling pathway for proliferation.

Water chickweed (Myosoton aquaticum (L.)) is a dicot broadleaf weed that is widespread in winter fields in China, and has evolved serious resistance to acetolactate synthase (ALS) inhibiting herbicides.

Patients with rheumatoid arthritis (RA) may be more susceptible to infection by coronavirus disease-19 (COVID-19) due to immune system dysfunction. However, there are still insufficient treatment strategies for patients with RA and COVID-19. Since Jingulian is a traditional Chinese medicine (TCM) with anti-viral and immune regulatory functions, our study aims to explore the detailed mechanisms of Jingulian in treating patients with RA and COVID-19.

WNT5a is a mainly "non-canonical" WNT ligand whose dysregulation is observed in lung diseases such as idiopathic pulmonary fibrosis (IPF), chronic obstructive pulmonary disease (COPD) and asthma. Germline deletion of Wnt5a disrupts embryonic lung development. However, the temporal-specific function of WNT5a remains unknown. In this study, we generated a conditional loss-of-function mouse model (Wnt5aCAG) and examined the specific role of Wnt5a during the saccular and alveolar phases of lung development. The lack of Wnt5a in the saccular phase blocked distal airway expansion and attenuated differentiation of endothelial and alveolar epithelial type I (AT1) cells and myofibroblasts. Postnatal Wnt5a inactivation disrupted alveologenesis, producing a phenotype resembling human bronchopulmonary dysplasia (BPD). Mutant lungs showed hypoalveolization, but endothelial and epithelial differentiation was unaffected. The major impact of Wnt5a inactivation on alveologenesis was on myofibroblast differentiation and migration, with reduced expression of key regulatory genes. These findings were validated in vitro using isolated lung fibroblasts. Conditional inactivation of the WNT5a receptors Ror1 and Ror2 in alveolar myofibroblasts recapitulated the Wnt5aCAG phenotype, demonstrating that myofibroblast defects are the major cause of arrested alveologenesis in Wnt5aCAG lungs. Finally, we show that WNT5a is reduced in human BPD lung samples, indicating the clinical relevance and potential role for WNT5a in pathogenesis of BPD.

Although millions of transcription factor binding sites, or cistromes, have been identified across the human genome, defining which of these sites is functional in a given condition remains challenging. Using CRISPR/Cas9 knockout screens and gene essentiality or fitness as the readout, we systematically investigated the essentiality of over 10,000 FOXA1 and CTCF binding sites in breast and prostate cancer cells. We found that essential FOXA1 binding sites act as enhancers to orchestrate the expression of nearby essential genes through the binding of lineage-specific transcription factors. In contrast, CRISPR screens of the CTCF cistrome revealed 2 classes of essential binding sites. The first class of essential CTCF binding sites act like FOXA1 sites as enhancers to regulate the expression of nearby essential genes, while a second class of essential CTCF binding sites was identified at topologically associated domain (TAD) boundaries and display distinct characteristics. Using regression methods trained on our screening data and public epigenetic profiles, we developed a model to predict essential cis-elements with high accuracy. The model for FOXA1 essentiality correctly predicts noncoding variants associated with cancer risk and progression. Taken together, CRISPR screens of cis-regulatory elements can define the essential cistrome of a given factor and can inform the development of predictive models of cistrome function.

Trimethylated histone H3 lysine 27 (H3K27me3) is linked to gene silencing, whereas H3K4me3 is associated with gene activation. These two marks frequently co-occupy gene promoters, forming bivalent domains. Bivalency signifies repressed but activatable states of gene expression and can be resolved to active, H3K4me3-prevalent states during multiple cellular processes, including differentiation, development and epithelial mesenchymal transition. However, the molecular mechanism underlying bivalency resolution remains largely unknown. Here, we show that the H3K27 demethylase UTX (also called KDM6A) is required for the resolution and activation of numerous retinoic acid (RA)-inducible bivalent genes during the RA-driven differentiation of mouse embryonic stem cells (ESCs). Notably, UTX loss in mouse ESCs inhibited the RA-driven bivalency resolution and activation of most developmentally critical homeobox (Hox) a-d genes. The UTX-mediated resolution and activation of many bivalent Hox genes during mouse ESC differentiation were recapitulated during RA-driven differentiation of human NT2/D1 embryonal carcinoma cells. In support of the importance of UTX in bivalency resolution, Utx-null mouse ESCs and UTX-depleted NT2/D1 cells displayed defects in RA-driven cellular differentiation. Our results define UTX as a bivalency-resolving histone modifier necessary for stem cell differentiation.

The four-subunit negative elongation factor (NELF) complex mediates RNA polymerase II (Pol II) pausing at promoter-proximal regions. Ablation of individual NELF subunits destabilizes the NELF complex and causes cell lethality, leading to the prevailing concept that NELF-mediated Pol II pausing is essential for cell proliferation. Using separation-of-function mutations, we show here that NELFB function in cell proliferation can be uncoupled from that in Pol II pausing. NELFB mutants sequestered in the cytoplasm and deprived of NELF nuclear function still support cell proliferation and part of the NELFB-dependent transcriptome. Mechanistically, cytoplasmic NELFB physically and functionally interacts with prosurvival signaling kinases, most notably phosphatidylinositol-3-kinase/AKT. Ectopic expression of membrane-tethered phosphatidylinositol-3-kinase/AKT partially bypasses the role of NELFB in cell proliferation, but not Pol II occupancy. Together, these data expand the current understanding of the physiological impact of Pol II pausing and underscore the multiplicity of the biological functions of individual NELF subunits.

Alveolar type II (ATII) cells produce pulmonary surfactant (PS) essential for maintaining lung function. The aberration or depletion of PS can cause alveolar collapse, a hallmark of acute respiratory distress syndrome (ARDS). However, the intricacies underlying these changes remain unclear. This study aimed to elucidate the mechanisms underlying PS perturbations in ATII cells using transcriptional RNA-seq, offering insights into the pathogenesis of ARDS.

Intellectual disability (ID) is one of the most prevalent chronic developmental brain disorders or phenotype of syndromic ID, affecting nearly 1-2% of the general population worldwide. Over recent decades, tremendous effort and high-throughput platforms have been devised to explore the complex heterogeneity, numerous genes and variants have been associated with the ID, especially de novo mutations and copy number variants. An organized resource containing the increasing genetic data is imperative to assist ID research. In this study, the integrative and annotated intellectual disability database has been developed, named 'IDGenetics', which contains known information about ID, including 815 genes and 17102 variants associated with 918 clinical diseases (3001 clinical phenotype) collected from 3822 publications and ID-related databases. Furthermore, in-depth data mining was performed to obtain an understanding of each entry, including functional annotation, gene/disease/phenotype network establishment and overlap analysis focusing on comorbidity. 1478 candidate genes (483 high-confidence and 995 low-confidence) were collected and prioritized by adopting the annotations of 12 functional prediction tools and algorithm. In addition, IDGenetics database provides concise search methods, convenient browsing functions, intuitive graphical displays and constantly updated features. IDGenetics will be a valuable and integrative resource for deciphering the genetic and functional architecture of ID and the improvement of clinical diagnosis, intervention and treatment.

Alternative polyadenylation (APA) plays an essential role in brain development; however, current transcriptome-wide association studies (TWAS) largely overlook APA in nominating susceptibility genes. Here, we performed a 3' untranslated region (3'UTR) APA TWAS (3'aTWAS) for 11 brain disorders by combining their genome-wide association studies data with 17,300 RNA-seq samples across 2,937 individuals. We identified 354 3'aTWAS-significant genes, including known APA-linked risk genes, such as SNCA in Parkinson's disease. Among these 354 genes, ~57% are not significant in traditional expression- and splicing-TWAS studies, since APA may regulate the translation, localization and protein-protein interaction of the target genes independent of mRNA level expression or splicing. Furthermore, we discovered ATXN3 as a 3'aTWAS-significant gene for amyotrophic lateral sclerosis, and its modulation substantially impacted pathological hallmarks of amyotrophic lateral sclerosis in vitro. Together, 3'aTWAS is a powerful strategy to nominate important APA-linked brain disorder susceptibility genes, most of which are largely overlooked by conventional expression and splicing analyses.

Hepatocellular carcinoma (HCC) is one of the most common heterogeneous tumors that occurs after chronic liver diseases and hepatitis virus infection. Immune-related genes (IRGs) and their ligands regulate the homeostasis of tumor microenvironment, which is essential for the treatment of HCC and its prognosis. This study aimed to investigate the clinical value of IRGs in predicting the prognosis of HCC.

Parkinson's disease (PD) is characterized by α-synuclein aggregation in dopaminergic (DA) neurons, which are sensitive to oxidative stress. Mitochondria aconitase 2 (ACO2) is an essential enzyme in the tricarboxylic acid cycle that orchestrates mitochondrial and autophagic functions to energy metabolism. Though widely linked to diseases, its relation to PD has not been fully clarified. Here we revealed that the peripheral ACO2 activity was significantly decreased in PD patients and associated with their onset age and disease durations. The knock-in mouse and Drosophila models with the A252T variant displayed aggravated motor deficits and DA neuron degeneration after 6-OHDA and rotenone-induction, and the ACO2 knockdown or blockade cells showed features of mitochondrial and autophagic dysfunction. Moreover, the transcription of autophagy-related genes LC3 and Atg5 was significantly downregulated via inhibited histone acetylation at the H3K9 and H4K5 sites. These data provided multi-dimensional evidences supporting the essential roles of ACO2, and as a potential early biomarker to be used in clinical trials for assessing the effects of antioxidants in PD. Moreover, ameliorating energy metabolism by targeting ACO2 could be considered as a potential therapeutic strategy for PD and other neurodegenerative disorders.

Calcium and integrin-binding protein 2 (CIB2) belongs to a protein family with four known members, CIB1 through CIB4, which are characterized by multiple calcium-binding EF-hand domains. Among the family members, the Cib1 and Cib2 genes are expressed in mouse cochlear hair cells, and mutations in the human CIB2 gene have been associated with nonsyndromic deafness DFNB48 and syndromic deafness USH1J. To further explore the function of CIB1 and CIB2 in hearing, we established Cib1 and Cib2 knockout mice using the clustered regularly interspaced short palindromic repeat (CRISPR)-associated Cas9 nuclease (CRISPR/Cas9) genome editing technique. We found that loss of CIB1 protein does not affect auditory function, whereas loss of CIB2 protein causes profound hearing loss in mice. Further investigation revealed that hair cell stereocilia development is affected in Cib2 knockout mice. Noticeably, loss of CIB2 abolishes mechanoelectrical transduction (MET) currents in auditory hair cells. In conclusion, we show here that although both CIB1 and CIB2 are readily detected in the cochlea, only loss of CIB2 results in profound hearing loss, and that CIB2 is essential for auditory hair cell MET.

Characterizing Mycobacterium abscessus complex (MABSC) biofilms under host-relevant conditions is essential to the design of informed therapeutic strategies targeted to this persistent, drug-tolerant, population of extracellular bacilli. Using synthetic cystic fibrosis medium (SCFM) which we previously reported to closely mimic the conditions encountered by MABSC in actual cystic fibrosis (CF) sputum and a new model of biofilm formation, we show that MABSC biofilms formed under these conditions are substantially different from previously reported biofilms grown in standard laboratory media in terms of their composition, gene expression profile and stress response. Extracellular DNA (eDNA), mannose-and glucose-containing glycans and phospholipids, rather than proteins and mycolic acids, were revealed as key extracellular matrix (ECM) constituents holding clusters of bacilli together. None of the environmental cues previously reported to impact biofilm development had any significant effect on SCFM-grown biofilms, most likely reflecting the fact that SCFM is a nutrient-rich environment in which MABSC finds a variety of ways of coping with stresses. Finally, molecular determinants were identified that may represent attractive new targets for the development of adjunct therapeutics targeting MABSC biofilms in persons with CF.

How tumor suppressor p53 selectively responds to specific signals, especially in normal cells, is poorly understood. We performed genome-wide profiling of p53 chromatin interactions and target gene expression in human embryonic stem cells (hESCs) in response to early differentiation, induced by retinoic acid, versus DNA damage, caused by adriamycin. Most p53-binding sites are unique to each state and define stimulus-specific p53 responses in hESCs. Differentiation-activated p53 targets include many developmental transcription factors and, in pluripotent hESCs, are bound by OCT4 and NANOG at chromatin enriched in both H3K27me3 and H3K4me3. Activation of these genes occurs with recruitment of p53 and H3K27me3-specific demethylases, UTX and JMJD3, to chromatin. In contrast, genes associated with cell migration and motility are bound by p53 specifically after DNA damage. Surveillance functions of p53 in cell death and cell cycle regulation are conserved in both DNA damage and differentiation. Comparative genomic analysis of p53-targets in mouse and human ESCs supports an inter-species divergence in p53 regulatory functions during evolution. Our findings expand the registry of p53-regulated genes to define p53-regulated opposition to pluripotency during early differentiation, a process highly distinct from stress-induced p53 response in hESCs.

Embryonic stem cells are maintained in a self-renewing and pluripotent state by multiple regulatory pathways. Pluripotent-specific transcriptional networks are sequentially reactivated as somatic cells reprogram to achieve pluripotency. How epigenetic regulators modulate this process and contribute to somatic cell reprogramming is not clear. Here we performed a functional RNAi screen to identify the earliest epigenetic regulators required for reprogramming. We identified components of the SAGA histone acetyltransferase complex, in particular Gcn5, as critical regulators of reprogramming initiation. Furthermore, we showed in mouse pluripotent stem cells that Gcn5 strongly associates with Myc and that, upon initiation of somatic reprogramming, Gcn5 and Myc form a positive feed-forward loop that activates a distinct alternative splicing network and the early acquisition of pluripotency-associated splicing events. These studies expose a Myc-SAGA pathway that drives expression of an essential alternative splicing regulatory network during somatic cell reprogramming.

Histone acetylation is associated with active transcription in eukaryotic cells. It helps to open up the chromatin by neutralizing the positive charge of histone lysine residues and providing binding platforms for "reader" proteins. The bromodomain (BRD) has long been thought to be the sole protein module that recognizes acetylated histones. Recently, we identified the YEATS domain of AF9 (ALL1 fused gene from chromosome 9) as a novel acetyl-lysine-binding module and showed that the ENL (eleven-nineteen leukemia) YEATS domain is an essential acetyl-histone reader in acute myeloid leukemias. The human genome encodes four YEATS domain proteins, including GAS41, a component of chromatin remodelers responsible for H2A.Z deposition onto chromatin; however, the importance of the GAS41 YEATS domain in human cancer remains largely unknown. Here we report that GAS41 is frequently amplified in human non-small cell lung cancer (NSCLC) and is required for cancer cell proliferation, survival, and transformation. Biochemical and crystal structural studies demonstrate that GAS41 binds to histone H3 acetylated on H3K27 and H3K14, a specificity that is distinct from that of AF9 or ENL. ChIP-seq (chromatin immunoprecipitation [ChIP] followed by high-throughput sequencing) analyses in lung cancer cells reveal that GAS41 colocalizes with H3K27ac and H3K14ac on the promoters of actively transcribed genes. Depletion of GAS41 or disruption of the interaction between its YEATS domain and acetylated histones impairs the association of histone variant H2A.Z with chromatin and consequently suppresses cancer cell growth and survival both in vitro and in vivo. Overall, our study identifies GAS41 as a histone acetylation reader that promotes histone H2A.Z deposition in NSCLC.

N 6-methyladenosine (m6A) is one of the common forms of RNA modifications. METTL3 is the essential factor that has methyltransferase activity. One important role of m6A is to regulate translation of mRNAs via reader YTHDF1. However, whether m6A could indirectly affect the translation of mRNA through other approaches remains unclear.

Approximately ten percent of male infertility is caused by non-obstructive azoospermia (NOA), but the etiologies of many NOA remain elusive. Recently, a genome-wide association study (GWAS) of NOA in Han Chinese men was conducted, and only a few genetic variants associated with NOA were found, which might have resulted from genetic heterogeneity. However, those variants that lack genome-wide significance might still be essential for fertility. Functional analysis of genes surrounding these variants in Drosophila identified some spermatogenesis-essential genes. As a complementary method of Drosophila screening, SK1 background Saccharomvces cerevisiae was used as a model to screen meiosis-related genes from the NOA GWAS data in this study. After functional screening, GDA1 (orthologous to humanENTPD6) was found to be a novel meiosis-related gene. The deletion of GDA1 resulted in the failure of yeast sporulation. Further investigations showed that Gda1p was important for pre-meiotic S phase entry. Interestingly, the meiotic role of Gda1p was dependent on its guanosine diphosphatase activity, but not it's cytoplasmic, transmembrane or stem domains. These yeast data suggest that ENTPD6 may be a novel meiosis-associated NOA-related gene, and the yeast model provides a good approach to analyze GWAS results of NOA.

Rett syndrome (RTT) is a neurodevelopmental disorder caused by mutations in methyl-CpG-binding protein-2 (MECP2), a transcriptional regulator of many genes, including brain-derived neurotrophic factor (BDNF). BDNF levels are reduced in RTT autopsy brains and in multiple brain areas of Mecp2-deficient mice. Furthermore, experimental interventions that increase BDNF levels improve RTT-like phenotypes in Mecp2 mutant mice. Here, we characterized the actions of a small-molecule ligand of the BDNF receptor TrkB in hippocampal function in Mecp2 mutant mice. Systemic treatment of female Mecp2 heterozygous (HET) mice with LM22A-4 for 4 weeks improved hippocampal-dependent object location memory and restored hippocampal long-term potentiation (LTP). Mechanistically, LM22A-4 acts to dampen hyperactive hippocampal network activity, reduce the frequency and amplitude of miniature excitatory postsynaptic currents (mEPSCs), and reduce the frequency of spontaneous tetrodotoxin-resistant Ca2+ signals in Mecp2 mutant hippocampal neurons, making them comparable to those features observed in wild-type neurons. Together, these observations indicate that LM22A-4 is a promising therapeutic candidate for the treatment of hippocampal dysfunction in RTT.