Genome sequencing and comprehensive rare-variant analysis of 465 families with neurodevelopmental disorders.

Despite significant progress in unraveling the genetic causes of neurodevelopmental disorders (NDDs), a substantial proportion of individuals with NDDs remain without a genetic diagnosis after microarray and/or exome sequencing. Here, we aimed to assess the power of short-read genome sequencing (GS), complemented with long-read GS, to identify causal variants in participants with NDD from the National Institute for Health and Care Research (NIHR) BioResource project. Short-read GS was conducted on 692 individuals (489 affected and 203 unaffected relatives) from 465 families. Additionally, long-read GS was performed on five affected individuals who had structural variants (SVs) in technically challenging regions, had complex SVs, or required distal variant phasing. Causal variants were identified in 36% of affected individuals (177/489), and a further 23% (112/489) had a variant of uncertain significance after multiple rounds of re-analysis. Among all reported variants, 88% (333/380) were coding nuclear SNVs or insertions and deletions (indels), and the remainder were SVs, non-coding variants, and mitochondrial variants. Furthermore, long-read GS facilitated the resolution of challenging SVs and invalidated variants of difficult interpretation from short-read GS. This study demonstrates the value of short-read GS, complemented with long-read GS, in investigating the genetic causes of NDDs. GS provides a comprehensive and unbiased method of identifying all types of variants throughout the nuclear and mitochondrial genomes in individuals with NDD.

Elucidating the molecular mechanisms associated with TARS2-related mitochondrial disease.

TARS2 encodes human mitochondrial threonyl tRNA-synthetase that is responsible for generating mitochondrial Thr-tRNAThr and clearing mischarged Ser-tRNAThr during mitochondrial translation. Pathogenic variants in TARS2 have hitherto been reported in a pair of siblings and an unrelated patient with an early onset mitochondrial encephalomyopathy and a combined respiratory chain enzyme deficiency in muscle. We here report five additional unrelated patients with TARS2-related mitochondrial diseases, expanding the clinical phenotype to also include epilepsy, dystonia, hyperhidrosis and severe hearing impairment. In addition, we document seven novel TARS2 variants-one nonsense variant and six missense variants-that we demonstrate are pathogenic and causal of the disease presentation based on population frequency, homology modeling and functional studies that show the effects of the pathogenic variants on TARS2 stability and/or function.

De Novo VPS4A Mutations Cause Multisystem Disease with Abnormal Neurodevelopment.

The endosomal sorting complexes required for transport (ESCRTs) are essential for multiple membrane modeling and membrane-independent cellular processes. Here we describe six unrelated individuals with de novo missense variants affecting the ATPase domain of VPS4A, a critical enzyme regulating ESCRT function. Probands had structural brain abnormalities, severe neurodevelopmental delay, cataracts, growth impairment, and anemia. In cultured cells, overexpression of VPS4A mutants caused enlarged endosomal vacuoles resembling those induced by expression of known dominant-negative ATPase-defective forms of VPS4A. Proband-derived fibroblasts had enlarged endosomal structures with abnormal accumulation of the ESCRT protein IST1 on the limiting membrane. VPS4A function was also required for normal endosomal morphology and IST1 localization in iPSC-derived human neurons. Mutations affected other ESCRT-dependent cellular processes, including regulation of centrosome number, primary cilium morphology, nuclear membrane morphology, chromosome segregation, mitotic spindle formation, and cell cycle progression. We thus characterize a distinct multisystem disorder caused by mutations affecting VPS4A and demonstrate that its normal function is required for multiple human developmental and cellular processes.

Evaluation of copy number variants for genetic hearing loss: a review of current approaches and recent findings.

Structural variation includes a change in copy number, orientation, or location of a part of the genome. Copy number variants (CNVs) are a common cause of genetic hearing loss, comprising nearly 20% of diagnosed cases. While large deletions involving the gene STRC are the most common pathogenic CNVs, a significant proportion of known hearing loss genes also contain pathogenic CNVs. In this review, we provide an overview of currently used methods for detection of CNVs in genes known to cause hearing loss including molecular techniques such as multiplex ligation probe amplification (MLPA) and digital droplet polymerase chain reaction (ddPCR), array-CGH and single-nucleotide polymorphism (SNP) arrays, as well as techniques for detection of CNVs using next-generation sequencing data analysis including targeted gene panel, exome, and genome sequencing data. In addition, in this review, we compile published data on pathogenic hearing loss CNVs to provide an up-to-date overview. We show that CNVs have been identified in 29 different non-syndromic hearing loss genes. An understanding of the contribution of CNVs to genetic hearing loss is critical to the current diagnosis of hearing loss and is crucial for future gene therapies. Thus, evaluation for CNVs is required in any modern pipeline for genetic diagnosis of hearing loss.

