Whole genome sequencing of 91 multiplex schizophrenia families reveals increased burden of rare, exonic copy number variation in schizophrenia probands and genetic heterogeneity

Investigating copy number variants in schizophrenia pedigrees using a new consensus pipeline called PECAN

Copy number variants (CNVs) have been implicated in many human diseases, including psychiatric disorders. Whole genome sequencing offers advantages in CNV calling compared to previous array-based methods. Here we present a robust and transparent CNV calling pipeline, PECAN (PEdigree Copy number vAriaNt calling), for short-read, whole genome sequencing data, comprised of a novel combination of four calling methods and structural variant genotyping. This method is scalable and can incorporate pedigree information to retain lower-confidence CNVs that would otherwise be discarded. We have robustly benchmarked PECAN using gold-standard CNV calls for two well-established evaluation samples, NA12878 and HG002, showing that PECAN performs with high precision and recall on both datasets, outperforming another pedigree-based CNV calling pipeline. As part of this work, we provide a list of high-confidence gold standard CNVs for the NA12878 reference sample, curated from multiple studies. We applied PECAN to a collection of pedigrees multiply affected with schizophrenia and identified a rare deletion that perfectly co-segregates with schizophrenia in one of the pedigrees. The CNV overlaps the gene PITRM1, which has been implicated in a complex phenotype including ataxia, developmental delay, and schizophrenia-like episodes in affected adults.

Third-generation genome sequencing implicates medium-sized structural variants in chronic schizophrenia

Background

Schizophrenia (SCZ) is a heterogeneous psychiatric disorder, with significant contribution from genetic factors particularly for chronic cases with negative symptoms and cognitive deficits. To date, Genome Wide Association Studies (GWAS) and exome sequencing have associated SCZ with a number of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs), but there is still missing heritability. Medium-sized structural variants (SVs) are difficult to detect using SNP arrays or second generation sequencing, and may account for part of the missing heritability of SCZ.

Aims and objectives

To identify SVs associated with severe chronic SCZ across the whole genome.

Study design

10 multiplex families with probands suffering from chronic SCZ with negative symptoms and cognitive deficits were recruited, with all their affected members demonstrating uni-lineal inheritance. Control subjects comprised one affected member from the affected lineage, and unaffected members from each paternal and maternal lineage.

Methods

Third generation sequencing was applied to peripheral blood samples from 10 probands and 5 unaffected controls. Bioinformatic tools were used to identify SVs from the long sequencing reads, with confirmation of findings in probands by short-read Illumina sequencing, Sanger sequencing and visual manual validation with Integrated Genome Browser.

Results

In the 10 probands, we identified and validated 88 SVs (mostly in introns and medium-sized), within 79 genes, which were absent in the 5 unaffected control subjects. These 79 genes were enriched in 20 biological pathways which were related to brain development, neuronal migration, neurogenesis, neuronal/synaptic function, learning/memory, and hearing. These identified SVs also showed evidence for enrichment of genes that are highly expressed in the adolescent striatum.

Conclusion

A substantial part of the missing heritability in SCZ may be explained by medium-sized SVs detectable only by third generation sequencing. We have identified a number of such SVs potentially conferring risk for SCZ, which implicate multiple brain-related genes and pathways. In addition to previously-identified pathways involved in SCZ such as neurodevelopment and neuronal/synaptic functioning, we also found novel evidence for enrichment in hearing-related pathways and genes expressed in the adolescent striatum.

Identification of copy number variations in the genome of Dairy Gir cattle

Studying structural variants that can control complex traits is relevant for dairy cattle production, especially for animals that are tolerant to breeding conditions in the tropics, such as the Dairy Gir cattle. This study identified and characterized high confidence copy number variation regions (CNVR) in the Gir breed genome. A total of 38 animals were whole-genome sequenced, and 566 individuals were genotyped with a high-density SNP panel, among which 36 animals had both sequencing and SNP genotyping data available. Two sets of high confidence CNVR were established: one based on common CNV identified in the studied population (CNVR_POP), and another with CNV identified in sires with both sequence and SNP genotyping data available (CNVR_ANI). We found 10 CNVR_POP and 45 CNVR_ANI, which covered 1.05 Mb and 4.4 Mb of the bovine genome, respectively. Merging these CNV sets for functional analysis resulted in 48 unique high confidence CNVR. The overlapping genes were previously related to embryonic mortality, environmental adaptation, evolutionary process, immune response, longevity, mammary gland, resistance to gastrointestinal parasites, and stimuli recognition, among others. Our results contribute to a better understanding of the Gir breed genome. Moreover, the CNV identified in this study can potentially affect genes related to complex traits, such as production, health, and reproduction.

