Analyses of non-coding somatic drivers in 2,658 cancer whole genomes.

The discovery of drivers of cancer has traditionally focused on protein-coding genes1-4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5' region of TP53, in the 3' untranslated regions of NFKBIZ and TOB1, focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

Genomic features of renal cell carcinoma developed during end-stage renal disease and dialysis.

Patients with end-stage renal disease (ESRD) or receiving dialysis have a much higher risk for renal cell carcinoma (RCC), but carcinogenic mechanisms and genomic features remain little explored and undefined. This study's goal was to identify the genomic features of ESRD RCC and characterize them for associations with tumor histology and dialysis exposure. In this study, we obtained 33 RCCs, with various histological subtypes, that developed in ESRD patients receiving dialysis and performed whole-genome sequencing and transcriptome analyses. Driver events, copy-number alteration (CNA) analysis and mutational signature profiling were performed using an analysis pipeline that integrated data from germline and somatic SNVs, Indels and structural variants as well as CNAs, while transcriptome data were analyzed for differentially expressed genes and through gene set enrichment analysis. ESRD related clear cell RCCs' driver genes and mutations mirrored those in sporadic ccRCCs. Longer dialysis periods significantly correlated with a rare mutational signature SBS23, whose etiology is unknown, and increased mitochondrial copy number. All acquired cystic disease (ACD)-RCCs, which developed specifically in ESRD patients, showed chromosome 16q amplification. Gene expression analysis suggests similarity between certain ACD-RCCs and papillary RCCs and in TCGA papillary RCCs with chromosome 16 gain identified enrichment for genes related to DNA repair, as well as pathways related to reactive oxygen species, oxidative phosphorylation and targets of Myc. This analysis suggests that ESRD or dialysis could induce types of cellular stress that impact some specific types of genomic damage leading to oncogenesis.

Variant spectrum of von Hippel-Lindau disease and its genomic heterogeneity in Japan.

Von Hippel-Lindau (VHL) disease is an autosomal dominant, inherited syndrome with variants in the VHL gene, causing predisposition to multi-organ neoplasms with vessel abnormality. Germline variants in VHL can be detected in 80-90% of patients clinically diagnosed with VHL disease. Here, we summarize the results of genetic tests for 206 Japanese VHL families, and elucidate the molecular mechanisms of VHL disease, especially in variant-negative unsolved cases. Of the 206 families, genetic diagnosis was positive in 175 families (85%), including 134 families (65%) diagnosed by exon sequencing (15 novel variants) and 41 (20%) diagnosed by multiplex ligation-dependent probe amplification (MLPA) (one novel variant). The deleterious variants were significantly enriched in VHL disease Type 1. Interestingly, five synonymous or non-synonymous variants within exon 2 caused exon 2 skipping, which is the first report of exon 2 skipping caused by several missense variants. Whole genome and target deep sequencing analysis were performed for 22 unsolved cases with no variant identified and found three cases with VHL mosaicism (variant allele frequency: 2.5-22%), one with mobile element insertion in the VHL promoter region, and two with a pathogenic variant of BAP1 or SDHB. The variants associated with VHL disease are heterogeneous, and for more accuracy of the genetic diagnosis of VHL disease, comprehensive genome and DNA/RNA analyses are required to detect VHL mosaicism, complicated structure variants and other related gene variants.

Hereditary cancer variants and homologous recombination deficiency in biliary tract cancer.

