Epilepsies of presumed genetic etiology show enrichment of rare variants that occur in the general population.

Previous studies suggested that severe epilepsies, e.g., developmental and epileptic encephalopathies (DEEs), are mainly caused by ultra-rare de novo genetic variants. For milder disease, rare genetic variants could contribute to the phenotype. To determine the importance of rare variants for different epilepsy types, we analyzed a whole-exome sequencing cohort of 9,170 epilepsy-affected individuals and 8,436 control individuals. Here, we separately analyzed three different groups of epilepsies: severe DEEs, genetic generalized epilepsy (GGE), and non-acquired focal epilepsy (NAFE). We required qualifying rare variants (QRVs) to occur in control individuals with an allele count ≥ 1 and a minor allele frequency ≤ 1:1,000, to be predicted as deleterious (CADD ≥ 20), and to have an odds ratio in individuals with epilepsy ≥ 2. We identified genes enriched with QRVs primarily in NAFE (n = 72), followed by GGE (n = 32) and DEE (n = 21). This suggests that rare variants may play a more important role for causality of NAFE than for DEE. Moreover, we found that genes harboring QRVs, e.g., HSGP2, FLNA, or TNC, encode proteins that are involved in structuring the brain extracellular matrix. The present study confirms an involvement of rare variants for NAFE that occur also in the general population, while in DEE and GGE, the contribution of such variants appears more limited.

Diagnostic Yield of Targeted Hearing Loss Gene Panel Sequencing in a Large German Cohort With a Balanced Age Distribution from a Single Diagnostic Center: An Eight-year Study.

OBJECTIVES: Hereditary hearing loss exhibits high degrees of genetic and clinical heterogeneity. To elucidate the population-specific and age-related genetic and clinical spectra of hereditary hearing loss, we investigated the sequencing data of causally associated hearing loss genes in a large cohort of hearing-impaired probands with a balanced age distribution from a single center in Southwest Germany. DESIGN: Genetic testing was applied to 305 hearing-impaired probands/families with a suspected genetic hearing loss etiology and a balanced age distribution over a period of 8 years (2011-2018). These individuals were representative of the regional population according to age and sex distributions. The genetic testing workflow consisted of single-gene screening (n = 21) and custom-designed hearing loss gene panel sequencing (n = 284) targeting known nonsyndromic and syndromic hearing loss genes in a diagnostic setup. Retrospective reanalysis of sequencing data was conducted by applying the current American College of Medical Genetics and Genomics/Association for Molecular Pathology guidelines. RESULTS: A genetic diagnosis was established for 75 (25%) of the probands that involved 75 causal variants in 35 genes, including 16 novel causal variants and 9 medically significant variant reclassifications. Nearly half of the solved cases (47%; n = 35) were related to variants in the five most frequently affected genes: GJB2 (25%), MYO15A, WFS1, SLC26A4, and COL11A1 (all 5%). Nearly one-quarter of the cases (23%; n = 17) were associated with variants in seven additional genes (TMPRSS3, COL4A3, LOXHD1, EDNRB, MYO6, TECTA, and USH2A). The remaining one-third of single cases (33%; n = 25) were linked to variants in 25 distinct genes. Diagnostic rates and gene distribution were highly dependent on phenotypic characteristics. A positive family history of autosomal-recessive inheritance in combination with early onset and higher grades of hearing loss significantly increased the solve rate up to 60%, while late onset and lower grades of hearing loss yielded significantly fewer diagnoses. Regarding genetic diagnoses, autosomal-dominant genes accounted for 37%, autosomal-recessive genes for 60%, and X-linked genes for 3% of the solved cases. Syndromic/nonsyndromic hearing loss mimic genes were affected in 27% of the genetic diagnoses. CONCLUSIONS: The genetic epidemiology of the largest German cohort subjected to comprehensive targeted sequencing for hereditary hearing loss to date revealed broad causal gene and variant spectra in this population. Targeted hearing loss gene panel analysis proved to be an effective tool for ensuring an appropriate diagnostic yield in a routine clinical setting including the identification of novel variants and medically significant reclassifications. Solve rates were highly sensitive to phenotypic characteristics. The unique population-adapted and balanced age distribution of the cohort favoring late hearing loss onset uncovered a markedly large contribution of autosomal-dominant genes to the diagnoses which may be a representative for other age balanced cohorts in other populations.

