Rare coding variants in ten genes confer substantial risk for schizophrenia.

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

Robertsonian Fusion and Centromere Repositioning Contributed to the Formation of Satellite-free Centromeres During the Evolution of Zebras.

Centromeres are epigenetically specified by the histone H3 variant CENP-A and typically associated with highly repetitive satellite DNA. We previously discovered natural satellite-free neocentromeres in Equus caballus and Equus asinus. Here, through ChIP-seq with an anti-CENP-A antibody, we found an extraordinarily high number of centromeres lacking satellite DNA in the zebras Equus burchelli (15 of 22) and Equus grevyi (13 of 23), demonstrating that the absence of satellite DNA at the majority of centromeres is compatible with genome stability and species survival and challenging the role of satellite DNA in centromere function. Nine satellite-free centromeres are shared between the two species in agreement with their recent separation. We assembled all centromeric regions and improved the reference genome of E. burchelli. Sequence analysis of the CENP-A binding domains revealed that they are LINE-1 and AT-rich with four of them showing DNA amplification. In the two zebras, satellite-free centromeres emerged from centromere repositioning or following Robertsonian fusion. In five chromosomes, the centromeric function arose near the fusion points, which are located within regions marked by traces of ancestral pericentromeric sequences. Therefore, besides centromere repositioning, Robertsonian fusions are an important source of satellite-free centromeres during evolution. Finally, in one case, a satellite-free centromere was seeded on an inversion breakpoint. At 11 chromosomes, whose primary constrictions seemed to be associated with satellite repeats by cytogenetic analysis, satellite-free neocentromeres were instead located near the ancestral inactivated satellite-based centromeres; therefore, the centromeric function has shifted away from a satellite repeat containing locus to a satellite-free new position.

Sequencing and de novo assembly of 150 genomes from Denmark as a population reference.

Hundreds of thousands of human genomes are now being sequenced to characterize genetic variation and use this information to augment association mapping studies of complex disorders and other phenotypic traits. Genetic variation is identified mainly by mapping short reads to the reference genome or by performing local assembly. However, these approaches are biased against discovery of structural variants and variation in the more complex parts of the genome. Hence, large-scale de novo assembly is needed. Here we show that it is possible to construct excellent de novo assemblies from high-coverage sequencing with mate-pair libraries extending up to 20 kilobases. We report de novo assemblies of 150 individuals (50 trios) from the GenomeDenmark project. The quality of these assemblies is similar to those obtained using the more expensive long-read technology. We use the assemblies to identify a rich set of structural variants including many novel insertions and demonstrate how this variant catalogue enables further deciphering of known association mapping signals. We leverage the assemblies to provide 100 completely resolved major histocompatibility complex haplotypes and to resolve major parts of the Y chromosome. Our study provides a regional reference genome that we expect will improve the power of future association mapping studies and hence pave the way for precision medicine initiatives, which now are being launched in many countries including Denmark.

hgtseq: A Standard Pipeline to Study Horizontal Gene Transfer.

Horizontal gene transfer (HGT) is well described in prokaryotes: it plays a crucial role in evolution, and has functional consequences in insects and plants. However, less is known about HGT in humans. Studies have reported bacterial integrations in cancer patients, and microbial sequences have been detected in data from well-known human sequencing projects. Few of the existing tools for investigating HGT are highly automated. Thanks to the adoption of Nextflow for life sciences workflows, and to the standards and best practices curated by communities such as nf-core, fully automated, portable, and scalable pipelines can now be developed. Here we present nf-core/hgtseq to facilitate the analysis of HGT from sequencing data in different organisms. We showcase its performance by analysing six exome datasets from five mammals. Hgtseq can be run seamlessly in any computing environment and accepts data generated by existing exome and whole-genome sequencing projects; this will enable researchers to expand their analyses into this area. Fundamental questions are still open about the mechanisms and the extent or role of horizontal gene transfer: by releasing hgtseq we provide a standardised tool which will enable a systematic investigation of this phenomenon, thus paving the way for a better understanding of HGT.

