Synexpression groups in eukaryotes

An "R-spondin code" for multimodal signaling ON-OFF states

R-spondins (RSPOs) are a family of secreted proteins and stem cell growth factors that are potent co-activators of Wnt signaling. Recently, RSPO2 and RSPO3 were shown to be multifunctional, not only amplifying Wnt- but also binding BMP- and FGF receptors to downregulate signaling. The common mechanism underlying these diverse functions is that RSPO2 and RSPO3 act as "endocytosers" that link transmembrane proteins to ZNRF3/RNF43 E3 ligases and trigger target internalization. Thus, RSPOs are natural protein targeting chimeras for cell surface proteins. Conducting data mining and cell surface binding assays we report additional candidate RSPO targets, including SMO, PTC1,2, LGI1, ROBO4, and PTPR(F/S). We propose that there is an "R-spondin code" that imparts combinatorial signaling ON-OFF states of multiple growth factors. This code involves the modular RSPO domains, notably distinct motifs in the divergent RSPO-TSP1 domains to mediate target interaction and internalization. The RSPO code offers a novel framework for the understanding how diverse signaling pathways may be coordinately regulated in development and disease.

CLOCK evolved in cnidaria to synchronize internal rhythms with diel environmental cues

The circadian clock enables anticipation of the day/night cycle in animals ranging from cnidarians to mammals. Circadian rhythms are generated through a transcription-translation feedback loop (TTFL or pacemaker) with CLOCK as a conserved positive factor in animals. However, CLOCK's functional evolutionary origin and mechanism of action in basal animals are unknown. In the cnidarian Nematostella vectensis, pacemaker gene transcript levels, including NvClk (the Clock ortholog), appear arrhythmic under constant darkness, questioning the role of NvCLK. Utilizing CRISPR/Cas9, we generated a NvClk allele mutant (NvClkΔ), revealing circadian behavior loss under constant dark (DD) or light (LL), while maintaining a 24 hr rhythm under light-dark condition (LD). Transcriptomics analysis revealed distinct rhythmic genes in wild-type (WT) polypsunder LD compared to DD conditions. In LD, NvClkΔ/Δ polyps exhibited comparable numbers of rhythmic genes, but were reduced in DD. Furthermore, under LD, the NvClkΔ/Δ polyps showed alterations in temporal pacemaker gene expression, impacting their potential interactions. Additionally, differential expression of non-rhythmic genes associated with cell division and neuronal differentiation was observed. These findings revealed that a light-responsive pathway can partially compensate for circadian clock disruption, and that the Clock gene has evolved in cnidarians to synchronize rhythmic physiology and behavior with the diel rhythm of the earth's biosphere.

The pelvic organs receive no parasympathetic innervation

The pelvic organs (bladder, rectum, and sex organs) have been represented for a century as receiving autonomic innervation from two pathways - lumbar sympathetic and sacral parasympathetic - by way of a shared relay, the pelvic ganglion, conceived as an assemblage of sympathetic and parasympathetic neurons. Using single-cell RNA sequencing, we find that the mouse pelvic ganglion is made of four classes of neurons, distinct from both sympathetic and parasympathetic ones, albeit with a kinship to the former, but not the latter, through a complex genetic signature. We also show that spinal lumbar preganglionic neurons synapse in the pelvic ganglion onto equal numbers of noradrenergic and cholinergic cells, both of which therefore serve as sympathetic relays. Thus, the pelvic viscera receive no innervation from parasympathetic or typical sympathetic neurons, but instead from a divergent tail end of the sympathetic chains, in charge of its idiosyncratic functions.

Genomic investigation of duplication, functional conservation, and divergence in the LRR-RLK Family of Saccharum

BACKGROUND

Sugarcane (Saccharum spp.) holds exceptional global significance as a vital crop, serving as a primary source of sucrose, bioenergy, and various by-products. The optimization of sugarcane breeding by fine-tuning essential traits has become crucial for enhancing crop productivity and stress resilience. Leucine-rich repeat receptor-like kinases (LRR-RLK) genes present promising targets for this purpose, as they are involved in various aspects of plant development and defense processes.

RESULTS

Here, we present a detailed overview of phylogeny and expression of 288 (495 alleles) and 312 (1365 alleles) LRR-RLK genes from two founding Saccharum species, respectively. Phylogenetic analysis categorized these genes into 15 subfamilies, revealing considerable expansion or reduction in certain LRR-type subfamilies. Compared to other plant species, both Saccharum species had more significant LRR-RLK genes. Examination of cis-acting elements demonstrated that SsLRR-RLK and SoLRR-RLK genes exhibited no significant difference in the types of elements included, primarily involved in four physiological processes. This suggests a broad conservation of LRR-RLK gene function during Saccharum evolution. Synteny analysis indicated that all LRR-RLK genes in both Saccharum species underwent gene duplication, primarily through whole-genome duplication (WGD) or segmental duplication. We identified 28 LRR-RLK genes exhibiting novel expression patterns in response to different tissues, gradient development leaves, and circadian rhythm in the two Saccharum species. Additionally, SoLRR-RLK104, SoLRR-RLK7, SoLRR-RLK113, and SsLRR-RLK134 were identified as candidate genes for sugarcane disease defense response regulators through transcriptome data analysis of two disease stresses. This suggests LRR-RLK genes of sugarcane involvement in regulating various biological processes, including leaf development, plant morphology, photosynthesis, maintenance of circadian rhythm stability, and defense against sugarcane diseases.

CONCLUSIONS

This investigation into gene duplication, functional conservation, and divergence of LRR-RLK genes in two founding Saccharum species lays the groundwork for a comprehensive genomic analysis of the entire LRR-RLK gene family in Saccharum. The results reveal LRR-RLK gene played a critical role in Saccharum adaptation to diverse conditions, offering valuable insights for targeted breeding and precise phenotypic adjustments.

Evolution and expression analysis of the caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) gene family in jute (Corchorus L.)

BACKGROUND

Jute is considered one of the most important crops for fiber production and multipurpose usages. Caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) is a crucial enzyme involved in lignin biosynthesis in plants. The potential functions of CCoAOMT in lignin biosynthesis of jute have been reported in several studies. However, little is known about the evolution of the CCoAOMT gene family, and either their expression level at different developing stages in different jute cultivars, as well as under abiotic stresses including salt and drought stress.

RESULTS

In the present study, 66 CCoAOMT genes from 12 species including 12 and eight CCoAOMTs in Corchorus olitorius and C. capsularis were identified. Phylogenetic analysis revealed that CCoAOMTs could be divided into six groups, and gene expansion was observed in C. olitorius. Furthermore, gene expression analysis of developing jute fibers was conducted at different developmental stages (15, 30, 45, 60, and 90 days after sowing [DAS]) in six varieties (Jute-179 [J179], Lubinyuanguo [LB], and Qiongyueqing [QY] for C. capsularis; Funong No.5 [F5], Kuanyechangguo [KY], and Cvlv [CL] for C. olitorius). The results showed that CCoAOMT1 and CCoAOMT2 were the dominant genes in the CCoAOMT family. Of these two dominant CCoAOMTs, CCoAOMT2 showed a constitutive expression level during the entire growth stages, while CCoAOMT1 exhibited differential expression patterns. These two genes showed higher expression levels in C. olitorius than in C. capsularis. The correlation between lignin content and CCoAOMT gene expression levels indicated that this gene family influences the lignin content of jute. Using real-time quantitative reverse transcription PCR (qRT-PCR), a substantial up-regulation of CCoAOMTs was detected in stem tissues of jute 24 h after drought treatment, with an up to 17-fold increase in expression compared to that of untreated plants.

CONCLUSIONS

This study provides a basis for comprehensive genomic studies of the entire CCoAOMT gene family in C. capsularis and C. olitorius. Comparative genomics analysis among the CCoAOMT gene families of 12 species revealed the close evolutionary relationship among Corchorus, Theobroma cacao and Gossypium raimondii. This study also shows that CCoAOMTs are not only involved in lignin biosynthesis, but also are associated with the abiotic stress response in jute, and suggests the potential use of these lignin-related genes to genetically improve the fiber quality of jute.

Insulin secretion deficits in a Prader-Willi syndrome β-cell model are associated with a concerted downregulation of multiple endoplasmic reticulum chaperones

