Short tandem repeat expansions in cortical layer-specific genes implicate in phenotypic severity and adaptability of autism spectrum disorder.

AIM: Short tandem repeats (STRs) are repetitive DNA sequences and highly mutable in various human disorders. While the involvement of STRs in various genetic disorders has been extensively studied, their role in autism spectrum disorder (ASD) remains largely unexplored. In this study, we aimed to investigate genetic association of STR expansions with ASD using whole genome sequencing (WGS) and identify risk loci associated with ASD phenotypes. METHODS: We analyzed WGS data of 634 ASD families and performed genome-wide evaluation for 12,929 STR loci. We found rare STR expansions that exceeded normal repeat lengths in autism cases compared to unaffected controls. By integrating single cell RNA and ATAC sequencing datasets of human postmortem brains, we prioritized STR loci in genes specifically expressed in cortical development stages. A deep learning method was used to predict functionality of ASD-associated STR loci. RESULTS: In ASD cases, rare STR expansions predominantly occurred in early cortical layer-specific genes involved in neurodevelopment, highlighting the cellular specificity of STR-associated genes in ASD risk. Leveraging deep learning prediction models, we demonstrated that these STR expansions disrupted the regulatory activity of enhancers and promoters, suggesting a potential mechanism through which they contribute to ASD pathogenesis. We found that individuals with ASD-associated STR expansions exhibited more severe ASD phenotypes and diminished adaptability compared to non-carriers. CONCLUSION: Short tandem repeat expansions in cortical layer-specific genes are associated with ASD and could potentially be a risk genetic factor for ASD. Our study is the first to show evidence of STR expansion associated with ASD in an under-investigated population.

Widespread somatic L1 retrotransposition in normal colorectal epithelium.

Throughout an individual's lifetime, genomic alterations accumulate in somatic cells1-11. However, the mutational landscape induced by retrotransposition of long interspersed nuclear element-1 (L1), a widespread mobile element in the human genome12-14, is poorly understood in normal cells. Here we explored the whole-genome sequences of 899 single-cell clones established from three different cell types collected from 28 individuals. We identified 1,708 somatic L1 retrotransposition events that were enriched in colorectal epithelium and showed a positive relationship with age. Fingerprinting of source elements showed 34 retrotransposition-competent L1s. Multidimensional analysis demonstrated that (1) somatic L1 retrotranspositions occur from early embryogenesis at a substantial rate, (2) epigenetic on/off of a source element is preferentially determined in the early organogenesis stage, (3) retrotransposition-competent L1s with a lower population allele frequency have higher retrotransposition activity and (4) only a small fraction of L1 transcripts in the cytoplasm are finally retrotransposed in somatic cells. Analysis of matched cancers further suggested that somatic L1 retrotransposition rate is substantially increased during colorectal tumourigenesis. In summary, this study illustrates L1 retrotransposition-induced somatic mosaicism in normal cells and provides insights into the genomic and epigenomic regulation of transposable elements over the human lifetime.

Ultra-Low Level Somatic Mutations and Structural Variations in Focal Cortical Dysplasia Type II.

OBJECTIVE: Brain somatic mutations in mTOR pathway genes are a major genetic etiology of focal cortical dysplasia type II (FCDII). Despite a greater ability to detect low-level somatic mutations in the brain by deep sequencing and analytics, about 40% of cases remain genetically unexplained. METHODS: We included 2 independent cohorts consisting of 21 patients with mutation-negative FCDII without apparent mutations on conventional deep sequencing of bulk brain. To find ultra-low level somatic variants or structural variants, we isolated cells exhibiting phosphorylation of the S6 ribosomal protein (p-S6) in frozen brain tissues using fluorescence-activated cell sorting (FACS). We then performed deep whole-genome sequencing (WGS; >90×) in p-S6+ cells in a cohort of 11 patients with mutation-negative. Then, we simplified the method to whole-genome amplification and target gene sequencing of p-S6+ cells in independent cohort of 10 patients with mutation-negative followed by low-read depth WGS (10×). RESULTS: We found that 28.6% (6 of 21) of mutation-negative FCDII carries ultra-low level somatic mutations (less than 0.2% of variant allele frequency [VAF]) in mTOR pathway genes. Our method showed ~34 times increase of the average mutational burden in FACS mediated enrichment of p-S6+ cells (average VAF = 5.84%) than in bulky brain tissues (average VAF = 0.17%). We found that 19% (4 of 21) carried germline structural variations in GATOR1 complex undetectable in whole exome or targeted gene sequencing. CONCLUSIONS: Our method facilitates the detection of ultra-low level somatic mutations, in specifically p-S6+ cells, and germline structural variations and increases the genetic diagnostic rate up to ~80% for the entire FCDII cohort. ANN NEUROL 2023;93:1082-1093.

