Alu element in the RNA binding motif protein, X-linked 2 (RBMX2) gene found to be linked to bipolar disorder.

OBJECTIVE: We have used long-read single molecule, real-time (SMRT) sequencing to fully characterize a ~12Mb genomic region on chromosome Xq24-q27, significantly linked to bipolar disorder (BD) in an extended family from a genetic sub-isolate. This family segregates BD in at least four generations with 24 affected individuals. METHODS: We selected 16 family members for targeted sequencing. The selected individuals either carried the disease haplotype, were non-carriers of the disease haplotype, or served as married-in controls. We designed hybrid capture probes enriching for 5-9Kb fragments spanning the entire 12Mb region that were then sequenced to screen for candidate structural variants (SVs) that could explain the increased risk for BD in this extended family. RESULTS: Altogether, 201 variants were detected in the critically linked region. Although most of these represented common variants, three variants emerged that showed near-perfect segregation among all BD type I affected individuals. Two of the SVs were identified in or near genes belonging to the RNA Binding Motif Protein, X-Linked (RBMX) gene family-a 330bp Alu (subfamily AluYa5) deletion in intron 3 of the RBMX2 gene and an intergenic 27bp tandem repeat deletion between the RBMX and G protein-coupled receptor 101 (GPR101) genes. The third SV was a 50bp tandem repeat insertion in intron 1 of the Coagulation Factor IX (F9) gene. CONCLUSIONS: Among the three genetically linked SVs, additional evidence supported the Alu element deletion in RBMX2 as the leading candidate for contributing directly to the disease development of BD type I in this extended family.

Genomic targets for high-resolution inference of kinship, ancestry and disease susceptibility in orang-utans (genus: Pongo).

BACKGROUND: Orang-utans comprise three critically endangered species endemic to the islands of Borneo and Sumatra. Though whole-genome sequencing has recently accelerated our understanding of their evolutionary history, the costs of implementing routine genome screening and diagnostics remain prohibitive. Capitalizing on a tri-fold locus discovery approach, combining data from published whole-genome sequences, novel whole-exome sequencing, and microarray-derived genotype data, we aimed to develop a highly informative gene-focused panel of targets that can be used to address a broad range of research questions. RESULTS: We identified and present genomic co-ordinates for 175,186 SNPs and 2315 Y-chromosomal targets, plus 185 genes either known or presumed to be pathogenic in cardiovascular (N = 109) or respiratory (N = 43) diseases in humans - the primary and secondary causes of captive orang-utan mortality - or a majority of other human diseases (N = 33). As proof of concept, we designed and synthesized 'SeqCap' hybrid capture probes for these targets, demonstrating cost-effective target enrichment and reduced-representation sequencing. CONCLUSIONS: Our targets are of broad utility in studies of orang-utan ancestry, admixture and disease susceptibility and aetiology, and thus are of value in addressing questions key to the survival of these species. To facilitate comparative analyses, these targets could now be standardized for future orang-utan population genomic studies. The targets are broadly compatible with commercial target enrichment platforms and can be utilized as published here to synthesize applicable probes.

Representative Sequencing: Unbiased Sampling of Solid Tumor Tissue.

Although thousands of solid tumors have been sequenced to date, a fundamental under-sampling bias is inherent in current methodologies. This is caused by a tissue sample input of fixed dimensions (e.g., 6 mm biopsy), which becomes grossly under-powered as tumor volume scales. Here, we demonstrate representative sequencing (Rep-Seq) as a new method to achieve unbiased tumor tissue sampling. Rep-Seq uses fixed residual tumor material, which is homogenized and subjected to next-generation sequencing. Analysis of intratumor tumor mutation burden (TMB) variability shows a high level of misclassification using current single-biopsy methods, with 20% of lung and 52% of bladder tumors having at least one biopsy with high TMB but low clonal TMB overall. Misclassification rates by contrast are reduced to 2% (lung) and 4% (bladder) when a more representative sampling methodology is used. Rep-Seq offers an improved sampling protocol for tumor profiling, with significant potential for improved clinical utility and more accurate deconvolution of clonal structure.

Structural variation and its potential impact on genome instability: Novel discoveries in the EGFR landscape by long-read sequencing.

