Precision medicine screening using whole-genome sequencing and advanced imaging to identify disease risk in adults.

Reducing premature mortality associated with age-related chronic diseases, such as cancer and cardiovascular disease, is an urgent priority. We report early results using genomics in combination with advanced imaging and other clinical testing to proactively screen for age-related chronic disease risk among adults. We enrolled active, symptom-free adults in a study of screening for age-related chronic diseases associated with premature mortality. In addition to personal and family medical history and other clinical testing, we obtained whole-genome sequencing (WGS), noncontrast whole-body MRI, dual-energy X-ray absorptiometry (DXA), global metabolomics, a new blood test for prediabetes (Quantose IR), echocardiography (ECHO), ECG, and cardiac rhythm monitoring to identify age-related chronic disease risks. Precision medicine screening using WGS and advanced imaging along with other testing among active, symptom-free adults identified a broad set of complementary age-related chronic disease risks associated with premature mortality and strengthened WGS variant interpretation. This and other similarly designed screening approaches anchored by WGS and advanced imaging may have the potential to extend healthy life among active adults through improved prevention and early detection of age-related chronic diseases (and their risk factors) associated with premature mortality.

The human noncoding genome defined by genetic diversity.

Understanding the significance of genetic variants in the noncoding genome is emerging as the next challenge in human genomics. We used the power of 11,257 whole-genome sequences and 16,384 heptamers (7-nt motifs) to build a map of sequence constraint for the human species. This build differed substantially from traditional maps of interspecies conservation and identified regulatory elements among the most constrained regions of the genome. Using new Hi-C experimental data, we describe a strong pattern of coordination over 2 Mb where the most constrained regulatory elements associate with the most essential genes. Constrained regions of the noncoding genome are up to 52-fold enriched for known pathogenic variants as compared to unconstrained regions (21-fold when compared to the genome average). This map of sequence constraint across thousands of individuals is an asset to help interpret noncoding elements in the human genome, prioritize variants and reconsider gene units at a larger scale.

Profiling of Short-Tandem-Repeat Disease Alleles in 12,632 Human Whole Genomes.

Short tandem repeats (STRs) are hyper-mutable sequences in the human genome. They are often used in forensics and population genetics and are also the underlying cause of many genetic diseases. There are challenges associated with accurately determining the length polymorphism of STR loci in the genome by next-generation sequencing (NGS). In particular, accurate detection of pathological STR expansion is limited by the sequence read length during whole-genome analysis. We developed TREDPARSE, a software package that incorporates various cues from read alignment and paired-end distance distribution, as well as a sequence stutter model, in a probabilistic framework to infer repeat sizes for genetic loci, and we used this software to infer repeat sizes for 30 known disease loci. Using simulated data, we show that TREDPARSE outperforms other available software. We sampled the full genome sequences of 12,632 individuals to an average read depth of approximately 30× to 40× with Illumina HiSeq X. We identified 138 individuals with risk alleles at 15 STR disease loci. We validated a representative subset of the samples (n = 19) by Sanger and by Oxford Nanopore sequencing. Additionally, we validated the STR calls against known allele sizes in a set of GeT-RM reference cell-line materials (n = 6). Several STR loci that are entirely guanine or cytosines (G or C) have insufficient read evidence for inference and therefore could not be assayed precisely by TREDPARSE. TREDPARSE extends the limit of STR size detection beyond the physical sequence read length. This extension is critical because many of the disease risk cutoffs are close to or beyond the short sequence read length of 100 to 150 bases.

Fast and accurate HLA typing from short-read next-generation sequence data with xHLA.

The HLA gene complex on human chromosome 6 is one of the most polymorphic regions in the human genome and contributes in large part to the diversity of the immune system. Accurate typing of HLA genes with short-read sequencing data has historically been difficult due to the sequence similarity between the polymorphic alleles. Here, we introduce an algorithm, xHLA, that iteratively refines the mapping results at the amino acid level to achieve 99-100% four-digit typing accuracy for both class I and II HLA genes, taking only [Formula: see text]3 min to process a 30× whole-genome BAM file on a desktop computer.

