Knowledge of the structure of gene expression is essential for mammalian transcriptomics research. We analyzed a collection of more than one million porcine expressed sequence tags (ESTs), of which two-thirds were generated in the Sino-Danish Pig Genome Project and one-third are from public databases. The Sino-Danish ESTs were generated from one normalized and 97 non-normalized cDNA libraries representing 35 different tissues and three developmental stages.

OX40 and OX40 ligand (OX40L), as essential immune checkpoint (IC) modulators, are highly correlated with cancer immunity regulation as well as tumor microenvironment (TME). Immunotherapy showed outstanding advantages in small-cell lung cancer (SCLC) therapy. However, functions and clinical significance of OX40 and OX40L in SCLC were not clear yet.

Knowledge of immune cell phenotypes, function, and developmental trajectory in acute myeloid leukemia (AML) microenvironment is essential for understanding mechanisms of evading immune surveillance and immunotherapy response of targeting special microenvironment components.

Magnesium-protoporphyrin IX monomethyl ester cyclase (MPEC) catalyzes the conversion of MPME to divinyl protochlorophyllide (DVpchlide). This is an essential enzyme during chlorophyll (Chl) biosynthesis but details of its function in rice are still lacking. Here, we identified a novel rice mutant yellow-leaf 1 (yl-1), which showed decreased Chl accumulation, abnormal chloroplast ultrastructure and attenuated photosynthetic activity. Map-based cloning and over-expression analysis suggested that YL-1 encodes a subunit of MPEC. The YL-1 protein localizes in chloroplasts, and it is mainly expressed in green tissues, with greatest abundance in leaves and young panicles. Results of qRT-PCR showed that Chl biosynthesis upstream genes were highly expressed in the yl-1 mutant, while downstream genes were compromised, indicating that YL-1 plays a pivotal role in the Chl biosynthesis. Furthermore, the expression levels of photosynthesis and chloroplast development genes were also affected. RNA-seq results futher proved that numerous membrane-associated genes, including many plastid membrane-associated genes, have altered expression pattern in the yl-1 mutant, implying that YL-1 is required for plastid membrane stability. Thus, our study confirms a putative MPME cyclase as a novel key enzyme essential for Chl biosynthesis and chloroplast membrane stability in rice.

One of the challenging problems in the etiology of diseases is to explore the relationships between initiation and progression of diseases and abnormalities in local regions of metabolic pathways. To gain insight into such relationships, we applied the "k-clique" subpathway identification method to all disease-related gene sets. For each disease, the disease risk regions of metabolic pathways were then identified and considered as subpathways associated with the disease. We finally built a disease-metabolic subpathway network (DMSPN). Through analyses based on network biology, we found that a few subpathways, such as that of cytochrome P450, were highly connected with many diseases, and most belonged to fundamental metabolisms, suggesting that abnormalities of fundamental metabolic processes tend to cause more types of diseases. According to the categories of diseases and subpathways, we tested the clustering phenomenon of diseases and metabolic subpathways in the DMSPN. The results showed that both disease nodes and subpathway nodes displayed slight clustering phenomenon. We also tested correlations between network topology and genes within disease-related metabolic subpathways, and found that within a disease-related subpathway in the DMSPN, the ratio of disease genes and the ratio of tissue-specific genes significantly increased as the number of diseases caused by the subpathway increased. Surprisingly, the ratio of essential genes significantly decreased and the ratio of housekeeping genes remained relatively unchanged. Furthermore, the coexpression levels between disease genes and other types of genes were calculated for each subpathway in the DMSPN. The results indicated that those genes intensely influenced by disease genes, including essential genes and tissue-specific genes, might be significantly associated with the disease diversity of subpathways, suggesting that different kinds of genes within a disease-related subpathway may play significantly differential roles on the diversity of diseases caused by the corresponding subpathway.