BACKGROUND & AIMS: The heritability and actionability of variants in homologous recombination-related genes in biliary tract cancers (BTCs) are uncertain. Although associations between BTC and BRCA germline variants have been reported, homologous recombination deficiency has not been investigated in BTCs. METHODS: We sequenced germline variants in 27 cancer-predisposing genes in 1,292 BTC cases and 37,583 controls without a personal nor family history of cancer. We compared pathogenic germline variant frequencies between cases and controls and documented the demographic and clinical characteristics of carriers. In addition, whole-genome sequencing of 45 BTC tissues was performed to evaluate homologous recombination deficiency status. RESULTS: Targeted sequencing identified 5,018 germline variants, which were classified into 317 pathogenic, 3,611 variants of uncertain significance, and 1,090 benign variants. Seventy-one BTC cases (5.5%) had at least one pathogenic variant among 27 cancer-predisposing genes. Pathogenic germline variants enriched in BTCs were present in BRCA1, BRCA2, APC, and MSH6 (p <0.00185). PALB2 variants were marginally associated with BTC (p = 0.01). APC variants were predominantly found in ampulla of Vater carcinomas. Whole-genome sequencing demonstrated that three BTCs with pathogenic germline variants in BRCA2 and PALB2, accompanied by loss of heterozygosity, displayed homologous recombination deficiency. Conversely, pathogenic germline variants without a second hit or variants of other homologous recombination-related genes such as ATM and BRIP1 showed homologous recombination-proficient phenotypes. CONCLUSIONS: In this study, we describe the heritability and actionability of variants in homologous recombination-related genes, which could be used to guide screening and therapeutic strategies for BTCs. IMPACT AND IMPLICATIONS: We found that 5.5% of biliary tract cancers (BTCs) in a Japanese population possessed hereditary cancer-predisposing gene alterations, including in BRCA and genes associated with colorectal cancer. Two hits in homologous recombination-related genes were required to confer a homologous recombination-deficient phenotype. PARP inhibitors and DNA-damaging regimens may be effective strategies against BTCs exhibiting homologous recombination deficiency. Hence, in this study, genome-wide sequencing has revealed a potential new therapeutic strategy that could be applied to a subset of BTCs.

Mitochondrial DNA Variants at Low-Level Heteroplasmy and Decreased Copy Numbers in Chronic Kidney Disease (CKD) Tissues with Kidney Cancer.

The human mitochondrial genome (mtDNA) is a circular DNA molecule with a length of 16.6 kb, which contains a total of 37 genes. Somatic mtDNA mutations accumulate with age and environmental exposure, and some types of mtDNA variants may play a role in carcinogenesis. Recent studies observed mtDNA variants not only in kidney tumors but also in adjacent kidney tissues, and mtDNA dysfunction results in kidney injury, including chronic kidney disease (CKD). To investigate whether a relationship exists between heteroplasmic mtDNA variants and kidney function, we performed ultra-deep sequencing (30,000×) based on long-range PCR of DNA from 77 non-tumor kidney tissues of kidney cancer patients with CKD (stages G1 to G5). In total, this analysis detected 697 single-nucleotide variants (SNVs) and 504 indels as heteroplasmic (0.5% ≤ variant allele frequency (VAF) < 95%), and the total number of detected SNVs/indels did not differ between CKD stages. However, the number of deleterious low-level heteroplasmic variants (pathogenic missense, nonsense, frameshift and tRNA) significantly increased with CKD progression (p < 0.01). In addition, mtDNA copy numbers (mtDNA-CNs) decreased with CKD progression (p < 0.001). This study demonstrates that mtDNA damage, which affects mitochondrial genes, may be involved in reductions in mitochondrial mass and associated with CKD progression and kidney dysfunction.

Association of an IGHV3-66 gene variant with Kawasaki disease.

In a meta-analysis of three GWAS for susceptibility to Kawasaki disease (KD) conducted in Japan, Korea, and Taiwan and follow-up studies with a total of 11,265 subjects (3428 cases and 7837 controls), a significantly associated SNV in the immunoglobulin heavy variable gene (IGHV) cluster in 14q33.32 was identified (rs4774175; OR = 1.20, P = 6.0 × 10-9). Investigation of nonsynonymous SNVs of the IGHV cluster in 9335 Japanese subjects identified the C allele of rs6423677, located in IGHV3-66, as the most significant reproducible association (OR = 1.25, P = 6.8 × 10-10 in 3603 cases and 5731 controls). We observed highly skewed allelic usage of IGHV3-66, wherein the rs6423677 A allele was nearly abolished in the transcripts in peripheral blood mononuclear cells of both KD patients and healthy adults. Association of the high-expression allele with KD strongly indicates some active roles of B-cells or endogenous immunoglobulins in the disease pathogenesis. Considering that significant association of SNVs in the IGHV region with disease susceptibility was previously known only for rheumatic heart disease (RHD), a complication of acute rheumatic fever (ARF), these observations suggest that common B-cell related mechanisms may mediate the symptomology of KD and ARF as well as RHD.

Genome-wide meta-analysis and replication studies in multiple ethnicities identify novel adolescent idiopathic scoliosis susceptibility loci.