Prader-Willi syndrome (PWS) is a multisystem disorder with neurobehavioral, metabolic, and hormonal phenotypes, caused by loss of expression of a paternally-expressed imprinted gene cluster. Prior evidence from a PWS mouse model identified abnormal pancreatic islet development with retention of aged insulin and deficient insulin secretion. To determine the collective roles of PWS genes in β-cell biology, we used genome-editing to generate isogenic, clonal INS-1 insulinoma lines having 3.16 Mb deletions of the silent, maternal- (control) and active, paternal-allele (PWS). PWS β-cells demonstrated a significant cell autonomous reduction in basal and glucose-stimulated insulin secretion. Further, proteomic analyses revealed reduced levels of cellular and secreted hormones, including all insulin peptides and amylin, concomitant with reduction of at least ten endoplasmic reticulum (ER) chaperones, including GRP78 and GRP94. Critically, differentially expressed genes identified by whole transcriptome studies included reductions in levels of mRNAs encoding these secreted peptides and the group of ER chaperones. In contrast to the dosage compensation previously seen for ER chaperones in Grp78 or Grp94 gene knockouts or knockdown, compensation is precluded by the stress-independent deficiency of ER chaperones in PWS β-cells. Consistent with reduced ER chaperones levels, PWS INS-1 β-cells are more sensitive to ER stress, leading to earlier activation of all three arms of the unfolded protein response. Combined, the findings suggest that a chronic shortage of ER chaperones in PWS β-cells leads to a deficiency of protein folding and/or delay in ER transit of insulin and other cargo. In summary, our results illuminate the pathophysiological basis of pancreatic β-cell hormone deficits in PWS, with evolutionary implications for the multigenic PWS-domain, and indicate that PWS-imprinted genes coordinate concerted regulation of ER chaperone biosynthesis and β-cell secretory pathway function.

MCC is a centrosomal protein that relocalizes to non-centrosomal apical sites during intestinal cell differentiation

The gene mutated in colorectal cancer (MCC) encodes a coiled-coil protein implicated, as its name suggests, in the pathogenesis of hereditary human colon cancer. To date, however, the contributions of MCC to intestinal homeostasis and disease remain unclear. Here, we examine the subcellular localization of MCC, both at the mRNA and protein levels, in the adult intestinal epithelium. Our findings reveal that Mcc transcripts are restricted to proliferating crypt cells, including Lgr5+ stem cells, where the Mcc protein is distinctly associated with the centrosome. Upon intestinal cellular differentiation, Mcc is redeployed to the apical domain of polarized villus cells where non-centrosomal microtubule organizing centers (ncMTOCs) are positioned. Using intestinal organoids, we show that the shuttling of the Mcc protein depends on phosphorylation by casein kinases 1δ and ε, which are critical modulators of WNT signaling. Together, our findings support a role for MCC in establishing and maintaining the cellular architecture of the intestinal epithelium as a component of both the centrosome and ncMTOC.

Expression Profiling and MicroRNA Regulatory Networks of Homeobox Family Genes in Sugarcane Saccharum spontaneum L

Homeobox (HB) genes play important roles in plant growth and development processes, particularly in the formation of lateral organs. Thus, they could influence leaf morphogenesis and biomass formation in plants. However, little is known about HBs in sugarcane, a crucial sugar crop, due to its complex genetic background. Here, 302 allelic sequences for 104 HBs were identified and divided into 13 subfamilies in sugarcane Saccharum spontaneum. Comparative genomics revealed that whole-genome duplication (WGD)/segmental duplication significantly promoted the expansion of the HB family in S. spontaneum, with SsHB26, SsHB63, SsHB64, SsHB65, SsHB67, SsHB95, and SsHB96 being retained from the evolutionary event before the divergence of dicots and monocots. Based on the analysis of transcriptome and degradome data, we speculated that SsHB15 and SsHB97 might play important roles in regulating sugarcane leaf morphogenesis, with miR166 and SsAGO10 being involved in the regulation of SsHB15 expression. Moreover, subcellular localization and transcriptional activity detection assays demonstrated that these two genes, SsHB15 and SsHB97, were functional transcription factors. This study demonstrated the evolutionary relationship and potential functions of SsHB genes and will enable the further investigation of the functional characterization and the regulatory mechanisms of SsHBs.

Gene function prediction in five model eukaryotes exclusively based on gene relative location through machine learning

The function of most genes is unknown. The best results in automated function prediction are obtained with machine learning-based methods that combine multiple data sources, typically sequence derived features, protein structure and interaction data. Even though there is ample evidence showing that a gene's function is not independent of its location, the few available examples of gene function prediction based on gene location rely on sequence identity between genes of different organisms and are thus subjected to the limitations of the relationship between sequence and function. Here we predict thousands of gene functions in five model eukaryotes (Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster, Mus musculus and Homo sapiens) using machine learning models exclusively trained with features derived from the location of genes in the genomes to which they belong. Our aim was not to obtain the best performing method to automated function prediction but to explore the extent to which a gene's location can predict its function in eukaryotes. We found that our models outperform BLAST when predicting terms from Biological Process and Cellular Component Ontologies, showing that, at least in some cases, gene location alone can be more useful than sequence to infer gene function.

Evolutionary and Gene Expression Analyses Reveal New Insights into the Role of LSU Gene-Family in Plant Responses to Sulfate-Deficiency

LSU proteins belong to a plant-specific gene family initially characterized by their strong induction in response to sulfate (S) deficiency. In the last few years, LSUs have arisen as relevant hubs in protein–protein interaction networks, in which they play relevant roles in the response to abiotic and biotic stresses. Most of our knowledge on LSU genomic organization, expression and function comes from studies in Arabidopsis and tobacco, while little is known about the LSU gene repertoire and evolution of this family in land plants. In this work, a total of 270 LSU family members were identified using 134 land plant species with whole-genome sequences available. Phylogenetic analysis revealed that LSU genes belong to a Spermatophyta-specific gene family, and their homologs are distributed in three major groups, two for dicotyledons and one group for monocotyledons. Protein sequence analyses showed four new motifs that further support the subgroup classification by phylogenetic analyses. Moreover, we analyzed the expression of LSU genes in one representative species of each phylogenetic group (wheat, tomato and Arabidopsis) and found a conserved response to S deficiency, suggesting that these genes might play a key role in S stress responses. In summary, our results indicate that LSU genes belong to the Spermatophyta-specific gene family and their response to S deficiency is conserved in angiosperms.

Functional Modules in Gametogenesis

Gametogenesis, the production of eggs and sperm, is a fundamental process in sexually reproducing animals. Following gametogenesis commitment and sexual fate decision, germ cells undergo several developmental processes to halve their genomic size and acquire sex-specific characteristics of gametes, including cellular size, motility, and cell polarity. However, it remains unclear how different gametogenesis processes are initially integrated. With the advantages of the teleost fish medaka (Oryzias latipes), in which germline stem cells continuously produce eggs and sperm in mature gonads and a sexual switch gene in germ cells is identified, we found that distinct pathways initiate gametogenesis cooperatively after commitment to gametogenesis. This evokes the concept of functional modules, in which functionally interlocked genes are grouped to yield distinct gamete characteristics. The various combinations of modules may allow us to explain the evolution of diverse reproductive systems, such as parthenogenesis and hermaphroditism.

Major Vault Protein (MVP) Associated With BRAF V600E Mutation Is an Immune Microenvironment-Related Biomarker Promoting the Progression of Papillary Thyroid Cancer via MAPK/ERK and PI3K/AKT Pathways

Papillary thyroid cancer (PTC) is the most common malignancy of the endocrine system, with an increase in incidence frequency. Major vault protein (MVP) is the main structural protein of the vault complex that has already been investigated in specific cancers. Yet the underlying biological functions and molecular mechanisms of MVP in PTC still remain considerably uncharacterized. Comprehensive analyses are predicated on several public datasets and local RNA-Seq cohort. Clinically, we found that MVP was upregulated in human PTC than in non-cancerous thyroid tissue and was correlated with vital clinicopathological parameters in PTC patients. MVP expression was associated with BRAF V600E, RAS, TERT, and RET status, and it was correlated with worse progression-free survival in PTC patients. Functionally, enrichment analysis provided new clues for the close relationship between MVP with cancer-related signaling pathways and the immune microenvironment in PTC. In PTC with high MVP expression, we found CD8+ T cells, regulatory T cells, and follicular helper T cells have a higher infiltration level. Intriguingly, MVP expression was positively correlated with multiple distinct phases of the anti-cancer immunity cycle. MVP knockdown significantly suppressed cell viability and colony formation, and promoted apoptosis. In addition, downregulated MVP markedly inhibited the migration and invasion potential of PTC cells. The rescue experiments showed that MVP could reverse the level of cell survival and migration. Mechanistically, MVP exerts its oncogenic function in PTC cells through activating PI3K/AKT/mTOR and MAPK/ERK pathways. These results point out that MVP is a reliable biomarker related to the immune microenvironment and provide a basis for elucidating the oncogenic roles of MVP in PTC progression.

