Exome sequencing identified rare recurrent copy number variants and hereditary breast cancer susceptibility.

Copy number variants (CNVs) are a major source of genetic variation and can disrupt genes or affect gene dosage. They are known to be causal or underlie predisposition to various diseases. However, the role of CNVs in inherited breast cancer susceptibility has not been thoroughly investigated. To address this, we performed whole-exome sequencing based analysis of rare CNVs in 98 high-risk Northern Finnish breast cancer cases. After filtering, selected candidate alleles were validated and characterized with a combination of orthogonal methods, including PCR-based approaches, optical genome mapping and long-read sequencing. This revealed three recurrent alterations: a 31 kb deletion co-occurring with a retrotransposon insertion (delins) in RAD52, a 13.4 kb deletion in HSD17B14 and a 64 kb partial duplication of RAD51C. Notably, all these genes encode proteins involved in pathways previously identified as essential for breast cancer development. Variants were genotyped in geographically matched cases and controls (altogether 278 hereditary and 1983 unselected breast cancer cases, and 1229 controls). The RAD52 delins and HSD17B14 deletion both showed significant enrichment among cases with indications of hereditary disease susceptibility. RAD52 delins was identified in 7/278 cases (2.5%, P = 0.034, OR = 2.86, 95% CI = 1.10-7.45) and HSD17B14 deletion in 8/278 cases (2.9%, P = 0.014, OR = 3.28, 95% CI = 1.31-8.23), the frequency of both variants in the controls being 11/1229 (0.9%). This suggests a role for RAD52 and HSD17B14 in hereditary breast cancer susceptibility. The RAD51C duplication was very rare, identified only in 2/278 of hereditary cases and 2/1229 controls (P = 0.157, OR = 4.45, 95% CI = 0.62-31.70). The identification of recurrent CNVs in these genes, and especially the relatively high frequency of RAD52 and HSD17B14 alterations in the Finnish population, highlights the importance of studying CNVs alongside single nucleotide variants when searching for genetic factors underlying hereditary disease predisposition.

Enteroids to Study Pediatric Intestinal Drug Transport.

Intestinal maturational changes after birth affect the pharmacokinetics (PK) of drugs, having major implications for drug safety and efficacy. However, little is known about ontogeny-related PK patterns in the intestine. To explore the accuracy of human enteroid monolayers for studying drug transport in the pediatric intestine, we compared the drug transporter functionality and expression in enteroid monolayers and tissue from pediatrics and adults. Enteroid monolayers were cultured of 14 pediatric [median (range) age: 44 weeks (2 days-13 years)] and 5 adult donors, in which bidirectional drug transport experiments were performed. In parallel, we performed similar experiments with tissue explants in Ussing chamber using 11 pediatric [median (range) age: 54 weeks (15 weeks-10 years)] and 6 adult tissues. Enalaprilat, propranolol, talinolol, and rosuvastatin were used to test paracellular, transcellular, and transporter-mediated efflux by P-gp and breast cancer resistance protein (BCRP), respectively. In addition, we compared the expression patterns of ADME-related genes in pediatric and adult enteroid monolayers with tissues using RNA sequencing. Efflux transport by P-gp and BCRP was comparable between the enteroids and tissue. Efflux ratios (ERs) of talinolol and rosuvastatin by P-gp and BCRP, respectively, were higher in enteroid monolayers compared to Ussing chamber, likely caused by experimental differences in model setup and cellular layers present. Explorative statistics on the correlation with age showed trends of increasing ER with age for P-gp in enteroid monolayers; however, it was not significant. In the Ussing chamber setup, lower enalaprilat and propranolol transport was observed with age. Importantly, the RNA sequencing pathway analysis revealed that age-related variation in drug metabolism between neonates and adults was present in both enteroids and intestinal tissue. Age-related differences between 0 and 6 months old and adults were observed in tissue as well as in enteroid monolayers, although to a lesser extent. This study provides the first data for the further development of pediatric enteroids as an in vitro model to study age-related variation in drug transport. Overall, drug transport in enteroids was in line with data obtained from ex vivo tissue (using chamber) experiments. Additionally, pathway analysis showed similar PK-related differences between neonates and adults in both tissue and enteroid monolayers. Given the challenge to elucidate the effect of developmental changes in the pediatric age range in human tissue, intestinal enteroids derived from pediatric patients could provide a versatile experimental platform to study pediatric phenotypes.

