Glycosylphosphatidylinositol biosynthesis and remodeling are required for neural tube closure, heart development, and cranial neural crest cell survival

Adaptive introgression reveals the genetic basis of a sexually selected syndrome in wall lizards

The joint expression of particular colors, morphologies, and behaviors is a common feature of adaptation, but the genetic basis for such "phenotypic syndromes" remains poorly understood. Here, we identified a complex genetic architecture associated with a sexually selected syndrome in common wall lizards, by capitalizing on the adaptive introgression of coloration and morphology into a distantly related lineage. Consistent with the hypothesis that the evolution of phenotypic syndromes in vertebrates is facilitated by developmental linkage through neural crest cells, most of the genes associated with the syndrome are involved in neural crest cell regulation. A major locus was a ~400-kb region, characterized by standing structural genetic variation and previously implied in the evolutionary innovation of coloration and beak size in birds. We conclude that features of the developmental and genetic architecture contribute to maintaining trait integration, facilitating the extensive and rapid introgressive spread of suites of sexually selected characters.

Protein lipidation in health and disease: molecular basis, physiological function and pathological implication

Posttranslational modifications increase the complexity and functional diversity of proteins in response to complex external stimuli and internal changes. Among these, protein lipidations which refer to lipid attachment to proteins are prominent, which primarily encompassing five types including S-palmitoylation, N-myristoylation, S-prenylation, glycosylphosphatidylinositol (GPI) anchor and cholesterylation. Lipid attachment to proteins plays an essential role in the regulation of protein trafficking, localisation, stability, conformation, interactions and signal transduction by enhancing hydrophobicity. Accumulating evidence from genetic, structural, and biomedical studies has consistently shown that protein lipidation is pivotal in the regulation of broad physiological functions and is inextricably linked to a variety of diseases. Decades of dedicated research have driven the development of a wide range of drugs targeting protein lipidation, and several agents have been developed and tested in preclinical and clinical studies, some of which, such as asciminib and lonafarnib are FDA-approved for therapeutic use, indicating that targeting protein lipidations represents a promising therapeutic strategy. Here, we comprehensively review the known regulatory enzymes and catalytic mechanisms of various protein lipidation types, outline the impact of protein lipidations on physiology and disease, and highlight potential therapeutic targets and clinical research progress, aiming to provide a comprehensive reference for future protein lipidation research.

Drosophila models of phosphatidylinositol glycan biosynthesis class A congenital disorder of glycosylation (PIGA-CDG) mirror patient phenotypes

Mutations in the phosphatidylinositol glycan biosynthesis class A (PIGA) gene cause a rare, X-linked recessive congenital disorder of glycosylation. Phosphatidylinositol glycan biosynthesis class A congenital disorder of glycosylation (PIGA-CDG) is characterized by seizures, intellectual and developmental delay, and congenital malformations. The PIGA gene encodes an enzyme involved in the first step of glycosylphosphatidylinositol (GPI) anchor biosynthesis. There are over 100 GPI-anchored proteins that attach to the cell surface and are involved in cell signaling, immunity, and adhesion. Little is known about the pathophysiology of PIGA-CDG. Here, we describe the first Drosophila model of PIGA-CDG and demonstrate that loss of PIG-A function in Drosophila accurately models the human disease. As expected, complete loss of PIG-A function is larval lethal. Heterozygous null animals appear healthy but, when challenged, have a seizure phenotype similar to what is observed in patients. To identify the cell-type specific contributions to disease, we generated neuron- and glia-specific knockdown of PIG-A. Neuron-specific knockdown resulted in reduced lifespan and a number of neurological phenotypes but no seizure phenotype. Glia-knockdown also reduced lifespan and, notably, resulted in a very strong seizure phenotype. RNA sequencing analyses demonstrated that there are fundamentally different molecular processes that are disrupted when PIG-A function is eliminated in different cell types. In particular, loss of PIG-A in neurons resulted in upregulation of glycolysis, but loss of PIG-A in glia resulted in upregulation of protein translation machinery. Here, we demonstrate that Drosophila is a good model of PIGA-CDG and provide new data resources for future study of PIGA-CDG and other GPI anchor disorders.

Neuronal Wiring Receptors Dprs and DIPs Are GPI Anchored and This Modification Contributes to Their Cell Surface Organization

