A targeted sequencing study of glutamatergic candidate genes in suicide attempters with bipolar disorder

Suicidal behavior has been shown to have a heritable component that is partly driven by psychiatric disorders [Brent and Mann, 2005]. However, there is also an independent factor contributing to the heritability of suicidal behavior. We previously conducted a genome-wide association study (GWAS) of bipolar suicide attempters and bipolar non-attempters to assess this independent factor [Willour et al., 2012]. This GWAS implicated glutamatergic neurotransmission in attempted suicide. In the current study, we have conducted a targeted next-generation sequencing study of the glutamatergic N-methyl-D-aspartate (NMDA) receptor, neurexin, and neuroligin gene families in 476 bipolar suicide attempters and 473 bipolar non-attempters. The goal of this study was to gather sequence information from coding and regulatory regions of these glutamatergic genes to identify variants associated with attempted suicide. We identified 186 coding variants and 4,298 regulatory variants predicted to be functional in these genes. No individual variants were overrepresented in cases or controls to a degree that was statistically significant after correction for multiple testing. Additionally, none of the gene-level results were statistically significant following correction. While this study provides no direct support for a role of the examined glutamatergic candidate genes, further sequencing in expanded gene sets and datasets will be required to ultimately determine whether genetic variation in glutamatergic signaling influences suicidal behavior. © 2016 Wiley Periodicals, Inc.

Transposon mutagenesis identifies candidate genes that cooperate with loss of transforming growth factor-beta signaling in mouse intestinal neoplasms

Colorectal cancer (CRC) results from the accumulation of gene mutations and epigenetic alterations in colon epithelial cells, which promotes CRC formation through deregulating signaling pathways. One of the most commonly deregulated signaling pathways in CRC is the transforming growth factor β (TGF-β) pathway. Importantly, the effects of TGF-β signaling inactivation in CRC are modified by concurrent mutations in the tumor cell, and these concurrent mutations determine the ultimate biological effects of impaired TGF-β signaling in the tumor. However, many of the mutations that cooperate with the deregulated TGF-β signaling pathway in CRC remain unknown. Therefore, we sought to identify candidate driver genes that promote the formation of CRC in the setting of TGF-β signaling inactivation. We performed a forward genetic screen in mice carrying conditionally inactivated alleles of the TGF-β receptor, type II (Tgfbr2) using Sleeping Beauty (SB) transposon mediated mutagenesis. We used TAPDANCE and Gene-centric statistical methods to identify common insertion sites (CIS) and, thus, candidate tumor suppressor genes and oncogenes within the tumor genome. CIS analysis of multiple neoplasms from these mice identified many candidate Tgfbr2 cooperating genes and the Wnt/β-catenin, Hippo and MAPK pathways as the most commonly affected pathways. Importantly, the majority of candidate genes were also found to be mutated in human CRC. The SB transposon system provides an unbiased method to identify Tgfbr2 cooperating genes in mouse CRC that are functionally relevant and that may provide further insight into the pathogenesis of human CRC.

Glaucoma Risk Alleles in the Ocular Hypertension Treatment Study

PURPOSE

Primary open-angle glaucoma (POAG) is a major cause of blindness and visual disability. Several genetic risk factors for POAG and optic nerve features have been identified. We measured the relative risk for glaucoma that these factors contribute to participants in the Ocular Hypertension Treatment Study (OHTS).

DESIGN

Comparative series.

PARTICIPANTS

One thousand fifty-seven of 1636 participants (65%) of the OHTS were enrolled in this genetics ancillary study.

METHODS

Samples of DNA were available from 1057 OHTS participants. Of these, 209 developed POAG (cases) and 848 did not develop glaucoma (controls) between 1994 and 2009. The frequencies of 13 risk alleles previously associated with POAG or with optic disc features in other cohorts were compared between POAG cases and controls in the OHTS cohort using analyses of variance. The 2 largest subgroups, non-Hispanic whites (n = 752; 70.7%) and blacks (n = 249, 23.7%), also were analyzed separately. The probability of glaucoma developing over the course of the OHTS was compared between participants stratified for transmembrane and coiled-coil domains 1 (TMCO1) risk alleles using Kaplan-Meier and Cox proportional hazards analyses.

MAIN OUTCOME MEASURES

Association of POAG with known genetic factors.

RESULTS