Exploring the missing heritability in subjects with hearing loss, enlarged vestibular aqueducts, and a single or no pathogenic SLC26A4 variant.

Pathogenic variants in SLC26A4 have been associated with autosomal recessive hearing loss (arHL) and a unilateral or bilateral enlarged vestibular aqueduct (EVA). SLC26A4 is the second most frequently mutated gene in arHL. Despite the strong genotype-phenotype correlation, a significant part of cases remains genetically unresolved. In this study, we investigated a cohort of 28 Dutch index cases diagnosed with HL in combination with an EVA but without (M0) or with a single (M1) pathogenic variant in SLC26A4. To explore the missing heritability, we first determined the presence of the previously described EVA-associated haplotype (Caucasian EVA (CEVA)), characterized by 12 single nucleotide variants located upstream of SLC26A4. We found this haplotype and a delimited V1-CEVA haplotype to be significantly enriched in our M1 patient cohort (10/16 cases). The CEVA haplotype was also present in two M0 cases (2/12). Short- and long-read whole genome sequencing and optical genome mapping could not prioritize any of the variants present within the CEVA haplotype as the likely pathogenic defect. Short-read whole-genome sequencing of the six M1 cases without this haplotype and the two M0/CEVA cases only revealed previously overlooked or misinterpreted splice-altering SLC26A4 variants in two cases, who are now genetically explained. No deep-intronic or structural variants were identified in any of the M1 subjects. With this study, we have provided important insights that will pave the way for elucidating the missing heritability in M0 and M1 SLC26A4 cases. For pinpointing the pathogenic effect of the CEVA haplotype, additional analyses are required addressing defect(s) at the RNA, protein, or epigenetic level.

Optical genome mapping identifies a germline retrotransposon insertion in SMARCB1 in two siblings with atypical teratoid rhabdoid tumors.

In a subset of pediatric cancers, a germline cancer predisposition is highly suspected based on clinical and pathological findings, but genetic evidence is lacking, which hampers genetic counseling and predictive testing in the families involved. We describe a family with two siblings born from healthy parents who were both neonatally diagnosed with atypical teratoid rhabdoid tumor (ATRT). This rare and aggressive pediatric tumor is associated with biallelic inactivation of SMARCB1, and in 30% of the cases, a predisposing germline mutation is involved. Whereas the tumors of both siblings showed loss of expression of SMARCB1 and acquired homozygosity of the locus, whole exome and whole genome sequencing failed to identify germline or somatic SMARCB1 pathogenic mutations. We therefore hypothesized that the insertion of a pathogenic repeat-rich structure might hamper its detection, and we performed optical genome mapping (OGM) as an alternative strategy to identify structural variation in this locus. Using this approach, an insertion of ~2.8 kb within intron 2 of SMARCB1 was detected. Long-range PCR covering this region remained unsuccessful, but PacBio HiFi genome sequencing identified this insertion to be a SINE-VNTR-Alu, subfamily E (SVA-E) retrotransposon element, which was present in a mosaic state in the mother. This SVA-E insertion disrupts correct splicing of the gene, resulting in loss of a functional allele. This case demonstrates the power of OGM and long-read sequencing to identify genomic variations in high-risk cancer-predisposing genes that are refractory to detection with standard techniques, thereby completing the clinical and molecular diagnosis of such complex cases and greatly improving counseling and surveillance of the families involved. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Human enteroid monolayers as a potential alternative for Ussing chamber and Caco-2 monolayers to study passive permeability and drug efflux.

