Abstract
Defective left–right (LR) organization involving abnormalities in cilia ultrastructure causes laterality disorders including situs inversus (SI) and heterotaxy (Htx) with the prevalence approximately 1/10,000 births. In this study, we describe two unrelated family trios with abnormal cardiac LR patterning. Through whole-exome sequencing (WES), we identified compound heterozygous mutations (c.805-1G >C; p. Ile269GlnfsTer8/c.1117dupA; p.Thr373AsnfsTer19) (c.29T>C; p.Ile10Thr/c.356A>G; p.His119Arg) of NEK3, encoding a NIMA (never in mitosis A)-related kinase, in two affected individuals, respectively. Protein levels of NEK3 were abrogated in Patient-1 with biallelic loss-of function (LoF) NEK3 mutations that causes premature stop codon. Subsequence transcriptome analysis revealed that NNMT (nicotinamide N-methyltransferase) and SIRT2 (sirtuin2) was upregulated by NEK3 knockdown in human retinal pigment epithelial (RPE) cells in vitro, which associates α-tubulin deacetylation by western blot and immunofluorescence. Transmission electron microscopy (TEM) analysis further identified defective ciliary ultrastructure in Patient-1. Furthermore, inner ring components of nuclear pore complex (NPC) including nucleoporin (NUP)205, NUP188, and NUP155 were significantly downregulated in NEK3-silenced cells. In conclusion, we identified biallelic mutations of NEK3 predispose individual to abnormal cardiac left–right patterning via SIRT2-mediated α-tubulin deacetylation and downregulation of inner ring nucleoporins. Our study suggested that NEK3 could be a candidate gene for human ciliopathies.
Similar content being viewed by others
Introduction
Left–right (LR) asymmetry is an essential aspect of embryonic development in most bilateral vertebrates and invertebrates1. Accurate positioning of asymmetric organs along the LR axis is crucial for their proper function, and defects in LR organization cause a range of malformations including situs inversus (SI) totalis (SIT) and heterotaxy (Htx)2. SIT is a congenital condition in which the organs in the chest and abdomen are arranged in a complete mirror-image reversal of the usual positions; the prevalence of SIT is estimated to range from 1/25,000 to 1/8000 (ref. 3). Up to 20% of patients with SIT have Kartagener syndrome (KS) which is a trilogy of symptoms (nasal polyps, bronchiectasis, and SIT). KS is a subgroup of primary ciliary dyskinesia (PCD)4. Approximately 50% of patients with PCD present with SIT, designated as KS5. KS and PCD are associated with the same respiratory symptoms that are secondary to ultrastructural anomalies of cilia6. Different from SIT, Htx is a condition that involves the internal organs being abnormally arranged but not in a mirror image within the chest and abdomen7. Over 80% of indiciduals with Htx present with complex congenital heart disease (CHD) including transposition of the great arteries and double outlet right ventricle8. SIT patients have a lower risk of CHD (3–9%) than patients with Htx. However, the risk of CDH in SIT is still much higher than in the normal situation, situs solitus (0.6–0.8%)9.
Previous studies have established the fundamental roles of embryonic cilia-driven fluid flow in LR organization in mammals10,11. Classified by microtubular structure, two types of cilia are found: primary/non-motile cilia and regular motile cilia. Most regular motile cilia consist of a ring of nine peripheral microtubule doublets that surrounds a central pair of single microtubules (9 + 2 structure). In contrast to the 9 + 2 pattern of regular cilia, primary or sensory cilia are solitary, immotile organelles that have 9 + 0 microtubule configuration and are present on most cell types12. Primary cilia are microtubule-based organelles involved in a variety of signaling cascades13. Of the many signaling pathways associated with primary cilia, NIMA (never in mitosis A)-related kinases (Nek1-Nek11), a family of Ser-Thr kinases, have been extensively illustrated in many cellular processes14. Although the functions of mammalian Nek proteins remain largely unknown, common themes in Nek kinases are cilia, centrosomes, and microtubule organization15. Nek2 is crucial for proper mitosis through its effects on centrosomes16, and Nek1 and Nek8 are involved in primary cilium formation17,18. Nek6 can phosphorylate the kinesin Eg5, which regulates microtubule movement19. Nek8 phosphorylates Bicd2, a protein that is important for microtubule maintenance and axonal stability20. A recent study identified NEK10 as a ciliated cell-specific kinase whose activity regulates the motile ciliary proteome to promote ciliary length and mucociliary transport21. More importantly, several Nek kinases have been implicated in embryonic LR determination. Rare genic copy number variations in NEK2 were identified in patients with Htx22. Nek2 is a key switch balancing cilia biogenesis that be required for normal LR patterning in Xenopus23. Nek8 homozygous null mice exhibit randomization of LR asymmetry24.
Different from most NIMA family members that closely associated with cell cycle regulation, inhibition of Nek3 activity did not affect cell cycle progression25. Two haplotypes have been identified in NEK3 that differ by a 1 nt indel (rs3837575, 0.541 in GnomAD) at the end of exon 10 (supplementary Fig. S1). Both haplotypes appear to encode functional proteins via alternative in-frame splicing of exon 11. In the full-length protein (NM_002498), the length of exon 10 is 73 bp. However, in the haplotype that lacking exon 11 (NM_001146099), the length of exon 10 is 72 bp. Expression of a phospho-defective mutant or truncated Nek3 leads to distorted neuronal morphology and deacetylation of microtubules via HDAC6 (ref. 26). In addition, NEK3 was reported to play important roles in cancer progression. NEK3 activates Rac1 and contributes to prolactin-mediated breast cancer motility through paxillin phosphorylation27. Phosphorylation of NEK3 at T165 modulates focal adhesion remodeling necessary for breast cancer cell migration28. In the following, by performing trio-based whole-exome sequencing (WES), we identified biallelic mutations in NEK3 in two unrelated patients with abnormal cardiac LR patterning and found NEK3 deficiency promotes NNMT expression and downregulates expression of inner ring components of NPC, which potentially associates with defective ultrastructure of cilia and abnormal cardiac LR organization.
Materials and Methods
Subjects
In this study, we collected two unrelated family trios with abnormal cardiac LR patterning. The affected individual in Family-1 is a 6-year-old boy who had been diagnosed as isolated dextrocardia at the age of 2 months by X-ray and color ultrasonic diagnosis at Pediatric Cardiovascular Center of the Children’s Hospital affiliated to Fudan University, Shanghai, China. In the Family-2, a 29-year-old female was diagnosed as SI totalis with mild scoliosis by X-ray at Pukou District Central Hospital, Nanjing, China. The fresh peripheral blood samples of patients and unaffected family members were collected for genomic DNA and RNA isolation, WES, and real-time PCR. For studies of affected individual and their families, written informed consent was obtained from all participants prior to the start of the study. All procedures in the study were approved by the Medical Ethics Committee of Children’s Hospital of Fudan University (2016-079) (Shanghai, China).
Genetic analysis
WES was performed by utilizing the SureSelect human all exon platform (v.6; Agilent Technologies, Santa Clara, CA, USA). The Genome Analysis Toolkit software package (GATK4.1.2.0) was used for detection of single-nucleotide variants and indels. Variant filtering was performed according to minor allele frequencies (MAFs) under 0.01 in the GnomAD database (http://gnomad.broadinstitute.org/), given the rarity of situs abnormality and associations with protein-altering variants (missense, nonsense, and splice-site variants as well as coding indels). Recessive and homozygous variants were then maintained for subsequent analysis. The effects of the identified variants were assessed using SIFT (http://sift.jcvi.org), PolyPhen-2 (genetics.bwh.harvard.edu/pph2), and Mutationtaster (http://mutationtaster.org). Conservation analysis was performed using Multiple Sequence Alignment ClustalW2 (ebi.ac.uk/Tools/msa/clustalw2). Candidate causal variants of NEK3 were then confirmed by sanger sequencing and the sequences of primers are provided in Supplementary Table S1.
cDNA analysis
Total RNA was isolated from fresh peripheral blood samples of the Patient-1 and his unaffected parents by using the RNAprep pure Blood Kit (DP433, TIANGEN) followed by first-strand cDNA synthesis (KR106, TIANGEN). cDNA was then amplified by polymerase chain reaction (PCR) with specific NEK3 primers followed by gel extraction (DP209, TIANGEN) and Sanger sequencing. The primers used in PCR reactions are listed in Supplementary Table S1.
Cell culture
Human TERT retinal pigment epithelial (RPE) cells (ATCC®CRL-4000) were grown in Dulbecco’s modified Eagle’s medium (DMEM)/ F12 1:1: Liquid (HyClone) supplemented with 10% FBS (fetal bovine serum) (Gibco), 100 IU/ml penicillin, 100 μg/ml streptomycin at 37 °C in a 5% CO2 atmosphere.
siRNA silencing
SiRNA reagent kits and the negative control were ordered from RiboBio Co. Ltd (Guangzhou, China). Three pairs of siRNAs (#1-3) against human NEK3 (NM_002498) and two pairs of siRNAs (#1–2) against human NUP205 (NM_015135) were designed for gene silencing and siRNA sequences are provided in Supplementary Table S2. Briefly, human RPE cells were grown in the presence of 10% FBS under identical culture conditions and were then transfected with Ctrl or each NEK3 or NUP205 siRNA (75 nM) using the LipofectAMINE 3000 reagent (Invitrogen, CA, USA) according to the standard protocols. At 60-h post-transfection, cells were collected to evaluate knockdown efficiency by western blot.
RT-qPCR
Total RNA was isolated from Ctrl or NEK3-siRNA-treated human RPE cells using the RNAprep pure Cell Kit (DP430, Tiangen, Beijing, China) followed by first-strand cDNA synthesis. RT-qPCR was performed using the FastFire qPCR PreMix (SYBR Green) (FP207, TIANGEN) on an ABI StepOnePlus instrument. The delta-delta-Ct (ddCt) algorithm was utilized to analyze the relative changes in gene expression by using GAPDH as the internal reference (housekeeping gene). The data are presented as the mean ± SD of four independent experiments; each sample was assayed in triplicate in each experiment. RT-qPCR primer sequences are available in Supplementary Table S1 and were designed using Primer3Plus software to span exon–exon junctions.
Western blot
siRNA-treated human RPE cells or nasal mucosa tissues of the patient were lysed on ice for 30 min in RIPA lysis buffer (Sangon Biotech) and lysates were clarified by centrifugation at 12,000 r.p.m. for 10 min at 4 °C. Western blotting was performed following standard procedures. Membranes were blocked with 5% non-fat milk for 1 h at room temperature and then incubated with different primary antibodies in 0.5% bovine serum albumin (BSA) or 5% milk. The bands on immunoblots were quantified using Image J software.
Immunofluorescence
Human RPE cells were seeded on glass cover slips and grown in presence of 10% FBS under identical culture conditions. At 60-h post-transfection, cells were fixed and permeabilized for 10 min using 4% PFA and 0.8% Triton-X100. After blocking, cells were incubated with the primary antibody as aforementioned (anti-AcTub; Sigma-Aldrich, 1:500 dilution) overnight at 4 °C. Cells were then incubated with the secondary antibody (Alexa 488 anti-mouse; Thermo Fisher Scientific, 1:1000 dilution) for 2 h at room temperature, followed by staining with DAPI (Thermo Fisher Scientific)/PBS for 6 min. Confocal imaging was performed using an SP8 system (Leica), and images were processed with using the Leica AF software suite.
High-speed video microscopy and TEM analysis
Nasal tissues of Patient-1 and the healthy individual were collected and maintained in L-15 medium (Invitrogen, CA) for video microscopy using a Leica inverted microscope (Leica DMI300B) as described previously29,30. Movies were recorded at 200 frames/s at room temperature by use of a 680 PROSILICA GE camera (Allied Vision, PA). Transmission electron microscopy (TEM) analysis was performed at the Department of Electron Microscopy, Fudan University as previously29,30. Briefly, the samples were fixed with glutaraldehyde (2.5% w/v) in 0.1 M sodium cacodylate buffer. Following fixation in osmium tetroxide (1% w/v), samples were dehydrated through graded ethanol series and embedded in epoxy resin. Sections were cut at 50–70 nm thickness and stained with 2% methanolic uranyl acetate and Reynold’s lead citrate for analysis by using a TEM (JEM-1400; Jeol, Tokyo, Japan).
RNA sequencing
Transcriptome libraries were prepared using a TruSeq RNA Sample Preparation kit (Illumina) and then sequenced using the Illumina HiSeq PE platform at 2 × 151-bp read lengths. The expression level of each transcript was calculated according to fragments per kilobase of exon per million mapped reads (FPKM). Differentially expressed genes (DEGs) between the different groups were selected based on a log fold-change greater than 2 and a false-discovery rate of less than 0.05. Gene ontology (GO) functional-enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed using GOATOLLS and KOBAS. DEGs were considered significantly enriched when the p value was less than 0.05 following Bonferroni correction.
Antibodies
For immunoblots, proteins were probed with one of the following primary antibodies in either 5% non-fat milk or 5% bovine serum albumin: anti-HEC1 (#3613, 1:1000 dilution) and anti-NEK3 (#37636, 1:800 dilution) antibodies from Abcam; anti-XPO1 (#46249, 1:1000 dilution), anti-phospho-Akt (Ser473) (#4051, 1:1000 dilution), and anti-GAPDH (#5174, 1:3000 dilution) from Cell Signaling Technology; anti-NUP205 (Santa Cruz Biotechnology, 1:200 dilution), and anti-tubulin (acetylated) (T6793, Sigma-Aldrich, 1:1000 dilution for WB).
Statistics
The data are presented as the mean ± standard. Comparisons between groups were analyzed with two-tailed nonparametric Mann–Whitney test by using statistical software SPSS version 16 (SPSS Inc., Chicago, USA). A p value of less than 0.05 was considered statistically significant.
Results
Case presentation
In this study, we collected two unrelated family trios with abnormal cardiac LR patterning. The affected individual in Family-1 is a 6-year-old boy of Chinese origin; both his parents were unaffected (Fig. 1a). Chest X-ray showed isolated dextrocardia and raised left side of the diaphragm in this boy (Fig. 1b). Echocardiography did not reveal other detectable abnormalities in heart structure or motion. Typical PCD phenotypes including bronchitis, sinusitis, and otitis media were temporarily excluded by physical examination as previously29,30. In the Family-2, the patient is a 29-year-old female; her husband and 3-year-old son were healthy (Fig. 1a). Chest X-ray showed the mirror-image arrangement of the abdominal organs including heart and liver in the patient (Fig. 1b). She exhibited normal growth parameters without mental retardation and could withstand some high-intensity exercises like mountain-climbing. Other abnormalities in the heart, lung, and liver were excluded by physical examination. The Patient-2 had no clinical symptoms of PCD including bronchitis, sinusitis, and otitis media according to her complaint in life history. Chromosomal microarray analysis did not detect potential pathogenic copy number variations in these two affected individuals.
Biallelic mutations in NEK3 in two patients with defective cardiac LR patterning
Through trio-based WES analysis, we identified compound heterozygous non-synonymous mutations (c.805-1G>C/c.1117dupA; c.29T>C/c.356A>G) of NEK3 (NM_002498) (Fig. 1c), encoding a NIMA (never in mitosis A)-related kinase, in two patients. In the patient-1, the MAF (MAF) were 0.00068 for paternally derived c.805-1G>C and 0.00047 for maternally derived c.1117dupA (Table 1). According to GnomAD database, c.805-1G>C and c.1117dupA (p. Thr373fs) were annotated as splicing acceptor site and frameshift variants, respectively. In the patient-2, the MAF were 0.00009 for c.29T>C (p.Ile10Thr) and 0.00048 for c.356A> G (p.His119Arg) (Table 1). These two missense variants p.Ile10Thr and p.His119R were highly conserved in different species including Zebrafish (Danio rerio) and Tropical clawed frog (Xenopus tropicalis) (Fig. 1d), and were uniformly predicted to be damaging by Poly-phen2, SIFT and Mutationtaster (Table 1).
c.805-1G>C leads to exon 10 skipping and generates premature translation termination
To study the effects of splicing acceptor site variant c.805-1G>C on mRNA splicing, we performed real-time PCR using cDNAs derived from fresh peripheral blood samples of patient and his unaffected parents of Family-1. Two different pairs of primers were designed for PCR amplification as indicated (Fig. 2a). We did not identify the difference on the size of PCR products when using Primer-1. However, a smaller size of PCR band that uniquely exists in the cDNA samples of patient and his father that carrying c.805-1G>C was identified when using Primer-2 (Fig. 2b). Sanger sequencing further demonstrated that exon 10 (73-bp length) of NEK3 was entirely skipped in this smaller band, which would induce an out-of-frame alteration. Thus, both c.805-1G>C (p.Ile269GlnfsTer8) and c.1117dupA (p.Thr373AsnfsTer19) could generate prematurely terminated protein translation (Fig. 2c). As a consequence, protein level of NEK3 were almost abrogated in the nasal mucosa of Patient-1 when comparing with a health male (Fig. 2d).
NEK3 deficiency promotes expression of NNMT and SIRT2, and deacetylates α-tubulin
To mimic the biological consequences of NEK3 biallelic loss-of-function (LoF) mutations on cellular biological processes, we knocked down the expression of NEK3 by using small interfering RNAs (siRNAs). Human TERT RPE cells were transfected with control (Ctrl) or each NEK3 siRNA#1-3. Silencing efficiency was tested by western blots showing that siRNAs #1 and #2 largely depleted NEK3 protein levels (Fig. 3a), whereas siRNA #3 and Ctrl siRNA had no significant effect on NEK3 silencing. Subsequent RNA-seq based transcriptome analysis identified 826 and 701 DEGs in NEK3-silenced hRPE cells mediated by siRNAs #1 and #2, respectively, whereas 36 DEGs were identified in cells transfected with siRNA #3 compared with Ctrl siRNA (Fig. 3b, c). After exclusion of unspecific DEGs generated by Ctrl siRNA or ineffective siRNA#3, a total of 342 DEGs were specially shared by siRNA #1 and #2 mediated NEK3 silencing (Fig. 3d). GO enrichment analysis of these 342 DEGs revealed that sister chromatid segregation (p = 6.19E−8), DNA replication (p = 1.01E−7), and spindle microtubule (p = 1.38E−7) were the top 3 enriched terms (Fig. 3e). This finding suggested NEK3 deficiency might interfere microtubule assembly.
In this study, mRNA levels of both NNMT (nicotinamide N-methyltransferase) and Sirt2 (sirtuin2) were significantly upregulated in the NEK3-silenced cells by transcriptome analysis followed by qPCR validation (Fig. 4a). Previous study found expression of NNMT in human neuroblastoma cells significantly increased the expression of Sirt1/2/3 (ref. 31). Sirt2 activity could be inhibited by the substrate of NNMT, nicotinamide. Sirt2 is the major effector that involving the α-tubulin deacetylation and regulates cilia length32,33. Acetylated microtubules have been considered to be stable, long-lived microtubules and α-tubulin deacetylation is associated with cilia resorption or disassembly34. In addition, a well-defined ciliary protein KIF7 that controls cilium architecture35,36,37 was significantly downregulated by NEK3 knockdown according to transcriptome analysis (Supplementary Table S2). Therefore, we performed immunofluorescence staining for acetylated α-tubulin. By using a confocal microscope, we observed that cilia biogenesis was disrupted by NEK3 knockdown mediated by both siRNA #1 and #2 whereas siRNA#3 had no effect on this process when compared with Ctrl siRNA in human RPE cells. Representative images of Ctrl siRNA, NEK3 siRNA#1, and siRNA#3 were displayed in Fig. 4b. Western blot further confirmed the markedly downregulation of acetylated α-tubulin and HEC1, the main component of the NDC80 kinetochore complex, in NEK3-silenced RPE cells whereas XPO1 and phosphor-Akt (Ser473) were not significantly affected (Fig. 4c).
Defective ciliary ultrastructure in Patient-1 with NEK3 mutations
To study the nature of NEK3-associated structural defects, we established culture of primary human nasal respiratory epithelial cells and identified restricted ciliary motion with shorten forward and recovery stroke in Patient-1 as previously (29,306) (Fig. 5a) We next assessed by TEM the ultrastructural effect of NEK3 mutations on cilia from the nasal mucosa of Patient-1 (Fig. 5b). Overall, we examined >50 cilia cross-sections without bias toward measurements in any particular cilia region. the percentage of abnormal cilia is 43.8% (25/57) in Patient-1 and 1.92% (1/52) in the healthy individual (Fig. 5c). Ultrastructural analyses identified multiple ultrastructural abnormalities including the absence of the central-pair microtubules and nine peripheral microtubule doublets when compared to a healthy male (Fig. 5d). These findings indicated that NEK3 mutations are associated with the defective microtubule assembly and ultrastructural abnormalities in the Patient-1.
NEK3 deficiency downregulates mRNA levels of main components of inner rings of NPC
Inner ring complex of the NPC was constituted by Nup93, Nup205, Nup188, and Nup155 (ref. 38). Inner ring components of the NPC were essential for nuclear assembly and nuclear envelope formation39,40. A recent study found inner ring components of the NPC, but not other nucleoporins, specially localize and function at the cilium base41. Depletion of Nup188 and its binding partner Nup93 leads to a global defect in LR patterning in Xenopus41. Moreover, we previously identified biallelic NUP205 mutations in patients with abnormal cardiac LR patterning and depletion of Nup205 cause abnormal cardiac LR organization in Zebrafish42. In addition, inner ring components of the NPC were essential for nuclear assembly and nuclear envelope formation43,44.
In this study, transcriptome analysis revealed that three components of inner ring of NPC, including NUP205, NUP188, and NUP155 were significantly downregulated in the NEK3-silenced cells (Supplementary Table S3), which explains that chromatin assembly or disassembly pathway was enriched in GO analysis (p = 2.34E−7). Downregulation of NUP205 is further confirmed by western blot (Fig. 6a). To assess the effects of NUP205 deficiency on cellular biological processes, we explored siRNA-mediated NUP205 knockdown in hRPE cells and performed transcriptome analysis. Silencing efficiency was tested by real-time qPCR showing that both two siRNAs #1 and #2 effectively knockdown NUP205 expression (Fig. 6b). Very interestingly, RNA-seq based transcriptome analysis revealed that NEK3 is significantly downregulated in NUP205-silenced cells (Fig. 6c) followed by western blot validation (Fig. 6d). This finding implicates the potential reciprocal regulatory interactions between the NEK3 and NUP205 on gene expression.
Discussion
Abnormalities in cilia ultrastructure lead to a range of human phenotypes termed as ciliopathies, which include nephronophthisis, respiratory infection, infertility, retinal degeneration, situs abnormality, and skeletal dysplasia45. Nek kinases were initially studied for their roles in regulating cilia, centrosomes, and microtubule dynamics during cell cycle. Many NEKs including NEK1/2/8/9/10 have found to be mutated in human ciliopathies21,22,46,47,48. Here, we firstly identified biallelic LoF mutations of NEK3 in two unrelated patients with cilia-related abnormal cardiac LR patterning. Previous study found that NEK3 was potentially involved in actin polymerization by regulating the levels of F-action (polymerized form), which is essential for prolactin-mediated cytoskeletal reorganization27.
NEK3 have two functional alleles in humans via alternative in-frame splicing of exon 11. In the full-length allele-1 (NM_002498), c.805-1G>C would induce an out-of-frame deletion as shown in the Patient-1. However, in the shortened allele-2 (NM_001146099) in which the length of exon 10 is 72 bp, c.805-1G>C have no effect on the reading frame and just induce a non-frameshift deletion, which possibly still expresses a functional protein. This finding suggests the genetic consequence of the c.805-1G>C of NEK3 on disease occurrence could be altered in different allele.
HDAC6 and SIRT2 were two major effectors that deacetylates α-tubulin32,49. Previous studies have found that Nek3 regulates microtubule via HDAC6 acetylation in neurons involving axonal degeneration26. Given the similarities between cilia and neuronal axons (both are microtubule-rich structures extending from the cell body), we initiated studies and found cilia biogenesis in human RPE cell was disrupted by NEK3 deficiency via SIRT2 deacetylation. In our study, mRNA levels of HDAC6 and other HDAC members were not significantly changed in NEK3-silenced cells. Intriguingly, SIRT2 ablation was reported to have no effect on tubulin acetylation in brain or the progression of Huntington’s disease phenotypes50. These findings further suggested SIRT2 and HDAC6 exhibited potential tissue specificity on microtubule deacetylation. Moreover, overexpression of SIRT2 in renal epithelial cells was reported to cause decreased cilia length and morpholino knockdown of SIRT2 increased cilia length in Zebrafish51,52, which implicates the regulatory role of SIRT2 on cilia formation.
In this study, NNMT was the one of most upregulated DEGs in NEK3-silenced cells. NNMT catalyzes the reaction between nicotinamide and S-adenosylmethionine to produce 1-methylnicotinamide and S-adenosylhomocysteine. NNMT could promote the expression and stability of SIRTs31,53, and be inhibited by substrate of NNMT, nicotinamide51. These findings suggested upregulated NNMT might promote the expression of SIRT2 in the NEK3-silenced cells. In addition, NNMT is overexpressed in various human tumors and involved in the development and progression of several carcinomas54.
Both mRNA and protein levels of Hec1 were markedly downregulated in the NEK3-depleted cells. The outer kinetochore complex component Ndc80/Hec1 plays a critical role in regulating microtubule attachment to the spindle for accurate sister chromatid segregation55,56. GO analysis demonstrated that sister chromatid segregation is the most enriched term by DEGs of transcriptome analysis in our study. Very interestingly, a controversial model was recently proposed that biased chromatid segregation was involved in the LR symmetry breaking in mice without directly involving the Nodal gene57.
We identified the potential reciprocal regulatory interactions between the NEK3 and NUP205. Previous study found phosphorylation of NUP98 by multiple kinases, mainly NEK2/6/7, is crucial for NPC disassembly during mitotic entry58. Thus, investigating the effect of NEK3 site-specific point mutation on phosphorylation of inner ring components of the NPC like NUP205 and NUP188 are a potential avenue for further mechanistic understanding of these observations. Overall, our study shows that recessive variants in NEK3 are associated with abnormal cardiac LR patterning. Meanwhile, NEK3 could be a candidate gene for human ciliopathies.
Data availability
Raw sequencing data have been deposited in the NCBI Sequence Read Archive (SRA) database under accession number SRP281163. Other data that support the findings of this study are available from the corresponding authors upon reasonable request.
References
Hamada, H., Meno, C., Watanabe, D. & Saijoh, Y. Establishment of vertebrate left-right asymmetry. Nat. Rev. Genet. 3, 103–113 (2002).
Sutherland, M. J. & Ware, S. M. Disorders of left-right asymmetry: heterotaxy and situs inversus. Am. J. Med. Genet. C Semin. Med. Genet. 151C, 307–317 (2009).
Bartoloni, L. et al. Mutations in the DNAH11 (axonemal heavy chain dynein type 11) gene cause one form of situs inversus totalis and most likely primary ciliary dyskinesia. Proc. Natl Acad. Sci. USA 99, 10282–10286 (2002).
Leigh, M. W. et al. Clinical and genetic aspects of primary ciliary dyskinesia/Kartagener syndrome. Genet. Med. 11, 473–487 (2009).
Pereira, R. et al. Clinical and genetic analysis of children with kartagener syndrome. Cells 8, 900 (2019).
Vetrini, F. et al. Bi-allelic mutations in PKD1L1 are associated with laterality defects in humans. Am. J. Hum. Genet. 99, 886–893 (2016).
Lin, A. E., Ticho, B. S., Houde, K., Westgate, M. N. & Holmes, L. B. Heterotaxy: associated conditions and hospital-based prevalence in newborns. Genet. Med. 2, 157–172 (2000).
Taketazu, M., Lougheed, J., Yoo, S. J., Lim, J. S. & Hornberger, L. K. Spectrum of cardiovascular disease, accuracy of diagnosis, and outcome in fetal heterotaxy syndrome. Am. J. Cardiol. 97, 720–724 (2006).
Peeters, H. & Devriendt, K. Human laterality disorders. Eur. J. Med. Genet. 49, 349–362 (2006).
Lowe, L. A. et al. Conserved left-right asymmetry of nodal expression and alterations in murine situs inversus. Nature 381, 158–161 (1996).
Hirokawa, N., Tanaka, Y., Okada, Y. & Takeda, S. Nodal flow and the generation of left-right asymmetry. Cell 125, 33–45 (2006).
Horani, A. & Ferkol, T. W. Advances in the genetics of primary ciliary dyskinesia: clinical implications. Chest 154, 645–652 (2018).
Fliegauf, M., Benzing, T. & Omran, H. When cilia go bad: cilia defects and ciliopathies. Nat. Rev. Mol. Cell Biol. 8, 880–893 (2007).
Fry, A. M., O’Regan, L., Sabir, S. R. & Bayliss, R. Cell cycle regulation by the NEK family of protein kinases. J. Cell Sci. 125, 4423–4433 (2012).
Peres, de et al. Checking NEKs: overcoming a bottleneck in human diseases. Molecules 25, 1778 (2020).
Fry, A. M., Meraldi, P. & Nigg, E. A. A centrosomal function for the human Nek2 protein kinase, a member of the NIMA family of cell cycle regulators. EMBO J. 17, 470–481 (1998).
Mahjoub, M. R., Trapp, M. L. & Quarmby, L. M. NIMA-related kinases defective in murine models of polycystic kidney diseases localize to primary cilia and centrosomes. J. Am. Soc. Nephrol. 16, 3485–3489 (2005).
Shalom, O., Shalva, N., Altschuler, Y. & Motro, B. The mammalian Nek1 kinase is involved in primary cilium formation. FEBS Lett. 582, 1465–1470 (2008).
Rapley, J. et al. The NIMA-family kinase Nek6 phosphorylates the kinesin Eg5 at a novel site necessary for mitotic spindle formation. J. Cell Sci. 121, 3912–3921 (2008).
Holland, P. M. et al. Purification, cloning, and characterization of Nek8, a novel NIMA-related kinase, and its candidate substrate Bicd2. J. Biol. Chem. 277, 16229–16240 (2002).
Chivukula, R. R. et al. A human ciliopathy reveals essential functions for NEK10 in airway mucociliary clearance. Nat. Med. 26, 244–251 (2020).
Fakhro, K. A. et al. Rare copy number variations in congenital heart disease patients identify unique genes in left-right patterning. Proc. Natl Acad. Sci. USA 108, 2915–2920 (2011).
Endicott, S. J., Basu, B., Khokha, M. & Brueckner, M. The NIMA-like kinase Nek2 is a key switch balancing cilia biogenesis and resorption in the development of left-right asymmetry. Development 142, 4068–4079 (2015).
Manning, D. K. et al. Loss of the ciliary kinase Nek8 causes left-right asymmetry defects. J. Am. Soc. Nephrol. 24, 100–112 (2013).
Tanaka, K. & Nigg, E. A. Cloning and characterization of the murine Nek3 protein kinase, a novel member of the NIMA family of putative cell cycle regulators. J. Biol. Chem. 274, 13491–13497 (1999).
Chang, J., Baloh, R. H. & Milbrandt, J. The NIMA-family kinase Nek3 regulates microtubule acetylation in neurons. J. Cell Sci. 122, 2274–2282 (2009).
Miller, S. L., Antico, G., Raghunath, P. N., Tomaszewski, J. E. & Clevenger, C. V. Nek3 kinase regulates prolactin-mediated cytoskeletal reorganization and motility of breast cancer cells. Oncogene 26, 4668–4678 (2007).
Harrington, K. M. & Clevenger, C. V. Identification of NEK3 kinase threonine 165 as a novel regulatory phosphorylation site that modulates focal adhesion remodeling necessary for breast cancer cell migration. J. Biol. Chem. 291, 21388–21406 (2016).
Guo, Z., Chen, W., Huang, J., Wang, L. & Qian, L. Clinical and genetic analysis of patients with primary ciliary dyskinesia caused by novel DNAAF3 mutations. J. Hum. Genet. 64, 711–719 (2019).
Guo, Z., Chen, W., Wang, L. & Qian, L. Clinical and genetic spectrum of children with primary ciliary dyskinesia in China. J. Pediatr. 225, 157–165 (2020).
Liu, K. Y. et al. Nicotinamide N-methyltransferase increases complex I activity in SH-SY5Y cells via sirtuin 3. Biochem. Biophys. Res. Commun. 467, 491–496 (2015).
North, B. J., Marshall, B. L., Borra, M. T., Denu, J. M. & Verdin, E. The human Sir2 ortholog, SIRT2, is an NAD+-dependent tubulin deacetylase. Mol. Cell 11, 437–444 (2003).
Skoge, R. H., Dölle, C. & Ziegler, M. Regulation of SIRT2-dependent α-tubulin deacetylation by cellular NAD levels. DNA Repair (Amst.) 23, 33–38 (2014).
Eshun-Wilson, L. et al. Effects of α-tubulin acetylation on microtubule structure and stability. Proc. Natl Acad. Sci. USA 116, 10366–10371 (2019).
Liem, K. F. Jr., He, M., Ocbina, P. J. & Anderson, K. V. Mouse Kif7/Costal2 is a cilia-associated protein that regulates Sonic hedgehog signaling. Proc. Natl Acad. Sci. USA 106, 13377–13382 (2009).
He, M. et al. The kinesin-4 protein Kif7 regulates mammalian Hedgehog signalling by organizing the cilium tip compartment. Nat. Cell Biol. 16, 663–672 (2014).
Jiang, S. et al. Interplay between the kinesin and tubulin mechanochemical cycles underlies microtubule tip tracking by the non-motile ciliary kinesin kif7. Dev. Cell 49, 711–730.e8 (2019).
Burdine, R. D., Preston, C. C., Leonard, R. J., Bradley, T. A. & Faustino, R. S. Nucleoporins in cardiovascular disease. J. Mol. Cell. Cardiol. 141, 43–52 (2020).
Ródenas, E., Klerkx, E. P., Ayuso, C., Audhya, A. & Askjaer, P. Early embryonic requirement for nucleoporin Nup35/NPP-19 in nuclear assembly. Dev. Biol. 327, 399–409 (2009).
Franz, C. et al. Nup155 regulates nuclear envelope and nuclear pore complex formation in nematodes and vertebrates. EMBO J. 24, 3519–3531 (2005).
Del, Viso, F. et al. Congenital heart disease genetics uncovers context-dependent organization and function of nucleoporins at cilia. Dev. Cell 38, 478–492 (2016).
Chen, W. et al. Bi-allelic mutations in NUP205 and NUP210 are associated with abnormal cardiac left-right patterning. Circ. Genom. Precis. Med. 12, e002492 (2019).
De-Magistris, P., Tatarek-Nossol, M., Dewor, M. & Antonin, W. A self-inhibitory interaction within Nup155 and membrane binding are required for nuclear pore complex formation. J. Cell. Sci. 131, jcs208538 (2018).
Hawryluk-Gara, L. A., Shibuya, E. K. & Wozniak, R. W. Vertebrate Nup53 interacts with the nuclear lamina and is required for the assembly of a Nup93-containing complex. Mol. Biol. Cell. 16, 2382–2394 (2005).
Badano, J. L., Mitsuma, N., Beales, P. L. & Katsanis, N. The ciliopathies: an emerging class of human genetic disorders. Annu. Rev. Genomics Hum. Genet. 7, 125–148 (2006).
Thiel, C. et al. NEK1 mutations cause short-rib polydactyly syndrome type majewski. Am. J. Hum. Genet 88, 106–114 (2011).
Casey, J. P. et al. Recessive NEK9 mutation causes a lethal skeletal dysplasia with evidence of cell cycle and ciliary defects. Hum. Mol. Genet. 25, 1824–1835 (2016).
Grampa, V. et al. Novel nek8 mutations cause severe syndromic renal cystic dysplasia through yap dysregulation. PLoS Genet. 12, e1005894 (2016).
Hubbert, C. et al. HDAC6 is a microtubule-associated deacetylase. Nature 417, 455–458 (2002).
Bobrowska, A., Donmez, G., Weiss, A., Guarente, L. & Bates, G. SIRT2 ablation has no effect on tubulin acetylation in brain, cholesterol biosynthesis or the progression of Huntington’s disease phenotypes in vivo. PLoS ONE 7, e34805 (2012).
Zhou, X. et al. SIRT2 regulates ciliogenesis and contributes to abnormal centrosome amplification caused by loss of polycystin-1. Hum. Mol. Genet. 23, 1644–1655 (2014).
Lim, J., Son, J., Ryu, J. & Kim, J. E. Sirt2 affects primary cilia formation by regulating mtor signaling in retinal pigmented epithelial cells. Int. J. Mol. Sci. 21, 2240 (2020).
Hong, S. et al. Nicotinamide N-methyltransferase regulates hepatic nutrient metabolism through Sirt1 protein stabilization. Nat. Med. 21, 887–894 (2015).
Eckert, M. A. et al. Proteomics reveals NNMT as a master metabolic regulator of cancer-associated fibroblasts. Nature 569, 723–728 (2019).
Martin-Lluesma, S., Stucke, V. M. & Nigg, E. A. Role of Hec1 in spindle checkpoint signaling and kinetochore recruitment of Mad1/Mad2. Science 297, 2267–2270 (2002).
Tauchman, E. C., Boehm, F. J. & DeLuca, J. G. Stable kinetochore-microtubule attachment is sufficient to silence the spindle assembly checkpoint in human cells. Nat. Commun. 6, 10036 (2015).
Sauer, S. & Klar, A. J. Left-right symmetry breaking in mice by left-right dynein may occur via a biased chromatid segregation mechanism, without directly involving the Nodal gene. Front. Oncol. 2, 166 (2012).
Laurell, E. et al. Phosphorylation of Nup98 by multiple kinases is crucial for NPC disassembly during mitotic entry. Cell 144, 539–550 (2011).
Acknowledgements
We are deeply grateful to all the patients and their families for their participation. This work was supported by grants from National Key Research and Development Project (grant no. 2018YFA0801300), Shanghai Rising-Star Program (grant no. 18QA1403500), the National Natural Science Foundation of China (grant no. 81870285, 81741097, and 81601283), the Shanghai Science and Technology International Exchange and Cooperation Fund (grant no.18410741900) to X.Z.
Author information
Authors and Affiliations
Contributions
X.Z. designed the project, analyzed the data, and wrote the manuscript. Y.Z. performed PCR amplification, immunofluorescence assay, and western blots. W.C. performed TEM analysis and interpreted the clinical information. W.Z. and Z.L. performed real-time qPCR and Sanger sequencing, respectively. All authors read and approved the final version of the manuscript.
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Edited by S. Lavandero
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zhang, Y., Chen, W., Zeng, W. et al. Biallelic loss of function NEK3 mutations deacetylate α-tubulin and downregulate NUP205 that predispose individuals to cilia-related abnormal cardiac left–right patterning. Cell Death Dis 11, 1005 (2020). https://doi.org/10.1038/s41419-020-03214-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41419-020-03214-1
This article is cited by
-
Biallelic mutations of TTC12 and TTC21B were identified in Chinese patients with multisystem ciliopathy syndromes
Human Genomics (2022)
-
Mutation in XPO5 causes adult-onset autosomal dominant familial focal segmental glomerulosclerosis
Human Genomics (2022)
-
Acox2 is a regulator of lysine crotonylation that mediates hepatic metabolic homeostasis in mice
Cell Death & Disease (2022)
-
Biallelic DNAH9 mutations are identified in Chinese patients with defective left–right patterning and cilia-related complex congenital heart disease
Human Genetics (2022)