No association was detected between the known POAG risk alleles when the OHTS cohort was examined as a whole. However, in the subgroup of non-Hispanic whites, allele frequencies at the TMCO1 locus were statistically different between cases and controls (P = 0.00028). By 13 years, non-Hispanic white participants with TMCO1 risk alleles had a 12% higher cumulative frequency of glaucoma developing than participants with no TMCO1 risk alleles. Moreover, the Cox proportional hazard analysis demonstrated that TMCO1 alleles increased relative risk comparable with that of some previously analyzed clinical measures (i.e., intraocular pressure).

CONCLUSIONS

The size of the OHTS cohort and its composition of 2 large racial subgroups may limit its power to detect some glaucoma risk factors. However, TMCO1 genotype was found to increase the risk of glaucoma developing among non-Hispanic whites, the largest racial subgroup in the OHTS cohort, at a magnitude similar to clinical predictors of disease that long have been associated with glaucoma.

Tank-Binding Kinase 1 (TBK1) Gene and Open-Angle Glaucomas (An American Ophthalmological Society Thesis)

PURPOSE

To investigate the role of TANK-binding kinase 1 (TBK1) gene copy-number variations (ie, gene duplications and triplications) in the pathophysiology of various open-angle glaucomas.

METHODS

In previous studies, we discovered that copy-number variations in the TBK1 gene are associated with normal-tension glaucoma. Here, we investigated the prevalence of copy-number variations in cohorts of patients with other open-angle glaucomas-juvenile-onset open-angle glaucoma (n=30), pigmentary glaucoma (n=209), exfoliation glaucoma (n=225), and steroid-induced glaucoma (n=79)-using a quantitative polymerase chain reaction assay.

RESULTS

No TBK1 gene copy-number variations were detected in patients with juvenile-onset open-angle glaucoma, pigmentary glaucoma, or steroid-induced glaucoma. A TBK1 gene duplication was detected in one (0.44%) of the 225 exfoliation glaucoma patients.

CONCLUSIONS

TBK1 gene copy-number variations (gene duplications and triplications) have been previously associated with normal-tension glaucoma. An exploration of other open-angle glaucomas detected a TBK1 copy-number variation in a patient with exfoliation glaucoma, which is the first example of a TBK1 mutation in a glaucoma patient with a diagnosis other than normal-tension glaucoma. A broader phenotypic range may be associated with TBK1 copy-number variations, although mutations in this gene are most often detected in patients with normal-tension glaucoma.

SQSTM1 Mutations and Glaucoma

Glaucoma is the most common cause of irreversible blindness worldwide. One subset of glaucoma, normal tension glaucoma (NTG) occurs in the absence of high intraocular pressure. Mutations in two genes, optineurin (OPTN) and TANK binding kinase 1 (TBK1), cause familial NTG and have known roles in the catabolic cellular process autophagy. TKB1 encodes a kinase that phosphorylates OPTN, an autophagy receptor, which ultimately activates autophagy. The sequestosome (SQSTM1) gene also encodes an autophagy receptor and also is a target of TBK1 phosphorylation. Consequently, we hypothesized that mutations in SQSTM1 may also cause NTG. We tested this hypothesis by searching for glaucoma-causing mutations in a cohort of NTG patients (n = 308) and matched controls (n = 157) using Sanger sequencing. An additional 1098 population control samples were also analyzed using whole exome sequencing. A total of 17 non-synonymous mutations were detected which were not significantly skewed between cases and controls when analyzed separately, or as a group (p > 0.05). These data suggest that SQSTM1 mutations are not a common cause of NTG.

Automated Axon Counting in Rodent Optic Nerve Sections with AxonJ

We have developed a publicly available tool, AxonJ, which quantifies the axons in optic nerve sections of rodents stained with paraphenylenediamine (PPD). In this study, we compare AxonJ's performance to human experts on 100x and 40x images of optic nerve sections obtained from multiple strains of mice, including mice with defects relevant to glaucoma. AxonJ produced reliable axon counts with high sensitivity of 0.959 and high precision of 0.907, high repeatability of 0.95 when compared to a gold-standard of manual assessments and high correlation of 0.882 to the glaucoma damage staging of a previously published dataset. AxonJ allows analyses that are quantitative, consistent, fully-automated, parameter-free, and rapid on whole optic nerve sections at 40x. As a freely available ImageJ plugin that requires no highly specialized equipment to utilize, AxonJ represents a powerful new community resource augmenting studies of the optic nerve using mice.