The Drosophila Dpr and DIP proteins belong to the immunoglobulin superfamily of cell surface proteins (CSPs). Their hetero- and homophilic interactions have been implicated in a variety of neuronal functions, including synaptic connectivity, cell survival, and axon fasciculation. However, the signaling pathways underlying these diverse functions are unknown. To gain insight into Dpr-DIP signaling, we sought to examine how these CSPs are associated with the membrane. Specifically, we asked whether Dprs and DIPs are integral membrane proteins or membrane anchored through the addition of glycosylphosphatidylinositol (GPI) linkage. We demonstrate that most Dprs and DIPs are GPI anchored to the membrane of insect cells and validate these findings for some family members in vivo using Drosophila larvae, where GPI anchor cleavage results in loss of surface labeling. Additionally, we show that GPI cleavage abrogates aggregation of insect cells expressing cognate Dpr-DIP partners. To test if the GPI anchor affects Dpr and DIP localization, we replaced it with a transmembrane domain and observed perturbation of subcellular localization on motor neurons and muscles. These data suggest that membrane anchoring of Dprs and DIPs through GPI linkage is required for localization and that Dpr-DIP intracellular signaling likely requires transmembrane coreceptors.

SMPD4 mediated sphingolipid metabolism regulates brain and primary cilia development

Genetic variants in multiple sphingolipid biosynthesis genes cause human brain disorders. A recent study collected patients from twelve unrelated families with variants in the gene SMPD4 , a neutral sphingomyelinase which metabolizes sphingomyelin into ceramide at an early stage of the biosynthesis pathway. These patients have severe developmental brain malformations including microcephaly and cerebellar hypoplasia. However, the mechanism of SMPD4 was not known and we pursued a new mouse model. We hypothesized that the role of SMPD4 in producing ceramide is important for making primary cilia, a crucial organelle mediating cellular signaling. We found that the mouse model has cerebellar hypoplasia due to failure of Purkinje cell development. Human induced pluripotent stem cells exhibit neural progenitor cell death and have shortened primary cilia which is rescued by adding exogenous ceramide. SMPD4 production of ceramide is crucial for human brain development.

GPI-anchoring disorders and the heart: Is cardiomyopathy an overlooked feature?

Glycosylphosphatidylinositol anchoring disorders (GPI-ADs) are a subgroup of congenital disorders of glycosylation. GPI biosynthesis requires proteins encoded by over 30 genes of which 24 genes are linked to neurodevelopmental disorders. Patients, especially those with PIGA-encephalopathy, have a high risk of premature mortality which sometimes is attributed to cardiomyopathy. We aimed to explore the occurrence of cardiomyopathy among patients with GPI-ADs and to raise awareness about this potentially lethal feature. Unpublished patients with genetically proven GPI-ADs and cardiomyopathy were identified through an international collaboration and recruited through the respective clinicians. We also reviewed the literature for published patients with cardiomyopathy and GPI-AD and contacted the corresponding authors for additional information. We identified four novel and unrelated patients with GPI-AD and cardiomyopathy. Cardiomyopathy was diagnosed before adulthood and was the cause of early demise in two patients. Only one patients underwent cardiac workup after being diagnosed with a GPI-AD. All were diagnosed with PIGA-encephalopathy and three had a disease-causing variant at the same residue. The literature reports five additional children with GPI-AD related cardiomyopathy, three of which died before adulthood. We have shown that patients with GPI-ADs are at risk of developing cardiomyopathy and that regular cardiac workup with echocardiography is necessary.

Drosophila models of PIGA-CDG mirror patient phenotypes

Mutations in the phosphatidylinositol glycan biosynthesis class A (PIGA) gene cause a rare, X-linked recessive congenital disorder of glycosylation (CDG). PIGA-CDG is characterized by seizures, intellectual and developmental delay, and congenital malformations. The PIGA gene encodes an enzyme involved in the first step of GPI anchor biosynthesis. There are over 100 GPI anchored proteins that attach to the cell surface and are involved in cell signaling, immunity, and adhesion. Little is known about the pathophysiology of PIGA-CDG. Here we describe the first Drosophila model of PIGA-CDG and demonstrate that loss of PIG-A function in Drosophila accurately models the human disease. As expected, complete loss of PIG-A function is larval lethal. Heterozygous null animals appear healthy, but when challenged, have a seizure phenotype similar to what is observed in patients. To identify the cell-type specific contributions to disease, we generated neuron- and glia-specific knockdown of PIG-A. Neuron-specific knockdown resulted in reduced lifespan and a number of neurological phenotypes, but no seizure phenotype. Glia-knockdown also reduced lifespan and, notably, resulted in a very strong seizure phenotype. RNAseq analyses demonstrated that there are fundamentally different molecular processes that are disrupted when PIG-A function is eliminated in different cell types. In particular, loss of PIG-A in neurons resulted in upregulation of glycolysis, but loss of PIG-A in glia resulted in upregulation of protein translation machinery. Here we demonstrate that Drosophila is a good model of PIGA-CDG and provide new data resources for future study of PIGA-CDG and other GPI anchor disorders.