After oral administration, the intestine is the first site of drug absorption, making it a key determinant of the bioavailability of a drug, and hence drug efficacy and safety. Existing non-clinical models of the intestinal barrier in vitro often fail to mimic the barrier and absorption of the human intestine. We explore if human enteroid monolayers are a suitable tool for intestinal absorption studies compared to primary tissue (Ussing chamber) and Caco-2 cells. Bidirectional drug transport was determined in enteroid monolayers, fresh tissue (Ussing chamber methodology) and Caco-2 cells. Apparent permeability (Papp) and efflux ratios for enalaprilat (paracellular), propranolol (transcellular), talinolol (P-glycoprotein (P-gp)) and rosuvastatin (Breast cancer resistance protein (BCRP)) were determined and compared between all three methodologies and across intestinal regions. Bulk RNA sequencing was performed to compare gene expression between enteroid monolayers and primary tissue. All three models showed functional efflux transport by P-gp and BCRP with higher basolateral to apical (B-to-A) transport compared to apical-to-basolateral (A-to-B). B-to-A Papp values were similar for talinolol and rosuvastatin in tissue and enteroids. Paracellular transport of enalaprilat was lower and transcellular transport of propranolol was higher in enteroids compared to tissue. Enteroids appeared show more region- specific gene expression compared to tissue. Fresh tissue and enteroid monolayers both show active efflux by P-gp and BCRP in jejunum and ileum. Hence, the use of enteroid monolayers represents a promising and versatile experimental platform to complement current in vitro models.

Mobile element insertions in rare diseases: a comparative benchmark and reanalysis of 60,000 exome samples.

Mobile element insertions (MEIs) are a known cause of genetic disease but have been underexplored due to technical limitations of genetic testing methods. Various bioinformatic tools have been developed to identify MEIs in Next Generation Sequencing data. However, most tools have been developed specifically for genome sequencing (GS) data rather than exome sequencing (ES) data, which remains more widely used for routine diagnostic testing. In this study, we benchmarked six MEI detection tools (ERVcaller, MELT, Mobster, SCRAMble, TEMP2 and xTea) on ES data and on GS data from publicly available genomic samples (HG002, NA12878). For all the tools we evaluated sensitivity and precision of different filtering strategies. Results show that there were substantial differences in tool performance between ES and GS data. MELT performed best with ES data and its combination with SCRAMble increased substantially the detection rate of MEIs. By applying both tools to 10,890 ES samples from Solve-RD and 52,624 samples from Radboudumc we were able to diagnose 10 patients who had remained undiagnosed by conventional ES analysis until now. Our study shows that MELT and SCRAMble can be used reliably to identify clinically relevant MEIs in ES data. This may lead to an additional diagnosis for 1 in 3000 to 4000 patients in routine clinical ES.

Unravelling undiagnosed rare disease cases by HiFi long-read genome sequencing.

Solve-RD is a pan-European rare disease (RD) research program that aims to identify disease-causing genetic variants in previously undiagnosed RD families. We utilised 10-fold coverage HiFi long-read sequencing (LRS) for detecting causative structural variants (SVs), single nucleotide variants (SNVs), insertion-deletions (InDels), and short tandem repeat (STR) expansions in extensively studied RD families without clear molecular diagnoses. Our cohort includes 293 individuals from 114 genetically undiagnosed RD families selected by European Rare Disease Network (ERN) experts. Of these, 21 families were affected by so-called 'unsolvable' syndromes for which genetic causes remain unknown, and 93 families with at least one individual affected by a rare neurological, neuromuscular, or epilepsy disorder without genetic diagnosis despite extensive prior testing. Clinical interpretation and orthogonal validation of variants in known disease genes yielded thirteen novel genetic diagnoses due to de novo and rare inherited SNVs, InDels, SVs, and STR expansions. In an additional four families, we identified a candidate disease-causing SV affecting several genes including an MCF2 / FGF13 fusion and PSMA3 deletion. However, no common genetic cause was identified in any of the 'unsolvable' syndromes. Taken together, we found (likely) disease-causing genetic variants in 13.0% of previously unsolved families and additional candidate disease-causing SVs in another 4.3% of these families. In conclusion, our results demonstrate the added value of HiFi long-read genome sequencing in undiagnosed rare diseases.