A Homozygous Nme7 Mutation Is Associated with Situs Inversus Totalis

We investigated the cause of situs inversus totalis (SIT) in two siblings from a consanguineous family. Genotyping and whole-exome analysis revealed a homozygous change in NME7, resulting in deletion of an exon causing an in-frame deletion of 34 amino acids located in the second NDK domain of the protein and segregated with the defective lateralization in the family. NME7 is an important developmental gene, and NME7 protein is a component of the γ-tubulin ring complex. This mutation is predicted to affect the interaction of NME7 protein with this complex as it deletes the amino acids crucial for the binding. SIT associated with homozygous deletion in our patients is in line with Nme7(-/-) mutant mice phenotypes consisting of congenital hydrocephalus and SIT, indicating a novel human laterality patterning role for NME7. Further cases are required to elaborate the full human phenotype associated with NME7 mutations.

A Mutation in LTBP2 Causes Congenital Glaucoma in Domestic Cats (Felis catus)

The glaucomas are a group of diseases characterized by optic nerve damage that together represent a leading cause of blindness in the human population and in domestic animals. Here we report a mutation in LTBP2 that causes primary congenital glaucoma (PCG) in domestic cats. We identified a spontaneous form of PCG in cats and established a breeding colony segregating for PCG consistent with fully penetrant, autosomal recessive inheritance of the trait. Elevated intraocular pressure, globe enlargement and elongated ciliary processes were consistently observed in all affected cats by 8 weeks of age. Varying degrees of optic nerve damage resulted by 6 months of age. Although subtle lens zonular instability was a common feature in this cohort, pronounced ectopia lentis was identified in less than 10% of cats examined. Thus, glaucoma in this pedigree is attributed to histologically confirmed arrest in the early post-natal development of the aqueous humor outflow pathways in the anterior segment of the eyes of affected animals. Using a candidate gene approach, significant linkage was established on cat chromosome B3 (LOD 18.38, θ = 0.00) using tightly linked short tandem repeat (STR) loci to the candidate gene, LTBP2. A 4 base-pair insertion was identified in exon 8 of LTBP2 in affected individuals that generates a frame shift that completely alters the downstream open reading frame and eliminates functional domains. Thus, we describe the first spontaneous and highly penetrant non-rodent model of PCG identifying a valuable animal model for primary glaucoma that closely resembles the human disease, providing valuable insights into mechanisms underlying the disease and a valuable animal model for testing therapies.

Quantitative measurement of retinal ganglion cell populations via histology-based random forest classification

The inner surface of the retina contains a complex mixture of neurons, glia, and vasculature, including retinal ganglion cells (RGCs), the final output neurons of the retina and primary neurons that are damaged in several blinding diseases. The goal of the current work was two-fold: to assess the feasibility of using computer-assisted detection of nuclei and random forest classification to automate the quantification of RGCs in hematoxylin/eosin (H&E)-stained retinal whole-mounts; and if possible, to use the approach to examine how nuclear size influences disease susceptibility among RGC populations. To achieve this, data from RetFM-J, a semi-automated ImageJ-based module that detects, counts, and collects quantitative data on nuclei of H&E-stained whole-mounted retinas, were used in conjunction with a manually curated set of images to train a random forest classifier. To test performance, computer-derived outputs were compared to previously published features of several well-characterized mouse models of ophthalmic disease and their controls: normal C57BL/6J mice; Jun-sufficient and Jun-deficient mice subjected to controlled optic nerve crush (CONC); and DBA/2J mice with naturally occurring glaucoma. The result of these efforts was development of RetFM-Class, a command-line-based tool that uses data output from RetFM-J to perform random forest classification of cell type. Comparative testing revealed that manual and automated classifications by RetFM-Class correlated well, with 83.2% classification accuracy for RGCs. Automated characterization of C57BL/6J retinas predicted 54,642 RGCs per normal retina, and identified a 48.3% Jun-dependent loss of cells at 35 days post CONC and a 71.2% loss of RGCs among 16-month-old DBA/2J mice with glaucoma. Output from automated analyses was used to compare nuclear area among large numbers of RGCs from DBA/2J mice (n = 127,361). In aged DBA/2J mice with glaucoma, RetFM-Class detected a decrease in median and mean nucleus size of cells classified into the RGC category, as did an independent confirmation study using manual measurements of nuclear area demarcated by BRN3A-immunoreactivity. In conclusion, we have demonstrated that histology-based random forest classification is feasible and can be utilized to study RGCs in a high-throughput fashion. Despite having some limitations, this approach demonstrated a significant association between the size of the RGC nucleus and the DBA/2J form of glaucoma.

