FAM13A as a Novel Hypoxia-Induced Gene in Non-Small Cell Lung Cancer

Breathing new life into the study of COPD with genes identified from genome-wide association studies

COPD is a major cause of morbidity and mortality globally. While the significance of environmental exposures in disease pathogenesis is well established, the functional contribution of genetic factors has only in recent years drawn attention. Notably, many genes associated with COPD risk are also linked with lung function. Because reduced lung function precedes COPD onset, this association is consistent with the possibility that derangements leading to COPD could arise during lung development. In this review, we summarise the role of leading genes (HHIP, FAM13A, DSP, AGER and TGFB2) identified by genome-wide association studies in lung development and COPD. Because many COPD genome-wide association study genes are enriched in lung epithelial cells, we focus on the role of these genes in the lung epithelium in development, homeostasis and injury.

Chromosome-level genome assembly and population genomics of Mongolian racerunner (Eremias argus) provide insights into high-altitude adaptation in lizards

BACKGROUND

Although the extreme environmental adaptation of organisms is a hot topic in evolutionary biology, genetic adaptation to high-altitude environment remains poorly characterized in ectothermic animals. Squamates are among the most diverse terrestrial vertebrates, with tremendous ecological plasticity and karyotype diversity, and are a unique model system to investigate the genetic footprints of adaptation.

RESULTS

We report the first chromosome-level assembly of the Mongolian racerunner (Eremias argus) and our comparative genomics analyses found that multiple chromosome fissions/fusions events are unique to lizards. We further sequenced the genomes of 61 Mongolian racerunner individuals that were collected from altitudes ranging from ~ 80 to ~ 2600 m above sea level (m.a.s.l.). Population genomic analyses revealed many novel genomic regions under strong selective sweeps in populations endemic to high altitudes. Genes embedded in those genomic regions are mainly associated with energy metabolism and DNA damage repair pathways. Moreover, we identified and validated two substitutions of PHF14 that may enhance the lizards' tolerance to hypoxia at high altitudes.

CONCLUSIONS

Our study reveals the molecular mechanism of high-altitude adaptation in ectothermic animal using lizard as a research subject and provides a high-quality lizard genomic resource for future research.

FAM13A promotes proliferation of bovine preadipocytes by targeting Hypoxia-Inducible factor-1 signaling pathway

The family with sequence similarity 13 member A (FAM13A) gene has been discovered in recent years and is related to metabolism. In this study, the function of FAM13A in precursor adipocyte proliferation in Qinchuan cattle was investigated using fluorescence quantitative polymerase chain reaction (PCR), western blotting, 5-ethynyl-2'-deoxyuridine staining, and other tests. FAM13A promoted precursor adipocyte proliferation. To determine the pathway FAM13A was involved in, transcriptome sequencing, fluorescence quantitative PCR, western blotting, and other tests were used, which identified the hypoxia inducible factor-1 (HIF-1) signalling pathway. Finally, cobalt chloride, a chemical mimic of hypoxia, was used to treat precursor adipocytes. mRNA and protein levels of FAM13A were significantly increased after hypoxia. Thus, FAM13A promoted bovine precursor adipocyte proliferation by inhibiting the HIF-1 signalling pathway, whereas chemically induced hypoxia negatively regulated FAM13A expression, regulating cell proliferation.

FAM13A as potential therapeutic target in modulating TGF-β-induced airway tissue remodeling in COPD

Genome-wide association studies have shown that a gene variant in the Family with sequence similarity 13, member A (FAM13A) is strongly associated with reduced lung function and the appearance of respiratory symptoms in patients with chronic obstructive pulmonary disease (COPD). A key player in smoking-induced tissue injury and airway remodeling is the transforming growth factor-β1 (TGF-β1). To determine the role of FAM13A in TGF-β1 signaling, FAM13A^-/- airway epithelial cells were generated using CRISPR-Cas9, whereas overexpression of FAM13A was achieved using lipid nanoparticles. Wild-type (WT) and FAM13A^-/- cells were treated with TGF-β1, followed by gene and/or protein expression analyses. FAM13A^-/- cells augmented TGF-β1-induced increase in collagen type 1 (COL1A1), matrix metalloproteinase 2 (MMP2), expression compared with WT cells. This effect was mediated by an increase in β-catenin (CTNNB1) expression in FAM13A^-/- cells compared with WT cells after TGF-β1 treatment. FAM13A overexpression was partially protective from TGF-β1-induced COL1A1 expression. Finally, we showed that airway epithelial-specific FAM13A protein expression is significantly increased in patients with severe COPD compared with control nonsmokers, and negatively correlated with lung function. In contrast, β-catenin (CTNNB1), which has previously been linked to be regulated by FAM13A, is decreased in the airway epithelium of smokers with COPD compared with non-COPD subjects. Together, our data showed that FAM13A may be protective from TGF-β1-induced fibrotic response in the airway epithelium via sequestering CTNNB1 from its regulation on downstream targets. Therapeutic increase in FAM13A expression in the airway epithelium of smokers at risk for COPD, and those with mild COPD, may reduce the extent of airway tissue remodeling.