Ultrafast prediction of somatic structural variations by filtering out reads matched to pan-genome k-mer sets.

Variant callers typically produce massive numbers of false positives for structural variations, such as cancer-relevant copy-number alterations and fusion genes resulting from genome rearrangements. Here we describe an ultrafast and accurate detector of somatic structural variations that reduces read-mapping costs by filtering out reads matched to pan-genome k-mer sets. The detector, which we named ETCHING (for efficient detection of chromosomal rearrangements and fusion genes), reduces the number of false positives by leveraging machine-learning classifiers trained with six breakend-related features (clipped-read count, split-reads count, supporting paired-end read count, average mapping quality, depth difference and total length of clipped bases). When benchmarked against six callers on reference cell-free DNA, validated biomarkers of structural variants, matched tumour and normal whole genomes, and tumour-only targeted sequencing datasets, ETCHING was 11-fold faster than the second-fastest structural-variant caller at comparable performance and memory use. The speed and accuracy of ETCHING may aid large-scale genome projects and facilitate practical implementations in precision medicine.

Identification of a pleiotropic effect of ADIPOQ on cardiac dysfunction and Alzheimer's disease based on genetic evidence and health care records.

Observations of comorbidity in heart diseases, including cardiac dysfunction (CD) are increasing, including and cognitive impairment, such as Alzheimer's disease and dementia (AD/D). This comorbidity might be due to a pleiotropic effect of genetic variants shared between CD and AD/D. Here, we validated comorbidity of CD and AD/D based on diagnostic records from millions of patients in Korea and the University of California, San Francisco Medical Center (odds ratio 11.5 [8.5-15.5, 95% Confidence Interval (CI)]). By integrating a comprehensive human disease-SNP association database (VARIMED, VARiants Informing MEDicine) and whole-exome sequencing of 50 brains from individuals with and without Alzheimer's disease (AD), we identified missense variants in coding regions including APOB, a known risk factor for CD and AD/D, which potentially have a pleiotropic role in both diseases. Of the identified variants, site-directed mutation of ADIPOQ (268 G > A; Gly90Ser) in neurons produced abnormal aggregation of tau proteins (p = 0.02), suggesting a functional impact for AD/D. The association of CD and ADIPOQ variants was confirmed based on domain deletion in cardiac cells. Using the UK Biobank including data from over 500000 individuals, we examined a pleiotropic effect of the ADIPOQ variant by comparing CD- and AD/D-associated phenotypic evidence, including cardiac hypertrophy and cognitive degeneration. These results indicate that convergence of health care records and genetic evidences may help to dissect the molecular underpinnings of heart disease and associated cognitive impairment, and could potentially serve a prognostic function. Validation of disease-disease associations through health care records and genomic evidence can determine whether health conditions share risk factors based on pleiotropy.

Analysis of low-level somatic mosaicism reveals stage and tissue-specific mutational features in human development.

Most somatic mutations that arise during normal development are present at low levels in single or multiple tissues depending on the developmental stage and affected organs. However, the effect of human developmental stages or mutations of different organs on the features of somatic mutations is still unclear. Here, we performed a systemic and comprehensive analysis of low-level somatic mutations using deep whole-exome sequencing (average read depth ~500×) of 498 multiple organ tissues with matched controls from 190 individuals. Our results showed that early clone-forming mutations shared between multiple organs were lower in number but showed higher allele frequencies than late clone-forming mutations [0.54 vs. 5.83 variants per individual; 6.17% vs. 1.5% variant allele frequency (VAF)] along with less nonsynonymous mutations and lower functional impacts. Additionally, early and late clone-forming mutations had unique mutational signatures that were distinct from mutations that originated from tumors. Compared with early clone-forming mutations that showed a clock-like signature across all organs or tissues studied, late clone-forming mutations showed organ, tissue, and cell-type specificity in the mutation counts, VAFs, and mutational signatures. In particular, analysis of brain somatic mutations showed a bimodal occurrence and temporal-lobe-specific signature. These findings provide new insights into the features of somatic mosaicism that are dependent on developmental stage and brain regions.

Non-coding de novo mutations in chromatin interactions are implicated in autism spectrum disorder.