Trio exome sequencing is highly relevant in prenatal diagnostics.

OBJECTIVE: About 3% of newborns show malformations, with about 20% of the affected having genetic causes. Clarification of genetic diseases in postnatal diagnostics was significantly improved with high-throughput sequencing, in particular through whole exome sequencing covering all protein-coding regions. Here, we aim to extend the use of this technology to prenatal diagnostics. METHOD: Between 07/2018 and 10/2020, 500 pregnancies with fetal ultrasound abnormalities were analyzed after genetic counseling as part of prenatal diagnostics using WES of the fetus and parents. RESULTS: Molecular genetic findings could explain ultrasound abnormalities in 38% of affected fetuses. In 47% of these, disease-causing de novo variants were found. Pathogenic variants in genes with autosomal recessive or X-linked inheritance were detected in more than one-third (70/189 = 37%). The latter are associated with increased probability of recurrence, making their detection important for further pregnancies. Average time from sample receipt to report was 12 days in the recent cases. CONCLUSION: Trio exome sequencing is a useful addition to prenatal diagnostics due to its high diagnostic yield and short processing time (comparable to chromosome analysis). It covers a wide spectrum of genetic changes. Comprehensive interdisciplinary counseling before and after diagnostics is indispensable.

Genetic architecture of inherited retinal degeneration in Germany: A large cohort study from a single diagnostic center over a 9-year period.

We aimed to unravel the molecular genetic basis of inherited retinal degeneration (IRD) in a comprehensive cohort of patients diagnosed in the largest center for IRD in Germany. A cohort of 2,158 affected patients from 1,785 families diagnosed with IRD was analyzed by targeted next-generation sequencing (NGS). Patients with single-gene disorders (i.e., choroideremia and retinoschisis) were analyzed by Sanger sequencing and multiplex ligation-dependent probe amplification. Our study cohort accounts for ∼7% of the estimated 30,000 patients with IRD in Germany, thereby providing representative data for both the prevalence of IRDs and the mutation spectrum of IRD genes for the population in Germany. We achieved a molecular diagnostic rate of 35-95%, depending on the clinical entities, with a high detection rate for achromatopsia, retinoschisis, and choroideremia, and a low detection rate for central areolar choroidal dystrophy and macular dystrophy. A total of 1,161 distinct variants were identified, including 501 novel variants, reaffirming the known vast genetic heterogeneity of IRD in a mainly outbred European population. This study demonstrates the clinical utility of panel-based NGS in a large and highly heterogeneous cohort from an outbred population and for the first time gives a comprehensive representation of the genetic landscape of IRDs in Germany. The data are valuable and crucial for the scientific community and healthcare providers, but also for the pharmaceutical industry in the progressing field of personalized medicine and gene therapy.

Copy number variation profiling in pharmacogenes using panel-based exome resequencing and correlation to human liver expression.

Structural variants including copy number variations (CNV) have gained widespread attention, especially in pharmacogenomics but for several genes functional relevance and clinical evidence are still lacking. Detection of CNVs in next-generation sequencing data is challenging but offers widespread applications. We developed a cohort-based CNV detection workflow to extract CNVs from read counts of targeted NGS of 340 genes involved in absorption, distribution, metabolism and excretion (ADME) of drugs. We applied our method to 150 human liver tissue samples and correlated identified CNVs to mRNA expression levels. In total, we identified 445 deletions (73%) and 167 duplications (27%) in 36 pharmacogenes including all well-known CNVs of CYPs, GSTs, SULTs, UGTs, numerous described rare CNVs of CYP2E1, SLC16A3 or UGT2B15 as well as novel observations, e.g., for SLC22A12, SLC22A17 and GPS2 (G Protein Pathway Suppressor 2). We were able to fine-map complex CNVs of CYP2A6 and CYP2D6 with exon resolution. Correlation analysis confirmed known expression patterns for common CNVs and suggested an influence on expression variability for some rare CNVs. Our straightforward CNV detection workflow can be easily applied to any NGS coverage data and helped to analyze CNVs in an ADME-NGS panel of 340 pharmacogenes to improve genotype-phenotype correlations.