Whole-genome sequencing identifies common-to-rare variants associated with human blood metabolites.

Genetic factors modifying the blood metabolome have been investigated through genome-wide association studies (GWAS) of common genetic variants and through exome sequencing. We conducted a whole-genome sequencing study of common, low-frequency and rare variants to associate genetic variations with blood metabolite levels using comprehensive metabolite profiling in 1,960 adults. We focused the analysis on 644 metabolites with consistent levels across three longitudinal data collections. Genetic sequence variations at 101 loci were associated with the levels of 246 (38%) metabolites (P ≤ 1.9 × 10-11). We identified 113 (10.7%) among 1,054 unrelated individuals in the cohort who carried heterozygous rare variants likely influencing the function of 17 genes. Thirteen of the 17 genes are associated with inborn errors of metabolism or other pediatric genetic conditions. This study extends the map of loci influencing the metabolome and highlights the importance of heterozygous rare variants in determining abnormal blood metabolic phenotypes in adults.

Deep sequencing of 10,000 human genomes.

We report on the sequencing of 10,545 human genomes at 30×-40× coverage with an emphasis on quality metrics and novel variant and sequence discovery. We find that 84% of an individual human genome can be sequenced confidently. This high-confidence region includes 91.5% of exon sequence and 95.2% of known pathogenic variant positions. We present the distribution of over 150 million single-nucleotide variants in the coding and noncoding genome. Each newly sequenced genome contributes an average of 8,579 novel variants. In addition, each genome carries on average 0.7 Mb of sequence that is not found in the main build of the hg38 reference genome. The density of this catalog of variation allowed us to construct high-resolution profiles that define genomic sites that are highly intolerant of genetic variation. These results indicate that the data generated by deep genome sequencing is of the quality necessary for clinical use.

Analytical performance validation of a coronary heart disease risk assessment multi-analyte proteomic test.

BACKGROUND: Coronary heart disease (CHD) remains prevalent despite efforts to improve CHD risk assessment. The authors developed a multi-analyte immunoassay-based CHD risk assessment (CHDRA) algorithm, clinically validated in a multicenter study, to improve CHDRA in intermediate risk individuals. OBJECTIVE: Clinical laboratory validation of the CHDRA biomarker assays' analytical performance. METHODS: Multiplexed immunoassay panels developed for the seven CHDRA assays were evaluated with donor sera in a clinical laboratory. Specificity, sensitivity, interfering substances and reproducibility of the CHDRA assays, along with the effects of pre-analytical specimen processing, were evaluated. RESULTS: Analytical measurements of the CHDRA panel proteins (CTACK, Eotaxin, Fas Ligand, HGF, IL-16, MCP-3 and sFas) exhibited acceptable accuracy (80 - 120%), cross-reactivity (< 1%), interference (< 30% at high concentrations of bilirubin, lipids, hemoglobin and HAMA), sensitivity and reproducibility (< 20% CV across multiple runs, operators and instruments). Recoveries from donor sera subjected to typical clinical laboratory pre-analytical conditions were within 80 - 120%. The pre-analytical variables did not substantively impact the CHDRA scores. CONCLUSIONS: The CHDRA panel analytical validation in a clinical laboratory meets or exceeds the specifications established during the clinical utility studies. Risk score reproducibility across multiple test scenarios suggests the assays are not susceptible to clinical laboratory pre-analytical and analytical variation.

Company profile: Aviir, Inc.

Aviir, Inc. is a venture-funded biotechnology company developing and commercializing laboratory tests to provide personalized information to physicians and patients, with the goal of preventing cardiovascular disease and metabolic syndromes. Leveraging advanced research, Aviir developed and launched MIRISK VP™, a risk assessment test to better identify individuals at risk of a heart attack. Aviir also offers an extensive menu of other cardiovascular and metabolic tests through its Clinical Laboratory Improvement Amendments (CLIA)-certified laboratory. Efforts are likewise focused on expanding genomics testing capability to address sudden cardiac death attributed to inherited cardiovascular diseases. This completes their integrated precision diagnostics approach that combines biomarker immunoassays with genomic and transcription analysis, along with core clinical chemistry to deliver a comprehensive personal health solution.

