Rare Single Nucleotide Variants Research Articles

Using coding and non-coding rare variants to target candidate genes in patients with severe tinnitus

Tinnitus is the phantom percept of an internal non-verbal set of noises and tones. It is reported by 15% of the population and it is usually associated with hearing and/or brain disorders. The role of structural variants (SVs) in coding and non-coding regions has not been investigated in patients with severe tinnitus. In this study, we performed whole-genome sequencing in 97 unrelated Swedish individuals with chronic tinnitus (TIGER cohort). Rare single nucleotide variants (SNV), large structural variants (LSV), and copy number variations (CNV) were retrieved to perform a gene enrichment analysis in TIGER and in a subgroup of patients with severe tinnitus (SEVTIN, n = 34), according to the tinnitus handicap inventory (THI) scores. An independent exome sequencing dataset of 147 Swedish tinnitus patients was used as a replication cohort (JAGUAR cohort) and population-specific datasets from Sweden (SweGen) and Non-Finish Europeans (NFE) from gnomAD were used as control groups. SEVTIN patients showed a higher prevalence of hyperacusis, hearing loss, and anxiety when they were compared to individuals in the TIGER cohort. We found an enrichment of rare missense variants in 6 and 8 high-constraint genes in SEVTIN and TIGER cohorts, respectively. Of note, an enrichment of missense variants was found in the CACNA1E gene in both SEVTIN and TIGER. We replicated the burden of missense variants in 9 high-constrained genes in the JAGUAR cohort, including the gene NAV2, when data were compared with NFE. Moreover, LSVs in constrained regions overlapping CACNA1E, NAV2, and TMEM132D genes were observed in TIGER and SEVTIN.

MetaRNN: differentiating rare pathogenic and rare benign missense SNVs and InDels using deep learning

Multiple computational approaches have been developed to improve our understanding of genetic variants. However, their ability to identify rare pathogenic variants from rare benign ones is still lacking. Using context annotations and deep learning methods, we present pathogenicity prediction models, MetaRNN and MetaRNN-indel, to help identify and prioritize rare nonsynonymous single nucleotide variants (nsSNVs) and non-frameshift insertion/deletions (nfINDELs). We use independent test sets to demonstrate that these new models outperform state-of-the-art competitors and achieve a more interpretable score distribution. Importantly, prediction scores from both models are comparable, enabling easy adoption of integrated genotype-phenotype association analysis methods. All pre-computed nsSNV scores are available at http://www.liulab.science/MetaRNN. The stand-alone program is also available at https://github.com/Chang-Li2019/MetaRNN.

High-coverage whole-genome sequencing of the expanded 1000 Genomes Project cohort including 602 trios

SummaryThe 1000 Genomes Project (1kGP) is the largest fully open resource of whole-genome sequencing (WGS) data consented for public distribution without access or use restrictions. The final, phase 3 release of the 1kGP included 2,504 unrelated samples from 26 populations and was based primarily on low-coverage WGS. Here, we present a high-coverage 3,202-sample WGS 1kGP resource, which now includes 602 complete trios, sequenced to a depth of 30X using Illumina. We performed single-nucleotide variant (SNV) and short insertion and deletion (INDEL) discovery and generated a comprehensive set of structural variants (SVs) by integrating multiple analytic methods through a machine learning model. We show gains in sensitivity and precision of variant calls compared to phase 3, especially among rare SNVs as well as INDELs and SVs spanning frequency spectrum. We also generated an improved reference imputation panel, making variants discovered here accessible for association studies.

Extreme purifying selection against point mutations in the human genome