Novel phosphopantothenoylcysteine synthetase (PPCS) mutations with prominent neuromuscular features: Expanding the phenotypical spectrum of PPCS-related disorders.

Biallelic pathogenic variants in phosphopantothenoylcysteine synthetase, PPCS, are a rare cause of a severe early-onset dilated cardiomyopathy with high morbidity and mortality. To date, only five individuals with PPCS-mutations have been reported. Here, we report a female infant who presented in the neonatal period with hypotonia, a necrotizing myopathy with intermittent rhabdomyolysis and other extracardiac manifestations before developing a progressive and ultimately fatal dilated cardiomyopathy. Gene agnostic trio genome sequencing revealed two rare variants in the PPCS [MIM: 609853] in trans, a previously reported pathogenic c.320_334del p. (Pro107_Ala111del) variant, and a c.613-3C>G intronic variant of uncertain significance. Functional studies confirmed the likely pathogenicity of this variant. Our case provides clinical and histopathological evidence for an associated neuromuscular phenotype not previously recognized and expands the evolving phenotypic spectrum of PPCS-related disorders. We also performed a literature search of all previously published cases and summarize the common features.

Bi-allelic Mutations in NDUFA6 Establish Its Role in Early-Onset Isolated Mitochondrial Complex I Deficiency.

Isolated complex I deficiency is a common biochemical phenotype observed in pediatric mitochondrial disease and often arises as a consequence of pathogenic variants affecting one of the ∼65 genes encoding the complex I structural subunits or assembly factors. Such genetic heterogeneity means that application of next-generation sequencing technologies to undiagnosed cohorts has been a catalyst for genetic diagnosis and gene-disease associations. We describe the clinical and molecular genetic investigations of four unrelated children who presented with neuroradiological findings and/or elevated lactate levels, highly suggestive of an underlying mitochondrial diagnosis. Next-generation sequencing identified bi-allelic variants in NDUFA6, encoding a 15 kDa LYR-motif-containing complex I subunit that forms part of the Q-module. Functional investigations using subjects' fibroblast cell lines demonstrated complex I assembly defects, which were characterized in detail by mass-spectrometry-based complexome profiling. This confirmed a marked reduction in incorporated NDUFA6 and a concomitant reduction in other Q-module subunits, including NDUFAB1, NDUFA7, and NDUFA12. Lentiviral transduction of subjects' fibroblasts showed normalization of complex I. These data also support supercomplex formation, whereby the ∼830 kDa complex I intermediate (consisting of the P- and Q-modules) is in complex with assembled complex III and IV holoenzymes despite lacking the N-module. Interestingly, RNA-sequencing data provided evidence that the consensus RefSeq accession number does not correspond to the predominant transcript in clinically relevant tissues, prompting revision of the NDUFA6 RefSeq transcript and highlighting not only the importance of thorough variant interpretation but also the assessment of appropriate transcripts for analysis.

ExpansionHunter: a sequence-graph-based tool to analyze variation in short tandem repeat regions.

SUMMARY: We describe a novel computational method for genotyping repeats using sequence graphs. This method addresses the long-standing need to accurately genotype medically important loci containing repeats adjacent to other variants or imperfect DNA repeats such as polyalanine repeats. Here we introduce a new version of our repeat genotyping software, ExpansionHunter, that uses this method to perform targeted genotyping of a broad class of such loci. AVAILABILITY AND IMPLEMENTATION: ExpansionHunter is implemented in C++ and is available under the Apache License Version 2.0. The source code, documentation, and Linux/macOS binaries are available at https://github.com/Illumina/ExpansionHunter/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Spinal muscular atrophy diagnosis and carrier screening from genome sequencing data.

