Artificial intelligence enables comprehensive genome interpretation and nomination of candidate diagnoses for rare genetic diseases.

BACKGROUND: Clinical interpretation of genetic variants in the context of the patient's phenotype is becoming the largest component of cost and time expenditure for genome-based diagnosis of rare genetic diseases. Artificial intelligence (AI) holds promise to greatly simplify and speed genome interpretation by integrating predictive methods with the growing knowledge of genetic disease. Here we assess the diagnostic performance of Fabric GEM, a new, AI-based, clinical decision support tool for expediting genome interpretation. METHODS: We benchmarked GEM in a retrospective cohort of 119 probands, mostly NICU infants, diagnosed with rare genetic diseases, who received whole-genome or whole-exome sequencing (WGS, WES). We replicated our analyses in a separate cohort of 60 cases collected from five academic medical centers. For comparison, we also analyzed these cases with current state-of-the-art variant prioritization tools. Included in the comparisons were trio, duo, and singleton cases. Variants underpinning diagnoses spanned diverse modes of inheritance and types, including structural variants (SVs). Patient phenotypes were extracted from clinical notes by two means: manually and using an automated clinical natural language processing (CNLP) tool. Finally, 14 previously unsolved cases were reanalyzed. RESULTS: GEM ranked over 90% of the causal genes among the top or second candidate and prioritized for review a median of 3 candidate genes per case, using either manually curated or CNLP-derived phenotype descriptions. Ranking of trios and duos was unchanged when analyzed as singletons. In 17 of 20 cases with diagnostic SVs, GEM identified the causal SVs as the top candidate and in 19/20 within the top five, irrespective of whether SV calls were provided or inferred ab initio by GEM using its own internal SV detection algorithm. GEM showed similar performance in absence of parental genotypes. Analysis of 14 previously unsolved cases resulted in a novel finding for one case, candidates ultimately not advanced upon manual review for 3 cases, and no new findings for 10 cases. CONCLUSIONS: GEM enabled diagnostic interpretation inclusive of all variant types through automated nomination of a very short list of candidate genes and disorders for final review and reporting. In combination with deep phenotyping by CNLP, GEM enables substantial automation of genetic disease diagnosis, potentially decreasing cost and expediting case review.

Targeted Identification of Protein Interactions in Eukaryotic mRNA Translation.

To identify protein-protein interactions and phosphorylated amino acid sites in eukaryotic mRNA translation, replicate TAP-MudPIT and control experiments are performed targeting Saccharomyces cerevisiae genes previously implicated in eukaryotic mRNA translation by their genetic and/or functional roles in translation initiation, elongation, termination, or interactions with ribosomal complexes. Replicate tandem affinity purifications of each targeted yeast TAP-tagged mRNA translation protein coupled with multidimensional liquid chromatography and tandem mass spectrometry analysis are used to identify and quantify copurifying proteins. To improve sensitivity and minimize spurious, nonspecific interactions, a novel cross-validation approach is employed to identify the most statistically significant protein-protein interactions. Using experimental and computational strategies discussed herein, the previously described protein composition of the canonical eukaryotic mRNA translation initiation, elongation, and termination complexes is calculated. In addition, statistically significant unpublished protein interactions and phosphorylation sites for S. cerevisiae's mRNA translation proteins and complexes are identified.

Next-Generation Sequencing Strategies.

More than a decade ago, the term "next-generation" sequencing was coined to describe what was, at the time, revolutionary new methods to sequence RNA and DNA at a faster pace and cheaper cost than could be performed by standard bench-top protocols. Since then, the field of DNA sequencing has evolved at a rapid pace, with new breakthroughs allowing capacity to exponentially increase and cost to dramatically decrease. As genome-scale sequencing has become routine, a paradigm shift is occurring in genomics, which uses the power of high-throughput, rapid sequencing power with large-scale studies. These new approaches to genetic discovery will provide direct impact to fields such as personalized medicine, evolution, and biodiversity. This work reviews recent technology advances and methods in next-generation sequencing and highlights current large-scale sequencing efforts driving the evolution of the genomics space.