Structural variation (SV) is typically defined as variation within the human genome that exceeds 50 base pairs (bp). SV may be copy number neutral or it may involve duplications, deletions, and complex rearrangements. Recent studies have shown SV to be associated with many human diseases. However, studies of SV have been challenging due to technological constraints. With the advent of third generation (long-read) sequencing technology, exploration of longer stretches of DNA not easily examined previously has been made possible. In the present study, we utilized third generation (long-read) sequencing techniques to examine SV in the EGFR landscape of four haplotypes derived from two human samples. We analyzed the EGFR gene and its landscape (+/- 500,000 base pairs) using this approach and were able to identify a region of non-coding DNA with over 90% similarity to the most common activating EGFR mutation in non-small cell lung cancer. Based on previously published Alu-element genome instability algorithms, we propose a molecular mechanism to explain how this non-coding region of DNA may be interacting with and impacting the stability of the EGFR gene and potentially generating this cancer-driver gene. By these techniques, we were also able to identify previously hidden structural variation in the four haplotypes and in the human reference genome (hg38). We applied previously published algorithms to compare the relative stabilities of these five different EGFR gene landscape haplotypes to estimate their relative potentials to generate the EGFR exon 19, 15 bp canonical deletion. To our knowledge, the present study is the first to use the differences in genomic architecture between targeted cancer-linked phased haplotypes to estimate their relative potentials to form a common cancer-linked driver mutation.

Tracking Cancer Evolution Reveals Constrained Routes to Metastases: TRACERx Renal.

Clear-cell renal cell carcinoma (ccRCC) exhibits a broad range of metastatic phenotypes that have not been systematically studied to date. Here, we analyzed 575 primary and 335 metastatic biopsies across 100 patients with metastatic ccRCC, including two cases sampledat post-mortem. Metastatic competence was afforded by chromosome complexity, and we identify 9p loss as a highly selected event driving metastasis and ccRCC-related mortality (p = 0.0014). Distinct patterns of metastatic dissemination were observed, including rapid progression to multiple tissue sites seeded by primary tumors of monoclonal structure. By contrast, we observed attenuated progression in cases characterized by high primary tumor heterogeneity, with metastatic competence acquired gradually and initial progression to solitary metastasis. Finally, we observed early divergence of primitive ancestral clones and protracted latency of up to two decades as a feature of pancreatic metastases.

Targeted Bisulfite Sequencing Using the SeqCap Epi Enrichment System.

Cytosine methylation has been shown to have a role in a host of biological processes. In mammalian biology these include stem cell differentiation, embryonic development, genomic imprinting, inflammation, and silencing of transposable elements. Given the central importance of these processes, it is not surprising to find aberrant cytosine methylation patterns associated with many disorders in humans, including cancer, cardiovascular disease, and neurological disease. While whole genome shotgun bisulfite sequencing (WGBS) has recently become feasible, generating high sequence coverage data for the entire genome is expensive, both in terms of money and analysis time, when generally only a small subset of the genome is of interest to most researchers. This report details a procedure for the targeted enrichment of bisulfite treated DNA via SeqCap Epi, allowing high resolution focus of next generation sequencing onto a subset of the genome for high resolution cytosine methylation analysis. Regions ranging in size from only a few kb up to over 200 Mb may be targeted, including the use of the SeqCap Epi CpGiant design which is designed to target 5.5 million CpGs in the human genome. Finally, multiple samples may be multiplexed and sequenced together to provide an inexpensive method of generating methylation data for a large number of samples in a high throughput fashion.

LINE-1 Mediated Insertion into Poc1a (Protein of Centriole 1 A) Causes Growth Insufficiency and Male Infertility in Mice.

Skeletal dysplasias are a common, genetically heterogeneous cause of short stature that can result from disruptions in many cellular processes. We report the identification of the lesion responsible for skeletal dysplasia and male infertility in the spontaneous, recessive mouse mutant chagun. We determined that Poc1a, encoding protein of the centriole 1a, is disrupted by the insertion of a processed Cenpw cDNA, which is flanked by target site duplications, suggestive of a LINE-1 retrotransposon-mediated event. Mutant fibroblasts have impaired cilia formation and multipolar spindles. Male infertility is caused by defective spermatogenesis early in meiosis and progressive germ cell loss. Spermatogonial stem cell transplantation studies revealed that Poc1a is essential for normal function of both Sertoli cells and germ cells. The proliferative zone of the growth plate is small and disorganized because chondrocytes fail to re-align after cell division and undergo increased apoptosis. Poc1a and several other genes associated with centrosome function can affect the skeleton and lead to skeletal dysplasias and primordial dwarfisms. This mouse mutant reveals how centrosome dysfunction contributes to defects in skeletal growth and male infertility.