Tuning a three-component reaction for trapping kinase substrate complexes.

The upstream protein kinases responsible for thousands of phosphorylation events in the phosphoproteome remain to be discovered. We developed a three-component chemical reaction which converts the transient noncovalent substrate-kinase complex into a covalently cross-linked product by utilizing a dialdehyde-based cross-linker, 1. Unfortunately, the reaction of 1 with a lysine in the kinase active site and an engineered cysteine on the substrate to form an isoindole cross-linked product could not be performed in the presence of competing cellular proteins due to nonspecific side reactions. In order to more selectively target the cross-linker to protein kinases in cell lysates, we replaced the weak, kinase-binding adenosine moiety of 1 with a potent protein kinase inhibitor scaffold. In addition, we replaced the o-phthaldialdehyde moiety in 1 with a less-reactive thiophene-2,3-dicarboxaldehyde moiety. The combination of these two structural modifications provides for cross-linking of a cysteine-containing substrate to its corresponding kinase in the presence of competing cellular proteins.

Structure of the kinase domain of an imatinib-resistant Abl mutant in complex with the Aurora kinase inhibitor VX-680.

We present a high-resolution (2.0 A) crystal structure of the catalytic domain of a mutant form of the Abl tyrosine kinase (H396P; Abl-1a numbering) that is resistant to the Abl inhibitor imatinib. The structure is determined in complex with the small-molecule inhibitor VX-680 (Vertex Pharmaceuticals, Cambridge, MA), which blocks the activity of various imatinib-resistant mutant forms of Abl, including one (T315I) that is resistant to both imatinib and BMS-354825 (dasatinib), a dual Src/Abl inhibitor that seems to be clinically effective against all other imatinib-resistant forms of BCR-Abl. VX-680 is shown to have significant inhibitory activity against BCR-Abl bearing the T315I mutation in patient-derived samples. The Abl kinase domain bound to VX-680 is not phosphorylated on the activation loop in the crystal structure but is nevertheless in an active conformation, previously unobserved for Abl and inconsistent with the binding of imatinib. The adoption of an active conformation is most likely the result of synergy between the His(396)Pro mutation, which destabilizes the inactive conformation required for imatinib binding, and the binding of VX-680, which favors the active conformation through hydrogen bonding and steric effects. VX-680 is bound to Abl in a mode that accommodates the substitution of isoleucine for threonine at residue 315 (the "gatekeeper" position). The avoidance of the innermost cavity of the Abl kinase domain by VX-680 and the specific recognition of the active conformation explain the effectiveness of this compound against mutant forms of BCR-Abl, including those with mutations at the gatekeeper position.

Inhibition of drug-resistant mutants of ABL, KIT, and EGF receptor kinases.

To realize the full potential of targeted protein kinase inhibitors for the treatment of cancer, it is important to address the emergence of drug resistance in treated patients. Mutant forms of BCR-ABL, KIT, and the EGF receptor (EGFR) have been found that confer resistance to the drugs imatinib, gefitinib, and erlotinib. The mutations weaken or prevent drug binding, and interestingly, one of the most common sites of mutation in all three kinases is a highly conserved "gatekeeper" threonine residue near the kinase active site. We have identified existing clinical compounds that bind and inhibit drug-resistant mutant variants of ABL, KIT, and EGFR. We found that the Aurora kinase inhibitor VX-680 and the p38 inhibitor BIRB-796 inhibit the imatinib- and BMS-354825-resistant ABL(T315I) kinase. The KIT/FLT3 inhibitor SU-11248 potently inhibits the imatinib-resistant KIT(V559D/T670I) kinase, consistent with the clinical efficacy of SU-11248 against imatinib-resistant gastrointestinal tumors, and the EGFR inhibitors EKB-569 and CI-1033, but not GW-572016 and ZD-6474, potently inhibit the gefitinib- and erlotinib-resistant EGFR(L858R/T790M) kinase. EKB-569 and CI-1033 are already in clinical trials, and our results suggest that they should be considered for testing in the treatment of gefitinib/erlotinib-resistant non-small cell lung cancer. The results highlight the strategy of screening existing clinical compounds against newly identified drug-resistant mutant variants to find compounds that may serve as starting points for the development of next-generation drugs, or that could be used directly to treat patients that have acquired resistance to first-generation targeted therapy.