Large-scale genome sequencing has enabled the measurement of strong purifying selection in protein-coding genes. Here we describe a new method, called ExtRaINSIGHT, for measuring such selection in noncoding as well as coding regions of the human genome. ExtRaINSIGHT estimates the prevalence of “ultraselection” by the fractional depletion of rare single-nucleotide variants, after controlling for variation in mutation rates. Applying ExtRaINSIGHT to 71,702 whole genome sequences from gnomAD v3, we find abundant ultraselection in evolutionarily ancient miRNAs and neuronal protein-coding genes, as well as at splice sites. By contrast, we find much less ultraselection in other noncoding RNAs and transcription factor binding sites, and only modest levels in ultraconserved elements. We estimate that ~0.4–0.7% of the human genome is ultraselected, implying ~ 0.26–0.51 strongly deleterious mutations per generation. Overall, our study sheds new light on the genome-wide distribution of fitness effects by combining deep sequencing data and classical theory from population genetics.

The genetic architecture of schizophrenia: review of large-scale genetic studies.

Schizophrenia is a complex and often chronic psychiatric disorder with high heritability. Diagnosis of schizophrenia is still made clinically based on psychiatric symptoms; no diagnostic tests or biomarkers are available. Pathophysiology-based diagnostic scheme and treatments are also not available. Elucidation of the pathogenesis is needed for development of pathology-based diagnostics and treatments. In the past few decades, genetic research has made substantial advances in our understanding of the genetic architecture of schizophrenia. Rare copy number variations (CNVs) and rare single-nucleotide variants (SNVs) detected by whole-genome CNV analysis and whole-genome/-exome sequencing analysis have provided the great advances. Common single-nucleotide polymorphisms (SNPs) detected by large-scale genome-wide association studies have also provided important information. Large-scale genetic studies have been revealed that both rare and common genetic variants play crucial roles in this disorder. In this review, we focused on CNVs, SNVs, and SNPs, and discuss the latest research findings on the pathogenesis of schizophrenia based on these genetic variants. Rare variants with large effect sizes can provide mechanistic hypotheses. CRISPR-based genetics approaches and induced pluripotent stem cell technology can facilitate the functional analysis of these variants detected in patients with schizophrenia. Recent advances in long-read sequence technology are expected to detect variants that cannot be detected by short-read sequence technology. Various studies that bring together data from common variant and transcriptomic datasets provide biological insight. These new approaches will provide additional insight into the pathophysiology of schizophrenia and facilitate the development of pathology-based therapeutics.

Targeted Next-Generation Sequencing of Thymic Epithelial Tumours Revealed Pathogenic Variants in KIT, ERBB2, KRAS, and TP53 in 30% of Thymic Carcinomas.

Simple SummaryThe biology of thymic epithelial tumours (TETs), including thymomas and thymic carcinomas, and particularly the extent of molecular dysregulation, is poorly understood. Through next-generation sequencing of 15 genes implicated in common solid tumours in 53 TETs, we found a larger number of single nucleotide variants (SNVs) in thymic carcinomas than thymomas. About 30% of thymic carcinomas had at least one somatic pathogenic gene variant in TP53, ERBB2, KIT, or KRAS, whereas variants of uncertain clinical significance in KIT, ERBB2, and FOXL2 were found exclusively in thymomas. The presence of somatic pathogenic variants was non-significantly associated with shorter disease-free survival in thymic carcinoma patients. No somatic pathogenic or likely pathogenic SNV was found in thymomas. Importantly, we also evaluated germline SNVs, adding to the number of known genetic alterations in TETs and thereby enhancing our molecular understanding of these neoplasms.A better understanding of the molecular pathogenesis of thymic epithelial tumours (TETs) could revolutionise their treatment. We evaluated thymomas and thymic carcinomas by next-generation sequencing (NGS) of somatic or germline single nucleotide variants (SNVs) in genes commonly mutated in solid tumours. In total, 19 thymomas and 34 thymic carcinomas were analysed for nonsynonymous SNVs in 15 genes by targeted NGS (reference genome: hg19/GRCh37). Ten SNVs in TP53 (G154V, R158P, L194H, R267fs, R273C, R306 *, Q317 *), ERBB2 (V773M), KIT (L576P), and KRAS (Q61L) considered somatic and pathogenic/likely pathogenic were detected in 10 of 34 (29.4%) thymic carcinomas. No somatic SNVs confirmed as pathogenic/likely pathogenic were found in thymomas. Rare SNVs of uncertain or unknown functional and clinical significance, to our knowledge not reported previously in TETs, were found in ERBB2 (S703R), KIT (I690V), and FOXL2 (P157S) in 3 of 19 (16%) thymomas. The most frequent germline SNVs were TP53 P72R (94% TETs), ERBB2 I655V (40% TETs), and KIT M541L (9% TETs). No significant difference in median disease-free survival (DFS) was found between thymic carcinoma patients with and without pathogenic SNVs (p = 0.190); however, a trend toward a longer DFS was observed in the latter (16.0 vs. 30.0 months, respectively). In summary, NGS analysis of TETs revealed several SNVs in genes related to the p53, AKT, MAPK, and K-Ras signalling pathways. Thymic carcinomas showed greater genetic dysregulation than thymomas. The germline and rare SNVs of uncertain clinical significance reported in this study add to the number of known genetic alterations in TETs, thus extending our molecular understanding of these neoplasms. Druggable KIT alterations in thymic carcinomas have potential as therapeutic targets.