Targeted gene sequencing in 6994 individuals with neurodevelopmental disorder with epilepsy.

PURPOSE: We aimed to gain insight into frequencies of genetic variants in genes implicated in neurodevelopmental disorder with epilepsy (NDD+E) by investigating large cohorts of patients in a diagnostic setting. METHODS: We analyzed variants in NDD+E using epilepsy gene panel sequencing performed between 2013 and 2017 by two large diagnostic companies. We compared variant frequencies in 6994 panels with another 8588 recently published panels as well as exome-wide de novo variants in 1942 individuals with NDD+E and 10,937 controls. RESULTS: Genes with highest frequencies of ultrarare variants in NDD+E comprised SCN1A, KCNQ2, SCN2A, CDKL5, SCN8A, and STXBP1, concordant with the two other epilepsy cohorts we investigated. In only 46% of the analyzed 262 dominant and X-linked panel genes ultrarare variants in patients were reported. Among genes with contradictory evidence of association with epilepsy, CACNB4, CLCN2, EFHC1, GABRD, MAGI2, and SRPX2 showed equal frequencies in cases and controls. CONCLUSION: We show that improvement of panel design increased diagnostic yield over time, but panels still display genes with low or no diagnostic yield. With our data, we hope to improve current diagnostic NDD+E panel design and provide a resource of ultrarare variants in individuals with NDD+E to the community.

Molecular tumor analysis and liquid biopsy: a feasibility investigation analyzing circulating tumor DNA in patients with central nervous system lymphomas.

BACKGROUND: Central nervous system lymphomas (CNSL) is a devastating disease. Currently, a confirmatory biopsy is required prior to treatment. OBJECTIVE: Our investigation aims to prove the feasibility of a minimally-invasive diagnostic approach for the molecular characterization of CNSL. METHODS: Tissue biopsies from 6 patients with suspected CNSL were analyzed using a 649gene next-generation sequencing (NGS) tumor panel (tumor vs. reference tissue (EDTA-blood)). The individual somatic mutation pattern was used as a basis for the digital PCR analyzing circulating tumor DNA (ctDNA) from plasma and cerebrospinal fluid (CSF) samples, identifying one selected tumor mutation during this first step of the feasibility investigation. RESULTS: NGS-analysis of biopsy tissue revealed a specific somatic mutation pattern in all confirmed lymphoma samples (n = 5, NGS-sensitivity 100%) and none in the sample identified as normal brain tissue (NGS-specificity 100%). cfDNA-extraction was dependent on the extraction-kit used and feasible in 3 samples, in all of which somatic mutations were detectable (100%). Analysis of CSF-derived cfDNA was superior to plasma-derived cfDNA and routine microscopic analysis (lymphoma cells: n = 2, 40%). One patient showed a divergent molecular pattern, typical of Burkitt-Lymphoma (HIV+, serologic evidence of EBV-infection). Lumbar puncture was tolerated without complications, whereas biopsy caused 3 hemorrhages. CONCLUSIONS: Our investigation provides evidence that analysis of cfDNA in central nervous system tumors is feasible using the described protocol. Molecular characterization of CNSL could be achieved by analysis of CSF-derived cfDNA. Knowledge of a tumor's specific mutation pattern may allow initiation of targeted therapies, treatment surveillance and could lead to minimally-invasive diagnostics in the future.

Immune monitoring and TCR sequencing of CD4 T cells in a long term responsive patient with metastasized pancreatic ductal carcinoma treated with individualized, neoepitope-derived multipeptide vaccines: a case report.