Reversible deficits in apical transporter trafficking associated with deficiency in diacylglycerol acyltransferase.

Deficiency in diacylglycerol acyltransferase (DGAT1) is a rare cause of neonatal diarrhea, without a known mechanism or in vitro model. A patient presenting at our institution at 7 weeks of life with failure to thrive and diarrhea was found by whole-exome sequencing to have a homozygous DGAT1 truncation mutation. Duodenal biopsies showed loss of DGAT1 and deficits in apical membrane transporters and junctional proteins in enterocytes. When placed on a very low-fat diet, the patient's diarrhea resolved with normalization of brush border transporter localization in endoscopic biopsies. DGAT1 knockdown in Caco2-BBe cells modeled the deficits in apical trafficking, with loss of apical DPPIV and junctional occludin. Elevation in cellular lipid levels, including diacylglycerol (DAG) and phospholipid metabolites of DAG, was documented by lipid analysis in DGAT1 knockdown cells. Culture of the DGAT1 knockdown cells in lipid-depleted media led to re-establishment of occludin and return of apical DPPIV. DGAT1 loss appears to elicit global changes in enterocyte polarized trafficking that could account for deficits in absorption seen in the patient. The in vitro modeling of this disease should allow for investigation of possible therapeutic targets.

Genetic Risk for Subsequent Neoplasms Among Long-Term Survivors of Childhood Cancer.

Purpose Childhood cancer survivors are at increased risk of subsequent neoplasms (SNs), but the germline genetic contribution is largely unknown. We assessed the contribution of pathogenic/likely pathogenic (P/LP) mutations in cancer predisposition genes to their SN risk. Patients and Methods Whole-genome sequencing (30-fold) was performed on samples from childhood cancer survivors who were ≥ 5 years since initial cancer diagnosis and participants in the St Jude Lifetime Cohort Study, a retrospective hospital-based study with prospective clinical follow-up. Germline mutations in 60 genes known to be associated with autosomal dominant cancer predisposition syndromes with moderate to high penetrance were classified by their pathogenicity according to the American College of Medical Genetics and Genomics guidelines. Relative rates (RRs) and 95% CIs of SN occurrence by mutation status were estimated using multivariable piecewise exponential regression stratified by radiation exposure. Results Participants were 3,006 survivors (53% male; median age, 35.8 years [range, 7.1 to 69.8 years]; 56% received radiotherapy), 1,120 SNs were diagnosed among 439 survivors (14.6%), and 175 P/LP mutations were identified in 5.8% (95% CI, 5.0% to 6.7%) of survivors. Mutations were associated with significantly increased rates of breast cancer (RR, 13.9; 95% CI, 6.0 to 32.2) and sarcoma (RR, 10.6; 95% CI, 4.3 to 26.3) among irradiated survivors and with increased rates of developing any SN (RR, 4.7; 95% CI, 2.4 to 9.3), breast cancer (RR, 7.7; 95% CI, 2.4 to 24.4), nonmelanoma skin cancer (RR, 11.0; 95% CI, 2.9 to 41.4), and two or more histologically distinct SNs (RR, 18.6; 95% CI, 3.5 to 99.3) among nonirradiated survivors. Conclusion The findings support referral of all survivors for genetic counseling for potential clinical genetic testing, which should be prioritized for nonirradiated survivors with any SN and for those with breast cancer or sarcoma in the field of prior irradiation.

Ecophysiological Examination of the Lake Erie Microcystis Bloom in 2014: Linkages between Biology and the Water Supply Shutdown of Toledo, OH.