Novel KEL allele associated with loss of Kpb identified in a white blood donor.

The importance of identifying variant alleles among blood donors is significant to the safety of transfusion for recipients. Molecular methods have become more prominent in the routine process of antigen typing donor units. Some variant antigens cannot be detected using only serologic methods. Molecular testing allows the determination of nucleotide sequences that are used to predict a phenotype. Antigens of the Kell blood group system are known for being highly immunogenic and causing adverse reactions upon antibody formation. A female white blood donor who typed Kp(b-) using serologic methods on multiple donations since 2005 was the subject of a typing discrepancy investigation. Routine genotyping using a commercial genotyping kit (HemoID DQS Panel; Agena Bioscience, San Diego, CA) predicted the donor to type Kp(a+b+). Investigation of the discrepancy between these two results identified a rare single nucleotide variant in the KEL gene at nucleotide position c.948G>T that alters amino acid residue 316 from tryptophan (Trp) to cysteine (Cys). After discovery of the novel allele, adsorption and elution studies were performed to see if there was weakened Kpb expression. The elution studies yielded negative results, which indicated that Kpb is not expressed. The KEL transcripts expressed by the donor were determined using cDNA analysis, and the predicted amino acid sequence of the novel allele was modeled to investigate the impact of the amino acid sequence on the structure of the KEL polypeptide. Both SWISS-MODEL and Robetta software were used to evaluate the impact of the p.Trp316Cys on the three-dimensional protein structure. There was no conformational change noted with SWISS-MODEL, whereas the Robetta software showed a significant conformational change compared with the normal Kp(b+) reference sequence. Because the donor is homozygous for variants associated with k and Jsb expression, it was not possible to determine whether the novel allele is associated with loss of Kpb only or loss of all Kell antigens.

Two ovarian candidate enhancers, identified by time series enhancer RNA analyses, harbor rare genetic variations identified in ovarian insufficiency.

The genetic regulation of ovarian development remains largely unclear. Indeed, in most cases of impaired ovarian development-such as 46,XX disorders of sex development (DSD) without SRY, and premature ovarian insufficiency (POI)-the genetic causes have not been identified, and the vast majority of disease-associated sequence variants could lie within non-coding regulatory sequences. In this study, we aimed to identify enhancers of five ovarian genes known to play key roles in early ovarian development, basing our analysis on the expression of enhancer derived transcripts (eRNAs), which are considered to characterize active enhancers. Temporal expression profile changes in mouse WT1-positive ovarian cells were obtained from cap analysis of gene expression at E13.5, E16.5 and P0. We compared the chronological expression profiles of ovarian-specific eRNA with expression profiles for each of the ovarian-specific genes, yielding two candidate sequences for enhancers of Wnt4 and Rspo1. Both sequences are conserved between mouse and human, and we confirmed their enhancer activities using transient expression assays in murine granulosa cells. Furthermore, by sequencing the region in patients with impaired ovarian development in 24 patients, such as POI, gonadal dysgenesis and 46,XX DSD, we identified rare single nucleotide variants in both sequences. Our results demonstrate that combined analysis of the temporal expression profiles of eRNA and mRNA of target genes presents a powerful tool for locating cis-element enhancers, and a means of identifying disease-associated sequence variants that lie within non-coding regulatory sequences, thus advancing an important unmet need in forward human genetics.