Hypoxia-Induced FAM13A Regulates the Proliferation and Metastasis of Non-Small Cell Lung Cancer Cells

Hypoxia in non-small cell lung cancer (NSCLC) affects cancer progression, metastasis and metabolism. We previously showed that FAM13A was induced by hypoxia in NSCLC but the biological function of this gene has not been fully elucidated. This study aimed to investigate the role of hypoxia-induced FAM13A in NSCLC progression and metastasis. Lentiviral shRNAs were used for FAM13A gene silencing in NSCLC cell lines (A549, CORL-105). MTS assay, cell tracking VPD540 dye, wound healing assay, invasion assay, BrdU assay and APC Annexin V staining assays were performed to examine cell proliferation ability, migration, invasion and apoptosis rate in NSCLC cells. The results of VPD540 dye and MTS assays showed a significant reduction in cell proliferation after FAM13A knockdown in A549 cells cultured under normal and hypoxia (1% O₂) conditions (p < 0.05), while the effect of FAM13A downregulation on CORL-105 cells was observed after 96 h exposition to hypoxia. Moreover, FAM13A inhibition induced S phase cell cycle arrest in A549 cells under hypoxia conditions. Silencing of FAM13A significantly suppressed migration of A549 and CORL-105 cells in both oxygen conditions, especially after 72 and 96 h (p < 0.001 in normoxia, p < 0.01 after hypoxia). It was showed that FAM13A reduction resulted in disruption of the F-actin cytoskeleton altering A549 cell migration. Cell invasion rates were significantly decreased in A549 FAM13A depleted cells compared to controls (p < 0.05), mostly under hypoxia. FAM13A silencing had no effect on apoptosis induction in NSCLC cells. In the present study, we found that FAM13A silencing has a negative effect on proliferation, migration and invasion activity in NSCLC cells in normal and hypoxic conditions. Our data demonstrated that FAM13A depleted post-hypoxic cells have a decreased cell proliferation ability and metastatic potential, which indicates FAM13A as a potential therapeutic target in lung cancer.

A Novel Promoter CpG-Based Signature for Long-Term Survival Prediction of Breast Cancer Patients

DNA methylation has been reported as one of the most critical epigenetic aberrations during the tumorigenesis and development of breast cancer (BC). This study explored a novel promoter CpG-based signature for long-term survival prediction of BC patients. We used The Cancer Genome Atlas (TCGA) data as training set, and results were validated in an independent dataset from Gene Expression Omnibus (GEO). First, the differential methylation CpG sites were screened in TCGA dataset, of which the candidate promoter CpG sites were preliminarily identified with the univariate Cox regression analysis and the least absolute shrinkage and selection operator regression analysis. Second, the signature was constructed with stepwise regression analysis and multivariate Cox proportional hazards model, which was validated with the survival analysis of two cohorts each from TCGA and GEO databases. The 10-year receiver operating characteristic curves of risk score presented an area under the curve of over 0.7 for both cohorts. A nomogram was also constructed and released. Moreover, Gene Set Enrichment Analysis was performed to identify the more active pathways in high-risk patients. The CpG sites-target gene correlations and differential methylation regions were further explored. In conclusion, the promoter CpG-based signature exhibited good prognostic prediction efficacy in the long-term overall survival of BC patients.

Hypoxia Activates SOX5/Wnt/β-Catenin Signaling by Suppressing MiR-338-3p in Gastric Cancer

