Automated prediction of the clinical impact of structural copy number variations

Copy number variants (CNVs) play an important role in many biological processes, including the development of genetic diseases, making them attractive targets for genetic analyses. The interpretation of the effect of these structural variants is a challenging problem due to highly variable numbers of gene, regulatory, or other genomic elements affected by the CNV. This led to the demand for the interpretation tools that would relieve researchers, laboratory diagnosticians, genetic counselors, and clinical geneticists from the laborious process of annotation and classification of CNVs. We designed and validated a prediction method (ISV; Interpretation of Structural Variants) that is based on boosted trees which takes into account annotations of CNVs from several publicly available databases. The presented approach achieved more than 98% prediction accuracy on both copy number loss and copy number gain variants while also allowing CNVs being assigned “uncertain” significance in predictions. We believe that ISV’s prediction capability and explainability have a great potential to guide users to more precise interpretations and classifications of CNVs.

Loss-of-function of OTUD7A in the schizophrenia-associated 15q13.3 deletion impairs synapse development and function in human neurons

Identifying causative gene(s) within disease-associated large genomic regions of copy number variants (CNVs) is challenging. Here, by targeted sequencing of genes within schizophrenia (SZ)-associated CNVs in 1,779 SZ cases and 1,418 controls, we identified three rare putative loss-of-function (LoF) mutations in OTU deubiquitinase 7A (OTUD7A) within the 15q13.3 deletion in cases, but none in controls. To tie OTUD7A LoF with any SZ-relevant cellular phenotypes, we modeled the OTUD7A LoF mutation, rs757148409, in human induced pluripotent stem cell (hiPSC)-derived induced excitatory neurons (iNs) by CRISPR/Cas9 engineering. The mutant iNs showed a ~50% decrease in OTUD7A expression without undergoing nonsense-mediated mRNA decay. The mutant iNs also exhibited marked reduction of dendritic complexity, density of synaptic proteins GluA1 and PSD-95, and neuronal network activity. Congruent with the neuronal phenotypes in mutant iNs, our transcriptomic analysis showed that the set of OTUD7A LoF-downregulated genes was enriched for those relating to synapse development and function, and was associated with SZ and other neuropsychiatric disorders. These results suggest that OTUD7A LoF impairs synapse development and neuronal function in human neurons, providing mechanistic insight into the possible role of OTUD7A in driving neuropsychiatric phenotypes associated with the 15q13.3 deletion.

Schizophrenia-associated somatic copy number variants from 12,834 cases reveal contribution to risk and recurrent, isoform-specific NRXN1 disruptions

While inherited and de novo copy number variants (CNV) have been implicated in the genetic architecture of schizophrenia (SCZ), the contribution of somatic CNVs (sCNVs), present in some but not all cells of the body, remains unknown. Here we explore the role of sCNVs in SCZ by analyzing blood-derived genotype arrays from 12,834 SCZ cases and 11,648 controls. sCNVs were more common in cases (0.91%) than in controls (0.51%, p = 2.68e-4). We observed recurrent somatic deletions of exons 1-5 of the NRXN1 gene in 5 SCZ cases. Allele-specific Hi-C maps revealed ectopic, allele-specific loops forming between a potential novel cryptic promoter and non-coding cis regulatory elements upon deletions in the 5' region of NRXN1. We also observed recurrent intragenic deletions of ABCB11, a gene associated with anti-psychotic response, in 5 treatment-resistant SCZ cases. Taken together our results indicate an important role of sCNVs to SCZ risk and treatment-responsiveness.

Somatic mutation involving diverse genes leads to a spectrum of focal cortical malformations

Post-zygotically acquired genetic variants, or somatic variants, that arise during cortical development have emerged as important causes of focal epilepsies, particularly those due to malformations of cortical development. Pathogenic somatic variants have been identified in many genes within the PI3K-AKT3-mTOR-signaling pathway in individuals with hemimegalencephaly and focal cortical dysplasia (type II), and more recently in SLC35A2 in individuals with focal cortical dysplasia (type I) or non-dysplastic epileptic cortex. Given the expanding role of somatic variants across different brain malformations, we sought to delineate the landscape of somatic variants in a large cohort of patients who underwent epilepsy surgery with hemimegalencephaly or focal cortical dysplasia. We evaluated samples from 123 children with hemimegalencephaly (n=16), focal cortical dysplasia type I and related phenotypes (n=48), focal cortical dysplasia type II (n=44), or focal cortical dysplasia type III (n=15) classified using imaging and pathological findings. We performed high-depth exome sequencing in brain tissue-derived DNA from each case and identified somatic single nucleotide, indel, and large copy number variants. In 75% of individuals with hemimegalencephaly and 29% with focal cortical dysplasia type II, we identified pathogenic variants in PI3K-AKT-mTOR pathway genes. Four of 48 cases with focal cortical dysplasia type I (8%) had a likely pathogenic variant in SLC35A2. While no other gene had multiple disease-causing somatic variants across the focal cortical dysplasia type I cohort, four individuals in this group had a single pathogenic or likely pathogenic somatic variant in CASK, KRAS, NF1, and NIPBL, genes associated with neurodevelopmental disorders. No rare pathogenic or likely pathogenic somatic variants in any neurological disease genes like those identified in the focal cortical dysplasia type I cohort were found in 63 neurologically normal controls (P = 0.017), suggesting a role for these novel variants. We also identified a somatic loss-of-function variant in the known epilepsy gene, PCDH19, present in a very small number of alleles in the dysplastic tissue from a female patient with focal cortical dysplasia IIIa with hippocampal sclerosis. In contrast to focal cortical dysplasia type II, neither focal cortical dysplasia type I nor III had somatic variants in genes that converge on a unifying biological pathway, suggesting greater genetic heterogeneity compared to type II. Importantly, we demonstrate that FCD types I, II, and III, are associated with somatic gene variants across a broad range of genes, many associated with epilepsy in clinical syndromes caused by germline variants, as well as including some not previously associated with radiographically evident cortical brain malformations.

