Growth arrest-specific 2 protein family: Structure and function

A novel variant in GAS2 is associated with autosomal dominant nonsyndromic hearing impairment in a Chinese family

Knockout of GAS2 (growth arrest-specific protein 2), causes disorganization and destabilization of microtubule bundles in supporting cells of the cochlear duct, leading to hearing loss in vivo. However, the molecular mechanism through which GAS2 variant results in hearing loss remains unknown. By Whole-exome sequencing, we identified a novel heterozygous splicing variant in GAS2 (c.616-2 A > G) as the only candidate mutation segregating with late-onset and progressive nonsyndromic hearing loss (NSHL) in a large dominant family. This splicing mutation causes an intron retention and produces a C-terminal truncated protein (named GAS2mu). Mechanistically, the degradation of GAS2mu via the ubiquitin-proteasome pathway is enhanced, and cells expressing GAS2mu exhibit disorganized microtubule bundles. Additionally, GAS2mu further promotes apoptosis by increasing the Bcl-xS/Bcl-xL ratio instead of through the p53-dependent pathway as wild-type GAS2 does, indicating that GAS2mu acts as a toxic molecule to exacerbate apoptosis. Our findings demonstrate that this novel variant of GAS2 promotes its own protein degradation, microtubule disorganization and cellular apoptosis, leading to hearing loss in carriers. This study expands the spectrum of GAS2 variants and elucidates the underlying pathogenic mechanisms, providing a foundation for future investigations of new therapeutic strategies to prevent GAS2-associated progressive hearing loss.

A comprehensive analysis of GAS2 family members identifies that GAS2L1 is a novel biomarker and promotes the proliferation of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a common primary liver cancer with a high incidence and mortality. Members of the growth-arresting-specific 2 (GAS2) family are involved in various biological processes in human malignancies. To date, there is only a limited amount of information available about the expression profile and clinical importance of GAS2 family in HCC. In this study, we found that GAS2L1 and GAS2L3 were distinctly upregulated in HCC specimens compared to non-tumor specimens. Pan-cancer assays indicated that GAS2L1 and GAS2L3 were highly expressed in most cancers. The Pearson's correlation revealed that the expressions of GAS2, GAS2L1 and GAS2L2 were negatively associated with methylation levels. Survival assays indicated that GAS2L1 and GAS2L3 were independent prognostic factors for HCC patients. Immune cell infiltration analysis revealed that GAS2, GAS2L1 and GAS2L3 were associated with several immune cells. Finally, we confirmed that GAS2L1 was highly expressed in HCC cells and its knockdown suppressed the proliferation of HCC cells. Taken together, our findings suggested the expression patterns and prognostic values of GAS2 members in HCC, providing insights for further study of the GAS2 family as sensitive diagnostic and prognostic markers for HCC.

1-L Transcription of SARS-CoV-2 Spike Protein S1 Subunit

The COVID-19 pandemic prompted rapid research on SARS-CoV-2 pathogenicity. Consequently, new data can be used to advance the molecular understanding of SARS-CoV-2 infection. The present bioinformatics study discusses the "spikeopathy" at the molecular level and focuses on the possible post-transcriptional regulation of the SARS-CoV-2 spike protein S1 subunit in the host cell/tissue. A theoretical protein-RNA recognition code was used to check the compatibility of the SARS-CoV-2 spike protein S1 subunit with mRNAs in the human transcriptome (1-L transcription). The principle for this method is elucidated on the defined RNA binding protein GEMIN5 (gem nuclear organelle-associated protein 5) and RNU2-1 (U2 spliceosomal RNA). Using the method described here, it was shown that 45% of the genes/proteins identified by 1-L transcription of the SARS-CoV-2 spike protein S1 subunit are directly linked to COVID-19, 39% are indirectly linked to COVID-19, and 16% cannot currently be associated with COVID-19. The identified genes/proteins are associated with stroke, diabetes, and cardiac injury.

Development of the RF-GSEA Method for Identifying Disulfidptosis-Related Genes and Application in Hepatocellular Carcinoma