Role of ADAMTS13, VWF and F8 genes in deep vein thrombosis.

BackgroundWe previously described the association between rare ADAMTS13 single nucleotide variants (SNVs) and deep vein thrombosis (DVT). Moreover, DVT patients with at least one rare ADAMTS13 SNV had a lower ADAMTS13 activity than non-carriers.AimsTo confirm ADAMTS13 variants association with DVT and reduced plasma ADAMTS13 activity levels in a larger population. To investigate the role of VWF and F8 variants.MethodsADAMTS13, VWF and F8 were sequenced using next-generation sequencing in 594 Italian DVT patients and 571 controls. Genetic association testing was performed using logistic regression and gene-based tests. The association between rare ADAMTS13 variants and the respective plasmatic activity, available for 365 cases and 292 controls, was determined using linear regression. All analyses were age-, sex- adjusted.ResultsWe identified 48 low-frequency/common and 272 rare variants. Nine low-frequency/common variants had a P<0.05, but a false discovery rate between 0.06 and 0.24. Of them, 7 were found in ADAMTS13 (rs28641026, rs28503257, rs685523, rs3124768, rs3118667, rs739469, rs3124767; all protective) and 2 in VWF (rs1800382 [risk], rs7962217 [protective]). Rare ADAMTS13 variants were significantly associated with DVT using the burden, variable threshold (VT) and UNIQ (P<0.05), but not with C-ALPHA, SKAT and SKAT-O tests. Rare VWF and F8 variants were not associated with DVT. Carriers of rare ADAMTS13 variants had lower ADAMTS13 activity than non-carriers (ß -6.2, 95%CI -11,-1.5). This association was stronger for DVT patients than controls (ß -7.5, 95%CI -13.5,-1.5 vs. ß -2.9, 95%CI -10.4,4.5).ConclusionsADAMTS13 and VWF low-frequency/common variants mainly showed a protective effect, although their association with DVT was not confirmed. DVT patients carrying a rare ADAMTS13 variants had slightly reduced ADAMTS13 activity levels, but a higher DVT risk. Rare VWF and FVIII variants were not associated with DVT suggesting that other mechanisms are responsible for the high VWF and FVIII levels measured in DVT patients.

Rare single nucleotide variants in COL5A1 promoter do not play a major role in keratoconus susceptibility associated with rs1536482

BackgroundKeratoconus is a chronic degenerative disorder of the cornea characterized by thinning and cone-shaped protrusions. Although genetic factors play a key role in keratoconus development, the etiology is still under investigation. The occurrence of single-nucleotide polymorphisms (SNPs) associated with keratoconus in Russian patients is poorly studied. The purpose of this study was to validate whether three reported keratoconus-associated SNPs (rs1536482 near the COL5A1 gene, rs2721051 near the FOXO1 gene, rs1324183 near the MPDZ gene) are also actual for a Russian cohort of patients. Additionally, we investigated the COL5A1 promoter sequence for single-nucleotide variants (SNVs) in a subgroup of keratoconus patients with at least one rs1536482 minor allele (rs1536482+) to assess the role of these SNVs in keratoconus susceptibility associated with rs1536482.MethodsThis case-control study included 150 keratoconus patients and two control groups (main and additional, 205 and 474 participants, respectively). We performed PCR targeting regions flanking SNVs and the COL5A1 promoter, followed by Sanger sequencing of amplicons. The additional control group was genotyped using an SNP array.ResultsThe minor allele frequency was significantly different between the keratoconus and control cohorts (main and combined) for rs1536482, rs2721051, and rs1324183 (p-value < 0.05). The rare variants rs1043208782 and rs569248712 were found in the COL5A1 promoter in two out of 94 rs1536482+ keratoconus patients.Conclusionrs1536482, rs2721051, and rs1324183 were associated with keratoconus in a Russian cohort. SNVs in the COL5A1 promoter do not play a major role in keratoconus susceptibility associated with rs1536482.