Three-dimensional chromatin interactions regulate gene expressions. The significance of de novo mutations (DNMs) in chromatin interactions remains poorly understood for autism spectrum disorder (ASD). We generated 813 whole-genome sequences from 242 Korean simplex families to detect DNMs, and identified target genes which were putatively affected by non-coding DNMs in chromatin interactions. Non-coding DNMs in chromatin interactions were significantly involved in transcriptional dysregulations related to ASD risk. Correspondingly, target genes showed spatiotemporal expressions relevant to ASD in developing brains and enrichment in biological pathways implicated in ASD, such as histone modification. Regarding clinical features of ASD, non-coding DNMs in chromatin interactions particularly contributed to low intelligence quotient levels in ASD probands. We further validated our findings using two replication cohorts, Simons Simplex Collection (SSC) and MSSNG, and showed the consistent enrichment of non-coding DNM-disrupted chromatin interactions in ASD probands. Generating human induced pluripotent stem cells in two ASD families, we were able to demonstrate that non-coding DNMs in chromatin interactions alter the expression of target genes at the stage of early neural development. Taken together, our findings indicate that non-coding DNMs in ASD probands lead to early neurodevelopmental disruption implicated in ASD risk via chromatin interactions.

Estimation of intrafamilial DNA contamination in family trio genome sequencing using deviation from Mendelian inheritance.

With the increasing number of sequencing projects involving families, quality control tools optimized for family genome sequencing are needed. However, accurately quantifying contamination in a DNA mixture is particularly difficult when genetically related family members are the sources. We developed TrioMix, a maximum likelihood estimation (MLE) framework based on Mendel's law of inheritance, to quantify DNA mixture between family members in genome sequencing data of parent-offspring trios. TrioMix can accurately deconvolute any intrafamilial DNA contamination, including parent-offspring, sibling-sibling, parent-parent, and even multiple familial sources. In addition, TrioMix can be applied to detect genomic abnormalities that deviate from Mendelian inheritance patterns, such as uniparental disomy (UPD) and chimerism. A genome-wide depth and variant allele frequency plot generated by TrioMix facilitates tracing the origin of Mendelian inheritance deviations. We showed that TrioMix could accurately deconvolute genomes in both simulated and real data sets.

Clonal dynamics in early human embryogenesis inferred from somatic mutation.

Cellular dynamics and fate decision in early human embryogenesis remain largely unknown owing to the challenges of performing studies in human embryos1. Here, we explored whole-genomes of 334 single-cell colonies and targeted deep sequences of 379 bulk tissues obtained from various anatomical locations of seven recently deceased adult human donors. Using somatic mutations as an intrinsic barcode, we reconstructed early cellular phylogenies that demonstrate (1) an endogenous mutational rate that is higher in the first cell division but decreases to approximately one per cell per cell division later in life; (2) universal unequal contribution of early cells to embryo proper, resulting from early cellular bottlenecks that stochastically set aside epiblast cells within the embryo; (3) examples of varying degrees of early clonal imbalances between tissues on the left and right sides of the body, different germ layers and specific anatomical parts and organs; (4) emergence of a few ancestral cells that will substantially contribute to adult cell pools in blood and liver; and (5) presence of mitochondrial DNA heteroplasmy in the fertilized egg. Our approach also provides insights into the age-related mutational processes and loss of sex chromosomes in normal somatic cells. In sum, this study provides a foundation for future studies to complete cellular phylogenies in human embryogenesis.

Novel Alzheimer's disease risk variants identified based on whole-genome sequencing of APOE ε4 carriers.

Alzheimer's disease (AD) is a progressive neurodegenerative disease associated with a complex genetic etiology. Besides the apolipoprotein E ε4 (APOE ε4) allele, a few dozen other genetic loci associated with AD have been identified through genome-wide association studies (GWAS) conducted mainly in individuals of European ancestry. Recently, several GWAS performed in other ethnic groups have shown the importance of replicating studies that identify previously established risk loci and searching for novel risk loci. APOE-stratified GWAS have yielded novel AD risk loci that might be masked by, or be dependent on, APOE alleles. We performed whole-genome sequencing (WGS) on DNA from blood samples of 331 AD patients and 169 elderly controls of Korean ethnicity who were APOE ε4 carriers. Based on WGS data, we designed a customized AD chip (cAD chip) for further analysis on an independent set of 543 AD patients and 894 elderly controls of the same ethnicity, regardless of their APOE ε4 allele status. Combined analysis of WGS and cAD chip data revealed that SNPs rs1890078 (P = 6.64E-07) and rs12594991 (P = 2.03E-07) in SORCS1 and CHD2 genes, respectively, are novel genetic variants among APOE ε4 carriers in the Korean population. In addition, nine possible novel variants that were rare in individuals of European ancestry but common in East Asia were identified. This study demonstrates that APOE-stratified analysis is important for understanding the genetic background of AD in different populations.