Genetic perturbation of the maize methylome.

DNA methylation can play important roles in the regulation of transposable elements and genes. A collection of mutant alleles for 11 maize (Zea mays) genes predicted to play roles in controlling DNA methylation were isolated through forward- or reverse-genetic approaches. Low-coverage whole-genome bisulfite sequencing and high-coverage sequence-capture bisulfite sequencing were applied to mutant lines to determine context- and locus-specific effects of these mutations on DNA methylation profiles. Plants containing mutant alleles for components of the RNA-directed DNA methylation pathway exhibit loss of CHH methylation at many loci as well as CG and CHG methylation at a small number of loci. Plants containing loss-of-function alleles for chromomethylase (CMT) genes exhibit strong genome-wide reductions in CHG methylation and some locus-specific loss of CHH methylation. In an attempt to identify stocks with stronger reductions in DNA methylation levels than provided by single gene mutations, we performed crosses to create double mutants for the maize CMT3 orthologs, Zmet2 and Zmet5, and for the maize DDM1 orthologs, Chr101 and Chr106. While loss-of-function alleles are viable as single gene mutants, the double mutants were not recovered, suggesting that severe perturbations of the maize methylome may have stronger deleterious phenotypic effects than in Arabidopsis thaliana.

Targeted enrichment beyond the consensus coding DNA sequence exome reveals exons with higher variant densities.

BACKGROUND: Enrichment of loci by DNA hybridization-capture, followed by high-throughput sequencing, is an important tool in modern genetics. Currently, the most common targets for enrichment are the protein coding exons represented by the consensus coding DNA sequence (CCDS). The CCDS, however, excludes many actual or computationally predicted coding exons present in other databases, such as RefSeq and Vega, and non-coding functional elements such as untranslated and regulatory regions. The number of variants per base pair (variant density) and our ability to interrogate regions outside of the CCDS regions is consequently less well understood. RESULTS: We examine capture sequence data from outside of the CCDS regions and find that extremes of GC content that are present in different subregions of the genome can reduce the local capture sequence coverage to less than 50% relative to the CCDS. This effect is due to biases inherent in both the Illumina and SOLiD sequencing platforms that are exacerbated by the capture process. Interestingly, for two subregion types, microRNA and predicted exons, the capture process yields higher than expected coverage when compared to whole genome sequencing. Lastly, we examine the variation present in non-CCDS regions and find that predicted exons, as well as exonic regions specific to RefSeq and Vega, show much higher variant densities than the CCDS. CONCLUSIONS: We show that regions outside of the CCDS perform less efficiently in capture sequence experiments. Further, we show that the variant density in computationally predicted exons is more than 2.5-times higher than that observed in the CCDS.

Whole exome capture in solution with 3 Gbp of data.

We have developed a solution-based method for targeted DNA capture-sequencing that is directed to the complete human exome. Using this approach allows the discovery of greater than 95% of all expected heterozygous singe base variants, requires as little as 3 Gbp of raw sequence data and constitutes an effective tool for identifying rare coding alleles in large scale genomic studies.

A triplet repeat expansion genetic mouse model of infantile spasms syndrome, Arx(GCG)10+7, with interneuronopathy, spasms in infancy, persistent seizures, and adult cognitive and behavioral impairment.

Infantile spasms syndrome (ISS) is a catastrophic pediatric epilepsy with motor spasms, persistent seizures, mental retardation, and in some cases, autism. One of its monogenic causes is an insertion mutation [c.304ins (GCG)(7)] on the X chromosome, expanding the first polyalanine tract of the interneuron-specific transcription factor Aristaless-related homeobox (ARX) from 16 to 23 alanine codons. Null mutation of the Arx gene impairs GABA and cholinergic interneuronal migration but results in a neonatal lethal phenotype. We developed the first viable genetic mouse model of ISS that spontaneously recapitulates salient phenotypic features of the human triplet repeat expansion mutation. Arx((GCG)10+7) ("Arx plus 7") pups display abnormal spasm-like myoclonus and other key EEG features, including multifocal spikes, electrodecremental episodes, and spontaneous seizures persisting into maturity. The neurobehavioral profile of Arx mutants was remarkable for lowered anxiety, impaired associative learning, and abnormal social interaction. Laminar decreases of Arx+ cortical interneurons and a selective reduction of calbindin-, but not parvalbumin- or calretinin-expressing interneurons in neocortical layers and hippocampus indicate that specific classes of synaptic inhibition are missing from the adult forebrain, providing a basis for the seizures and cognitive disorder. A significant reduction of calbindin-, NPY (neuropeptide Y)-expressing, and cholinergic interneurons in the mutant striatum suggest that dysinhibition within this network may contribute to the dyskinetic motor spasms. This mouse model narrows the range of critical pathogenic elements within brain inhibitory networks essential to recreate this complex neurodevelopmental syndrome.