Disulfidptosis is a newly discovered cellular programmed cell death mode. Presently, a considerable number of genes related to disulfidptosis remain undiscovered, and its significance in hepatocellular carcinoma remains unrevealed. We have developed a powerful analytical method called RF-GSEA for identifying potential genes associated with disulfidptosis. This method draws inspiration from gene regulation networks and graph theory, and it is implemented through a combination of random forest regression model and Gene Set Enrichment Analysis. Subsequently, to validate the practical application value of this method, we applied it to hepatocellular carcinoma. Based on the RF-GSEA method, we developed a disulfidptosis-related signature. Lastly, we looked into how the disulfidptosis-related signature is connected to HCC prognosis, the tumor microenvironment, the effectiveness of immunotherapy, and the sensitivity of chemotherapy drugs. The RF-GSEA method identified a total of 220 disulfidptosis-related genes, from which 7 were selected to construct the disulfidptosis-related signature. The high-disulfidptosis-related score group had a worse prognosis compared to the low-disulfidptosis-related score group and showed lower infiltration levels of immune-promoting cells. The high-disulfidptosis-related score group had a higher likelihood of benefiting from immunotherapy compared to the low-disulfidptosis-related score group. The RF-GSEA method is a powerful tool for identifying disulfidptosis-related genes. The disulfidptosis-related signature effectively predicts HCC prognosis, immunotherapy response, and drug sensitivity.

Transcriptome RNA Sequencing Reveals That Circular RNAs Are Abundantly Expressed in Embryonic Breast Muscle of Duck

Circular RNAs are widespread in various species and have important roles in myogenesis. However, the circular RNAs involved in breast muscle development in ducks have not yet been studied. Here, to identify circular RNAs during duck skeletal muscle development, three pectorales from Shan Ma ducks at E13 and E19, which represent undifferentiated and differentiated myoblasts, respectively, were collected and subjected to RNA sequencing. A total of 16,622 circular RNAs were identified, of which approximately 80% were exonic circular RNAs and 260 were markedly differentially expressed between E19 and E13. The parental genes of the differentially expressed circular RNAs were significantly enriched in muscle-related biological processes. Moreover, we found that the overexpression of circGAS2-2 promoted cell cycle progression and increased the proliferation viability of duck primary myoblasts; conversely, knockdown of circGAS2-2 retarded the cell cycle and reduced the proliferation viability of myoblasts. Taken together, our results demonstrate that circular RNAs are widespread and variously expressed during the development of duck skeletal muscle and that circGAS2-2 is involved in the regulation of myogenesis.

Genome-wide association study in a rat model of temperament identifies multiple loci for exploratory locomotion and anxiety-like traits

Common genetic factors likely contribute to multiple psychiatric diseases including mood and substance use disorders. Certain stable, heritable traits reflecting temperament, termed externalizing or internalizing, play a large role in modulating vulnerability to these disorders. To model these heritable tendencies, we selectively bred rats for high and low exploration in a novel environment [bred High Responders (bHR) vs. Low Responders (bLR)]. To identify genes underlying the response to selection, we phenotyped and genotyped 538 rats from an F₂ cross between bHR and bLR. Several behavioral traits show high heritability, including the selection trait: exploratory locomotion (EL) in a novel environment. There were significant phenotypic and genetic correlations between tests that capture facets of EL and anxiety. There were also correlations with Pavlovian conditioned approach (PavCA) behavior despite the lower heritability of that trait. Ten significant and conditionally independent loci for six behavioral traits were identified. Five of the six traits reflect different facets of EL that were captured by three behavioral tests. Distance traveled measures from the open field and the elevated plus maze map onto different loci, thus may represent different aspects of novelty-induced locomotor activity. The sixth behavioral trait, number of fecal boli, is the only anxiety-related trait mapping to a significant locus on chromosome 18 within which the Pik3c3 gene is located. There were no significant loci for PavCA. We identified a missense variant in the Plekhf1 gene on the chromosome 1:95 Mb QTL and Fancf and Gas2 as potential candidate genes that may drive the chromosome 1:107 Mb QTL for EL traits. The identification of a locomotor activity-related QTL on chromosome 7 encompassing the Pkhd1l1 and Trhr genes is consistent with our previous finding of these genes being differentially expressed in the hippocampus of bHR vs. bLR rats. The strong heritability coupled with identification of several loci associated with exploratory locomotion and emotionality provide compelling support for this selectively bred rat model in discovering relatively large effect causal variants tied to elements of internalizing and externalizing behaviors inherent to psychiatric and substance use disorders.

Genome-wide transcriptional profiling of pulmonary functional sequelae in ARDS- secondary to SARS-CoV-2 infection

BACKGROUND

Up to 80% of patients surviving acute respiratory distress syndrome (ARDS) secondary to SARS-CoV-2 infection present persistent anomalies in pulmonary function after hospital discharge. There is a limited understanding of the mechanistic pathways linked to post-acute pulmonary sequelae.

AIM

To identify the molecular underpinnings associated with severe lung diffusion involvement in survivors of SARS-CoV-2-induced ARDS.

METHODS

