Low-coverage Whole-genome Sequencing Research Articles

Correlation of Immunological and Molecular Profiles with Response to Crizotinib in Alveolar Soft Part Sarcoma: An Exploratory Study Related to the EORTC 90101 "CREATE" Trial.

Alveolar soft part sarcoma (ASPS) is a rare subtype of soft tissue sarcoma characterized by an unbalanced translocation, resulting in ASPSCR1-TFE3 fusion that transcriptionally upregulates MET expression. The European Organization for Research and Treatment of Cancer (EORTC) 90101 “CREATE” phase II trial evaluated the MET inhibitor crizotinib in ASPS patients, achieving only limited antitumor activity. We performed a comprehensive molecular analysis of ASPS tissue samples collected in this trial to identify potential biomarkers correlating with treatment outcome. A tissue microarray containing 47 ASPS cases was used for the characterization of the tumor microenvironment using multiplex immunofluorescence. DNA isolated from 34 available tumor samples was analyzed to detect recurrent gene copy number alterations (CNAs) and mutations by low-coverage whole-genome sequencing and whole-exome sequencing. Pathway enrichment analysis was used to identify diseased-associated pathways in ASPS sarcomagenesis. Kaplan–Meier estimates, Cox regression, and the Fisher’s exact test were used to correlate histopathological and molecular findings with clinical data related to crizotinib treatment, aiming to identify potential factors associated with patient outcome. Tumor microenvironment characterization showed the presence of PD-L1 and CTLA-4 in 10 and 2 tumors, respectively, and the absence of PD-1 in all specimens. Apart from CD68, other immunological markers were rarely expressed, suggesting a low level of tumor-infiltrating lymphocytes in ASPS. By CNA analysis, we detected a number of broad and focal alterations. The most common alteration was the loss of chromosomal region 1p36.32 in 44% of cases. The loss of chromosomal regions 1p36.32, 1p33, 1p22.2, and 8p was associated with shorter progression-free survival. Using whole-exome sequencing, 13 cancer-associated genes were found to be mutated in at least three cases. Pathway enrichment analysis identified genetic alterations in NOTCH signaling, chromatin organization, and SUMOylation pathways. NOTCH4 intracellular domain dysregulation was associated with poor outcome, while inactivation of the beta-catenin/TCF complex correlated with improved outcome in patients receiving crizotinib. ASPS is characterized by molecular heterogeneity. We identify genetic aberrations potentially predictive of treatment outcome during crizotinib therapy and provide additional insights into the biology of ASPS, paving the way to improve treatment approaches for this extremely rare malignancy.

Cerebrospinal Fluid Cell-Free DNA-Based Detection of High Level of Genomic Instability Is Associated With Poor Prognosis in NSCLC Patients With Leptomeningeal Metastases.

IntroductionLeptomeningeal metastasis (LM) commonly occurs in non-small cell lung cancer (NSCLC) patients and has a poor prognosis. Due to limited access to leptomeningeal lesions, the genetic characteristics of LM have not been explored to date. Cerebrospinal fluid (CSF) may be the most representative liquid biopsy medium to obtain genomic information from LM in NSCLC.MethodsCSF biopsies and matched peripheral blood biopsies were collected from 33 NSCLC patients with LM. We profiled genetic alterations from LM by comparing CSF cell-free DNA (cfDNA) with plasma cfDNA. Somatic mutations were examined using targeted sequencing. Genomic instability was analyzed by low-coverage whole-genome sequencing (WGS).ResultsDriver mutations were detected in 100% of CSF cfDNA with much higher variant allele frequency than that in matched plasma cfDNA (57.5%). Furthermore, we found that the proportions of CSF cfDNA fragments below 150 bp were significantly higher than those in plasma cfDNA. These findings indicate enrichment of circulating tumor DNA (ctDNA) in CSF and explain the high sensitivity of mutation detection in the CSF. The absence of some mutations in CSF cfDNA—especially the first-/second-generation mutation T790M, which confers resistance to epidermal growth factor receptor (EGFR)-Tyrosine kinase inhibitors (TKIs)—that were present in plasma cfDNA samples indicates different mechanisms of cancer evolution between LM and extracranial lesions. In addition, 86.6% of CSF ctDNA samples revealed high levels of genomic instability compared with 2.5% in plasma cfDNA samples. A higher number of large-scale state transitions (LSTs) in CSF cfDNA were associated with a shorter overall survival (OS).ConclusionOur results suggest that LM and extracranial lesions develop independently. Both CSF cfDNA genetic profiling and plasma cfDNA genetic profiling are necessary for clinical decision-making for NSCLC patients with LM. Through CSF-based low-coverage WGS, a high level of LSTs was identified as a potential biomarker of poor prognosis.

