Analysis of 200 unrelated individuals with a constitutional NF1 deep intronic pathogenic variant reveals that variants flanking the alternatively spliced NF1 exon 31 [23a] cause a classical neurofibromatosis type 1 phenotype while altering predominantly NF1 isoform type II

Neurofibromatosis type 1 results from loss-of-function NF1 pathogenic variants (PVs). Up to 30% of all NF1 PVs disrupt mRNA splicing, including deep intronic variants. Here, we retrospectively investigated the spectrum of NF1 deep intronic PVs in a cohort of 8,090 unrelated individuals from the University of Alabama at Birmingham (UAB) dataset with a molecularly confirmed neurofibromatosis type 1. All variants were identified through their effect on the NF1 transcript, followed by variant characterization at the DNA-level. A total of 68 distinct variants, which were ≥ 20 nucleotides away from the closest exon-intron junction, were identified in 2.5% unrelated individuals with NF1 (200/8,090). Nine different pathogenic splice variants, identified in 20 probands, led to exonization of different parts of intron 30 [23.2] or 31 [23a]. The two major NF1 transcript isoforms, distinguished by the absence (type I) or presence (type II) of the alternatively spliced cassette exon 31 [23a], are equally expressed in blood in control individuals without NF1 or NF1-affected individuals carrying their PV not in the introns flanking exon 31 [23a]. By fragment and cloning analysis we demonstrated that the exonization of intron 31 [23a] sequences due to deep intronic PV predominantly affects the NF1 isoform II. Seven additional (likely) pathogenic NF1 deep intronic variants not observed in the UAB dataset were found by classification of 36 variants identified by a literature search. Hence, the unique list of these 75 deep intronic (likely) PVs should be included in any comprehensive NF1 testing strategy.

Size matters: the impact of nucleus size on results from spatial transcriptomics

BACKGROUND

Visium Spatial Gene Expression (ST) is a method combining histological spatial information with transcriptomics profiles directly from tissue sections. The use of spatial information has made it possible to discover new modes of gene expression regulations. However, in the ST experiment, the nucleus size of cells may exceed the thickness of a tissue slice. This may, in turn, negatively affect comprehensive capturing the transcriptomics profile in a single slice, especially for tissues having large differences in the size of nuclei.

METHODS

Here, we defined the effect of Consecutive Slices Data Integration (CSDI) on unveiling accurate spot clustering and deconvolution of spatial transcriptomic spots in human postmortem brains. By considering the histological information as reference, we assessed the improvement of unsupervised clustering and single nuclei RNA-seq and ST data integration before and after CSDI.

RESULTS

Apart from the escalated number of defined clusters representing neuronal layers, the pattern of clusters in consecutive sections was concordant only after CSDI. Besides, the assigned cell labels to spots matches the histological pattern of tissue sections after CSDI.

CONCLUSION

CSDI can be applied to investigate consecutive sections studied with ST in the human cerebral cortex, avoiding misinterpretation of spot clustering and annotation, increasing accuracy of cell recognition as well as improvement in uncovering the layers of grey matter in the human brain.

Analysis of BRCA1 and BRCA2 alternative splicing in predisposition to ovarian cancer

BACKGROUND

The mRNA splicing is regulated on multiple levels, resulting in the proper distribution of genes' transcripts in each cell and maintaining cell homeostasis. At the same time, the expression of alternative transcripts can change in response to underlying genetic variants, often missed during routine diagnostics.

AIM

The main aim of this study was to define the frequency of aberrant splicing in BRCA1 and BRCA2 genes in blood RNA extracted from ovarian cancer patients who were previously found negative for the presence of pathogenic alterations in the 25 most commonly analysed ovarian cancer genes, including BRCA1 and BRCA2.

MATERIAL AND METHODS

Frequency and spectrum of splicing alterations in BRCA1 and BRCA2 genes were analysed in blood RNA from 101 ovarian cancer patients and healthy controls (80 healthy women) using PCR followed by gel electrophoresis and Sanger sequencing. The expression of splicing events was examined using RT-qPCR.

RESULTS