Molecular response of chorioretinal endothelial cells to complement injury: implications for macular degeneration

Age-related macular degeneration (AMD) is a common, blinding disease of the elderly in which macular photoreceptor cells, retinal pigment epithelium and choriocapillaris endothelial cells ultimately degenerate. Recent studies have found that degeneration of the choriocapillaris occurs early in this disease and that endothelial cell drop-out is concomitant with increased deposition of the complement membrane attack complex (MAC) at the choroidal endothelium. However, the impact of MAC injury to choroidal endothelial cells is poorly understood. To model this event in vitro, and to study the downstream consequences of MAC injury, endothelial cells were exposed to complement from human serum, compared to heat-inactivated serum, which lacks complement components. Cells exposed to complement components in human serum showed increased labelling with antibodies directed against the MAC, time- and dose-dependent cell death, as assessed by lactate dehydrogenase assay and increased permeability. RNA-Seq analysis following complement injury revealed increased expression of genes associated with angiogenesis including matrix metalloproteinase (MMP)-3 and -9, and VEGF-A. The MAC-induced increase in MMP9 RNA expression was validated using C5-depleted serum compared to C5-reconstituted serum. Increased levels of MMP9 were also established, using western blot and zymography. These data suggest that, in addition to cell lysis, complement attack on choroidal endothelial cells promotes an angiogenic phenotype in surviving cells.

Heterozygous triplication of upstream regulatory sequences leads to dysregulation of matrix metalloproteinase 19 in patients with cavitary optic disc anomaly

Patients with a congenital optic nerve disease, cavitary optic disc anomaly (CODA), are born with profound excavation of the optic nerve resembling glaucoma. We previously mapped the gene that causes autosomal-dominant CODA in a large pedigree to a chromosome 12q locus. Using comparative genomic hybridization and quantitative PCR analysis of this pedigree, we report identifying a 6-Kbp heterozygous triplication upstream of the matrix metalloproteinase 19 (MMP19) gene, present in all 17 affected family members and no normal members. Moreover, the triplication was not detected in 78 control subjects or in the Database of Genomic Variants. We further detected the same 6-Kbp triplication in one of 24 unrelated CODA patients and in none of 172 glaucoma patients. Analysis with a Luciferase assay showed that the 6-Kbp sequence has transcription enhancer activity. A 773-bp fragment of the 6-Kbp DNA segment increased downstream gene expression eightfold, suggesting that triplication of this sequence may lead to dysregulation of the downstream gene, MMP19, in CODA patients. Lastly, immunohistochemical analysis of human donor eyes revealed strong expression of MMP19 in optic nerve head. These data strongly suggest that triplication of an enhancer may lead to overexpression of MMP19 in the optic nerve that causes CODA.

Audioprofile Surfaces: The 21st Century Audiogram

OBJECTIVE

To present audiometric data in 3 dimensions by considering age as an addition dimension.

METHODS

Audioprofile surfaces (APSs) were fitted to a set of audiograms by plotting each measurement of an audiogram as an independent point in 3 dimensions with the x, y, and z axes representing frequency, hearing loss in dB, and age, respectively.

RESULTS

Using the Java-based APS viewer as a standalone application, APSs were pre-computed for 34 loci. By selecting APSs for the appropriate genetic locus, a clinician can compare this APS-generated average surface to a specific patient's audiogram.

CONCLUSION

Audioprofile surfaces provide an easily interpreted visual representation of a person's hearing acuity relative to others with the same genetic cause of hearing loss. Audioprofile surfaces will support the generation and testing of sophisticated hypotheses to further refine our understanding of the biology of hearing.

Stereo Photo Measured ONH Shape Predicts Development of POAG in Subjects With Ocular Hypertension

PURPOSE

To identify objective, quantitative optic nerve head (ONH) structural features and model the contributions of glaucoma.