Genome sequencing identifies rare tandem repeat expansions and copy number variants in Lennox-Gastaut syndrome.

Epilepsies are a group of common neurological disorders with a substantial genetic basis. Despite this, the molecular diagnosis of epilepsies remains challenging due to its heterogeneity. Studies utilizing whole-genome sequencing may provide additional insights into genetic causes of epilepsies of unknown aetiology. Whole-genome sequencing was used to evaluate a cohort of adults with unexplained developmental and epileptic encephalopathies (n = 30), for whom prior genetic tests, including whole-exome sequencing in some cases, were negative or inconclusive. Rare single nucleotide variants, insertions/deletions, copy number variants and tandem repeat expansions were analysed. Seven pathogenic or likely pathogenic single nucleotide variants, and two pathogenic deleterious copy number variants were identified in nine patients (32.1% of the cohort). One of the copy number variants, identified in a patient with Lennox–Gastaut syndrome, was too small to be detected by chromosomal microarray techniques. We also identified two tandem repeat expansions with clinical implications in two other patients with Lennox–Gastaut syndrome: a CGG repeat expansion in the 5′untranslated region of DIP2B, and a CTG expansion in ATXN8OS (previously implicated in spinocerebellar ataxia type 8). Three patients had KCNA2 pathogenic variants. One of them died of sudden unexpected death in epilepsy. The other two patients had, in addition to a KCNA2 variant, a second de novo variant impacting potential epilepsy-relevant genes (KCNIP4 and UBR5). Overall, whole-genome sequencing provided a genetic explanation in 32.1% of the total cohort. This is also the first report of coding and non-coding tandem repeat expansions identified in patients with Lennox–Gastaut syndrome. This study demonstrates that using whole-genome sequencing, the examination of multiple types of rare genetic variation, including those found in the non-coding region of the genome, can help resolve unexplained epilepsies.

Personalized Genetic Diagnosis of Congenital Heart Defects in Newborns.

Congenital heart disease is a group of pathologies characterized by structural malformations of the heart or great vessels. These alterations occur during the embryonic period and are the most frequently observed severe congenital malformations, the main cause of neonatal mortality due to malformation, and the second most frequent congenital malformations overall after malformations of the central nervous system. The severity of different types of congenital heart disease varies depending on the combination of associated anatomical defects. The causes of these malformations are usually considered multifactorial, but genetic variants play a key role. Currently, use of high-throughput genetic technologies allows identification of pathogenic aneuploidies, deletions/duplications of large segments, as well as rare single nucleotide variants. The high incidence of congenital heart disease as well as the associated complications makes it necessary to establish a diagnosis as early as possible to adopt the most appropriate measures in a personalized approach. In this review, we provide an exhaustive update of the genetic bases of the most frequent congenital heart diseases as well as other syndromes associated with congenital heart defects, and how genetic data can be translated to clinical practice in a personalized approach.

Whole-exome sequencing reveals insights into genetic susceptibility to Congenital Zika Syndrome.

