A 3-year-old female patient presenting with an unknown syndrome of a neonatal progeroid appearance, lipodystrophy, pulmonary hypertension, cutis marmorata, feeding disorder and failure to thrive was investigated by whole-genome sequencing. This revealed a de novo, heterozygous, frame-shift mutation in the Caveolin1 gene (CAV1) (p.Phe160X). Mutations in CAV1, encoding the main component of the caveolae in plasma membranes, cause Berardinelli-Seip congenital lipodystrophy type 3 (BSCL). Although BSCL is recessive, heterozygous carriers either show a reduced phenotype of partial lipodystrophy, pulmonary hypertension, or no phenotype. To investigate the pathogenic mechanisms underlying this syndrome in more depth, we performed next generation RNA sequencing of peripheral blood, which showed several dysregulated pathways in the patient that might be related to the phenotypic progeroid features (apoptosis, DNA repair/replication, mitochondrial). Secondly, we found a significant down-regulation of known Cav1 interaction partners, verifying the dysfunction of CAV1. Other known progeroid genes and lipodystrophy genes were also dysregulated. Next, western blotting of lysates of cultured fibroblasts showed that the patient shows a significantly decreased expression of wild-type CAV1 protein, demonstrating a loss-of-function mutation, though her phenotype is more severe that other heterozygotes with similar mutations. This phenotypic variety could be explained by differences in genetic background. Indications for this are supported by additional rare variants we found in AGPAT2 and LPIN1 lipodystrophy genes. CAV1, AGPAT2 and LPIN1 all play an important role in triacylglycerol (TAG) biosynthesis in adipose tissue, and the defective function in different parts of this pathway, though not all to the same extend, could contribute to a more severe lipoatrophic phenotype in this patient. In conclusion, we report, for the first time, an association of CAV1 dysfunction with a syndrome of severe premature aging and lipodystrophy. This may contribute to a better understanding of the aging process and pathogenic mechanisms that contribute to premature aging.

Matching molecularly targeted therapies with cancer subtype-specific gene mutations is revolutionizing oncology care. However, for rare cancers this approach is problematic due to the often poor understanding of the disease's natural history and phenotypic heterogeneity, making treatment of these cancers a particularly unmet medical need in clinical oncology. Advanced Sézary syndrome (SS), an aggressive, exceedingly rare variant of cutaneous T-cell lymphoma (CTCL) is a prototypical example of a rare cancer. Through whole genome and RNA sequencing (RNA-seq) of a SS patient's tumor we discovered a highly expressed gene fusion between CTLA4 (cytotoxic T lymphocyte antigen 4) and CD28 (cluster of differentiation 28), predicting a novel stimulatory molecule on the surface of tumor T cells. Treatment with the CTLA4 inhibitor ipilimumab resulted in a rapid clinical response. Our findings suggest a novel driver mechanism for SS, and cancer in general, and exemplify an emerging model of cancer treatment using exploratory genomic analysis to identify a personally targeted treatment option when conventional therapies are exhausted.

Alzheimer's disease (AD) is characterized by deficits in cerebral metabolic rates of glucose in the posterior cingulate (PC) and precuneus in AD subjects, and in APOEε4 carriers, decades before the onset of measureable cognitive deficits. However, the cellular and molecular basis of this phenotype remains to be clarified. Given the roles of astrocytes in energy storage and brain immunity, we sought to characterize the transcriptome of AD PC astrocytes. Cells were laser capture microdissected from AD (n = 10) and healthy elderly control (n = 10) subjects for RNA sequencing. We generated >5.22 billion reads and compared sequencing data between controls and AD patients. We identified differentially expressed mitochondria-related genes including TRMT61B, FASTKD2, and NDUFA4L2, and using pathway and weighted gene coexpression analyses, we identified differentially expressed immune response genes. A number of these genes, including CLU, C3, and CD74, have been implicated in beta amyloid generation or clearance. These data provide key insights into astrocyte-specific contributions to AD, and we present this data set as a publicly available resource.

To determine the phenotypic spectrum caused by mutations in GRIN1 encoding the NMDA receptor subunit GluN1 and to investigate their underlying functional pathophysiology.