Stargazin and other transmembrane AMPA receptor regulating proteins interact with synaptic scaffolding protein MAGI-2 in brain.

The spatial coordination of neurotransmitter receptors with other postsynaptic signaling and structural molecules is regulated by a diverse array of cell-specific scaffolding proteins. The synaptic trafficking of AMPA receptors by the stargazin protein in some neurons, for example, depends on specific interactions between the C terminus of stargazin and the PDZ [postsynaptic density-95 (PSD-95)/Discs large/zona occludens-1] domains of membrane-associated guanylate kinase scaffolding proteins PSD-93 or PSD-95. Stargazin [Cacng2 (Ca2+ channel gamma2 subunit)] is one of four closely related proteins recently categorized as transmembrane AMPA receptor regulating proteins (TARPs) that appear to share similar functions but exhibit distinct expression patterns in the CNS. We used yeast two-hybrid screening to identify MAGI-2 (membrane associated guanylate kinase, WW and PDZ domain containing 2) as a novel candidate interactor with the cytoplasmic C termini of the TARPs. MAGI-2 [also known as S-SCAM (synaptic scaffolding molecule)] is a multi-PDZ domain scaffolding protein that interacts with several different ligands in brain, including PTEN (phosphatase and tensin homolog), dasm1 (dendrite arborization and synapse maturation 1), dendrin, axin, beta- and delta-catenin, neuroligin, hyperpolarization-activated cation channels, beta1-adrenergic receptors, and NMDA receptors. We confirmed that MAGI-2 coimmunoprecipitated with stargazin in vivo from mouse cerebral cortex and used in vitro assays to localize the interaction to the C-terminal -TTPV amino acid motif of stargazin and the PDZ1, PDZ3, and PDZ5 domains of MAGI-2. Expression of stargazin recruited MAGI-2 to cell membranes and cell-cell contact sites in transfected HEK-293T cells dependent on the presence of the stargazin -TTPV motif. These experiments identify MAGI-2 as a strong candidate for linking TARP/AMPA receptor complexes to a wide range of other postsynaptic molecules and pathways and advance our knowledge of protein interactions at mammalian CNS synapses.

Neonatal epilepsy syndromes and GEFS+: mechanistic considerations.

Genetic analyses of familial epilepsies over the past decade have identified mutations in several different ion channel genes that result in neonatal or early-onset seizure disorders, including benign familial neonatal convulsions (BFNC), generalized epilepsy with febrile seizures plus (GEFS+), and severe myoclonic epilepsy of infancy (SMEI). These genes encode voltage-gated Na+ channel subunits (SCN1A, SCN2A, SCN1B), voltage-gated K+ channel subunits (KCNQ2, KCNQ3), and a ligand-gated neurotransmitter receptor subunit (GABRG2). While the opportunity to genotype patients for mutations in these genes can have an immediate and significant impact on our ability to diagnose and provide genetic counseling to patients, the ultimate goal is to use this molecular knowledge to develop effective treatments and cures for each disorder. This will necessitate elucidation of the molecular, cellular, and network mechanisms that translate ion channel defects into specific epilepsy phenotypes. The functional analysis of epileptogenic channel mutations in vitro and in vivo has already provided a vast amount of raw biophysical data, but attempts to interpret these data to explain clinical phenotypes so far appear to raise as many questions as they answer. Nevertheless, patterns are beginning to emerge from these early studies that will help define the full scope of the challenges ahead while simultaneously providing the foundation of future efforts to overcome them. Here, I discuss some of the potential mechanisms that have been uncovered recently linking mutant ion channel genes to neonatal epilepsy syndromes and GEFS+.

