Clonal and "Intrinsic" Heterogeneity of Somatic Variants in Microsatellite-Stable Colorectal Carcinomas and Their Metastases.

To test the traditional model of tumor progression, Darwinian-type evolution, against the more recent Big Bang model, we selected 6 microsatellite-stable colorectal standard-type adenocarcinomas and their synchronous lymph node and liver metastases. Somatic genomic variants were identified by whole-exome sequencing (WES) of large tumor fragments from the primaries and 1 liver metastasis each, and used to design targeted resequencing next-generation sequencing (NGS) panels, 1 per case. Targeted deep resequencing (mean coverage, 2725; median, 2222) was performed with DNA from punch samples (1-mm tissue microarrayer needles) obtained from different regions of the primaries and their metastases. In total, 255 genomic variants were interrogated in 108 punch samples. Clonal heterogeneity was infrequent: a pattern of clonal heterogeneity consistent with a role in metastasis formation was observed only in 1 case in a single gene (p. Asp604Tyr of the PTPRT gene). However, when comparing variant allele frequencies (VAFs) of genomic variants in adjacent positions on chromosomes ("matched genomic variant loci") across punch samples, differences that exceeded 2 SD of the NGS assay variations (ad hoc dubbed VAF dysbalance) were observed in 7.1% of the punch samples (2.6%-12.0% per case), which indicates an intricate intermixing of mutated and nonmutated tumor cells ("intrinsic heterogeneity"). Additional OncoScan array analyses on a subset of the punch samples (31 in total) showed gross genomic aberrations as a possible explanation in only some (39.2%) of the matched genomic variant loci with VAF dysbalance. Our study provides a fairly direct (statistical model-free) view of the genomic states of microsatellite-stable colorectal carcinomas and their metastases, and suggests that Darwinian-type tumor evolution is not the key pathway of the metastasizing disease; instead, we recorded an "intrinsic" genomic heterogeneity, which may echo an initial Big Bang-like event.

A disorder clinically resembling cystic fibrosis caused by biallelic variants in the AGR2 gene.

PURPOSE: We sought to describe a disorder clinically mimicking cystic fibrosis (CF) and to elucidate its genetic cause. METHODS: Exome/genome sequencing and human phenotype ontology data of nearly 40 000 patients from our Bio/Databank were analysed. RNA sequencing of samples from the nasal mucosa from patients, carriers and controls followed by transcriptome analysis was performed. RESULTS: We identified 13 patients from 9 families with a CF-like phenotype consisting of recurrent lower respiratory infections (13/13), failure to thrive (13/13) and chronic diarrhoea (8/13), with high morbidity and mortality. All patients had biallelic variants in AGR2, (1) two splice-site variants, (2) gene deletion and (3) three missense variants. We confirmed aberrant AGR2 transcripts caused by an intronic variant and complete absence of AGR2 transcripts caused by the large gene deletion, resulting in loss of function (LoF). Furthermore, transcriptome analysis identified significant downregulation of components of the mucociliary machinery (intraciliary transport, cilium organisation), as well as upregulation of immune processes. CONCLUSION: We describe a previously unrecognised autosomal recessive disorder caused by AGR2 variants. AGR2-related disease should be considered as a differential diagnosis in patients presenting a CF-like phenotype. This has implications for the molecular diagnosis and management of these patients. AGR2 LoF is likely the disease mechanism, with consequent impairment of the mucociliary defence machinery. Future studies should aim to establish a better understanding of the disease pathophysiology and to identify potential drug targets.

Bi-allelic variants in HOPS complex subunit VPS41 cause cerebellar ataxia and abnormal membrane trafficking.