We did not identify any novel, potentially pathogenic splicing alterations. Nevertheless, we detected six naturally occurring transcripts, named BRCA1ΔE9-10, BRCA1ΔE11, BRCA1ΔE11q, and BRCA2ΔE3, BRCA2ΔE12 and BRCA2ΔE17-18 of which three (BRCA1ΔE11q, BRCA1ΔE11 and BRCA2ΔE3) were significantly higher expressed in the ovarian cancer cohort than in healthy controls (p ≤ 0.0001).

CONCLUSIONS

This observation indicates that the upregulation of selected naturally occurring transcripts can be stimulated by non-genetic mechanisms and be a potential systemic response to disease progression and/or treatment. However, this hypothesis requires further examination.

High prevalence of somatic PIK3CA and TP53 pathogenic variants in the normal mammary gland tissue of sporadic breast cancer patients revealed by duplex sequencing

The mammary gland undergoes hormonally stimulated cycles of proliferation, lactation, and involution. We hypothesized that these factors increase the mutational burden in glandular tissue and may explain high cancer incidence rate in the general population, and recurrent disease. Hence, we investigated the DNA sequence variants in the normal mammary gland, tumor, and peripheral blood from 52 reportedly sporadic breast cancer patients. Targeted resequencing of 542 cancer-associated genes revealed subclonal somatic pathogenic variants of: PIK3CA, TP53, AKT1, MAP3K1, CDH1, RB1, NCOR1, MED12, CBFB, TBX3, and TSHR in the normal mammary gland at considerable allelic frequencies (9 × 10-2- 5.2 × 10-1), indicating clonal expansion. Further evaluation of the frequently damaged PIK3CA and TP53 genes by ultra-sensitive duplex sequencing demonstrated a diversified picture of multiple low-level subclonal (in 10-2-10-4 alleles) hotspot pathogenic variants. Our results raise a question about the oncogenic potential in non-tumorous mammary gland tissue of breast-conserving surgery patients.

Targeted massively parallel sequencing of candidate regions on chromosome 22q predisposing to multiple schwannomas: An analysis of 51 individuals in a single-center experience

Constitutional LZTR1 or SMARCB1 pathogenic variants (PVs) have been found in ∼86% of familial and ∼40% of sporadic schwannomatosis cases. Hence, we performed massively parallel sequencing of the entire LZTR1, SMARCB1, and NF2 genomic loci in 35 individuals with schwannomas negative for constitutional first-hit PVs in the LZTR1/SMARCB1/NF2 coding sequences; however, with 22q deletion and/or a different NF2 PV in each tumor, including six cases with only one tumor available. Furthermore, we verified whether any other LZTR1/SMARCB1/NF2 (likely) PVs could be found in 16 cases carrying a SMARCB1 constitutional variant in the 3'-untranslated region (3'-UTR) c.*17C>T, c.*70C>T, or c.*82C>T. As no additional variants were found, functional studies were performed to clarify the effect of these 3'-UTR variants on the transcript. The 3'-UTR variants c.*17C>T and c.*82C>T showed pathogenicity by negatively affecting the SMARCB1 transcript level. Two novel deep intronic SMARCB1 variants, c.500+883T>G and c.500+887G>A, resulting in out-of-frame missplicing of intron 4, were identified in two unrelated individuals. Further resequencing of the entire repeat-masked genomics sequences of chromosome 22q in individuals negative for PVs in the SMARCB1/LZTR1/NF2 coding- and noncoding regions revealed five potential schwannomatosis-predisposing candidate genes, that is, MYO18B, NEFH, SGSM1, SGSM3, and SBF1, pending further verification.

Clinical spectrum of individuals with pathogenic NF1 missense variants affecting p.Met1149, p.Arg1276, and p.Lys1423: genotype-phenotype study in neurofibromatosis type 1