Quantifying the Impact of Rare and Ultra-rare Coding Variation across the Phenotypic Spectrum.

There is a limited understanding about the impact of rare protein-truncating variants across multiple phenotypes. We explore the impact of this class of variants on 13 quantitative traits and 10 diseases using whole-exome sequencing data from 100,296 individuals. Protein-truncating variants in genes intolerant to this class of mutations increased risk of autism, schizophrenia, bipolar disorder, intellectual disability, and ADHD. In individuals without these disorders, there was an association with shorter height, lower education, increased hospitalization, and reduced age at enrollment. Gene sets implicated from GWASs did not show a significant protein-truncating variants burden beyond what was captured by established Mendelian genes. In conclusion, we provide a thorough investigation of the impact of rare deleterious coding variants on complex traits, suggesting widespread pleiotropic risk.

Mutations in SNX14 cause a distinctive autosomal-recessive cerebellar ataxia and intellectual disability syndrome.

Intellectual disability and cerebellar atrophy occur together in a large number of genetic conditions and are frequently associated with microcephaly and/or epilepsy. Here we report the identification of causal mutations in Sorting Nexin 14 (SNX14) found in seven affected individuals from three unrelated consanguineous families who presented with recessively inherited moderate-severe intellectual disability, cerebellar ataxia, early-onset cerebellar atrophy, sensorineural hearing loss, and the distinctive association of progressively coarsening facial features, relative macrocephaly, and the absence of seizures. We used homozygosity mapping and whole-exome sequencing to identify a homozygous nonsense mutation and an in-frame multiexon deletion in two families. A homozygous splice site mutation was identified by Sanger sequencing of SNX14 in a third family, selected purely by phenotypic similarity. This discovery confirms that these characteristic features represent a distinct and recognizable syndrome. SNX14 encodes a cellular protein containing Phox (PX) and regulator of G protein signaling (RGS) domains. Weighted gene coexpression network analysis predicts that SNX14 is highly coexpressed with genes involved in cellular protein metabolism and vesicle-mediated transport. All three mutations either directly affected the PX domain or diminished SNX14 levels, implicating a loss of normal cellular function. This manifested as increased cytoplasmic vacuolation as observed in cultured fibroblasts. Our findings indicate an essential role for SNX14 in neural development and function, particularly in development and maturation of the cerebellum.

Mutation of SALL2 causes recessive ocular coloboma in humans and mice.

Ocular coloboma is a congenital defect resulting from failure of normal closure of the optic fissure during embryonic eye development. This birth defect causes childhood blindness worldwide, yet the genetic etiology is poorly understood. Here, we identified a novel homozygous mutation in the SALL2 gene in members of a consanguineous family affected with non-syndromic ocular coloboma variably affecting the iris and retina. This mutation, c.85G>T, introduces a premature termination codon (p.Glu29*) predicted to truncate the SALL2 protein so that it lacks three clusters of zinc-finger motifs that are essential for DNA-binding activity. This discovery identifies SALL2 as the third member of the Drosophila homeotic Spalt-like family of developmental transcription factor genes implicated in human disease. SALL2 is expressed in the developing human retina at the time of, and subsequent to, optic fissure closure. Analysis of Sall2-deficient mouse embryos revealed delayed apposition of the optic fissure margins and the persistence of an anterior retinal coloboma phenotype after birth. Sall2-deficient embryos displayed correct posterior closure toward the optic nerve head, and upon contact of the fissure margins, dissolution of the basal lamina occurred and PAX2, known to be critical for this process, was expressed normally. Anterior closure was disrupted with the fissure margins failing to meet, or in some cases misaligning leading to a retinal lesion. These observations demonstrate, for the first time, a role for SALL2 in eye morphogenesis and that loss of function of the gene causes ocular coloboma in humans and mice.

