Deep intronic mutations and human disease

Polymorphisms in JAK2 Gene are Associated with Production Traits and Mastitis Resistance in Dairy Cattle

The present study was designed to investigate the effects of single nucleotide polymorphisms (SNPs) in the JAK2 gene on the production and mastitis related traits in dairy cattle. Blood and milk samples were collected from 201 lactating dairy cattle of three breeds, i.e. Holstein Friesian (HF), Jersey (J) and Achai (A) and their crosses maintained at well-established dairy farms in Khyber Pakhtunkhwa, Pakistan. Generalized linear model was used to evaluate the association between genotypes and the studied traits. A DNA pool was made from randomly selected 30 samples which revealed three SNPs, i.e. SNP 1 in 5’ upstream region (G>A, rs379754157), SNP 2 in intron 15 (A>G, rs134192265), and SNP 3 in exon 20 (A>G, rs110298451) that were further validated in the population under study using SNaPshot technique. Of the three SNPs, SNP 1 did not obey Hardy-Weinberg equilibrium (P<0.05). SNP 2 and SNP 3 were found to be in strong linkage disequilibrium and allele G was highly prevalent compared to allele A in these SNPs. in SNP 1, the GG genotype was associated with significantly (P<0.01) higher SCC, whereas SNP 2 and SNP 3 were significantly (P<0.01) associated with higher lactose percentage compared to the other geno-types. The haplogroups association analysis revealed that H1H2 (GG GG AG) has significantly lower SCC than H2H2 (GG GG GG). The results infer that JAK2 could be an important candidate gene and the studied SNPs might be useful genetic markers for production and mastitis related traits.

A combined RNA-seq and whole genome sequencing approach for identification of non-coding pathogenic variants in single families

Inherited retinal degenerations (IRDs) are at the focus of current genetic therapeutic advancements. For a genetic treatment such as gene therapy to be successful, an accurate genetic diagnostic is required. Genetic diagnostics relies on the assessment of the probability that a given DNA variant is pathogenic. Non-coding variants present a unique challenge for such assessments as compared to coding variants. For one, non-coding variants are present at much higher number in the genome than coding variants. In addition, our understanding of the rules that govern the non-coding regions of the genome is less complete than our understanding of the coding regions. Methods that allow for both the identification of candidate non-coding pathogenic variants and their functional validation may help overcome these caveats allowing for a greater number of patients to benefit from advancements in genetic therapeutics. We present here an unbiased approach combining whole genome sequencing (WGS) with patient-induced pluripotent stem cell (iPSC)-derived retinal organoids (ROs) transcriptome analysis. With this approach, we identified and functionally validated a novel pathogenic non-coding variant in a small family with a previously unresolved genetic diagnosis.

A novel 23 bp indel mutation in PRL gene is associated with growth traits in Luxi Blackhead sheep

Prolactin is a highly versatile pituitary hormone with multiple biological functions. PRL expression is regulated by POU1F1 and the prophet of POU1F1 (PROP1). The aim of this study was to investigate the indel variations in ovine PRL and the directly related (PROP1 and POU1F1) genes, and their associations with growth traits in Luxi Blackhead (LXBH) sheep. A monomorphism in PROP1 and POU1F1 genes, and one novel 23-bp insertion mutation in the PRL gene, were identified in LXBH sheep. The 23 bp insertion mutation within PRL gene was significantly associated with several body measurements (e.g., body weight, body height) in sheep of different ages (p < 0.05). Ram lambs (p = 0.036) of genotype insertion/insertion (II) had significantly higher body weights. Weaners (p = 0.018) of genotypes insertion/insertion (II) and insertion/deletion (ID) also had significantly higher body weights compared with male sheep of deletion/deletion (DD) genotype. Moreover, among ewe lambs, individuals of genotype insertion/insertion (II) had a higher paunch girth compared to those with other genotypes (p = 0.044). These findings indicate that a 23 bp indel variant of the ovine PRL gene is correlated with body measurements in LXBH sheep. The findings have potential utility for sheep breeding programs based on marker-assisted selection.

Increasing the Genetic Diagnosis Yield in Inherited Retinal Dystrophies: Assigning Pathogenicity to Novel Non-canonical Splice Site Variants

Aims: We aimed to validate the pathogenicity of genetic variants identified in inherited retinal dystrophy (IRD) patients, which were located in non-canonical splice sites (NCSS). Methods: After next generation sequencing (NGS) analysis (target gene panels or whole exome sequencing (WES)), NCSS variants were prioritized according to in silico predictions. In vivo and in vitro functional tests were used to validate their pathogenicity. Results: Four novel NCSS variants have been identified. They are located in intron 33 and 34 of ABCA4 (c.4774-9G>A and c.4849-8C>G, respectively), intron 2 of POC1B (c.101-3T>G) and intron 3 of RP2 (c.884-14G>A). Functional analysis detected different aberrant splicing events, including intron retention, exon skipping and intronic nucleotide addition, whose molecular effect was either the disruption or the elongation of the open reading frame of the corresponding gene. Conclusions: Our data increase the genetic diagnostic yield of IRD patients and expand the landscape of pathogenic variants, which will have an impact on the genotype–phenotype correlations and allow patients to opt for the emerging gene and cell therapies.

The Experimentally Obtained Functional Impact Assessments of 5' Splice Site GT'GC Variants Differ Markedly from Those Predicted

Introduction:

5' splice site GT>GC or +2T>C variants have been frequently reported to cause human genetic disease and are routinely scored as pathogenic splicing mutations. However, we have recently demonstrated that such variants in human disease genes may not invariably be pathogenic. Moreover, we found that no splicing prediction tools appear to be capable of reliably distinguishing those +2T>C variants that generate wild-type transcripts from those that do not.

Methodology

Herein, we evaluated the performance of a novel deep learning-based tool, SpliceAI, in the context of three datasets of +2T>C variants, all of which had been characterized functionally in terms of their impact on pre-mRNA splicing. The first two datasets refer to our recently described “in vivo” dataset of 45 known disease-causing +2T>C variants and the “in vitro” dataset of 103 +2T>C substitutions subjected to full-length gene splicing assay. The third dataset comprised 12 BRCA1 +2T>C variants that were recently analyzed by saturation genome editing.

Results

Comparison of the SpliceAI-predicted and experimentally obtained functional impact assessments of these variants (and smaller datasets of +2T>A and +2T>G variants) revealed that although SpliceAI performed rather better than other prediction tools, it was still far from perfect. A key issue was that the impact of those +2T>C (and +2T>A) variants that generated wild-type transcripts represents a quantitative change that can vary from barely detectable to an almost full expression of wild-type transcripts, with wild-type transcripts often co-existing with aberrantly spliced transcripts.

Conclusion

Our findings highlight the challenges that we still face in attempting to accurately identify splice-altering variants.

Recent advances in primary immunodeficiency: from molecular diagnosis to treatment

The technological advances in diagnostics and therapy of primary immunodeficiency are progressing at a fast pace. This review examines recent developments in the field of inborn errors of immunity, from their definition to their treatment. We will summarize the challenges posed by the growth of next-generation sequencing in the clinical setting, touch briefly on the expansion of the concept of inborn errors of immunity beyond the classic immune system realm, and finally review current developments in targeted therapies, stem cell transplantation, and gene therapy.