Congenital Zika Syndrome (CZS) is a critical illness with a wide range of severity caused by Zika virus (ZIKV) infection during pregnancy. Life-threatening neurodevelopmental dysfunctions are among the most common phenotypes observed in affected newborns. Risk factors that contribute to susceptibility and response to ZIKV infection may be related to the virus itself, the environment, and maternal genetic background. Nevertheless, the newborn’s genetic contribution to the critical illness is still not elucidated. Here, we aimed to identify possible genetic variants as well as relevant biological pathways that might be associated with CZS phenotypes. For this purpose, we performed a whole-exome sequencing in 40 children born to women with confirmed exposure to ZIKV during pregnancy. We investigated the occurrence of rare harmful single-nucleotide variants (SNVs) possibly associated with inborn errors in genes ontologically related to CZS phenotypes. Moreover, an exome-wide association analysis was also performed using a case-control design (29 CZS cases and 11 controls), for both common and rare variants. Five out of the 29 CZS patients harbored known pathogenic variants likely to contribute to mild to severe manifestations observed. Approximately, 30% of affected individuals carried at least one pathogenic or likely pathogenic SNV in genes candidates to play a role in CZS. Our common variant association analysis detected a suggestive protective effect of the rs2076469 in DISP3 gene (p-value: 1.39 x 10−5). The IL12RB2 gene (p-value: 2.18x10-11) also showed an unusual distribution of nonsynonymous rare SNVs in control samples. Finally, genes harboring harmful variants are involved in processes related to CZS phenotypes such as neurological development and immunity. Therefore, both rare and common variations may be likely to contribute as the underlying genetic cause of CZS susceptibility. The variations and pathways identified in this study may also have implications for the development of therapeutic strategies in the future.

Missense RHD single nucleotide variants induce weakened D antigen expression by altering splicing and/or protein expression.

Although D variant phenotype is known to be due to genetic defects, including rare missense single nucleotide variants (SNVs), within the RHD gene, few studies have addressed the molecular and cellular mechanisms driving this altered expression. We and others showed previously that splicing is commonly disrupted by SNVs in constitutive splice sites and their vicinity. We thus sought to investigate whether rare missense SNVs located in "deep" exonic regions could also impair this mechanism. Forty-six missense SNVs reported within exons 6 and 7 were first selected from the Human RhesusBase. Their respective effect on splicing was assessed by using an in vitro assay. An RhD-negative cell model was further generated by using the CRISPR-Cas9 approach. RhD-mutated proteins were overexpressed in the newly created model, and cell membrane expression of the D antigen was measured by flow cytometry. Minigene splicing assay showed that 14 of 46 (30.4%) missense SNVs alter splicing. Very interestingly, further investigation of two missense SNVs, which both affect codon 338 and confer a weak D phenotype, showed various mechanisms: c.1012C>G (p.Leu338Val) disrupts splicing only, while c.1013T>C (p.Leu338Pro) alters only the protein structure, in agreement with in silico prediction tools and 3D protein structure visualization. Our functional data set suggests that missense SNVs damage quantitatively D antigen expression by, at least, two different mechanisms (splicing alteration and protein destabilization) that may act independently. These data thereby contribute to extend the current knowledge of the molecular mechanisms governing weakened D expression.

Whole-genome sequencing identifies functional noncoding variation in SEMA3C that cosegregates with dyslexia in a multigenerational family

Dyslexia is a common heritable developmental disorder involving impaired reading abilities. Its genetic underpinnings are thought to be complex and heterogeneous, involving common and rare genetic variation. Multigenerational families segregating apparent monogenic forms of language-related disorders can provide useful entrypoints into biological pathways. In the present study, we performed a genome-wide linkage scan in a three-generational family in which dyslexia affects 14 of its 30 members and seems to be transmitted with an autosomal dominant pattern of inheritance. We identified a locus on chromosome 7q21.11 which cosegregated with dyslexia status, with the exception of two cases of phenocopy (LOD = 2.83). Whole-genome sequencing of key individuals enabled the assessment of coding and noncoding variation in the family. Two rare single-nucleotide variants (rs144517871 and rs143835534) within the first intron of the SEMA3C gene cosegregated with the 7q21.11 risk haplotype. In silico characterization of these two variants predicted effects on gene regulation, which we functionally validated for rs144517871 in human cell lines using luciferase reporter assays. SEMA3C encodes a secreted protein that acts as a guidance cue in several processes, including cortical neuronal migration and cellular polarization. We hypothesize that these intronic variants could have a cis-regulatory effect on SEMA3C expression, making a contribution to dyslexia susceptibility in this family.

Contribution of Copy Number Variation in Idiopathic Hypogonadotropic Hypogonadism