Nyssorhynchus darlingi genome-wide studies related to microgeographic dispersion and blood-seeking behavior

BackgroundIn Brazil, malaria is concentrated in the Amazon Basin, where more than 99% of the annual cases are reported. The main goal of this study was to investigate the population structure and genetic association of the biting behavior of Nyssorhynchus (also known as Anopheles) darlingi, the major malaria vector in the Amazon region of Brazil, using low-coverage genomic sequencing data.MethodsSamples were collected in the municipality of Mâncio Lima, Acre state, Brazil between 2016 and 2017. Different approaches using genotype imputation and no gene imputation for data treatment and low-coverage sequencing genotyping were performed. After the samples were genotyped, population stratification analysis was performed.ResultsWeak but statistically significant stratification signatures were identified between subpopulations separated by distances of approximately 2–3 km. Genome-wide association studies (GWAS) were performed to compare indoor/outdoor biting behavior and blood-seeking at dusk/dawn. A statistically significant association was observed between biting behavior and single nucleotide polymorphism (SNP) markers adjacent to the gene associated with cytochrome P450 (CYP) 4H14, which is associated with insecticide resistance. A statistically significant association between blood-seeking periodicity and SNP markers adjacent to genes associated with the circadian cycle was also observed.ConclusionThe data presented here suggest that low-coverage whole-genome sequencing with adequate processing is a powerful tool to genetically characterize vector populations at a microgeographic scale in malaria transmission areas, as well as for use in GWAS. Female mosquitoes entering houses to take a blood meal may be related to a specific CYP4H14 allele, and female timing of blood-seeking is related to circadian rhythm genes.Graphical

Cost-effective low-coverage whole-genome sequencing assay for the risk stratification of gastric cancer.

BACKGROUNDGastric cancer (GC), a multifactorial disease, is caused by pathogens, such as Helicobacter pylori (H. pylori) and Epstein-Barr virus (EBV), and genetic components. AIMTo investigate microbiomes and host genome instability by cost-effective, low-coverage whole-genome sequencing, as biomarkers for GC subtyping. METHODSSamples from 40 GC patients were collected from Taizhou Hospital, Zhejiang Province, affiliated with Wenzhou Medical University. DNA from the samples was subjected to low-coverage whole-genome sequencing with a median genome coverage of 1.86 × (range: 1.03 × to 3.17 ×) by Illumina × 10, followed by copy number analyses using a customized bioinformatics workflow ultrasensitive chromosomal aneuploidy detector. RESULTSOf the 40 GC samples, 20 (50%) were found to be enriched with microbiomes. EBV DNA was detected in 5 GC patients (12.5%). H. pylori DNA was found in 15 (37.5%) patients. The other 20 (50%) patients were found to have relatively higher genomic instability. Copy number amplifications of the oncogenes, ERBB2 and KRAS, were found in 9 (22.5%) and 7 (17.5%) of the GC samples, respectively. EBV enrichment was found to be associated with tumors in the gastric cardia and fundus. H. pylori enrichment was found to be associated with tumors in the pylorus and antrum. Tumors with elevated genomic instability showed no localization and could be observed in any location. Additionally, H. pylori-enriched GC was found to be associated with the Borrmann type II/III and gastritis history. EBV-enriched GC was not associated with gastritis. No statistically significant correlation was observed between genomic instability and gastritis. Furthermore, these three different molecular subtypes showed distinct survival outcomes (P = 0.019). EBV-positive tumors had the best prognosis, whereas patients with high genomic instability (CIN+) showed the worst survival. Patients with H. pylori infection showed intermediate prognosis compared with the other two subtypes.CONCLUSIONThus, using low-coverage whole-genome sequencing, GC can be classified into three categories based on disease etiology; this classification may prove useful for GC diagnosis and precision medicine.

Paroxysmal Kinesigenic Dyskinesia Caused by 16p11.2 Microdeletion and Related Clinical Features.