Exonic deletions in IMMP2L in schizophrenia with enhanced glycation stress subtype

We previously identified a subtype of schizophrenia (SCZ) characterized by increased plasma pentosidine, a marker of glycation and oxidative stress (PEN-SCZ). However, the genetic factors associated with PEN-SCZ have not been fully clarified. We performed a genome-wide copy number variation (CNV) analysis to identify CNVs associated with PEN-SCZ to provide an insight into the novel therapeutic targets for PEN-SCZ. Plasma pentosidine was measured by high-performance liquid chromatography in 185 patients with SCZ harboring rare CNVs detected by array comparative genomic hybridization. In three patients with PEN-SCZ showing additional autistic features, we detected a novel deletion at 7q31.1 within exons 2 and 3 of IMMP2L, which encodes the inner mitochondrial membrane peptidase subunit 2. The deletion was neither observed in non-PEN-SCZ nor in public database of control subjects. IMMP2L is one of the SCZ risk loci genes identified in a previous SCZ genome-wide association study, and its trans-populational association was recently described. Interestingly, deletions in IMMP2L have been previously linked with autism spectrum disorder. Disrupted IMMP2L function has been shown to cause glycation/oxidative stress in neuronal cells in an age-dependent manner. To our knowledge, this is the first genome-wide CNV study to suggest the involvement of IMMP2L exons 2 and 3 in the etiology of PEN-SCZ. The combination of genomic information with plasma pentosidine levels may contribute to the classification of biological SCZ subtypes that show additional autistic features. Modifying IMMP2L functions may be useful for treating PEN-SCZ if the underlying biological mechanism can be clarified in further studies.

Exome sequencing of Japanese schizophrenia multiplex families supports the involvement of calcium ion channels

Background

Most sequencing studies of schizophrenia (SCZ) have focused on de novo genetic variants due to interpretability. However, investigating shared rare variants among patients in the same multiplex family is also important. Relatively large-scale analyses of SCZ multiplex families have been done in Caucasian populations, but whether detected variants are also pathogenic in the Japanese population is unclear because of ethnic differences in rare variants.

Materials and methods

We performed whole-exome sequencing (WES) of 14 Japanese SCZ multiplex families. After quality control and filtering, we identified rare variants shared among affected persons within the same family. A gene ontology (GO) analysis was performed to identify gene categories possibly affected by these candidate variants.

Results

We found 530 variants in 486 genes as potential candidate variants from the 14 SCZ multiplex families examined. The GO analysis demonstrated significant enrichment in calcium channel activity.

Conclusion

This study provides supporting evidence that calcium ion channel activity is involved in SCZ. WES of multiplex families is a potential means of identifying disease-associated rare variants for SCZ.

Novel variant in HHAT as a cause of different sex development with partial gonadal dysgenesis associated with microcephaly, eye defects, and distal phalangeal hypoplasia of both thumbs: Case report

The palmitoylation of the Hedgehog (Hh) family of morphogens, named sonic hedgehog (SHH), desert hedgehog (DHH), and Indian hedgehog (IHH), is crucial for effective short- and long-range signaling. The hedgehog acyltransferase (HHAT) attaches the palmitate molecule to the Hh; therefore, variants in HHAT cause a broad spectrum of phenotypes. A missense HHAT novel variant c.1001T>A/p.(Met334Lys) was described in a patient first referred for a 46,XY different sexual development with partial gonadal dysgenesis but with microcephaly, eye defects, and distal phalangeal hypoplasia of both thumbs. The in silico analysis of the variant predicted an affectation of the nearest splicing site. Thus, in vitro minigene studies were carried out, which demonstrated that the variant does not affect the splicing. Subsequent protein in silico studies supported the pathogenicity of the variant, and, in conclusion, this was considered the cause of the patient's phenotype.

The conserved ASTN2/BRINP1 locus at 9q33.1-33.2 is associated with major psychiatric disorders in a large pedigree from Southern Spain

We investigated the genetic causes of major mental disorders (MMDs) including schizophrenia, bipolar disorder I, major depressive disorder and attention deficit hyperactive disorder, in a large family pedigree from Alpujarras, South of Spain, a region with high prevalence of psychotic disorders. We applied a systematic genomic approach based on karyotyping (n = 4), genotyping by genome-wide SNP array (n = 34) and whole-genome sequencing (n = 12). We performed genome-wide linkage analysis, family-based association analysis and polygenic risk score estimates. Significant linkage was obtained at chromosome 9 (9q33.1–33.2, LOD score = 4.11), a suggestive region that contains five candidate genes ASTN2, BRINP1, C5, TLR4 and TRIM32, previously associated with MMDs. Comprehensive analysis associated the MMD phenotype with genes of the immune system with dual brain functions. Moreover, the psychotic phenotype was enriched for genes involved in synapsis. These results should be considered once studying the genetics of psychiatric disorders in other families, especially the ones from the same region, since founder effects may be related to the high prevalence.