Low-Level Brain Somatic Mutations Are Implicated in Schizophrenia.

BACKGROUND: Somatic mutations arising from the brain have recently emerged as significant contributors to neurodevelopmental disorders, including childhood intractable epilepsy and cortical malformations. However, whether brain somatic mutations are implicated in schizophrenia (SCZ) is not well established. METHODS: We performed deep whole exome sequencing (average read depth > 550×) of matched dorsolateral prefrontal cortex and peripheral tissues from 27 patients with SCZ and 31 age-matched control individuals, followed by comprehensive and strict analysis of somatic mutations, including mutagenesis signature, substitution patterns, and involved pathways. In particular, we explored the impact of deleterious mutations in GRIN2B through primary neural culture. RESULTS: We identified an average of 4.9 and 5.6 somatic mutations per exome per brain in patients with SCZ and control individuals, respectively. These mutations presented with average variant allele frequencies of 8.0% in patients with SCZ and 7.6% in control individuals. Although mutational profiles, such as the number and type of mutations, showed no significant difference between patients with SCZ and control individuals, somatic mutations in SCZ brains were significantly enriched for SCZ-related pathways, including dopamine receptor, glutamate receptor, and long-term potentiation pathways. Furthermore, we showed that brain somatic mutations in GRIN2B (encoding glutamate ionotropic NMDA receptor subunit 2B), which were found in two patients with SCZ, disrupted the location of GRIN2B across the surface of dendrites among primary cultured neurons. CONCLUSIONS: Taken together, this study shows that brain somatic mutations are associated with the pathogenesis of SCZ.

Brain somatic mutations observed in Alzheimer's disease associated with aging and dysregulation of tau phosphorylation.

The role of brain somatic mutations in Alzheimer's disease (AD) is not well understood. Here, we perform deep whole-exome sequencing (average read depth 584×) in 111 postmortem hippocampal formation and matched blood samples from 52 patients with AD and 11 individuals not affected by AD. The number of somatic single nucleotide variations (SNVs) in AD brain specimens increases significantly with aging, and the rate of mutation accumulation in the brain is 4.8-fold slower than that in AD blood. The putatively pathogenic brain somatic mutations identified in 26.9% (14 of 52) of AD individuals are enriched in PI3K-AKT, MAPK, and AMPK pathway genes known to contribute to hyperphosphorylation of tau. We show that a pathogenic brain somatic mutation in PIN1 leads to a loss-of-function mutation. In vitro mimicking of haploinsufficiency of PIN1 aberrantly increases tau phosphorylation and aggregation. This study provides new insights into the genetic architecture underlying the pathogenesis of AD.

Tracing Oncogene Rearrangements in the Mutational History of Lung Adenocarcinoma.

Mutational processes giving rise to lung adenocarcinomas (LADCs) in non-smokers remain elusive. We analyzed 138 LADC whole genomes, including 83 cases with minimal contribution of smoking-associated mutational signature. Genomic rearrangements were not correlated with smoking-associated mutations and frequently served as driver events of smoking-signature-low LADCs. Complex genomic rearrangements, including chromothripsis and chromoplexy, generated 74% of known fusion oncogenes, including EML4-ALK, CD74-ROS1, and KIF5B-RET. Unlike other collateral rearrangements, these fusion-oncogene-associated rearrangements were frequently copy-number-balanced, representing a genomic signature of early oncogenesis. Analysis of mutation timing revealed that fusions and point mutations of canonical oncogenes were often acquired in the early decades of life. During a long latency, cancer-related genes were disrupted or amplified by complex rearrangements. The genomic landscape was different between subgroups-EGFR-mutant LADCs had frequent whole-genome duplications with p53 mutations, whereas fusion-oncogene-driven LADCs had frequent SETD2 mutations. Our study highlights LADC oncogenesis driven by endogenous mutational processes.

Brain Somatic Mutations in MTOR Disrupt Neuronal Ciliogenesis, Leading to Focal Cortical Dyslamination.