The kinetochore protein, CENPF, is mutated in human ciliopathy and microcephaly phenotypes.

BACKGROUND: Mutations in microtubule-regulating genes are associated with disorders of neuronal migration and microcephaly. Regulation of centriole length has been shown to underlie the pathogenesis of certain ciliopathy phenotypes. Using a next-generation sequencing approach, we identified mutations in a novel centriolar disease gene in a kindred with an embryonic lethal ciliopathy phenotype and in a patient with primary microcephaly. METHODS AND RESULTS: Whole exome sequencing data from a non-consanguineous Caucasian kindred exhibiting mid-gestation lethality and ciliopathic malformations revealed two novel non-synonymous variants in CENPF, a microtubule-regulating gene. All four affected fetuses showed segregation for two mutated alleles [IVS5-2A>C, predicted to abolish the consensus splice-acceptor site from exon 6; c.1744G>T, p.E582X]. In a second unrelated patient exhibiting microcephaly, we identified two CENPF mutations [c.1744G>T, p.E582X; c.8692 C>T, p.R2898X] by whole exome sequencing. We found that CENP-F colocalised with Ninein at the subdistal appendages of the mother centriole in mouse inner medullary collecting duct cells. Intraflagellar transport protein-88 (IFT-88) colocalised with CENP-F along the ciliary axonemes of renal epithelial cells in age-matched control human fetuses but did not in truncated cilia of mutant CENPF kidneys. Pairwise co-immunoprecipitation assays of mitotic and serum-starved HEKT293 cells confirmed that IFT88 precipitates with endogenous CENP-F. CONCLUSIONS: Our data identify CENPF as a new centriolar disease gene implicated in severe human ciliopathy and microcephaly related phenotypes. CENP-F has a novel putative function in ciliogenesis and cortical neurogenesis.

Whole-exome sequencing implicates DGKH as a risk gene for panic disorder in the Faroese population.

The demographic history of the isolated population of the Faroe Islands may have induced enrichment of variants rarely seen in outbred European populations, including enrichment of risk variants for panic disorder (PD). PD is a common mental disorder, characterized by recurring and unprovoked panic attacks, and genetic factors have been estimated to explain around 40% of the risk. In this study the potential enrichment of PD risk variants was explored based on whole-exome sequencing of 54 patients with PD and 211 control individuals from the Faroese population. No genome-wide significant associations were found, however several single variants and genes showed strong association with PD, where DGKH was found to be the strongest PD associated gene. Interestingly DGKH has previously demonstrated genome-wide significant association with bipolar disorder as well as evidence of association to other mental disorders. Additionally, we found an enrichment of PD risk variants in the Faroese population; variants with otherwise low frequency in more outbreed European populations. © 2016 Wiley Periodicals, Inc.

Experimental validation of methods for differential gene expression analysis and sample pooling in RNA-seq.

BACKGROUND: Massively parallel cDNA sequencing (RNA-seq) experiments are gradually superseding microarrays in quantitative gene expression profiling. However, many biologists are uncertain about the choice of differentially expressed gene (DEG) analysis methods and the validity of cost-saving sample pooling strategies for their RNA-seq experiments. Hence, we performed experimental validation of DEGs identified by Cuffdiff2, edgeR, DESeq2 and Two-stage Poisson Model (TSPM) in a RNA-seq experiment involving mice amygdalae micro-punches, using high-throughput qPCR on independent biological replicate samples. Moreover, we sequenced RNA-pools and compared their results with sequencing corresponding individual RNA samples. RESULTS: False-positivity rate of Cuffdiff2 and false-negativity rates of DESeq2 and TSPM were high. Among the four investigated DEG analysis methods, sensitivity and specificity of edgeR was relatively high. We documented the pooling bias and that the DEGs identified in pooled samples suffered low positive predictive values. CONCLUSIONS: Our results highlighted the need for combined use of more sensitive DEG analysis methods and high-throughput validation of identified DEGs in future RNA-seq experiments. They indicated limited utility of sample pooling strategies for RNA-seq in similar setups and supported increasing the number of biological replicate samples.