The critical role of cardiolipin in metazoan differentiation, development, and maturation

Cardiolipins are phospholipids that are central to proper mitochondrial functioning. Because mitochondria play crucial roles in differentiation, development, and maturation, we would also expect cardiolipin to play major roles in these processes. Indeed, cardiolipin has been implicated in the mechanism of three human diseases that affect young infants, implying developmental abnormalities. In this review, we will: (1) Review the biology of cardiolipin; (2) Outline the evidence for essential roles of cardiolipin during organismal development, including embryogenesis and cell maturation in vertebrate organisms; (3) Place the role(s) of cardiolipin during embryogenesis within the larger context of the roles of mitochondria in development; and (4) Suggest avenues for future research.

SUMOylation of Nuclear γ-Actin by SUMO2 supports DNA Damage Repair against Myocardial Ischemia-Reperfusion Injury

Myocardial infarction triggers oxidative DNA damage, apoptosis and adverse cardiac remodeling in the heart. Small ubiquitin-like modifier (SUMO) proteins mediate post-translational SUMOylation of the cardiac proteins in response to oxidative stress signals. Upregulation of isoform SUMO2 could attenuate myocardial injury via increasing protein SUMOylation. The present study aimed to discover the identity and cardioprotective activities of SUMOylated proteins. A plasmid vector for expressing N-Strep-SUMO2 protein was generated and introduced into H9c2 rat cardiomyocytes. The SUMOylated proteins were isolated with Strep-Tactin^® agarose beads and identified by MALDI-TOF-MS technology. As a result, γ-actin was identified from a predominant protein band of ~42 kDa and verified by Western blotting. The roles of SUMO2 and γ-actin SUMOylation were subsequently determined in a mouse model of myocardial infarction induced by ligating left anterior descending coronary artery and H9c2 cells challenged by hypoxia-reoxygenation. In vitro lentiviral-mediated SUMO2 expression in H9c2 cells were used to explore the role of SUMOylation of γ-actin. SUMOylation of γ-actin by SUMO2 was proven to be a new cardioprotective mechanism from the following aspects: 1) SUMO2 overexpression reduced the number of TUNEL positive cells, the levels of 8-OHdG and p-γ-H2ax while promoted the nuclear deposition of γ-actin in mouse model and H9c2 cell model of myocardial infarction; 2) SUMO-2 silencing decreased the levels of nuclear γ-actin and SUMOylation while exacerbated DNA damage; 3) Mutated γ-actin (K⁶⁸R/K²⁸⁴R) void of SUMOylation sites failed to protect cardiomyocytes against hypoxia-reoxygenation challenge. The present study suggested that SUMO2 upregulation promoted DNA damage repair and attenuated myocardial injury via increasing SUMOylation of γ-actin in the cell nucleus.

Genomic, transcriptomic, and metabolomic profiles of hiPSC-derived dopamine neurons from clinically discordant brothers with identical PRKN deletions

