Building a precision medicine infrastructure at a national level: The Swedish experience.

Precision medicine has the potential to transform healthcare by moving from one-size-fits-all to personalised treatment and care. This transition has been greatly facilitated through new high-throughput sequencing technologies that can provide the unique molecular profile of each individual patient, along with the rapid development of targeted therapies directed to the Achilles heels of each disease. To implement precision medicine approaches in healthcare, many countries have adopted national strategies and initiated genomic/precision medicine initiatives to provide equal access to all citizens. In other countries, such as Sweden, this has proven more difficult due to regionally organised healthcare. Using a bottom-up approach, key stakeholders from academia, healthcare, industry and patient organisations joined forces and formed Genomic Medicine Sweden (GMS), a national infrastructure for the implementation of precision medicine across the country. To achieve this, Genomic Medicine Centres have been established to provide regionally distributed genomic services, and a national informatics infrastructure has been built to allow secure data handling and sharing. GMS has a broad scope focusing on rare diseases, cancer, pharmacogenomics, infectious diseases and complex diseases, while also providing expertise in informatics, ethical and legal issues, health economy, industry collaboration and education. In this review, we summarise our experience in building a national infrastructure for precision medicine. We also provide key examples how precision medicine already has been successfully implemented within our focus areas. Finally, we bring up challenges and opportunities associated with precision medicine implementation, the importance of international collaboration, as well as the future perspective in the field of precision medicine.

The value of genome-wide analysis in craniosynostosis.

Background: This study assessed the diagnostic yield of high-throughput sequencing methods in a cohort of craniosynostosis (CS) patients not presenting causal variants identified through previous targeted analysis. Methods: Whole-genome or whole-exome sequencing (WGS/WES) was performed in a cohort of 59 patients (from 57 families) assessed by retrospective phenotyping as having syndromic or nonsyndromic CS. Results: A syndromic form was identified in 51% of the unrelated cases. A genetic cause was identified in 38% of syndromic cases, with novel variants detected in FGFR2 (a rare Alu insertion), TWIST1, TCF12, KIAA0586, HDAC9, FOXP1, and NSD2. Additionally, we report two patients with rare recurrent variants in KAT6A and YY1 as well as two patients with structural genomic aberrations: one with a 22q13 duplication and one with a complex rearrangement involving chromosome 2 (2p25 duplication including SOX11 and deletion of 2q22). Moreover, we identified potentially relevant variants in 87% of the remaining families with no previously detected causal variants, including novel variants in ADAMTSL4, ASH1L, ATRX, C2CD3, CHD5, ERF, H4C5, IFT122, IFT140, KDM6B, KMT2D, LTBP1, MAP3K7, NOTCH2, NSD1, SOS1, SPRY1, POLR2A, PRRX1, RECQL4, TAB2, TAOK1, TET3, TGFBR1, TCF20, and ZBTB20. Conclusion: These results confirm WGS/WES as a powerful diagnostic tool capable of either targeted in silico or broad genomic analysis depending on phenotypic presentation (e.g., classical or unusual forms of syndromic CS).

The outcome of targeted NGS screening in patients with syndromic forms of sagittal and pansynostosis - IL11RA is an emerging core-gene for pansynostosis.