Protein phosphatase 2A (PP2A) plays an important role in the control of the cell cycle. We previously reported that the PP2A inhibitors, cantharidin and okadaic acid (OA), efficiently repressed the growth of cancer cells. In the present study, we found that PP2A inhibitors arrested the cell cycle at the G2 phase through a mechanism that was dependent on the JNK pathway. Microarrays further showed that PP2A inhibitors induced expression changes in multiple genes that participate in cell cycle transition. To verify whether these expression changes were executed in a PP2A-dependent manner, we targeted the PP2A catalytic subunit (PP2Ac) using siRNA and evaluated gene expression with a microarray. After the cross comparison of these microarray data, we identified that CDK1 was potentially the same target when treated with either PP2A inhibitors or PP2Ac siRNA. In addition, we found that the down-regulation of CDK1 occurred in a JNK-dependent manner. Luciferase reporter gene assays demonstrated that repression of the transcription of CDK1 was executed through the JNK-dependent activation of the Sp1 transcription factor. By constructing deletion mutants of the CDK1 promoter and by using ChIP assays, we identified an element in the CDK1 promoter that responded to the JNK/Sp1 pathway after stimulation with PP2A inhibitors. Cantharidin and OA also up-regulated the expression of p21, an inhibitor of CDK1, via autophagy rather than PP2A/JNK pathway. Thus, this present study found that the PP2A/JNK/Sp1/CDK1 pathway and the autophagy/p21 pathway participated in G2/M cell cycle arrest triggered by PP2A inhibitors.

Follicle rupture, the final step in ovulation, utilizes conserved molecular mechanisms including matrix metalloproteinases (Mmps), steroid signaling, and adrenergic signaling. It is still unknown how follicles become competent for follicle rupture/ovulation. Here, we identify a zinc-finger transcription factor Hindsight (Hnt) as the first transcription factor regulating follicle's competency for ovulation in Drosophila. Hnt is not expressed in immature stage-13 follicle cells but is upregulated in mature stage-14 follicle cells, which is essential for follicle rupture/ovulation. Hnt upregulates Mmp2 expression in posterior follicle cells (essential for the breakdown of the follicle wall) and Oamb expression in all follicle cells (the receptor for receiving adrenergic signaling and inducing Mmp2 activation). Hnt's role in regulating Mmp2 and Oamb can be replaced by its human homolog Ras-responsive element-binding protein 1 (RREB-1). Our data suggest that Hnt/RREB-1 plays conserved role in regulating follicle maturation and competency for ovulation.

A major challenge in the application of the CRISPR-Cas13d system is to accurately predict its guide-dependent on-target and off-target effect. Here, we perform CRISPR-Cas13d proliferation screens and design a deep learning model, named DeepCas13, to predict the on-target activity from guide sequences and secondary structures. DeepCas13 outperforms existing methods to predict the efficiency of guides targeting both protein-coding and non-coding RNAs. Guides targeting non-essential genes display off-target viability effects, which are closely related to their on-target efficiencies. Choosing proper negative control guides during normalization mitigates the associated false positives in proliferation screens. We apply DeepCas13 to the guides targeting lncRNAs, and identify lncRNAs that affect cell viability and proliferation in multiple cell lines. The higher prediction accuracy of DeepCas13 over existing methods is extensively confirmed via a secondary CRISPR-Cas13d screen and quantitative RT-PCR experiments. DeepCas13 is freely accessible via http://deepcas13.weililab.org .

Recognition of modified histones by "reader" proteins constitutes a key mechanism regulating diverse chromatin-associated processes important for normal and neoplastic development. We recently identified the YEATS domain as a novel acetyllysine-binding module; however, the functional importance of YEATS domain-containing proteins in human cancer remains largely unknown. Here, we show that the YEATS2 gene is highly amplified in human non-small cell lung cancer (NSCLC) and is required for cancer cell growth and survival. YEATS2 binds to acetylated histone H3 via its YEATS domain. The YEATS2-containing ATAC complex co-localizes with H3K27 acetylation (H3K27ac) on the promoters of actively transcribed genes. Depletion of YEATS2 or disruption of the interaction between its YEATS domain and acetylated histones reduces the ATAC complex-dependent promoter H3K9ac levels and deactivates the expression of essential genes. Taken together, our study identifies YEATS2 as a histone H3K27ac reader that regulates a transcriptional program essential for NSCLC tumorigenesis.