Survivors attended to a complete pulmonary evaluation 3 months after hospital discharge. RNA sequencing (RNA-seq) was performed using Illumina technology in whole-blood samples from 50 patients with moderate to severe diffusion impairment (D_LCO<60%) and age- and sex-matched individuals with mild-normal lung function (D_LCO≥60%). A transcriptomic signature for optimal classification was constructed using random forest. Transcriptomic data were analyzed for biological pathway enrichment, cellular deconvolution, cell/tissue-specific gene expression and candidate drugs.

RESULTS

RNA-seq identified 1357 differentially expressed transcripts. A model composed of 14 mRNAs allowed the optimal discrimination of survivors with severe diffusion impairment (AUC=0.979). Hallmarks of lung sequelae involved cell death signaling, cytoskeleton reorganization, cell growth and differentiation and the immune response. Resting natural killer (NK) cells were the most important immune cell subtype for the prediction of severe diffusion impairment. Components of the signature correlated with neutrophil, lymphocyte and monocyte counts. A variable expression profile of the transcripts was observed in lung cell subtypes and bodily tissues. One upregulated gene, TUBB4A, constitutes a target for FDA-approved drugs.

CONCLUSIONS

This work defines the transcriptional programme associated with post-acute pulmonary sequelae and provides novel insights for targeted interventions and biomarker development.

Growth arrest‐specific protein 2 ( GAS2 ) interacts with CXCR4 to promote T‐cell leukemogenesis partially via c‐MYC

Although growth arrest‐specific protein 2 (GAS2) promotes the growth of T‐cell acute lymphoblastic leukemia (T‐ALL) cells in culture, the effect of GAS2 on T‐cell leukemogenesis has not been studied, and the mechanism remains unclear. In the present study, xenograft studies showed that GAS2 silencing impaired T‐cell leukemogenesis and decreased leukemic cell infiltration. Mechanistically, GAS2 regulated the protein expression of C‐X‐C chemokine receptor type 4 (CXCR4) rather than its transcript expression. Immunoprecipitation revealed that GAS2 interacted with CXCR4, and confocal analysis showed that GAS2 was partially co‐expressed with CXCR4, which provided a strong molecular basis for GAS2 to regulate CXCR4 expression. Importantly, CXCR4 overexpression alleviated the inhibitory effect of GAS2 silencing on the growth and migration of T‐ALL cells. Moreover, GAS2 or CXCR4 silencing inhibited the expression of NOTCH1 and c‐MYC. Forced expression of c‐MYC rescued the growth suppression induced by GAS2 or CXCR4 silencing. Meanwhile, GAS2 deficiency, specifically in blood cells, had a mild effect on normal hematopoiesis, including T‐cell development, and GAS2 silencing did not affect the growth of normal human CD3⁺ or CD34⁺ cells. Overall, our data indicate that GAS2 promotes T‐cell leukemogenesis through its interaction with CXCR4 to activate NOTCH1/c‐MYC, whereas impaired GAS2 expression has a mild effect on normal hematopoiesis. Therefore, our study suggests that targeting the GAS2/CXCR4 axis is a potential therapeutic strategy for T‐ALL.

Determining protein polarization proteome-wide using physical dissection of individual Stentor coeruleus cells

Cellular components are non-randomly arranged with respect to the shape and polarity of the whole cell.^1-4 Patterning within cells can extend down to the level of individual proteins and mRNA.^5,6 But how much of the proteome is actually localized with respect to cell polarity axes? Proteomics combined with cellular fractionation^7-11 has shown that most proteins localize to one or more organelles but does not tell us how many proteins have a polarized localization with respect to the large-scale polarity axes of the intact cell. Genome-wide localization studies in yeast^12-15 found that only a few percent of proteins have a localized position relative to the cell polarity axis defined by sites of polarized cell growth. Here, we describe an approach for analyzing protein distribution within a cell with a visibly obvious global patterning-the giant ciliate Stentor coeruleus.^16,17 Ciliates, including Stentor, have highly polarized cell shapes with visible surface patterning.^1,18 A Stentor cell is roughly 2 mm long, allowing a "proteomic dissection" in which microsurgery is used to separate cellular fragments along the anterior-posterior axis, followed by comparative proteomic analysis. In our analysis, 25% of the proteome, including signaling proteins, centrin/SFI proteins, and GAS2 orthologs, shows a polarized location along the cell's anterior-posterior axis. We conclude that a large proportion of all proteins are polarized with respect to global cell polarity axes and that proteomic dissection provides a simple and effective approach for spatial proteomics.