Introduction: While the role of single nucleotide variants (SNVs) in causal genes for Idiopathic Hypogonadotropic hypogonadism (IHH) is known, the contribution of copy number variants (CNVs) to IHH has not been systematically studied. Here, we examined the prevalence of CNVs in a large IHH cohort and their associated phenotypic spectrum. Methods: Exome sequencing (ES) from 1,441 IHH probands was analyzed using the GATK-gCNV pipeline to identify CNVs that spanned known IHH genes with a site frequency of <1%. ES data was also analyzed for rare SNVs in all IHH genes. IHH subjects were evaluated for reproductive and non-reproductive phenotypes (kidney, eye, cardiovascular, midline/ facial, bone and neurological abnormalities). Results: (i) CNV prevalence in IHH: Three percent of the IHH probands (46/1441 probands, 38 males and 8 females) harbored CNVs in 19/36 high confidence IHH genes (26 deletions 20 duplications). The vast majority of CNVs disrupted either the ANOS1 (26%) or FGFR1 (17%) genes. Intriguingly, CHD7 (a gene that carries SNVs in ~12% of our IHH subjects) did not harbor any CNVs. (ii) Phenotypic analysis: We found that CNVs were more common in IHH subjects with anosmia (Kallmann syndrome, N=33) compared to normosmic IHH (N=13). More than half of the subjects with CNVs affecting IHH genes carried at least 1 non-reproductive feature (26/46 IHH subjects, 56%). A syndromic presentation with multiple non-reproductive phenotypes was more common in IHH subjects harboring multigenic CNVs (IHH gene & additional genes) compared to IHH subjects with single IHH gene CNV or SNV in a single IHH gene (61% vs. 12% or 18%, respectively; p 0.0006). Conclusions: CNVs in known IHH genes contribute to 3% of the IHH genetic architecture in our cohort. IHH subjects with larger multigenic CNVs displayed phenotypes consistent with contiguous gene syndromes. The absence of genic CNVs in genes frequently found to carry SNVs, such as CHD7, requires additional analysis to establish the biologic or mechanistic reasons for their rarity.

A comprehensive system for detecting rare single nucleotide variants based on competitive DNA probe and duplex-specific nuclease

Single nucleotide variants (SNVs) have emerged as increasingly important biomarkers, particularly in the diagnosis and prognosis of cancers. However, most SNVs are rarely detected in blood samples from cancer patients as they are surrounded by abundant concomitant wild-type nucleic acids. Herein, we design a system that features a combination of competitive DNA probe system (CDPS) and duplex-specific nuclease (DSN) that we referred to as CAD. A theoretical model was established for the CAD system based on reaction networks. Guided by the theoretical model, we found that a minor loss in sensitivity significantly improved the specificity of the system, thus creating a theoretical discrimination factor (DF) > 100 for most conditions. This non-equivalent tradeoff between sensitivity and specificity provides a new concept for the analysis of rare DNA-sequence variants. As a demonstration of practicality, we applied as-proposed CAD system to identify low variant allele frequency (VAF) in a synthetic template (0.1% VAF) and human genomic DNA (1% VAF). This work promises complete guidance for the design of enzyme-based nucleic acid analysis.

Familial Psychosis Associated With a Missense Mutation at MACF1 Gene Combined With the Rare Duplications DUP3p26.3 and DUP16q23.3, Affecting the CNTN6 and CDH13 Genes.

Psychosis is a highly heritable and heterogeneous psychiatric condition. Its genetic architecture is thought to be the result of the joint effect of common and rare variants. Families with high prevalence are an interesting approach to shed light on the rare variant’s contribution without the need of collecting large cohorts. To unravel the genomic architecture of a family enriched for psychosis, with four affected individuals, we applied a system genomic approach based on karyotyping, genotyping by whole-exome sequencing to search for rare single nucleotide variants (SNVs) and SNP array to search for copy-number variants (CNVs). We identified a rare non-synonymous variant, g.39914279 C > G, in the MACF1 gene, segregating with psychosis. Rare variants in the MACF1 gene have been previously detected in SCZ patients. Besides, two rare CNVs, DUP3p26.3 and DUP16q23.3, were also identified in the family affecting relevant genes (CNTN6 and CDH13, respectively). We hypothesize that the co-segregation of these duplications with the rare variant g.39914279 C > G of MACF1 gene precipitated with schizophrenia and schizoaffective disorder.