Membrane trafficking is a complex, essential process in eukaryotic cells responsible for protein transport and processing. Deficiencies in vacuolar protein sorting (VPS) proteins, key regulators of trafficking, cause abnormal intracellular segregation of macromolecules and organelles and are linked to human disease. VPS proteins function as part of complexes such as the homotypic fusion and vacuole protein sorting (HOPS) tethering complex, composed of VPS11, VPS16, VPS18, VPS33A, VPS39 and VPS41. The HOPS-specific subunit VPS41 has been reported to promote viability of dopaminergic neurons in Parkinson's disease but to date has not been linked to human disease. Here, we describe five unrelated families with nine affected individuals, all carrying homozygous variants in VPS41 that we show impact protein function. All affected individuals presented with a progressive neurodevelopmental disorder consisting of cognitive impairment, cerebellar atrophy/hypoplasia, motor dysfunction with ataxia and dystonia, and nystagmus. Zebrafish disease modelling supports the involvement of VPS41 dysfunction in the disorder, indicating lysosomal dysregulation throughout the brain and providing support for cerebellar and microglial abnormalities when vps41 was mutated. This provides the first example of human disease linked to the HOPS-specific subunit VPS41 and suggests the importance of HOPS complex activity for cerebellar function.

The Inhibitory Response to PI3K/AKT Pathway Inhibitors MK-2206 and Buparlisib Is Related to Genetic Differences in Pancreatic Ductal Adenocarcinoma Cell Lines.

The aberrant activation of the phosphoinositide 3-kinase (PI3K)/ protein kinase B (AKT) pathway is common in pancreatic ductal adenocarcinomas (PDAC). The application of inhibitors against PI3K and AKT has been considered as a therapeutic option. We investigated PDAC cell lines exposed to increasing concentrations of MK-2206 (an AKT1/2/3 inhibitor) and Buparlisib (a pan-PI3K inhibitor). Cell proliferation, metabolic activity, biomass, and apoptosis/necrosis were evaluated. Further, whole-exome sequencing (WES) and RNA sequencing (RNA-seq) were performed to analyze the recurrent aberrations and expression profiles of the inhibitor target genes and the genes frequently mutated in PDAC (Kirsten rat sarcoma virus (KRAS), Tumor protein p53 (TP53)). MK-2206 and Buparlisib demonstrated pronounced cytotoxic effects and limited cell-line-specific effects in cell death induction. WES revealed two sequence variants within the direct target genes (PIK3CA c.1143C > G in Colo357 and PIK3CD c.2480C > G in Capan-1), but a direct link to the Buparlisib response was not observed. RNA-seq demonstrated that the expression level of the inhibitor target genes did not affect the efficacy of the corresponding inhibitors. Moreover, increased resistance to MK-2206 was observed in the analyzed cell lines carrying a KRAS variant. Further, increased resistance to both inhibitors was observed in SU.86.86 carrying two TP53 missense variants. Additionally, the presence of the PIK3CA c.1143C > G in KRAS-variant-carrying cell lines was observed to correlate with increased sensitivity to Buparlisib. In conclusion, the present study reveals the distinct antitumor effects of PI3K/AKT pathway inhibitors against PDAC cell lines. Aberrations in specific target genes, as well as KRAS and TP53, individually or together, affect the efficacy of the two PI3K/AKT pathway inhibitors.

Inhibitory Response to CK II Inhibitor Silmitasertib and CDKs Inhibitor Dinaciclib Is Related to Genetic Differences in Pancreatic Ductal Adenocarcinoma Cell Lines.

Casein kinase II (CK2) and cyclin-dependent kinases (CDKs) frequently interact within multiple pathways in pancreatic ductal adenocarcinoma (PDAC). Application of CK2- and CDK-inhibitors have been considered as a therapeutic option, but are currently not part of routine chemotherapy regimens. We investigated ten PDAC cell lines exposed to increasing concentrations of silmitasertib and dinaciclib. Cell proliferation, metabolic activity, biomass, and apoptosis/necrosis were evaluated, and bioinformatic clustering was used to classify cell lines into sensitive groups based on their response to inhibitors. Furthermore, whole exome sequencing (WES) and RNA sequencing (RNA-Seq) was conducted to assess recurrent mutations and the expression profile of inhibitor targets and genes frequently mutated in PDAC, respectively. Dinaciclib and silmitasertib demonstrated pronounced and limited cell line specific effects in cell death induction, respectively. WES revealed no genomic variants causing changes in the primary structure of the corresponding inhibitor target proteins. RNA-Seq demonstrated that the expression of all inhibitor target genes was higher in the PDAC cell lines compared to non-neoplastic pancreatic tissue. The observed differences in PDAC cell line sensitivity to silmitasertib or dinaciclib did not depend on target gene expression or the identified gene variants. For the PDAC hotspot genes kirsten rat sarcoma virus (KRAS) and tumor protein p53 (TP53), three and eight variants were identified, respectively. In conclusion, both inhibitors demonstrated in vitro efficacy on the PDAC cell lines. However, aberrations and expression of inhibitor target genes did not appear to affect the efficacy of the corresponding inhibitors. In addition, specific aberrations in TP53 and KRAS affected the efficacy of both inhibitors.