Background and ObjectivesIsolated paroxysmal kinesigenic dyskinesia (PKD) is mainly caused by PRRT2 variants and TMEM151A variants. Patients with proximal 16p11.2 microdeletion (16p11.2MD) (including PRRT2) often have neurodevelopmental phenotypes, whereas a few patients have PKD. Here, we aimed to identify 16p11.2MD in patients with PKD and describe the related phenotypes.MethodsWhole-exome sequencing and bioinformatics analysis of copy number variant (CNV) were performed in patients with PKD carrying neither PRRT2 nor TMEM151A variant. Quantitative PCR and low-coverage whole-genome sequencing verified the CNV.ResultsWe identified 9 sporadic patients with PKD and 16p11.2MD (∼535 kb), accounting for 9.6% (9/94) of our patients. Together with 9 previously reported patients with PKD and 16p11.2MD, we found that 16p11.2MD was de novo in 11 of 12 tested patients and inherited from a parent in the other patient. And 80% (12/15) of these patients had a mild language delay, 64.3% (9/14) had compromised learning ability, 42.9% (6/14) had a mild motor delay, and 50% (6/12) had abnormal neuroimaging findings. No severe autism disorders were observed.DiscussionMild developmental problems may be overlooked. A detailed inquiry of developmental history and CNV testing are necessary to distinguish patients with 16p11.2MD from isolated PKD.

Abstract PD14-05: Portraying tumor evolution of lobular breast cancer through phylogenetic analysis