Inherited Renal Tubulopathies-Challenges and Controversies

Electrolyte homeostasis is maintained by the kidney through a complex transport function mostly performed by specialized proteins distributed along the renal tubules. Pathogenic variants in the genes encoding these proteins impair this function and have consequences on the whole organism. Establishing a genetic diagnosis in patients with renal tubular dysfunction is a challenging task given the genetic and phenotypic heterogeneity, functional characteristics of the genes involved and the number of yet unknown causes. Part of these difficulties can be overcome by gathering large patient cohorts and applying high-throughput sequencing techniques combined with experimental work to prove functional impact. This approach has led to the identification of a number of genes but also generated controversies about proper interpretation of variants. In this article, we will highlight these challenges and controversies.

RNA Splicing Defects in Hypertrophic Cardiomyopathy: Implications for Diagnosis and Therapy

Hypertrophic cardiomyopathy (HCM), the most common inherited heart disease, is predominantly caused by mutations in genes that encode sarcomere-associated proteins. Effective gene-based diagnosis is critical for the accurate clinical management of patients and their family members. However, the introduction of high-throughput DNA sequencing approaches for clinical diagnostics has vastly expanded the number of variants of uncertain significance, leading to many inconclusive results that limit the clinical utility of genetic testing. More recently, developments in RNA analysis have been improving diagnostic outcomes by identifying new variants that interfere with splicing. This review summarizes recent discoveries of RNA mis-splicing in HCM and provides an overview of research that aims to apply the concept of RNA therapeutics to HCM.

Fundoscopy-directed genetic testing to re-evaluate negative whole exome sequencing results

Background

Whole exome sequencing (WES) allows for an unbiased search of the genetic cause of a disease. Employing it as a first-tier genetic testing can be favored due to the associated lower incremental cost per diagnosis compared to when using it later in the diagnostic pathway. However, there are technical limitations of WES that can lead to inaccurate negative variant callings. Our study presents these limitations through a re-evaluation of negative WES results using subsequent tests primarily driven by fundoscopic findings. These tests included targeted gene testing, inherited retinal gene panels, whole genome sequencing (WGS), and array comparative genomic hybridization.

Results

Subsequent genetic testing guided by fundoscopy findings identified the following variant types causing retinitis pigmentosa that were not detected by WES: frameshift deletion and nonsense variants in the RPGR gene, 353-bp Alu repeat insertions in the MAK gene, and large exonic deletion variants in the EYS and PRPF31 genes. Deep intronic variants in the ABCA4 gene causing Stargardt disease and the GUCY2D gene causing Leber congenital amaurosis were also identified.

Conclusions

Negative WES analyses inconsistent with the phenotype should raise clinical suspicion. Subsequent genetic testing may detect genetic variants missed by WES and can make patients eligible for gene replacement therapy and upcoming clinical trials. When phenotypic findings support a genetic etiology, negative WES results should be followed by targeted gene sequencing, array based approach or whole genome sequencing.

Analysis of the relationship between MIR155HG variants and gastric Cancer susceptibility

Background

Gastric cancer is one of the most common cancers in the world and a major cause of cancer-related death. This study aims to determine whether genetic variations in MIR155HG could be associated with gastric cancer risk.

Materials & methods

A total of 506 gastric cancer patients and 500 healthy controls were enrolled in this study. Genotypes were examined with the MassARRAY platform and data management and analysis were conducted with the Typer Software. Odds ratios (OR) and 95% confidence intervals (CIs) were calculated with logistic regression adjusting for age and gender to evaluate the associations between SNPs with gastric cancer in genetic model analysis.

Results

The “CC” genotype of rs4143370 decreased the risk of gastric cancer in genotype model (p = 0.020) and recessive model (p = 0.018). Inversely, the “CC” genotype of rs1893650 increased the risk of gastric cancer in genotype model (p = 0.023) and recessive model (p = 0.014). Stratified analysis showed that rs11911469 was associated with an increased risk of gastric cancer only among the male group in the dominant model (p = 0.039) and additive model (p = 0.030). The haplotype analysis showed a strong linkage disequilibrium among these six SNPs (rs4143370, rs77699734, rs11911469, rs1893650, rs34904192 and rs928883).

Conclusion

This study confirmed the relationship between SNPs of MIR155HG and the gastric cancer risk among the Chinese Han population. Our data may provide a new perspective to understand the aetiology of gastric cancer.

The diagnosis of inborn errors of metabolism by an integrative "multi-omics" approach: A perspective encompassing genomics, transcriptomics, and proteomics

Given the rapidly decreasing cost and increasing speed and accessibility of massively parallel technologies, the integration of comprehensive genomic, transcriptomic, and proteomic data into a "multi-omics" diagnostic pipeline is within reach. Even though genomic analysis has the capability to reveal all possible perturbations in our genetic code, analysis typically reaches a diagnosis in just 35% of cases, with a diagnostic gap arising due to limitations in prioritization and interpretation of detected variants. Here we review the utility of complementing genetic data with transcriptomic data and give a perspective for the introduction of proteomics into the diagnostic pipeline. Together these methodologies enable comprehensive capture of the functional consequence of variants, unobtainable by the analysis of each methodology in isolation. This facilitates functional annotation and reprioritization of candidate genes and variants-a promising approach to shed light on the underlying molecular cause of a patient's disease, increasing diagnostic rate, and allowing actionability in clinical practice.

Practical guide to genetic screening for inherited eye diseases

Genetic eye diseases affect around one in 1000 people worldwide for which the molecular aetiology remains unknown in the majority. The identification of disease-causing gene variant(s) allows a better understanding of the disorder and its inheritance. There is now an approved retinal gene therapy for autosomal recessive RPE65-retinopathy, and numerous ocular gene/mutation-targeted clinical trials underway, highlighting the importance of establishing a genetic diagnosis so patients can fully access the latest research developments and treatment options. In this review, we will provide a practical guide to managing patients with these conditions including an overview of inheritance patterns, required pre- and post-test genetic counselling, different types of cytogenetic and genetic testing available, with a focus on next generation sequencing using targeted gene panels, whole exome and genome sequencing. We will expand on the pros and cons of each modality, variant interpretation and options for family planning for the patient and their family. With the advent of genomic medicine, genetic screening will soon become mainstream within all ophthalmology subspecialties for prevention of disease and provision of precision therapeutics.

A deep intronic splice variant advises reexamination of presumably dominant SPG7 Cases

OBJECTIVE

To identify causative mutations in a patient affected by ataxia and spastic paraplegia.

METHODS

Whole-exome sequencing (WES) and whole-genome sequencing (WGS) were performed using patient's DNA sample. RT-PCR and cDNA Sanger sequencing were performed on RNA extracted from patient's fibroblasts, as well as western blot.

RESULTS

A novel missense variant in SPG7 (c.2195T> C; p.Leu732Pro) was first found by whole-exome sequencing (WES), while the second, also unreported, deep intronic variant (c.286 + 853A>G) was identified by whole-genome sequencing (WGS). RT-PCR confirmed the in silico predictions showing that this variant activated a cryptic splice site, inducing the inclusion of a pseudoexon into the mRNA sequence, which encoded a premature stop codon. Western blot showed decreased SPG7 levels in patient's fibroblasts.

INTERPRETATION