Adolescent idiopathic scoliosis (AIS) is the most common musculoskeletal disorder of childhood development. The genetic architecture of AIS is complex, and the great majority of risk factors are undiscovered. To identify new AIS susceptibility loci, we conducted the first genome-wide meta-analysis of AIS genome-wide association studies, including 7956 cases and 88 459 controls from 3 ancestral groups. Three novel loci that surpassed genome-wide significance were uncovered in intragenic regions of the CDH13 (P-value_rs4513093 = 1.7E-15), ABO (P-value_ rs687621 = 7.3E-10) and SOX6 (P-value_rs1455114 = 2.98E-08) genes. Restricting the analysis to females improved the associations at multiple loci, most notably with variants within CDH13 despite the reduction in sample size. Genome-wide gene-functional enrichment analysis identified significant perturbation of pathways involving cartilage and connective tissue development. Expression of both SOX6 and CDH13 was detected in cartilage chondrocytes and chromatin immunoprecipitation sequencing experiments in that tissue revealed multiple HeK27ac-positive peaks overlapping associated loci. Our results further define the genetic architecture of AIS and highlight the importance of vertebral cartilage development in its pathogenesis.

Association analyses of East Asian individuals and trans-ancestry analyses with European individuals reveal new loci associated with cholesterol and triglyceride levels.

Large-scale meta-analyses of genome-wide association studies (GWAS) have identified >175 loci associated with fasting cholesterol levels, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). With differences in linkage disequilibrium (LD) structure and allele frequencies between ancestry groups, studies in additional large samples may detect new associations. We conducted staged GWAS meta-analyses in up to 69,414 East Asian individuals from 24 studies with participants from Japan, the Philippines, Korea, China, Singapore, and Taiwan. These meta-analyses identified (P < 5 × 10-8) three novel loci associated with HDL-C near CD163-APOBEC1 (P = 7.4 × 10-9), NCOA2 (P = 1.6 × 10-8), and NID2-PTGDR (P = 4.2 × 10-8), and one novel locus associated with TG near WDR11-FGFR2 (P = 2.7 × 10-10). Conditional analyses identified a second signal near CD163-APOBEC1. We then combined results from the East Asian meta-analysis with association results from up to 187,365 European individuals from the Global Lipids Genetics Consortium in a trans-ancestry meta-analysis. This analysis identified (log10Bayes Factor ≥6.1) eight additional novel lipid loci. Among the twelve total loci identified, the index variants at eight loci have demonstrated at least nominal significance with other metabolic traits in prior studies, and two loci exhibited coincident eQTLs (P < 1 × 10-5) in subcutaneous adipose tissue for BPTF and PDGFC. Taken together, these analyses identified multiple novel lipid loci, providing new potential therapeutic targets.

Combined Genetic and Genealogic Studies Uncover a Large BAP1 Cancer Syndrome Kindred Tracing Back Nine Generations to a Common Ancestor from the 1700s.

We recently discovered an inherited cancer syndrome caused by BRCA1-Associated Protein 1 (BAP1) germline mutations, with high incidence of mesothelioma, uveal melanoma and other cancers and very high penetrance by age 55. To identify families with the BAP1 cancer syndrome, we screened patients with family histories of multiple mesotheliomas and melanomas and/or multiple cancers. We identified four families that shared an identical BAP1 mutation: they lived across the US and did not appear to be related. By combining family histories, molecular genetics, and genealogical approaches, we uncovered a BAP1 cancer syndrome kindred of ~80,000 descendants with a core of 106 individuals, whose members descend from a couple born in Germany in the early 1700s who immigrated to North America. Their descendants spread throughout the country with mutation carriers affected by multiple malignancies. Our data show that, once a proband is identified, extended analyses of these kindreds, using genomic and genealogical studies to identify the most recent common ancestor, allow investigators to uncover additional branches of the family that may carry BAP1 mutations. Using this knowledge, we have identified new branches of this family carrying BAP1 mutations. We have also implemented early-detection strategies that help identify cancers at early-stage, when they can be cured (melanomas) or are more susceptible to therapy (MM and other malignancies).

Meta-analysis of genome-wide association studies of adult height in East Asians identifies 17 novel loci.