PURPOSE: Spinal muscular atrophy (SMA), caused by loss of the SMN1 gene, is a leading cause of early childhood death. Due to the near identical sequences of SMN1 and SMN2, analysis of this region is challenging. Population-wide SMA screening to quantify the SMN1 copy number (CN) is recommended by the American College of Medical Genetics and Genomics. METHODS: We developed a method that accurately identifies the CN of SMN1 and SMN2 using genome sequencing (GS) data by analyzing read depth and eight informative reference genome differences between SMN1/2. RESULTS: We characterized SMN1/2 in 12,747 genomes, identified 1568 samples with SMN1 gains or losses and 6615 samples with SMN2 gains or losses, and calculated a pan-ethnic carrier frequency of 2%, consistent with previous studies. Additionally, 99.8% of our SMN1 and 99.7% of SMN2 CN calls agreed with orthogonal methods, with a recall of 100% for SMA and 97.8% for carriers, and a precision of 100% for both SMA and carriers. CONCLUSION: This SMN copy-number caller can be used to identify both carrier and affected status of SMA, enabling SMA testing to be offered as a comprehensive test in neonatal care and an accurate carrier screening tool in GS sequencing projects.

Hawaiian monk seals exhibit behavioral flexibility when targeting prey of different size and shape.

Animals use diverse feeding strategies to capture and consume prey, with many species switching between strategies to accommodate different prey. Many marine animals exhibit behavioral flexibility when feeding to deal with spatial and temporal heterogeneity in prey resources. However, little is known about flexibility in the feeding behavior of many large marine predators. Here, we documented the feeding behavior and kinematics of the endangered Hawaiian monk seal (Neomonachus schauinslandi, n=7) through controlled feeding trials. Seals were fed multiple prey types (e.g. night smelt, capelin, squid and herring) that varied in size and shape to examine behavioral flexibility in feeding. Hawaiian monk seals primarily used suction feeding (91% of all feeding trials) across all prey types, but biting, specifically pierce feeding, was also observed (9% of all feeding trials). Suction feeding was characterized by shorter temporal events, a smaller maximum gape and gape angle, and a fewer number of jaw motions than pierce feeding; suction feeding kinematic performance was also more variable compared with pierce feeding. Seals showed behavioral flexibility in their use of the two strategies. Suction feeding was used most frequently when targeting small to medium sized prey and biting was used with increasing frequency on larger prey. The feeding kinematics differed between feeding strategies and prey types, showing that Hawaiian monk seals adjusted their behaviors to particular feeding contexts. Hawaiian monk seals are opportunistic marine predators and their ability to adapt their feeding strategy and behavior to specific foraging scenarios allows them to target diverse prey resources.

Swimming and diving energetics in dolphins: a stroke-by-stroke analysis for predicting the cost of flight responses in wild odontocetes.

Exponential increases in hydrodynamic drag and physical exertion occur when swimmers move quickly through water, and underlie the preference for relatively slow routine speeds by marine mammals regardless of body size. Because of this and the need to balance limited oxygen stores when submerged, flight (escape) responses may be especially challenging for this group. To examine this, we used open-flow respirometry to measure the energetic cost of producing a swimming stroke during different levels of exercise in bottlenose dolphins (Tursiops truncatus). These data were then used to model the energetic cost of high-speed escape responses by other odontocetes ranging in mass from 42 to 2738 kg. The total cost per stroke during routine swimming by dolphins, 3.31±0.20 J kg-1 stroke-1, was doubled during maximal aerobic performance. A comparative analysis of locomotor costs (LC; in J kg-1 stroke-1), representing the cost of moving the flukes, revealed that LC during routine swimming increased with body mass (M) for odontocetes according to LC=1.46±0.0005M; a separate relationship described LC during high-speed stroking. Using these relationships, we found that continuous stroking coupled with reduced glide time in response to oceanic noise resulted in a 30.5% increase in metabolic rate in the beaked whale, a deep-diving odontocete considered especially sensitive to disturbance. By integrating energetics with swimming behavior and dive characteristics, this study demonstrates the physiological consequences of oceanic noise on diving mammals, and provides a powerful tool for predicting the biological significance of escape responses by cetaceans facing anthropogenic disturbances.

Arboretum: reconstruction and analysis of the evolutionary history of condition-specific transcriptional modules.

Comparative functional genomics studies the evolution of biological processes by analyzing functional data, such as gene expression profiles, across species. A major challenge is to compare profiles collected in a complex phylogeny. Here, we present Arboretum, a novel scalable computational algorithm that integrates expression data from multiple species with species and gene phylogenies to infer modules of coexpressed genes in extant species and their evolutionary histories. We also develop new, generally applicable measures of conservation and divergence in gene regulatory modules to assess the impact of changes in gene content and expression on module evolution. We used Arboretum to study the evolution of the transcriptional response to heat shock in eight species of Ascomycota fungi and to reconstruct modules of the ancestral environmental stress response (ESR). We found substantial conservation in the stress response across species and in the reconstructed components of the ancestral ESR modules. The greatest divergence was in the most induced stress, primarily through module expansion. The divergence of the heat stress response exceeds that observed in the response to glucose depletion in the same species. Arboretum and its associated analyses provide a comprehensive framework to systematically study regulatory evolution of condition-specific responses.

Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals.