The Full-Length Transcriptome Provides New Insights Into the Transcript Complexity of Abdominal Adipose and Subcutaneous Adipose in Pekin Ducks

Adipose tissues have a central role in organisms, and adipose content is a crucial economic trait of poultry. Pekin duck is an ideal model to study the mechanism of abdominal and subcutaneous adipose deposition for its high ability of adipose synthesis and deposition. Alternative splicing contributes to functional diversity in abdominal and subcutaneous adipose. However, there has been no systematic analysis of the dynamics of differential alternative splicing of abdominal and subcutaneous adipose in Pekin duck. In our study, the Pacific Biosciences (PacBio) Iso-Seq technology was applied to explore the transcriptional complexity of abdominal and subcutaneous adipose in Pekin ducks. In total, 143,931 and 111,337 full-length non-chimeric transcriptome sequences of abdominal and subcutaneous adipocytes were obtained from 41.78 GB raw data, respectively. These data led us to identify 19,212 long non-coding RNAs (lncRNAs) and 74,571 alternative splicing events. In addition, combined with the next-generation sequencing technology, we correlated the structure and function annotation with the differential expression profiles of abdominal and subcutaneous adipose transcripts. This study identified lots of novel alternative splicing events and major transcripts of transcription factors related to adipose synthesis. STAT3 was reported as a vital gene for adipogenesis, and we found that its major transcript is STAT3-1, which may play a considerable role in the process of adipose synthesis in Pekin duck. This study greatly increases our understanding of the gene models, genome annotations, genome structures, and the complexity and diversity of abdominal and subcutaneous adipose in Pekin duck. These data provide insights into the regulation of alternative splicing events, which form an essential part of transcript diversity during adipogenesis in poultry. The results of this study provide an invaluable resource for studying alternative splicing and tissue-specific expression.

A Novel Toxoplasma Inner Membrane Complex Suture-Associated Protein Regulates Suture Protein Targeting and Colocalizes with Membrane Trafficking Machinery

The cytoskeleton of Toxoplasma gondii is composed of the inner membrane complex (IMC) and an array of underlying microtubules that provide support at the periphery of the parasite. Specific subregions of the IMC carry out distinct roles in replication, motility, and host cell invasion. Building on our previous in vivo biotinylation (BioID) experiments of the IMC, we identified here a novel protein that localizes to discrete puncta that are embedded in the parasite’s cytoskeleton along the IMC sutures. Gene knockout analysis showed that loss of the protein results in defects in cytoskeletal suture protein targeting, cytoskeletal integrity, parasite morphology, and host cell invasion. We then used deletion analyses to identify a domain in the N terminus of the protein that is critical for both localization and function. Finally, we used the protein as bait for in vivo biotinylation, which identified several other proteins that colocalize in similar spot-like patterns. These putative interactors include several proteins that are implicated in membrane trafficking and are also associated with the cytoskeleton. Together, these data reveal an unexpected link between the IMC sutures and membrane trafficking elements of the parasite and suggest that the suture puncta are likely a portal for trafficking cargo across the IMC.

Genetic expression and mutational profile analysis in different pathologic stages of hepatocellular carcinoma patients

BACKGROUND

The clinical pathologic stages (stage I, II, III-IV) of hepatocellular carcinoma (HCC) are closely linked to the clinical prognosis of patients. This study aims at investigating the gene expression and mutational profile in different clinical pathologic stages of HCC.

METHODS

Based on the TCGA-LIHC cohort, we utilized a series of analytical approaches, such as statistical analysis, random forest, decision tree, principal component analysis (PCA), to identify the differential gene expression and mutational profiles. The expression patterns of several targeting genes were also verified by analyzing the Chinese HLivH060PG02 HCC cohort, several GEO datasets, HPA database, and diethylnitrosamine-induced HCC mouse model.

RESULTS

We identified a series of targeting genes with copy number variation, which is statistically associated with gene expression. Non-synonymous mutations mainly existed in some genes (e.g.,TTN, TP53, CTNNB1). Nevertheless, no association between gene mutation frequency and pathologic stage distribution was detected. The random forest and decision tree modeling analysis data showed a group of genes related to different HCC pathologic stages, including GAS2L3 and SEMA3F. Additionally, our PCA data indicated several genes associated with different pathologic stages, including SNRPA and SNRPD2. Compared with adjacent normal tissues, we observed a highly expressed level of GAS2L3, SNRPA, and SNRPD2 (P = 0.002) genes in HCC tissues of our HLivH060PG02 cohort. We also detected the high expression pattern of GAS2L3, SEMA3F, SNRPA, and SNRPD2 in the datasets of GSE102079, GSE76427, GSE64041, GSE121248, GSE84005, and the qPCR assay using diethylnitrosamine-induced HCC mouse model. Moreover, SEMA3F and SNRPD2 protein were highly stained in the HCC tissues of the HPA database. The high expression level of these four genes was associated with the poor survival prognosis of HCC cases.