Human height is associated with risk of multiple diseases and is profoundly determined by an individual's genetic makeup and shows a high degree of ethnic heterogeneity. Large-scale genome-wide association (GWA) analyses of adult height in Europeans have identified nearly 180 genetic loci. A recent study showed high replicability of results from Europeans-based GWA studies in Asians; however, population-specific loci may exist due to distinct linkage disequilibrium patterns. We carried out a GWA meta-analysis in 93 926 individuals from East Asia. We identified 98 loci, including 17 novel and 81 previously reported loci, associated with height at P < 5 × 10(-8), together explaining 8.89% of phenotypic variance. Among the newly identified variants, 10 are commonly distributed (minor allele frequency, MAF > 5%) in Europeans, with comparable frequencies with in Asians, and 7 single-nucleotide polymorphisms are with low frequency (MAF < 5%) in Europeans. In addition, our data suggest that novel biological pathway such as the protein tyrosine phosphatase family is involved in regulation of height. The findings from this study considerably expand our knowledge of the genetic architecture of human height in Asians.

A genome-wide association analysis identifies NMNAT2 and HCP5 as susceptibility loci for Kawasaki disease.

Kawasaki disease (KD), a systemic vasculitis of infants and children, manifests as fever and mucocutaneous inflammation. Although its etiology is largely unknown, the epidemiological data suggest that genetic factors are important in KD susceptibility. To identify genetic variants influencing KD susceptibility, we performed a genome-wide association study (GWAS) and replication study using a total of 915 children with KD and 4553 controls in the Korean population. Six single-nucleotide polymorphisms (SNPs) in three loci were associated significantly with KD susceptibility (P<1.0 × 10-5), including the previously reported BLK locus (rs6993775, odds ratio (OR)=1.52, P=2.52 × 10-11). The other two loci were newly identified: NMNAT2 on chromosome 1q25.3 (rs2078087, OR=1.33, P=1.15 × 10-6) and the human leukocyte antigen (HLA) region on chromosome 6p21.3 (HLA-C, HLA-B, MICA and HCP5) (rs9380242, rs9378199, rs9266669 and rs6938467; OR=1.33-1.51, P=8.93 × 10-6 to 5.24 × 10-8). Additionally, SNP rs17280682 in NLRP14 was associated significantly with KD with a family history (18 cases vs 4553 controls, OR=6.76, P=5.46 × 10-6). These results provide new insights into the pathogenesis and pathophysiology of KD.

Genome-wide association study identifies three novel loci for type 2 diabetes.

Although over 60 loci for type 2 diabetes (T2D) have been identified, there still remains a large genetic component to be clarified. To explore unidentified loci for T2D, we performed a genome-wide association study (GWAS) of 6 209 637 single-nucleotide polymorphisms (SNPs), which were directly genotyped or imputed using East Asian references from the 1000 Genomes Project (June 2011 release) in 5976 Japanese patients with T2D and 20 829 nondiabetic individuals. Nineteen unreported loci were selected and taken forward to follow-up analyses. Combined discovery and follow-up analyses (30 392 cases and 34 814 controls) identified three new loci with genome-wide significance, which were MIR129-LEP [rs791595; risk allele = A; risk allele frequency (RAF) = 0.080; P = 2.55 × 10(-13); odds ratio (OR) = 1.17], GPSM1 [rs11787792; risk allele = A; RAF = 0.874; P = 1.74 × 10(-10); OR = 1.15] and SLC16A13 (rs312457; risk allele = G; RAF = 0.078; P = 7.69 × 10(-13); OR = 1.20). This study demonstrates that GWASs based on the imputation of genotypes using modern reference haplotypes such as that from the 1000 Genomes Project data can assist in identification of new loci for common diseases.

A meta-analysis identifies adolescent idiopathic scoliosis association with LBX1 locus in multiple ethnic groups.