We previously reported on two brothers who carry identical compound heterozygous PRKN mutations yet present with significantly different Parkinson's Disease (PD) clinical phenotypes. Juvenile cases demonstrate that PD is not necessarily an aging-associated disease. Indeed, evidence for a developmental component to PD pathogenesis is accumulating. Thus, we hypothesized that the presence of additional genetic modifiers, including genetic loci relevant to mesencephalic dopamine neuron development, could potentially contribute to the different clinical manifestations of the two brothers. We differentiated human-induced pluripotent stem cells (hiPSCs) derived from the two brothers into mesencephalic neural precursor cells and early postmitotic dopaminergic neurons and performed wholeexome sequencing and transcriptomic and metabolomic analyses. No significant differences in the expression of canonical dopamine neuron differentiation markers were observed. Yet our transcriptomic analysis revealed a significant downregulation of the expression of three neurodevelopmentally relevant cell adhesion molecules, CNTN6, CNTN4 and CHL1, in the cultures of the more severely affected brother. In addition, several HLA genes, known to play a role in neurodevelopment, were differentially regulated. The expression of EN2, a transcription factor crucial for mesencephalic dopamine neuron development, was also differentially regulated. We further identified differences in cellular processes relevant to dopamine metabolism. Lastly, wholeexome sequencing, transcriptomics and metabolomics data all revealed differences in glutathione (GSH) homeostasis, the dysregulation of which has been previously associated with PD. In summary, we identified genetic differences which could potentially, at least partially, contribute to the discordant clinical PD presentation of the two brothers.

Loss of the N-acetylgalactosamine side chain of the GPI-anchor impairs bone formation and brain functions and accelerates the prion disease pathology

Glycosylphosphatidylinositol (GPI) is a posttranslational glycolipid modification of proteins that anchors proteins in lipid rafts on the cell surface. Although some GPI-anchored proteins (GPI-APs), including the prion protein PrP^C, have a glycan side chain composed of N-acetylgalactosamine (GalNAc)−galactose−sialic acid on the core structure of GPI glycolipid, in vivo functions of this GPI-GalNAc side chain are largely unresolved. Here, we investigated the physiological and pathological roles of the GPI-GalNAc side chain in vivo by knocking out its initiation enzyme, PGAP4, in mice. We show that Pgap4 mRNA is highly expressed in the brain, particularly in neurons, and mass spectrometry analysis confirmed the loss of the GalNAc side chain in PrP^C GPI in PGAP4-KO mouse brains. Furthermore, PGAP4-KO mice exhibited various phenotypes, including an elevated blood alkaline phosphatase level, impaired bone formation, decreased locomotor activity, and impaired memory, despite normal expression levels and lipid raft association of various GPI-APs. Thus, we conclude that the GPI-GalNAc side chain is required for in vivo functions of GPI-APs in mammals, especially in bone and the brain. Moreover, PGAP4-KO mice were more vulnerable to prion diseases and died earlier after intracerebral inoculation of the pathogenic prion strains than wildtype mice, highlighting the protective roles of the GalNAc side chain against prion diseases.

CNS glycosylphosphatidylinositol deficiency results in delayed white matter development, ataxia and premature death in a novel mouse model

The glycosylphosphatidylinositol (GPI) anchor is a post-translational modification added to approximately 150 different proteins to facilitate proper membrane anchoring and trafficking to lipid rafts. Biosynthesis and remodeling of the GPI anchor requires the activity of over 20 distinct genes. Defects in the biosynthesis of GPI anchors in humans lead to inherited glycosylphosphatidylinositol deficiency (IGD). IGD patients display a wide range of phenotypes though the central nervous system (CNS) appears to be the most commonly affected tissue. A full understanding of the etiology of these phenotypes has been hampered by the lack of animal models due to embryonic lethality of GPI biosynthesis gene null mutants. Here we model IGD by genetically ablating GPI production in the CNS with a conditional mouse allele of phosphatidylinositol glycan anchor biosynthesis, class A (Piga) and Nestin-Cre. We find that the mutants do not have structural brain defects but do not survive past weaning. The mutants show progressive decline with severe ataxia consistent with defects in cerebellar development. We show that the mutants have reduced myelination and defective Purkinje cell development. Surprisingly, we found that Piga was expressed in a fairly restricted pattern in the early postnatal brain consistent with the defects we observed in our model. Thus, we have generated a novel mouse model of the neurological defects of IGD which demonstrates a critical role for GPI biosynthesis in cerebellar and white matter development.

A CRISPR-Cas9-engineered mouse model for GPI-anchor deficiency mirrors human phenotypes and exhibits hippocampal synaptic dysfunctions

Inherited GPI-anchor biosynthesis deficiencies (IGDs) explain many cases of syndromic intellectual disability. Although diagnostic methods are improving, the pathophysiology underlying the disease remains unclear. Furthermore, we lack rodent models suitable for characterizing cognitive and social disabilities. To address this issue, we generated a viable mouse model for an IGD that mirrors the condition in human patients with a behavioral phenotype and susceptibility to epilepsy. Using this model, we obtained neurological insights such as deficits in synaptic transmission that will facilitate understanding of the pathophysiology of IGDs.