Accuracy and efficiency of germline variant calling pipelines for human genome data

Advances in next-generation sequencing technology have enabled whole genome sequencing (WGS) to be widely used for identification of causal variants in a spectrum of genetic-related disorders, and provided new insight into how genetic polymorphisms affect disease phenotypes. The development of different bioinformatics pipelines has continuously improved the variant analysis of WGS data. However, there is a necessity for a systematic performance comparison of these pipelines to provide guidance on the application of WGS-based scientific and clinical genomics. In this study, we evaluated the performance of three variant calling pipelines (GATK, DRAGEN and DeepVariant) using the Genome in a Bottle Consortium, "synthetic-diploid" and simulated WGS datasets. DRAGEN and DeepVariant show better accuracy in SNP and indel calling, with no significant differences in their F1-score. DRAGEN platform offers accuracy, flexibility and a highly-efficient execution speed, and therefore superior performance in the analysis of WGS data on a large scale. The combination of DRAGEN and DeepVariant also suggests a good balance of accuracy and efficiency as an alternative solution for germline variant detection in further applications. Our results facilitate the standardization of benchmarking analysis of bioinformatics pipelines for reliable variant detection, which is critical in genetics-based medical research and clinical applications.

Accuracy and efficiency of germline variant calling pipelines for human genome data

Advances in next-generation sequencing technology have enabled whole genome sequencing (WGS) to be widely used for identification of causal variants in a spectrum of genetic-related disorders, and provided new insight into how genetic polymorphisms affect disease phenotypes. The development of different bioinformatics pipelines has continuously improved the variant analysis of WGS data. However, there is a necessity for a systematic performance comparison of these pipelines to provide guidance on the application of WGS-based scientific and clinical genomics. In this study, we evaluated the performance of three variant calling pipelines (GATK, DRAGEN and DeepVariant) using the Genome in a Bottle Consortium, "synthetic-diploid" and simulated WGS datasets. DRAGEN and DeepVariant show better accuracy in SNP and indel calling, with no significant differences in their F1-score. DRAGEN platform offers accuracy, flexibility and a highly-efficient execution speed, and therefore superior performance in the analysis of WGS data on a large scale. The combination of DRAGEN and DeepVariant also suggests a good balance of accuracy and efficiency as an alternative solution for germline variant detection in further applications. Our results facilitate the standardization of benchmarking analysis of bioinformatics pipelines for reliable variant detection, which is critical in genetics-based medical research and clinical applications.

Discovery of rare variants implicated in schizophrenia using next-generation sequencing

Schizophrenia is a highly heritable psychiatric disorder that affects 1% of the population. Genome-wide association studies have identified common variants in candidate genes associated with schizophrenia, but the genetics mechanisms of this disorder have not yet been elucidated. The discovery of rare genetic variants that contribute to schizophrenia symptoms promises to help explain the missing heritability of the disease. Next generation sequencing techniques are revolutionizing the field of psychiatric genetics. Various statistical approaches have been developed for rare variant association testing in case-control and family studies. Targeted resequencing, whole exome sequencing and whole genome sequencing combined with these computational tools are used for the discovery of rare genetic variations in schizophrenia. The findings provide useful information for characterizing the rare mutations and elucidating the genetic mechanisms by which the variants cause schizophrenia.

Comparison of three variant callers for human whole genome sequencing

Testing of patients with genetics-related disorders is in progress of shifting from single gene assays to gene panel sequencing, whole-exome sequencing (WES) and whole-genome sequencing (WGS). Since WGS is unquestionably becoming a new foundation for molecular analyses, we decided to compare three currently used tools for variant calling of human whole genome sequencing data. We tested DeepVariant, a new TensorFlow machine learning-based variant caller, and compared this tool to GATK 4.0 and SpeedSeq, using 30×, 15× and 10× WGS data of the well-known NA12878 DNA reference sample. According to our comparison, the performance on SNV calling was almost similar in 30× data, with all three variant callers reaching F-Scores (i.e. harmonic mean of recall and precision) equal to 0.98. In contrast, DeepVariant was more precise in indel calling than GATK and SpeedSeq, as demonstrated by F-Scores of 0.94, 0.90 and 0.84, respectively. We conclude that the DeepVariant tool has great potential and usefulness for analysis of WGS data in medical genetics.