METHODS

Baseline stereoscopic optic disc images of 1635 glaucoma-free participants at risk for developing primary open-angle glaucoma (POAG) were collected as part of the Ocular Hypertension Treatment Study. A stereo correspondence algorithm designed for fundus images was applied to extract the three-dimensional (3D) information about the ONH. Principal component analysis was used to identify ONH 3D structural features and the contributions of demographic features, clinical variables, and disease were modeled using linear regression and linear component analysis. The computationally identified features were evaluated based on associations with glaucoma and ability to predict which participants would develop POAG.

RESULTS

The computationally identified features were significantly associated with future POAG, POAG-related demographics (age, ethnicity), and clinical measurements (horizontal and vertical cup-to-disc ratio, central corneal thickness, and refraction). Models predicting future POAG development using the OHTS baseline data and STEP features achieved an AUC of 0.722 in cross-validation testing. This was a significant improvement over using only demographics (age, sex, and ethnicity), which had an AUC of 0.599.

CONCLUSIONS

Methods for identifying objective, quantitative measurements of 3D ONH structure were developed using a large dataset. The identified features were significantly associated with POAG and POAG-related variables. Further, these features increased predictive model accuracy in predicting future POAG. The results indicate that the computationally identified features might be useful in POAG early screening programs or as endophenotypes to investigate POAG genetics.

North Carolina Macular Dystrophy Is Caused by Dysregulation of the Retinal Transcription Factor PRDM13

PURPOSE

To identify specific mutations causing North Carolina macular dystrophy (NCMD).

DESIGN

Whole-genome sequencing coupled with reverse transcription polymerase chain reaction (RT-PCR) analysis of gene expression in human retinal cells.

PARTICIPANTS

A total of 141 members of 12 families with NCMD and 261 unrelated control individuals.

METHODS

Genome sequencing was performed on 8 affected individuals from 3 families affected with chromosome 6-linked NCMD (MCDR1) and 2 individuals affected with chromosome 5-linked NCMD (MCDR3). Variants observed in the MCDR1 locus with frequencies <1% in published databases were confirmed using Sanger sequencing. Confirmed variants absent from all published databases were sought in 8 additional MCDR1 families and 261 controls. The RT-PCR analysis of selected genes was performed in stem cell-derived human retinal cells.

MAIN OUTCOME MEASURES

Co-segregation of rare genetic variants with disease phenotype.

RESULTS

Five sequenced individuals with MCDR1-linked NCMD shared a haplotype of 14 rare variants spanning 1 Mb of the disease-causing allele. One of these variants (V1) was absent from all published databases and all 261 controls, but was found in 5 additional NCMD kindreds. This variant lies in a DNase 1 hypersensitivity site (DHS) upstream of both the PRDM13 and CCNC genes. Sanger sequencing of 1 kb centered on V1 was performed in the remaining 4 NCMD probands, and 2 additional novel single nucleotide variants (V2 in 3 families and V3 in 1 family) were identified in the DHS within 134 bp of the location of V1. A complete duplication of the PRDM13 gene was also discovered in a single family (V4). The RT-PCR analysis of PRDM13 expression in developing retinal cells revealed marked developmental regulation. Next-generation sequencing of 2 individuals with MCDR3-linked NCMD revealed a 900-kb duplication that included the entire IRX1 gene (V5). The 5 mutations V1 to V5 segregated perfectly in the 102 affected and 39 unaffected members of the 12 NCMD families.

CONCLUSIONS

We identified 5 rare mutations, each capable of arresting human macular development. Four of these strongly implicate the involvement of PRDM13 in macular development, whereas the pathophysiologic mechanism of the fifth remains unknown but may involve the developmental dysregulation of IRX1.

Hypomorphic mutations in TRNT1 cause retinitis pigmentosa with erythrocytic microcytosis