MicroRNAs are known to be important in a variety of cancer types. The specific expression and roles of miR-338-3p in the context of gastric cancer, however, remains largely unknown. In this study, we found that miR-338-3p was expressed significantly lower in established/primary human gastric cancer cells than that in human gastric epithelial cells; miR-338-3p is also decreased in human gastric cancer tissues and was positively associated with the worse prognosis of patients with gastric cancer. Enforced expression of miR-338-3p could inhibit cell growth, survival, and proliferation, while inducing cell apoptosis. In addition, miR-338-3p negatively regulated SOX5 expression through directly binding to the 3′-untranslated region of SOX5, and an inverse correlation was found between miR-338-3p and SOX5 messenger RNA expression in gastric cancer tissues. Furthermore, miR-338-3p-induced inactivation of Wnt/β-catenin signaling was greatly abrogated by SOX5 upregulation. Finally, we found that hypoxic conditions were linked with reduced miR-338-3p expression in the context of gastric cancer. In conclusion, miR-338-3p acts as a tumor suppressor in gastric cancer, possibly by directly targeting SOX5 and blocking Wnt/β-catenin signaling. These findings might provide novel therapeutic targets for gastric cancer.

Hypoxia causes sex-specific hepatic toxicity at the transcriptome level in marine medaka (Oryzias melastigma)

Hypoxia, a low environmental oxygen level, is a common problem in the ocean globally. Hypoxia has been known to cause disruption to the endocrine system of marine organisms in both laboratory and field studies. Our previous studies have demonstrated the sex-specific response to hypoxia in the neural and reproductive systems of marine fish. In the current report, we aim to study the sex-specific hepatic response of fish at the transcriptome level to hypoxic stress. By using a comparative transcriptome analysis, followed by a systematic bioinformatics analysis including Database for Annotation, Visualization and Integrated Discovery (DAVID) and Ingenuity Pathway Analysis (IPA), we found that hypoxia altered expression of genes related to cell proliferation and apoptosis of hepatocytes, which are associated with human pathologies, such as liver inflammation hepatic steatosis and steatohepatitis. Furthermore, we observed sex-specific responses in the livers of fish through different cell signaling pathways. In female fish, hypoxia causes dysregulation of expression of genes related to impairment in endoplasmic reticulum structure and liver metabolism. In male fish, genes associated with redox homeostasis and fatty acid metabolism were altered by hypoxic stress. The findings of this study support the notion that hypoxia could cause sex-specific changes (hepatic toxicity and changes) in marine fish.

Functional Variant in 3'UTR of FAM13A Is Potentially Associated with Susceptibility and Survival of Lung Squamous Carcinoma

FAM13A is associated with aging lung disease (primarily chronic obstructive pulmonary disorder and pulmonary fibrosis) and shows stable expression throughout lung development. However, a few systematic studies of FAM13A have been conducted to assess the pathogenesis of lung cancer, particularly susceptibility. We predicted that single-nucleotide polymorphisms (SNPs) in FAM13A may be associated with lung cancer development. We systematically selected five functional SNPs (rs2602120, rs3017895, rs9224, rs7657817, and rs3756050) and genotyped them with the Genesky proprietary improved Multiligase Detection Reaction multiplex SNP genotyping system in a case-control study of 626 lung cancer cases and 667 cancer-free controls. The functional effects of FAM13A and specific miRNAs (miRNA-22-5p and miRNA-1301-3p) were evaluated based on The Cancer Genome Atlas database. We found that rs9224 in the 3' untranslated region (UTR) of FAM13A was potentially associated with an increased risk of lung squamous carcinoma (LUSQ) (additive model: odds ratio = 1.47, 95% confidence interval = 1.04-2.07, p = 0.028). In addition, the results of expression quantitative trait loci analysis suggested that the rs9224 polymorphism affects the expression of FAM13A (p = 0.050) and miRNA-22-5p (p = 0.031) in LUSQ. Further, survival analysis indicated decreased overall survival in the presence of the variant alleles of rs9224 (p = 0.048). The present results indicate that variant genotypes of rs9224 in the FAM13A 3'UTR may modify LUSQ susceptibility by affecting the binding of miRNA-22-5p and predict a poor prognosis of patients with LUSQ.

What lies beneath? Molecular evolution during the radiation of caecilian amphibians

Background