CONCLUSIONS

Our study provides evidence regarding the gene expression and mutational profile in different clinical pathologic stages of TCGA HCC cases. Identifying four targeting genes, including GAS2L3, SNRPA, SNRPD2, and SEMA3F, offers insight into the molecular mechanisms associated with different prognoses of HCC.

Novel computational analysis of large transcriptome datasets identifies sets of genes distinguishing chronic obstructive pulmonary disease from healthy lung samples

Chronic obstructive pulmonary disease (COPD) kills over three million people worldwide every year. Despite its high global impact, the knowledge about the underlying molecular mechanisms is still limited. In this study, we aimed to extend the available knowledge by identifying a small set of COPD-associated genes. We analysed different publicly available gene expression datasets containing whole lung tissue (WLT) and airway epithelium (AE) samples from over 400 human subjects for differentially expressed genes (DEGs). We reduced the resulting sets of 436 and 663 DEGs using a novel computational approach that utilises a random depth-first search to identify genes which improve the distinction between COPD patients and controls along the first principle component of the data. Our method identified small sets of 10 and 15 genes in the WLT and AE, respectively. These sets of genes significantly (p < 10^–20) distinguish COPD patients from controls with high fidelity. The final sets revealed novel genes like cysteine rich protein 1 (CRIP1) or secretoglobin family 3A member 2 (SCGB3A2) that may underlie fundamental molecular mechanisms of COPD in these tissues.

Integrative analysis regarding the correlation between GAS2 family genes and human glioma prognosis

Background

Emerging oncogenes were reportedly linked to the complicated subtypes and pathogenesis of clinical gliomas. Herein, we first comprehensively explored the potential correlation between growth‐arrest‐specific two family genes (GAS2, GAS2L1, GAS2L2, GAS2L3) and gliomas by bioinformatics analysis and cellular experiments.

Methods

Based on the available datasets of TCGA (The Cancer Genome Atlas), CGGA (Chinese Glioma Genome Atlas), and Oncomine databases, we performed a series of analyses, such as gene expression, survival prognosis, DNA methylation, immune infiltration, and partner enrichment. We also utilized two glioma cell lines to conduct the colony formation and wound‐healing assay.

Results

GAS2L3 gene was highly expressed in glioma tissues compared to normal brain tissues (p < 0.05). We further observed the relationship between the high expressed GAS2L3 and poor clinical prognosis of brain low‐grade glioma (LGG) cases in our Cox proportional hazard model (hazard ratio [HR] = 0.1715, p < 0.001). Moreover, DNA hypomethylation status of GAS2L3 was correlated with the high expression of GAS2L3 in LGG tissues and the poor clinical prognosis of primary glioma cases (p < 0.05). We also found that the high expression of GAS2L3 was associated with the infiltration level of immune cells, especially the T cells (p < 0.0001). Functional enrichment analysis of GAS2L3‐correlated genes and interaction partners further indicated that GAS2L3 might take part in the occurrence of glioma by influencing a series of biological behaviors, such as cell division, cytoskeleton binding, and cell adhesion. Additionally, our cellular experiment data suggested that a highly expressed GAS2L3 gene contributes to the enhanced proliferation and migration of glioma cells.

Conclusion

This study first analyzed the potential role of GAS2 family genes, especially GAS2L3, in the clinical prognosis and possible functional mechanisms of glioma, which gives a novel insight into the relationship between GAS2L3 and LGG.

Growth arrest-specific 2 protein family: Structure and function

Members of the growth arrest–specific 2 (GAS2) protein family consist of a putative actin‐binding (CH) domain and a microtubule‐binding (GAR) domain and are considered miniversions of spectraplakins. There are four members in the GAS2 family, viz. GAS2, GAS2L1, GAS2L2 and GAS2L3. Although GAS2 is defined as a family of growth arrest–specific proteins, the significant differences in the expression patterns, interaction characteristics and biological issues or diseases among the different GAS2 family members have not been systemically reviewed to date. Therefore, we summarized the available evidence on the structures and functions of GAS2 family members. This review facilitates a comprehensive molecular understanding of the involvement of the GAS2 family members in an array of biological processes, including cytoskeleton reorganization, cell cycle, apoptosis and cancer development.