Haplotype-resolved inversion landscape reveals hotspots of mutational recurrence associated with genomic disorders

Unlike copy number variants (CNVs), inversions remain an underexplored genetic variation class. By integrating multiple genomic technologies, we discover 729 inversions in 41 human genomes. Approximately 85% of inversions <2 kbp form by twin-priming during L1-retrotransposition; 80% of the larger inversions are balanced and affect twice as many base pairs as CNVs. Balanced inversions show an excess of common variants, and 72% are flanked by segmental duplications (SDs) or mobile elements. Since this suggests recurrence due to non-allelic homologous recombination, we developed complementary approaches to identify recurrent inversion formation. We describe 40 recurrent inversions encompassing 0.6% of the genome, showing inversion rates up to 2.7*10-4 per locus and generation. Recurrent inversions exhibit a sex-chromosomal bias, and significantly co-localize to the critical regions of genomic disorders. We propose that inversion recurrence results in an elevated number of heterozygous carriers and structural SD diversity, which increases mutability in the population and predisposes to disease-causing CNVs.

CONGA: Copy number variation genotyping in ancient genomes and low-coverage sequencing data

To date, ancient genome analyses have been largely confined to the study of single nucleotide polymorphisms (SNPs). Copy number variants (CNVs) are a major contributor of disease and of evolutionary adaptation, but identifying CNVs in ancient shotgun-sequenced genomes is hampered by (a) most published genomes being <1x coverage, (ii) ancient DNA fragments being typically <80 bps. These characteristics preclude state-of-the-art CNV detection software to be effectively applied to ancient genomes. Here we present CONGA, an algorithm tailored for genotyping deletion and duplication events in genomes with low depths of coverage. Simulations show that CONGA can genotype deletions and duplications >1 Kbps with F-scores >0.77 and >0.82, respectively at >=0.5x. Further, down-sampling experiments using published ancient BAM files reveal that >1 Kbps deletions could be genotyped at F-score >0.75 at >=1x coverage. Using CONGA, we analyse deletion events at 10,018 loci in 56 ancient human genomes spanning the last 50,000 years, with coverages 0.4x-26x. We find inter-individual genetic diversity measured using deletions and SNPs to be highly correlated, suggesting that deletion frequencies broadly reflect demographic history. We also identify signatures of purifying selection on deletions, such as an excess of singletons compared to those in SNPs. CONGA paves the way for systematic studies of drift, mutation load, and adaptation in ancient and modern-day gene pools through the lens of CNVs.

Genomic and transcriptomic somatic alterations of hepatocellular carcinoma in non-cirrhotic livers

Background: Liver cancer is the second leading cause of cancer-related deaths worldwide. Hepatocellular carcinoma (HCC) risk factors include chronic hepatitis, cirrhosis, and alcohol abuse, whereby tumorigenesis is induced through inflammation and subsequent fibrotic response. However, a subset of HCC arises in non-cirrhotic livers. We characterized the genomic and transcriptomic landscape of non-cirrhotic HCC to identify features underlying the disease's development and progression. Methods: Whole genome and transcriptome sequencing was performed on 30 surgically resectable tumors comprised of primarily of non-cirrhotic HCC and adjacent normal tissue. Using somatic variants, capture reagents were created and employed on an additional 87 cases of mixed cirrhotic/non-cirrhotic HCC. Cases were analyzed to identify viral integrations, single nucleotide variants (SNVs), insertions and deletions (INDELS), copy number variants, loss of heterozygosity, gene fusions, structural variants, and differential gene expression. Results: We detected 3,750 SNVs/INDELS and extensive CNVs and expression changes. Recurrent TERT promoter mutations occurred in >52% of non-cirrhotic discovery samples. Frequently mutated genes included TP53, CTNNB1, and APOB. Cytochrome P450 mediated metabolism was significantly downregulated. Structural variants were observed at MACROD2, WDPCP, and NCKAP5 in >20% of samples. Furthermore, NR1H4 fusions involving gene partners EWSR1, GNPTAB, and FNIP1 were detected and validated in 2 non-cirrhotic samples. Conclusion: Genomic analysis can elucidate mechanisms that may contribute to non-cirrhotic HCC tumorigenesis. The comparable mutational landscape between cirrhotic and non-cirrhotic HCC supports previous work suggesting a convergence at the genomic level during disease progression. It is therefore possible genomic-based treatments can be applied to both HCC subtypes with progressed disease.