BACKGROUND: Cancer vaccines can effectively establish clinically relevant tumor immunity. Novel sequencing approaches rapidly identify the mutational fingerprint of tumors, thus allowing to generate personalized tumor vaccines within a few weeks from diagnosis. Here, we report the case of a 62-year-old patient receiving a four-peptide-vaccine targeting the two sole mutations of his pancreatic tumor, identified via exome sequencing. METHODS: Vaccination started during chemotherapy in second complete remission and continued monthly thereafter. We tracked IFN-γ+ T cell responses against vaccine peptides in peripheral blood after 12, 17 and 34 vaccinations by analyzing T-cell receptor (TCR) repertoire diversity and epitope-binding regions of peptide-reactive T-cell lines and clones. By restricting analysis to sorted IFN-γ-producing T cells we could assure epitope-specificity, functionality, and TH1 polarization. RESULTS: A peptide-specific T-cell response against three of the four vaccine peptides could be detected sequentially. Molecular TCR analysis revealed a broad vaccine-reactive TCR repertoire with clones of discernible specificity. Four identical or convergent TCR sequences could be identified at more than one time-point, indicating timely persistence of vaccine-reactive T cells. One dominant TCR expressing a dual TCRVα chain could be found in three T-cell clones. The observed T-cell responses possibly contributed to clinical outcome: The patient is alive 6 years after initial diagnosis and in complete remission for 4 years now. CONCLUSIONS: Therapeutic vaccination with a neoantigen-derived four-peptide vaccine resulted in a diverse and long-lasting immune response against these targets which was associated with prolonged clinical remission. These data warrant confirmation in a larger proof-of concept clinical trial.

Frequent genes in rare diseases: panel-based next generation sequencing to disclose causal mutations in hereditary neuropathies.

Hereditary neuropathies comprise a wide variety of chronic diseases associated to more than 80 genes identified to date. We herein examined 612 index patients with either a Charcot-Marie-Tooth phenotype, hereditary sensory neuropathy, familial amyloid neuropathy, or small fiber neuropathy using a customized multigene panel based on the next generation sequencing technique. In 121 cases (19.8%), we identified at least one putative pathogenic mutation. Of these, 54.4% showed an autosomal dominant, 33.9% an autosomal recessive, and 11.6% an X-linked inheritance. The most frequently affected genes were PMP22 (16.4%), GJB1 (10.7%), MPZ, and SH3TC2 (both 9.9%), and MFN2 (8.3%). We further detected likely or known pathogenic variants in HINT1, HSPB1, NEFL, PRX, IGHMBP2, NDRG1, TTR, EGR2, FIG4, GDAP1, LMNA, LRSAM1, POLG, TRPV4, AARS, BIC2, DHTKD1, FGD4, HK1, INF2, KIF5A, PDK3, REEP1, SBF1, SBF2, SCN9A, and SPTLC2 with a declining frequency. Thirty-four novel variants were considered likely pathogenic not having previously been described in association with any disorder in the literature. In one patient, two homozygous mutations in HK1 were detected in the multigene panel, but not by whole exome sequencing. A novel missense mutation in KIF5A was considered pathogenic because of the highly compatible phenotype. In one patient, the plasma sphingolipid profile could functionally prove the pathogenicity of a mutation in SPTLC2. One pathogenic mutation in MPZ was identified after being previously missed by Sanger sequencing. We conclude that panel based next generation sequencing is a useful, time- and cost-effective approach to assist clinicians in identifying the correct diagnosis and enable causative treatment considerations.

The second report of a new hypomyelinating disease due to a defect in the VPS11 gene discloses a massive lysosomal involvement.

Vesicular protein sorting-associated proteins (VPS, including VPS11) are indispensable in the endocytic network, in particular the endosome-lysosome biogenesis. Exome sequencing revealed the homozygous variant p.Leu387_ Gly395del in the VPS11 gene in two siblings. On immunoblotting, the mutant VPS11 protein showed a distinctly reduced immunostaining intensity. The children presented with primary and severe developmental delay associated with myoclonic seizures, spastic tetraplegia, trunk and neck hypotonia, blindness, hearing loss, and microcephaly. Neuro-imaging showed severe hypomyelination affecting cerebral and cerebellar white matter and corpus callosum, in the absence of a peripheral neuropathy. Electron microscopy of a skin biopsy revealed clusters of membranous cytoplasmic bodies in dermal unmyelinated nerve axons, and numbers of vacuoles in eccrine sweat glands, similar to what is seen in a classic lysosomal storage disease (LSD). Bone marrow cytology showed a high number of storage macrophages with a micro-vacuolated cytoplasm. Biochemically, changes in urinary glycosphingolipids were reminiscent of those in prosaposin deficiency (another LSD). The clinical and neuro-imaged features in our patients were almost identical to those in some recently reported patients with another variant in the VPS11 gene, p.Cys846Gly; underlining the presumed pathogenic potential of VPS11 defects. A new feature was the morphological evidence for lysosomal storage in VPS11 deficiency: This newly characterised disease can be viewed as belonging to the complex field of LSD.