There is growing recognition that mammalian cells produce many thousands of large intergenic transcripts. However, the functional significance of these transcripts has been particularly controversial. Although there are some well-characterized examples, most (>95%) show little evidence of evolutionary conservation and have been suggested to represent transcriptional noise. Here we report a new approach to identifying large non-coding RNAs using chromatin-state maps to discover discrete transcriptional units intervening known protein-coding loci. Our approach identified approximately 1,600 large multi-exonic RNAs across four mouse cell types. In sharp contrast to previous collections, these large intervening non-coding RNAs (lincRNAs) show strong purifying selection in their genomic loci, exonic sequences and promoter regions, with greater than 95% showing clear evolutionary conservation. We also developed a functional genomics approach that assigns putative functions to each lincRNA, demonstrating a diverse range of roles for lincRNAs in processes from embryonic stem cell pluripotency to cell proliferation. We obtained independent functional validation for the predictions for over 100 lincRNAs, using cell-based assays. In particular, we demonstrate that specific lincRNAs are transcriptionally regulated by key transcription factors in these processes such as p53, NFkappaB, Sox2, Oct4 (also known as Pou5f1) and Nanog. Together, these results define a unique collection of functional lincRNAs that are highly conserved and implicated in diverse biological processes.

Gene duplication and the evolution of ribosomal protein gene regulation in yeast.

Coexpression of genes within a functional module can be conserved at great evolutionary distances, whereas the associated regulatory mechanisms can substantially diverge. For example, ribosomal protein (RP) genes are tightly coexpressed in Saccharomyces cerevisiae, but the cis and trans factors associated with them are surprisingly diverged across Ascomycota fungi. Little is known, however, about the functional impact of such changes on actual expression levels or about the selective pressures that affect them. Here, we address this question in the context of the evolution of the regulation of RP gene expression by using a comparative genomics approach together with cross-species functional assays. We show that an activator (Ifh1) and a repressor (Crf1) that control RP gene regulation in normal and stress conditions in S. cerevisiae are derived from the duplication and subsequent specialization of a single ancestral protein. We provide evidence that this regulatory innovation coincides with the duplication of RP genes in a whole-genome duplication (WGD) event and may have been important for tighter control of higher levels of RP transcripts. We find that subsequent loss of the derived repressor led to the loss of a stress-dependent repression of RPs in the fungal pathogen Candida glabrata. Our comparative computational and experimental approach shows how gene duplication can constrain and drive regulatory evolution and provides a general strategy for reconstructing the evolutionary trajectory of gene regulation across species.

Utility of Exome Sequencing for Diagnosis in Unexplained Pediatric-Onset Epilepsy.

Importance: Genomic advances inform our understanding of epilepsy and can be translated to patients as precision diagnoses that influence clinical treatment, prognosis, and counseling. Objective: To delineate the genetic landscape of pediatric epilepsy and clinical utility of genetic diagnoses for patients with epilepsy. Design, Setting, and Participants: This cohort study used phenotypic data from medical records and treating clinicians at a pediatric hospital to identify patients with unexplained pediatric-onset epilepsy. Exome sequencing was performed for 522 patients and available biological parents, and sequencing data were analyzed for single nucleotide variants (SNVs) and copy number variants (CNVs). Variant pathogenicity was assessed, patients were provided with their diagnostic results, and clinical utility was evaluated. Patients were enrolled from August 2018 to October 2021, and data were analyzed through December 2022. Exposures: Phenotypic features associated with diagnostic genetic results. Main Outcomes and Measures: Main outcomes included diagnostic yield and clinical utility. Diagnostic findings included variants curated as pathogenic, likely pathogenic (PLP), or diagnostic variants of uncertain significance (VUS) with clinical features consistent with the involved gene's associated phenotype. The proportion of the cohort with diagnostic findings, the genes involved, and their clinical utility, defined as impact on clinical treatment, prognosis, or surveillance, are reported. Results: A total of 522 children (269 [51.5%] male; mean [SD] age at seizure onset, 1.2 [1.4] years) were enrolled, including 142 children (27%) with developmental epileptic encephalopathy and 263 children (50.4%) with intellectual disability. Of these, 100 participants (19.2%) had identifiable genetic explanations for their seizures: 89 participants had SNVs (87 germline, 2 somatic mosaic) involving 69 genes, and 11 participants had CNVs. The likelihood of identifying a genetic diagnosis was highest in patients with intellectual disability (adjusted odds ratio [aOR], 2.44; 95% CI, 1.40-4.26), early onset seizures (aOR, 0.93; 95% CI, 0.88-0.98), and motor impairment (aOR, 2.19; 95% CI 1.34-3.58). Among 43 patients with apparently de novo variants, 2 were subsequently determined to have asymptomatic parents harboring mosaic variants. Of 71 patients who received diagnostic results and were followed clinically, 29 (41%) had documented clinical utility resulting from their genetic diagnoses. Conclusions and Relevance: These findings suggest that pediatric-onset epilepsy is genetically heterogeneous and that some patients with previously unexplained pediatric-onset epilepsy had genetic diagnoses with direct clinical implications.