CaBPs are a family of Ca(2+)-binding proteins related to calmodulin and are localized in the brain and sensory organs, including the retina and cochlea. Although their physiological roles are not yet fully elucidated, CaBPs modulate Ca(2+) signaling through effectors such as voltage-gated Ca(v) Ca(2+) channels. In this study, we identified a splice-site mutation (c.637+1G>T) in Ca(2+)-binding protein 2 (CABP2) in three consanguineous Iranian families affected by moderate-to-severe hearing loss. This mutation, most likely a founder mutation, probably leads to skipping of exon 6 and premature truncation of the protein (p.Phe164Serfs(∗)4). Compared with wild-type CaBP2, the truncated CaBP2 showed altered Ca(2+) binding in isothermal titration calorimetry and less potent regulation of Ca(v)1.3 Ca(2+) channels. We show that genetic defects in CABP2 cause moderate-to-severe sensorineural hearing impairment. The mutation might cause a hypofunctional CaBP2 defective in Ca(2+) sensing and effector regulation in the inner ear.

The pathogenesis of hypothalamic hamartoma (HH) associated with epilepsy is unknown. We have identified an individual with HH and refractory epilepsy exhibiting subtle dysmorphic features. High-resolution karyotype identified a duplication of the terminal end of 6p (6p25.1-25.3), confirmed by fluorescent in situ-hybridization (FISH). Copy number analysis with high-density (250K) single nucleotide polymorphism (SNP) genotyping microarrays characterized the abnormality as a series of amplified regions between 1.4 Mb and 10.2 Mb, with a small tandem deletion from 8.8 Mb to 9.7 Mb. There are 38 RefSeq genes within the duplicated regions, and no known coding sequences within the deletion. This unique patient helps identify 6p25.1-25.3 as a possible susceptibility locus for sporadic HH.

We report the development and validation of experimental methods, study designs, and analysis software for pooling-based genomewide association (GWA) studies that use high-throughput single-nucleotide-polymorphism (SNP) genotyping microarrays. We first describe a theoretical framework for establishing the effectiveness of pooling genomic DNA as a low-cost alternative to individually genotyping thousands of samples on high-density SNP microarrays. Next, we describe software called "GenePool," which directly analyzes SNP microarray probe intensity data and ranks SNPs by increased likelihood of being genetically associated with a trait or disorder. Finally, we apply these methods to experimental case-control data and demonstrate successful identification of published genetic susceptibility loci for a rare monogenic disease (sudden infant death with dysgenesis of the testes syndrome), a rare complex disease (progressive supranuclear palsy), and a common complex disease (Alzheimer disease) across multiple SNP genotyping platforms. On the basis of these theoretical calculations and their experimental validation, our results suggest that pooling-based GWA studies are a logical first step for determining whether major genetic associations exist in diseases with high heritability.

Hearing impairment (HI) is a common disorder of sensorineural function with a highly heterogeneous genetic background. Although substantial progress has been made in the understanding of the genetic etiology of hereditary HI, many genes implicated in HI remain undiscovered. Via exome and Sanger sequencing of DNA samples obtained from consanguineous Pakistani families that segregate profound prelingual sensorineural HI, we identified rare homozygous missense variants in four genes (ADAMTS1, MPDZ, MVD, and SEZ6) that are likely the underlying cause of HI. Linkage analysis provided statistical evidence that these variants are associated with autosomal recessive nonsyndromic HI. In silico analysis of the mutant proteins encoded by these genes predicted structural, conformational or interaction changes. RNAseq data analysis revealed expression of these genes in the sensory epithelium of the mouse inner ear during embryonic, postnatal, and adult stages. Immunohistochemistry of the mouse cochlear tissue, further confirmed the expression of ADAMTS1, SEZ6, and MPDZ in the neurosensory hair cells of the organ of Corti, while MVD expression was more prominent in the spiral ganglion cells. Overall, supported by in silico mutant protein analysis, animal models, linkage analysis, and spatiotemporal expression profiling in the mouse inner ear, we propose four new candidate genes for HI and expand our understanding of the etiology of HI.