Identification of a deep intronic variant in SPG7, which could only have been detected by performing WGS, led to a diagnosis in this HSP patient. This case challenges the notion of an autosomal dominant inheritance for SPG7, and illustrates the importance of performing WGS subsequently or alternatively to WES to find additional mutations, especially in patients carrying one variant in a gene causing a predominantly autosomal recessive disease.

Identification of splice defects due to noncanonical splice site or deep-intronic variants in ABCA4

Pathogenic variants in the ATP-binding cassette transporter A4 (ABCA4) gene cause a continuum of retinal disease phenotypes, including Stargardt disease. Noncanonical splice site (NCSS) and deep-intronic variants constitute a large fraction of disease-causing alleles, defining the functional consequences of which remains a challenge. We aimed to determine the effect on splicing of nine previously reported or unpublished NCSS variants, one near exon splice variant and nine deep-intronic variants in ABCA4, using in vitro splice assays in human embryonic kidney 293T cells. Reverse transcription-polymerase chain reaction and Sanger sequence analysis revealed splicing defects for 12 out of 19 variants. Four deep-intronic variants create pseudoexons or elongate the upstream exon. Furthermore, eight NCSS variants cause a partial deletion or skipping of one or more exons in messenger RNAs. Among the 12 variants, nine lead to premature stop codons and predicted truncated ABCA4 proteins. At least two deep-intronic variants affect splice enhancer and silencer motifs and, therefore, these conserved sequences should be carefully evaluated when predicting the outcome of NCSS and deep-intronic variants.

A large intragenic deletion in the CLCN1 gene causes Hereditary Myotonia in pigs

Mutations in the CLCN1 gene are the primary cause of non-dystrophic Hereditary Myotonia in several animal species. However, there are no reports of Hereditary Myotonia in pigs to date. Therefore, the objective of the present study was to characterize the clinical and molecular findings of Hereditary Myotonia in an inbred pedigree. The clinical, electromyographic, histopathological, and molecular findings were evaluated. Clinically affected pigs presented non-dystrophic recessive Hereditary Myotonia. Nucleotide sequence analysis of the entire coding region of the CLCN1 gene revealed the absence of the exons 15 and 16 in myotonic animals. Analysis of the genomic region flanking the deletion unveiled a large intragenic deletion of 4,165 nucleotides. Interestingly, non-related, non-myotonic pigs expressed transcriptional levels of an alternate transcript (i.e., X2) that was identical to the deleted X1 transcript of myotonic pigs. All myotonic pigs and their progenitors were homozygous recessive and heterozygous, respectively, for the 4,165-nucleotide deletion. This is the first study reporting Hereditary Myotonia in pigs and characterizing its clinical and molecular findings. Moreover, to the best of our knowledge, Hereditary Myotonia has never been associated with a genomic deletion in the CLCN1 gene in any other species.

A Deep Intronic Variant in LDLR in Familial Hypercholesterolemia

BACKGROUND

Familial hypercholesterolemia (FH) is an inherited disorder characterized by high plasma LDL-C (low-density lipoprotein-cholesterol) levels. The vast majority of FH patients carry a mutation in the coding region of LDLR, APOB, or PCSK9. We set out to identify the culprit genetic defect in a large family with clinical FH, in whom no mutations were identified in the coding regions of these FH genes.

METHODS

Whole genome sequencing was performed in 5 affected and 4 unaffected individuals from a family with an unexplained autosomal dominant FH trait. The effect on splicing of the identified novel intronic LDLR mutation was ascertained by cDNA sequencing. The prevalence of the novel variant was assessed in 1 245 FH patients without an FH causing mutation identified by Sanger sequencing and in 2 154 patients referred for FH analysis by next-generation sequencing (covering the intronic region).

RESULTS

A novel deep intronic variant in LDLR (c.2140+103G>T) was found to cosegregate with high LDL-C in 5 patients, but was not present in 4 unaffected family members. The variant was shown to result in a 97 nucleotides insertion leading to a frameshift and premature stop codon in exon 15 of LDLR. The prevalence of the intronic variant was 0.24% (3/1245) in a cohort of FH patients without a known FH causing mutation and 0.23% (5/2154) in a population of FH patients referred for analysis by next-generation sequencing. Cosegregation analysis of a second family showed full penetrance of the novel variant with the FH phenotype over 3 generations.

CONCLUSIONS

The c.2140+103G>T mutation in LDLR is a novel intronic variant identified in FH that cosegregates with the FH phenotype. Our findings underline the need to analyze the intronic regions of LDLR in patients with FH, especially those in whom no mutation is found in the coding regions of LDLR, APOB, or PCSK9.

Deep-intronic variants in CNGB3 cause achromatopsia by pseudoexon activation

Our comprehensive cohort of 1100 unrelated achromatopsia (ACHM) patients comprises a considerable number of cases (~5%) harboring only a single pathogenic variant in the major ACHM gene CNGB3. We sequenced the entire CNGB3 locus in 33 of these patients to find a second variant which eventually explained the patients' phenotype. Forty-seven intronic CNGB3 variants were identified in 28 subjects after a filtering step based on frequency and the exclusion of variants found in cis with pathogenic alleles. In a second step, in silico prediction tools were used to filter out those variants with little odds of being deleterious. This left three variants that were analyzed using heterologous splicing assays. Variant c.1663-1205G>A, found in 14 subjects, and variant c.1663-2137C>T, found in two subjects, were indeed shown to exert a splicing defect by causing pseudoexon insertion into the transcript. Subsequent screening of further unsolved CNGB3 subjects identified four additional cases harboring the c.1663-1205G>A variant which makes it the eighth most frequent CNGB3 variant in our cohort. Compound heterozygosity could be validated in ten cases. Our study demonstrates that whole gene sequencing can be a powerful approach to identify the second pathogenic allele in patients apparently harboring only one disease-causing variant.

Pathogenic Abnormal Splicing Due to Intronic Deletions that Induce Biophysical Space Constraint for Spliceosome Assembly