Identification of RHEX as a novel biomarker related to progression and immunity of non-small cell lung carcinoma

Background

The therapeutic response and prognosis of patients with non-small cell lung carcinoma (NSCLC) are widely related to immunity. To improve the prognosis of patients and provide reliable information to guide appropriate personalized treatment strategies, it is necessary to identify reliable prognostic or predictive indicators closely related to tumor phenotype and immune traits in NSCLC.

Methods

Based on The Cancer Genome Atlas (TCGA)-NSCLC mRNA expression profile data, a novel approach combining differential gene expression analysis, single-sample gene set enrichment analysis (ssGSEA), and weighted gene co-expression network analysis (WGCNA) was used to screen hub genes. Subsequently, the regulator of hemoglobinization and erythroid cell expansion (RHEX) was identified as a key gene using the log-rank test and confirmed in the ArrayExpress database. The relationship between RHEX and clinicopathological parameters was analyzed using the Wilcoxon rank-sum test. More importantly, through gene set enrichment analysis (GSEA) and cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) algorithms, and with reference to the Tumor IMmune Estimation Resource (TIMER) database, we explored the relevant pathways of RHEX and its relationship with tumor-infiltrating immune cells (TICs). Finally, we depicted the association between RHEX and immunomodulators in the TCGA and a web portal TISIDB.

Results

The RHEX mRNA expression levels in tumor tissues were lower than those in normal tissues and declined with the progression of NSCLC. Meanwhile, RHEX overexpression was associated with high immune infiltration levels and a favorable clinical prognosis. RHEX may participate in tumor microenvironment (TME) regulation through multiple tumor-immune related pathways, especially the JAK-STAT signaling pathway. Furthermore, RHEX expression affected the infiltrating abundance of multiple TICs and positively correlated with most of the immunomodulators in NSCLC.

Conclusions

Our study is the first to propose that RHEX is an immune-related gene with prognostic value in NSCLC and reveals the underlying mechanism between RHEX and tumor-immune system interactions. These results ultimately provide guidance for prognosis and immunotherapy for NSCLC patients.

Detection and identification of cis-regulatory elements using change-point and classification algorithms

BACKGROUND

Transcriptional regulation is primarily mediated by the binding of factors to non-coding regions in DNA. Identification of these binding regions enhances understanding of tissue formation and potentially facilitates the development of gene therapies. However, successful identification of binding regions is made difficult by the lack of a universal biological code for their characterisation.

RESULTS

We extend an alignment-based method, changept, and identify clusters of biological significance, through ontology and de novo motif analysis. Further, we apply a Bayesian method to estimate and combine binary classifiers on the clusters we identify to produce a better performing composite.

CONCLUSIONS

The analysis we describe provides a computational method for identification of conserved binding sites in the human genome and facilitates an alternative interrogation of combinations of existing data sets with alignment data.

Discovery of a genetic module essential for assigning left-right asymmetry in humans and ancestral vertebrates

The vertebrate left-right axis is specified during embryogenesis by a transient organ: the left-right organizer (LRO). Species including fish, amphibians, rodents and humans deploy motile cilia in the LRO to break bilateral symmetry, while reptiles, birds, even-toed mammals and cetaceans are believed to have LROs without motile cilia. We searched for genes whose loss during vertebrate evolution follows this pattern and identified five genes encoding extracellular proteins, including a putative protease with hitherto unknown functions that we named ciliated left-right organizer metallopeptide (CIROP). Here, we show that CIROP is specifically expressed in ciliated LROs. In zebrafish and Xenopus, CIROP is required solely on the left side, downstream of the leftward flow, but upstream of DAND5, the first asymmetrically expressed gene. We further ascertained 21 human patients with loss-of-function CIROP mutations presenting with recessive situs anomalies. Our findings posit the existence of an ancestral genetic module that has twice disappeared during vertebrate evolution but remains essential for distinguishing left from right in humans.

NTF3 Correlates With Prognosis and Immune Infiltration in Hepatocellular Carcinoma

Background: The potential role of Neurotrophic factor-3(NTF3) in liver cancer is unknown. Therefore, we aimed to explore the clinical value of NTF3 in hepatocellular carcinoma (HCC).

Methods: We used a variety of databases to analyze the expression, relationship with prognosis and immune significance of NTF3 in liver cancer through bioinformatics.

Results: NTF3 was low expressed in HCC and was an independent prognostic factor in patients with HCC. CIBERSORT analysis indicated that NTF3 expression was positively correlated with CD4+ cells, mast cells, NK cells, macrophages and B cells in the tumor microenvironment. Furthermore, we found that NTF3 expression was negatively correlated with the immune checkpoints PD-L1, TIGIT and TIM-3. Functional network analysis revealed that NTF3 regulates HCC progression through a variety of cancer-related kinases, transcription factors and signaling pathways.

Conclusions: We demonstrate that NTF3 correlates with prognosis and immune infiltration in HCC.

Deletion of pancreas-specific miR-216a reduces beta-cell mass and inhibits pancreatic cancer progression in mice

miRNAs have crucial functions in many biological processes and are candidate biomarkers of disease. Here, we show that miR-216a is a conserved, pancreas-specific miRNA with important roles in pancreatic islet and acinar cells. Deletion of miR-216a in mice leads to a reduction in islet size, β-cell mass, and insulin levels. Single-cell RNA sequencing reveals a subpopulation of β-cells with upregulated acinar cell markers under a high-fat diet. miR-216a is induced by TGF-β signaling, and inhibition of miR-216a increases apoptosis and decreases cell proliferation in pancreatic cells. Deletion of miR-216a in the pancreatic cancer-prone mouse line KrasG12D;Ptf1aCreER reduces the propensity of pancreatic cancer precursor lesions. Notably, circulating miR-216a levels are elevated in both mice and humans with pancreatic cancer. Collectively, our study gives insights into how β-cell mass and acinar cell growth are modulated by a pancreas-specific miRNA and also suggests miR-216a as a potential biomarker for diagnosis of pancreatic diseases.

Clustering of Sex-Biased Genes and Transposable Elements in the Genome of the Medaka Fish Oryzias latipes

Although genes with similar expression patterns are sometimes found in the same genomic regions, almost nothing is known about the relative organization in genomes of genes and transposable elements (TEs), which might influence each other at the regulatory level. In this study, we used transcriptomic data from male and female gonads of the Japanese medaka Oryzias latipes to define sexually biased genes and TEs and analyze their relative genomic localization. We identified 20,588 genes expressed in the adult gonads of O. latipes. Around 39% of these genes are differentially expressed between male and female gonads. We further analyzed the expression of TEs using the program SQuIRE and showed that more TE copies are overexpressed in testis than in ovaries (36% vs. 10%, respectively). We then developed a method to detect genomic regions enriched in testis- or ovary-biased genes. This revealed that sex-biased genes and TEs are not randomly distributed in the genome and a part of them form clusters with the same expression bias. We also found a correlation of expression between TE copies and their closest genes, which increases with decreasing intervening distance. Such a genomic organization suggests either that TEs hijack the regulatory sequences of neighboring sexual genes, allowing their expression in germ line cells and consequently new insertions to be transmitted to the next generation, or that TEs are involved in the regulation of sexual genes, and might therefore through their mobility participate in the rewiring of sex regulatory networks.

The Organizer and Its Signaling in Embryonic Development

Germ layer specification and axis formation are crucial events in embryonic development. The Spemann organizer regulates the early developmental processes by multiple regulatory mechanisms. This review focuses on the responsive signaling in organizer formation and how the organizer orchestrates the germ layer specification in vertebrates. Accumulated evidence indicates that the organizer influences embryonic development by dual signaling. Two parallel processes, the migration of the organizer’s cells, followed by the transcriptional activation/deactivation of target genes, and the diffusion of secreting molecules, collectively direct the early development. Finally, we take an in-depth look at active signaling that originates from the organizer and involves germ layer specification and patterning.