Focal malformations of cortical development (FMCDs), including focal cortical dysplasia (FCD) and hemimegalencephaly (HME), are major etiologies of pediatric intractable epilepsies exhibiting cortical dyslamination. Brain somatic mutations in MTOR have recently been identified as a major genetic cause of FMCDs. However, the molecular mechanism by which these mutations lead to cortical dyslamination remains poorly understood. Here, using patient tissue, genome-edited cells, and mouse models with brain somatic mutations in MTOR, we discovered that disruption of neuronal ciliogenesis by the mutations underlies cortical dyslamination in FMCDs. We found that abnormal accumulation of OFD1 at centriolar satellites due to perturbed autophagy was responsible for the defective neuronal ciliogenesis. Additionally, we found that disrupted neuronal ciliogenesis accounted for cortical dyslamination in FMCDs by compromising Wnt signals essential for neuronal polarization. Altogether, this study describes a molecular mechanism by which brain somatic mutations in MTOR contribute to the pathogenesis of cortical dyslamination in FMCDs.

Clonal History and Genetic Predictors of Transformation Into Small-Cell Carcinomas From Lung Adenocarcinomas.

Purpose Histologic transformation of EGFR mutant lung adenocarcinoma (LADC) into small-cell lung cancer (SCLC) has been described as one of the major resistant mechanisms for epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors (TKIs). However, the molecular pathogenesis is still unclear. Methods We investigated 21 patients with advanced EGFR-mutant LADCs that were transformed into EGFR TKI-resistant SCLCs. Among them, whole genome sequencing was applied for nine tumors acquired at various time points from four patients to reconstruct their clonal evolutionary history and to detect genetic predictors for small-cell transformation. The findings were validated by immunohistochemistry in 210 lung cancer tissues. Results We identified that EGFR TKI-resistant LADCs and SCLCs share a common clonal origin and undergo branched evolutionary trajectories. The clonal divergence of SCLC ancestors from the LADC cells occurred before the first EGFR TKI treatments, and the complete inactivation of both RB1 and TP53 were observed from the early LADC stages in sequenced tumors. We extended the findings by immunohistochemistry in the early-stage LADC tissues of 75 patients treated with EGFR TKIs; inactivation of both Rb and p53 was strikingly more frequent in the small-cell-transformed group than in the nontransformed group (82% v 3%; odds ratio, 131; 95% CI, 19.9 to 859). Among patients registered in a predefined cohort (n = 65), an EGFR mutant LADC that harbored completely inactivated Rb and p53 had a 43× greater risk of small-cell transformation (relative risk, 42.8; 95% CI, 5.88 to 311). Branch-specific mutational signature analysis revealed that apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC)-induced hypermutation was frequent in the branches toward small-cell transformation. Conclusion EGFR TKI-resistant SCLCs are branched out early from the LADC clones that harbor completely inactivated RB1 and TP53. The evaluation of RB1 and TP53 status in EGFR TKI-treated LADCs is informative in predicting small-cell transformation.

Somatic Mutations in TSC1 and TSC2 Cause Focal Cortical Dysplasia.

Focal cortical dysplasia (FCD) is a major cause of the sporadic form of intractable focal epilepsies that require surgical treatment. It has recently been reported that brain somatic mutations in MTOR account for 15%-25% of FCD type II (FCDII), characterized by cortical dyslamination and dysmorphic neurons. However, the genetic etiologies of FCDII-affected individuals who lack the MTOR mutation remain unclear. Here, we performed deep hybrid capture and amplicon sequencing (read depth of 100×-20,012×) of five important mTOR pathway genes-PIK3CA, PIK3R2, AKT3, TSC1, and TSC2-by using paired brain and saliva samples from 40 FCDII individuals negative for MTOR mutations. We found that 5 of 40 individuals (12.5%) had brain somatic mutations in TSC1 (c.64C>T [p.Arg22Trp] and c.610C>T [p.Arg204Cys]) and TSC2 (c.4639G>A [p.Val1547Ile]), and these results were reproducible on two different sequencing platforms. All identified mutations induced hyperactivation of the mTOR pathway by disrupting the formation or function of the TSC1-TSC2 complex. Furthermore, in utero CRISPR-Cas9-mediated genome editing of Tsc1 or Tsc2 induced the development of spontaneous behavioral seizures, as well as cytomegalic neurons and cortical dyslamination. These results show that brain somatic mutations in TSC1 and TSC2 cause FCD and that in utero application of the CRISPR-Cas9 system is useful for generating neurodevelopmental disease models of somatic mutations in the brain.