A precise genetic diagnosis is the single most important step for families with genetic disorders to enable personalized and preventative medicine. In addition to genetic variants in coding regions (exons) that can change a protein sequence, abnormal pre-mRNA splicing can be devastating for the encoded protein, inducing a frameshift or in-frame deletion/insertion of multiple residues. Non-coding variants that disrupt splicing are extremely challenging to identify. Stemming from an initial clinical discovery in two index Australian families, we define 25 families with genetic disorders caused by a class of pathogenic non-coding splice variant due to intronic deletions. These pathogenic intronic deletions spare all consensus splice motifs, though they critically shorten the minimal distance between the 5' splice-site (5'SS) and branchpoint. The mechanistic basis for abnormal splicing is due to biophysical constraint precluding U1/U2 spliceosome assembly, which stalls in A-complexes (that bridge the 5'SS and branchpoint). Substitution of deleted nucleotides with non-specific sequences restores spliceosome assembly and normal splicing, arguing against loss of an intronic element as the primary causal basis. Incremental lengthening of 5'SS-branchpoint length in our index EMD case subject defines 45-47 nt as the critical elongation enabling (inefficient) spliceosome assembly for EMD intron 5. The 5'SS-branchpoint space constraint mechanism, not currently factored by genomic informatics pipelines, is relevant to diagnosis and precision medicine across the breadth of Mendelian disorders and cancer genomics.

COSSMO: predicting competitive alternative splice site selection using deep learning

Motivation

Alternative splice site selection is inherently competitive and the probability of a given splice site to be used also depends on the strength of neighboring sites. Here, we present a new model named the competitive splice site model (COSSMO), which explicitly accounts for these competitive effects and predicts the percent selected index (PSI) distribution over any number of putative splice sites. We model an alternative splicing event as the choice of a 3′ acceptor site conditional on a fixed upstream 5′ donor site or the choice of a 5′ donor site conditional on a fixed 3′ acceptor site. We build four different architectures that use convolutional layers, communication layers, long short-term memory and residual networks, respectively, to learn relevant motifs from sequence alone. We also construct a new dataset from genome annotations and RNA-Seq read data that we use to train our model.

Results

COSSMO is able to predict the most frequently used splice site with an accuracy of 70% on unseen test data, and achieve an R² of 0.6 in modeling the PSI distribution. We visualize the motifs that COSSMO learns from sequence and show that COSSMO recognizes the consensus splice site sequences and many known splicing factors with high specificity.

Availability and implementation

Model predictions, our training dataset, and code are available from http://cossmo.genes.toronto.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Association of lncRNA CCAT2 and CASC8 Gene Polymorphisms with Hepatocellular Carcinoma

The worldwide incidence of hepatocellular carcinoma (HCC), the major histological type of primary liver cancer, is heterogeneous due to the variable prevalence of etiological factors, indicating a correlation of HCC risk with genetic variations among individuals. Among long non-coding RNAs (lncRNAs) located in the chromosome 8q24 loci and involved in the carcinogenesis are colon cancer associated transcript 2 (CCAT2) and cancer susceptibility candidate 8 (CASC8). In this study, the association of CCAT2 and CASC8 gene polymorphisms with the occurrence of HCC was explored between 397 HCC patients and 1195 controls. We found that carriers of rs6983267 GG in CCAT2 were more susceptible to HCC, with the odds ratio (OR) and adjusted odds ratio (AOR) being 1.532 (95% CI, 1.103–2.129; p = 0.011) and 1.627 (95% CI, 1.120–2.265; p = 0.033), respectively. Moreover, for patients stratified by age (under 65), gender (male only), or status of drinking (habitual drinkers), a protective effect of CASC8 rs3843549 on presenting high Child–Pugh scores, metastatic vascular invasion, or large-size tumors was observed in a dominant model. Collectively, our data reveal association of CCAT2 and CASC8 gene polymorphisms with the occurrence and progression of HCC.

Analysis of MIR155HG variants and colorectal cancer susceptibility in Han Chinese population

BACKGROUND

MIR155HG plays an important role in malignant tumors, but it is rarely reported in the occurrence and development of colorectal cancer (CRC). This study investigated the effects of MIR155HG polymorphisms on CRC susceptibility from the perspective of molecular genetics.

METHODS

Eight SNPs in MIR155HG were selected and genotyped among 514 CRC cases and 510 healthy controls using the Agena MassARRAY platform. The associations between these SNPs and the CRC risk were evaluated under genetic models using conditional logistic regression analysis. The HaploReg v4.1 database was used for SNPs functional prediction.

RESULTS

The allele "C" of rs12482371 (p = 0.047), allele "C" of rs1893650 (p = 0.025), and the allele "A" of rs928883 (p = 0.037) in MIR155HG were significantly associated with CRC risk. Genetic model analysis revealed that rs12482371 and rs1893650 increased CRC risk; whereas rs928883 was associated with reduced CRC risk. Stratification analysis showed that rs9383938 was a protective factor in CRC patients under 60 years old. Rs12482371 and rs1893650 were associated with the CRC risk in females. Rs11911469 and rs34904192 may affect the clinical stage and lymph node metastasis. Moreover, the haplotypes CTT and GTC of LD block rs4143370|rs77218221|rs12482371, and the haplotypes CATGA and CACGG of LD block rs77699734|rs11911469|rs1893650|rs34904192|rs928883 were significantly associated with CRC risk.

CONCLUSION

This study revealed that MIR155HG SNPs were associated with CRC susceptibility and could be predictive biomarkers for CRC risk.

A deep intronic splice mutation of STAT3 underlies hyper IgE syndrome by negative dominance

Heterozygous in-frame mutations in coding regions of human STAT3 underlie the only known autosomal dominant form of hyper IgE syndrome (AD HIES). About 5% of familial cases remain unexplained. The mutant proteins are loss-of-function and dominant-negative when tested following overproduction in recipient cells. However, the production of mutant proteins has not been detected and quantified in the cells of heterozygous patients. We report a deep intronic heterozygous STAT3 mutation, c.1282-89C>T, in 7 relatives with AD HIES. This mutation creates a new exon in the STAT3 complementary DNA, which, when overexpressed, generates a mutant STAT3 protein (D427ins17) that is loss-of-function and dominant-negative in terms of tyrosine phosphorylation, DNA binding, and transcriptional activity. In immortalized B cells from these patients, the D427ins17 protein was 2 kDa larger and 4-fold less abundant than wild-type STAT3, on mass spectrometry. The patients' primary B and T lymphocytes responded poorly to STAT3-dependent cytokines. These findings are reminiscent of the impaired responses of leukocytes from other patients with AD HIES due to typical STAT3 coding mutations, providing further evidence for the dominance of the mutant intronic allele. These findings highlight the importance of sequencing STAT3 introns in patients with HIES without candidate variants in coding regions and essential splice sites. They also show that AD HIES-causing STAT3 mutant alleles can be dominant-negative even if the encoded protein is produced in significantly smaller amounts than wild-type STAT3.

Relation between DNA ionization potentials, single base substitutions and pathogenic variants

Background

It is nowadays clear that single base substitutions that occur in the human genome, of which some lead to pathogenic conditions, are non-random and influenced by their flanking nucleobase sequences. However, despite recent progress, the understanding of these "non-local" effects is still far from being achieved.

Results

To advance this problem, we analyzed the relationship between the base mutability in specific gene regions and the electron hole transport along the DNA base stacks, as it is one of the mechanisms that have been suggested to contribute to these effects. More precisely, we studied the connection between the normalized frequency of single base substitutions and the vertical ionization potential of the base and its flanking sequence, estimated using MP2/6-31G* ab initio quantum chemistry calculations. We found a statistically significant overall anticorrelation between these two quantities: the lower the vIP value, the more probable the substitution. Moreover, the slope of the regression lines varies. It is larger for introns than for exons and untranslated regions, and for synonymous than for missense substitutions. Interestingly, the correlation appears to be more pronounced when considering the flanking sequence of the substituted base in the 3’ rather than in the 5’ direction, which corresponds to the preferred direction of charge migration. A weaker but still statistically significant correlation is found between the ionization potentials and the pathogenicity of the base substitutions. Moreover, pathogenicity is also preferentially associated with larger changes in ionization potentials upon base substitution.

Conclusions

With this analysis we gained new insights into the complex biophysical mechanisms that are at the basis of mutagenesis and pathogenicity, and supported the role of electron-hole transport in these matters.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5867-y) contains supplementary material, which is available to authorized users.

Association of TLR4 and TLR9 polymorphisms and haplotypes with cervical cancer susceptibility