Multicenter Consensus Approach to Evaluation of Neonatal Hypotonia in the Genomic Era: A Review.

IMPORTANCE: Infants with hypotonia can present with a variety of potentially severe clinical signs and symptoms and often require invasive testing and multiple procedures. The wide range of clinical presentations and potential etiologies leaves diagnosis and prognosis uncertain, underscoring the need for rapid elucidation of the underlying genetic cause of disease. OBSERVATIONS: The clinical application of exome sequencing or genome sequencing has dramatically improved the timely yield of diagnostic testing for neonatal hypotonia, with diagnostic rates of greater than 50% in academic neonatal intensive care units (NICUs) across Australia, Canada, the UK, and the US, which compose the International Precision Child Health Partnership (IPCHiP). A total of 74% (17 of 23) of patients had a change in clinical care in response to genetic diagnosis, including 2 patients who received targeted therapy. This narrative review discusses the common causes of neonatal hypotonia, the relative benefits and limitations of available testing modalities used in NICUs, and hypotonia management recommendations. CONCLUSIONS AND RELEVANCE: This narrative review summarizes the causes of neonatal hypotonia and the benefits of prompt genetic diagnosis, including improved prognostication and identification of targeted treatments which can improve the short-term and long-term outcomes. Institutional resources can vary among different NICUs; as a result, consideration should be given to rule out a small number of relatively unique conditions for which rapid targeted genetic testing is available. Nevertheless, the consensus recommendation is to use rapid genome or exome sequencing as a first-line testing option for NICU patients with unexplained hypotonia. As part of the IPCHiP, this diagnostic experience will be collected in a central database with the goal of advancing knowledge of neonatal hypotonia and improving evidence-based practice.

Refinements and considerations for trio whole-genome sequence analysis when investigating Mendelian diseases presenting in early childhood.

To facilitate early deployment of whole-genome sequencing (WGS) for severely ill children, a standardized pipeline for WGS analysis with timely turnaround and primary care pediatric uptake is needed. We developed a bioinformatics pipeline for comprehensive gene-agnostic trio WGS analysis of children suspected of having an undiagnosed monogenic disease that included detection and interpretation of primary genetic mechanisms of disease, including SNVs/indels, CNVs/SVs, uniparental disomy (UPD), imprinted genes, short tandem repeat expansions, mobile element insertions, SMN1/2 copy number calling, and mitochondrial genome variants. We assessed primary care practitioner experience and competence in a large cohort of 521 families (comprising 90% WGS trios). Children were identified by primary practitioners for recruitment, and we used the UK index of multiple deprivation to confirm lack of patient socio-economic status ascertainment bias. Of the 521 children sequenced, 176 (34%) received molecular diagnoses, with rates as high as 45% for neurology clinics. Twenty-three of the diagnosed cases (13%) required bespoke methods beyond routine SNV/CNV analysis. In our multidisciplinary clinician user experience assessment, both pediatricians and clinical geneticists expressed strong support for rapid WGS early in the care pathway, but requested further training in determining patient selection, consenting, and variant interpretation. Rapid trio WGS provides an efficacious single-pass screening test for children when deployed by primary practitioners in clinical settings that carry high a priori risk for rare pediatric disease presentations.

Heterozygous frameshift variants in HNRNPA2B1 cause early-onset oculopharyngeal muscular dystrophy.