Retinitis pigmentosa (RP) is a highly heterogeneous group of disorders characterized by degeneration of the retinal photoreceptor cells and progressive loss of vision. While hundreds of mutations in more than 100 genes have been reported to cause RP, discovering the causative mutations in many patients remains a significant challenge. Exome sequencing in an individual affected with non-syndromic RP revealed two plausibly disease-causing variants in TRNT1, a gene encoding a nucleotidyltransferase critical for tRNA processing. A total of 727 additional unrelated individuals with molecularly uncharacterized RP were completely screened for TRNT1 coding sequence variants, and a second family was identified with two members who exhibited a phenotype that was remarkably similar to the index patient. Inactivating mutations in TRNT1 have been previously shown to cause a severe congenital syndrome of sideroblastic anemia, B-cell immunodeficiency, recurrent fevers and developmental delay (SIFD). Complete blood counts of all three of our patients revealed red blood cell microcytosis and anisocytosis with only mild anemia. Characterization of TRNT1 in patient-derived cell lines revealed reduced but detectable TRNT1 protein, consistent with partial function. Suppression of trnt1 expression in zebrafish recapitulated several features of the human SIFD syndrome, including anemia and sensory organ defects. When levels of trnt1 were titrated, visual dysfunction was found in the absence of other phenotypes. The visual defects in the trnt1-knockdown zebrafish were ameliorated by the addition of exogenous human TRNT1 RNA. Our findings indicate that hypomorphic TRNT1 mutations can cause a recessive disease that is almost entirely limited to the retina.

Novel TMEM98 mutations in pedigrees with autosomal dominant nanophthalmos

PURPOSE

Autosomal dominant nanophthalmos is an inherited eye disorder characterized by a structurally normal but smaller eye. Patients with nanophthalmos have high hyperopia (far-sightedness), a greater incidence of angle-closure glaucoma, and increased risk of surgical complications. In this study, the clinical features and the genetic basis of nanophthalmos were investigated in two large autosomal dominant nanophthalmos pedigrees.

METHODS

Fourteen members of a Caucasian pedigree from the United States and 15 members of a pedigree from the Mariana Islands enrolled in a genetic study of nanophthalmos and contributed DNA samples. Twenty of 29 family members underwent eye examinations that included measurement of axial eye length and/or refractive error. The genetic basis of nanophthalmos in the pedigrees was studied with linkage analysis, whole exome sequencing, and candidate gene (i.e., TMEM98) sequencing to identify the nanophthalmos-causing gene.

RESULTS

Nine members of the pedigree from the United States and 11 members of the pedigree from the Mariana Islands were diagnosed with nanophthalmos that is transmitted as an autosomal dominant trait. The patients with nanophthalmos had abnormally short axial eye lengths, which ranged from 15.9 to 18.4 mm. Linkage analysis of the nanophthalmos pedigree from the United States identified nine large regions of the genome (greater than 10 Mbp) that were coinherited with disease in this family. Genes within these "linked regions" were examined for disease-causing mutations using exome sequencing, and a His196Pro mutation was detected in the TMEM98 gene, which was recently reported to be a nanophthalmos gene. Sanger sequencing subsequently showed that all other members of this pedigree with nanophthalmos also carry the His196Pro TMEM98 mutation. Testing the Mariana Islands pedigree for TMEM98 mutations identified a 34 bp heterozygous deletion that spans the 3' end of exon 4 in all affected family members. Neither TMEM98 mutation was detected in public exome sequence databases.

CONCLUSIONS

A recent report identified a single TMEM98 missense mutation in a nanophthalmos pedigree. Our discovery of two additional TMEM98 mutations confirms the important role of the gene in the pathogenesis of autosomal dominant nanophthalmos.

RetFM-J, an ImageJ-based module for automated counting and quantifying features of nuclei in retinal whole-mounts

The present article introduces RetFM-J, a semi-automated ImageJ-based module that detects, counts, and collects quantitative data on nuclei of the inner retina from H&E-stained whole-mounted retinas. To illustrate performance, computer-derived outputs were analyzed in inbred C57BL/6J mice. Automated characterization yielded computer-derived outputs that closely matched manual counts. As a method using open-source software that is freely available, inexpensive staining reagents that are robust, and imaging equipment that is routine to most laboratories, RetFM-J could be utilized in a wide variety of experiments benefiting from high-throughput, quantitative, uniform analyses of total cellularity in the inner retina.

Sequencing methods and datasets to improve functional interpretation of sleeping beauty mutagenesis screens

BACKGROUND