ARNT2 mutation causes hypopituitarism, post-natal microcephaly, visual and renal anomalies.

We describe a previously unreported syndrome characterized by secondary (post-natal) microcephaly with fronto-temporal lobe hypoplasia, multiple pituitary hormone deficiency, seizures, severe visual impairment and abnormalities of the kidneys and urinary tract in a highly consanguineous family with six affected children. Homozygosity mapping and exome sequencing revealed a novel homozygous frameshift mutation in the basic helix-loop-helix transcription factor gene ARNT2 (c.1373_1374dupTC) in affected individuals. This mutation results in absence of detectable levels of ARNT2 transcript and protein from patient fibroblasts compared with controls, consistent with nonsense-mediated decay of the mutant transcript and loss of ARNT2 function. We also show expression of ARNT2 within the central nervous system, including the hypothalamus, as well as the renal tract during human embryonic development. The progressive neurological abnormalities, congenital hypopituitarism and post-retinal visual pathway dysfunction in affected individuals demonstrates for the first time the essential role of ARNT2 in the development of the hypothalamo-pituitary axis, post-natal brain growth, and visual and renal function in humans.

Scalable and efficient DNA sequencing analysis on different compute infrastructures aiding variant discovery.

DNA variation analysis has become indispensable in many aspects of modern biomedicine, most prominently in the comparison of normal and tumor samples. Thousands of samples are collected in local sequencing efforts and public databases requiring highly scalable, portable, and automated workflows for streamlined processing. Here, we present nf-core/sarek 3, a well-established, comprehensive variant calling and annotation pipeline for germline and somatic samples. It is suitable for any genome with a known reference. We present a full rewrite of the original pipeline showing a significant reduction of storage requirements by using the CRAM format and runtime by increasing intra-sample parallelization. Both are leading to a 70% cost reduction in commercial clouds enabling users to do large-scale and cross-platform data analysis while keeping costs and CO2 emissions low. The code is available at https://nf-co.re/sarek.

The promise of explainable deep learning for omics data analysis: Adding new discovery tools to AI.

Deep learning has already revolutionised the way a wide range of data is processed in many areas of daily life. The ability to learn abstractions and relationships from heterogeneous data has provided impressively accurate prediction and classification tools to handle increasingly big datasets. This has a significant impact on the growing wealth of omics datasets, with the unprecedented opportunity for a better understanding of the complexity of living organisms. While this revolution is transforming the way these data are analyzed, explainable deep learning is emerging as an additional tool with the potential to change the way biological data is interpreted. Explainability addresses critical issues such as transparency, so important when computational tools are introduced especially in clinical environments. Moreover, it empowers artificial intelligence with the capability to provide new insights into the input data, thus adding an element of discovery to these already powerful resources. In this review, we provide an overview of the transformative effects explainable deep learning is having on multiple sectors, ranging from genome engineering and genomics, from radiomics to drug design and clinical trials. We offer a perspective to life scientists, to better understand the potential of these tools, and a motivation to implement them in their research, by suggesting learning resources they can use to move their first steps in this field.

Omics in a Digital World: The Role of Bioinformatics in Providing New Insights Into Human Aging.