Single nucleotide polymorphisms (SNPs) in TLR genes may serve as a crucial marker for early susceptibility of various cancers including cervical cancer. The present study was therefore designed to ascertain the role of TLR4 and TLR9 SNPs and haplotypes to hrHPV infection and cervical cancer susceptibility. The study included 110 cervical cancer biopsies and 141 cervical smears from age-matched healthy controls of Gujarati ethnicity of Western India. hrHPV 16 and 18 were detected using Real-time PCR. Eight SNPs, four each in TLR4 and TLR9 were analyzed using Polymerase Chain Reaction-Restriction Fragment Length Polymorphism and Allele-Specific PCR. HPV 16 and 18 were detected in 68% cervical cancer cases. TLR4 rs4986790, rs1927911 and TLR9 rs187084 showed association with HPV 16/18 infection. CC and CT genotypes of TLR4 rs11536889 and rs1927911 respectively, and TC, CC genotypes of TLR9 rs187084, as well as minor alleles of TLR4 rs4986790 and TLR9 rs187084, were associated with the increased risk of cervical cancer. Stage-wise analysis revealed TLR9 rs187084 and rs352140 to be associated with early-stage cancer. TLR4 haplotype GTAC and TLR9 haplotype GATC were associated with the increased risk of cervical cancer while TLR4 haplotype GCAG was associated with the decreased risk. TLR4 haplotype GCAG and TLR9 haplotype GATC showed association with increased susceptibility to hrHPV infection. In conclusion, the present study revealed association of TLR4 and TLR9 polymorphisms and haplotypes with hrHPV infection and cervical cancer risk. Further evaluation of a larger sample size covering diverse ethnic populations globally is warranted.

In vivo measurements reveal a single 5'-intron is sufficient to increase protein expression level in Caenorhabditis elegans

Introns can increase gene expression levels using a variety of mechanisms collectively referred to as Intron Mediated Enhancement (IME). IME has been measured in cell culture and plant models by quantifying expression of intronless and intron-bearing reporter genes in vitro. We developed hardware and software to implement microfluidic chip-based gene expression quantification in vivo. We altered position, number and sequence of introns in reporter genes controlled by the hsp-90 promoter. Consistent with plant and mammalian studies, we determined a single, natural or synthetic, 5'-intron is sufficient for the full IME effect conferred by three synthetic introns, while a 3'-intron is not. We found coding sequence can affect IME; the same three synthetic introns that increase mcherry protein concentration by approximately 50%, increase mEGFP by 80%. We determined IME effect size is not greatly affected by the stronger vit-2 promoter. Our microfluidic imaging approach should facilitate screens for factors affecting IME and other intron-dependent processes.

Two Insertion/Deletion Variants within SPAG17 Gene Are Associated with Goat Body Measurement Traits

Simple Summary

Sperm-associated antigen 17 (SPAG17) is a reproduction and skeletal development related gene. This study aimed to identify crucial insertion-deletion (indel) variations, which influence the body measurement traits of goats. Two intronic indels (14 bp and 17 bp indels) were identified by sequencing. In Shaanbei white cashmere goat (SBWC), the different genotypes of the 14 bp indel were markedly associated with goat body height, chest width, body length, and chest depth. The genotypes of the 17 bp indel were significantly associated with body height and chest width. The different combined genotypes were significantly associated with body height and chest width of SBWC and ten traits of Hainan black goat. These results suggested that the 14 and 17 bp indels within SPAG17 can be used in goat growth related traits marker-assisted selection breeding, especially body height.

Abstract

Sperm-associated antigen 17 (SPAG17) gene encodes a multifunctional cytoplasmic protein, which influences not only reproduction but also skeletal development related body measurement traits, especially body height. Thus, this study aimed to identify crucial insertion-deletion (indel) variations, which influence the body measurement traits of goats in large goat populations (n = 1725). As a result, two intronic indels (14 bp and 17 bp indel) were identified by sequencing. For the two indel loci, the distributions of genotypes and alleles were significantly different between the Shaanbei white cashmere goat (SBWC) and the Hainan black goat (HNBG). In SBWC goats, the different genotypes of the 14 bp indel were markedly associated with goat body height, chest width, body length and chest depth. The genotypes of the 17 bp indel were significantly related to body height and chest width. At the two loci, for all seven analyzed traits of SBWC goat, the growth data of DD homozygotes were the worst, which means that the 14 bp insertion and the 17 bp deletion were beneficial and detrimental variations, respectively. Moreover, the combined genotypes were significantly related to body height and chest width of SBWC goats and ten traits of HNBG. These results suggested that the 14 and 17 bp indels within SPAG17 can be used in goat growth related traits marker-assisted selection breeding, especially body height.

Beyond the Big Five: Investigating Myostatin Structure, Polymorphism and Expression in Camelus dromedarius

Myostatin, a negative regulator of skeletal muscle mass in animals, has been shown to play a role in determining muscular hypertrophy in several livestock species, and a high degree of polymorphism has been previously reported for this gene in humans and cattle. In this study, we provide a characterization of the myostatin gene in the dromedary (Camelus dromedarius) at the genomic, transcript and protein level. The gene was found to share high structural and sequence similarity with other mammals, notably Old World camelids. 3D modeling highlighted several non-conservative SNP variants compared to the bovine, as well as putative functional variants involved in the stability of the myostatin dimer. NGS data for nine dromedaries from various countries revealed 66 novel SNPs, all of them falling either upstream or downstream the coding region. The analysis also confirmed the presence of three previously described SNPs in intron 1, predicted here to alter both splicing and transcription factor binding sites (TFBS), thus possibly impacting myostatin processing and/or regulation. Several putative TFBS were identified in the myostatin upstream region, some of them belonging to the myogenic regulatory factor family. Patterns of SNP distribution across countries, as suggested by Bayesian clustering of the nine dromedaries using the 69 SNPs, pointed to weak geographic differentiation, in line with known recurrent gene flow at ancient trading centers along caravan routes. Myostatin expression was investigated in a set of 8 skeletal muscles, both at transcript and protein level, via Digital Droplet PCR and Western Blotting, respectively. No significant differences were observed at the transcript level, while, at the protein level, the only significant differences concerned the promyostatin dimer (75 kDa), in four pair-wise comparisons, all involving the tensor fasciae latae muscle. Beside the mentioned band at 75 kDa, additional bands were observed at around 40 and 25 kDa, corresponding to the promyostatin monomer and the active C-terminal myostatin dimer, respectively. Their weaker intensity suggests that the unprocessed myostatin dimers could act as important reservoirs of slowly available myostatin forms. Under this assumption, the sequential cleavage steps may contribute additional layers of control within an already complex regulatory framework.

Functional mRNA analysis reveals aberrant splicing caused by novel intronic mutation in WDR45 in NBIA patient

WDR45 gene-associated neurodegeneration with brain iron accumulation (NBIA), referred to as beta-propeller protein-associated neurodegeneration (BPAN), is a rare disorder that presents with a very nonspecific clinical phenotype in children constituting global developmental delay. This case report illustrates the power of a combination of trio exome sequencing, in silico splicing analysis, and mRNA analysis to provide sufficient evidence for pathogenicity of a relatively intronic variant in WDR45, and in so doing, find a genetic diagnosis for a 6-year-old patient with developmental delay and seizures, a diagnosis which may otherwise have only been found once the characteristic MRI patterns of the disease became more obvious in young adulthood.