PLXNA2 and LRRC40 as candidate genes in autism spectrum disorder.

Autism spectrum disorder (ASD) is a neurodevelopmental disability with high heritability yet the genetic etiology remains elusive. Therefore, it is necessary to elucidate new genotype-phenotype relationships for ASD to improve both the etiological knowledge and diagnosis. In this work, a copy-number variant and whole-exome sequencing analysis were performed in an ASD patient with a complex neurobehavioral phenotype with epilepsy and attention deficit hyperactivity disorder. We identified rare recessive single nucleotide variants in the two genes, PLXNA2 encoding Plexin A2 that participates in neurodevelopment, and LRRC40, which encodes Leucine-rich repeat containing protein 40, a protein of unknown function. PLXNA2 showed the heterozygous missense variants c.614G>A (p.Arg205Gln) and c.4904G>A (p.Arg1635Gln) while LRRC40 presented the homozygous missense variant c.1461G>T (p.Leu487Phe). In silico analysis predicted that these variants could be pathogenic. We studied PLXNA2 and LRRC40 mRNA and proteins in fibroblasts from the patient and controls. We observed a significant PlxnA2 subcellular delocalization and very low levels of LRRC40 in the patient. Moreover, we found a novel interaction between PlxnA2 and LRRC40 suggesting that participate in a common neural pathway. This interaction was significant decreased in the patient's fibroblasts. In conclusion, our results identified PLXNA2 and LRRC40 genes as candidates in ASD providing novel clues for the pathogenesis. Further attention to these genes is warranted in genetic studies of patients with neurodevelopmental disorders, particularly ASD. LAY SUMMARY: Genomics is improving the knowledge and diagnosis of patients with autism spectrum disorder (ASD) yet the genetic etiology remains elusive. Here, using genomic analysis together with experimental functional studies, we identified in an ASD complex patient the PLXNA2 and LRRC40 recessive genes as ASD candidates. Furthermore, we found that the proteins of these genes interact in a common neural network. Therefore, more attention to these genes is warranted in genetic studies of patients with neurodevelopmental disorders, particularly ASD.

Allelic Dropout Is a Common Phenomenon That Reduces the Diagnostic Yield of PCR-Based Sequencing of Targeted Gene Panels.

Primary cardiomyopathies (CMPs) are monogenic but multi-allelic disorders with dozens of genes involved in pathogenesis. The implementation of next-generation sequencing (NGS) approaches has resulted in more time- and cost-efficient DNA diagnostics of cardiomyopathies. However, the diagnostic yield of genetic testing for each subtype of CMP fails to exceed 60%. The aim of this study was to demonstrate that allelic dropout (ADO) is a common phenomenon that reduces the diagnostic yield in primary cardiomyopathy genetic testing based on targeted gene panels assayed on the Ion Torrent platform. We performed mutational screening with three custom targeted gene panels based on sets of oligoprimers designed automatically using AmpliSeq Designer® containing 1049 primer pairs for 37 genes with a total length of 153 kb. DNA samples from 232 patients were screened with at least one of these targeted gene panels. We detected six ADO events in both IonTorrent PGM (three cases) and capillary Sanger sequencing (three cases) data, identifying ADO-causing variants in all cases. All ADO events occurred due to common or rare single nucleotide variants (SNVs) in the oligoprimer binding sites and were detected because of the presence of “marker” SNVs in the target DNA fragment. We ultimately identified that PCR-based NGS involves a risk of ADO that necessitates the use of Sanger sequencing to validate NGS results. We assume that oligoprimer design without ADO data affects the amplification efficiency of up to 0.77% of amplicons.