A small molecule-kinase interaction map for clinical kinase inhibitors.

Kinase inhibitors show great promise as a new class of therapeutics. Here we describe an efficient way to determine kinase inhibitor specificity by measuring binding of small molecules to the ATP site of kinases. We have profiled 20 kinase inhibitors, including 16 that are approved drugs or in clinical development, against a panel of 119 protein kinases. We find that specificity varies widely and is not strongly correlated with chemical structure or the identity of the intended target. Many novel interactions were identified, including tight binding of the p38 inhibitor BIRB-796 to an imatinib-resistant variant of the ABL kinase, and binding of imatinib to the SRC-family kinase LCK. We also show that mutations in the epidermal growth factor receptor (EGFR) found in gefitinib-responsive patients do not affect the binding affinity of gefitinib or erlotinib. Our results represent a systematic small molecule-protein interaction map for clinical compounds across a large number of related proteins.

Disruption of forkhead transcription factor (FOXO) family members in mice reveals their functional diversification.

Genetic analysis in Caenorhabditis elegans has uncovered essential roles for DAF-16 in longevity, metabolism, and reproduction. The mammalian orthologs of DAF-16, the closely-related FOXO subclass of forkhead transcription factors (FKHR/FOXO1, FKHRL1/FOXO3a, and AFX/FOXO4), also have important roles in cell cycle arrest, apoptosis and stress responses in vitro, but their in vivo physiological roles are largely unknown. To elucidate their role in normal development and physiology, we disrupted each of the Foxo genes in mice. Foxo1-null embryos died on embryonic day 10.5 as a consequence of incomplete vascular development. Foxo1-null embryonic and yolk sac vessels were not well developed at embryonic day 9.5, and Foxo1 expression was found in a variety of embryonic vessels, suggesting a crucial role of this transcription factor in vascular formation. On the other hand, both Foxo3a- and Foxo4-null mice were viable and grossly indistinguishable from their littermate controls, indicating dispensability of these two members of the Foxo transcription factor family for normal vascular development. Foxo3a-null females showed age-dependent infertility and had abnormal ovarian follicular development. In contrast, histological analyses of Foxo4-null mice did not identify any consistent abnormalities. These results demonstrate that the physiological roles of Foxo genes are functionally diverse in mammals.

Identification and characterization of mouse SSX genes: a multigene family on the X chromosome with restricted cancer/testis expression.

Human SSX was first identified as the gene involved in the t(X;18) translocation in synovial sarcoma. SSX is a multigene family, with 9 complete genes on chromosome Xp11. Normally expressed almost exclusively in testis, SSX mRNA is expressed in various human tumors, defining SSX as a cancer/testis antigen. We have now cloned the mouse ortholog of SSX. Mouse SSX genes can be divided into Ssxa and Ssxb subfamilies based on sequence homology. Ssxa has only one member, whereas 12 Ssxb genes, Ssxb1 to Ssxb12, were identified by cDNA cloning from mouse testis and mouse tumors. Both Ssxa and Ssxb are located on chromosome X and show tissue-restricted mRNA expression to testis among normal tissues. All putative human and mouse SSX proteins share conserved KRAB and SSX-RD domains. Mouse tumors were found to express some, but not all, Ssxb genes, similar to the SSX activation in human tumors.

The forkhead transcription factor Foxo1 regulates adipocyte differentiation.

An outstanding question in adipocyte biology is how hormonal cues are relayed to the nucleus to activate the transcriptional program that promotes adipogenesis. The forkhead transcription factor Foxo1 is regulated by insulin via Akt-dependent phosphorylation and nuclear exclusion. We show that Foxo1 is induced in the early stages of adipocyte differentiation but that its activation is delayed until the end of the clonal expansion phase. Constitutively active Foxo1 prevents the differentiation of preadipocytes, while dominant-negative Foxo1 restores adipocyte differentiation of fibroblasts from insulin receptor-deficient mice. Further, Foxo1 haploinsufficiency protects from diet-induced diabetes in mice. We propose that Foxo1 plays an important role in the integration of hormone-activated signaling pathways with the complex transcriptional cascade that promotes adipocyte differentiation.