Hepatocyte growth factor-overexpressing mice that harbor a deletion of the Ink4a/p16 locus (HP mice) form melanomas with low metastatic potential in response to UV irradiation. Here we report that these tumors become highly metastatic following hemizygous deletion of the Nme1 and Nme2 metastasis suppressor genes (HPN mice). Whole-genome sequencing of melanomas from HPN mice revealed a striking increase in lung metastatic activity that is associated with missense mutations in eight signature genes (Arhgap35, Atp8b4, Brca1, Ift172, Kif21b, Nckap5, Pcdha2, and Zfp869). RNA-seq analysis of transcriptomes from HP and HPN primary melanomas identified a 32-gene signature (HPN lung metastasis signature) for which decreased expression is strongly associated with lung metastatic potential. Analysis of transcriptome data from The Cancer Genome Atlas revealed expression profiles of these genes that predict improved survival of patients with cutaneous or uveal melanoma. Silencing of three representative HPN lung metastasis signature genes (ARRDC3, NYNRIN, RND3) in human melanoma cells resulted in increased invasive activity, consistent with roles for these genes as mediators of the metastasis suppressor function of NME1 and NME2. In conclusion, our studies have identified a family of genes that mediate suppression of melanoma lung metastasis, and which may serve as prognostic markers and/or therapeutic targets for clinical management of metastatic melanoma.

Although variant alleles of hundreds of genes are associated with sensorineural deafness in children, the genes and alleles involved remain largely unknown in the Sub-Saharan regions of Africa. We ascertained 56 small families mainly of Yoruba ethno-lingual ancestry in or near Ibadan, Nigeria, that had at least one individual with nonsyndromic, severe-to-profound, prelingual-onset, bilateral hearing loss not attributed to nongenetic factors. We performed a combination of exome and Sanger sequencing analyses to evaluate both nuclear and mitochondrial genomes. No biallelic pathogenic variants were identified in GJB2, a common cause of deafness in many populations. Potential causative variants were identified in genes associated with nonsyndromic hearing loss (CIB2, COL11A1, ILDR1, MYO15A, TMPRSS3, and WFS1), nonsyndromic hearing loss or Usher syndrome (CDH23, MYO7A, PCDH15, and USH2A), and other syndromic forms of hearing loss (CHD7, OPA1, and SPTLC1). Several rare mitochondrial variants, including m.1555A>G, were detected in the gene MT-RNR1 but not in control Yoruba samples. Overall, 20 (33%) of 60 independent cases of hearing loss in this cohort of families were associated with likely causal variants in genes reported to underlie deafness in other populations. None of these likely causal variants were present in more than one family, most were detected as compound heterozygotes, and 77% had not been previously associated with hearing loss. These results indicate an unusually high level of genetic heterogeneity of hearing loss in Ibadan, Nigeria and point to challenges for molecular genetic screening, counseling, and early intervention in this population.

Cardiomyopathies are diseases of the heart muscle and are important causes of heart failure. Dilated cardiomyopathy (DCM) is a common form of cardiomyopathy that can be acquired, syndromic or non-syndromic. The current study was conducted to explore the genetic defects in a Pakistani family with cardiac disease and features of Marfan's syndrome (MFS).

Personalized medicine promises individualized disease prediction and treatment. The convergence of machine learning (ML) and available multimodal data is key moving forward. We build upon previous work to deliver multimodal predictions of Parkinson's disease (PD) risk and systematically develop a model using GenoML, an automated ML package, to make improved multi-omic predictions of PD, validated in an external cohort. We investigated top features, constructed hypothesis-free disease-relevant networks, and investigated drug-gene interactions. We performed automated ML on multimodal data from the Parkinson's progression marker initiative (PPMI). After selecting the best performing algorithm, all PPMI data was used to tune the selected model. The model was validated in the Parkinson's Disease Biomarker Program (PDBP) dataset. Our initial model showed an area under the curve (AUC) of 89.72% for the diagnosis of PD. The tuned model was then tested for validation on external data (PDBP, AUC 85.03%). Optimizing thresholds for classification increased the diagnosis prediction accuracy and other metrics. Finally, networks were built to identify gene communities specific to PD. Combining data modalities outperforms the single biomarker paradigm. UPSIT and PRS contributed most to the predictive power of the model, but the accuracy of these are supplemented by many smaller effect transcripts and risk SNPs. Our model is best suited to identifying large groups of individuals to monitor within a health registry or biobank to prioritize for further testing. This approach allows complex predictive models to be reproducible and accessible to the community, with the package, code, and results publicly available.