A recurrent, homozygous EMC10 frameshift variant is associated with a syndrome of developmental delay with variable seizures and dysmorphic features.

PURPOSE: The endoplasmic reticulum membrane complex (EMC) is a highly conserved, multifunctional 10-protein complex related to membrane protein biology. In seven families, we identified 13 individuals with highly overlapping phenotypes who harbor a single identical homozygous frameshift variant in EMC10. METHODS: Using exome, genome, and Sanger sequencing, a recurrent frameshift EMC10 variant was identified in affected individuals in an international cohort of consanguineous families. Multiple families were independently identified and connected via Matchmaker Exchange and internal databases. We assessed the effect of the frameshift variant on EMC10 RNA and protein expression and evaluated EMC10 expression in normal human brain tissue using immunohistochemistry. RESULTS: A homozygous variant EMC10 c.287delG (Refseq NM_206538.3, p.Gly96Alafs*9) segregated with affected individuals in each family, who exhibited a phenotypic spectrum of intellectual disability (ID) and global developmental delay (GDD), variable seizures and variable dysmorphic features (elongated face, curly hair, cubitus valgus, and arachnodactyly). The variant arose on two founder haplotypes and results in significantly reduced EMC10 RNA expression and an unstable truncated EMC10 protein. CONCLUSION: We propose that a homozygous loss-of-function variant in EMC10 causes a novel syndromic neurodevelopmental phenotype. Remarkably, the recurrent variant is likely the result of a hypermutable site and arose on distinct founder haplotypes.

Combining exome/genome sequencing with data repository analysis reveals novel gene-disease associations for a wide range of genetic disorders.

PURPOSE: Within this study, we aimed to discover novel gene-disease associations in patients with no genetic diagnosis after exome/genome sequencing (ES/GS). METHODS: We followed two approaches: (1) a patient-centered approach, which after routine diagnostic analysis systematically interrogates variants in genes not yet associated to human diseases; and (2) a gene variant centered approach. For the latter, we focused on de novo variants in patients that presented with neurodevelopmental delay (NDD) and/or intellectual disability (ID), which are the most common reasons for genetic testing referrals. Gene-disease association was assessed using our data repository that combines ES/GS data and Human Phenotype Ontology terms from over 33,000 patients. RESULTS: We propose six novel gene-disease associations based on 38 patients with variants in the BLOC1S1, IPO8, MMP15, PLK1, RAP1GDS1, and ZNF699 genes. Furthermore, our results support causality of 31 additional candidate genes that had little published evidence and no registered OMIM phenotype (56 patients). The phenotypes included syndromic/nonsyndromic NDD/ID, oral-facial-digital syndrome, cardiomyopathies, malformation syndrome, short stature, skeletal dysplasia, and ciliary dyskinesia. CONCLUSION: Our results demonstrate the value of data repositories which combine clinical and genetic data for discovering and confirming gene-disease associations. Genetic laboratories should be encouraged to pursue such analyses for the benefit of undiagnosed patients and their families.

Biallelic ZNF335 mutations cause basal ganglia abnormality with progressive cerebral/cerebellar atrophy.