Pathogenic germline mutations in PIGV lead to glycosylphosphatidylinositol biosynthesis deficiency (GPIBD). Individuals with pathogenic biallelic mutations in genes of the glycosylphosphatidylinositol (GPI)-anchor pathway exhibit cognitive impairments, motor delay, and often epilepsy. Thus far, the pathophysiology underlying the disease remains unclear, and suitable rodent models that mirror all symptoms observed in human patients have not been available. Therefore, we used CRISPR-Cas9 to introduce the most prevalent hypomorphic missense mutation in European patients, Pigv:c.1022C > A (p.A341E), at a site that is conserved in mice. Mirroring the human pathology, mutant Pigv^341E mice exhibited deficits in motor coordination, cognitive impairments, and alterations in sociability and sleep patterns, as well as increased seizure susceptibility. Furthermore, immunohistochemistry revealed reduced synaptophysin immunoreactivity in Pigv^341E mice, and electrophysiology recordings showed decreased hippocampal synaptic transmission that could underlie impaired memory formation. In single-cell RNA sequencing, Pigv^341E-hippocampal cells exhibited changes in gene expression, most prominently in a subtype of microglia and subicular neurons. A significant reduction in Abl1 transcript levels in several cell clusters suggested a link to the signaling pathway of GPI-anchored ephrins. We also observed elevated levels of Hdc transcripts, which might affect histamine metabolism with consequences for circadian rhythm. This mouse model will not only open the doors to further investigation into the pathophysiology of GPIBD, but will also deepen our understanding of the role of GPI-anchor–related pathways in brain development.

Limb-clasping, cognitive deficit and increased vulnerability to kainic acid-induced seizures in neuronal glycosylphosphatidylinositol deficiency mouse models

Posttranslational modification of a protein with glycosylphosphatidylinositol (GPI) is a conserved mechanism exists in all eukaryotes. Thus far, >150 human GPI-anchored proteins have been discovered and ~30 enzymes have been reported to be involved in the biosynthesis and maturation of mammalian GPI. Phosphatidylinositol glycan biosynthesis class A protein (PIGA) catalyzes the very first step of GPI anchor biosynthesis. Patients carrying a mutation of the PIGA gene usually suffer from inherited glycosylphosphatidylinositol deficiency (IGD) with intractable epilepsy and intellectual developmental disorder. We generated three mouse models with PIGA deficits specifically in telencephalon excitatory neurons (Ex-M-cko), inhibitory neurons (In-M-cko) or thalamic neurons (Th-H-cko), respectively. Both Ex-M-cko and In-M-cko mice showed impaired long-term fear memory and were more susceptible to kainic acid-induced seizures. In addition, In-M-cko demonstrated a severe limb-clasping phenotype. Hippocampal synapse changes were observed in Ex-M-cko mice. Our Piga conditional knockout mouse models provide powerful tools to understand the cell-type specific mechanisms underlying inherited GPI deficiency and to test different therapeutic modalities.

Using an aquatic model, Xenopus laevis, to uncover the role of chromodomain 1 in craniofacial disorders

The chromodomain family member chromodomain 1 (CHD1) has been shown to have numerous critical molecular functions including transcriptional regulation, splicing, and DNA repair. Complete loss of function of this gene is not compatible with life. On the other hand, missense and copy number variants of CHD1 can result in intellectual disabilities and craniofacial malformations in human patients including cleft palate and Pilarowski-Bjornsson Syndrome. We have used the aquatic developmental model organism Xenopus laevis, to determine a specific role for Chd1 in such cranioafcial disorders. Protein and gene knockdown techniques in Xenopus, including antisense oligos and mosaic Crispr/Cas9-mediated mutagenesis, recapitulated the craniofacial defects observed in humans. Further analysis indicated that embryos deficient in Chd1 had defects in cranial neural crest development and jaw cartilage morphology. Additionally, flow cytometry and immunohistochemistry revealed that decreased Chd1 resulted in increased in apoptosis in the developing head. Together, these experiments demonstrate that Chd1 is critical for fundamental processes and cell survival in craniofacial development. We also presented evidence that Chd1 is regulated by retinoic acid signaling during craniofacial development. Expression levels of chd1 mRNA, specifically in the head, were increased by RAR agonist exposure and decreased upon antagonist treatment. Subphenotypic levels of an RAR antagonist and Chd1 morpholinos synergized to result in orofacial defects. Further, RAR DNA binding sequences (RAREs) were detected in chd1 regulatory regions by bioinformatic analysis. In summary, by combining human genetics and experiments in an aquatic model we now have a better understanding of the role of CHD1 in craniofacial disorders.