Orofaciodigital syndrome (OFD) is clinically heterogeneous and is characterized by abnormalities in the oral cavity, facial features, digits, and central nervous system. At least 18 subtypes of the condition have been described in the literature. OFD is caused by variants in several genes with overlapping phenotypes. We studied a consanguineous Pakistani family with two affected siblings with an atypical form of OFD type 4 (OFD4). In addition to the typical features of OFD4 that include limb defects and growth retardation, the siblings displayed rare features of scaphocephaly and seizures. Exome sequencing analysis revealed a novel homozygous splice site variant c.257-1G>A in TCTN3 that segregated with disease. This homozygous splice site variant in TCTN3 is most likely the underlying cause of the atypical form of OFD4 observed in this family. Our results contribute to the phenotypic spectrum of TCTN3 associated ciliopathies and will facilitate better clinical diagnosis.

The Foundational Data Initiative for Parkinson Disease (FOUNDIN-PD) is an international collaboration producing fundamental resources for Parkinson disease (PD). FOUNDIN-PD generated a multi-layered molecular dataset in a cohort of induced pluripotent stem cell (iPSC) lines differentiated to dopaminergic (DA) neurons, a major affected cell type in PD. The lines were derived from the Parkinson's Progression Markers Initiative study, which included participants with PD carrying monogenic PD variants, variants with intermediate effects, and variants identified by genome-wide association studies and unaffected individuals. We generated genetic, epigenetic, regulatory, transcriptomic, and longitudinal cellular imaging data from iPSC-derived DA neurons to understand molecular relationships between disease-associated genetic variation and proximate molecular events. These data reveal that iPSC-derived DA neurons provide a valuable cellular context and foundational atlas for modeling PD genetic risk. We have integrated these data into a FOUNDIN-PD data browser as a resource for understanding the molecular pathogenesis of PD.

A genome-wide, whole brain approach to investigate genetic effects on neuroimaging phenotypes for identifying quantitative trait loci is described. The Alzheimer's Disease Neuroimaging Initiative 1.5 T MRI and genetic dataset was investigated using voxel-based morphometry (VBM) and FreeSurfer parcellation followed by genome-wide association studies (GWAS). One hundred forty-two measures of grey matter (GM) density, volume, and cortical thickness were extracted from baseline scans. GWAS, using PLINK, were performed on each phenotype using quality-controlled genotype and scan data including 530,992 of 620,903 single nucleotide polymorphisms (SNPs) and 733 of 818 participants (175 AD, 354 amnestic mild cognitive impairment, MCI, and 204 healthy controls, HC). Hierarchical clustering and heat maps were used to analyze the GWAS results and associations are reported at two significance thresholds (p<10(-7) and p<10(-6)). As expected, SNPs in the APOE and TOMM40 genes were confirmed as markers strongly associated with multiple brain regions. Other top SNPs were proximal to the EPHA4, TP63 and NXPH1 genes. Detailed image analyses of rs6463843 (flanking NXPH1) revealed reduced global and regional GM density across diagnostic groups in TT relative to GG homozygotes. Interaction analysis indicated that AD patients homozygous for the T allele showed differential vulnerability to right hippocampal GM density loss. NXPH1 codes for a protein implicated in promotion of adhesion between dendrites and axons, a key factor in synaptic integrity, the loss of which is a hallmark of AD. A genome-wide, whole brain search strategy has the potential to reveal novel candidate genes and loci warranting further investigation and replication.

Recent advances in the treatment of cancer have focused on targeting genomic aberrations with selective therapeutic agents. In rare tumors, where large-scale clinical trials are daunting, this targeted genomic approach offers a new perspective and hope for improved treatments. Cancers of the ampulla of Vater are rare tumors that comprise only about 0.2% of gastrointestinal cancers. Consequently, they are often treated as either distal common bile duct or pancreatic cancers.