Annual cyanobacterial blooms dominated by Microcystis have occurred in western Lake Erie (U.S./Canada) during summer months since 1995. The production of toxins by bloom-forming cyanobacteria can lead to drinking water crises, such as the one experienced by the city of Toledo in August of 2014, when the city was rendered without drinking water for >2 days. It is important to understand the conditions and environmental cues that were driving this specific bloom to provide a scientific framework for management of future bloom events. To this end, samples were collected and metatranscriptomes generated coincident with the collection of environmental metrics for eight sites located in the western basin of Lake Erie, including a station proximal to the water intake for the city of Toledo. These data were used to generate a basin-wide ecophysiological fingerprint of Lake Erie Microcystis populations in August 2014 for comparison to previous bloom communities. Our observations and analyses indicate that, at the time of sample collection, Microcystis populations were under dual nitrogen (N) and phosphorus (P) stress, as genes involved in scavenging of these nutrients were being actively transcribed. Targeted analysis of urea transport and hydrolysis suggests a potentially important role for exogenous urea as a nitrogen source during the 2014 event. Finally, simulation data suggest a wind event caused microcystin-rich water from Maumee Bay to be transported east along the southern shoreline past the Toledo water intake. Coupled with a significant cyanophage infection, these results reveal that a combination of biological and environmental factors led to the disruption of the Toledo water supply. This scenario was not atypical of reoccurring Lake Erie blooms and thus may reoccur in the future.

A high-throughput molecular data resource for cutaneous neurofibromas.

Neurofibromatosis type 1 (NF1) is a genetic disorder with a range of clinical manifestations such as widespread growth of benign tumours called neurofibromas, pain, learning disorders, bone deformities, vascular abnormalities and even malignant tumours. With the establishment of the Children's Tumour Foundation biobank, neurofibroma samples can now be collected directly from patients to be analysed by the larger scientific community. This work describes a pilot study to characterize one class of neurofibroma, cutaneous neurofibromas, by molecularly profiling of ~40 cutaneous neurofibromas collected from 11 individual patients. Data collected from each tumour includes (1) SNP Arrays, (2) Whole genome sequencing (WGS) and (3) RNA-Sequencing. These data are now freely available for further analysis at http://www.synapse.org/cutaneousNF.

Genomic regulation of invasion by STAT3 in triple negative breast cancer.

Breast cancer is a heterogeneous disease comprised of four molecular subtypes defined by whether the tumor-originating cells are luminal or basal epithelial cells. Breast cancers arising from the luminal mammary duct often express estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth receptor 2 (HER2). Tumors expressing ER and/or PR are treated with anti-hormonal therapies, while tumors overexpressing HER2 are targeted with monoclonal antibodies. Immunohistochemical detection of ER, PR, and HER2 receptors/proteins is a critical step in breast cancer diagnosis and guided treatment. Breast tumors that do not express these proteins are known as "triple negative breast cancer" (TNBC) and are typically basal-like. TNBCs are the most aggressive subtype, with the highest mortality rates and no targeted therapy, so there is a pressing need to identify important TNBC tumor regulators. The signal transducer and activator of transcription 3 (STAT3) transcription factor has been previously implicated as a constitutively active oncogene in TNBC. However, its direct regulatory gene targets and tumorigenic properties have not been well characterized. By integrating RNA-seq and ChIP-seq data from 2 TNBC tumors and 5 cell lines, we discovered novel gene signatures directly regulated by STAT3 that were enriched for processes involving inflammation, immunity, and invasion in TNBC. Functional analysis revealed that STAT3 has a key role regulating invasion and metastasis, a characteristic often associated with TNBC. Our findings suggest therapies targeting STAT3 may be important for preventing TNBC metastasis.

Geospatial Resolution of Human and Bacterial Diversity with City-Scale Metagenomics.