Animal models of cancer are useful to generate complementary datasets for comparison to human tumor data. Insertional mutagenesis screens, such as those utilizing the Sleeping Beauty (SB) transposon system, provide a model that recapitulates the spontaneous development and progression of human disease. This approach has been widely used to model a variety of cancers in mice. Comprehensive mutation profiles are generated for individual tumors through amplification of transposon insertion sites followed by high-throughput sequencing. Subsequent statistical analyses identify common insertion sites (CISs), which are predicted to be functionally involved in tumorigenesis. Current methods utilized for SB insertion site analysis have some significant limitations. For one, they do not account for transposon footprints - a class of mutation generated following transposon remobilization. Existing methods also discard quantitative sequence data due to uncertainty regarding the extent to which it accurately reflects mutation abundance within a heterogeneous tumor. Additionally, computational analyses generally assume that all potential insertion sites have an equal probability of being detected under non-selective conditions, an assumption without sufficient relevant data. The goal of our study was to address these potential confounding factors in order to enhance functional interpretation of insertion site data from tumors.

RESULTS

We describe here a novel method to detect footprints generated by transposon remobilization, which revealed minimal evidence of positive selection in tumors. We also present extensive characterization data demonstrating an ability to reproducibly assign semi-quantitative information to individual insertion sites within a tumor sample. Finally, we identify apparent biases for detection of inserted transposons in several genomic regions that may lead to the identification of false positive CISs.

CONCLUSION

The information we provide can be used to refine analyses of data from insertional mutagenesis screens, improving functional interpretation of results and facilitating the identification of genes important in cancer development and progression.

Characterization of Cav1.4 complexes (α11.4, β2, and α2δ4) in HEK293T cells and in the retina

In photoreceptor synaptic terminals, voltage-gated Cav1.4 channels mediate Ca(2+) signals required for transmission of visual stimuli. Like other high voltage-activated Cav channels, Cav1.4 channels are composed of a main pore-forming Cav1.4 α1 subunit and auxiliary β and α2δ subunits. Of the four distinct classes of β and α2δ, β2 and α2δ4 are thought to co-assemble with Cav1.4 α1 subunits in photoreceptors. However, an understanding of the functional properties of this combination of Cav subunits is lacking. Here, we provide evidence that Cav1.4 α1, β2, and α2δ4 contribute to Cav1.4 channel complexes in the retina and describe their properties in electrophysiological recordings. In addition, we identified a variant of β2, named here β2X13, which, along with β2a, is present in photoreceptor terminals. Cav1.4 α1, β2, and α2δ4 were coimmunoprecipitated from lysates of transfected HEK293 cells and mouse retina and were found to interact in the outer plexiform layer of the retina containing the photoreceptor synaptic terminals, by proximity ligation assays. In whole-cell patch clamp recordings of transfected HEK293T cells, channels (Cav1.4 α1 + β2X13) containing α2δ4 exhibited weaker voltage-dependent activation than those with α2δ1. Moreover, compared with channels (Cav1.4 α1 + α2δ4) with β2a, β2X13-containing channels exhibited greater voltage-dependent inactivation. The latter effect was specific to Cav1.4 because it was not seen for Cav1.2 channels. Our results provide the first detailed functional analysis of the Cav1.4 subunits that form native photoreceptor Cav1.4 channels and indicate potential heterogeneity in these channels conferred by β2a and β2X13 variants.

Cordova: web-based management of genetic variation data

Cordova is an out-of-the-box solution for building and maintaining an online database of genetic variations integrated with pathogenicity prediction results from popular algorithms. Our primary motivation for developing this system is to aid researchers and clinician-scientists in determining the clinical significance of genetic variations. To achieve this goal, Cordova provides an interface to review and manually or computationally curate genetic variation data as well as share it for clinical diagnostics and the advancement of research.

AVAILABILITY AND IMPLEMENTATION

Cordova is open source under the MIT license and is freely available for download at https://github.com/clcg/cordova.

An international effort towards developing standards for best practices in analysis, interpretation and reporting of clinical genome sequencing results in the CLARITY Challenge

Background

There is tremendous potential for genome sequencing to improve clinical diagnosis and care once it becomes routinely accessible, but this will require formalizing research methods into clinical best practices in the areas of sequence data generation, analysis, interpretation and reporting. The CLARITY Challenge was designed to spur convergence in methods for diagnosing genetic disease starting from clinical case history and genome sequencing data. DNA samples were obtained from three families with heritable genetic disorders and genomic sequence data were donated by sequencing platform vendors. The challenge was to analyze and interpret these data with the goals of identifying disease-causing variants and reporting the findings in a clinically useful format. Participating contestant groups were solicited broadly, and an independent panel of judges evaluated their performance.