Risk stratification of cervical lesions using capture sequencing and machine learning method based on HPV and human integrated genomic profiles

From initial human papillomavirus (HPV) infection and precursor stages, the development of cervical cancer takes decades. High-sensitivity HPV DNA testing is currently recommended as primary screening method for cervical cancer, whereas better triage methodologies are encouraged to provide accurate risk management for HPV-positive women. Given that virus-driven genomic variation accumulates during cervical carcinogenesis, we designed a 39 Mb custom capture panel targeting 17 HPV types and 522 mutant genes related to cervical cancer. Using capture-based next-generation sequencing, HPV integration status, somatic mutation and copy number variation were analyzed on 34 paired samples, including 10 cases of HPV infection (HPV+), 10 cases of cervical intraepithelial neoplasia (CIN) grade and 14 cases of CIN2+ (CIN2: n = 1; CIN2-3: n = 3; CIN3: n = 9; squamous cell carcinoma: n = 1). Finally, the machine learning algorithm (Random Forest) was applied to build the risk stratification model for cervical precursor lesions based on CIN2+ enriched biomarkers. Generally, HPV integration events (11 in HPV+, 25 in CIN1 and 56 in CIN2+), non-synonymous mutations (2 in CIN1, 12 in CIN2+) and copy number variations (19.1 in HPV+, 29.4 in CIN1 and 127 in CIN2+) increased from HPV+ to CIN2+. Interestingly, ‘common’ deletion of mitochondrial chromosome was significantly observed in CIN2+ (P = 0.009). Together, CIN2+ enriched biomarkers, classified as HPV information, mutation, amplification, deletion and mitochondrial change, successfully predicted CIN2+ with average accuracy probability score of 0.814, and amplification and deletion ranked as the most important features. Our custom capture sequencing combined with machine learning method effectively stratified the risk of cervical lesions and provided valuable integrated triage strategies.

Clinical, pathophysiological and genetic features of motor symptoms in autosomal dominant Alzheimer's disease

Owing to an early and marked deposition of amyloid-β in the basal ganglia, autosomal dominant Alzheimer's disease could distinctly involve motor symptoms. Therefore, we aimed to assess the prevalence and characteristics of motor signs in autosomal dominant Alzheimer's disease. Baseline Unified Parkinson Disease Rating Scale part three scores (UPDRS-III) from 433 participants of the Dominantly Inherited Alzheimer's Network observational study were analysed. Motor symptoms were scrutinized with respect to associations with mutation carrier status, mutation site within PSEN1, basal ganglia amyloid-β as measured by Pittsburgh compound B PET, estimated years to symptom onset and Clinical Dementia Rating Scale-Sum of Boxes. Motor findings in mutation carriers were compared to patients with sporadic Alzheimer's disease using data of the National Alzheimer's Coordination Center. Mutation carriers showed motor findings at a higher frequency (28.4% versus 12.8%; P < 0.001) and severity (mean UPDRS-III scores 2.0 versus 0.4; P < 0.001) compared to non-carriers. Eleven of the 27 UPDRS-III items were statistically more frequently affected in mutation carriers after adjustment for multiple comparisons. Ten of these 11 items were subscale components of bradykinesia. In cognitively asymptomatic mutation carriers, dysdiadochokinesia was more frequent compared to non-carriers (right hand: 3.8% versus 0%; adjusted P = 0.023; left: 4.4% versus 0.6%; adjusted P = 0.031). In this cohort, the positive predictive value for mutation carrier status in cognitively asymptomatic participants (50% a priori risk) of dysdiadochokinesia was 100% for the right and 87.5% for the left side. Mutation carriers with motor findings more frequently were basal ganglia amyloid-β positive (84% versus 63.3%; P = 0.006) and showed more basal ganglia amyloid-β deposition (Pittsburgh compound B-standardized uptake value ratio 2.472 versus 1.928; P = 0.002) than those without. Frequency and severity of motor findings were greater in post-codon 200 PSEN1 mutations (36%; mean UPDRS-III score 3.03) compared to mutations pre-codon 200 PSEN1 (19.3%, P = 0.022; 0.91, P = 0.013). In mutation carriers, motor symptom severity was significantly positively correlated with basal ganglia amyloid-β deposition, Clinical Dementia Rating scores and estimated years to symptom onset. Mutation carriers with a Clinical Dementia Rating global score of 2 exhibited more pronounced motor symptoms than sporadic Alzheimer's disease patients with the same Clinical Dementia Rating global score (mean UPDRS-III scores 20.71 versus 5.96; P < 0.001). With a prevalence of approximately 30% and increasing severity with progression of dementia, motor symptoms are proven as a clinically relevant finding in autosomal dominant Alzheimer's disease, in particular in advanced dementia stages, that correlates with deposition of amyloid-β in the basal ganglia. In a very small per cent of cognitively asymptomatic members of families with autosomal dominant Alzheimer's disease, dysdiadochokinesia may increase the chance of an individual's status as mutation carrier.

Aberrant splicing contributes to severe α-spectrin-linked congenital hemolytic anemia

The etiology of severe hemolytic anemia in most patients with recessive hereditary spherocytosis (rHS) and the related disorder hereditary pyropoikilocytosis (HPP) is unknown. Whole exome sequencing of DNA from probands of 24 rHS or HPP kindreds identified numerous mutations in erythrocyte membrane α-spectrin (SPTA1). Twenty-eight mutations were novel, with null alleles frequently found in trans to missense mutations. No mutations were identified in a third of SPTA1 alleles (17/48). Whole genome sequencing revealed linkage disequilibrium between the common rHS-linked α-spectrinBug Hill polymorphism and a rare intron 30 variant in all 17 mutation-negative alleles. In vitro minigene studies and in vivo splicing analyses revealed the intron 30 variant changes a weak alternate branch point (BP) to a strong BP. This change leads to increased utilization of an alternate 3' splice acceptor site, perturbing normal α-spectrin mRNA splicing and creating an elongated mRNA transcript. In vivo mRNA stability studies revealed the newly created termination codon in the elongated transcript activates nonsense mediated decay leading to spectrin deficiency. These results demonstrate a unique mechanism of human genetic disease contributes to the etiology of a third of cases of rHS, facilitating diagnosis and treatment of severe anemia, and identifying a new target for therapeutic manipulation.

Case for genome sequencing in infants and children with rare, undiagnosed or genetic diseases

Up to 350 million people worldwide suffer from a rare disease, and while the individual diseases are rare, in aggregate they represent a substantial challenge to global health systems. The majority of rare disorders are genetic in origin, with children under the age of five disproportionately affected. As these conditions are difficult to identify clinically, genetic and genomic testing have become the backbone of diagnostic testing in this population. In the last 10 years, next-generation sequencing technologies have enabled testing of multiple disease genes simultaneously, ranging from targeted gene panels to exome sequencing (ES) and genome sequencing (GS). GS is quickly becoming a practical first-tier test, as cost decreases and performance improves. A growing number of studies demonstrate that GS can detect an unparalleled range of pathogenic abnormalities in a single laboratory workflow. GS has the potential to deliver unbiased, rapid and accurate molecular diagnoses to patients across diverse clinical indications and complex presentations. In this paper, we discuss clinical indications for testing and historical testing paradigms. Evidence supporting GS as a diagnostic tool is supported by superior genomic coverage, types of pathogenic variants detected, simpler laboratory workflow enabling shorter turnaround times, diagnostic and reanalysis yield, and impact on healthcare.