Analytical Performance Evaluation of a 523-Gene Circulating Tumor DNA Assay for Next-Generation Sequencing-Based Comprehensive Tumor Profiling in Liquid Biopsy Samples.

Next-generation sequencing (NGS)-based comprehensive tumor profiling from liquid biopsy samples can significantly improve diagnosis and monitoring of tumors when high-quality tissue material is difficult to obtain. In addition, it offers the potential to capture the entire complexity of the tumor, which is particularly important for highly heterogeneous or metastatic tumors. Here, we report the findings of an analytical performance evaluation of the TruSight Oncology 500 circulating tumor DNA (ctDNA) assay, a 523-gene NGS panel developed for ctDNA-based comprehensive genomic profiling of tumors, using reference and patient samples. Using 30 ng cell-free DNA, the assay showed high sensitivity and low variant detection variability for single-nucleotide variants, insertions and deletions, and fusions down to a variant allele frequency (VAF) of 0.5% in the reference samples and VAFs that were highly concordant with previous digital droplet PCR results in the patient samples. At reduced input amounts (20, 15, and 5 ng) and below VAFs of 0.5%, sensitivity was considerably lower and variant detection variability increased. Covering 523 tumor-associated genes, the assay demonstrated a convincing performance comparable to NGS-based ctDNA assays with smaller gene panels, highlighting its value to screen large numbers of different genes.

Case Report: Long-Term Survival of a Patient with Cerebral Metastasized Ovarian Carcinoma Treated with a Personalized Peptide Vaccine and Anti-PD-1 Therapy.

Ovarian cancer is one of the most common cancers among women and the most lethal malignancy of all gynecological cancers. Surgery is promising in the early stages; however, most patients are first diagnosed in the advanced stages, where treatment options are limited. Here, we present a 49-year-old patient who was first diagnosed with stage III ovarian cancer. After the tumor progressed several times under guideline therapies with no more treatment options available at that time, the patient received a fully individualized neoantigen-derived peptide vaccine in the setting of an individual healing attempt. The tumor was analyzed for somatic mutations via whole exome sequencing and potential neoepitopes were vaccinated over a period of 50 months. During vaccination, the patient additionally received anti-PD-1 therapy to prevent further disease progression. Vaccine-induced T-cell responses were detected using intracellular cytokine staining. After eleven days of in vitro expansion, four T-cell activation markers (namely IFN-É£, TNF-α, IL-2, and CD154) were measured. The proliferation capacity of neoantigen-specific T-cells was determined using a CFSE proliferation assay. Immune monitoring revealed a very strong CD4+ T-cell response against one of the vaccinated peptides. The vaccine-induced T-cells simultaneously expressed CD154, TNF, IL-2, and IFN-É£ and showed a strong proliferation capacity upon neoantigen stimulation. Next-generation sequencing, as well as immunohistochemical analysis, revealed a loss of Beta-2 microglobulin (B2M), which is essential for MHC class I presentation. The results presented here implicate that the application of neoantigen-derived peptide vaccines might be considered for those cancer stages, where promising therapeutic options are lacking. Furthermore, we provide more data that endorse the intensive investigation of B2M loss as a tumor escape mechanism in clinical trials using anti-cancer vaccines together with immune-checkpoint inhibitors.

REVEAL--visual eQTL analytics.