Vecuum: identification and filtration of false somatic variants caused by recombinant vector contamination.

MOTIVATION: Advances in sequencing technologies have remarkably lowered the detection limit of somatic variants to a low frequency. However, calling mutations at this range is still confounded by many factors including environmental contamination. Vector contamination is a continuously occurring issue and is especially problematic since vector inserts are hardly distinguishable from the sample sequences. Such inserts, which may harbor polymorphisms and engineered functional mutations, can result in calling false variants at corresponding sites. Numerous vector-screening methods have been developed, but none could handle contamination from inserts because they are focusing on vector backbone sequences alone. RESULTS: We developed a novel method-Vecuum-that identifies vector-originated reads and resultant false variants. Since vector inserts are generally constructed from intron-less cDNAs, Vecuum identifies vector-originated reads by inspecting the clipping patterns at exon junctions. False variant calls are further detected based on the biased distribution of mutant alleles to vector-originated reads. Tests on simulated and spike-in experimental data validated that Vecuum could detect 93% of vector contaminants and could remove up to 87% of variant-like false calls with 100% precision. Application to public sequence datasets demonstrated the utility of Vecuum in detecting false variants resulting from various types of external contamination. AVAILABILITY AND IMPLEMENTATION: Java-based implementation of the method is available at http://vecuum.sourceforge.net/ CONTACT: swkim@yuhs.acSupplementary information: Supplementary data are available at Bioinformatics online.

Multilevel Genomics-Based Taxonomy of Renal Cell Carcinoma.

On the basis of multidimensional and comprehensive molecular characterization (including DNA methalylation and copy number, RNA, and protein expression), we classified 894 renal cell carcinomas (RCCs) of various histologic types into nine major genomic subtypes. Site of origin within the nephron was one major determinant in the classification, reflecting differences among clear cell, chromophobe, and papillary RCC. Widespread molecular changes associated with TFE3 gene fusion or chromatin modifier genes were present within a specific subtype and spanned multiple subtypes. Differences in patient survival and in alteration of specific pathways (including hypoxia, metabolism, MAP kinase, NRF2-ARE, Hippo, immune checkpoint, and PI3K/AKT/mTOR) could further distinguish the subtypes. Immune checkpoint markers and molecular signatures of T cell infiltrates were both highest in the subtype associated with aggressive clear cell RCC. Differences between the genomic subtypes suggest that therapeutic strategies could be tailored to each RCC disease subset.

Association analysis of the perturbation of interactions in biological pathways and anticancer drug activity.

Understanding how different genomic mutational landscapes in patients with cancer lead to different responses to anticancer drugs is an important challenge for realizing precision medicine for cancer. Many studies have analyzed the comprehensive anticancer drug-response profiles and genomic profiles of cancer cell lines to identify the relationship between the anticancer drug response and genomic alternations. However, few studies have focused on interpreting these profiles with a network perspective. In this work, we analyzed genomic alterations in cancer cell lines by considering which interactions in the signaling pathway were perturbed by mutations. With our interaction-centric approach, we identified novel interaction/drug response associations for two drugs (afatinib and ixabepilone) for which no gene-centric association could be found. When we compared the performance of classifiers for predicting the responses to 164 drugs, the classifiers trained with interaction-centric features outperformed the classifiers trained with gene-centric features, despite the smaller number of features (p-value = 2.0 × 10(-3)). By incorporating the interaction information from signaling pathways, we revealed associations between genomic alterations and drug responses that could be missed when using a gene-centric approach.

Sensitivity Analysis of the Parameters of Korea's Pandemic Influenza Preparedness Plan.

OBJECTIVES: Our aim was to evaluate Korea's Pandemic Influenza Preparedness Plan. METHODS: We conducted a sensitivity analysis on the expected number of outpatients and hospital bed occupancy, with 1,000,000 parameter combinations, in a situation of pandemic influenza, using the mathematical simulation program InfluSim. RESULTS: Given the available resources in Korea, antiviral treatment and social distancing must be combined to reduce the number of outpatients and hospitalizations sufficiently; any single intervention is not enough. The antiviral stockpile of 4-6% is sufficient for the expected eligible number of cases to be treated. However, the eligible number assumed (30% for severe cases and 26% for extremely severe cases) is very low compared to the corresponding number in European countries, where up to 90% of the population are assumed to be eligible for antiviral treatment. CONCLUSIONS: A combination of antiviral treatment and social distancing can mitigate a pandemic, but will only bring it under control for the most optimistic parameter combinations.