ABL1 Is a Prognostic Marker and Associated with Immune Infiltration in Hepatocellular Carcinoma

Background

The role of ABL1 in hepatocellular carcinoma (HCC) is still unclear. Therefore, this study aims to explore the potential role of ABL1 in the progression of HCC using bioinformatics methods.

Methods

We analyzed the expression, prognostic potential, and immune cell effect of ABL1 in HCC by using a variety of datasets.

Results

ABL1 is highly expressed in HCC and associated with unfavorable overall survival (OS) and disease-free survival (DFS). Functional network analysis revealed that ABL1 plays an important role in mitochondrial activity, ATP metabolism, protein translation and metabolism, various neurological diseases, nonalcoholic fatty liver disease, and notch signaling pathway. In addition, we found that ABL1 expression was closely correlated with B cells, CD8 + T cells, CD4 + T cells, macrophages, neutrophils, and dendritic cells. Furthermore, ABL1 expression was positively associated with the expression levels of immune checkpoint genes, such as PD-1L, TIM3, TIGIT, and CTLA4.

Conclusion

ABL1 is associated with immune infiltration and prognosis of HCC.

Rab7 is required for mesoderm patterning and gastrulation in Xenopus

Early embryogenesis requires tightly controlled temporal and spatial coordination of cellular behavior and signaling. Modulations are achieved at multiple levels, from cellular transcription to tissue-scale behavior. Intracellularly, the endolysosomal system emerges as an important regulator at different levels, but in vivo studies are rare. In the frog Xenopus, little is known about the developmental roles of endosomal regulators, or their potential involvement in signaling, especially for late endosomes. Here, we analyzed a hypothesized role of Rab7 in this context, a small GTPase known for its role as a late endosomal regulator. First, rab7 showed strong maternal expression. Following localized zygotic transcript enrichment in the mesodermal ring and neural plate, it was found in tailbud-stage neural ectoderm, notochord, pronephros, eyes and neural crest tissues. Inhibition resulted in strong axis defects caused by a requirement of rab7 for mesodermal patterning and correct gastrulation movements. To test a potential involvement in growth factor signaling, we analyzed early Wnt-dependent processes in the mesoderm. Our results suggest a selective requirement for ligand-induced Wnt activation, implicating a context-dependent role of Rab7.

Summary: The late endosomal regulator Rab7 is required for gastrulation movements and axis elongation in Xenopus by regulating early mesoderm patterning.

Genome-wide identification and expression profiling of DREB genes in Saccharum spontaneum

Background

The dehydration-responsive element-binding proteins (DREBs) are important transcription factors that interact with a DRE/CRT (C-repeat) sequence and involve in response to multiple abiotic stresses in plants. Modern sugarcane are hybrids from the cross between Saccharum spontaneum and Saccharum officinarum, and the high sugar content is considered to the attribution of S. officinaurm, while the stress tolerance is attributed to S. spontaneum. To understand the molecular and evolutionary characterization and gene functions of the DREBs in sugarcane, based on the recent availability of the whole genome information, the present study performed a genome-wide in silico analysis of DREB genes and transcriptome analysis in the polyploidy S. spontaneum.

Results

Twelve DREB1 genes and six DREB2 genes were identified in S. spontaneum genome and all proteins contained a conserved AP2/ERF domain. Eleven SsDREB1 allele genes were assumed to be originated from tandem duplications, and two of them may be derived after the split of S. spontaneum and the proximal diploid species sorghum, suggesting tandem duplication contributed to the expansion of DREB1-type genes in sugarcane. Phylogenetic analysis revealed that one DREB2 gene was lost during the evolution of sugarcane. Expression profiling showed different SsDREB genes with variable expression levels in the different tissues, indicating seven SsDREB genes were likely involved in the development and photosynthesis of S. spontaneum. Furthermore, SsDREB1F, SsDREB1L, SsDREB2D, and SsDREB2F were up-regulated under drought and cold condition, suggesting that these four genes may be involved in both dehydration and cold response in sugarcane.

Conclusions

These findings demonstrated the important role of DREBs not only in the stress response, but also in the development and photosynthesis of S. spontaneum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07799-5.

A novel essential protein identification method based on PPI networks and gene expression data

Background

Some proposed methods for identifying essential proteins have better results by using biological information. Gene expression data is generally used to identify essential proteins. However, gene expression data is prone to fluctuations, which may affect the accuracy of essential protein identification. Therefore, we propose an essential protein identification method based on gene expression and the PPI network data to calculate the similarity of "active" and "inactive" state of gene expression in a cluster of the PPI network. Our experiments show that the method can improve the accuracy in predicting essential proteins.

Results

In this paper, we propose a new measure named JDC, which is based on the PPI network data and gene expression data. The JDC method offers a dynamic threshold method to binarize gene expression data. After that, it combines the degree centrality and Jaccard similarity index to calculate the JDC score for each protein in the PPI network. We benchmark the JDC method on four organisms respectively, and evaluate our method by using ROC analysis, modular analysis, jackknife analysis, overlapping analysis, top analysis, and accuracy analysis. The results show that the performance of JDC is better than DC, IC, EC, SC, BC, CC, NC, PeC, and WDC. We compare JDC with both NF-PIN and TS-PIN methods, which predict essential proteins through active PPI networks constructed from dynamic gene expression.

Conclusions

We demonstrate that the new centrality measure, JDC, is more efficient than state-of-the-art prediction methods with same input. The main ideas behind JDC are as follows: (1) Essential proteins are generally densely connected clusters in the PPI network. (2) Binarizing gene expression data can screen out fluctuations in gene expression profiles. (3) The essentiality of the protein depends on the similarity of "active" and "inactive" state of gene expression in a cluster of the PPI network.

Up-regulation of L Antigen Family Member 3 Associates With Aggressive Progression of Breast Cancer

The role of L Antigen Family Member 3 (LAGE3) in breast cancer (BC) has not been sufficiently studied. In this study, we explored the clinical value and biological functions of LAGE3 in BC. Comprehensive analysis of LAGE3 was carried out on The Cancer Genome Atlas, Molecular Taxonomy of Breast Cancer International Consortium and Gene Expression Omnibus datasets. Results showed that LAGE3 expression was higher in BC tissues than in normal breast tissues of public datasets and our local cohort. Moreover, its expression was higher in BC patients with larger tumor size, significant lymph node metastasis, higher tumor grade, and more advanced disease stage. High expression of LAGE3 was correlated with poor prognosis, and LAGE3 could independently predict survival of BC patients. Functional enrichment analysis revealed a correlation between LAGE3 expression and biochemical metabolism and immune-related terms and cancer-related pathways. Analysis of tumor microenvironment indicated that LAGE3 expression was associated with the immune cell infiltration and anti-cancer immunity cycle. LAGE3 expression was higher in triple-negative breast cancer (TNBC) compared to hormone receptor-positive BC, but not HER2-positive subtype. Suppression of LAGE3 expression inhibited the proliferation and induced apoptosis of TNBC cell lines. Besides, the down-regulation of LAGE3 attenuated the migration and invasion but reduced the expression level of epithelial-mesenchymal-transition related proteins in TNBC cell lines. In conclusion, this study demonstrated for the first time that LAGE3 promotes the progression of BC. Therefore, it may be a potential diagnostic and prognostic biomarker, as well as a treatment target for BC.

Isthmin1, a secreted signaling protein, acts downstream of diverse embryonic patterning centers in development