MOTIVATION: The analysis of expression quantitative trait locus (eQTL) data is a challenging scientific endeavor, involving the processing of very large, heterogeneous and complex data. Typical eQTL analyses involve three types of data: sequence-based data reflecting the genotypic variations, gene expression data and meta-data describing the phenotype. Based on these, certain genotypes can be connected with specific phenotypic outcomes to infer causal associations of genetic variation, expression and disease. To this end, statistical methods are used to find significant associations between single nucleotide polymorphisms (SNPs) or pairs of SNPs and gene expression. A major challenge lies in summarizing the large amount of data as well as statistical results and to generate informative, interactive visualizations. RESULTS: We present Reveal, our visual analytics approach to this challenge. We introduce a graph-based visualization of associations between SNPs and gene expression and a detailed genotype view relating summarized patient cohort genotypes with data from individual patients and statistical analyses. AVAILABILITY: Reveal is included in Mayday, our framework for visual exploration and analysis. It is available at http://it.inf.uni-tuebingen.de/software/reveal/. CONTACT: guenter.jaeger@uni-tuebingen.de.

The role of rare genetic variants enrichment in epilepsies of presumed genetic etiology.

Previous studies suggested that severe epilepsies e.g., developmental and epileptic encephalopathies (DEE) are mainly caused by ultra-rare de novo genetic variants. For milder phenotypes, rare genetic variants could contribute to the phenotype. To determine the importance of rare variants for different epilepsy types, we analyzed a whole-exome sequencing cohort of 9,170 epilepsy-affected individuals and 8,436 controls. Here, we separately analyzed three different groups of epilepsies : severe DEEs, genetic generalized epilepsy (GGE), and non-acquired focal epilepsy (NAFE). We required qualifying rare variants (QRVs) to occur in controls at a minor allele frequency ≤ 1:1,000, to be predicted as deleterious (CADD≥20), and to have an odds ratio in epilepsy cases ≥2. We identified genes enriched with QRVs in DEE (n=21), NAFE (n=72), and GGE (n=32) - the number of enriched genes are found greatest in NAFE and least in DEE. This suggests that rare variants may play a more important role for causality of NAFE than in DEE. Moreover, we found that QRV-carrying genes e.g., HSGP2, FLNA or TNC are involved in structuring the brain extracellular matrix. The present study confirms an involvement of rare variants for NAFE, while in DEE and GGE, the contribution of such variants appears more limited.

Tübingen model study: large-scale introduction of rapid antigen testing in the population and the viral dynamics of SARS-CoV-2.

Despite of contact restrictions, population mobility remains the main reason for the spread of SARS-CoV-2. The state of Baden-Württemberg (BW), Germany, approved a model study in Tübingen (TÜMOD) to evaluate how mandatory rapid diagnostic tests (RDT) could reduce transmission. Between 16 March and 24 April 2021, approximately 165,000 residents and visitors to the city were screened for SARS CoV-2 infection using Abbott Panbio™ COVID-19 Antigen rapid test device. We assessed incidences and recorded epidemiological characteristics in a subset of 4,118 participants recruited at three of the nine testing stations. PCR tests were performed in RDT-positives to determine the positive predictive value (PPV), and circulating variants of SARS-CoV-2 were identified by whole-genome sequencing. 2,282 RDT-negative samples were tested by pooled PCR to calculate the false negative rate (FNR). Viral load was compared between variants. 116 (3%) participants were positive by RDT, and of these, 57 (49%) were positive by PCR, 55 (47%) were negative. This resulted in a PPV of 51%. Of the 57 positives, 52 SARS-CoV-2 genomes were successfully sequenced. Of these, 50 belonged to the B.1.1.7 lineage, which had a high viral load (average Ct = 19). Of the 2,282 RDT negatives tested, all were PCR negative (FNR 0%). At the end of TÜMOD, the incidence in Tübingen, which was initially lower, had reached the incidence in the state of BW. While it is difficult to assess the impact of TÜMOD on incidence independent of confounding factors, further studies are needed to identify the effect of close-meshed testing on infection rates.

Case Report: Targeting of individual somatic tumor mutations by multipeptide vaccination tailored for HLA class I and II presentation induces strong CD4 and CD8 T-cell responses in a patient with metastatic castration sensitive prostate cancer.