Here, we have studied the prevalence and spectrum of genetic alterations in syndromic forms of sagittal and pansynostosis. Eighteen patients with sagittal synostosis (isolated or combined with other synostoses, except coronal) or pansynostosis were phenotypically assessed by retrospective analysis of medical records, three-dimensional computed tomography skull reconstructions, and registered photos. Patient DNAs were analyzed using a targeted next-generation sequencing (NGS) panel including 63 craniosynostosis (CS) related genes. Pathogenic and likely pathogenic variants were found in 72% of the cases, mainly affecting FGFR2, TWIST1, IL11RA, and SKI. Two patients that were negative at NGS screening - one with a supernumerary marker chromosome with duplication of 15q25.2q26.3 and one with a pathogenic PHEX variant - were identified using microarray and single gene analysis, respectively. The overall diagnostic rate in the cohort was thus 83%. We identified two novel likely pathogenic variants in FGFR2 (NM_022970.3: c.811_812delGGinsCC, p.Gly271Pro) and TWIST1 (NM_000474.3: c.476T > A, p.Leu159His), and a novel variant of unclear phenotypic significance in RUNX2 (NM_001024630.3: c.340G > A, p.Val114Ile) which could suggest a modulatory effect. Notably, we also identified three new patients with pansynostosis and a Crouzon-like phenotype with IL11RA mutation. Targeted NGS using a broad panel of CS-related genes is a simple and powerful tool for detecting pathogenic mutations in patients with syndromic forms of CS and multiple suture involvement, in particular pansynostosis. Our results provide additional evidence of an association between pansynostosis and IL11RA, an emerging core gene for autosomal recessive CS.

[The utility of whole genome sequencing in rare disease diagnostics]. / Helgenomanalys vid sällsynta diagnoser ger stor patientnytta.

If a disease affects fewer than 1 in 2 000, the European Union defines it as a rare disease. Globally, about 300 million people live with a rare disease, and in Sweden about 400 000. There are approximately 7 000 different rare diseases. The clinical picture varies from a single symptom to complex patterns with multiple organs affected, often combined with cognitive and motor impairment. At least 72 % of all rare diseases are genetic and 70% have childhood onset. Many patients are undiagnosed and do not receive optimal treatment. Today, only 5% of rare diseases have an approved treatment option. With modern genetic high throughput techniques, many disease-causing mutations are identified, increasing the possibility of personalized treatment and prevention strategies, designed by the individual's genetic conditions, i.e. precision medicine.

Cytogenetically visible inversions are formed by multiple molecular mechanisms.

Cytogenetically detected inversions are generally assumed to be copy number and phenotypically neutral events. While nonallelic homologous recombination is thought to play a major role, recent data suggest the involvement of other molecular mechanisms in inversion formation. Using a combination of short-read whole-genome sequencing (WGS), 10X Genomics Chromium WGS, droplet digital polymerase chain reaction and array comparative genomic hybridization we investigated the genomic structure of 18 large unique cytogenetically detected chromosomal inversions and achieved nucleotide resolution of at least one chromosomal inversion junction for 13/18 (72%). Surprisingly, we observed that seemingly copy number neutral inversions can be accompanied by a copy-number gain of up to 350 kb and local genomic complexities (3/18, 17%). In the resolved inversions, the mutational signatures are consistent with nonhomologous end-joining (8/13, 62%) or microhomology-mediated break-induced replication (5/13, 38%). Our study indicates that short-read 30x coverage WGS can detect a substantial fraction of chromosomal inversions. Moreover, replication-based mechanisms are responsible for approximately 38% of those events leading to a significant proportion of inversions that are actually accompanied by additional copy-number variation potentially contributing to the overall phenotypic presentation of those patients.

NGS targeted screening of 100 Scandinavian patients with coronal synostosis.

Craniosynostosis (CS), the premature closure of one or more cranial sutures, occurs both as part of a syndrome or in isolation (nonsyndromic form). Here, we have studied the prevalence and spectrum of genetic alterations associated with coronal suture closure in 100 Scandinavian patients treated at a single craniofacial unit. All patients were phenotypically assessed and analyzed with a custom-designed 63 gene NGS-panel. Most cases (78%) were syndromic forms of CS. Pathogenic and likely pathogenic variants explaining the phenotype were found in 80% of the families with syndromic CS and in 14% of those with nonsyndromic CS. Sixty-five percent of the families had mutations in the CS core genes FGFR2, TWIST1, FGFR3, TCF12, EFNB1, FGFR1, and POR. Five novel pathogenic/likely pathogenic variants in TWIST1, TCF12, and EFNB1 were identified. We also found novel variants in SPECC1L, IGF1R, and CYP26B1 with a possible modulator phenotypic effect. Our findings demonstrate that NGS targeted sequencing is a powerful tool to detect pathogenic mutations in patients with coronal CS and further emphasize the importance of thorough assessment of the patient's phenotype for reliable interpretation of the molecular findings. This is particularly important in patients with complex phenotypes and rare forms of CS.