We report 281 individuals carrying a pathogenic recurrent NF1 missense variant at p.Met1149, p.Arg1276, or p.Lys1423, representing three nontruncating NF1 hotspots in the University of Alabama at Birmingham (UAB) cohort, together identified in 1.8% of unrelated NF1 individuals. About 25% (95% confidence interval: 20.5-31.2%) of individuals heterozygous for a pathogenic NF1 p.Met1149, p.Arg1276, or p.Lys1423 missense variant had a Noonan-like phenotype, which is significantly more compared with the "classic" NF1-affected cohorts (all p < .0001). Furthermore, p.Arg1276 and p.Lys1423 pathogenic missense variants were associated with a high prevalence of cardiovascular abnormalities, including pulmonic stenosis (all p < .0001), while p.Arg1276 variants had a high prevalence of symptomatic spinal neurofibromas (p < .0001) compared with "classic" NF1-affected cohorts. However, p.Met1149-positive individuals had a mild phenotype, characterized mainly by pigmentary manifestations without externally visible plexiform neurofibromas, symptomatic spinal neurofibromas or symptomatic optic pathway gliomas. As up to 0.4% of unrelated individuals in the UAB cohort carries a p.Met1149 missense variant, this finding will contribute to more accurate stratification of a significant number of NF1 individuals. Although clinically relevant genotype-phenotype correlations are rare in NF1, each affecting only a small percentage of individuals, together they impact counseling and management of a significant number of the NF1 population.

Bayesian multilevel model of micro RNA levels in ovarian-cancer and healthy subjects

In transcriptomics, micro RNAs (miRNAs) has gained much interest especially as potential disease indicators. However, apart from holding a great promise related to their clinical application, a lot of inconsistent results have been published. Our aim was to compare the miRNA expression levels in ovarian cancer and healthy subjects using the Bayesian multilevel model and to assess their potential usefulness in diagnosis. We have analyzed a case-control observational data on expression profiling of 49 preselected miRNA-based ovarian cancer indicators in 119 controls and 59 patients. A Bayesian multilevel model was used to characterize the effect of disease on miRNA levels controlling for differences in age and body weight. The difference between the miRNA level and health status of the patient on the scale of the data variability were discussed in the context of their potential usefulness in diagnosis. Additionally, the cross-validated area under the ROC curve (AUC) was used to assess the expected out-of-sample discrimination index of a different sets of miRNAs. The proposed model allowed us to describe the set of miRNA levels in patients and controls. Three highly correlated miRNAs: miR-101-3p, miR-142-5p, miR-148a-3p rank the highest with almost identical effect sizes that ranges from 0.45 to 1.0. For those miRNAs the credible interval for AUC ranged from 0.63 to 0.67 indicating their limited discrimination potential. A little benefit in adding information from other miRNAs was observed. There were several miRNAs in the dataset (miR-604, hsa-miR-221-5p) for which inferences were uncertain. For those miRNAs more experimental effort is needed to fully assess their effect in the context of new hits discovery and usefulness as disease indicators. The proposed multilevel Bayesian model can be used to characterize the panel of miRNA profile and to assess the difference in expression levels between healthy and cancer individuals.

Expanding the clinical phenotype of individuals with a 3-bp in-frame deletion of the NF1 gene (c.2970_2972del): an update of genotype-phenotype correlation

PURPOSE

Neurofibromatosis type 1 (NF1) is characterized by a highly variable clinical presentation, but almost all NF1-affected adults present with cutaneous and/or subcutaneous neurofibromas. Exceptions are individuals heterozygous for the NF1 in-frame deletion, c.2970_2972del (p.Met992del), associated with a mild phenotype without any externally visible tumors.

METHODS

A total of 135 individuals from 103 unrelated families, all carrying the constitutional NF1 p.Met992del pathogenic variant and clinically assessed using the same standardized phenotypic checklist form, were included in this study.

RESULTS

None of the individuals had externally visible plexiform or histopathologically confirmed cutaneous or subcutaneous neurofibromas. We did not identify any complications, such as symptomatic optic pathway gliomas (OPGs) or symptomatic spinal neurofibromas; however, 4.8% of individuals had nonoptic brain tumors, mostly low-grade and asymptomatic, and 38.8% had cognitive impairment/learning disabilities. In an individual with the NF1 constitutional c.2970_2972del and three astrocytomas, we provided proof that all were NF1-associated tumors given loss of heterozygosity at three intragenic NF1 microsatellite markers and c.2970_2972del.

CONCLUSION

We demonstrate that individuals with the NF1 p.Met992del pathogenic variant have a mild NF1 phenotype lacking clinically suspected plexiform, cutaneous, or subcutaneous neurofibromas. However, learning difficulties are clearly part of the phenotypic presentation in these individuals and will require specialized care.