Evolution leaves an imprint in species through genetic change. At the molecular level, evolutionary changes can be explored by studying ratios of nucleotide substitutions. The interplay among molecular evolution, derived phenotypes, and ecological ranges can provide insights into adaptive radiations. Caecilians (order Gymnophiona), probably the least known of the major lineages of vertebrates, are limbless tropical amphibians, with adults of most species burrowing in soils (fossoriality). This enigmatic order of amphibians are very distinct phenotypically from other extant amphibians and likely from the ancestor of Lissamphibia, but little to nothing is known about the molecular changes underpinning their radiation. We hypothesised that colonization of various depths of tropical soils and of freshwater habitats presented new ecological opportunities to caecilians.

Results

A total of 8540 candidate groups of orthologous genes from transcriptomic data of five species of caecilian amphibians and the genome of the frog Xenopus tropicalis were analysed in order to investigate the genetic machinery behind caecilian diversification. We found a total of 168 protein-coding genes with signatures of positive selection at different evolutionary times during the radiation of caecilians. The majority of these genes were related to functional elements of the cell membrane and extracellular matrix with expression in several different tissues. The first colonization of the tropical soils was connected to the largest number of protein-coding genes under positive selection in our analysis. From the results of our study, we highlighted molecular changes in genes involved in perception, reduction-oxidation processes, and aging that likely were involved in the adaptation to different soil strata.

Conclusions

The genes inferred to have been under positive selection provide valuable insights into caecilian evolution, potentially underpin adaptations of caecilians to their extreme environments, and contribute to a better understanding of fossorial adaptations and molecular evolution in vertebrates.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5694-1) contains supplementary material, which is available to authorized users.

Trade-offs in aging lung diseases: a review on shared but opposite genetic risk variants in idiopathic pulmonary fibrosis, lung cancer and chronic obstructive pulmonary disease

PURPOSE OF REVIEW

The process of aging involves biological changes that increases susceptibility for disease. In the aging lung disease IPF, GWAS studies identified genes associated with risk for disease. Recently, several of these genes were also found to be involved in risk for COPD or lung cancer. This review describes GWAS-derived risk genes for IPF that overlap with risk genes for lung cancer or COPD.

RECENT FINDINGS

Risk genes that overlap between aging lung diseases, include FAM13A, DSP and TERT. Most interestingly, disease predisposing alleles for IPF are opposite to those for COPD or lung cancer. Studies show that the alleles are associated with differential gene expression and with physiological traits in the general population. The opposite allelic effect sizes suggest the presence of trade-offs in the aging lung. For TERT, the trade-off involves cellular senescence versus proliferation and repair. For FAM13A and DSP, trade-offs may involve protection from noxious gases or tissue integrity.

SUMMARY

The overlap in risk genes in aging lung diseases provides evidence that processes associated with FAM13A, DSP and TERT are important for healthy aging. The opposite effect size of the disease risk alleles may represent trade-offs, for which a model involving an apicobasal gene expression gradient is presented.

Integration of Molecular Interactome and Targeted Interaction Analysis to Identify a COPD Disease Network Module

The polygenic nature of complex diseases offers potential opportunities to utilize network-based approaches that leverage the comprehensive set of protein-protein interactions (the human interactome) to identify new genes of interest and relevant biological pathways. However, the incompleteness of the current human interactome prevents it from reaching its full potential to extract network-based knowledge from gene discovery efforts, such as genome-wide association studies, for complex diseases like chronic obstructive pulmonary disease (COPD). Here, we provide a framework that integrates the existing human interactome information with experimental protein-protein interaction data for FAM13A, one of the most highly associated genetic loci to COPD, to find a more comprehensive disease network module. We identified an initial disease network neighborhood by applying a random-walk method. Next, we developed a network-based closeness approach (C_AB) that revealed 9 out of 96 FAM13A interacting partners identified by affinity purification assays were significantly close to the initial network neighborhood. Moreover, compared to a similar method (local radiality), the C_AB approach predicts low-degree genes as potential candidates. The candidates identified by the network-based closeness approach were combined with the initial network neighborhood to build a comprehensive disease network module (163 genes) that was enriched with genes differentially expressed between controls and COPD subjects in alveolar macrophages, lung tissue, sputum, blood, and bronchial brushing datasets. Overall, we demonstrate an approach to find disease-related network components using new laboratory data to overcome incompleteness of the current interactome.