Correction of pseudoexon splicing caused by a novel intronic dysferlin mutation

OBJECTIVE

Dysferlin is a large transmembrane protein that functions in critical processes of membrane repair and vesicle fusion. Dysferlin-deficiency due to mutations in the dysferlin gene leads to muscular dystrophy (Miyoshi myopathy (MM), limb girdle muscular dystrophy type 2B (LGMD2B), distal myopathy with anterior tibial onset (DMAT)), typically with early adult onset. At least 416 pathogenic dysferlin mutations are known, but for approximately 17% of patients, one or both of their pathogenic variants remain undefined following standard exon sequencing methods that interrogate exons and nearby flanking intronic regions but not the majority of intronic regions.

METHODS

We sequenced RNA from myogenic cells to identify a novel dysferlin pathogenic variant in two affected siblings that previously had only one disease-causing variant identified. We designed antisense oligonucleotides (AONs) to bypass the effects of this mutation on RNA splicing.

RESULTS

We identified a new pathogenic point mutation deep within dysferlin intron 50i. This intronic variant causes aberrant mRNA splicing and inclusion of an additional pseudoexon (PE, we term PE50.1) within the mature dysferlin mRNA. PE50.1 inclusion alters the protein sequence, causing premature translation termination. We identified this mutation in 23 dysferlinopathy patients (seventeen families), revealing it to be one of the more prevalent dysferlin mutations. We used AON-mediated exon skipping to correct the aberrant PE50.1 splicing events in vitro, which increased normal mRNA production and significantly restored dysferlin protein expression.

INTERPRETATION

Deep intronic mutations can be a common underlying cause of dysferlinopathy, and importantly, could be treatable with AON-based exon-skipping strategies.

Genomic Analysis in the Age of Human Genome Sequencing

Affordable genome sequencing technologies promise to revolutionize the field of human genetics by enabling comprehensive studies that interrogate all classes of genome variation, genome-wide, across the entire allele frequency spectrum. Ongoing projects worldwide are sequencing many thousands-and soon millions-of human genomes as part of various gene mapping studies, biobanking efforts, and clinical programs. However, while genome sequencing data production has become routine, genome analysis and interpretation remain challenging endeavors with many limitations and caveats. Here, we review the current state of technologies for genetic variant discovery, genotyping, and functional interpretation and discuss the prospects for future advances. We focus on germline variants discovered by whole-genome sequencing, genome-wide functional genomic approaches for predicting and measuring variant functional effects, and implications for studies of common and rare human disease.

New Concepts About Familial Isolated Hyperparathyroidism

CONTEXT

Familial isolated hyperparathyroidism (FIHP) is defined as familial primary hyperparathyroidism (FH) without a characteristic extra-parathyroidal feature of a more complex hyperparathyroid syndrome. During 80 years, new concepts of FIHP have been developed within this definition. FIHP has been difficult to study due to small kindreds and mildly symptomatic cases.

EVIDENCE ACQUISITION

Searches were through PubMed for FIHP, other FH syndromes, and the gene(s) mutated in each.

EVIDENCE SYNTHESIS

Within its definition, the current concept of FIHP has clinical and mutational components. It can include incomplete expressions of MEN1 FHH, or HPT-JT or their mutations. Newest concepts of FIHP focus upon kindreds without mutation of either the MEN1, CASR, or CDC73 gene; 17% have germline activating mutation of the gene for the GCM2 transcription factor. Other genes for FIHP will probably be identified shortly. The FIHP kindreds with or without GCM2 mutation contain a median of only 2 cases of PHPT. The small kindred size in both subgroups of FIHP is probably caused low rate of screening among relatives. PHPT in FIHP with GCM2 mutation seems similar to PHPT in MEN1. Persons with FIHP and GCM2 mutation present as adults with mild hypercalcemia and multiple parathyroid tumors.

CONCLUSIONS

The current concept of FIHP led to a focus on small kindreds without mutation of MEN1, CASR, or CDC73. These assisted in the identification of germline activating GCM2 mutations in 17%. There is a need for clinical and mutational characterization in more cases to determine any unique clinical features of FIHP, either with or without mutation of GCM2.

Expanding the Boundaries of RNA Sequencing as a Diagnostic Tool for Rare Mendelian Disease

Gene-panel and whole-exome analyses are now standard methodologies for mutation detection in Mendelian disease. However, the diagnostic yield achieved is at best 50%, leaving the genetic basis for disease unsolved in many individuals. New approaches are thus needed to narrow the diagnostic gap. Whole-genome sequencing is one potential strategy, but it currently has variant-interpretation challenges, particularly for non-coding changes. In this study we focus on transcriptome analysis, specifically total RNA sequencing (RNA-seq), by using monogenetic neuromuscular disorders as proof of principle. We examined a cohort of 25 exome and/or panel "negative" cases and provided genetic resolution in 36% (9/25). Causative mutations were identified in coding and non-coding exons, as well as in intronic regions, and the mutational pathomechanisms included transcriptional repression, exon skipping, and intron inclusion. We address a key barrier of transcriptome-based diagnostics: the need for source material with disease-representative expression patterns. We establish that blood-based RNA-seq is not adequate for neuromuscular diagnostics, whereas myotubes generated by transdifferentiation from an individual's fibroblasts accurately reflect the muscle transcriptome and faithfully reveal disease-causing mutations. Our work confirms that RNA-seq can greatly improve diagnostic yield in genetically unresolved cases of Mendelian disease, defines strengths and challenges of the technology, and demonstrates the suitability of cell models for RNA-based diagnostics. Our data set the stage for development of RNA-seq as a powerful clinical diagnostic tool that can be applied to the large population of individuals with undiagnosed, rare diseases and provide a framework for establishing minimally invasive strategies for doing so.

Large-scale meta-analysis of mutations identified in panels of breast/ovarian cancer-related genes - Providing evidence of cancer predisposition genes

OBJECTIVE

Germline mutations occurring in the highly penetrant genes BRCA1 and BRCA2 are responsible for only certain cases of familial breast cancer (BC) and ovarian cancer (OC). Thus, the use of NGS multi-gene panel (MGP) testing has recently become very popular.

METHODS

To estimate a reliable BC and OC risk associated with pathogenic variants in the selected candidate BC/OC predisposition genes, a comprehensive meta-analysis of 48 MGP-based studies analyzing BC/OC patients was conducted. The role of 37 genes was evaluated, comparing, in total, the mutation frequency in ~120,000 BC/OC cases and ~120,000 controls, which guaranteed strong statistical support with high confidence for most analyzed genes.

RESULTS

We characterized the strategies of MGP analyses and the types and localizations of the identified mutations and showed that 13 and 11 of the analyzed genes were significantly associated with an increased BC and OC risk, respectively. The risk attributed to some of these genes (e.g., CDKN2A and PALB2 for BC) was similar to that observed for BRCA2. The analysis also showed a substantial difference in the profile of genes contributing to either BC or OC risk, including genes specifically associated with a high risk of OC but not BC (e.g., RAD51C, and RAD51D).