Correction: Expanding the clinical phenotype of individuals with a 3-bp in-frame deletion of the NF1 gene (c.2970_2972del): an update of genotype-phenotype correlation

A correction has been published to this Article. The PDF and HTML have been updated accordingly.

Genotype-Phenotype Correlation in NF1: Evidence for a More Severe Phenotype Associated with Missense Mutations Affecting NF1 Codons 844-848

Neurofibromatosis type 1 (NF1), a common genetic disorder with a birth incidence of 1:2,000-3,000, is characterized by a highly variable clinical presentation. To date, only two clinically relevant intragenic genotype-phenotype correlations have been reported for NF1 missense mutations affecting p.Arg1809 and a single amino acid deletion p.Met922del. Both variants predispose to a distinct mild NF1 phenotype with neither externally visible cutaneous/plexiform neurofibromas nor other tumors. Here, we report 162 individuals (129 unrelated probands and 33 affected relatives) heterozygous for a constitutional missense mutation affecting one of five neighboring NF1 codons-Leu844, Cys845, Ala846, Leu847, and Gly848-located in the cysteine-serine-rich domain (CSRD). Collectively, these recurrent missense mutations affect ∼0.8% of unrelated NF1 mutation-positive probands in the University of Alabama at Birmingham (UAB) cohort. Major superficial plexiform neurofibromas and symptomatic spinal neurofibromas were more prevalent in these individuals compared with classic NF1-affected cohorts (both p < 0.0001). Nearly half of the individuals had symptomatic or asymptomatic optic pathway gliomas and/or skeletal abnormalities. Additionally, variants in this region seem to confer a high predisposition to develop malignancies compared with the general NF1-affected population (p = 0.0061). Our results demonstrate that these NF1 missense mutations, although located outside the GAP-related domain, may be an important risk factor for a severe presentation. A genotype-phenotype correlation at the NF1 region 844-848 exists and will be valuable in the management and genetic counseling of a significant number of individuals.

Application of high-resolution genomic profiling in the differential diagnosis of liposarcoma

Background

Rarity and heterogeneity of liposarcomas (LPS) make their diagnosis difficult even for sarcoma-experts pathologists. The molecular mechanism underlying the development and progression of liposarcomas (LPS) remains only partially known. In order to identify and compare the genomic profiles, we analyzed array-based comparative genomic hybridization (array-CGH) profiles of 66 liposarcomas, including well-differentiated (WDLPS), dedifferentiated (DDLPS) and myxoid (MLPS) subtypes.

Results

Copy number aberrations (CNAs) were identified in 98% of WDLPS and DDLPS and in 95% of MLPS cases. The minimal common region of amplification at 12q14.1q21.1 was observed in 96% of WDLPS and DDLPS cases. Four regions of CNAs, including losses of chromosome 6, 11 and 13 and gains of chromosome 14 were classified as recurrent in DDLPS; at least one was identified in 74% of DDLPS tumors. The DDLPS-associated losses were much more common in tumors with increased genomic complexity. In MLPS, the most frequent CNAs were losses of chromosome 6 (40%) and gains of chromosome 1 (30%), with the minimal overlapping regions 6q14.1q22.31 and 1q25.1q32.2, respectively.

Conclusions

Our findings show that the application of array-CGH allows to delineate clearly the genomic profiles of WDLPS, DDLPS and MLPS that reflect biological differences between these tumors. Although CNAs varied widely, the subtypes of tumors have characteristic genomic profiles that could facilitate the differential diagnosis of LPS subtypes, especially between WDLPS and DDLPS.

Electronic supplementary material

The online version of this article (doi:10.1186/s13039-017-0309-5) contains supplementary material, which is available to authorized users.

Genomic findings in patients with clinical suspicion of 22q11.2 deletion syndrome