Replicative and non-replicative mechanisms in the formation of clustered CNVs are indicated by whole genome characterization.

Clustered copy number variants (CNVs) as detected by chromosomal microarray analysis (CMA) are often reported as germline chromothripsis. However, such cases might need further investigations by massive parallel whole genome sequencing (WGS) in order to accurately define the underlying complex rearrangement, predict the occurrence mechanisms and identify additional complexities. Here, we utilized WGS to delineate the rearrangement structure of 21 clustered CNV carriers first investigated by CMA and identified a total of 83 breakpoint junctions (BPJs). The rearrangements were further sub-classified depending on the patterns observed: I) Cases with only deletions (n = 8) often had additional structural rearrangements, such as insertions and inversions typical to chromothripsis; II) cases with only duplications (n = 7) or III) combinations of deletions and duplications (n = 6) demonstrated mostly interspersed duplications and BPJs enriched with microhomology. In two cases the rearrangement mutational signatures indicated both a breakage-fusion-bridge cycle process and haltered formation of a ring chromosome. Finally, we observed two cases with Alu- and LINE-mediated rearrangements as well as two unrelated individuals with seemingly identical clustered CNVs on 2p25.3, possibly a rare European founder rearrangement. In conclusion, through detailed characterization of the derivative chromosomes we show that multiple mechanisms are likely involved in the formation of clustered CNVs and add further evidence for chromoanagenesis mechanisms in both "simple" and highly complex chromosomal rearrangements. Finally, WGS characterization adds positional information, important for a correct clinical interpretation and deciphering mechanisms involved in the formation of these rearrangements.

Hyperestrogenism Affects Adult Height Outcome in Growth Hormone Treated Boys With Silver-Russell Syndrome.

Background: Intrauterine growth retardation and short stature are common features in Silver-Russell syndrome (SRS). Despite recombinant growth hormone (rGH) treatment, poor pubertal height gain, affecting adult height (AH), is common. This study investigated whether growth patterns and estrogen concentrations are associated with AH outcome in rGH treated SRS males. Methods: In this retrospective longitudinal single-center study, 11 males with SRS were classified as non-responders (NR = 6) or responders (R = 5), depending on AH adjusted for midparental height. Epigenetic analysis and longitudinal growth measures, including bone age, rGH related parameters, pubertal development, gonadotropins and estrogen concentrations, were analyzed until AH. Results: Pubarche before 9 years was only observed in one NR. At 10 years of age, there was no difference in gonadotropins between NR and R. However, estradiol (E2) concentrations at 10 years of age showed a strong association to AH adjusted for MPH (r = -0.78, p < 0.001). Serum E2 (pmol/L) was significantly higher in NR at ages 10 years [median (range) 2 (<2-5) vs. <2 (<2)], 12 years [23 (10-57) vs. 2 (<2-2)] and 14 years [77 (54-87) vs. 24 (<2-38)] but not at 16 years. Birth weight standard deviation score (SDS) was lower in NR [-4.1 (-4.7 to -2.1) vs. -2.7 (-3.3 to -1.7)]. Weight gain (SDS) until pubertal onset was greater in NR [2.4 (1.4-3.5) vs. 0.8 (-0.4 to 1.7)] and pubertal height gain (SDS) was lower in NR [-1.0 (-2.7-0.4) vs. 0.1 (-0.1 to 1.1)]. At AH, a number of NR and R had high E2 concentrations and small testes. Conclusion: Increased E2 concentrations at age 10, 12, and 14 years were associated to less pubertal height gain, thus affecting AH. Due to the small number of patients, the results need to be confirmed in larger cohorts. The finding of impaired testicular development stresses the need of hormonal evaluation as a complement to clinical and radiological assessment when predicting AH in males with SRS.