Results

A total of 30 international groups were engaged. The entries reveal a general convergence of practices on most elements of the analysis and interpretation process. However, even given this commonality of approach, only two groups identified the consensus candidate variants in all disease cases, demonstrating a need for consistent fine-tuning of the generally accepted methods. There was greater diversity of the final clinical report content and in the patient consenting process, demonstrating that these areas require additional exploration and standardization.

Conclusions

The CLARITY Challenge provides a comprehensive assessment of current practices for using genome sequencing to diagnose and report genetic diseases. There is remarkable convergence in bioinformatic techniques, but medical interpretation and reporting are areas that require further development by many groups.

Altered gene expression in dry age-related macular degeneration suggests early loss of choroidal endothelial cells

PURPOSE

Age-related macular degeneration (AMD) is a major cause of blindness in developed countries. The molecular pathogenesis of early events in AMD is poorly understood. We investigated differential gene expression in samples of human retinal pigment epithelium (RPE) and choroid from early AMD and control maculas with exon-based arrays.

METHODS

Gene expression levels in nine human donor eyes with early AMD and nine control human donor eyes were assessed using Affymetrix Human Exon ST 1.0 arrays. Two controls did not pass quality control and were removed. Differentially expressed genes were annotated using the Database for Annotation, Visualization and Integrated Discovery (DAVID), and gene set enrichment analysis (GSEA) was performed on RPE-specific and endothelium-associated gene sets. The complement factor H (CFH) genotype was also assessed, and differential expression was analyzed regarding high AMD risk (YH/HH) and low AMD risk (YY) genotypes.

RESULTS

Seventy-five genes were identified as differentially expressed (raw p value <0.01; ≥50% fold change, mean log2 expression level in AMD or control ≥ median of all average gene expression values); however, no genes were significant (adj. p value <0.01) after correction for multiple hypothesis testing. Of 52 genes with decreased expression in AMD (fold change <0.5; raw p value <0.01), 18 genes were identified by DAVID analysis as associated with vision or neurologic processes. The GSEA of the RPE-associated and endothelium-associated genes revealed a significant decrease in genes typically expressed by endothelial cells in the early AMD group compared to controls, consistent with previous histologic and proteomic studies. Analysis of the CFH genotype indicated decreased expression of ADAMTS9 in eyes with high-risk genotypes (fold change = -2.61; raw p value=0.0008).

CONCLUSIONS

GSEA results suggest that RPE transcripts are preserved or elevated in early AMD, concomitant with loss of endothelial cell marker expression. These results are consistent with the notion that choroidal endothelial cell dropout or dedifferentiation occurs early in the pathogenesis of AMD.

Advancing genetic testing for deafness with genomic technology

BACKGROUND

Non-syndromic hearing loss (NSHL) is the most common sensory impairment in humans. Until recently its extreme genetic heterogeneity precluded comprehensive genetic testing. Using a platform that couples targeted genomic enrichment (TGE) and massively parallel sequencing (MPS) to sequence all exons of all genes implicated in NSHL, we tested 100 persons with presumed genetic NSHL and in so doing established sequencing requirements for maximum sensitivity and defined MPS quality score metrics that obviate Sanger validation of variants.

METHODS

We examined DNA from 100 sequentially collected probands with presumed genetic NSHL without exclusions due to inheritance, previous genetic testing, or type of hearing loss. We performed TGE using post-capture multiplexing in variable pool sizes followed by Illumina sequencing. We developed a local Galaxy installation on a high performance computing cluster for bioinformatics analysis.

RESULTS

To obtain maximum variant sensitivity with this platform 3.2-6.3 million total mapped sequencing reads per sample were required. Quality score analysis showed that Sanger validation was not required for 95% of variants. Our overall diagnostic rate was 42%, but this varied by clinical features from 0% for persons with asymmetric hearing loss to 56% for persons with bilateral autosomal recessive NSHL.

CONCLUSIONS

These findings will direct the use of TGE and MPS strategies for genetic diagnosis for NSHL. Our diagnostic rate highlights the need for further research on genetic deafness focused on novel gene identification and an improved understanding of the role of non-exonic mutations. The unsolved families we have identified provide a valuable resource to address these areas.