Search

Abstract Chordoma is a rare bone cancer with variable clinical outcomes. Here, we recruited 184 sporadic chordoma patients from the US and Canada and collected their clinical and treatment data. The average age at diagnosis was 45.5 years (Range 5–78) and the chordoma site distribution was 49.2% clivus, 26.2% spinal, and 24.0% sacral. Most patients (97.5%) received surgery as the primary treatment, among whom 85.3% also received additional treatment. Except for the most prevalent cancers like prostate, lung, breast, and skin cancer, there was no discernible enrichment for any specific cancer type among patients or their family members. Among a subset of patients (N = 70) with tumor materials, we conducted omics analyses and obtained targeted panel sequencing and SNP array genotyping data for 51 and 49 patients, respectively. The most recurrent somatic driver mutations included PIK3CA (12%), followed by chromatin remodeling genes PBRM1 and SETD2. Amplification of the 6q27 region, containing the chordoma susceptibility gene TBXT, was detected in eight patients (16.3%). Clival patients appeared to be less likely to carry driver gene mutations, chromosome arm level deletion events (e.g., 5p, 5p, and 9p), or 6q27 amplification compared to sacral patients. After adjusting for age, sex, tumor site, and additional treatment, patients with somatic deletions of 14q (OR = 13.73, 95% CI 1.96–96.02, P = 0.008) and 18p (OR = 13.68, 95% CI 1.77–105.89, P = 0.012) were more likely to have persistent chordoma. The study highlights genomic heterogeneity in chordoma, potentially linked to location and clinical progression.

Abstract Burkitt lymphoma (BL) is responsible for many childhood cancers in sub-Saharan Africa, where it is linked to recurrent or chronic infection by Epstein-Barr virus or Plasmodium falciparum. However, whether human leukocyte antigen (HLA) polymorphisms, which regulate immune response, are associated with BL has not been well investigated, which limits our understanding of BL etiology. Here we investigate this association among 4,645 children aged 0-15 years, 800 with BL, enrolled in Uganda, Tanzania, Kenya, and Malawi. HLA alleles are imputed with accuracy >90% for HLA class I and 85-89% for class II alleles. BL risk is elevated with HLA-DQA1*04:01 (adjusted odds ratio [OR] = 1.61, 95% confidence interval [CI] = 1.32-1.97, P = 3.71 × 10−6), with rs2040406(G) in HLA-DQA1 region (OR = 1.43, 95% CI = 1.26-1.63, P = 4.62 × 10−8), and with amino acid Gln at position 53 versus other variants in HLA-DQA1 (OR = 1.36, P = 2.06 × 10−6). The associations with HLA-DQA1*04:01 (OR = 1.29, P = 0.03) and rs2040406(G) (OR = 1.68, P = 0.019) persist in mutually adjusted models. The higher risk rs2040406(G) variant for BL is associated with decreased HLA-DQB1 expression in eQTLs in EBV transformed lymphocytes. Our results support the role of HLA variation in the etiology of BL and suggest that a promising area of research might be understanding the link between HLA variation and EBV control.

Abstract In high-income countries, mosaic chromosomal alterations in peripheral blood leukocytes are associated with an elevated risk of adverse health outcomes, including hematologic malignancies. We investigate mosaic chromosomal alterations in sub-Saharan Africa among 931 children with Burkitt lymphoma, an aggressive lymphoma commonly characterized by immunoglobulin-MYC chromosomal rearrangements, 3822 Burkitt lymphoma-free children, and 674 cancer-free men from Ghana. We find autosomal and X chromosome mosaic chromosomal alterations in 3.4% and 1.7% of Burkitt lymphoma-free children, and 8.4% and 3.7% of children with Burkitt lymphoma (P-values = 5.7×10−11 and 3.74×10−2, respectively). Autosomal mosaic chromosomal alterations are detected in 14.0% of Ghanaian men and increase with age. Mosaic chromosomal alterations in Burkitt lymphoma cases include gains on chromosomes 1q and 8, the latter spanning MYC, while mosaic chromosomal alterations in Burkitt lymphoma-free children include copy-neutral loss of heterozygosity on chromosomes 10, 14, and 16. Our results highlight mosaic chromosomal alterations in sub-Saharan African populations as a promising area of research.