The alpha-amino-3-hydroxyl-5-methyl-4-isoxazolepropionate receptor trafficking regulator "stargazin" is related to the claudin family of proteins by Its ability to mediate cell-cell adhesion.

Mutations in the Cacng2 gene encoding the neuronal transmembrane protein stargazin result in recessively inherited epilepsy and ataxia in "stargazer" mice. Functional studies suggest a dual role for stargazin, both as a modulatory gamma subunit for voltage-dependent calcium channels and as a regulator of post-synaptic membrane targeting for alpha-amino-3-hydroxyl-5-methyl-4-isoxazolepropionate (AMPA)-type glutamate receptors. Co-immunoprecipitation experiments demonstrate that stargazin can bind proteins of either complex in vivo, but it remains unclear whether it can associate with both complexes simultaneously. Cacng2 is one of eight closely related genes (Cacng1-8) encoding proteins with four transmembrane segments, cytoplasmic termini, and molecular masses between 25 and 44 kDa. This group of Cacng genes constitutes only one branch of a larger monophyletic assembly dominated by over 20 genes encoding proteins known as claudins. Claudins regulate cell adhesion and paracellular permeability as fundamental components of non-neuronal tight junctions. Because stargazin is structurally similar to claudins, we hypothesized that it might also have retained claudin-like functions inherited from a common ancestor. Here, we report that expression of stargazin in mouse L-fibroblasts results in cell aggregation comparable with that produced by claudins, and present evidence that the interaction is heterotypic and calcium dependent. The data suggest that the cell adhesion function of stargazin preceded its current role in neurons as a regulator of either voltage-dependent calcium channels or AMPA receptors. We speculate these complexes may have co-opted the established presence of stargazin at sites of close cell-cell contact to facilitate their own evolving intercellular signaling functions.

Allelic mutations of the sodium channel SCN8A reveal multiple cellular and physiological functions.

Allelic mutations of Scn8a in the mouse have revealed the range of neurological disorders that can result from alternations of one neuronal sodium channel. Null mutations produce the most severe phenotype, with motor neuron failure leading to paralysis and juvenile lethality. Two less severe mutations cause ataxia, tremor, muscle weakness, and dystonia. The electrophysiological effects have been studied at the cellular level by recording from neurons from the mutant mice. The data demonstrate that Scn8a is required for the complex spiking of cerebellar Purkinje cells and for persistent sodium current in several classes of neurons, including some with pacemaker roles. The mouse mutations of Scn8a have also provided insight into the mode of inheritance of channelopathies, and led to the identification of a modifier gene that affects transcript splicing. These mutations demonstrate the value of mouse models to elucidate the pathophysiology of human disease.

Mutations in high-voltage-activated calcium channel genes stimulate low-voltage-activated currents in mouse thalamic relay neurons.

Ca2+ currents, especially those activated at low voltages (LVA), influence burst generation in thalamocortical circuitry and enhance the abnormal rhythmicity associated with absence epilepsy. Mutations in several genes for high-voltage-activated (HVA) Ca2+ channel subunits are linked to spike-wave seizure phenotypes in mice; however, none of these mutations are predicted to increase intrinsic membrane excitability or directly enhance LVA currents. We examined biophysical properties of both LVA and HVA Ca2+ currents in thalamic cells of tottering (tg; Cav2.1/alpha1A subunit), lethargic (lh; beta4 subunit), and stargazer (stg; gamma2 subunit) brain slices. We observed 46, 51, and 45% increases in peak current densities of LVA Ca2+ currents evoked at -50 mV from -110 mV in tg, lh, and stg mice, respectively, compared with wild type. The half-maximal voltages for steady-state inactivation of LVA currents were shifted in a depolarized direction by 7.5-13.5 mV in all three mutants, although no alterations in the time-constant for recovery from inactivation of LVA currents were found. HVA peak current densities in tg and stg were increased by 22 and 45%, respectively, and a 5 mV depolarizing shift of the activation curve was observed in lh. Despite elevated LVA amplitudes, no alterations in mRNA expression of the genes mediating T-type subunits, Cav3.1/alpha1G, Cav3.2/alpha1H, or Cav3.3/alpha1I, were detected in the three mutants. Our data demonstrate that mutation of Cav2.1 or regulatory subunit genes increases intrinsic membrane excitability in thalamic neurons by potentiating LVA Ca2+ currents. These alterations increase the probability for abnormal thalamocortical synchronization and absence epilepsy in tg, lh, and stg mice.