To date, less than 10 pedigrees have been reported with ZNF335 mutations since it was discovered in 2012 and little is known about ZNF335-related clinical spectrum. We describe a 12 years old male patient who is only child of nonconsanguineous Turkish parents. Trio whole genome sequencing identified previously unreported compound heterozygous variants in ZNF335, namely, c.3889T > A p.(Ser1297Thr) and c.758G > A p.(Arg253Gln) where transmitted by his father and mother, respectively. Patient' magnetic resonance imaging findings were overlapping to those observed in the previous cases with ZNF335 mutations. Here we report the oldest patient with biallelic ZNF335 mutations. We recommend screening for ZNF335 defects in patients with basal ganglia anomaly, secondary white matter abnormalities and microcephaly.

The Rostock International Parkinson's Disease (ROPAD) Study: Protocol and Initial Findings.

BACKGROUND: Genetic stratification of Parkinson's disease (PD) patients facilitates gene-tailored research studies and clinical trials. The objective of this study was to describe the design of and the initial data from the Rostock International Parkinson's Disease (ROPAD) study, an epidemiological observational study aiming to genetically characterize ~10,000 participants. METHODS: Recruitment criteria included (1) clinical diagnosis of PD, (2) relative of participant with a reportable LRRK2 variant, or (3) North African Berber or Ashkenazi Jew. DNA analysis involved up to 3 successive steps: (1) variant (LRRK2) and gene (GBA) screening, (2) panel sequencing of 68 PD-linked genes, and (3) genome sequencing. RESULTS: Initial data based on the first 1360 participants indicated that the ROPAD enrollment strategy revealed a genetic diagnostic yield of ~14% among a PD cohort from tertiary referral centers. CONCLUSIONS: The ROPAD screening protocol is feasible for high-throughput genetic characterization of PD participants and subsequent prioritization for gene-focused research efforts and clinical trials. © 2020 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

VPS26C homozygous nonsense variant in two cousins with neurodevelopmental deficits, growth failure, skeletal abnormalities, and distinctive facial features.

In this report, we describe two cousins with cognitive impairment, growth failure, skeletal abnormalities, and distinctive facial features. Genome sequencing failed to identify variants in known disease-associated genes explaining the phenotype. Extended comprehensive analysis of the two affected cousins' genomes, however, revealed that both share the homozygous nonsense variant c.178G>T (p.Glu60*) in the VPS26C gene. This gene encodes VPS26C, a member of the retriever integral membrane protein recycling pathway. The potential vital biological role of VPS26C, the nature of the variant which is predicted to result in loss-of-function, expression studies revealing significant reduction in the mutant transcript, and the co-segregation of the homozygous variant with the phenotype in two affected individuals all support that VPS26C is a novel gene associated with a previously unrecognized syndrome characterized by neurodevelopmental deficits, growth failure, skeletal abnormalities, and distinctive facial features.

Novel FANCI mutations in Fanconi anemia with VACTERL association.

Fanconi anemia (FA) is an inherited bone marrow failure syndrome caused by mutations in DNA repair genes; some of these patients may have features of the VACTERL association. Autosomal recessive mutations in FANCI are a rare cause of FA. We identified FANCI mutations by next generation sequencing in three patients in our FA cohort among several whose mutated gene was unknown. Four of the six mutations are novel and all mutations are likely deleterious to protein function. There are now 16 reported cases of FA due to FANCI of whom 7 have at least 3 features of the VACTERL association (44%). This suggests that the VACTERL association in patients with FA may be seen in patients with FANCI mutations more often than previously recognized.

Fanconi anemia is associated with a defect in the BRCA2 partner PALB2.

The Fanconi anemia and BRCA networks are considered interconnected, as BRCA2 gene defects have been discovered in individuals with Fanconi anemia subtype D1. Here we show that a defect in the BRCA2-interacting protein PALB2 is associated with Fanconi anemia in an individual with a new subtype. PALB2-deficient cells showed hypersensitivity to cross-linking agents and lacked chromatin-bound BRCA2; these defects were corrected upon ectopic expression of PALB2 or by spontaneous reversion.

Loss-of-Function FANCL Mutations Associate with Severe Fanconi Anemia Overlapping the VACTERL Association.