BACKGROUND: Aging is a complex phenotype influenced by a combination of genetic and environmental factors. Although many studies addressed its cellular and physiological age-related changes, the molecular causes of aging remain undetermined. Considering the biological complexity and heterogeneity of the aging process, it is now clear that full understanding of mechanisms underlying aging can only be achieved through the integration of different data types and sources, and with new computational methods capable to achieve such integration. RECENT ADVANCES: In this review, we show that an omics vision of the age-dependent changes occurring as the individual ages can provide researchers with new opportunities to understand the mechanisms of aging. Combining results from single-cell analysis with systems biology tools would allow building interaction networks and investigate how these networks are perturbed during aging and disease. The development of high-throughput technologies such as next-generation sequencing, proteomics, metabolomics, able to investigate different biological markers and to monitor them simultaneously during the aging process with high accuracy and specificity, represents a unique opportunity offered to biogerontologists today. CRITICAL ISSUES: Although the capacity to produce big data drastically increased over the years, integration, interpretation and sharing of high-throughput data remain major challenges. In this paper we present a survey of the emerging omics approaches in aging research and provide a large collection of datasets and databases as a useful resource for the scientific community to identify causes of aging. We discuss their peculiarities, emphasizing the need for the development of methods focused on the integration of different data types. FUTURE DIRECTIONS: We critically review the contribution of bioinformatics into the omics of aging research, and we propose a few recommendations to boost collaborations and produce new insights. We believe that significant advancements can be achieved by following major developments in bioinformatics, investing in diversity, data sharing and community-driven portable bioinformatics methods. We also argue in favor of more engagement and participation, and we highlight the benefits of new collaborations along these lines. This review aims at being a useful resource for many researchers in the field, and a call for new partnerships in aging research.

Association between the interleukin-1beta polymorphisms and Alzheimer's disease: a systematic review and meta-analysis.

The pro-inflammatory cytokine interleukin(IL)-1beta is a main component in inflammatory pathways and is overexpressed in the brain of Alzheimer's disease (AD) patients. Several studies report associations between IL-1beta polymorphisms and AD, but findings from different studies are controversial. Our aim was to verify the correlation between the single nucleotide polymorphisms (SNPs) of the IL-1beta, at sites -511 and +3953, and AD by meta-analysis. Computerized bibliographic searches of PUBMED and AlzGene database (http://www.alzgene.org) were supplemented with manual searches of reference lists. There is evidence for association between IL-1beta +3953 SNP and AD, with an OR=1.60 (95% C.I.: 1.16-2.22; Z=2.83 p=0.005) for TT genotype. No significant difference in genotype distribution of the IL-1beta -511 SNP in AD was obtained, but high between-study heterogeneity was found. To reduce heterogeneity, subgroup analyses were performed using, as stratifying variables, characteristics of the population under study (age, gender, type of AD diagnosis, Mini Mental State Examination of the controls) and characteristics related to the study design (statistical power of individual studies). The frequency of the IL-1beta -511 TT genotype resulted significantly higher than other genotypes only when the Caucasian studies with the highest statistical power were included in the subgroup analysis (OR=1.32; 95% C.I.: 1.03-1.69; p=0.03), with no evidence of between-study heterogeneity. Our data support an association between the TT genotype of IL-1beta +3953 SNP and AD, and suggest a possible association of the -511 TT genotype. Unreplicability of the results seems to be due mainly to the lack of statistical power of the individual studies.

Autism spectrum disorder and attention deficit hyperactivity disorder have a similar burden of rare protein-truncating variants.

The exome sequences of approximately 8,000 children with autism spectrum disorder (ASD) and/or attention deficit hyperactivity disorder (ADHD) and 5,000 controls were analyzed, finding that individuals with ASD and individuals with ADHD had a similar burden of rare protein-truncating variants in evolutionarily constrained genes, both significantly higher than controls. This motivated a combined analysis across ASD and ADHD, identifying microtubule-associated protein 1A (MAP1A) as a new exome-wide significant gene conferring risk for childhood psychiatric disorders.

A genetic-demographic approach reveals male-specific association between survival and tumor necrosis factor (A/G)-308 polymorphism.