Comprehensive de novo mutation discovery with HiFi long-read sequencing.

BACKGROUND: Long-read sequencing (LRS) techniques have been very successful in identifying structural variants (SVs). However, the high error rate of LRS made the detection of small variants (substitutions and short indels < 20 bp) more challenging. The introduction of PacBio HiFi sequencing makes LRS also suited for detecting small variation. Here we evaluate the ability of HiFi reads to detect de novo mutations (DNMs) of all types, which are technically challenging variant types and a major cause of sporadic, severe, early-onset disease. METHODS: We sequenced the genomes of eight parent-child trios using high coverage PacBio HiFi LRS (~ 30-fold coverage) and Illumina short-read sequencing (SRS) (~ 50-fold coverage). De novo substitutions, small indels, short tandem repeats (STRs) and SVs were called in both datasets and compared to each other to assess the accuracy of HiFi LRS. In addition, we determined the parent-of-origin of the small DNMs using phasing. RESULTS: We identified a total of 672 and 859 de novo substitutions/indels, 28 and 126 de novo STRs, and 24 and 1 de novo SVs in LRS and SRS respectively. For the small variants, there was a 92 and 85% concordance between the platforms. For the STRs and SVs, the concordance was 3.6 and 0.8%, and 4 and 100% respectively. We successfully validated 27/54 LRS-unique small variants, of which 11 (41%) were confirmed as true de novo events. For the SRS-unique small variants, we validated 42/133 DNMs and 8 (19%) were confirmed as true de novo event. Validation of 18 LRS-unique de novo STR calls confirmed none of the repeat expansions as true DNM. Confirmation of the 23 LRS-unique SVs was possible for 19 candidate SVs of which 10 (52.6%) were true de novo events. Furthermore, we were able to assign 96% of DNMs to their parental allele with LRS data, as opposed to just 20% with SRS data. CONCLUSIONS: HiFi LRS can now produce the most comprehensive variant dataset obtainable by a single technology in a single laboratory, allowing accurate calling of substitutions, indels, STRs and SVs. The accuracy even allows sensitive calling of DNMs on all variant levels, and also allows for phasing, which helps to distinguish true positive from false positive DNMs.

Generation of a patient-derived induced pluripotent cell line (SCTCi016-A) carrying a homozygous variant in RPE65.

Leber congenital amaurosis (LCA) can be caused by mutations in more than 20 different genes. One of these, RPE65, encodes a protein essential for the visual cycle that is expressed in retinal pigment epithelium cells. In this work, we describe the generation and characterization of the human iPSC line SCTCi16-A. This hiPSC line was generated from peripheral blood mononuclear cells (PBMCs) from a patient affected with LCA caused by the homozygous c.11+5G>A variant in the RPE65 gene. Reprograming was conducted using episomal vectors containing OCT3/4, SOX2, KLF4, L-MYC, and LIN28.

The Predicted Splicing Variant c.11+5G>A in RPE65 Leads to a Reduction in mRNA Expression in a Cell-Specific Manner.