Extracellular signals play essential roles during embryonic patterning by providing positional information in a concentration-dependent manner, and many such signals, like Wnt, fibroblast growth factor (FGF), Hedgehog (Hh), and retinoic acid, act by being secreted into the extracellular space, thereby triggering receptor-mediated responses in other cells. Isthmin1 (ism1) is a secreted protein whose gene expression pattern coincides with that of early dorsal determinants, nodal ligand genes like sqt and cyc, and with fgf8 during various phases of zebrafish development. Ism1 functions in early embryonic patterning and development are poorly understood; however, it has recently been shown to interact with nodal pathway genes to control organ asymmetry in chicken. Here, we show that misexpression of ism1 deletion constructs disrupts embryonic patterning in zebrafish and exhibits genetic interactions with both Fgf and nodal signaling. Unlike Fgf and nodal pathway mutants, CRISPR/Cas9-engineered ism1 mutants did not show obvious developmental defects. Further, in vivo single molecule fluorescence correlation spectroscopy (FCCS) showed that Ism1 diffuses freely in the extra-cellular space, with a diffusion coefficient similar to that of Fgf8a; however, our measurements do not support direct molecular interactions between Ism1 and either nodal ligands or Fgf8a in the developing zebrafish embryo. Together, data from gain- and loss-of-function experiments suggest that zebrafish Ism1 plays a complex role in regulating extracellular signals during early embryonic development.

T-Box Transcription Factor 22 Is an Immune Microenvironment-Related Biomarker Associated With the BRAF V600E Mutation in Papillary Thyroid Carcinoma

Papillary thyroid cancer (PTC) is the most common malignant disease in endocrine systems. T-box transcription factor 22 (TBX22) is a phylogenetically conserved family member that has not been widely characterized in cancers. In this study, we explored the potential clinical significance and biological functions of TBX22 in PTC. Comprehensive analyses of TBX22 were based on the public databases and our local qRT-PCR cohort. We observed that TBX22 was significantly downregulated in PTC compared with normal tissues. TBX22 was associated with several clinicopathological factors in PTC. Low TBX22 expression correlated with BRAF ^V600E and TERT mutation. Functional enrichment analysis revealed that cancer-related pathways and immune progress were closely associated with TBX22 in PTC. In TBX22-low PTC, high immune infiltration levels with increased CD8⁺ T cells, natural killer, M1 macrophages, and T-regulatory cells were observed. TBX22 was negatively correlated with the activity of different steps of the anticancer immunity cycle. Functionally, overexpression of TBX22 inhibited the proliferation, invasion, and migration in PTC cells, while knocking down of TBX22 showed the opposite effects. The present findings disclose that TBX22, as an immune microenvironment-related biomarker, could be an important tumor suppresser gene and might inform the management of PTC patients better.

Identification and validation of L Antigen Family Member 3 as an immune-related biomarker associated with the progression of papillary thyroid cancer

BACKGROUND

Papillary thyroid cancer (PTC) is heterogeneous cancer with many different immune cells involved in its pathogenesis. L Antigen Family Member 3 (LAGE3) is an ESO/LAGE gene family member that has not been extensively studied in PTC.

METHODS

Comprehensive bioinformatics analyses of LAGE3 were based on The Cancer Genome Atlas, Gene Expression Omnibus, and Genomics of Drug Sensitivity in Cancer (GDSC) databases. We also performed RNA-sequencing on 78 paired samples from local PTC patients.

RESULTS

We observed that LAGE3 was significantly up-regulated in most solid tumor types, including PTC compared with corresponding normal tissues. The high level of LAGE3 was also significantly associated with advanced malignancy. LAGE3 expression was significantly associated with cancer-related pathways, biochemical metabolism, and immune-related terms. Further, tumor microenvironment analysis indicated LAGE3 was positively correlated with different immune cells infiltrating levels and the activity of different steps of the cancer-immunity cycle. Analyses based on the GDSC database revealed that low levels of LAGE3 might be resistant to WZ3105, I-BET-762, and PHA-793887. In addition, the experimental results validated that knocking down LAGE3 could affect proliferation, migration, and invasion in the PTC cell lines.

CONCLUSION

This study discloses that LAGE3 plays an oncogenic and cancer-immunological role, also providing novel PTC biological and clinical implications.

R-spondins are BMP receptor antagonists in Xenopus early embryonic development

BMP signaling plays key roles in development, stem cells, adult tissue homeostasis, and disease. How BMP receptors are extracellularly modulated and in which physiological context, is therefore of prime importance. R-spondins (RSPOs) are a small family of secreted proteins that co-activate WNT signaling and function as potent stem cell effectors and oncogenes. Evidence is mounting that RSPOs act WNT-independently but how and in which physiological processes remains enigmatic. Here we show that RSPO2 and RSPO3 also act as BMP antagonists. RSPO2 is a high affinity ligand for the type I BMP receptor BMPR1A/ALK3, and it engages ZNRF3 to trigger internalization and degradation of BMPR1A. In early Xenopus embryos, Rspo2 is a negative feedback inhibitor in the BMP4 synexpression group and regulates dorsoventral axis formation. We conclude that R-spondins are bifunctional ligands, which activate WNT- and inhibit BMP signaling via ZNRF3, with implications for development and cancer.

Identification of a peripheral blood gene signature predicting aortic valve calcification

Calcific aortic valve disease (CAVD) is a significant cause of illness and death worldwide. Identification of early predictive markers could help optimize patient management. RNA-sequencing was carried out on human fetal aortic valves at gestational weeks 9, 13, and 22 and on a case-control study with adult noncalcified and calcified bicuspid and tricuspid aortic valves. In dimension reduction and clustering analyses, diseased valves tended to cluster with fetal valves at week 9 rather than normal adult valves, suggesting that part of the disease program might be due to reiterated developmental processes. The analysis of groups of coregulated genes revealed predominant immune-metabolic signatures, including innate and adaptive immune responses involving lymphocyte T-cell metabolic adaptation. Cytokine and chemokine signaling, cell migration, and proliferation were all increased in CAVD, whereas oxidative phosphorylation and protein translation were decreased. Discrete immune-metabolic gene signatures were present at fetal stages and increased in adult controls, suggesting that these processes intensify throughout life and heighten in disease. Cellular stress response and neurodegeneration gene signatures were aberrantly expressed in CAVD, pointing to a mechanistic link between chronic inflammation and biological aging. Comparison of the valve RNA-sequencing data set with a case-control study of whole blood transcriptomes from asymptomatic individuals with early aortic valve calcification identified a highly predictive gene signature of CAVD and of moderate aortic valve calcification in overtly healthy individuals. These data deepen and broaden our understanding of the molecular basis of CAVD and identify a peripheral blood gene signature for the early detection of aortic valve calcification.

Genome-wide characterization and expression analysis of geranylgeranyl diphosphate synthase genes in cotton (Gossypium spp.) in plant development and abiotic stresses

Background

GGPP (geranylgeranyl diphosphate) is produced in the isoprenoid pathway and mediates the function of various plant metabolites, which is synthesized by GGPPS (GGPP synthases) in plants. GGPPS characterization has not been performed in any plant species except Arabidopsis thaliana. Here, we performed a complete computational and bioinformatics analysis of GGPPS and detected their transcription expression pattern in Gossypium hirsutum for the first time so that to explore their evolutionary relationship and potential functions. Finally, we unravelled evolutionary relationship, conserved sequence logos, gene duplication and potential involvement in plant development and abiotic stresses tolerance of GGPPS genes in G. hirsutum and other plant species.

Results

A total of 159 GGPPS genes from 18 plant species were identified and evolutionary analysis divided these GGPPS genes into five groups to indicate their divergence from a common ancestor. Further, GGPPS family genes were conserved during evolution and underwent segmental duplication. The identified 25 GhGGPPS genes showed diverse expression pattern particularly in ovule and fiber development indicating their vital and divers roles in the fiber development. Additionally, GhGGPPS genes exhibited wide range of responses when subjected to abiotic (heat, cold, NaCl and PEG) stresses and hormonal (BL, GA, IAA, SA and MeJA) treatments, indicating their potential roles in various biotic and abiotic stresses tolerance.

Conclusions

The GGPPS genes are evolutionary conserved and might be involve in different developmental stages and stress response. Some potential key genes (e.g. GhGGPP4, GhGGPP9, and GhGGPP15) were suggested for further study and provided valuable source for cotton breeding to improve fiber quality and resistant to various stresses.

Upregulation of LAGE3 correlates with prognosis and immune infiltrates in colorectal cancer: A bioinformatic analysis

BACKGROUND

Colorectal cancer (CRC) is the primary cause of cancer-related deaths worldwide. Identification of new CRC biomarkers is imperative to improve the prognosis and development of therapies against the disease. LAGE3 (L Antigen Family Member 3) functions as a tRNA modifier, although its potential role in CRC has not been fully elucidated.

METHODS