In bioinformatics, the rapid development of gene sequencing technology has produced an increasing amount of microarray data. This type of data shares the typical characteristics of small sample size and high feature dimensions. Searching for biomarkers from microarray data, which expression features of various diseases, is essential for the disease classification. feature selection has therefore became fundemental for the analysis of microarray data, which designs to remove irrelevant and redundant features. There are a large number of redundant features and irrelevant features in microarray data, which severely degrade the classification effectiveness. We propose an innovative feature selection method with the goal of obtaining feature dependencies from a priori knowledge and removing redundant features using spectral clustering. In this paper, the graph structure is firstly constructed by using the gene interaction network as a priori knowledge, and then a link prediction method based on graph neural network is proposed to enhance the graph structure data. Finally, a feature selection method based on spectral clustering is proposed to determine biomarkers. The classification accuracy on DLBCL and Prostate can be improved by 10.90% and 16.22% compared to traditional methods. Link prediction provides an average classification accuracy improvement of 1.96% and 1.31%, and is up to 16.98% higher than the published method. The results show that the proposed method can have full use of a priori knowledge to effectively select disease prediction biomarkers with high classification accuracy.

The MYC oncogene is often dysregulated in human cancer, including hepatocellular carcinoma (HCC). MYC is considered undruggable to date. Here, we comprehensively identify genes essential for survival of MYChigh but not MYClow cells by a CRISPR/Cas9 genome-wide screen in a MYC-conditional HCC model. Our screen uncovers novel MYC synthetic lethal (MYC-SL) interactions and identifies most MYC-SL genes described previously. In particular, the screen reveals nucleocytoplasmic transport to be a MYC-SL interaction. We show that the majority of MYC-SL nucleocytoplasmic transport genes are upregulated in MYChigh murine HCC and are associated with poor survival in HCC patients. Inhibiting Exportin-1 (XPO1) in vivo induces marked tumor regression in an autochthonous MYC-transgenic HCC model and inhibits tumor growth in HCC patient-derived xenografts. XPO1 expression is associated with poor prognosis only in HCC patients with high MYC activity. We infer that MYC may generally regulate and require altered expression of nucleocytoplasmic transport genes for tumorigenesis.

NSUN5, encoding a cytosine-5 RNA methyltransferase and located in the 7q11.23 locus, is a candidate gene for tetralogy of Fallot (TOF). Deletion of the 7q11.23 locus in humans is linked to cardiac outflow tract (OFT) disorders including TOF. We identified four potential pathogenic mutations in the coding region of NSUN5 and which were enriched in TOF patients by an association study of 132 TOF patients and 2,000 in-house controls (P = 1.44 × 10-5). We then generated a Nsun5 null (Nsun5 -/-) mouse model to validate the human findings by defining the functions of Nsun5 in OFT morphogenesis. The OFT did not develop properly in the Nsun5 deletion embryonic heart. We found a misalignment of the aorta and septum defects caused by the delayed fusion of the membraneous ventricular spetum as an OFT development delay. This caused OFT development delay in 27 of 64 (42.2%) Nsun5 -/- mice. Moreover, we also found OFT development delay in 8 of 51 (15.7%) Nsun5 +/- mice. Further functional experiments showed that the loss of Nsun5 function impaired the 5-methylcytosine (m5C) modification and translation efficiency of essential cardiac genes. Nsun5 is required for normal OFT morphogenesis and it regulates the m5C modification of essential cardiac genes. Our findings suggest the involvement of NSUN5 in the pathogenesis of TOF.

The sperm flagellum is essential for male fertility, and defects in flagellum biogenesis are associated with male infertility. Deficiency of coiled-coil domain-containing (CCDC) 42 (CCDC42) is specifically associated with malformation of mouse sperm flagella. Here, we find that the testis-specific protein CCDC38 interacts with CCDC42, localizing on the manchette and sperm tail during spermiogenesis. Inactivation of CCDC38 in male mice results in a distorted manchette, multiple morphological abnormalities of the flagella of spermatozoa and eventually male sterility. Furthermore, we find that CCDC38 interacts with intraflagellar transport protein 88 (IFT88), as well as outer dense fibrous 2 (ODF2), and the knockout of Ccdc38 reduces transport of ODF2 to the flagellum. Altogether, our results uncover the essential role of CCDC38 in sperm flagellum biogenesis, and suggest that some mutations of these genes might be associated with male infertility in humans.