The panoply of microorganisms and other species present in our environment influence human health and disease, especially in cities, but have not been profiled with metagenomics at a city-wide scale. We sequenced DNA from surfaces across the entire New York City (NYC) subway system, the Gowanus Canal, and public parks. Nearly half of the DNA (48%) does not match any known organism; identified organisms spanned 1,688 bacterial, viral, archaeal, and eukaryotic taxa, which were enriched for harmless genera associated with skin (e.g., Acinetobacter). Predicted ancestry of human DNA left on subway surfaces can recapitulate U.S. Census demographic data, and bacterial signatures can reveal a station's history, such as marine-associated bacteria in a hurricane-flooded station. Some evidence of pathogens was found (Bacillus anthracis), but a lack of reported cases in NYC suggests that the pathogens represent a normal, urban microbiome. This baseline metagenomic map of NYC could help long-term disease surveillance, bioterrorism threat mitigation, and health management in the built environment of cities.

Exome sequencing in amyotrophic lateral sclerosis identifies risk genes and pathways.

Amyotrophic lateral sclerosis (ALS) is a devastating neurological disease with no effective treatment. We report the results of a moderate-scale sequencing study aimed at increasing the number of genes known to contribute to predisposition for ALS. We performed whole-exome sequencing of 2869 ALS patients and 6405 controls. Several known ALS genes were found to be associated, and TBK1 (the gene encoding TANK-binding kinase 1) was identified as an ALS gene. TBK1 is known to bind to and phosphorylate a number of proteins involved in innate immunity and autophagy, including optineurin (OPTN) and p62 (SQSTM1/sequestosome), both of which have also been implicated in ALS. These observations reveal a key role of the autophagic pathway in ALS and suggest specific targets for therapeutic intervention.

Intragenic motifs regulate the transcriptional complexity of Pkhd1/PKHD1.

Autosomal recessive polycystic kidney disease (ARPKD) results from mutations in the human PKHD1 gene. Both this gene, and its mouse ortholog, Pkhd1, are primarily expressed in renal and biliary ductal structures. The mouse protein product, fibrocystin/polyductin complex (FPC), is a 445-kDa protein encoded by a 67-exon transcript that spans >500 kb of genomic DNA. In the current study, we observed multiple alternatively spliced Pkhd1 transcripts that varied in size and exon composition in embryonic mouse kidney, liver, and placenta samples, as well as among adult mouse pancreas, brain, heart, lung, testes, liver, and kidney. Using reverse transcription PCR and RNASeq, we identified 22 novel Pkhd1 kidney transcripts with unique exon junctions. Various mechanisms of alternative splicing were observed, including exon skipping, use of alternate acceptor/donor splice sites, and inclusion of novel exons. Bioinformatic analyses identified, and exon-trapping minigene experiments validated, consensus binding sites for serine/arginine-rich proteins that modulate alternative splicing. Using site-directed mutagenesis, we examined the functional importance of selected splice enhancers. In addition, we demonstrated that many of the novel transcripts were polysome bound, thus likely translated. Finally, we determined that the human PKHD1 R760H missense variant alters a splice enhancer motif that disrupts exon splicing in vitro and is predicted to truncate the protein. Taken together, these data provide evidence of the complex transcriptional regulation of Pkhd1/PKHD1 and identified motifs that regulate its splicing. Our studies indicate that Pkhd1/PKHD1 transcription is modulated, in part by intragenic factors, suggesting that aberrant PKHD1 splicing represents an unappreciated pathogenic mechanism in ARPKD. Key messages: Multiple mRNA transcripts are generated for Pkhd1 in renal tissues Pkhd1 transcription is modulated by standard splice elements and effectors Mutations in splice motifs may alter splicing to generate nonfunctional peptides.

De novo gene mutations highlight patterns of genetic and neural complexity in schizophrenia.