RNA-seq matrix and corresponding clinical information were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, then subjected to survival, enrichment, and tumor microenvironment analyses using packages implemented in R.

RESULTS

We found that LAGE3 was upregulated and significantly correlating with poor prognosis in multiple CRC cohorts. Additionally, multivariate Cox regression analysis revealed that LAGE3 was an independent prognostic factor in patients with CRC, whereas functional enrichment analysis indicated that it could regulate protein targeting, tRNA processing, and the PD-1/PD-L1 checkpoint pathway. Furthermore, CIBERSORT analysis indicated a negative relationship between LAGE3 and levels of infiltration for multiple immune cells, especially CD8 + T cells in CRC. Particularly, LAGE3 expression was inversely correlated with the expression of immune checkpoints as well as that of various immune cell types of signature genes.

CONCLUSION

Collectively, our results indicate that high LAGE3 expression correlates with adverse prognosis and poor immune infiltration in CRC patients.

Genome-Wide Identification and Expression Profile Analysis of WRKY Family Genes in the Autopolyploid Saccharum spontaneum

WRKY is one of the largest transcription factor families in plants and plays important roles in the regulation of developmental and physiological processes. To date, the WRKY gene family has not been identified in Saccharum species because of its complex polyploid genome. In this study, a total of 294 sequences for 154 SsWRKY genes were identified in the polyploid Saccharum spontaneum genome and then named on the basis of their chromosome locations, including 13 (8.4%) genes with four alleles, 29 (18.8%) genes with three alleles and 41 (26.6%) genes with two alleles. Among them, 73.8% and 16.0% of the SsWRKY genes originated from segmental duplications and tandem duplications, respectively. The WRKY members exhibited conserved gene structures and amino acid sequences among the allelic haplotypes, which were accompanied by variations in intron sizes. Phylogenetic and collinearity analyses revealed that 27 SsWRKYs originated after the split of sorghum and Saccharum, resulting in a significantly higher number of WRKYs in sugarcane than in the proximal diploid species sorghum. The analysis of RNA-seq data revealed that SsWRKYs' expression profiles in 46 different samples including different developmental stages revealed distinct temporal and spatial patterns with 52 genes expressed in all tissues, four genes not expressed in any tissues and 21 SsWRKY genes likely to be involved in photosynthesis. The comprehensive analysis of SsWRKYs' expression will provide an important and valuable foundation for further investigation of the regulatory mechanisms of WRKYs in physiological roles in sugarcane S. spontaneum.

OsRRM, an RNA-Binding Protein, Modulates Sugar Transport in Rice (Oryza sativa L.)

Sugar allocation between vegetative and reproductive tissues is vital to plant development, and sugar transporters play fundamental roles in this process. Although several transcription factors have been identified that control their transcription levels, the way in which the expression of sugar transporter genes is controlled at the posttranscriptional level is unknown. In this study, we showed that OsRRM, an RNA-binding protein, modulates sugar allocation in tissues on the source-to-sink route. The OsRRM expression pattern partly resembles that of several sugar transporter and transcription factor genes that specifically affect sugar transporter gene expression. The messenger RNA levels of almost all of the sugar transporter genes are severely reduced in the osrrm mutant, and this alters sugar metabolism and sugar signaling, which further affects plant height, flowering time, seed size, and starch synthesis. We further showed that OsRRM binds directly to messenger RNAs encoded by sugar transporter genes and thus may stabilize their transcripts. Therefore, we have uncovered the physiological function of OsRRM, which sheds new light on sugar metabolism and sugar signaling.

The BMP ligand Pinhead together with Admp supports the robustness of embryonic patterning

Vertebrate embryonic dorsoventral axis is robustly stable in the face of variations in bone morphogenetic protein (BMP) signaling. However, the molecular mechanism behind this robustness remains uncharacterized. In this study, we show that zebrafish Pinhead, together with Admp, plays an important compensatory role in ensuring the robustness of axial patterning through fine-tuning of BMP signaling. pinhead encodes a BMP-like ligand expressed in the ventrolateral margin of the early gastrula. Transcription of pinhead and admp is under opposing regulation, where pinhead depletion results in a compensatory increase in admp transcription and vice versa, leading to normal axis formation in pinhead or admp mutants. Expression of pinhead and admp is directly repressed by the BMP/Smad pathway. When BMP signals were inhibited or excessively activated, pinhead/admp expression changed accordingly, allowing for self-regulation. Thus, this study reveals a negative feedback loop between BMP signaling and pinhead/admp that effectively stabilizes embryonic patterning by buffering against fluctuations in BMP signaling.

Sequence environment of BMP-dependent activating elements controls transcriptional responses to Dpp signaling in Drosophila

Intercellular signaling pathways activate transcription factors, which, along with tissue-specific co-factors, regulate expression of target genes. Responses to TGFβ/BMP signals are mediated by Smad proteins, which form complexes and accumulate in the nucleus to directly bind and regulate enhancers of BMP targets upon signaling. In Drosophila, gene activation by BMP signaling often requires, in addition to direct input by Smads, the signal-dependent removal of the transcriptional repressor Brk. Previous studies on enhancers of BMP-activated genes have defined a BMP-responsive motif, the AE, which integrates activatory and repressive input by the Smad complex and Brk, respectively. Here, we address whether sequence variations within the core AE sequences might endow the motif with additional properties accounting for qualitative and quantitative differences in BMP responses, including tissue specificity of transcriptional activation and differential sensitivity to Smad and Brk inputs. By analyzing and cross-comparing three distinct BMP-responsive enhancers from the genes wit and D ad in two different epithelia, the wing imaginal disc and the follicular epithelium, we demonstrate that differences in the AEs contribute neither to the observed tissue-restriction of BMP responses nor to differences in the utilization of the Smad and Brk branches for transcriptional activation. Rather, our results suggest that the cis-environment of the BMP-response elements not only dictates tissue specificity but also differential sensitivity to the two BMP mediators.

LCE: an open web portal to explore gene expression and clinical associations in lung cancer

We constructed a lung cancer-specific database housing expression data and clinical data from over 6700 patients in 56 studies. Expression data from 23 genome-wide platforms were carefully processed and quality controlled, whereas clinical data were standardized and rigorously curated. Empowered by this lung cancer database, we created an open access web resource—the Lung Cancer Explorer (LCE), which enables researchers and clinicians to explore these data and perform analyses. Users can perform meta-analyses on LCE to gain a quick overview of the results on tumor vs non-malignant tissue (normal) differential gene expression and expression-survival association. Individual dataset-based survival analysis, comparative analysis, and correlation analysis are also provided with flexible options to allow for customized analyses from the user.

New insights into the evolution and functional divergence of the SWEET family in Saccharum based on comparative genomics

BACKGROUND

The SWEET (Sugars Will Eventually be Exported Transporters) gene family is a recently identified group of sugar transporters that play an indispensable role in sugar efflux, phloem loading, plant-pathogen interaction, nectar secretion, and reproductive tissue development. However, little information on Saccharum SWEET is available for this crop with a complex genetic background.

RESULTS

In this study, 22 SWEET genes were identified from Saccharum spontaneum Bacterial Artificial Chromosome libraries sequences. Phylogenetic analyses of SWEETs from 11 representative plant species showed that gene expansions of the SWEET family were mainly caused by the recent gene duplication in dicot plants, while these gene expansions were attributed to the ancient whole genome duplication (WGD) in monocot plant species. Gene expression profiles were obtained from RNA-seq analysis. SWEET1a and SWEET2s had higher expression levels in the transitional zone and maturing zone than in the other analyzed zones. SWEET1b was mainly expressed in the leaf tissues and the mature zone of the leaf of both S. spontaneum and S. officinarum, and displayed a peak in the morning and was undetectable in both sclerenchyma and parenchyma cells from the mature stalks of S. officinarum. SsSWEET4a\4b had higher expression levels than SWEET4c and were mainly expressed in the stems of seedlings and mature plants. SWEET13s are recently duplicated genes, and the expression of SWEET13s dramatically increased from the maturing to mature zones. SWEET16b's expression was not detected in S. officinarum, but displayed a rhythmic diurnal expression pattern.

CONCLUSIONS