Chromosome 22q11.2 deletion syndrome, one of the most common human genomic syndromes, has highly heterogeneous clinical presentation. Patients usually harbor a 1.5 to 3 Mb hemizygous deletion at chromosome 22q11.2, resulting in pathognomic TBX1, CRKL and/or MAPK1 haploinsufficiency. However, there are some individuals with clinical features resembling the syndrome who are eventually diagnosed with genomic disorders affecting other chromosomal regions. The objective of this study was to evaluate the additive value of high-resolution array-CGH testing in the cohort of 41 patients with clinical features of 22q11.2 deletion syndrome and negative results of standard cytogenetic diagnostic testing (karyotype and FISH for 22q11.2 locus). Array-CGH analysis revealed no aberrations at chromosomes 22 or 10 allegedly related to the syndrome. Five (12.2 %) patients were found to have other genomic imbalances, namely 17q21.31 microdeletion syndrome (MIM#610443), 1p36 deletion syndrome (MIM#607872), NF1 microduplication syndrome (MIM#613675), chromosome 6pter-p24 deletion syndrome (MIM#612582) and a novel interstitial deletion at 3q26.31 of 0.65 Mb encompassing a dosage-dependent gene NAALADL2. Our study demonstrates that the implementation of array-CGH into the panel of classic diagnostic procedures adds significantly to their efficacy. It allows for detection of constitutional genomic imbalances in 12 % of subjects with negative result of karyotype and FISH targeted for 22q11.2 region. Moreover, if used as first-tier genetic test, the method would provide immediate diagnosis in ∼40 % phenotypic 22q11.2 deletion subjects.

Cytokine profiling in exhaled breath condensate after exercise challenge in asthmatic children with post-exercise symptoms

INTRODUCTION

Markers of exhaled breath condensate (EBC) correlate with lung function impairment, airway remodeling and different aspects of the disease such as exercise-induced bronchoconstriction (EIB). Aim of the study was to determine the cytokine profile in EBC of children with asthma after an exercise treadmill challenge in order to obtain clinically useful information about mechanisms of EIB; also, to assess correlations between cytokine concentrations in EBC and clinical characteristics of the patients.

MATERIAL AND METHODS

The study population consisted of 25 randomly selected children, aged 8 to 19 years, with asthma and EIB symptoms despite the use of control medications. Patients on the day of the study visit underwent fractional exhaled nitric oxide measurement (FeNO) and baseline spirometry, performed an exercise treadmill challenge (ETC), and EBC samples were obtained at the end of the ETC.

RESULTS

In asthmatic children with positive ETC, monocyte hemotactic protein-1 (MCP-1) and IL-16 adjusted to pre-EBC forced expiratory volume in 1 s (FEV1) were significantly higher compared to children with negative ETC (p = 0.022 and p = 0.017 respectively). After adjustment to pre-EBC FEV1 other cytokines (IL-4, IL-5, IL-6, IL-8, MIG, TNF-α) were not related to post-exercise changes in FEV1.

CONCLUSIONS

We observed a specific inflammatory profile in the airways of asthmatic children with bronchoconstriction induced by exercise. The concentration of cytokines in EBC depended on the post-exercise decrease in FEV1, which was measured by the area under the curve (AUC). MCP-1 and IL-16, adjusted to pre-EBC FEV1, were significantly higher in children with a positive exercise challenge compared to those with a negative one.

Detection of somatic BRCA1/2 mutations in ovarian cancer - next-generation sequencing analysis of 100 cases

The overall prevalence of germline BRCA1/2 mutations is estimated between 11% and 15% of all ovarian cancers. Individuals with germline BRCA1/2 alterations treated with the PARP1 inhibitors (iPARP1) tend to respond better than patients with wild-type BRCA1/2. Additionally, also somatic BRCA1/2 alterations induce the sensitivity to iPARP1. Therefore, the detection of both germline and somatic BRCA1/2 mutations is required for effective iPARP1 treatment. The aim of this study was to identify the frequency and spectrum of germline and somatic BRCA1/2 alterations in a group of Polish patients with ovarian serous carcinoma. In total, 100 formalin-fixed paraffin-embedded (FFPE) ovarian serous carcinoma tissues were enrolled to the study. Mutational analysis of BRCA1/2 genes was performed by using next-generation sequencing. The presence of pathogenic variants was confirmed by Sanger sequencing. In addition, to confirm the germline or somatic status of the mutation, the nonneoplastic tissue was analyzed by bidirectional Sanger sequencing. In total, 27 (28% of patient samples) mutations (20 in BRCA1 and 7 in BRCA2) were identified. For 22 of 27 patients, nonneoplastic cells were available and sequencing revealed the somatic character of two BRCA1 (2/16; 12.5%) and two BRCA2 (2/6; 33%) mutations. Notably, we identified six novel frameshift or nonsense BRCA1/2 mutations. The heterogeneity of the detected mutations confirms the necessity of simultaneous analysis of BRCA1/2 genes in all patients diagnosed with serous ovarian carcinoma. Moreover, the use of tumor tissue for mutational analysis allowed the detection of both somatic and germline BRCA1/2 mutations.

Glomerulopathy in patients with distal duplication of chromosome 6p

Background

Duplication of the distal part of chromosome 6p is a rare genetic syndrome. Renal involvement has been reported in the majority of patients, including a wide range of congenital abnormalities of kidney and urinary tract and, occasionally, a proteinuric glomerulopathy.

Case presentation

Here, we report a 13-year-old girl with 6p25.3p22.1 duplication who presented with proteinuria in infancy, was later diagnosed as focal segmental glomerulosclerosis, progressed to end-stage renal disease and was successfully transplanted.

Conclusion

A systematic literature review suggests that 15–20 % of individuals with distal 6p duplication develop progressive proteinuric glomerulopathy. Monitoring of kidney function should be recommended in all cases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12882-016-0246-2) contains supplementary material, which is available to authorized users.