Monogenic hypertension is rare and caused by genetic mutations, but whether factors associated with mutations are disease-specific remains uncertain. Given two factors associated with high mutation rates, we tested how many previously known genes match with (i) proximity to telomeres or (ii) high adenine and thymine content in cardiovascular diseases (CVDs) related to vascular stiffening. We extracted genomic information using a genome data viewer. In human chromosomes, 64 of 79 genetic loci involving >25 rare mutations and single nucleotide polymorphisms satisfied (i) or (ii), resulting in an 81% matching rate. However, this high matching rate was no longer observed as we checked the two factors in genes associated with essential hypertension (EH), thoracic aortic aneurysm (TAA), and congenital heart disease (CHD), resulting in matching rates of 53%, 70%, and 75%, respectively. A matching of telomere proximity or high adenine and thymine content projects the list of loci involving rare mutations of monogenic hypertension better than those of other CVDs, likely due to adoption of rigorous criteria for true-positive signals. Our data suggest that the factor-disease matching rate is an accurate tool that can explain deleterious mutations of monogenic hypertension at a >80% match-unlike the relatively lower matching rates found in human genes of EH, TAA, CHD, and familial Parkinson's disease.

The histone variant H2A.Z is essential for maintaining embryonic stem cell (ESC) identity in part by keeping developmental genes in a poised bivalent state. However, how H2A.Z is deposited into the bivalent domains remains unknown. In mammals, two chromatin remodeling complexes, Tip60/p400 and SRCAP, exchange the canonical histone H2A for H2A.Z in the chromatin. Here we show that Glioma Amplified Sequence 41 (Gas41), a shared subunit of the two H2A.Z-depositing complexes, functions as a reader of histone lysine acetylation and recruits Tip60/p400 and SRCAP to deposit H2A.Z into specific chromatin regions including bivalent domains. The YEATS domain of Gas41 bound to acetylated histone H3K27 and H3K14 both in vitro and in cells. The crystal structure of the Gas41 YEATS domain in complex with the H3K27ac peptide revealed that, similar to the AF9 and ENL YEATS domains, Gas41 YEATS forms a serine-lined aromatic cage for acetyllysine recognition. Consistently, mutations in the aromatic residues of the Gas41 YEATS domain abrogated the interaction. In mouse ESCs, knockdown of Gas41 led to flattened morphology of ESC colonies, as the result of derepression of differentiation genes. Importantly, the abnormal morphology was rescued by expressing wild-type Gas41, but not the YEATS domain mutated counterpart that does not recognize histone acetylation. Mechanically, we found that Gas41 depletion led to reduction of H2A.Z levels and a concomitant reduction of H3K27me3 levels on bivalent domains. Together, our study reveals an essential role of the Gas41 YEATS domain in linking histone acetylation to H2A.Z deposition and maintenance of ESC identity.

N7-methylguanosine (m7G) modification, routinely occurring at mRNA 5' cap or within tRNAs/rRNAs, also exists internally in messenger RNAs (mRNAs). Although m7G-cap is essential for pre-mRNA processing and protein synthesis, the exact role of mRNA internal m7G modification remains elusive. Here, we report that mRNA internal m7G is selectively recognized by Quaking proteins (QKIs). By transcriptome-wide profiling/mapping of internal m7G methylome and QKI-binding sites, we identified more than 1,000 high-confidence m7G-modified and QKI-bound mRNA targets with a conserved "GANGAN (N = A/C/U/G)" motif. Strikingly, QKI7 interacts (via C terminus) with the stress granule (SG) core protein G3BP1 and shuttles internal m7G-modified transcripts into SGs to regulate mRNA stability and translation under stress conditions. Specifically, QKI7 attenuates the translation efficiency of essential genes in Hippo signaling pathways to sensitize cancer cells to chemotherapy. Collectively, we characterized QKIs as mRNA internal m7G-binding proteins that modulate target mRNA metabolism and cellular drug resistance.