Our study revealed the gene evolutionary history of SWEETs in Saccharum and SWEET1b was found to be a sucrose starvation-induced gene involved in the sugar transportation in the high photosynthetic zones. SWEET13c was identified as the key player in the efflux of sugar transportation in mature photosynthetic tissues. SWEET4a\4b were found to be mainly involved in sugar transportation in the stalk. SWEET1a\2a\4a\4b\13a\16b were suggested to be the genes contributing to the differences in sugar contents between S. spontaneum and S. officinarum. Our results are valuable for further functional analysis of SWEET genes and utilization of the SWEET genes for genetic improvement of Saccharum for biofuel production.

Moonlighting Phosphatase Activity of Klenow DNA Polymerase in the Presence of RNA

RNA is a key player in the cellular central dogma, including RNA transcription and protein synthesis. However, it is unknown whether RNA can directly interfere with DNA synthesis. Recently, we have found in vitro that while binding to DNA polymerase nonspecifically, RNA can transform DNA polymerase to display a moonlighting activity, dNTP phosphatase, in turn interfering with DNA synthesis. This phosphatase activity removes the γ-phosphate from dNTPs (generating dNDPs) and subsequently removes the β-phosphate from the formed dNDPs (generating dNMPs), confirmed by the noncleavable α,β-CH₂-dGTP and β,γ-CH₂-dGTP analogues. We also found that dGTP is the best substrate for the phosphatase, and the dNTP phosphatase activity is sensitive to the reaction medium. In addition, we have revealed that RNA can tune the activity of closely related proteins and give rise to new catalytic functions with subtle differences. Moreover, we have demonstrated in vitro that at the lower dNTP level, this phosphatase can directly inhibit DNA synthesis by dNTP depletion, though the phosphatase activity is 690-fold slower than the polymerase activity. Our observation in vitro suggests a plausible strategy for RNA to directly interfere with DNA polymerase and DNA synthesis in vivo.

NvERTx: a gene expression database to compare embryogenesis and regeneration in the sea anemone Nematostella vectensis

For over a century, researchers have been comparing embryogenesis and regeneration hoping that lessons learned from embryonic development will unlock hidden regenerative potential. This problem has historically been a difficult one to investigate because the best regenerative model systems are poor embryonic models and vice versa. Recently, however, there has been renewed interest in this question, as emerging models have allowed researchers to investigate these processes in the same organism. This interest has been further fueled by the advent of high-throughput transcriptomic analyses that provide virtual mountains of data. Here, we present Nematostella vectensis Embryogenesis and Regeneration Transcriptomics (NvERTx), a platform for comparing gene expression during embryogenesis and regeneration. NvERTx consists of close to 50 transcriptomic data sets spanning embryogenesis and regeneration in Nematostella These data were used to perform a robust de novo transcriptome assembly, with which users can search, conduct BLAST analyses, and plot the expression of multiple genes during these two developmental processes. The site is also home to the results of gene clustering analyses, to further mine the data and identify groups of co-expressed genes. The site can be accessed at http://nvertx.kahikai.org.

The E3 ubiquitin ligase UBR5 regulates centriolar satellite stability and primary cilia

Primary cilia are crucial for signal transduction in a variety of pathways, including hedgehog and Wnt. Disruption of primary cilia formation (ciliogenesis) is linked to numerous developmental disorders (known as ciliopathies) and diseases, including cancer. The ubiquitin-proteasome system (UPS) component UBR5 was previously identified as a putative positive regulator of ciliogenesis in a functional genomics screen. UBR5 is an E3 ubiquitin ligase that is frequently deregulated in tumors, but its biological role in cancer is largely uncharacterized, partly due to a lack of understanding of interacting proteins and pathways. We validated the effect of UBR5 depletion on primary cilia formation using a robust model of ciliogenesis, and identified CSPP1, a centrosomal and ciliary protein required for cilia formation, as a UBR5-interacting protein. We show that UBR5 ubiquitylates CSPP1, and that UBR5 is required for cytoplasmic organization of CSPP1-comprising centriolar satellites in centrosomal periphery, suggesting that UBR5-mediated ubiquitylation of CSPP1 or associated centriolar satellite constituents is one underlying requirement for cilia expression. Hence, we have established a key role for UBR5 in ciliogenesis that may have important implications in understanding cancer pathophysiology.

FSPP: A Tool for Genome-Wide Prediction of smORF-Encoded Peptides and Their Functions

smORFs are small open reading frames of less than 100 codons. Recent low throughput experiments showed a lot of smORF-encoded peptides (SEPs) played crucial rule in processes such as regulation of transcription or translation, transportation through membranes and the antimicrobial activity. In order to gather more functional SEPs, it is necessary to have access to genome-wide prediction tools to give profound directions for low throughput experiments. In this study, we put forward a functional smORF-encoded peptides predictor (FSPP) which tended to predict authentic SEPs and their functions in a high throughput method. FSPP used the overlap of detected SEPs from Ribo-seq and mass spectrometry as target objects. With the expression data on transcription and translation levels, FSPP built two co-expression networks. Combing co-location relations, FSPP constructed a compound network and then annotated SEPs with functions of adjacent nodes. Tested on 38 sequenced samples of 5 human cell lines, FSPP successfully predicted 856 out of 960 annotated proteins. Interestingly, FSPP also highlighted 568 functional SEPs from these samples. After comparison, the roles predicted by FSPP were consistent with known functions. These results suggest that FSPP is a reliable tool for the identification of functional small peptides. FSPP source code can be acquired at https://www.bioinfo.org/FSPP.

Identification of Isthmin 1 as a Novel Clefting and Craniofacial Patterning Gene in Humans

Orofacial clefts are one of the most common birth defects, affecting 1-2 per 1000 births, and have a complex etiology. High-resolution array-based comparative genomic hybridization has increased the ability to detect copy number variants (CNVs) that can be causative for complex diseases such as cleft lip and/or palate. Utilizing this technique on 97 nonsyndromic cleft lip and palate cases and 43 cases with cleft palate only, we identified a heterozygous deletion of Isthmin 1 in one affected case, as well as a deletion in a second case that removes putative 3' regulatory information. Isthmin 1 is a strong candidate for clefting, as it is expressed in orofacial structures derived from the first branchial arch and is also in the same "synexpression group" as fibroblast growth factor 8 and sprouty RTK signaling antagonist 1a and 2, all of which have been associated with clefting. CNVs affecting Isthmin 1 are exceedingly rare in control populations, and Isthmin 1 scores as a likely haploinsufficiency locus. Confirming its role in craniofacial development, knockdown or clustered randomly interspaced short palindromic repeats/Cas9-generated mutation of isthmin 1 in Xenopus laevis resulted in mild to severe craniofacial dysmorphologies, with several individuals presenting with median clefts. Moreover, knockdown of isthmin 1 produced decreased expression of LIM homeobox 8, itself a gene associated with clefting, in regions of the face that pattern the maxilla. Our study demonstrates a successful pipeline from CNV identification of a candidate gene to functional validation in a vertebrate model system, and reveals Isthmin 1 as both a new human clefting locus as well as a key craniofacial patterning gene.

Genomic regression analysis of coordinated expression

Co-expression analysis is widely used to predict gene function and to identify functionally related gene sets. However, co-expression analysis using human cancer transcriptomic data is confounded by somatic copy number alterations (SCNA), which produce co-expression signatures based on physical proximity rather than biological function. To better understand gene-gene co-expression based on biological regulation but not SCNA, we describe a method termed "Genomic Regression Analysis of Coordinated Expression" (GRACE) to adjust for the effect of SCNA in co-expression analysis. The results from analyses of TCGA, CCLE, and NCI60 data sets show that GRACE can improve our understanding of how a transcriptional network is re-wired in cancer. A user-friendly web database populated with data sets from The Cancer Genome Atlas (TCGA) is provided to allow customized query.

TGF-β Family Signaling in Tumor Suppression and Cancer Progression

Transforming growth factor-β (TGF-β) induces a pleiotropic pathway that is modulated by the cellular context and its integration with other signaling pathways. In cancer, the pleiotropic reaction to TGF-β leads to a diverse and varied set of gene responses that range from cytostatic and apoptotic tumor-suppressive ones in early stage tumors, to proliferative, invasive, angiogenic, and oncogenic ones in advanced cancer. Here, we review the knowledge accumulated about the molecular mechanisms involved in the dual response to TGF-β in cancer, and how tumor cells evolve to evade the tumor-suppressive responses of this signaling pathway and then hijack the signal, converting it into an oncogenic factor. Only through the detailed study of this complexity can the suitability of the TGF-β pathway as a therapeutic target against cancer be evaluated.