The (A/G)-308 polymorphism of the tumor necrosis factor alpha gene (TNF) is associated with age-related diseases, but its influence on longevity is controversial. We genotyped for this polymorphism 747 Italian volunteers (401 women and 346 men, age 19-110 years). By applying a genetic-demographic (GD) approach we found that, in men, the survival function of allele A carriers is lower than that of noncarriers at all the ages (p =.044). After defining (by exploiting again demographic information) three age classes, we found that the frequency of men carrying the A allele decreases with age (p =.019), thus confirming the GD analysis results. The same analyses gave negative results in women. Therefore, allele A has a detrimental effect on life expectancy, and this effect is specific to men. A haplotype analysis carried out in men by screening the TNFa, TNFc, and TNFe microsatellite polymorphisms (spanning about 20 kb) confirmed the association of the TNF region with life expectancy.

Role of mitochondrial DNA in longevity, aging and age-related diseases in humans: a reappraisal.

The genetic variability of H. sapiens mitochondrial DNA (mtDNA) can be either germ-line inherited or somatically acquired, and its effect on aging and longevity as well as on the pathogenesis of complex age-related diseases is a hot topic. Here we illustrate the complexity of such studies, related to the large genetic variability of mtDNA in different populations and the fact that the rate of the aging process is different in different cells, tissues and organs. As far as concern Alzheimer's disease, the accumulation of somatic mutations in several tissues have been investigated, as well as the inherited mtDNA variability. However, the issue is still controversial and further studies are needed to clarify the role of mtDNA variants in Alzheimer's disease. This review is aimed to summarize the most recent advances in this field. By high throughput mtDNA sequencing and the study of large cohorts of ethnically homogeneous subjects/patients, it is now possible to perform high dimensionality studies in order to clarify the genetic associations among several inherited mtDNA variants and longevity or age-associated diseases in humans.

The use of whole-exome sequencing to disentangle complex phenotypes.

The success of whole-exome sequencing to identify mutations causing single-gene disorders has been well documented. In contrast whole-exome sequencing has so far had limited success in the identification of variants causing more complex phenotypes that seem unlikely to be due to the disruption of a single gene. We describe a family where two male offspring of healthy first cousin parents present a complex phenotype consisting of peripheral neuropathy and bronchiectasis that has not been described previously in the literature. Due to the fact that both children had the same problems in the context of parental consanguinity we hypothesised illness resulted from either X-linked or autosomal recessive inheritance. Through the use of whole-exome sequencing we were able to simplify this complex phenotype and identified a causative mutation (p.R1070*) in the gene periaxin (PRX), a gene previously shown to cause peripheral neuropathy (Dejerine-Sottas syndrome) when this mutation is present. For the bronchiectasis phenotype we were unable to identify a causal single mutation or compound heterozygote, reflecting the heterogeneous nature of this phenotype. In conclusion, in this study we show that whole-exome sequencing has the power to disentangle complex phenotypes through the identification of causative genetic mutations for distinct clinical disorders that were previously masked.

STAG3 truncating variant as the cause of primary ovarian insufficiency.

Primary ovarian insufficiency (POI) is a distressing cause of infertility in young women. POI is heterogeneous with only a few causative genes having been discovered so far. Our objective was to determine the genetic cause of POI in a consanguineous Lebanese family with two affected sisters presenting with primary amenorrhoea and an absence of any pubertal development. Multipoint parametric linkage analysis was performed. Whole-exome sequencing was done on the proband. Linkage analysis identified a locus on chromosome 7 where exome sequencing successfully identified a homozygous two base pair duplication (c.1947_48dupCT), leading to a truncated protein p.(Y650Sfs*22) in the STAG3 gene, confirming it as the cause of POI in this family. Exome sequencing combined with linkage analyses offers a powerful tool to efficiently find novel genetic causes of rare, heterogeneous disorders, even in small single families. This is only the second report of a STAG3 variant; the first STAG3 variant was recently described in a phenotypically similar family with extreme POI. Identification of an additional family highlights the importance of STAG3 in POI pathogenesis and suggests it should be evaluated in families affected with POI.