Localized prostate cancer is curable, but metastatic castration sensitive prostate cancer has a low 5-year survival rate, while broad treatment options are lacking. Here we present an mCSPC patient under remission receiving individualized neoantigen-derived peptide vaccination as recurrence prophylaxis in the setting of an individual treatment attempt. The patient was initially analyzed for somatic tumor mutations and then consecutively treated with two different peptide vaccines over a period of 33 months. The first vaccine contained predicted HLA class I binding peptides only whereas the second vaccine contained both predicted HLA class I and II binding peptides. Intracellular cytokine staining after 12 day in-vitro expansion measuring four T-cell activation markers (IFNg, TNF-α, IL-2, CD154) was used to determine vaccine-induced T-cell responses. While the first vaccine induced only one robust CD4+ T-cell response after 21 vaccinations, co-vaccination of HLA class I and II peptides induced multiple strong and durable CD4+ and CD8+ T-cell responses already after sixth vaccinations. The vaccine-induced immune responses were robust and polyfunctional. PSA remained undetectable for 51 months. The results presented here implicate that neoantigen-targeting vaccines might be considered for those cancer subtypes where therapeutic options are limited. Furthermore, our findings suggest that both HLA class I and II restricted peptides should be considered for future peptide vaccination trials.

iHAT: interactive hierarchical aggregation table for genetic association data.

In the search for single-nucleotide polymorphisms which influence the observable phenotype, genome wide association studies have become an important technique for the identification of associations between genotype and phenotype of a diverse set of sequence-based data. We present a methodology for the visual assessment of single-nucleotide polymorphisms using interactive hierarchical aggregation techniques combined with methods known from traditional sequence browsers and cluster heatmaps. Our tool, the interactive Hierarchical Aggregation Table (iHAT), facilitates the visualization of multiple sequence alignments, associated metadata, and hierarchical clusterings. Different color maps and aggregation strategies as well as filtering options support the user in finding correlations between sequences and metadata. Similar to other visualizations such as parallel coordinates or heatmaps, iHAT relies on the human pattern-recognition ability for spotting patterns that might indicate correlation or anticorrelation. We demonstrate iHAT using artificial and real-world datasets for DNA and protein association studies as well as expression Quantitative Trait Locus data.

An eQTL biological data visualization challenge and approaches from the visualization community.

In 2011, the IEEE VisWeek conferences inaugurated a symposium on Biological Data Visualization. Like other domain-oriented Vis symposia, this symposium's purpose was to explore the unique characteristics and requirements of visualization within the domain, and to enhance both the Visualization and Bio/Life-Sciences communities by pushing Biological data sets and domain understanding into the Visualization community, and well-informed Visualization solutions back to the Biological community. Amongst several other activities, the BioVis symposium created a data analysis and visualization contest. Unlike many contests in other venues, where the purpose is primarily to allow entrants to demonstrate tour-de-force programming skills on sample problems with known solutions, the BioVis contest was intended to whet the participants' appetites for a tremendously challenging biological domain, and simultaneously produce viable tools for a biological grand challenge domain with no extant solutions. For this purpose expression Quantitative Trait Locus (eQTL) data analysis was selected. In the BioVis 2011 contest, we provided contestants with a synthetic eQTL data set containing real biological variation, as well as a spiked-in gene expression interaction network influenced by single nucleotide polymorphism (SNP) DNA variation and a hypothetical disease model. Contestants were asked to elucidate the pattern of SNPs and interactions that predicted an individual's disease state. 9 teams competed in the contest using a mixture of methods, some analytical and others through visual exploratory methods. Independent panels of visualization and biological experts judged entries. Awards were given for each panel's favorite entry, and an overall best entry agreed upon by both panels. Three special mention awards were given for particularly innovative and useful aspects of those entries. And further recognition was given to entries that correctly answered a bonus question about how a proposed "gene therapy" change to a SNP might change an individual's disease status, which served as a calibration for each approaches' applicability to a typical domain question. In the future, BioVis will continue the data analysis and visualization contest, maintaining the philosophy of providing new challenging questions in open-ended and dramatically underserved Bio/Life Sciences domains.