Missense variants in RNA-binding proteins (RBPs) underlie a spectrum of disease phenotypes, including amyotrophic lateral sclerosis, frontotemporal dementia, and inclusion body myopathy. Here, we present ten independent families with a severe, progressive muscular dystrophy, reminiscent of oculopharyngeal muscular dystrophy (OPMD) but of much earlier onset, caused by heterozygous frameshift variants in the RBP hnRNPA2/B1. All disease-causing frameshift mutations abolish the native stop codon and extend the reading frame, creating novel transcripts that escape nonsense-mediated decay and are translated to produce hnRNPA2/B1 protein with the same neomorphic C-terminal sequence. In contrast to previously reported disease-causing missense variants in HNRNPA2B1, these frameshift variants do not increase the propensity of hnRNPA2 protein to fibrillize. Rather, the frameshift variants have reduced affinity for the nuclear import receptor karyopherin ß2, resulting in cytoplasmic accumulation of hnRNPA2 protein in cells and in animal models that recapitulate the human pathology. Thus, we expand the phenotypes associated with HNRNPA2B1 to include an early-onset form of OPMD caused by frameshift variants that alter its nucleocytoplasmic transport dynamics.

Dual-vector gene therapy restores cochlear amplification and auditory sensitivity in a mouse model of DFNB16 hearing loss.

Hearing loss affects an estimated 466 million people worldwide, with a substantial fraction due to genetic causes. Approximately 16% of genetic hearing loss is caused by pathogenic mutations in STRC, a gene that encodes the protein stereocilin. To develop gene therapy strategies for patients with STRC hearing loss, we generated a mouse model with a targeted deletion in the Strc gene. We devised a novel dual-vector approach to circumvent the size limitation of AAV vectors and drive expression of full-length STRC protein. To target outer hair cells, which are difficult to transduce, we used synthetic AAV9-PHP.B vectors for efficient dual-vector transduction. We report robust recovery of exogenous STRC expression in outer hair cells of Strc-deficient mice, recovery of hair bundle morphology, substantially improved cochlear amplification, and enhanced auditory sensitivity. The data raise the prospect that our strategy could benefit ~2.3 million patients worldwide affected by STRC mutations.

Abundancy of polymorphic CGG repeats in the human genome suggest a broad involvement in neurological disease.

Expanded CGG-repeats have been linked to neurodevelopmental and neurodegenerative disorders, including the fragile X syndrome and fragile X-associated tremor/ataxia syndrome (FXTAS). We hypothesized that as of yet uncharacterised CGG-repeat expansions within the genome contribute to human disease. To catalogue the CGG-repeats, 544 human whole genomes were analyzed. In total, 6101 unique CGG-repeats were detected of which more than 93% were highly variable in repeat length. Repeats with a median size of 12 repeat units or more were always polymorphic but shorter repeats were often polymorphic, suggesting a potential intergenerational instability of the CGG region even for repeats units with a median length of four or less. 410 of the CGG repeats were associated with known neurodevelopmental disease genes or with strong candidate genes. Based on their frequency and genomic location, CGG repeats may thus be a currently overlooked cause of human disease.

Risk Factors Associated with a Second Primary Lung Cancer in Patients with an Initial Primary Lung Cancer.

BACKGROUND: Increased patient survivorship following initial primary lung cancer (IPLC) diagnosis and treatment has uncovered new clinical challenges as individuals post-IPLC are at growing subsequent risk of developing second primary lung cancer (SPLC). Proper SPLC surveillance guidelines aimed at monitoring IPLC survivors are crucial to enhancing health outcomes. This study aims to categorize risk factors associated with SPLC emergence in IPLC survivors for clinical use following IPLC treatment. MATERIALS AND METHODS: Using the Karmanos Cancer Institute Tumor Registry, patients diagnosed with IPLC from 2000 to 2017 were identified. Patients diagnosed with SPLC were matched to individuals who did not develop SPLC. Logistic and Cox regression analyses were performed to identify risk factors for SPLC emergence and overall survival (OS). RESULTS: One hundred twenty-one patients diagnosed with IPLC who later developed SPLC were identified and compared with 120 patients with IPLC who did not develop SPLC. Several factors such as stage at first diagnosis, histology, age, and smoking history were not associated with SPLC risk. The median time to SPLC was 1.79 years. Patients who were treated with surgical resection had a significantly higher probability of developing SPLC. After correcting for potential immortal time bias, the median OS was 3.63 years (95% confidence interval [CI], 3.05-5.00) and 7.31 years (95% CI, 4.62-10.90) for SPLC and no SPLC groups, respectively. CONCLUSION: This study uncovered notable associations and lack thereof between several competing SPLC risk factors, as well as mortality. Further characterization of SPLC risk factors is essential for enhancing surveillance recommendations.