To evaluate evidence for de novo etiologies in schizophrenia, we sequenced at high coverage the exomes of families recruited from two populations with distinct demographic structures and history. We sequenced a total of 795 exomes from 231 parent-proband trios enriched for sporadic schizophrenia cases, as well as 34 unaffected trios. We observed in cases an excess of de novo nonsynonymous single-nucleotide variants as well as a higher prevalence of gene-disruptive de novo mutations relative to controls. We found four genes (LAMA2, DPYD, TRRAP and VPS39) affected by recurrent de novo events within or across the two populations, which is unlikely to have occurred by chance. We show that de novo mutations affect genes with diverse functions and developmental profiles, but we also find a substantial contribution of mutations in genes with higher expression in early fetal life. Our results help define the genomic and neural architecture of schizophrenia.

Patterns and rates of exonic de novo mutations in autism spectrum disorders.

Autism spectrum disorders (ASD) are believed to have genetic and environmental origins, yet in only a modest fraction of individuals can specific causes be identified. To identify further genetic risk factors, here we assess the role of de novo mutations in ASD by sequencing the exomes of ASD cases and their parents (n = 175 trios). Fewer than half of the cases (46.3%) carry a missense or nonsense de novo variant, and the overall rate of mutation is only modestly higher than the expected rate. In contrast, the proteins encoded by genes that harboured de novo missense or nonsense mutations showed a higher degree of connectivity among themselves and to previous ASD genes as indexed by protein-protein interaction screens. The small increase in the rate of de novo events, when taken together with the protein interaction results, are consistent with an important but limited role for de novo point mutations in ASD, similar to that documented for de novo copy number variants. Genetic models incorporating these data indicate that most of the observed de novo events are unconnected to ASD; those that do confer risk are distributed across many genes and are incompletely penetrant (that is, not necessarily sufficient for disease). Our results support polygenic models in which spontaneous coding mutations in any of a large number of genes increases risk by 5- to 20-fold. Despite the challenge posed by such models, results from de novo events and a large parallel case-control study provide strong evidence in favour of CHD8 and KATNAL2 as genuine autism risk factors.

Targeted genome enrichment for efficient purification of endosymbiont DNA from host DNA.

Wolbachia endosymbionts are widespread in arthropods and are generally considered reproductive parasites, inducing various phenotypes including cytoplasmic incompatibility, parthenogenesis, feminization and male killing, which serve to promote their spread through populations. In contrast, Wolbachia infecting filarial nematodes that cause human diseases, including elephantiasis and river blindness, are obligate mutualists. DNA purification methods for efficient genomic sequencing of these unculturable bacteria have proven difficult using a variety of techniques. To efficiently capture endosymbiont DNA for studies that examine the biology of symbiosis, we devised a parallel strategy to an earlier array-based method by creating a set of SureSelect™ (Agilent) 120-mer target enrichment RNA oligonucleotides ("baits") for solution hybrid selection. These were designed from Wolbachia complete and partial genome sequences in GenBank and were tiled across each genomic sequence with 60 bp overlap. Baits were filtered for homology against host genomes containing Wolbachia using BLAT and sequences with significant host homology were removed from the bait pool. Filarial parasite Brugia malayi DNA was used as a test case, as the complete sequence of both Wolbachia and its host are known. DNA eluted from capture was size selected and sequencing samples were prepared using the NEBNext® Sample Preparation Kit. One-third of a 50 nt paired-end sequencing lane on the HiSeq™ 2000 (Illumina) yielded 53 million reads and the entirety of the Wolbachia genome was captured. We then used the baits to isolate more than 97.1 % of the genome of a distantly related Wolbachia strain from the crustacean Armadillidium vulgare, demonstrating that the method can be used to enrich target DNA from unculturable microbes over large evolutionary distances.

Exome sequencing supports a de novo mutational paradigm for schizophrenia.