On the significance of craniosynostosis in a case of Kabuki syndrome with a concomitant KMT2D mutation and 3.2 Mbp de novo 10q22.3q23.1 deletion.

Craniosynostosis has rarely been described in patients with Kabuki syndrome. We report here a boy with facial asymmetry due to combined premature synostosis of the right coronal and sagittal sutures as well as several symptoms reminiscent of Kabuki syndrome (KS). Our case supports previous observations and suggests that craniosynostosis is a part of the KS phenotype. The uniqueness of our case is the sporadic co-occurrence of two genetic disorders, that is, a de novo frameshift variant in the KMT2D gene and a de novo 3.2 Mbp 10q22.3q23.1 deletion. Our findings emphasize the importance of the initial clinical assessment of children with craniosynostosis and that genomic and monogenic disorders, such as Kabuki syndrome, should be considered among the differential diagnoses of syndromic forms of craniosynostosis.

Exon 2 duplication of the MID1 gene in a patient with a mild phenotype of Opitz G/BBB syndrome.

The X-linked form of Opitz G/BBB syndrome is a congenital midline malformation syndrome caused by MID1 loss-of-function mutations, including point mutations and small-sized duplications, insertions, and deletions. Three patients with an Opitz G/BBB syndrome phenotype and relatively large duplications of part of the MID1 gene have been described up to date. Here we report a 2-months-old boy with a very mild phenotype including craniofacial dysmorphism, swallowing difficulties, and a normal psychomotor development. Molecular karyotyping revealed a 57-kb duplication involving exon 2 of the MID1 gene. The in-frame tandem duplication was confirmed by MID1 transcript analysis. This alteration results likely in a mutant MID1 protein which contains 32 duplicated amino acids in the first part of the coiled-coil domain. The mild phenotype of the patient with the microduplication suggests that MID1 mutations can be found in patients with hypertelorism with or without other clinical signs and MID1 alterations might be missed in individuals not fulfilling the minimal criteria for diagnosis of X-linked Opitz G/BBB syndrome. This report further emphasizes the genotype-first approach in medical genetics in general and patients with unspecific clinical features in particular.

Estrogen receptor alpha gene polymorphism and endometrial cancer risk--a case-control study.

BACKGROUND: Estrogen is an established endometrial carcinogen. One of the most important mediators of estrogenic action is the estrogen receptor alpha. We have investigated whether polymorphic variation in the estrogen receptor alpha gene (ESR1) is associated with endometrial cancer risk. METHODS: In 702 cases with invasive endometrial cancer and 1563 controls, we genotyped five markers in ESR1 and used logistic regression models to estimate odds ratios (OR) and 95 percent confidence intervals (CI). RESULTS: We found an association between rs2234670, rs2234693, as well as rs9340799, markers in strong linkage disequilibrium (LD), and endometrial cancer risk. The association with rs9340799 was the strongest, OR 0.75 (CI 0.60-0.93) for heterozygous and OR 0.53 (CI 0.37-0.77) for homozygous rare compared to those homozygous for the most common allele. Haplotype models did not fit better to the data than single marker models. CONCLUSION: We found that intronic variation in ESR1 was associated with endometrial cancer risk.

Multiple displacement amplification to create a long-lasting source of DNA for genetic studies.

In many situations there may not be sufficient DNA collected from patient or population cohorts to meet the requirements of genome-wide analysis of SNPs, genomic copy number polymorphisms, or acquired copy number alternations. When the amount of available DNA for genotype analysis is limited, high performance whole-genome amplification (WGA) represents a new development in genetic analysis. It is especially useful for analysis of DNA extracted from stored histology slides, tissue samples, buccal swabs, or blood stains collected on filter paper. The multiple displacement amplification (MDA) method, which relies on isothermal amplification using the DNA polymerase of the bacteriophage phi29, is a recently developed technique for high performance WGA. This review addresses new trends in the technical performance of MDA and its applications to genetic analyses. The main challenge of WGA methods is to obtain balanced and faithful replication of all chromosomal regions without the loss of or preferential amplification of any genomic loci or allele. In multiple comparisons to other WGA methods, MDA appears to be most reliable for genotyping, with the most favorable call rates, best genomic coverage, and lowest amplification bias.