Pathogenic variants in RPE65 lead to retinal diseases, causing a vision impairment. In this work, we investigated the pathomechanism behind the frequent RPE65 variant, c.11+5G>A. Previous in silico predictions classified this change as a splice variant. Our prediction using novel software's suggested a 124-nt exon elongation containing a premature stop codon. This elongation was validated using midigenes-based approaches. Similar results were observed in patient-derived induced pluripotent stem cells (iPSC) and photoreceptor precursor cells. However, the splicing defect in all cases was detected at low levels and thereby does not fully explain the recessive condition of the resulting disease. Long-read sequencing discarded other rearrangements or variants that could explain the diseases. Subsequently, a more relevant model was employed: iPSC-derived retinal pigment epithelium (RPE) cells. In patient-derived iPSC-RPE cells, the expression of RPE65 was strongly reduced even after inhibiting a nonsense-mediated decay, contradicting the predicted splicing defect. Additional experiments demonstrated a cell-specific gene expression reduction due to the presence of the c.11+5G>A variant. This decrease also leads to the lack of the RPE65 protein, and differences in size and pigmentation between the patient and control iPSC-RPE. Altogether, our data suggest that the c.11+5G>A variant causes a cell-specific defect in the expression of RPE65 rather than the anticipated splicing defect which was predicted in silico.

Next Generation Exome Sequencing of Pediatric Asthma Identifies Rare and Novel Variants in Candidate Genes.

Multiple genes have been implicated to have a role in asthma predisposition by association studies. Pediatric patients often manifest a more extensive form of this disease and a particularly severe disease course. It is likely that genetic predisposition could play a more substantial role in this group. This study is aimed at identifying the spectrum of rare and novel variation in known pediatric asthma susceptibility genes using whole exome sequencing analysis in nine individual cases of childhood onset allergic asthma. DNA samples from the nine children with a history of bronchial asthma diagnosis underwent whole exome sequencing on Ion Proton. For each patient, the entire complement of rare variation within strongly associated candidate genes was catalogued. The analysis showed 21 variants in the subjects, 13 had been previously identified, and 8 were novel. Also, among of which, nineteen were nonsynonymous and 2 were nonsense. With regard to the novel variants, the 2 nonsynonymous variants in the PRKG1 gene (PRKG1: p.C519W and PRKG1: p.G520W) were presented in 4 cases, and a nonsynonymous variant in the MAVS gene (MAVS: p.A45V) was identified in 3 cases. The variants we found in this study will enrich the variant spectrum and build up the database in the Saudi population. Novel eight variants were identified in the study which provides more evidence in the genetic susceptibility in asthma among Saudi children, providing a genetic screening map for the molecular genetic determinants of allergic disease in Saudi children, with the goal of reducing the impact of chronic diseases on the health and the economy. We believe that the advanced specified statistical filtration/annotation programs used in this study succeeded to release such results in a preliminary study, exploring the genetic map of that disease in Saudi children.

Long-read technologies identify a hidden inverted duplication in a family with choroideremia.

The lack of molecular diagnoses in rare genetic diseases can be explained by limitations of current standard genomic technologies. Upcoming long-read techniques have complementary strengths to overcome these limitations, with a particular strength in identifying structural variants. By using optical genome mapping and long-read sequencing, we aimed to identify the pathogenic variant in a large family with X-linked choroideremia. In this family, aberrant splicing of exon 12 of the choroideremia gene CHM was detected in 2003, but the underlying genomic defect remained elusive. Optical genome mapping and long-read sequencing approaches now revealed an intragenic 1,752 bp inverted duplication including exon 12 and surrounding regions, located downstream of the wild-type copy of exon 12. Both breakpoint junctions were confirmed with Sanger sequencing and segregate with the X-linked inheritance in the family. The breakpoint junctions displayed sequence microhomology suggestive for an erroneous replication mechanism as the origin of the structural variant. The inverted duplication is predicted to result in a hairpin formation of the pre-mRNA with the wild-type exon 12, leading to exon skipping in the mature mRNA. The identified inverted duplication is deemed the hidden pathogenic cause of disease in this family. Our study shows that optical genome mapping and long-read sequencing have significant potential for the identification of (hidden) structural variants in rare genetic diseases.