Despite its high heritability, a large fraction of individuals with schizophrenia do not have a family history of the disease (sporadic cases). Here we examined the possibility that rare de novo protein-altering mutations contribute to the genetic component of schizophrenia by sequencing the exomes of 53 sporadic cases, 22 unaffected controls and their parents. We identified 40 de novo mutations in 27 cases affecting 40 genes, including a potentially disruptive mutation in DGCR2, a gene located in the schizophrenia-predisposing 22q11.2 microdeletion region. A comparison to rare inherited variants indicated that the identified de novo mutations show a large excess of non-synonymous changes in schizophrenia cases, as well as a greater potential to affect protein structure and function. Our analyses suggest a major role for de novo mutations in schizophrenia as well as a large mutational target, which together provide a plausible explanation for the high global incidence and persistence of the disease.

Genomic profiling of C/EBPß2 transformed mammary epithelial cells: a role for nuclear interleukin-1ß.

C/EBPß is essential for mammary gland growth and development and has been associated with poor prognosis in breast cancer. Overexpression of C/EBPß2 in MCF10A cells results in a variety of cancer phenotypes including EMT and ErbB independence. IL1ß is dramatically upregulated in MCF10A-C/EBPß2 cells but there is little, if any, processing to the mature 17 kD form. Although proIL1b has previously been considered to be biologically inactive, we demonstrate proIL1b is not only localized to the nucleus, but is also tightly associated with the chromatin. We show that proIL1ß is bound at specific locations in the genome and is positioned in such a way to play a role in the cancer phenotypes observed in MCF10A-C/EBPß2 cells. Moreover, nuclear IL1ß is detected in some human breast tumor samples. This study demonstrates the presence of nuclear proIL1ß in transformed mammary epithelial cells providing the first evidence that IL1ß may be a dual function cytokine.

Identification of genes and pathways involved in retinal neovascularization by microarray analysis of two animal models of retinal angiogenesis.

PURPOSE: Comparative retinal gene expression analysis in two rodent models of oxygen-induced retinopathy (OIR) was performed to identify the genes and pathways involved in retinal neovascularization. METHODS: Three independent experimental runs were conducted for each species, according to standard protocols for induction of OIR. Total retinal RNA was isolated at two time points, corresponding to the early response to relative hypoxia (P13 in mouse, P15 in rat) and to the later phase of maximum retinal neovascularization (P18 in mouse, P20 in rat) and was used to prepare labeled probes for hybridization. Gene expression was compared between normal and experimental conditions for each species at each time point. Probesets with a false-discovery rate of

Androgen regulated genes in human prostate xenografts in mice: relation to BPH and prostate cancer.

Benign prostatic hyperplasia (BPH) and prostate carcinoma (CaP) are linked to aging and the presence of androgens, suggesting that androgen regulated genes play a major role in these common diseases. Androgen regulation of prostate growth and development depends on the presence of intact epithelial-stromal interactions. Further, the prostatic stroma is implicated in BPH. This suggests that epithelial cell lines are inadequate to identify androgen regulated genes that could contribute to BPH and CaP and which could serve as potential clinical biomarkers. In this study, we used a human prostate xenograft model to define a profile of genes regulated in vivo by androgens, with an emphasis on identifying candidate biomarkers. Benign transition zone (TZ) human prostate tissue from radical prostatectomies was grafted to the sub-renal capsule site of intact or castrated male immunodeficient mice, followed by the removal or addition of androgens, respectively. Microarray analysis of RNA from these tissues was used to identify genes that were; 1) highly expressed in prostate, 2) had significant expression changes in response to androgens, and, 3) encode extracellular proteins. A total of 95 genes meeting these criteria were selected for analysis and validation of expression in patient prostate tissues using quantitative real-time PCR. Expression levels of these genes were measured in pooled RNAs from human prostate tissues with varying severity of BPH pathologic changes and CaP of varying Gleason score. A number of androgen regulated genes were identified. Additionally, a subset of these genes were over-expressed in RNA from clinical BPH tissues, and the levels of many were found to correlate with disease status. Our results demonstrate the feasibility, and some of the problems, of using a mouse xenograft model to characterize the androgen regulated expression profiles of intact human prostate tissues.