Chromosome 18q deletion syndrome with autoimmune diabetes mellitus: putative genomic loci for autoimmunity and immunodeficiency

A girl with 18q deletion syndrome was diagnosed with autoimmune diabetes mellitus and Hashimoto's thyroiditis at the age of 3 yr. In addition, the girl suffered from recurrent infections due to immunoglobulin A and IgG4 deficiency. She was also found to have CD3+CD4+FoxP3+, CD3+CD4+FoxP3+CD25+, and CD3+CD4+CD25+CD127 regulatory T cells deficiency. The exceptional coincidence of the two autoimmune disorders occurring at an early age, and associated with immune deficiency, implies that genes located on deleted 19.4 Mbp region at 18q21.32-q23 (chr18:58,660,699-78,012,870) might play a role in the pathogenesis of autoimmunity leading to β cell destruction and diabetes.

On the significance of germline cytogenetic rearrangements at MYCN locus in neuroblastoma

Background

MYCN oncogene amplification is the most important prognostic factor in neuroblastoma. 25% neuroblastoma tumors have somatic amplifications at this locus but little is known about its constitutional aberrations and their potential role in carcinogenesis. Here, we have performed an array-CGH and qPCR characterization of two patients with constitutional partial 2p trisomy including MYCN genomic region.

Results

One of the patients had congenital neuroblastoma and showed presence of minute areas of gains and losses within the common fragile site FRA2C at 2p24 encompassing MYCN. The link between 2p24 germline rearrangements and neuroblastoma development was reassessed by reviewing similar cases in the literature.

Conclusions

It appears that constitutional rearrangements involving chromosome 2p24 may play role in NB development.

Spectrum of NIPBL gene mutations in Polish patients with Cornelia de Lange syndrome

Cornelia de Lange syndrome (CdLS) is a rare multi-system genetic disorder characterised by growth and developmental delay, distinctive facial dysmorphism, limb malformations and multiple organ defects. The disease is caused by mutations in genes responsible for the formation and regulation of cohesin complex. About half of the cases result from mutations in the NIPBL gene coding delangin, a protein regulating the initialisation of cohesion. To date, approximately 250 point mutations have been identified in more than 300 CdLS patients worldwide. In the present study, conducted on a group of 64 unrelated Polish CdLS patients, 25 various NIPBL sequence variants, including 22 novel point mutations, were detected. Additionally, large genomic deletions on chromosome 5p13 encompassing the NIPBL gene locus were detected in two patients with the most severe CdLS phenotype. Taken together, 42 % of patients were found to have a deleterious alteration affecting the NIPBL gene, by and large private ones (89 %). The review of the types of mutations found so far in Polish patients, their frequency and correlation with the severity of the observed phenotype shows that Polish CdLS cases do not significantly differ from other populations.