The genus Streptomyces is a widespread genus within the phylum Actinobacteria and has been isolated from various environments worldwide. However, little is known about whether biogeography affects distributional pattern of Streptomyces in salty environments. Such information is essential for understanding the ecology of Streptomyces. Here we analyzed four house-keeping genes (16S rRNA, rpoB, recA and atpD) and salty-tolerance related genes (ectA-ectD) of 38 Streptomyces strains isolated from saline environments in Yunnan and Xinjiang Provinces of western China. The obtained Streptomyces strains were classified into three operational taxonomic units, each comprising habitat-specific geno- and ecotype STs. In combination with expressional variations of salty-tolerance related genes, the statistical analyses showed that spatial distance and environmental factors substantially influenced Streptomyces distribution in saline environments: the former had stronger influence at large spatial scales (>700 km), whereas the latter was influential at large (>700 km) and small spatial scales (<700 km). Plus, the quantitative analyses of salty-tolerence related genes (ectA-D) indicated that Streptomyces strains from salt lakes have higher expression of ectA-D genes and could accumulate larger quantities of ectoine and hydroxyectoine than strains from salt mines, which could help them resist to salinity in the hypersaline environments.

The nucleosome remodeling and deacetylase (NuRD) complex is essential for gene expression and cell fate determination, and missense mutations of NuRD caused neurodevelopmental diseases. However, the molecular pathogenesis of clinic NuRD variants is unknown. Here, we introduced a clinic CHD3 (L915F) variant into Caenorhabditis elegans homologue LET-418, impairing germline and vulva development and ultimately causing animal sterility. Our ATAC-seq and RNA-seq analyses revealed that this variant generated an abnormal open chromatin structure and disrupted the expression of developmental genes. Through genetic suppressor screens, we uncovered that intragenic mutations, likely renovating NuRD activity, restored animal viability. We also found that intergenic mutations in nucleosome remodeling factor NURF that counteracts NuRD rescued abnormal chromatin structure, gene expression, and animal sterility. We propose that two antagonistic chromatin-remodeling factors coordinate to establish the proper chromatin status and transcriptome and that inhibiting NURF may provide insights for treatment of NuRD mutation-related diseases.

Over the past decades, studies have reported that the combinatorial regulation of transcription factors (TFs) and microRNAs (miRNAs) is essential for the appropriate execution of biological events and developmental processes. Dysregulations of these regulators often cause diseases. However, there are no available resources on the regulatory cascades of TFs and miRNAs in the context of human diseases. To fulfill this vacancy, we established the TMREC database in this study. First, we integrated curated transcriptional and post-transcriptional regulations to construct the TF and miRNA regulatory network. Next, we identified all linear paths using the Breadth First Search traversal method. Finally, we used known disease-related genes and miRNAs to measure the strength of association between cascades and diseases. Currently, TMREC consists of 74,248 cascades and 25,194 cascade clusters, involving in 412 TFs, 266 miRNAs and 545 diseases. With the expanding of experimental support regulation data, we will regularly update the database. TMREC aims to help experimental biologists to comprehensively analyse gene expression regulation, to understand the aetiology and to predict novel therapeutic targets. TMREC is freely available at http://bioinfo.hrbmu.edu.cn/TMREC/.

Soybean (Glycine max) was domesticated from its wild relative Glycine soja. However, the genetic variations underlying soybean domestication are not well known. Comparative transcriptomics revealed that a small portion of the orthologous genes might have been fast evolving. In contrast, three gene expression clusters were identified as divergent by their expression patterns, which occupied 37.44% of the total genes, hinting at an essential role for gene expression alteration in soybean domestication. Moreover, the most divergent stage in gene expression between wild and cultivated soybeans occurred during seed development around the cotyledon stage (15 d after fertilization, G15). A module in which the co-expressed genes were significantly down-regulated at G15 of wild soybeans was identified. The divergent clusters and modules included substantial differentially expressed genes (DEGs) between wild and cultivated soybeans related to cell division, storage compound accumulation, hormone response, and seed maturation processes. Chromosomal-linked DEGs, quantitative trait loci controlling seed weight and oil content, and selection sweeps revealed candidate DEGs at G15 in the fruit-related divergence of G. max and G. soja. Our work establishes a transcriptomic selection mechanism for altering gene expression during soybean domestication, thus shedding light on the molecular networks underlying soybean seed development and breeding strategy.