The forkhead transcription factor Foxo1 links insulin signaling to Pdx1 regulation of pancreatic beta cell growth.

Diabetes is caused by an absolute (type 1) or relative (type 2) deficiency of insulin-producing beta cells. The mechanisms governing replication of terminally differentiated beta cells and neogenesis from progenitor cells are unclear. Mice lacking insulin receptor substrate-2 (Irs2) develop beta cell failure, suggesting that insulin signaling is required to maintain an adequate beta cell mass. We report that haploinsufficiency for the forkhead transcription factor Foxo1 reverses beta cell failure in Irs2(-/-) mice through partial restoration of beta cell proliferation and increased expression of the pancreatic transcription factor pancreas/duodenum homeobox gene-1 (Pdx1). Foxo1 and Pdx1 exhibit mutually exclusive patterns of nuclear localization in beta cells, and constitutive nuclear expression of a mutant Foxo1 is associated with lack of Pdx1 expression. We show that Foxo1 acts as a repressor of Foxa2-dependent (Hnf-3beta-dependent) expression from the Pdx1 promoter. We propose that insulin/IGFs regulate beta cell proliferation by relieving Foxo1 inhibition of Pdx1 expression in a subset of cells embedded within pancreatic ducts.

Regulation of insulin action and pancreatic beta-cell function by mutated alleles of the gene encoding forkhead transcription factor Foxo1.

Type 2 diabetes results from impaired action and secretion of insulin. It is not known whether the two defects share a common pathogenesis. We show that haploinsufficiency of the Foxo1 gene, encoding a forkhead transcription factor (forkhead box transcription factor O1), restores insulin sensitivity and rescues the diabetic phenotype in insulin-resistant mice by reducing hepatic expression of glucogenetic genes and increasing adipocyte expression of insulin-sensitizing genes. Conversely, a gain-of-function Foxo1 mutation targeted to liver and pancreatic beta-cells results in diabetes arising from a combination of increased hepatic glucose production and impaired beta-cell compensation due to decreased Pdx1 expression. These data indicate that Foxo1 is a negative regulator of insulin sensitivity in liver, adipocytes and pancreatic beta-cells. Impaired insulin signaling to Foxo1 provides a unifying mechanism for the common metabolic abnormalities of type 2 diabetes.NOTE: In the AOP version of this article, the name of the fourth author was misspelled as W K Cavanee rather than the correct spelling: W K Cavenee. This has been corrected in the full-text online version of the article. The name will appear correctly in the print version.

Analysis of the regulation of the A33 antigen gene reveals intestine-specific mechanisms of gene expression.

The A33 antigen is a transmembrane protein expressed almost exclusively by intestinal epithelial cells. The level of its expression is robust and uniform throughout the rostrocaudal axis of the human and mouse intestines. In the colon, strong expression is found in the basolateral membranes of both the proliferating cells in the lower regions of the crypts and the differentiating cells in the upper regions of crypts. Similarly, in the small intestine, the protein is highly expressed by all the epithelial cells in the crypts and by the differentiated cells migrating over the villi. Thus, the A33 antigen has emerged as a definitive marker for all intestinal epithelial cells, irrespective of cell lineage and differentiation status. To understand the molecular mechanisms mediating this rare tissue-specific expression pattern, we undertook a comprehensive analysis of the 5'-regulatory region of the human A33 antigen gene. This allowed us to point to positive cis-regulatory elements incorporating consensus Krüppel-like factor and caudal-related homeobox (CDX)-binding sites, located just upstream from the human A33 antigen transcription start site, as being important for the intestine-specific expression pattern of this gene. Further analysis provided evidence that the A33 antigen gene may be one of only a few target genes to be described thus far for the intestine-specific homeobox transcription factor, CDX1. Taken together, our data lead us to propose that the activity of CDX1 is pivotal in mediating the exquisite, intestine-specific expression pattern of the A33 antigen gene.