Nano-sulforaphane attenuates PhIP-induced early abnormal embryonic neuro-development

BACKGROUND

2-Amino-1-methyl-6-phenylimidazo[4,5-b]pyrimidine (PhIP), one of the most abundant heterocyclic aromatic amines (HAA) formed by cooking meat at high temperatures, may modify humans and rodents through the metabolic process prior to affecting nervous system development. In humans and rodents may be modified by metabolic processes and then affecting nervous system development.

METHODS

In this paper, PhIP was used to prepare a chicken embryo model with abnormal embryonic nervous system defects. Sulforaphane (SFN) is a derivative of a glucosinolate, which is abundant in cruciferous vegetables, and can pass through the placental barrier. Moreover, SFN has antioxidant and anti-apoptotic functions and is considered as a bioactive antioxidant with significant neuroprotective effects. Nano-sulforaphane (Nano-SFN, sulforaphane nanoparticles) was prepared by self-assembly using biocompatible, biodegradable methoxy polyethylene glycol 5000-b-polyglutamic acid 10,000 (mPEG5K-PGA10K) as the substrate, to explore the new application of Nano-SFN and its modified compounds as leading compounds in protecting against the abnormal development of the embryonic nervous system.

RESULTS

The results show that Nano-SFN could protect against PhIP-induced central nervous system (CNS, derived from neural tube) and peripheral nervous system (PNS, derived from neural crest cells, NCCs) defects and neural tube defects (NTDs), and increase the embryo survival rate.

CONCLUSIONS

This study indicates that Nano-SFN can effectively alleviate the developmental defects of embryonic nervous system induced by PhIP in the microenvironment and has a protective effect on embryonic development. It not only helps with expanding the application of SFN and improving its medicinal value, but also provides a possibility of SFN being developed as a novel drug for neuroprotection.

Association Study of Genetic Variants in Autophagy Pathway and Risk of Non-syndromic Cleft Lip With or Without Cleft Palate

Although genetic variants in autophagy pathway genes were associated with the risk of oral cancers and early development in embryos, their associations with non-syndromic cleft lip with or without cleft palate (NSCL/P) risk remained unclear. A two-stage case-control study (2,027 NSCL/P cases and 1,843 controls) was performed to investigate the associations between single nucleotide polymorphisms (SNPs) in 23 autophagy pathway genes and NSCL/P susceptibility. The logistic regression model was used to calculate effects of SNPs on NSCL/P susceptibility. Gene-based analysis was performed via the sequence kernel association test (SKAT) and multi-marker analysis of genomic annotation (MAGMA) methods. Expression quantitative trait loci (eQTL) analysis was conducted using NSCL/P lip tissue samples. Gene expression during embryonic development was evaluated using RNA-Seq. Functional roles were explored by luciferase activity assay, cell apoptosis, proliferation, and cycle in vitro. Rs2301104 in HIF1A was significantly associated with NSCL/P susceptibility in the combined analysis (OR: 1.29, 95% CI: 1.09-1.29, P = 3.39 × 10^-03), and showed strong evidence of association heterogeneity (P = 9.06 × 10^-03) with obvious association in the female (OR: 1.80; 95% CI: 1.32-2.45; P = 1.79 × 10^-04). The G allele of rs2301104 was associated with enhanced transcription activity and high expression of HIF1A compared with that of C allele. Moreover, rs2301104 exhibited an eQTL effect for HIF1A with its GC/CC genotypes associated with decreased HIF1A expression compared with those with GG genotypes (P = 3.1 × 10^-2). Knockdown of HIF1A induced cell apoptosis and inhibited cell proliferation in human embryonic palate mesenchyme (HEPM) and human oral epithelium cells (HOEC). This study demonstrated that rs2301104 in autophagy pathway gene HIF1A was associated with susceptibility of NSCL/P.