Abstract Background: Despite therapeutic advancements, a substantial number of patients with treated early-stage invasive lobular carcinoma (ILC) of the breast will still relapse as late as 20 years after the initial tumor diagnosis, thus making follow-up a challenging task. Here, using patient-matched primary and metastatic tumors from patients with ILC, we aimed to decipher the evolution of lobular breast cancer through time and explore metastatic dissemination events and underlying heterogeneity. Materials and Methods: We retrospectively identified 38 patients with metastatic ILC, of which 6 with de novo metastatic disease. In total, 99 formalin-fixed paraffin embedded tissue samples of primary tumor and metastasis, together with well-annotated clinicopathological data, were available. Low coverage whole-genome sequencing was performed to infer somatic copy number aberrations (CNA). Phylogenetic reconstruction was performed with Integer Linear Program for the Copy-Number Tree Problem algorithm. Results: Reconstruction of the phylogenetic trees showed different patterns of dissemination. For patients with multiple metastatic samples sequenced (4/4 patients), the metastases were found to be seeded from a common metastatic precursor. For a subset of patients with multiple primaries (5/10 patients), distant metastases were seeded from specific clones inside the primary tumor. Detailed analysis of known driver genes showed that alterations in common driver genes such as TP53, MYC and CCND1 were shared between primary tumor and metastases in a large proportion of the patients (47%, 32% and 21% of patients respectively) whilst CNA in other driver genes including FGFR2, PTEN or AKT2 among others, were frequently private in the metastases. Interestingly, we observed that in approximately 30% of the patients, the metastases harbored less CNA than the matched primary tumor, resulting in shorter metastatic branches and possibly indicating a decelerated tumor progression. Further analysis of the genomic and clinicopathological characteristics showed that these patients were mostly de novo metastatic (p-value=0.002), metastasizing primarily to the liver or bone and had significantly more aberrations in NF1, AKT1 and FGFR (p-value=0.04, 0.01 and 0.03 respectively), when compared to the rest of the cohort. Finally, most of the patients (63%) had a lower number of CNA accumulated after the metastatic lineage separation compared to the common trunk, implying a late dissemination event. Conclusions: In this study, we described different metastatic dissemination events taking place in ILC. Our analyses allowed us to identify a subset of tumors characterized by slow tumor progression, indicating potential tumor dormancy, as well as evidence of late dissemination evidence for most of the patients, thus highlighting the importance of early detection through cancer screening. These findings further shed light on how ILCs evolves through time and could potentially influence the clinical management of invasive lobular breast cancer. Citation Format: Danai Fimereli, David Venet, Mattia Rediti, David N Brown, Bram Boeckx, Marion Maetens, Samira Majjaj, Ghizlane Rouas, Maria Capra, Giuseppina Bonizzi, Federica Contaldo, Christine Galant, Martine Piccart, Giancarlo Pruneri, Denis Larsimont, Diether Lambrechts, Christine Desmedt, Françoise Rothé, Christos Sotiriou. Portraying tumor evolution of lobular breast cancer through phylogenetic analysis [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr PD14-05.

Ultra Low-Coverage Whole-Genome Sequencing as an Alternative to Genotyping Arrays in Genome-Wide Association Studies.

An array-based genotyping approach has been the standard practice for genome-wide association studies (GWASs); however, as sequencing costs plummet over the past years, ultra low-coverage whole-genome sequencing (ulcWGS <0.5× coverage) has emerged as a promising alternative that provides superior genomic coverage with substantial reduction of genotyping cost. To evaluate the potential utility of ulcWGS, we performed a whole-genome sequencing (WGS) of 72 European individuals to a target coverage of 0.4× and compared its performance with the widely used Infinium Global Screening Multi-Disease Array (GSA-MD). We showed that the number of variants captured by ulcWGS is comparable with imputed GSA-MD platform, particularly for low-frequency (95.5%) and common variants (99.9%), with high imputation R2 accuracy (mean 0.93 for SNPs and 0.86 for indels). Using deep-coverage 30× WGS as the “truth” genotypes, we found that ulcWGS has higher overall nonreference genotype concordance compared with imputed GSA-MD for both SNPs (0.90 vs. 0.88) and indels (0.86 vs. 0.83). In addition, ulcWGS proved to be as sensitive as the genotyping-based method in sex imputation and ancestry prediction producing similar principal component (PC) scores. Our findings provide important evidence that the cost efficient ulcWGS of <0.5× generates high genotype accuracy, outperforming the standard genotyping arrays, making it an attractive alternative to the array-based method in next-generation GWAS design.

Feasibility of Screening for Chromosome 15 Imprinting Disorders in 16 579 Newborns by Using a Novel Genomic Workflow

Newborn screening for Angelman syndrome (AS), Prader-Willi syndrome (PWS), and chromosome 15 duplication syndrome (Dup15q) may lead to benefit from early diagnosis and treatment. To examine the feasibility of newborn screening for these chromosome 15 imprinting disorders at population scale. In this diagnostic study, the validation data set for the first-tier SNRPN test, called methylation-specific quantitative melt analysis (MS-QMA), included 109 PWS, 48 AS, 9 Dup15q, and 1190 population control newborn blood spots (NBS) and peripheral tissue samples from participants recruited from January 2000 to December 2016. The test data set included NBS samples from 16 579 infants born in 2011. Infants with an NBS identified as positive for PWS, AS, or Dup15q by the first-tier test were referred for droplet digital polymerase chain reaction, real-time polymerase chain reaction, and low-coverage whole-genome sequencing for confirmatory testing. Data analyses were conducted between February 12, 2015, and August 15, 2020. In the validation data set, the median age for the 77 patients with PWS was 3.00 years (IQR, 0.01-44.50 years); for the 46 patients with AS, 2.76 years (IQR, 0.028 to 49.00 years); and for the 9 patients with Dup15q, 4.00 years (IQR, 1.00 to 28.00 years). Thirty-eight patients (51.4%) in the PWS group, 20 patients (45.5%) in the AS group, and 6 patients (66.7%) in the Dup15q group who had sex reported were male. The validation data set showed MS-QMA sensitivity of 99.0% for PWS, 93.8% for AS, and 77.8% for Dup15q; specificity of 100% for PWS, AS, and Dup15q; positive predictive and negative predictive values of 100% for PWS and AS; and a positive predictive value of 87.5% and negative predictive value of 100% for Dup15q. In the test data set of NBS samples from 16 579 infants, 92 had a positive test result using a methylation ratio cut-off of 3 standard deviations from the mean. Of these patients, 2 were confirmed to have PWS; 2, AS; and 1, maternal Dup15q. With the use of more conservative PWS- and AS-specific thresholds for positive calls from the validation data set, 9 positive NBS results were identified by MS-QMA in this cohort. The 2 PWS and 2 AS calls were confirmed by second-tier testing, but the 1 Dup15q case was not confirmed. Together, these results provided prevalence estimates of 1 in 8290 for both AS and PWS and 1 in 16 579 for maternal Dup15q, with positive predictive values for first-tier testing at 67.0% for AS, 33.0% for PWS, and 44.0% for combined detection of chromosome 15 imprinting disorders for the validation data set. The findings of this diagnostic study suggest that it is feasible to screen for all chromosome 15 imprinting disorders using SNRPN methylation analysis, with 5 individuals identified with these disorders out of 16 579 infants screened.

Comparison of Genotype Imputation for SNP Array and Low-Coverage Whole-Genome Sequencing Data.

Genotype imputation is the term used to describe the process of inferring unobserved genotypes in a sample of individuals. It is a key step prior to a genome-wide association study (GWAS) or genomic prediction. The imputation accuracy will directly influence the results from subsequent analyses. In this simulation-based study, we investigate the accuracy of genotype imputation in relation to some factors characterizing SNP chip or low-coverage whole-genome sequencing (LCWGS) data. The factors included the imputation reference population size, the proportion of target markers /SNP density, the genetic relationship (distance) between the target population and the reference population, and the imputation method. Simulations of genotypes were based on coalescence theory accounting for the demographic history of pigs. A population of simulated founders diverged to produce four separate but related populations of descendants. The genomic data of 20,000 individuals were simulated for a 10-Mb chromosome fragment. Our results showed that the proportion of target markers or SNP density was the most critical factor affecting imputation accuracy under all imputation situations. Compared with Minimac4, Beagle5.1 reproduced higher-accuracy imputed data in most cases, more notably when imputing from the LCWGS data. Compared with SNP chip data, LCWGS provided more accurate genotype imputation. Our findings provided a relatively comprehensive insight into the accuracy of genotype imputation in a realistic population of domestic animals.

Highlights from Recent Cancer Literature

Highlights from Recent Cancer Literature

Histopathological and Molecular Profiling of Clear Cell Sarcoma and Correlation with Response to Crizotinib: An Exploratory Study Related to EORTC 90101 "CREATE" Trial.

Simple SummaryClear cell sarcoma (CCSA) is a rare subtype of soft tissue sarcoma characterized by EWSR1 rearrangement and subsequent MET upregulation. The European Organisation for Research and Treatment of Cancer 90101 phase II trial evaluated the MET inhibitor crizotinib in CCSA but resulted in only sporadic responses. The aim of this exploratory study was to identify the molecular alterations potentially relevant for the treatment outcome by using archival CCSA samples and trial-related clinical data. We characterized MET signaling and revealed an infrequent activation of MET, which may explain the lack of response to crizotinib in the disease cohort. Based on sequencing analyses, we discovered copy number alterations, mutations and dysregulated pathways with potentially predictive or prognostic values for patients’ outcomes. This work describes the molecular heterogeneity in CCSA and provides deep insight into the biology of this ultra-rare malignancy, which may potentially lead to better therapeutic approaches.Clear cell sarcoma (CCSA) is characterized by a chromosomal translocation leading to EWSR1 rearrangement, resulting in aberrant transcription of multiple genes, including MET. The EORTC 90101 phase II trial evaluated the MET inhibitor crizotinib in CCSA but resulted in only sporadic responses. We performed an in-depth histopathological and molecular analysis of archival CCSA samples to identify alterations potentially relevant for the treatment outcome. Immunohistochemical characterization of MET signaling was performed using a tissue microarray constructed from 32 CCSA cases. The DNA from 24 available tumor specimens was analyzed by low-coverage whole-genome sequencing and whole-exome sequencing for the detection of recurrent copy number alterations (CNAs) and mutations. A pathway enrichment analysis was performed to identify the pathways relevant for CCSA tumorigenesis. Kaplan–Meier estimates and Fisher’s exact test were used to correlate the molecular findings with the clinical features related to crizotinib treatment, aiming to assess a potential association with the outcomes. The histopathological analysis showed the absence of a MET ligand and MET activation, with the presence of MET itself in most of cases. However, the expression/activation of MET downstream molecules was frequently observed, suggesting the role of other receptors in CCSA signal transduction. Using sequencing, we detected a number of CNAs at the chromosomal arm and region levels. The most common alteration was a gain of 8q24.21, observed in 83% of the cases. The loss of chromosomes 9q and 12q24 was associated with shorter survival. Based on exome sequencing, 40 cancer-associated genes were found to be mutated in more than one sample, with SRGAP3 and KMT2D as the most common alterations (each in four cases). The mutated genes encoded proteins were mainly involved in receptor tyrosine kinase signaling, polymerase-II transcription, DNA damage repair, SUMOylation and chromatin organization. Disruption in chromatin organization was correlated with longer progression-free survival in patients receiving crizotinib. Conclusions: The infrequent activation of MET may explain the lack of response to crizotinib observed in the majority of cases in the clinical trial. Our work describes the molecular heterogeneity in CCSA and provides further insight into the biology of this ultra-rare malignancy, which may potentially lead to better therapeutic approaches for CCSA.

Genome-Guided Discovery of Natural Products through Multiplexed Low-Coverage Whole-Genome Sequencing of Soil Actinomycetes on Oxford Nanopore Flongle.

ABSTRACTGenome mining is an important tool for discovery of new natural products; however, the number of publicly available genomes for natural product-rich microbes such as actinomycetes, relative to human pathogens with smaller genomes, is small. To obtain contiguous DNA assemblies and identify large (ca. 10 to greater than 100 kb) biosynthetic gene clusters (BGCs) with high GC (>70%) and high-repeat content, it is necessary to use long-read sequencing methods when sequencing actinomycete genomes. One of the hurdles to long-read sequencing is the higher cost. In the current study, we assessed Flongle, a recently launched platform by Oxford Nanopore Technologies, as a low-cost DNA sequencing option to obtain contiguous DNA assemblies and analyze BGCs. To make the workflow more cost-effective, we multiplexed up to four samples in a single Flongle sequencing experiment while expecting low-sequencing coverage per sample. We hypothesized that contiguous DNA assemblies might enable analysis of BGCs even at low sequencing depth. To assess the value of these assemblies, we collected high-resolution mass spectrometry data and conducted a multi-omics analysis to connect BGCs to secondary metabolites. In total, we assembled genomes for 20 distinct strains across seven sequencing experiments. In each experiment, 50% of the bases were in reads longer than 10 kb, which facilitated the assembly of reads into contigs with an average N50 value of 3.5 Mb. The programs antiSMASH and PRISM predicted 629 and 295 BGCs, respectively. We connected BGCs to metabolites for N,N-dimethyl cyclic-di-tryptophan, two novel lasso peptides, and three known actinomycete-associated siderophores, namely, mirubactin, heterobactin, and salinichelin.IMPORTANCE Short-read sequencing of GC-rich genomes such as those from actinomycetes results in a fragmented genome assembly and truncated biosynthetic gene clusters (often 10 to >100 kb long), which hinders our ability to understand the biosynthetic potential of a given strain and predict the molecules that can be produced. The current study demonstrates that contiguous DNA assemblies, suitable for analysis of BGCs, can be obtained through low-coverage, multiplexed sequencing on Flongle, which provides a new low-cost workflow ($30 to 40 per strain) for sequencing actinomycete strain libraries.

Chromosome-level assembly reveals a putative Y-autosomal fusion in the sex determination system of the Greenland Halibut (Reinhardtius hippoglossoides).

Despite the commercial importance of Greenland Halibut (Reinhardtius hippoglossoides), important gaps still persist in our knowledge of this species, including its reproductive biology and sex determination mechanism. Here, we combined single-molecule sequencing of long reads (Pacific Sciences) with chromatin conformation capture sequencing (Hi-C) data to assemble the first chromosome-level reference genome for this species. The high-quality assembly encompassed more than 598 Megabases (Mb) assigned to 1594 scaffolds (scaffold N50 = 25 Mb) with 96% of its total length distributed among 24 chromosomes. Investigation of the syntenic relationship with other economically important flatfish species revealed a high conservation of synteny blocks among members of this phylogenetic clade. Sex determination analysis revealed that similar to other teleost fishes, flatfishes also exhibit a high level of plasticity and turnover in sex determination mechanisms. A low-coverage whole-genome sequence analysis of 198 individuals revealed that Greenland Halibut possesses a male heterogametic XY system and several putative candidate genes implied in the sex determination of this species. Our study also suggests for the first time in flatfishes that a putative Y-autosomal fusion could be associated with a reduction of recombination typical of the early steps of sex chromosome evolution.

BIOM-36. SERIAL ASSESSMENT OF MEASURABLE RESIDUAL DISEASE IN MEDULLOBLASTOMA LIQUID BIOPSIES

Abstract Medulloblastoma is the most common malignant brain tumor in children. Despite the use of contemporary multi-modal therapy, up to one-third of patients will succumb to disease progression. Conventional response monitoring is based on magnetic resonance imaging and cerebrospinal fluid (CSF) cytology and a marker for measurable residual disease (MRD) is lacking. In this study, we show the clinical utility of profiling CSF-derived cell-free DNA (cfDNA) as a biomarker of MRD based on a sizable cohort of children with medulloblastoma enrolled on a prospective, multi-center, risk-adapted trial (SJMB03; NCT00085202). A total of 476 CSF samples were serially collected from 123 patients by lumbar puncture at post-operative baseline, during therapy, and at regular surveillance after completion of therapy. Using low-coverage whole-genome sequencing (lcWGS), tumor-associated copy-number alterations in CSF-derived cfDNA were investigated as an MRD surrogate for correlation with clinical features and outcomes. MRD was detected in post-operative baseline CSF samples for 85% and 54% of patients with metastatic and localized disease, respectively. The number of MRD-positive patients declined with therapy, yet those with persistent MRD detected at the end of therapy (20 of 68 CSFs at this timepoint) had significantly higher risk of disease progression (hazard ratio 8.94, 95%CI 4.10-19.49; log-rank p&amp;lt; 0.0001). Importantly, MRD-positivity from routine surveillance samples often preceded radiographic progression by more than 3 months. Comparative analyses of copy-number variations obtained from serial CSF samples enabled identification of relapse-specific alterations and early detection of relapse-dominant clones. Overall, our findings advocate for the routine incorporation of CSF-derived liquid biopsies in future medulloblastoma trials.

Haplotype-aware inference of human chromosome abnormalities

Extra or missing chromosomes-a phenomenon termed aneuploidy-frequently arise during human meiosis and embryonic mitosis and are the leading cause of pregnancy loss, including in the context of in vitro fertilization (IVF). While meiotic aneuploidies affect all cells and are deleterious, mitotic errors generate mosaicism, which may be compatible with healthy live birth. Large-scale abnormalities such as triploidy and haploidy also contribute to adverse pregnancy outcomes, but remain hidden from standard sequencing-based approaches to preimplantation genetic testing for aneuploidy (PGT-A). The ability to reliably distinguish meiotic and mitotic aneuploidies, as well as abnormalities in genome-wide ploidy, may thus prove valuable for enhancing IVF outcomes. Here, we describe a statistical method for distinguishing these forms of aneuploidy based on analysis of low-coverage whole-genome sequencing data, which is the current standard in the field. Our approach overcomes the sparse nature of the data by leveraging allele frequencies and linkage disequilibrium (LD) measured in a population reference panel. The method, which we term LD-informed PGT-A (LD-PGTA), retains high accuracy down to coverage as low as 0.05 × and at higher coverage can also distinguish between meiosis I and meiosis II errors based on signatures spanning the centromeres. LD-PGTA provides fundamental insight into the origins of human chromosome abnormalities, as well as a practical tool with the potential to improve genetic testing during IVF.

Low-Coverage Whole Genome Sequencing Using Laser Capture Microscopy with Combined Digital Droplet PCR: An Effective Tool to Study Copy Number and Kras Mutations in Early Lung Adenocarcinoma Development.

Defining detailed genomic characterization of early tumor progression is critical to identifying key regulators and pathways in carcinogenesis as potentially druggable targets. In human lung cancer, work to characterize early cancer development has mainly focused on squamous cancer, as the earliest lesions are more proximal in the airways and often accessible by repeated bronchoscopy. Adenocarcinomas are typically located distally in the lung, limiting accessibility for biopsy of pre-malignant and early stages. Mouse lung cancer models recapitulate many human genomic features and provide a model for tumorigenesis with pre-malignant atypical adenomatous hyperplasia and in situ adenocarcinomas often developing contemporaneously within the same animal. Here, we combined tissue characterization and collection by laser capture microscopy (LCM) with digital droplet PCR (ddPCR) and low-coverage whole genome sequencing (LC-WGS). ddPCR can be used to identify specific missense mutations in Kras (Kirsten rat sarcoma viral oncogene homolog, here focused on Kras Q61) and estimate the percentage of mutation predominance. LC-WGS is a cost-effective method to infer localized copy number alterations (CNAs) across the genome using low-input DNA. Combining these methods, the histological stage of lung cancer can be correlated with appearance of Kras mutations and CNAs. The utility of this approach is adaptable to other mouse models of human cancer.

Molecular genetic analysis and growth hormone response in patients with syndromic short stature

BackgroundSyndromic short stature is a genetic and phenotypic heterogeneous disorder with multiple causes. This study aims to identify genetic causes in patients with syndromic short stature of unknown cause and evaluate the efficacy of the growth hormone response.MethodsTrio-whole-exome sequencing was applied to identify pathogenic gene mutations in seven patents with short stature, multiple malformations, and/or intellectual disability. Whole-genome low-coverage sequencing was also performed to identify copy number variants in three patients with concurrent intellectual disability. Recombinant human growth hormone was administered to improve height in patients with an identified cause of syndromic short stature.ResultsOf the seven patients, three pathogenic/likely pathogenic gene mutations, including one FGFR3 mutation (c.1620C>A p.N540K), one novel GNAS mutation (c.2288C>T p.A763V), and one novel TRPS1 mutation (c.2527_c.2528dupTA p.S843fsX72), were identified in three patients. No copy number variants were identified in the three patients with concurrent intellectual disability. The proband with an FGFR3 mutation, a female 4 and 3/12 years of age, was diagnosed with hypochondroplasia. Long-acting growth hormone improved her height from 85.8 cm [− 5.05 standard deviation (SD)] to 100.4 cm (− 4.02 SD), and her increased height SD score (SDS) was 1.03 after 25 months of treatment. The proband with a GNAS mutation, a female 12 and 9/12 years of age, was diagnosed with pseudohypoparathyroidism Ia. After 14 months of treatment with short-acting growth hormone, her height improved from 139.3 cm (− 2.69 SD) to 145.0 cm (− 2.36 SD), and her increased height SDS was 0.33.ConclusionsTrio-whole-exome sequencing was an important approach to confirm genetic disorders in patients with syndromic short stature of unknown etiology. Short-term growth hormone was effective in improving height in patients with hypochondroplasia and pseudohypoparathyroidism Ia.

Serial assessment of measurable residual disease in medulloblastoma liquid biopsies.

Nearly one-third of children with medulloblastoma, a malignant embryonal tumor of the cerebellum, succumb to their disease. Conventional response monitoring by imaging and cerebrospinal fluid (CSF) cytology remains challenging, and a marker for measurable residual disease (MRD) is lacking. Here, we show the clinical utility of CSF-derived cell-free DNA (cfDNA) as a biomarker of MRD in serial samples collected from children with medulloblastoma (123 patients, 476 samples) enrolled on a prospective trial. Using low-coverage whole-genome sequencing, tumor-associated copy-number variations in CSF-derived cfDNA are investigated as an MRD surrogate. MRD is detected at baseline in 85% and 54% of patients with metastatic and localized disease, respectively. The number of MRD-positive patients declines with therapy, yet those with persistent MRD have significantly higher risk of progression. Importantly, MRD detection precedes radiographic progression in half who relapse. Our findings advocate for the prospective assessment of CSF-derived liquid biopsies in future trials for medulloblastoma.

A bioinformatic platform to integrate target capture and whole genome sequences of various read depths for phylogenomics.

The increasing availability of short‐read whole genome sequencing (WGS) provides unprecedented opportunities to study ecological and evolutionary processes. Although loci of interest can be extracted from WGS data and combined with target sequence data, this requires suitable bioinformatic workflows. Here, we test different assembly and locus extraction strategies and implement them into secapr, a pipeline that processes short‐read data into multilocus alignments for phylogenetics and molecular ecology analyses. We integrate the processing of data from low‐coverage WGS (<30×) and target sequence capture into a flexible framework, while optimizing de novo contig assembly and loci extraction. Specifically, we test different assembly strategies by contrasting their ability to recover loci from targeted butterfly protein‐coding genes, using four data sets: a WGS data set across different average coverages (10×, 5× and 2×) and a data set for which these loci were enriched prior to sequencing via target sequence capture. Using the resulting de novo contigs, we account for potential errors within contigs and infer phylogenetic trees to evaluate the ability of each assembly strategy to recover species relationships. We demonstrate that choosing multiple sizes of kmer simultaneously for assembly results in the highest yield of extracted loci from de novo assembled contigs, while data sets derived from sequencing read depths as low as 5× recovers the expected species relationships in phylogenetic trees. By making the tested assembly approaches available in the secapr pipeline, we hope to inspire future studies to incorporate complementary data and make an informed choice on the optimal assembly strategy.

Identification of a Rare Case With Nagashima-Type Palmoplantar Keratoderma and 18q Deletion Syndrome via Exome Sequencing and Low-Coverage Whole-Genome Sequencing.

Nagashima-type palmoplantar keratoderma (NPPK) is characterized by non-progressive, diffuse, and cross-gradient hyperkeratosis caused by mutations in the SERPINB7 gene on chromosome 18q21.33. Chromosome 18q deletion syndrome (18q- syndrome) is a terminal deletion or microdeletion syndrome characterized by intellectual disability and congenital malformations. This paper describes an 18-year-old man with palmoplantar keratoderma and diffuse white matter abnormalities in the brain. Trio-based exome sequencing (ES) revealed a suspected mosaic compound heterozygous mutation for c.796C>T (p.Arg266∗) in exon 8 inherited from the mother and a de novo exons 4–6 deletion of SERPINB7. Additional copy number variant (CNV) analysis of the ES data indicated a heterozygous gross deletion of 18q22.3-q23. The two SERPINB7 gene variants were verified by Sanger sequencing and quantitative real-time polymerase chain reaction (qRT-PCR). Finally, low-coverage whole-genome sequencing (WGS) confirmed the 18q22.3-q23 deletion and additionally detected a mosaic 18q21.33-q22.3 deletion, together explaining NPPK and the neurological phenotypes of the proband. The gross deletion of all exons of SERPINB7 was revealed for the first time. More rarely, c.796C>T (p.Arg266∗) was likely to be mosaic, while the exon deletion was mosaic. In conclusion, the combination of multiple molecular genetic testing methods provides comprehensive informative molecular findings and promotes the diagnosis of complex diseases, as in this case.