The diagnosis of VACTERL syndrome can be elusive, especially in the prenatal life, due to the presence of malformations that overlap those present in other genetic conditions, including the Fanconi anemia (FA). We report on three VACTERL cases within two families, where the two who arrived to be born died shortly after birth due to severe organs' malformations. The suspicion of VACTERL association was based on prenatal ultrasound assessment and postnatal features. Subsequent chromosome breakage analysis suggested the diagnosis of FA. Finally, by next-generation sequencing based on the analysis of the exome in one family and of a panel of Fanconi genes in the second one, we identified novel FANCL truncating mutations in both families. We used ectopic expression of wild-type FANCL to functionally correct the cellular FA phenotype for both mutations. Our study emphasizes that the diagnosis of FA should be considered when VACTERL association is suspected. Furthermore, we show that loss-of-function mutations in FANCL result in a severe clinical phenotype characterized by early postnatal death.

Microbial contamination and composition of oral samples subjected to clinical whole genome sequencing.

Biological material from the oral cavity is an excellent source of samples for genetic diagnostics. This is because collection is quick, easy-to-access, and non-invasive. We have set-up clinical whole genome sequence testing for patients with suspected hereditary disease. Beside the excellent quality of human DNA that can be isolated from such samples, we observed the presence of non-human DNA sequences at varying percentages. We investigated the proportion of non-human mapped reads (NHMR) sequenced from buccal swabs and saliva, the type of microbial genomes from which they were derived, and impact on molecular classification. Read sequences that did not map to the human reference genome were aligned to complete reference microbial reference sequences from the National Center for Biotechnology Information's (NCBI) RefSeq database using Kraken2. Out of 765 analyzed samples over 80% demonstrated more than 5% NHMRs. The majority of NHMRs were from bacterial genomes (average 69%, buccal swabs and 54% saliva), while the proportion of viruses was low, averaging 0.32% (buccal swabs) and 0.07% (saliva). We identified more than 30 different bacterial families of which Streptococcus mitis and Rothia mucilaginosa were the most common species. Importantly, the level of contamination did not impact the diagnostic yield.

Secondary findings in a large Pakistani cohort tested with whole genome sequencing.

Studies on genomic secondary findings (SFs) are diverse in participants' characteristics, sequencing methods, and versions of the ACMG SF list. Based on whole genome sequencing and the version 3.1 of the ACMG SF list, we studied SFs in 863 individuals from five different regions in Pakistan. We identified 24 ACMG SFs in 23 (2.7%) of 863 individuals: 18 of 24 were related to cardiovascular disease and four to cancer syndromes. In addition to ACMG SFs, we identified 16 (1.9%) participants with pathogenic and likely pathogenic variants in genes that were not related to the participants' clinical conditions but with clear medical actionability (non-ACMG SFs): 4 of 16 were related to eye diseases, two to metabolic disorders, and two to urinary system disorders. By testing a large Pakistani cohort with whole genome sequencing, we concluded that in countries such as Pakistan, the ACMG SF list could be expanded, and our non-ACMG SF list is one example.

BIRC2-BIRC3 amplification: a potentially druggable feature of a subset of head and neck cancers in patients with Fanconi anemia.

Head-and-neck squamous cell carcinomas (HNSCCs) are relatively common in patients with Fanconi anemia (FA), a hereditary chromosomal instability disorder. Standard chemo-radiation therapy is not tolerated in FA due to an overall somatic hypersensitivity to such treatment. The question is how to find a suitable alternative treatment. We used whole-exome and whole genome mRNA sequencing to identify major genomic and transcriptomic events associated with FA-HNSCC. CRISPR-engineered FA-knockout models were used to validate a number of top hits that were likely to be druggable. We identified deletion of 18q21.2 and amplification of 11q22.2 as prevailing copy-number alterations in FA HNSCCs, the latter of which was associated with strong overexpression of the cancer-related genes YAP1, BIRC2, BIRC3 (at 11q22.1-2). We then found the drug AZD5582, a known small molecule inhibitor of BIRC2-3, to selectively kill FA tumor cells that overexpressed BIRC2-3. This occurred at drug concentrations that did not affect the viability of untransformed FA cells. Our data indicate that 11q22.2 amplifications are relatively common oncogenic events in FA-HNSCCs, as holds for non FA-HNSCC. Therefore, chemotherapeutic inhibition of overexpressed BIRC2-3 may provide the basis for an approach to develop a clinically realistic treatment of FA-HNSCCs that carry 11q22.2 amplifications.