Comprehensive sequencing of the myocilin gene in a selected cohort of severe primary open-angle glaucoma patients.

Primary open-angle glaucoma (POAG) is the most common form of glaucoma, prevalent in approximately 1-2% of Caucasians in the UK over the age of 40. It is characterised by an open anterior chamber angle, raised intraocular pressure (IOP) and optic nerve damage leading to loss of sight. The myocilin gene (MYOC) is the most common glaucoma-causing gene, accounting for ~2% of British POAG cases. 358 patients were selected for next generation sequencing (NGS) with the following selection criteria: Caucasian ethnicity, intraocular pressure (IOP) 21-40 mm Hg, cup:disc ratio ≥0.6 and visual field mean deviation ≤-3. The entire MYOC gene (17,321 bp) was captured including the promoter, introns, UTRs and coding exons. We identify 12 exonic variants (one stop-gain, five missense and six synonymous variants), two promoter variants, 133 intronic variants, two 3' UTR variants and 23 intergenic variants. Four known or predicted pathogenic exonic variants (p.R126W, p.K216K, p.Q368* and p.T419A) were identified across 11 patients, which accounts for 3.07% of this POAG cohort. This is the first time that the entire region of MYOC has been sequenced and variants reported for a cohort of POAG patients.

A small gene sequencing panel realises a high diagnostic rate in patients with congenital nystagmus following basic phenotyping.

Nystagmus is a disorder of uncontrolled eye movement and can occur as an isolated trait (idiopathic INS, IINS) or as part of multisystem disorders such as albinism, significant visual disorders or neurological disease. Eighty-one unrelated patients with nystagmus underwent routine ocular phenotyping using commonly available phenotyping methods and were grouped into four sub-cohorts according to the level of phenotyping information gained and their findings. DNA was extracted and sequenced using a broad utility next generation sequencing (NGS) gene panel. A clinical subpanel of genes for nystagmus/albinism was utilised and likely causal variants were prioritised according to methods currently employed by clinical diagnostic laboratories. We determine the likely underlying genetic cause for 43.2% of participants with similar yields regardless of prior phenotyping. This study demonstrates that a diagnostic workflow combining basic ocular phenotyping and a clinically available targeted NGS panel, can provide a high diagnostic yield for patients with infantile nystagmus, enabling access to disease specific management at a young age and reducing the need for multiple costly, often invasive tests. By describing diagnostic yield for groups of patients with incomplete phenotyping data, it also permits the subsequent design of 'real-world' diagnostic workflows and illustrates the changing role of genetic testing in modern diagnostic workflows for heterogeneous ophthalmic disorders.

Identification of a functionally significant tri-allelic genotype in the Tyrosinase gene (TYR) causing hypomorphic oculocutaneous albinism (OCA1B).

Oculocutaneous albinism (OCA) and ocular albinism (OA) are inherited disorders of melanin biosynthesis, resulting in loss of pigment and severe visual deficits. OCA encompasses a range of subtypes with overlapping, often hypomorphic phenotypes. OCA1 is the most common cause of albinism in European populations and is inherited through autosomal recessive mutations in the Tyrosinase (TYR) gene. However, there is a high level of reported missing heritability, where only a single heterozygous mutation is found in TYR. This is also the case for other OCA subtypes including OCA2 caused by mutations in the OCA2 gene. Here we have interrogated the genetic cause of albinism in a well phenotyped, hypomorphic albinism population by sequencing a broad gene panel and performing segregation studies on phenotyped family members. Of eighteen probands we can confidently diagnose three with OA and OCA2, and one with a PAX6 mutation. Of six probands with only a single heterozygous mutation in TYR, all were found to have the two common variants S192Y and R402Q. Our results suggest that a combination of R402Q and S192Y with a deleterious mutation in a 'tri-allelic genotype' can account for missing heritability in some hypomorphic OCA1 albinism phenotypes.