BACKGROUND: Adolescent idiopathic scoliosis (AIS) is a common rotational deformity of the spine that presents in children worldwide, yet its etiology is poorly understood. Recent genome-wide association studies (GWAS) have identified a few candidate risk loci. One locus near the chromosome 10q24.31 LBX1 gene (OMIM #604255) was originally identified by a GWAS of Japanese subjects and replicated in additional Asian populations. To extend this result, and to create larger AIS cohorts for the purpose of large-scale meta-analyses in multiple ethnicities, we formed a collaborative group called the International Consortium for Scoliosis Genetics (ICSG). METHODS: Here, we report the first ICSG study, a meta-analysis of the LBX1 locus in six Asian and three non-Asian cohorts. RESULTS: We find significant evidence for association of this locus with AIS susceptibility in all nine cohorts. Results for seven cohorts containing both genders yielded P=1.22×10-43 for rs11190870, and P=2.94×10-48 for females in all nine cohorts. Comparing the regional haplotype structures for three populations, we refined the boundaries of association to a ∼25 kb block encompassing the LBX1 gene. The LBX1 protein, a homeobox transcription factor that is orthologous to the Drosophila ladybird late gene, is involved in proper migration of muscle precursor cells, specification of cardiac neural crest cells, and neuronal determination in developing neural tubes. CONCLUSIONS: Our results firmly establish the LBX1 region as the first major susceptibility locus for AIS in Asian and non-Hispanic white groups, and provide a platform for larger studies in additional ancestral groups.

Immuno-genomic analysis reveals eosinophilic feature and favorable prognosis of female non-smoking esophageal squamous cell carcinomas.

Most of esophageal squamous cell carcinoma (ESCC) develop in smoking males in Japan, but the genomic etiology and immunological characteristics of rare non-smoking female ECSS remain unclear. To elucidate the genomic and immunological features of ESCC in non-smoking females, we analyzed whole-genome or transcriptome sequencing data from 94 ESCCs, including 20 rare non-smoking female cases. In addition, 31,611 immune cells were extracted from four ESCC tissues and subject to single-cell RNA-seq. We compared their immuno-genomic and microbiome profiles between non-smoking female and smoking ESCCs. Non-smoking females showed much better prognosis. Whole-genome sequencing analysis showed no significant differences in driver genes or copy number alterations depending on smoking status. The mutational signatures specifically observed in non-smoking females ESCC could be attributed to aging. Immune profiling from RNA-seq revealed that ESCC in non-smoking females had high tumor microenvironment signatures and a high abundance of eosinophils with a favorable prognosis. Single-cell RNA-sequencing of intratumor immune cells revealed gender differences of eosinophils and their activation in female cases. ESCCs in non-smoking females have age-related mutational signatures and gender-specific tumor immune environment with eosinophils, which is likely to contribute to their favorable prognosis.

What types of nursing facilities are more likely to adopt a pharmacist's medication review recommendations?

BACKGROUND: Few studies in the literature have analyzed the determinants of pharmacist drug therapy recommendations in nursing facility settings, and those that have focus primarily on accepted/rejected recommendations by disease state. OBJECTIVE: To identify the set of nursing facility characteristics that are more likely to adopt a pharmacist's medication review recommendations. DESIGN: Cross-sectional, retrospective methods are used to examine 53 licensed nursing facilities receiving medication review services from a small independent consultant pharmacist practice with no ties to vendor pharmacy functions. SETTING: Nursing facilities in rural areas of central and western Minnesota in 2008. INTERVENTION: Medication review services. MAIN OUTCOME MEASURES: The number of recommendations made and accepted, which are aggregated to the level of the nursing facility. Poisson regression analysis is used to identify those nursing facility characteristics that predict total recommendations and total accepted recommendations. Data obtained from Medicare's Web site on each nursing facility's operating characteristics and quality indicators serve as covariates. RESULTS: At the 5% level, patient census (positively), greater certified nursing assistant staffing hours (positively), multisite facilities (positively), resident residency councils (negatively), and greater perceptions of registered nurse quality (negatively) predict a greater number of recommendations. Patient census (positively), greater licensed practical nurse staffing (negatively), having residency councils (negatively), and greater perceptions of registered nurse quality (negatively) significantly predict the number of accepted and implemented recommendations. CONCLUSION: Institutional specific factors, most notably, quality-of-care indicators, may affect a nursing facility's acceptance of a pharmacist's drug therapy review.

Genotype-phenotype mapping of a patient-derived lung cancer organoid biobank identifies NKX2-1-defined Wnt dependency in lung adenocarcinoma.

Human lung cancer is a constellation of tumors with various histological and molecular properties. To build a preclinical platform that covers this broad disease spectrum, we obtained lung cancer specimens from multiple sources, including sputum and circulating tumor cells, and generated a living biobank consisting of 43 lines of patient-derived lung cancer organoids. The organoids recapitulated the histological and molecular hallmarks of the original tumors. Phenotypic screening of niche factor dependency revealed that EGFR mutations in lung adenocarcinoma are associated with the independence from Wnt ligands. Gene engineering of alveolar organoids reveals that constitutive activation of EGFR-RAS signaling provides Wnt independence. Loss of the alveolar identity gene NKX2-1 confers Wnt dependency, regardless of EGFR signal mutation. Sensitivity to Wnt-targeting therapy can be stratified by the expression status of NKX2-1. Our results highlight the potential of phenotype-driven organoid screening and engineering for the fabrication of therapeutic strategies to combat cancer.

Comparative analyses define differences between BHD-associated renal tumour and sporadic chromophobe renal cell carcinoma.

BACKGROUND: Birt-Hogg-Dubé (BHD) syndrome, caused by germline alteration of folliculin (FLCN) gene, develops hybrid oncocytic/chromophobe tumour (HOCT) and chromophobe renal cell carcinoma (ChRCC), whereas sporadic ChRCC does not harbor FLCN alteration. To date, molecular characteristics of these similar histological types of tumours have been incompletely elucidated. METHODS: To elucidate renal tumourigenesis of BHD-associated renal tumours and sporadic renal tumours, we conducted whole genome sequencing (WGS) and RNA-sequencing (RNA-seq) of sixteen BHD-associated renal tumours from nine unrelated BHD patients, twenty-one sporadic ChRCCs and seven sporadic oncocytomas. We then compared somatic mutation profiles with FLCN variants and RNA expression profiles between BHD-associated renal tumours and sporadic renal tumours. FINDINGS: RNA-seq analysis revealed that BHD-associated renal tumours and sporadic renal tumours have totally different expression profiles. Sporadic ChRCCs were clustered into two distinct clusters characterized by L1CAM and FOXI1 expressions, molecular markers for renal tubule subclasses. Increased mitochondrial DNA (mtDNA) copy number with fewer variants was observed in BHD-associated renal tumours compared to sporadic ChRCCs. Cell-of-origin analysis using WGS data demonstrated that BHD-associated renal tumours and sporadic ChRCCs may arise from different cells of origin and second hit FLCN alterations may occur in early third decade of life in BHD patients. INTERPRETATION: These data further our understanding of renal tumourigenesis of these two different types of renal tumours with similar histology. FUNDING: This study was supported by JSPS KAKENHI Grants, RIKEN internal grant, and the Intramural Research Program of the National Institutes of Health (NIH), National Cancer Institute (NCI), Center for Cancer Research.

Immuno-genomic profiling of biopsy specimens predicts neoadjuvant chemotherapy response in esophageal squamous cell carcinoma.

Esophageal squamous cell carcinoma (ESCC) is one of the most aggressive cancers and is primarily treated with platinum-based neoadjuvant chemotherapy (NAC). Some ESCCs respond well to NAC. However, biomarkers to predict NAC sensitivity and their response mechanism in ESCC remain unclear. We perform whole-genome sequencing and RNA sequencing analysis of 141 ESCC biopsy specimens before NAC treatment to generate a machine-learning-based diagnostic model to predict NAC reactivity in ESCC and analyzed the association between immunogenomic features and NAC response. Neutrophil infiltration may play an important role in ESCC response to NAC. We also demonstrate that specific copy-number alterations and copy-number signatures in the ESCC genome are significantly associated with NAC response. The interactions between the tumor genome and immune features of ESCC are likely to be a good indicator of therapeutic capability and a therapeutic target for ESCC, and machine learning prediction for NAC response is useful.

High-power broadband laser source tunable from 3.0 µm to 4.4 µm based on a femtosecond Yb:fiber oscillator.

We describe a tunable broadband mid-IR laser source based on difference-frequency mixing of a 100 MHz femtosecond Yb:fiber laser oscillator and a Raman-shifted soliton generated with the same laser. The resulting light is tunable over 3.0 µm to 4.4 µm, with a FWHM bandwidth of 170 nm and maximum average output power up to 125 mW. The noise and coherence properties of this source are also investigated and described.