Immune dysregulation caused by homozygous mutations in CBLB.

CBL-B is an E3 ubiquitin ligase that ubiquitinates proteins downstream of immune receptors to downregulate positive signaling cascades. Distinct homozygous mutations in CBLB were identified in 3 unrelated children with early-onset autoimmunity, one of whom also had chronic urticaria. Patient T cells exhibited hyperproliferation in response to anti-CD3 cross-linking. One of the mutations, p.R496X, abolished CBL-B expression, and a second mutation, p.C464W, resulted in preserved CBL-B expression. The third mutation, p.H285L in the SH2 domain of CBL-B, was expressed at half the normal level in the patient's cells. Mice homozygous for the CBL-B p.H257L mutation, which corresponds to the patient's p.H285L mutation, had T and B cell hyperproliferation in response to antigen receptor cross-linking. CblbH257L mice had increased percentages of T regulatory cells (Tregs) that had normal in vitro suppressive function. However, T effector cells from the patient with the p.H285L mutation and CblbH257L mice were resistant to suppression by WT Tregs. Bone marrow-derived mast cells from CblbH257L mice were hyperactivated after FcεRI cross-linking, and CblbH257L mice demonstrated exaggerated IgE-mediated passive anaphylaxis. This study establishes CBL-B deficiency as a cause of immune dysregulation.

Disruption of the Fanconi anemia-BRCA pathway in cisplatin-sensitive ovarian tumors.

Ovarian tumor cells are often genomically unstable and hypersensitive to cisplatin. To understand the molecular basis for this phenotype, we examined the integrity of the Fanconi anemia-BRCA (FANC-BRCA) pathway in those cells. This pathway regulates cisplatin sensitivity and is governed by the coordinate activity of six genes associated with Fanconi anemia (FANCA, FANCC, FANCD2, FANCE, FANCF and FANCG) as well as BRCA1 and BRCA2 (FANCD1). Here we show that the FANC-BRCA pathway is disrupted in a subset of ovarian tumor lines. Mono-ubiquitination of FANCD2, a measure of the function of this pathway, and cisplatin resistance were restored by functional complementation with FANCF, a gene that is upstream in this pathway. FANCF inactivation in ovarian tumors resulted from methylation of its CpG island, and acquired cisplatin resistance correlated with demethylation of FANCF. We propose a model for ovarian tumor progression in which the initial methylation of FANCF is followed by FANCF demethylation and ultimately results in cisplatin resistance.

Patient-Derived Organotypic Epithelial Rafts Model Phenotypes in Juvenile-Onset Recurrent Respiratory Papillomatosis.

Juvenile-onset recurrent respiratory papillomatosis (JoRRP) is driven by human papillomavirus (HPV) low-risk strains and is associated with significant morbidity. While previous studies of 2D cultures have shed light on disease pathogenesis and demonstrated the utility of personalized medicine approaches, monolayer cultures lack the 3D tissue architecture and physiology of stratified, sequentially differentiated mucosal epithelium important in RRP disease pathogenesis. Herein we describe the establishment of JoRRP-derived primary cell populations that retain HPV genomes and viral gene expression in culture. These were directly compared to cells from matched adjacent non-diseased tissue, given the known RRP patient-to-patient variability. JoRRP papilloma versus control cells displayed decreased growth at subconfluency, with a switch to increased growth after reaching confluency, suggesting relative resistance to cell-cell contact and/or differentiation. The same papilloma cells grown as 3D organotypic rafts harbored hyperproliferation as compared to controls, with increased numbers of proliferating basal cells and inappropriately replicating suprabasal cells, mimicking phenotypes in the patient biopsies from which they were derived. These complementary model systems provide novel opportunities to elucidate disease mechanisms at distinct stages in JoRRP progression and to identify diagnostic, prognostic and therapeutic factors to personalize patient management and treatment.