A molecular atlas of the developing ectoderm defines neural, neural crest, placode, and nonneural progenitor identity in vertebrates

During vertebrate neurulation, the embryonic ectoderm is patterned into lineage progenitors for neural plate, neural crest, placodes and epidermis. Here, we use Xenopus laevis embryos to analyze the spatial and temporal transcriptome of distinct ectodermal domains in the course of neurulation, during the establishment of cell lineages. In order to define the transcriptome of small groups of cells from a single germ layer and to retain spatial information, dorsal and ventral ectoderm was subdivided along the anterior-posterior and medial-lateral axes by microdissections. Principal component analysis on the transcriptomes of these ectoderm fragments primarily identifies embryonic axes and temporal dynamics. This provides a genetic code to define positional information of any ectoderm sample along the anterior-posterior and dorsal-ventral axes directly from its transcriptome. In parallel, we use nonnegative matrix factorization to predict enhanced gene expression maps onto early and mid-neurula embryos, and specific signatures for each ectoderm area. The clustering of spatial and temporal datasets allowed detection of multiple biologically relevant groups (e.g., Wnt signaling, neural crest development, sensory placode specification, ciliogenesis, germ layer specification). We provide an interactive network interface, EctoMap, for exploring synexpression relationships among genes expressed in the neurula, and suggest several strategies to use this comprehensive dataset to address questions in developmental biology as well as stem cell or cancer research.

Vertebrate embryo germ layers become progressively regionalized by evolutionarily conserved molecular processes. Catching the early steps of this dynamic spatial cell diversification at the scale of the transcriptome was challenging, even with the advent of efficient RNA sequencing. We have microdissected complementary and defined areas of a single germ layer, the developing ectoderm, and explored how the transcriptome changes over time and space in the ectoderm during the differentiation of frog epidermis, neural plate, and neural crest. We have created EctoMap, a searchable interface using these regional transcriptomes, to predict the expression of the 31 thousand genes expressed in neurulae and their networks of co-expression, predictive of functional relationships. Through several examples, we illustrate how these data provide insights in development, cancer, evolution and stem cell biology.

Diurnal Cycling Transcription Factors of Pineapple Revealed by Genome-Wide Annotation and Global Transcriptomic Analysis

Circadian clock provides fitness advantage by coordinating internal metabolic and physiological processes to external cyclic environments. Core clock components exhibit daily rhythmic changes in gene expression, and the majority of them are transcription factors (TFs) and transcription coregulators (TCs). We annotated 1,398 TFs from 67 TF families and 80 TCs from 20 TC families in pineapple, and analyzed their tissue-specific and diurnal expression patterns. Approximately 42% of TFs and 45% of TCs displayed diel rhythmic expression, including 177 TF/TCs cycling only in the nonphotosynthetic leaf tissue, 247 cycling only in the photosynthetic leaf tissue, and 201 cycling in both. We identified 68 TF/TCs whose cycling expression was tightly coupled between the photosynthetic and nonphotosynthetic leaf tissues. These TF/TCs likely coordinate key biological processes in pineapple as we demonstrated that this group is enriched in homologous genes that form the core circadian clock in Arabidopsis and includes a STOP1 homolog. Two lines of evidence support the important role of the STOP1 homolog in regulating CAM photosynthesis in pineapple. First, STOP1 responds to acidic pH and regulates a malate channel in multiple plant species. Second, the cycling expression pattern of the pineapple STOP1 and the diurnal pattern of malate accumulation in pineapple leaf are correlated. We further examined duplicate-gene retention and loss in major known circadian genes and refined their evolutionary relationships between pineapple and other plants. Significant variations in duplicate-gene retention and loss were observed for most clock genes in both monocots and dicots.

Generation of animal form by the Chordin/Tolloid/BMP gradient: 100 years after D'Arcy Thompson

The classic book "On Growth and Form" by naturalist D'Arcy Thompson was published 100 years ago. To celebrate this landmark, we present experiments in the Xenopus embryo that provide a framework for understanding how simple, quantitative transformations of a morphogen gradient might have affected evolution and morphological diversity of organisms. D'Arcy Thompson proposed that different morphologies might be generated by modifying physical parameters in an underlying system of Cartesian coordinates that pre-existed in Nature and arose during evolutionary history. Chordin is a BMP antagonist secreted by the Spemann organizer located on the dorsal side of the gastrula. Chordin generates a morphogen gradient as first proposed by mathematician Alan Turing. The rate-limiting step of this dorsal-ventral (D-V) morphogen is the degradation of Chordin by the Tolloid metalloproteinase in the ventral side. Chordin is expressed at gastrula on the dorsal side where BMP signaling is low, while at the opposite side peak levels of BMP signaling are reached. In fishes, amphibians, reptiles and birds, high BMP signaling in the ventral region induces transcription of a secreted inhibitor of Tolloid called Sizzled. By depleting Sizzled exclusively in the ventral half of the embryo we were able to expand the ventro-posterior region in an otherwise normal embryo. Conversely, ventral depletion of Tolloid, which stabilizes Chordin, decreased ventral and tail structures, phenocopying the tolloid zebrafish mutation. We explain how historical constraints recorded in the language of DNA become subject to the universal laws of physics when an ancestral reaction-diffusion morphogen gradient dictates form.

Bud detachment in hydra requires activation of fibroblast growth factor receptor and a Rho-ROCK-myosin II signaling pathway to ensure formation of a basal constriction

Background:Hydra propagates asexually by exporting tissue into a bud, which detaches 4 days later as a fully differentiated young polyp. Prerequisite for detachment is activation of fibroblast growth factor receptor (FGFR) signaling. The mechanism which enables constriction and tissue separation within the monolayered ecto‐ and endodermal epithelia is unknown. Results: Histological sections and staining of F‐actin by phalloidin revealed conspicuous cell shape changes at the bud detachment site indicating a localized generation of mechanical forces and the potential enhancement of secretory functions in ectodermal cells. By gene expression analysis and pharmacological inhibition, we identified a candidate signaling pathway through Rho, ROCK, and myosin II, which controls bud base constriction and rearrangement of the actin cytoskeleton. Specific regional myosin phosphorylation suggests a crucial role of ectodermal cells at the detachment site. Inhibition of FGFR, Rho, ROCK, or myosin II kinase activity is permissive for budding, but represses myosin phosphorylation, rearrangement of F‐actin and constriction. The young polyp remains permanently connected to the parent by a broad tissue bridge. Conclusions: Our data suggest an essential role of FGFR and a Rho‐ROCK‐myosin II pathway in the control of cell shape changes required for bud detachment. Developmental Dynamics 246:502–516, 2017. © 2017 The Authors Developmental Dynamics published by Wiley Periodicals, Inc. on behalf of American Association of Anatomists

Hydra bud detachment involves the separation of two intact epithelia without cell death.

Remarkable cell shape changes and multicellular rosettes at the bud base indicate functional specification and strong mechanical forces.

mRNA colocalization, phospho‐myosin analysis and similar phenotypes obtained by pharmacological inhibition suggest a tight correlation between FGFR and a Rho‐ROCK‐Myosin II candidate signaling pathway.

Human fetal heart specific coexpression network involves congenital heart disease/defect candidate genes

Heart development is a complex process requiring dynamic transcriptional regulation. Disturbance of this process will lead to severe developmental defects such as congenital heart disease/defect (CHD). CHD is a group of complex disorder with high genetic heterogeneity, common pathways associated with CHD remains largely unknown. In the manuscript, we focused on the tissue specific genes in human fetal heart samples to explore such pathways. We used the RNA microarray dataset of human fetal tissues from ENCODE project to identify genes with heart tissue specific expression. A transcriptional network was constructed for these genes based on the Pearson correlation coefficients of their expression levels. Function, selective constraints and disease associations of these genes were then examined. Our analysis identified a network consisted of 316 genes with human fetal heart specific expression. The network was highly co-regulated and showed evolutionary conserved tissue expression pattern in tetrapod. Genes in this network are enriched in CHD specific genes and disease mutations. Using the transcriptomic data, we discovered a highly concerted gene network that might reflect a common pathway associated with the etiology of CHD. Such analysis should be helpful for disease associated gene identification in clinical studies.