We performed exome sequencing to evaluate the underlying molecular cause of a patient with bilateral conductive hearing loss due to multiple ossicular abnormalities as well as symphalangism of the fifth digits. This leads to the identification of a novel heterozygous start codon variant in the NOG gene (c.2T>C:p.Met1?) that hinders normal translation of the noggin protein. Variants in NOG lead to a spectrum of otologic, digit, and joint abnormalities, a combination suggested to be referred to as NOG-related-symphalangism spectrum disorder (NOG-SSD). Conductive hearing loss from such variants may stem from stapes footplate ankylosis, fixation of the malleoincudal joint, or fixation of the incus short process. In this case, the constellation of both stapes and incus fixation, an exceptionally tall stapes suprastructure, thickened long process of the incus, and enlarged incus body was encountered, leading to distinct challenges during otologic surgery to improve hearing thresholds. This case highlights multiple abnormalities to the ossicular chain in a patient with a start codon variant in NOG. We provide detailed imaging data on these malformations as well as surgical considerations and outcomes.

Bipolar disorder is a highly heritable psychiatric disorder. We performed a genome-wide association study (GWAS) including 20,352 cases and 31,358 controls of European descent, with follow-up analysis of 822 variants with P < 1 × 10-4 in an additional 9,412 cases and 137,760 controls. Eight of the 19 variants that were genome-wide significant (P < 5 × 10-8) in the discovery GWAS were not genome-wide significant in the combined analysis, consistent with small effect sizes and limited power but also with genetic heterogeneity. In the combined analysis, 30 loci were genome-wide significant, including 20 newly identified loci. The significant loci contain genes encoding ion channels, neurotransmitter transporters and synaptic components. Pathway analysis revealed nine significantly enriched gene sets, including regulation of insulin secretion and endocannabinoid signaling. Bipolar I disorder is strongly genetically correlated with schizophrenia, driven by psychosis, whereas bipolar II disorder is more strongly correlated with major depressive disorder. These findings address key clinical questions and provide potential biological mechanisms for bipolar disorder.

We performed shallow single-cell sequencing of genomic DNA across 1475 cells from a cell-line, COLO829, to resolve overall complexity and clonality. This melanoma tumor-line has been previously characterized by multiple technologies and is a benchmark for evaluating somatic alterations. In some of these studies, COLO829 has shown conflicting and/or indeterminate copy number and, thus, single-cell sequencing provides a tool for gaining insight. Following shallow single-cell sequencing, we first identified at least four major sub-clones by discriminant analysis of principal components of single-cell copy number data. Based on clustering, break-point and loss of heterozygosity analysis of aggregated data from sub-clones, we identified distinct hallmark events that were validated within bulk sequencing and spectral karyotyping. In summary, COLO829 exhibits a classical Dutrillaux's monosomic/trisomic pattern of karyotype evolution with endoreduplication, where consistent sub-clones emerge from the loss/gain of abnormal chromosomes. Overall, our results demonstrate how shallow copy number profiling can uncover hidden biological insights.

Congenital hearing impairment (HI) is genetically heterogeneous making its genetic diagnosis challenging. Investigation of novel HI genes and variants will enhance our understanding of the molecular mechanisms and to aid genetic diagnosis. We performed exome sequencing and analysis using DNA samples from affected members of two large families from Ghana and Pakistan, segregating autosomal-dominant (AD) non-syndromic HI (NSHI). Using in silico approaches, we modeled and evaluated the effect of the likely pathogenic variants on protein structure and function. We identified two likely pathogenic variants in SLC12A2, c.2935G>A:p.(E979K) and c.2939A>T:p.(E980V), which segregate with NSHI in a Ghanaian and Pakistani family, respectively. SLC12A2 encodes an ion transporter crucial in the homeostasis of the inner ear endolymph and has recently been reported to be implicated in syndromic and non-syndromic HI. Both variants were mapped to alternatively spliced exon 21 of the SLC12A2 gene. Exon 21 encodes for 17 residues in the cytoplasmatic tail of SLC12A2, is highly conserved between species, and preferentially expressed in cochlear tissues. A review of previous studies and our current data showed that out of ten families with either AD non-syndromic or syndromic HI, eight (80%) had variants within the 17 amino acid residue region of exon 21 (48 bp), suggesting that this alternate domain is critical to the transporter activity in the inner ear. The genotypic spectrum of SLC12A2 was expanded and the involvement of SLC12A2 in ADNSHI was confirmed. These results also demonstrate the role that SLC12A2 plays in ADNSHI in diverse populations including sub-Saharan Africans.