CONCLUSIONS

Our study provides strong statistical proof, defines the risk for many genes often considered candidates for BC/OC predisposition and excludes the role of other genes frequently analyzed in the MGPs. In the context of clinical diagnostics, the results support the knowledge-based interpretation of identified mutations.

Screening of BRCA1/2 deep intronic regions by targeted gene sequencing identifies the first germline BRCA1 variant causing pseudoexon activation in a patient with breast/ovarian cancer

BACKGROUND

Genetic analysis of BRCA1 and BRCA2 for the diagnosis of hereditary breast and ovarian cancer (HBOC) is commonly restricted to coding regions and exon-intron boundaries. Although germline pathogenic variants in these regions explain about ~20% of HBOC cases, there is still an important fraction that remains undiagnosed. We have screened BRCA1/2 deep intronic regions to identify potential spliceogenic variants that could explain part of the missing HBOC susceptibility.

METHODS

We analysed BRCA1/2 deep intronic regions by targeted gene sequencing in 192 high-risk HBOC families testing negative for BRCA1/2 during conventional analysis. Rare variants (MAF <0.005) predicted to create/activate splice sites were selected for further characterisation in patient RNA. The splicing outcome was analysed by RT-PCR and Sanger sequencing, and allelic imbalance was also determined when heterozygous exonic loci were present.

RESULTS

A novel transcript was detected in BRCA1 c.4185+4105C>T variant carrier. This variant promotes the inclusion of a pseudoexon in mature mRNA, generating an aberrant transcript predicted to encode for a non-functional protein. Quantitative and allele-specific assays determined haploinsufficiency in the variant carrier, supporting a pathogenic effect for this variant. Genotyping of 1030 HBOC cases and 327 controls did not identify additional carriers in Spanish population.

CONCLUSION

Screening of BRCA1/2 intronic regions has identified the first BRCA1 deep intronic variant associated with HBOC by pseudoexon activation. Although the frequency of deleterious variants in these regions appears to be low, our study highlights the importance of studying non-coding regions and performing comprehensive RNA assays to complement genetic diagnosis.

A rare deep intronic mutation of PKHD1 gene, c.8798-459 C > A, causes autosomal recessive polycystic kidney disease by pseudoexon activation

Autosomal recessive polycystic kidney disease (ARPKD), is a rare hepatorenal fibrocystic disorder primarily associated with progressive growth of multiple cysts in the kidneys causing progressive loss of renal function. The disease is linked to mutations in the PKHD1 gene. In this study, we describe the gene diagnosis and prenatal diagnosis for a consanguineous family with two fetuses diagnosed with polycystic kidney disease by fetal sonography during the pregnancy. Sequence analysis of cDNA synthesized from the PKHD1 mRNA of the second induced fetus identified a 111-nucleotide insert at the junction of exon 56 and 57 that originated from intervening sequence (IVS) 56. Further genomic sequencing of IVS 56 of the PKHD1 gene identified a rare homozygous deep intronic mutation (c.8798-459 C > A), which was inherited from the parents and not detectable in 100 unrelated control subjects. Moreover, we explored the pathogenicity of this deep intronic mutation by conducting a minigene splicing assay experiment, which demonstrated that the mutation causes a pseudoexon insertion, which results in a frameshift followed by a premature termination codon in exon 57. Eventually, the parents had a healthy baby by undergoing prenatal genetic diagnosis based on the targeted detection of the intron mutation. The newly identified deep intronic mutation is associated with a rare mechanism of abnormal splicing that expands the spectrum of known PKHD1 gene mutations. It can be used in evidence-based genetic and reproductive counseling for families with ARPKD.

Estimating contribution of rare non-coding variants to neuropsychiatric disorders

Owing to recent advances in DNA sequencing technology, a number of large-scale comprehensive analyses of genetic variations in protein-coding regions (i.e., whole-exome sequencing studies), have been conducted for neuropsychiatric and neurodevelopmental disorders, such as autism spectrum disorders, intellectual disability, and schizophrenia. These studies, especially those focusing on de novo (newly arising) mutations and extremely rare variants, have successfully identified previously unrecognized disease genes/mutations with a large effect size and deepen our understanding of the biology of neuropsychiatric diseases. Along with the continuously dropping sequencing cost, now the target of sequencing studies is expanding from the exome to the whole human genome. Several pioneering works have provided important insights into the contribution of rare non-coding variants to neuropsychiatric diseases. At the same time, these studies highlight need for further larger sample sizes and improvement in annotation of non-coding regulatory variants. In this review, key findings from recent studies as well as likely future directions are overviewed.

Intronic CNVs and gene expression variation in human populations

Introns can be extraordinarily large and they account for the majority of the DNA sequence in human genes. However, little is known about their population patterns of structural variation and their functional implication. By combining the most extensive maps of CNVs in human populations, we have found that intronic losses are the most frequent copy number variants (CNVs) in protein-coding genes in human, with 12,986 intronic deletions, affecting 4,147 genes (including 1,154 essential genes and 1,638 disease-related genes). This intronic length variation results in dozens of genes showing extreme population variability in size, with 40 genes with 10 or more different sizes and up to 150 allelic sizes. Intronic losses are frequent in evolutionarily ancient genes that are highly conserved at the protein sequence level. This result contrasts with losses overlapping exons, which are observed less often than expected by chance and almost exclusively affect primate-specific genes. An integrated analysis of CNVs and RNA-seq data showed that intronic loss can be associated with significant differences in gene expression levels in the population (CNV-eQTLs). These intronic CNV-eQTLs regions are enriched for intronic enhancers and can be associated with expression differences of other genes showing long distance intron-promoter 3D interactions. Our data suggests that intronic structural variation of protein-coding genes makes an important contribution to the variability of gene expression and splicing in human populations.

Cryptic intronic NBAS variant reveals the genetic basis of recurrent liver failure in a child

BACKGROUND

In almost half of patients with acute liver failure the cause is unknown, making targeted treatment and decisions about liver transplantation a challenge. Monogenic disorders may contribute to a significant proportion of these undiagnosed patients, and so the incorporation of technologies such as next generation sequencing (NGS) in the clinic could aid in providing a definitive diagnosis. However, this technology may present a major challenge in interpretation of sequence variants, particularly those in non-coding regions.

RESULTS

In this report we describe a case of Infantile liver failure syndrome 2 (ILFS2; MIM 616483) due to novel bi-allelic variants in the NBAS gene. A missense variant NM_015909.3(NBAS):c.2617C > T, NP_056993.2(NBAS):p.(Arg873Trp) was identified by whole genome sequencing (WGS). By combining WGS and reverse transcription-polymerase chain reaction (RT-PCR) we were able to identify a novel deep intronic variant, NM_015909.3(NBAS):c.2423 + 404G > C, leading to the inclusion of a pseudo-exon. This mechanism has not been described previously in this syndrome.

CONCLUSIONS

This study highlights the utility of analyzing NGS data in conjunction with investigating complementary DNA (cDNA) using techniques such as RT-PCR for detection of variants that otherwise would be likely to be missed in common NGS bioinformatic analysis pipelines. Combining these approaches, particularly when the phenotype match is strong, could lead to an increase in the diagnostic yield in acute liver failure and thus aid in targeted treatment, accurate genetic counseling and restoration of reproductive confidence.