Genotyping single-nucleotide polymorphisms by minisequencing using tag arrays.

The need for large-scale and high-throughput methods for SNP genotyping has rapidly increased during the last decade. Our system, presented here, combines the highly specific genotyping principle of minisequencing with the advantages of a microarray format that allows highly multiplexed and parallel analysis. Cyclic minisequencing reactions with fluorescently labeled dideoxynucleotides (ddNTPs) are performed in solution using multiplex PCR product as template and detection primers, designed to anneal immediately adjacent and upstream of the SNP site. The detection primers carry unique 5' tag sequences and oligonucleotides complementary to the tag sequence, cTags, are immobilized on a microarray. After extension, the tagged detection primers are allowed to hybridize to the cTags; then the fluorescent signals from the array are measured, and the genotypes are deduced according to the label incorporated. The "array of arrays" format of the system, accomplished by a silicon rubber grid giving separate reaction chambers, allows either 80 or 14 samples to be analyzed for up to 200 or 600 SNPs, respectively, on a single microscope slide.

Silhouette scores for assessment of SNP genotype clusters.

BACKGROUND: High-throughput genotyping of single nucleotide polymorphisms (SNPs) generates large amounts of data. In many SNP genotyping assays, the genotype assignment is based on scatter plots of signals corresponding to the two SNP alleles. In a robust assay the three clusters that define the genotypes are well separated and the distances between the data points within a cluster are short. "Silhouettes" is a graphical aid for interpretation and validation of data clusters that provides a measure of how well a data point was classified when it was assigned to a cluster. Thus "Silhouettes" can potentially be used as a quality measure for SNP genotyping results and for objective comparison of the performance of SNP assays at different circumstances. RESULTS: We created a program (ClusterA) for calculating "Silhouette scores", and applied it to assess the quality of SNP genotype clusters obtained by single nucleotide primer extension ("minisequencing") in the Tag-microarray format. A Silhouette score condenses the quality of the genotype assignment for each SNP assay into a single numeric value, which ranges from 1.0, when the genotype assignment is unequivocal, down to -1.0, when the genotype assignment has been arbitrary. In the present study we applied Silhouette scores to compare the performance of four DNA polymerases in our minisequencing system by analyzing 26 SNPs in both DNA polarities in 16 DNA samples. We found Silhouettes to provide a relevant measure for the quality of SNP assays at different reaction conditions, illustrated by the four DNA polymerases here. According to our result, the genotypes can be unequivocally assigned without manual inspection when the Silhouette score for a SNP assay is > 0.65. All four DNA polymerases performed satisfactorily in our Tag-array minisequencing system. CONCLUSION: "Silhouette scores" for assessing the quality of SNP genotyping clusters is convenient for evaluating the quality of SNP genotype assignment, and provides an objective, numeric measure for comparing the performance of SNP assays. The program we created for calculating Silhouette scores is freely available, and can be used for quality assessment of the results from all genotyping systems, where the genotypes are assigned by cluster analysis using scatter plots.

Oestrogen receptor alpha gene haplotype and postmenopausal breast cancer risk: a case control study.