Abstract Lung adenocarcinoma is the most common type of lung cancer. Known risk variants explain only a small fraction of lung adenocarcinoma heritability. Here, we conducted a two-stage genome-wide association study of lung adenocarcinoma of East Asian ancestry (21,658 cases and 150,676 controls; 54.5% never-smokers) and identified 12 novel susceptibility variants, bringing the total number to 28 at 25 independent loci. Transcriptome-wide association analyses together with colocalization studies using a Taiwanese lung expression quantitative trait loci dataset (n = 115) identified novel candidate genes, including FADS1 at 11q12 and ELF5 at 11p13. In a multi-ancestry meta-analysis of East Asian and European studies, four loci were identified at 2p11, 4q32, 16q23, and 18q12. At the same time, most of our findings in East Asian populations showed no evidence of association in European populations. In our studies drawn from East Asian populations, a polygenic risk score based on the 25 loci had a stronger association in never-smokers vs. individuals with a history of smoking (Pinteraction = 0.0058). These findings provide new insights into the etiology of lung adenocarcinoma in individuals from East Asian populations, which could be important in developing translational applications.

Abstract Background Few studies have examined epigenetic age acceleration (AA), the difference between DNA methylation (DNAm) predicted age and chronological age, in relation to somatic genomic features in paired cancer and normal tissue, with less work done in non-European populations. In this study, we aimed to examine DNAm age and its associations with breast cancer risk factors, subtypes, somatic genomic profiles including mutation and copy number alterations and other aging markers in breast tissue of Chinese breast cancer (BC) patients from Hong Kong. Methods We performed genome-wide DNA methylation profiling of 196 tumor and 188 paired adjacent normal tissue collected from Chinese BC patients in Hong Kong (HKBC) using Illumina MethylationEPIC array. The DNAm age was calculated using Horvath’s pan-tissue clock model. Somatic genomic features were based on data from RNA sequencing (RNASeq), whole-exome sequencing (WES), and whole-genome sequencing (WGS). Pearson’s correlation (r), Kruskal–Wallis test, and regression models were used to estimate associations of DNAm AA with somatic features and breast cancer risk factors. Results DNAm age showed a stronger correlation with chronological age in normal (Pearson r = 0.78, P

Abstract The Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial is a prospective cohort study of nearly 155,000 U.S. volunteers aged 55–74 at enrollment in 1993–2001. We developed the PLCO Atlas Project, a large resource for multi-trait genome-wide association studies (GWAS), by genotyping participants with available DNA and genomic consent. Genotyping on high-density arrays and imputation was performed, and GWAS were conducted using a custom semi-automated pipeline. Association summary statistics were generated from a total of 110,562 participants of European, African and Asian ancestry. Application programming interfaces (APIs) and open-source software development kits (SKDs) enable exploring, visualizing and open data access through the PLCO Atlas GWAS Explorer website, promoting Findable, Accessible, Interoperable, and Re-usable (FAIR) principles. Currently the GWAS Explorer hosts association data for 90 traits and >78,000,000 genomic markers, focusing on cancer and cancer-related phenotypes. New traits will be posted as association data becomes available. The PLCO Atlas is a FAIR resource of high-quality genetic and phenotypic data with many potential reuse opportunities for cancer research and genetic epidemiology.

Myelofibrosis is a risk factor for the development of Acute Myeloid Leukaemia. Here, the authors carry out an integrated genomic investigation of 933 myelofibrosis patients, and identified interactions between germline and somatic variation in patients who required haematopoietic cell transplantation.

Age-related DNA methylation is a potential mechanism contributing to breast cancer development. Studies of primarily Caucasian women have identified many CpG sites of age-related methylation in non-diseased breast tissue possibly driving cancer development over time. There is a paucity of studies involving Asian women whose ages at breast cancer onset are usually younger than Caucasians. We identified the 181 most consistent age-related methylation events in non-diseased breast tissue across published studies. Age-related methylation events were measured in adjacent normal and breast tumour tissue in an exclusively Asian population at the previously identified age-related methylation sites. Age-related methylation was found in 118 probes in adjacent normal breast tissue. Methylation of 99% of these sites was increased with age and predominantly located on CpG islands in promoter regions. To ascertain biological relevance to breast cancer, we focused on the 37 sites with overall higher methylation in tumour compared to adjacent normal samples. Some sites positively related to age, including AQP5 and CORO6, inversely correlated with gene expression. Several others have known involvement in suppression of carcinogenesis including GPC5 and SST, suggesting that perturbation of epigenetic regulation at these sites due to ageing may contribute to the progression of carcinogenesis. This study highlights an age-related methylation landscape in non-tumour tissue, consistent not just across studies, but also across different populations. We present candidate age-related methylation sites warranting further investigation as potential epigenetic drivers of breast cancer. They may serve as potential targets of site-specific demethylation intervention strategies for the prevention of age-related breast cancer.

Skull base chordomas are treated with surgery and chemotherapy but often recur due to incomplete resection, understanding the molecular underpinnings of the tumours may provide additional therapeutic strategies. Here, the authors carry out whole genome sequencing of 80 skull base chordoma tumours and identify the SWI/SNF component—PBRM1—as a frequently mutated gene.

Many tumours are heterogenous, which calls into question whether multiple biopsies are required to accurately assess the cancer. Here, the authors show that papillary renal cell carcinoma is clonal in nature, suggesting that in this cancer one biopsy is sufficient for diagnosis and molecular analyses.

Many tumors are known to be heterogeneous. Here, the authors examined multiple samples from 84 patients with lung adenocarcinoma and demonstrate that the intratumor heterogeneity of methylation and copy number associates with poor prognosis.

Abstract Myotonic dystrophy type I (DM1) is an autosomal dominant disease of which clinical manifestations resemble premature aging. We evaluated the contribution of telomere length in pathogenesis in 361 DM1 patients (12 with serial measurements) and 223 unaffected relative controls using qPCR assay. While no differences in baseline leukocyte relative telomere length (RTL) was noted, the data suggested an accelerated RTL attrition in DM1 (discovery cohort: T/S change/year = −0.013 in DM1 vs. −0.005 in controls, P = 0.04); similar trend was noted in validation cohort. Further investigations are needed to examine the role of TL in the pathophysiology of DM1.

Abstract Background Oral microbiota may be related to pancreatic cancer risk because periodontal disease, a condition linked to multiple specific microbes, has been associated with increased risk of pancreatic cancer. We evaluated the association between oral microbiota and pancreatic cancer in Iran. Methods A total of 273 pancreatic adenocarcinoma cases and 285 controls recruited from tertiary hospitals and a specialty clinic in Tehran, Iran provided saliva samples and filled out a questionnaire regarding demographics and lifestyle characteristics. DNA was extracted from saliva and the V4 region of the 16S rRNA gene was PCR amplified and sequenced on the MiSeq. The sequencing data were processed using the DADA2 plugin in QIIME 2 and taxonomy was assigned against the Human Oral Microbiome Database. Logistic regression and MiRKAT models were calculated with adjustment for potential confounders. Results No association was observed for alpha diversity with an average of 91.11 (standard deviation [SD] 2.59) sequence variants for cases and 89.42 (SD 2.58) for controls. However, there was evidence for an association between beta diversity and case status. The association between the Bray‐Curtis dissimilarity and pancreatic cancer was particularly strong with a MiRKAT P‐value of .000142 and specific principal coordinate vectors had strong associations with cancer risk. Several specific taxa were also associated with case status after adjustment for multiple comparisons. Conclusion The overall microbial community appeared to differ between pancreatic cancer cases and controls. Whether these reflect differences evident before development of pancreatic cancer will need to be evaluated in prospective studies.

Summary: Recent genomic studies suggest that Asian breast cancer (BC) may have distinct somatic features; however, most comparisons of BC genomic features across populations did not account for differences in age, subtype, and sequencing methods. In this study, we analyzed whole-exome sequencing (WES) data to characterize somatic copy number alterations (SCNAs) and mutation profiles in 98 Hong Kong BC (HKBC) patients and compared with those from The Cancer Genome Atlas of European ancestry (TCGA-EA, N = 686), which had similar distributions of age at diagnosis and PAM50 subtypes as in HKBC. We developed a two-sample Poisson model to compare driver gene selection pressure, which reflects the effect sizes of cancer driver genes, while accounting for differences in sample size, sequencing platforms, depths, and mutation calling methods. We found that somatic mutation and SCNA profiles were overall very similar between HKBC and TCGA-EA. The selection pressure for small insertions and deletions (indels) in GATA3 (false discovery rate (FDR) corrected p < 0.01) and single-nucleotide variants (SNVs) in TP53 (nominal p = 0.02, FDR corrected p = 0.28) was lower in HKBC than in TCGA-EA. Among the 13 signatures of single-base substitutions (SBS) that are common in BC, we found a suggestively higher contribution of SBS18 and a lower contribution of SBS1 in HKBC than in TCGA-EA, while the two APOBEC-induced signatures showed similar prevalence. Our results suggest that the genomic landscape of BC was largely very similar between HKBC and TCGA-EA, despite suggestive differences in some driver genes and mutational signatures that warrant future investigations in large and diverse Asian populations.

Abstract Background Heterogeneity of immune gene expression patterns of luminal breast cancer (BC), which is clinically heterogeneous and overall considered as low immunogenic, has not been well studied especially in non-European populations. Here, we aimed at characterizing the immune gene expression profile of luminal BC in an Asian population and associating it with patient characteristics and tumor genomic features. Methods We performed immune gene expression profiling of tumor and adjacent normal tissue in 92 luminal BC patients from Hong Kong using RNA-sequencing data and used unsupervised consensus clustering to stratify tumors. We then used luminal patients from The Cancer Genome Atlas (TCGA, N = 564) and a Korean breast cancer study (KBC, N = 112) as replication datasets. Results Based on the expression of 130 immune-related genes, luminal tumors were stratified into three distinct immune subtypes. Tumors in one subtype showed higher level of tumor-infiltrating lymphocytes (TILs), characterized by T cell gene activation, higher expression of immune checkpoint genes, higher nonsynonymous mutation burden, and higher APOBEC-signature mutations, compared with other luminal tumors. The high-TIL subtype was also associated with lower ESR1/ESR2 expression ratio and increasing body mass index. The comparison of the immune profile in tumor and matched normal tissue suggested a tumor-derived activation of specific immune responses, which was only seen in high-TIL patients. Tumors in a second subtype were characterized by increased expression of interferon-stimulated genes and enrichment for TP53 somatic mutations. The presence of three immune subtypes within luminal BC was replicated in TCGA and KBC, although the pattern was more similar in Asian populations. The germline APOBEC3B deletion polymorphism, which is prevalent in East Asian populations and was previously linked to immune activation, was not associated with immune subtypes in our study. This result does not support the hypothesis that the germline APOBEC3B deletion polymorphism is the driving force for immune activation in breast tumors in Asian populations. Conclusion Our findings suggest that immune gene expression and associated genomic features could be useful to further stratify luminal BC beyond the current luminal A/B classification and a subset of luminal BC patients may benefit from checkpoint immunotherapy, at least in Asian populations.

ABSTRACT Cigarettes and opium contain chemicals and particulate matter that may modify the oral microbiota. This study aimed to investigate the association between cigarette and opium use with the oral microbiota. A total of 558 participants were recruited from Iran between 2011 and 2015. Individuals were categorized as never cigarette nor opium users, ever cigarette-only smokers, ever opium-only users, and ever both cigarette and opium users. Participants provided saliva samples for 16S rRNA gene sequencing. Logistic regression, microbiome regression-based kernel association test (MiRKAT), and zero-inflated beta regression models were calculated. For every increase in 10 observed amplicon sequence variants (ASVs), the odds for being a cigarette-only smoker, opium-only user, and both user compared to never users decreased by 9% (odds ratio [OR] = 0.91; 95% confidence interval [95% CI] = 0.86 to 0.97), 13% (OR = 0.87; 95% CI = 0.75 to 1.01), and 12% (OR = 0.88; 95% CI = 0.80 to 0.96), respectively. The microbial communities differed by cigarette and opium use as indicated by MiRKAT models testing the three beta-diversity matrices (P < 0.05 for all). Three genera were less likely and one genus was more likely to be detected in cigarette-only smokers or opium-only users than in never users. The relative abundance of the phylum Actinobacteria (never, 14.78%; both, 21.20%) was higher and the phyla Bacteroidetes (never, 17.63%; both, 11.62%) and Proteobacteria (never, 9.06%; both, 3.70%) were lower in users of both cigarettes and opium, while the phylum Firmicutes (never, 54.29%; opium, 65.49%) was higher in opium-only users. Cigarette and opium use was associated with lower alpha-diversity, overall oral microbiota community composition, and both the presence and relative abundance of multiple taxa. IMPORTANCE Cigarette smoking and opium use are associated with periodontal disease caused by specific bacteria such as Porphyromonas gingivalis, which suggests a link between cigarette smoking and opium use and the oral microbiota. Alterations of the oral microbiota in cigarette smokers compared to nonsmokers have been reported, but this has not been studied across diverse populations. Additionally, the association of opium use with the oral microbiota has not been investigated to date. We conducted this study to investigate differences in the oral microbiota between ever users of cigarettes only, opium only, and both cigarettes and opium and never users of cigarettes and opium in Iran. Lower alpha-diversity, distinct overall oral microbial communities, and the presence and relative abundance of multiple taxa have been found for users of cigarettes and/or opium.

Abstract Background Environmental and genetic factors play an important role in the etiology of breast cancer. Several small blood-based DNA methylation studies have reported risk associations with methylation at individual CpGs and average methylation levels; however, these findings require validation in larger prospective cohort studies. To investigate the role of blood DNA methylation on breast cancer risk, we conducted a meta-analysis of four prospective cohort studies, including a total of 1663 incident cases and 1885 controls, the largest study of blood DNA methylation and breast cancer risk to date. Methods We assessed associations with methylation at 365,145 CpGs present in the HumanMethylation450 (HM450K) Beadchip, after excluding CpGs that did not pass quality controls in all studies. Each of the four cohorts estimated odds ratios (ORs) and 95% confidence intervals (CI) for the association between each individual CpG and breast cancer risk. In addition, each study assessed the association between average methylation measures and breast cancer risk, adjusted and unadjusted for cell-type composition. Study-specific ORs were combined using fixed-effect meta-analysis with inverse variance weights. Stratified analyses were conducted by age at diagnosis ( 10 years). The false discovery rate (q value) was used to account for multiple testing. Results The average age at blood draw ranged from 52.2 to 62.2 years across the four cohorts. Median follow-up time ranged from 6.6 to 8.4 years. The methylation measured at individual CpGs was not associated with breast cancer risk (q value > 0.59). In addition, higher average methylation level was not associated with risk of breast cancer (OR = 0.94, 95% CI = 0.85, 1.05; P = 0.26; P for study heterogeneity = 0.86). We found no evidence of modification of this association by age at diagnosis (P = 0.17), ER status (P = 0.88), time since blood collection (P = 0.98), or CpG location (P = 0.98). Conclusions Our data indicate that DNA methylation measured in the blood prior to breast cancer diagnosis in predominantly postmenopausal women is unlikely to be associated with substantial breast cancer risk on the HM450K array. Larger studies or with greater methylation coverage are needed to determine if associations exist between blood DNA methylation and breast cancer risk.

Abstract Background Prior research has established that the prevalence of pathogenic/likely pathogenic (P/LP) variants across all of the American College of Medical Genetics (ACMG) Secondary Findings (SF) genes is approximately 0.8–5%. We investigated the prevalence of P/LP variants in the 24 ACMG SF v2.0 cancer genes in a family-based cancer research cohort (n = 1173) and in cancer-free ethnicity-matched controls (n = 982). Methods We used InterVar to classify variants and subsequently conducted a manual review to further examine variants of unknown significance (VUS). Results In the 24 genes on the ACMG SF v2.0 list associated with a cancer phenotype, we observed 8 P/LP unique variants (8 individuals; 0.8%) in controls and 11 P/LP unique variants (14 individuals; 1.2%) in cases, a non-significant difference. We reviewed 115 VUS. The median estimated per-variant review time required was 30 min; the first variant within a gene took significantly (p = 0.0009) longer to review (median = 60 min) compared with subsequent variants (median = 30 min). The concordance rate was 83.3% for the variants examined by two reviewers. Conclusion The 115 VUS required database and literature review, a time- and labor-intensive process hampered by the difficulty in interpreting conflicting P/LP determinations. By rigorously investigating the 24 ACMG SF v2.0 cancer genes, our work establishes a benchmark P/LP variant prevalence rate in a familial cancer cohort and controls.

Abstract Background High-throughput sequencing has rapidly become an essential part of precision cancer medicine. But validating results obtained from analyzing and interpreting genomic data remains a rate-limiting factor. The gold standard, of course, remains manual validation by expert panels, which is not without its weaknesses, namely high costs in both funding and time as well as the necessarily selective nature of manual validation. But it may be possible to develop more economical, complementary means of validation. In this study we employed four synthetic data sets (variants with known mutations spiked into specific genomic locations) of increasing complexity to assess the sensitivity, specificity, and balanced accuracy of five open-source variant callers: FreeBayes v1.0, VarDict v11.5.1, MuTect v1.1.7, MuTect2, and MuSE v1.0rc. FreeBayes, VarDict, and MuTect were run in bcbio-next gen, and the results were integrated into a single Ensemble call set. The known mutations provided a level of “ground truth” against which we evaluated variant-caller performance. We further facilitated the comparison and evaluation by segmenting the whole genome into 10,000,000 base-pair fragments which yielded 316 segments. Results Differences among the numbers of true positives were small among the callers, but the numbers of false positives varied much more when the tools were used to analyze sets one through three. Both FreeBayes and VarDict produced strikingly more false positives than did the others, although VarDict, somewhat paradoxically also produced the highest number of true positives. The Ensemble approach yielded results characterized by higher specificity and balanced accuracy and fewer false positives than did any of the five tools used alone. Sensitivity and specificity, however, declined for all five callers as the complexity of the data sets increased, but we did not uncover anything more than limited, weak correlations between caller performance and certain DNA structural features: gene density and guanine-cytosine content. Altogether, MuTect2 performed the best among the callers tested, followed by MuSE and MuTect. Conclusions Spiking data sets with specific mutations –single-nucleotide variations (SNVs), single-nucleotide polymorphisms (SNPs), or structural variations (SVs) in this study—at known locations in the genome provides an effective and economical way to compare data analyzed by variant callers with ground truth. The method constitutes a viable alternative to the prolonged, expensive, and noncomprehensive assessment by expert panels. It should be further developed and refined, as should other comparatively “lightweight” methods of assessing accuracy. Given that the scientific community has not yet established gold standards for validating NGS-related technologies such as variant callers, developing multiple alternative means for verifying variant-caller accuracy will eventually lead to the establishment of higher-quality standards than could be achieved by prematurely limiting the range of innovative methods explored by members of the community.

It is unclear how often genetic mosaicism of chromosome X arises. Here, the authors examine women with cancer and cancer-free controls and show that X chromosome mosaicism occurs more frequently than on autosomes, especially on the inactive X chromosome, but is not linked to non-haematologic cancer risk

Telomeres, long nucleotide repeats and a protein complex at chromosome ends, shorten with each cell division and are susceptible to oxidative damage. Quantitative PCR (qPCR) is a widely-used technique to measure relative telomere length (RTL) in DNA samples but is challenging to optimize and significant lab-to-lab variability has been reported. In this study, we evaluated factors that may contribute to qPCR RTL measurement variability including DNA extraction methods, methods used for removing potential residual PCR inhibitors, sample storage conditions, and sample location in the PCR plate. Our results show that the DNA extraction and purification techniques, as well as sample storage conditions introduce significant variability in qPCR RTL results. We did not find significant differences in results based on sample location in the PCR plate or qPCR instrument used. These data suggest that lack of reproducibility in published association studies of RTL could be, in part, due to methodological inconsistencies. This study illustrates the importance of uniform sample handling, from DNA extraction through data generation and analysis, in using qPCR to determine RTL.

BACKGROUND:Lung adenocarcinoma (LUAD) is the most common histologic subtype of lung cancer and has a high risk of distant metastasis at every disease stage. We aimed to characterize the genomic landscape of LUAD and identify mutation signatures associated with tumor progression. METHODS AND FINDINGS:We performed an integrative genomic analysis, incorporating whole exome sequencing (WES), determination of DNA copy number and DNA methylation, and transcriptome sequencing for 101 LUAD samples from the Environment And Genetics in Lung cancer Etiology (EAGLE) study. We detected driver genes by testing whether the nonsynonymous mutation rate was significantly higher than the background mutation rate and replicated our findings in public datasets with 724 samples. We performed subclonality analysis for mutations based on mutant allele data and copy number alteration data. We also tested the association between mutation signatures and clinical outcomes, including distant metastasis, survival, and tumor grade. We identified and replicated two novel candidate driver genes, POU class 4 homeobox 2 (POU4F2) (mutated in 9 [8.9%] samples) and ZKSCAN1 (mutated in 6 [5.9%] samples), and characterized their major deleterious mutations. ZKSCAN1 was part of a mutually exclusive gene set that included the RTK/RAS/RAF pathway genes BRAF, EGFR, KRAS, MET, and NF1, indicating an important driver role for this gene. Moreover, we observed strong associations between methylation in specific genomic regions and somatic mutation patterns. In the tumor evolution analysis, four driver genes had a significantly lower fraction of subclonal mutations (FSM), including TP53 (p = 0.007), KEAP1 (p = 0.012), STK11 (p = 0.0076), and EGFR (p = 0.0078), suggesting a tumor initiation role for these genes. Subclonal mutations were significantly enriched in APOBEC-related signatures (p < 2.5×10-50). The total number of somatic mutations (p = 0.0039) and the fraction of transitions (p = 5.5×10-4) were associated with increased risk of distant metastasis. Our study's limitations include a small number of LUAD patients for subgroup analyses and a single-sample design for investigation of subclonality. CONCLUSIONS:These data provide a genomic characterization of LUAD pathogenesis and progression. The distinct clonal and subclonal mutation signatures suggest possible diverse carcinogenesis pathways for endogenous and exogenous exposures, and may serve as a foundation for more effective treatments for this lethal disease. LUAD's high heterogeneity emphasizes the need to further study this tumor type and to associate genomic findings with clinical outcomes.

Hodgkin lymphoma shows strong familial aggregation but no major susceptibility genes have been identified to date. The goal of this study was to identify high-penetrance variants using whole exome sequencing in 17 Hodgkin lymphoma prone families with three or more affected cases or obligate carriers (69 individuals), followed by targeted sequencing in an additional 48 smaller HL families (80 individuals). Alignment and variant calling were performed using standard methods. Dominantly segregating, rare, coding or potentially functional variants were further prioritized based on predicted deleteriousness, conservation, and potential importance in lymphoid malignancy pathways. We selected 23 genes for targeted sequencing. Only the p.A1065T variant in KDR (kinase insert domain receptor) also known as VEGFR2 (vascular endothelial growth factor receptor 2) was replicated in two independent Hodgkin lymphoma families. KDR is a type III receptor tyrosine kinase, the main mediator of vascular endothelial growth factor induced proliferation, survival, and migration. Its activity is associated with several diseases including lymphoma. Functional experiments have shown that p.A1065T, located in the activation loop, can promote constitutive autophosphorylation on tyrosine in the absence of vascular endothelial growth factor and that the kinase activity was abrogated after exposure to kinase inhibitors. A few other promising mutations were identified but appear to be “private”. In conclusion, in the largest sequenced cohort of Hodgkin lymphoma families to date, we identified a causal mutation in the KDR gene. While independent validation is needed, this mutation may increase downstream tumor cell proliferation activity and might be a candidate for targeted therapy.

Telomere length (TL) comparisons from different methods are challenging due to differences in laboratory techniques and data configuration. This study aimed to assess the validity of converting the quantitative polymerase chain reaction (qPCR) telomere/single copy gene (T/S) ratio to TL in kilobases (kb). We developed a linear regression equation to predict TL from qPCR T/S using flow cytometry with fluorescence in situ hybridization (flow FISH) TL data from 181 healthy donors (age range = 19⁻53) from the National Marrow Donor Program (NMDP) biorepository. TL measurements by qPCR and flow FISH were modestly correlated (R2 = 0.56, p < 0.0001). In Bland-Altman analyses, individuals with the shortest (≤10th percentile) or longest (≥90th) flow FISH TL had an over- or under-estimated qPCR TL (bias = 0.89 and −0.77 kb, respectively). Comparisons of calculated TL from the NMDP samples and 1810 age- and sex-matched individuals from the National Health and Nutrition Examination Survey showed significant differences (median = 7.1 versus 5.8 kb, respectively, p < 0.0001). Differences in annual TL attrition were also noted (31 versus 13 bp/year, respectively, p = 0.02). Our results demonstrate that TL calculated in kb from qPCR T/S may yield biased estimates for individuals with the shortest or longest TL, those often of high clinical interest. We also showed that calculated TL in kb from qPCR data are not comparable across populations and therefore are not necessarily useful.