Rationale: Pulmonary nontuberculous mycobacterial disease (PNTM) often affects white postmenopausal women, with a tall and lean body habitus and higher rates of scoliosis, pectus excavatum, mitral valve prolapse, and mutations in the CFTR gene. These clinical features and the familial clustering of the disease suggest an underlying genetic mechanism.
Objectives: To map the genes associated with PNTM, whole-exome sequencing was conducted in 12 PNTM families and 57 sporadic cases recruited at the National Institutes of Health Clinical Center during 2001–2013.
Methods: We performed a variant-level and a gene-level parametric linkage analysis on nine PNTM families (16 affected and 20 unaffected) as well as a gene-level association analysis on nine PNTM families and 55 sporadic cases.
Measurements and Main Results: The genome-wide variant-level linkage analysis using 4,328 independent common variants identified a 20-cM region on chromosome 6q12-6q16 (heterogeneity logarithm of odds score = 3.9), under a recessive disease model with 100% penetrance and a risk allele frequency of 5%. All genes on chromosome 6 were then tested in the gene-level linkage analysis, using the collapsed haplotype pattern method. The TTK protein kinase gene (TTK) on chromosome 6q14.1 was the most significant (heterogeneity logarithm of odds score = 3.38). In addition, the genes MAP2K4, RCOR3, KRT83, IFNLR1, and SLC29A1 were associated with PNTM in our gene-level association analysis.
Conclusions: The TTK gene encodes a protein kinase that is essential for mitotic checkpoints and the DNA damage response. TTK and other genetic loci identified in our study may contribute to the increased susceptibility to NTM infection and its progression to pulmonary disease.
Pulmonary nontuberculous mycobacterial disease (PNTM) has posed an increasing health burden globally in aging populations. The shared morphologic and clinical features among individuals affected with PNTM suggest a strong genetic influence. Several studies have connected the disease to multiple gene categories and pathways. However, no specific genetic regions or loci have been identified for PNTM.
We found evidence of genetic linkage on chromosome 6q12-q16 and further identified the gene TTK (TTK protein kinase) as a candidate gene for PNTM. TTK may contribute to the susceptibility of PNTM through its role in host DNA damage repair and cell survival.
Pulmonary nontuberculous mycobacterial (PNTM) disease is caused by infections with nontuberculous mycobacteria, which are ubiquitous in the environment and normal inhabitants of the drinking water system. Although exposure is nearly universal and infection is common in many geographic areas, only a proportion of infected people develop pulmonary disease (1–4). The disease is usually chronic, indolent, and difficult to diagnose and treat. According to the 2007 American Thoracic Society/Infectious Diseases Society of America statement (5), the diagnosis of PNTM is made from specific clinical, microbiological, and chest imaging findings. Population-based data from North America (3, 4, 6–8), Europe (9–13), Australia (14), and Asia (15, 16) have consistently shown a continued increase in both the prevalence and incidence of PNTM in the past two decades. Across these regions, PNTM is more common in individuals aged 50 years and older, with Mycobacterium avium complex as the dominant species of PNTM. In the United States, the majority of PNTM cases are white females (3, 4, 7, 8), whereas in Europe middle-aged men with chronic lung conditions predominate, possibly due to a historical excess of men smoking (9, 10).
Host and environmental factors interact to influence PNTM disease risk. It is associated with predisposing conditions including lung cancer, chronic obstructive pulmonary disease (COPD), and cystic fibrosis (17). In addition, a low body mass index has been associated with increased risk of PNTM disease (16, 18, 19). Climatic and population factors, such as greater evapotranspiration and greater population density, are thought to promote mycobacterial growth and disease transmission (17, 20). However, given that the organisms are so widespread yet disease is rare, host susceptibility likely plays a key role.
A distinct group of people affected by PNTM in the United States are postmenopausal women who are white, tall, thin, and nonsmokers, who have high rates of scoliosis, pectus excavatum, mitral valve prolapse, and cystic fibrosis transmembrane conductance regulator (CFTR) mutations (18). PNTM often presents among these women as nodular bronchiectatic PNTM, as opposed to the cavitary PNTM typically seen in men with a history of smoking and COPD (20). The shared morphologic features of women with PNTM, along with the observed familial clustering (19), support the hypothesis of genetic predisposition to the development of PNTM. Candidate gene studies (21–23) and pathway/network analyses (24, 25) have implicated certain genetic factors and biological pathways for PNTM.
In this study, we conducted a genome-wide linkage analysis among nine PNTM families to identify genetic regions that may harbor PNTM genes. With whole-exome sequencing data, we performed a gene-level linkage test, using a collapsed haplotype pattern (CHP) method, to further determine the gene(s) driving the detected linkage signals. A gene-based association analysis was also conducted on PNTM families and sporadic cases to identify rare variants associated with PNTM. Some of the results of this study have been previously reported in the form of an abstract (26).
A detailed description of the methods can be found in the online supplement.
Between 2001 and 2013, 12 PNTM families (21 PNTM-affected and 24 unaffected family members) and 57 sporadic PNTM cases were recruited at the National Institutes of Health Clinical Center (Bethesda, MD). PNTM diagnoses were made according to clinical guidelines (see the online supplement). Informed consent was obtained from all subjects, and this study was approved by the National Institutes of Health (protocols 01-I-0202, 07-I-0033, 07-I-0142, and 09-I-0172).
Whole-exome sequencing and variants calling have been previously reported (25). All variants were annotated by ANNOVAR (March 2015) (27). After quality controls, the final data set contained nine families (16 PNTM-affected and 20 unaffected individuals) and 55 sporadic cases with 249,843 single-nucleotide variants (SNVs) and 30,665 small insertion/deletion variants (INDELs). Our study population is primarily of European ancestry with a small number of individuals with East Asian ancestry (see the online supplement).
We conducted a parametric linkage analysis of the PNTM families at both the variant level and the gene level. The linkage analysis was limited to SNVs only, because the genetic distance between insertion/deletion variants in the human genome is not well characterized.
The variant-level parametric linkage analysis was performed in MERLIN (28), using a subset of common (minor allele frequency [MAF] > 0.10) and independent (pairwise r2 < 0.10) variants with call rate greater than 0.95. The genetic map GRCh37, from the HapMap project (https://www.genome.gov/10001688/international-hapmap-project/), was used for the analysis. We tested 60 disease models varying by mode of inheritance (dominant, recessive, and additive), penetrance (40%, 60%, 80%, and 100%), and risk allele frequency (0.01, 0.05, 0.10, 0.15, and 0.20). A heterogeneity logarithm of odds (HLOD) score greater than 3.3 was considered significant for genetic linkage (see the online supplement).
For chromosomal regions with the most significant linkage, a gene-level linkage analysis was performed to determine the gene(s) driving the linkage signals. We first applied the CHP method (29) to all sequenced variants to construct haplotypes for each gene and tested these genes for linkage, using MERLIN (see the online supplement). To infer the MAF range of the causal variants, we conducted a separate analysis to include only the rare variants with MAF not greater than 0.05 and the common variants with MAF greater than 0.05.
To search for rare variants associated with PNTM, we performed a genome-wide association test in our study cohort of PNTM families and sporadic cases. We applied a gene-based approach that extended the commonly used collapsing and kernel methods in unrelated case–control data to account for known pedigree relationships (30). East Asian individuals were excluded from the analysis, resulting in 66 case subjects and 15 control subjects in the analyzed cohort. Analysis was restricted to 16,026 genes with at least two rare variants defined as SNVs with call rate greater than 95% and MAF less than 0.05 in the 1000 Genomes Project European populations (see the online supplement). The widely accepted genome-wide threshold P < 5 × 10−8 was used to determine the significance.
A set of 4,328 independent SNVs with known genetic location was tested in the variant-level parametric linkage analysis for nine families (see Figure E1 in the online supplement). This analysis identified a 20-cM region on chromosome 6q12-q16 with evidence of linkage, and the evidence was strongest (HLOD score = 3.9) under a recessive disease model with 100% penetrance and a risk allele frequency of 0.05 (Figure 1; and see Figures E3 and E4). An empirical P value obtained from a permutation analysis under the null hypothesis of no linkage was less than 10−6. Under this disease model, family A (LOD score = 0.963), family B (LOD score = 0.977), and family D (LOD score = 0.711) each contributed the most to this linkage signal on chromosome 6 (see Figure E1). The information content across markers on chromosome 6 was consistently greater than 80%, suggesting an adequate power of the selected markers in extracting inheritance information (see Figure E3) (31).
In gene-level linkage analysis, considering all SNVs on chromosome 6, a total of 1,005 genes were collapsed using the CHP method and tested in a single-point parametric linkage analysis under various recessive disease models with 100% penetrance, varying by risk allele frequencies. The TTK protein kinase gene (TTK; OMIM [Online Mendelian Inheritance in Man] 604092) on chromosome 6q14.1 was identified as significant with an HLOD score of 3.38. The significance of the TTK gene was attenuated when common variants (MAF > 0.05) were excluded from the analysis (Table 1), suggesting that the risk alleles in the TTK gene associated with PNTM are likely to have a frequency between 0.05 and 0.10 in European populations.
Gene | Location† | Disease Model‡ | All* | Common* (>0.05) | Rare* (≤0.05) | |||
---|---|---|---|---|---|---|---|---|
No. of SNVs | HLOD§ | No. of SNVs | HLOD§ | No. SNVs | HLOD§ | |||
TTK | 6q14.1 | Recessive, 100% penetrance, and risk allele frequency of 0.05 | 27 | 3.38 | 13 | 3.35 | 14 | 1.30 |
TPBG | 6q14-6q15 | Recessive, 100% penetrance, and risk allele frequency of 0.05 | 6 | 3.08 | 1 | 2.49 | 5 | 0.41 |
ORC3 | 6q14.3-q16.1 | Recessive, 100% penetrance, and risk allele frequency of 0.20 | 18 | 3.10 | 10 | 3.18 | 8 | 0.33 |
ANKRD6 | 6q15 | Recessive, 100% penetrance, and risk allele frequency of 0.05 | 22 | 3.02 | 12 | 2.52 | 10 | 0.44 |
Among 27 SNVs and 4 INDELs sequenced in the TTK gene, 22 of them showed variations across the PNTM families (see Table E1). As suggested by the linkage results, the disease-causing variants may cosegregate in the PNTM families as either recessive homozygotes or compound heterozygotes. We applied a filtering strategy to search for the causal variants for PNTM (see Figure E5). Eight variants met the criteria of not being homozygotes of the risk allele in any unaffected individuals across all PNTM families. Variants were then considered recessive if they were homozygotes of the risk allele in all cases of a given family, allowing families to be affected by different variants. Variants were considered compound heterozygous if they were (1) heterozygotes in all cases of a given family, and (2) not compound heterozygotes in any control subjects of any PNTM family. However, we found that no variants in the TTK gene region met the above-mentioned criteria (data not shown).
Evidence of suggestive linkage was also found in the linkage region on chromosome 6q12-6q16 for the genes trophoblast glycoprotein (TPBG [OMIM 190920]; HLOD score = 3.08), origin recognition complex subunit 3 (ORC3 [OMIM 604972]; HLOD score = 3.10), and ankyrin repeat domain 6 (ANKRD6 [OMIM 610583]; HLOD = 3.02). For all of these genes, the linkage signals were likely to be driven by variants with common allele frequencies greater than 0.05 (Table 1).
A total of 16,026 genes across the genome were tested in the gene-based association analysis, with a mean coverage of seven SNVs in each gene region. Five genes reached genome-wide significance at P < 5 × 10−8 (Table 2).
Multipoint parametric linkage analysis identified a 20-cM region on chromosome 6q12-q16 under a recessive disease model with 100% penetrance and risk allele frequency of 0.05. The linkage was most significant at the TTK protein kinase gene (TTK), as suggested in our gene-level linkage analysis using the CHP method. In the same linkage region, suggestive linkage was also found with genes TPBG, ORC3, and ANKRD6.
The TTK gene on 6q14.1 encodes a dual-specificity protein kinase (hMps1), a key spindle assembly checkpoint protein that regulates the proper chromosomal alignment and segregation during mitosis (32–35). The expression level of hMps1 in several human tumors is correlated positively with tumor grade and negatively with survival (36–38), making it one of the most promising drug targets for cancer treatment (39). In addition to regulation of the spindle assembly checkpoint, hMps1 also participates in the DNA damage response. Depletion of hMps1 causes a defect in G2/M arrest, compromising DNA repair and cell survival (40–42). Dysfunction in DNA repair is a common finding in patients with COPD and tuberculosis (43). As studies on TTK/hMps1 have been exclusively focused on its tumorigenesis, its contribution to bronchiectasis; to bacterial, mycobacterial, or fungal infections; or overall progression to pulmonary disease is currently unknown. Impaired DNA damage repair in host cells induced by defects in the TTK gene may provide a beneficial environment for bacterial proliferation or may compromise local innate immune responses, allowing infections to cause disease.
Studies on the genes TPBG, ORC3, and ANKRD6 are limited and thus their relevance to PNTM is unclear. All three genes are expressed in lung tissues. The TPBG gene encodes an oncofetal antigen that is highly expressed in many types of cancer (44, 45). The protein encoded by the ORC3 gene is a subunit of the origin recognition complex (ORC), which is essential for the initiation of the DNA replication. Diversin, a planar cell polarity protein encoded by the ANKRD6 gene, is involved in the Wnt signaling pathways and has been implicated in human neural tube defects (46, 47).
Our gene-based association analysis identified five genes whose rare variants are associated with PNTM, of which two genes may affect PNTM through their involvement in host immune defense in bacterial or viral infection. The most significant association was found with the mitogen-activated protein kinase 4 gene (MAP2K4) on chromosome 17p12. The encoded protein is a member of the mitogen-activated protein kinase (MAPK) family, and is involved in a number of KEGG pathways including the Toll-like receptor signaling pathway; epithelial cell signaling in Helicobacter pylori infection; and hepatitis B, influenza A, and Epstein–Barr virus infection (48). The product of the IFN-λ receptor 1 gene (IFNLR1) on chromosome 1p36.11 forms a receptor complex with IL-10 receptor β (IL10RB) and interacts with IL-28A, IL-28B, and IL-29. Genetic variations in IFNLR1 are correlated with susceptibility to infection with hepatitis C virus and spontaneous viral clearance (49). Differential expression of IFNLR1 in low-risk and high-risk human papillomavirus–positive women suggests its potential contribution to lesion progression in high-risk human papillomavirus infection (50). In addition, the solute carrier family 29 member 1 gene (SLC29A1) encodes the human equilibrative nucleoside transporter 1 (hENT1) that mediates the transport of nucleosides and nucleoside analog drugs across the plasma membrane. Repression of hENT1 dampens pulmonary inflammation and improves lung function in animal models (51), which could potentially be relevant to progression to pulmonary disease in PNTM.
PNTM is an emerging public health issue, posing substantial health and economic burdens, particularly in aging populations. Our study systematically investigated the host genetics associated with nontuberculous mycobacterial infection and pulmonary disease. We combined traditional linkage analysis of individual variants with a novel gene-level approach based on collapsed haplotypes, and identified the TTK gene on chromosome 6q14.1 as a candidate gene for PNTM. However, our data suggest that the disease-causing variants of the TTK gene are likely common (MAF, 0.05–0.10) in the general population and may function as recessive homozygotes or compound heterozygotes, which is consistent with the observed low prevalence of the disease. Importantly, the causal variants may be located outside the exonic regions of the TTK gene and therefore not be included in this exome-sequencing data set. Furthermore, small insertion/deletion variants may also explain some of the familial aggregation of PNTM, but these currently cannot be analyzed in linkage and gene-based association analysis. Because the gene-based association analysis using pedigree data has slightly inflated type I errors (30), we applied a threshold (P < 5 × 10−8) more stringent than the Bonferroni-corrected P value (P < 3 × 10−6) to determine the statistical significance of the association results. This led to the identification of several genes that may be biologically relevant to PNTM.
Like other complex diseases, PNTM has a multifactorial etiology involving many genes and environmental exposures (25). Our family-based study identified the TTK gene region as linked to the familial form of the disease in our study population. Future studies with larger sample sizes are warranted to replicate our findings. The impact of this gene depends on the frequencies and effect sizes of the putative alleles in the general population, which requires further investigation with population sampling. Functional studies may also help to elucidate the role of TTK and the underlying disease mechanism of PNTM, and hence lead to the development of effective treatments of the disease. Findings from these genetic studies will provide valuable insights for the prevention, diagnosis, and treatment of PNTM.
1. | Winthrop KL, McNelley E, Kendall B, Marshall-Olson A, Morris C, Cassidy M, Saulson A, Hedberg K. Pulmonary nontuberculous mycobacterial disease prevalence and clinical features: an emerging public health disease. Am J Respir Crit Care Med 2010;182:977–982. |
2. | Cassidy PM, Hedberg K, Saulson A, McNelly E, Winthrop KL. Nontuberculous mycobacterial disease prevalence and risk factors: a changing epidemiology. Clin Infect Dis 2009;49:e124–e129. |
3. | Prevots DR, Shaw PA, Strickland D, Jackson LA, Raebel MA, Blosky MA, Montes de Oca R, Shea YR, Seitz AE, Holland SM, et al. Nontuberculous mycobacterial lung disease prevalence at four integrated health care delivery systems. Am J Respir Crit Care Med 2010;182:970–976. |
4. | Henkle E, Hedberg K, Schafer S, Novosad S, Winthrop KL. Population-based incidence of pulmonary nontuberculous mycobacterial disease in Oregon 2007 to 2012. Ann Am Thorac Soc 2015;12:642–647. |
5. | Griffith DE, Aksamit T, Brown-Elliott BA, Catanzaro A, Daley C, Gordin F, Holland SM, Horsburgh R, Huitt G, Iademarco MF, et al.; ATS Mycobacterial Diseases Subcommittee; American Thoracic Society; Infectious Diseases Society of America. An official ATS/IDSA statement: diagnosis, treatment, and prevention of nontuberculous mycobacterial diseases. Am J Respir Crit Care Med 2007;175:367–416. [Published erratum appears in Am J Respir Crit Care Med 2007;175:744–745.] |
6. | Marras TK, Mendelson D, Marchand-Austin A, May K, Jamieson FB. Pulmonary nontuberculous mycobacterial disease, Ontario, Canada, 1998–2010. Emerg Infect Dis 2013;19:1889–1891. |
7. | Mirsaeidi M, Machado RF, Garcia JGN, Schraufnagel DE. Nontuberculous mycobacterial disease mortality in the United States, 1999–2010: a population-based comparative study. PLoS One 2014;9:e91879. |
8. | Adjemian J, Olivier KN, Seitz AE, Holland SM, Prevots DR. Prevalence of nontuberculous mycobacterial lung disease in U.S. Medicare beneficiaries. Am J Respir Crit Care Med 2012;185:881–886. |
9. | Shah NM, Davidson JA, Anderson LF, Lalor MK, Kim J, Thomas HL, Lipman M, Abubakar I. Pulmonary Mycobacterium avium-intracellulare is the main driver of the rise in non-tuberculous mycobacteria incidence in England, Wales and Northern Ireland, 2007–2012. BMC Infect Dis 2016;16:195. |
10. | Russell CD, Claxton P, Doig C, Seagar A-L, Rayner A, Laurenson IF. Non-tuberculous mycobacteria: a retrospective review of Scottish isolates from 2000 to 2010. Thorax 2014;69:593–595. |
11. | Rindi L, Garzelli C. Increase in non-tuberculous mycobacteria isolated from humans in Tuscany, Italy, from 2004 to 2014. BMC Infect Dis 2016;16:44. |
12. | Ringshausen FC, Wagner D, de Roux A, Diel R, Hohmann D, Hickstein L, Welte T, Rademacher J. Prevalence of nontuberculous mycobacterial pulmonary disease, Germany, 2009–2014. Emerg Infect Dis 2016;22:1102–1105. |
13. | Ringshausen FC, Apel R-M, Bange F-C, de Roux A, Pletz MW, Rademacher J, Suhling H, Wagner D, Welte T. Burden and trends of hospitalisations associated with pulmonary non-tuberculous mycobacterial infections in Germany, 2005–2011. BMC Infect Dis 2013;13:231. |
14. | Thomson RM; NTM Working Group at Queensland TB Control Centre and Queensland Mycobacterial Reference Laboratory. Changing epidemiology of pulmonary nontuberculous mycobacteria infections. Emerg Infect Dis 2010;16:1576–1583. |
15. | Namkoong H, Kurashima A, Morimoto K, Hoshino Y, Hasegawa N, Ato M, Mitarai S. Epidemiology of pulmonary nontuberculous mycobacterial disease, Japan. Emerg Infect Dis 2016;22:1116–1117. |
16. | Ide S, Nakamura S, Yamamoto Y, Kohno Y, Fukuda Y, Ikeda H, Sasaki E, Yanagihara K, Higashiyama Y, Hashiguchi K, et al. Epidemiology and clinical features of pulmonary nontuberculous mycobacteriosis in Nagasaki, Japan. PLoS One 2015;10:e0128304. |
17. | Prevots DR, Marras TK. Epidemiology of human pulmonary infection with nontuberculous mycobacteria: a review. Clin Chest Med 2015;36:13–34. |
18. | Kim RD, Greenberg DE, Ehrmantraut ME, Guide SV, Ding L, Shea Y, Brown MR, Chernick M, Steagall WK, Glasgow CG, et al. Pulmonary nontuberculous mycobacterial disease: prospective study of a distinct preexisting syndrome. Am J Respir Crit Care Med 2008;178:1066–1074. |
19. | Colombo RE, Hill SC, Claypool RJ, Holland SM, Olivier KN. Familial clustering of pulmonary nontuberculous mycobacterial disease. Chest 2010;137:629–634. |
20. | McShane PJ, Glassroth J. Pulmonary disease due to nontuberculous mycobacteria: current state and new insights. Chest 2015;148:1517–1527. |
21. | Koh W-J, Kwon OJ, Kim EJ, Lee KS, Ki C-S, Kim JW. NRAMP1 gene polymorphism and susceptibility to nontuberculous mycobacterial lung diseases. Chest 2005;128:94–101. |
22. | Um S-W, Ki C-S, Kwon OJ, Koh W-J. HLA antigens and nontuberculous mycobacterial lung disease in Korean patients. Lung 2009;187:136–140. |
23. | Jang M-A, Kim S-Y, Jeong B-H, Park HY, Jeon K, Kim J-W, Ki C-S, Koh W-J. Association of CFTR gene variants with nontuberculous mycobacterial lung disease in a Korean population with a low prevalence of cystic fibrosis. J Hum Genet 2013;58:298–303. |
24. | Lipner EM, Garcia BJ, Strong M. Network analysis of human genes influencing susceptibility to mycobacterial infections. PLoS One 2016;11:e0146585. |
25. | Szymanski EP, Leung JM, Fowler CJ, Haney C, Hsu AP, Chen F, Duggal P, Oler AJ, McCormack R, Podack E, et al. Pulmonary nontuberculous mycobacterial infection: a multisystem, multigenic disease. Am J Respir Crit Care Med 2015;192:618–628. |
26. | Chen F, Szymanski EP, Olivier KN, Liu X, Tettelin H, Holland SM, Duggal P. Whole exome sequencing and linkage analysis identified TTK as a candidate gene for pulmonary nontuberculous mycobacterial lung disease. Presented at the 2016 Annual Meeting of the American Society of Human Genetics. October 18, 2016, Vancouver, BC, Canada. Abstract/Program #1318. |
27. | Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res 2010;38:e164. |
28. | Abecasis GR, Cherny SS, Cookson WO, Cardon LR. Merlin: rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet 2002;30:97–101. |
29. | Wang GT, Zhang D, Li B, Dai H, Leal SM. Collapsed haplotype pattern method for linkage analysis of next-generation sequence data. Eur J Hum Genet 2015;23:1739–1743. |
30. | Schaid DJ, McDonnell SK, Sinnwell JP, Thibodeau SN. Multiple genetic variant association testing by collapsing and kernel methods with pedigree or population structured data. Genet Epidemiol 2013;37:409–418. |
31. | Evans DM, Cardon LR. Guidelines for genotyping in genomewide linkage studies: single-nucleotide-polymorphism maps versus microsatellite maps. Am J Hum Genet 2004;75:687–692. |
32. | Lee S, Thebault P, Freschi L, Beaufils S, Blundell TL, Landry CR, Bolanos-Garcia VM, Elowe S. Characterization of spindle checkpoint kinase Mps1 reveals domain with functional and structural similarities to tetratricopeptide repeat motifs of Bub1 and BubR1 checkpoint kinases. J Biol Chem 2012;287:5988–6001. |
33. | Morin V, Prieto S, Melines S, Hem S, Rossignol M, Lorca T, Espeut J, Morin N, Abrieu A. CDK-dependent potentiation of MPS1 kinase activity is essential to the mitotic checkpoint. Curr Biol 2012;22:289–295. |
34. | Liu S-T, Chan GKT, Hittle JC, Fujii G, Lees E, Yen TJ. Human MPS1 kinase is required for mitotic arrest induced by the loss of CENP-E from kinetochores. Mol Biol Cell 2003;14:1638–1651. |
35. | von Schubert C, Cubizolles F, Bracher JM, Sliedrecht T, Kops GJPL, Nigg EA. Plk1 and Mps1 cooperatively regulate the spindle assembly checkpoint in human cells. Cell Reports 2015;12:66–78. |
36. | Slee RB, Grimes BR, Bansal R, Gore J, Blackburn C, Brown L, Gasaway R, Jeong J, Victorino J, March KL, et al. Selective inhibition of pancreatic ductal adenocarcinoma cell growth by the mitotic MPS1 kinase inhibitor NMS-P715. Mol Cancer Ther 2014;13:307–315. |
37. | Tannous BA, Kerami M, Van der Stoop PM, Kwiatkowski N, Wang J, Zhou W, Kessler AF, Lewandrowski G, Hiddingh L, Sol N, et al. Effects of the selective MPS1 inhibitor MPS1-IN-3 on glioblastoma sensitivity to antimitotic drugs. J Natl Cancer Inst 2013;105:1322–1331. |
38. | Carter SL, Eklund AC, Kohane IS, Harris LN, Szallasi Z. A signature of chromosomal instability inferred from gene expression profiles predicts clinical outcome in multiple human cancers. Nat Genet 2006;38:1043–1048. |
39. | Dominguez-Brauer C, Thu KL, Mason JM, Blaser H, Bray MR, Mak TW. Targeting mitosis in cancer: emerging strategies. Mol Cell 2015;60:524–536. |
40. | Wei J-H, Chou Y-F, Ou Y-H, Yeh Y-H, Tyan S-W, Sun T-P, Shen C-Y, Shieh S-Y. TTK/hMps1 participates in the regulation of DNA damage checkpoint response by phosphorylating CHK2 on threonine 68. J Biol Chem 2005;280:7748–7757. |
41. | Yeh Y-H, Huang Y-F, Lin T-Y, Shieh S-Y. The cell cycle checkpoint kinase CHK2 mediates DNA damage-induced stabilization of TTK/hMps1. Oncogene 2009;28:1366–1378. |
42. | Yu Z-C, Huang Y-F, Shieh S-Y. Requirement for human Mps1/TTK in oxidative DNA damage repair and cell survival through MDM2 phosphorylation. Nucleic Acids Res 2016;44:1133–1150. |
43. | da Silva ALG, Bresciani MJ, Karnopp TE, Weber AF, Ellwanger JH, Henriques JAP, Valim AR de M, Possuelo LG. DNA damage and cellular abnormalities in tuberculosis, lung cancer and chronic obstructive pulmonary disease. Multidiscip Respir Med 2015;10:38. |
44. | Kloudová K, Hromádková H, Partlová S, Brtnický T, Rob L, Bartůňková J, Hensler M, Halaška MJ, Špíšek R, Fialová A. Expression of tumor antigens on primary ovarian cancer cells compared to established ovarian cancer cell lines. Oncotarget 2016;7:46120–46126. |
45. | He P, Jiang S, Ma M, Wang Y, Li R, Fang F, Tian G, Zhang Z. Trophoblast glycoprotein promotes pancreatic ductal adenocarcinoma cell metastasis through Wnt/planar cell polarity signaling. Mol Med Rep 2015;12:503–509. |
46. | Allache R, Wang M, De Marco P, Merello E, Capra V, Kibar Z. Genetic studies of ANKRD6 as a molecular switch between Wnt signaling pathways in human neural tube defects. Birth Defects Res A Clin Mol Teratol 2015;103:20–26. |
47. | Wu G, Huang X, Hua Y, Mu D. Roles of planar cell polarity pathways in the development of neural [correction of neutral] tube defects. J Biomed Sci 2011;18:66. |
48. | Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res 2016;44:D457–D462. |
49. | Cui Q, Zhang Y-X, Su J, Chen X, Ding K, Lei N, Liu Y, Li J, Zhang Y, Yu R-B. Genetic variation in IL28RA is associated with the outcomes of HCV infection in a high-risk Chinese population. Infect Genet Evol 2011;11:1682–1689. |
50. | Cannella F, Scagnolari C, Selvaggi C, Stentella P, Recine N, Antonelli G, Pierangeli A. Interferon λ 1 expression in cervical cells differs between low-risk and high-risk human papillomavirus–positive women. Med Microbiol Immunol (Berl) 2014;203:177–184. |
51. | Morote-Garcia JC, Köhler D, Roth JM, Mirakaj V, Eldh T, Eltzschig HK, Rosenberger P. Repression of the equilibrative nucleoside transporters dampens inflammatory lung injury. Am J Respir Cell Mol Biol 2013;49:296–305. |
Supported in part with federal funds to the Institute for Genome Sciences from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services, under contract number HHSN272200900009C (S.M.H.; exome sequencing and analyses). F.C. is supported by the Maryland Genetics, Epidemiology, and Medicine Training Program (MD-GEM) sponsored by the Burroughs-Wellcome Fund.
Author Contributions: F.C. and P.D. designed and performed the statistical analysis and wrote the manuscript; E.P.S. assisted with the quality control of sequencing data; K.N.O. treated patients and helped conceive the study; X.L. and H.T. performed exome sequencing, alignment, and variant calling; S.M.H. conceived the project and revised the manuscript. All authors provided feedback on the manuscript.
This article has an online supplement, which is accessible from this issue’s table of contents at www.atsjournals.org
Originally Published in Press as DOI: 10.1164/rccm.201612-2479OC on August 4, 2017
Author disclosures are available with the text of this article at www.atsjournals.org.