INTRODUCTION: Oestrogen receptor alpha, which mediates the effect of oestrogen in target tissues, is genetically polymorphic. Because breast cancer development is dependent on oestrogenic influence, we have investigated whether polymorphisms in the oestrogen receptor alpha gene (ESR1) are associated with breast cancer risk. METHODS: We genotyped breast cancer cases and age-matched population controls for one microsatellite marker and four single-nucleotide polymorphisms (SNPs) in ESR1. The numbers of genotyped cases and controls for each marker were as follows: TAn, 1514 cases and 1514 controls; c.454-397C --> T, 1557 cases and 1512 controls; c.454-351A --> G, 1556 cases and 1512 controls; c.729C --> T, 1562 cases and 1513 controls; c.975C --> G, 1562 cases and 1513 controls. Using logistic regression models, we calculated odds ratios (ORs) and 95% confidence intervals (CIs). Haplotype effects were estimated in an exploratory analysis, using expectation-maximisation algorithms for case-control study data. RESULTS: There were no compelling associations between single polymorphic loci and breast cancer risk. In haplotype analyses, a common haplotype of the c.454-351A --> G or c.454-397C --> T and c.975C --> G SNPs appeared to be associated with an increased risk for ductal breast cancer: one copy of the c.454-351A --> G and c.975C --> G haplotype entailed an OR of 1.19 (95% CI 1.06-1.33) and two copies with an OR of 1.42 (95% CI 1.15-1.77), compared with no copies, under a model of multiplicative penetrance. The association with the c.454-397C --> T and c.975C --> G haplotypes was similar. Our data indicated that these haplotypes were more influential in women with a high body mass index. Adjustment for multiple comparisons rendered the associations statistically non-significant. CONCLUSION: We found suggestions of an association between common haplotypes in ESR1 and the risk for ductal breast cancer that is stronger in heavy women.

Microarrays for genotyping human group a rotavirus by multiplex capture and type-specific primer extension.

Human group A rotavirus (HRV) is the major cause of severe gastroenteritis in infants worldwide. HRV shares the feature of a high degree of genetic diversity with many other RNA viruses, and therefore, genotyping of this organism is more complicated than genotyping of more stable DNA viruses. We describe a novel microarray-based method that allows high-throughput genotyping of RNA viruses with a high degree of polymorphism by multiplex capture and type-specific extension on microarrays. Denatured reverse transcription (RT)-PCR products derived from two outer capsid genes of clinical isolates of HRV were hybridized to immobilized capture oligonucleotides representing the most commonly occurring P and G genotypes on a microarray. Specific primer extension of the type-specific capture oligonucleotides was applied to incorporate the fluorescent nucleotide analogue cyanine 5-labeled dUTP as a detectable label. Laser scanning and fluorescence detection of the microarrays was followed by visual or computer-assisted interpretation of the fluorescence patterns generated on the microarrays. Initially, the method detected HRV in all 40 samples and correctly determined both the G and the P genotypes of 35 of the 40 strains analyzed. After modification by inclusion of additional capture oligonucleotides specific for the initially unassigned genotypes, all genotypes could be correctly defined. The results of genotyping with the microarray fully agreed with the results obtained by nucleotide sequence analysis and sequence-specific multiplex RT-PCR. Owing to its robustness, simplicity, and general utility, the microarray-based method may gain wide applicability for the genotyping of microorganisms, including highly variable RNA and DNA viruses.

Quantitative evaluation by minisequencing and microarrays reveals accurate multiplexed SNP genotyping of whole genome amplified DNA.

Whole genome amplification (WGA) procedures such as primer extension preamplification (PEP) or multiple displacement amplification (MDA) have the potential to provide an unlimited source of DNA for large-scale genetic studies. We have performed a quantitative evaluation of PEP and MDA for genotyping single nucleotide polymorphisms (SNPs) using multiplex, four-color fluorescent minisequencing in a microarray format. Forty-five SNPs were genotyped and the WGA methods were evaluated with respect to genotyping success, signal-to-noise ratios, power of genotype discrimination, yield and imbalanced amplification of alleles in the MDA product. Both PEP and MDA products provided genotyping results with a high concordance to genomic DNA. For PEP products the power of genotype discrimination was lower than for MDA due to a 2-fold lower signal-to-noise ratio. MDA products were indistinguishable from genomic DNA in all aspects studied. To obtain faithful representation of the SNP alleles at least 0.3 ng DNA should be used per MDA reaction. We conclude that the use of WGA, and MDA in particular, is a highly promising procedure for producing DNA in sufficient amounts even for genome wide SNP mapping studies.