Histone H3K36 mutations promote sarcomagenesis through altered histone methylation landscape

Several types of pediatric cancers reportedly contain high-frequency missense mutations in histone H3, yet the underlying oncogenic mechanism remains poorly characterized. Here we report that the H3 lysine 36-to-methionine (H3K36M) mutation impairs the differentiation of mesenchymal progenitor cells and generates undifferentiated sarcoma in vivo. H3K36M mutant nucleosomes inhibit the enzymatic activities of several H3K36 methyltransferases. Depleting H3K36 methyltransferases, or expressing an H3K36I mutant that similarly inhibits H3K36 methylation, is sufficient to phenocopy the H3K36M mutation. After the loss of H3K36 methylation, a genome-wide gain in H3K27 methylation leads to a redistribution of polycomb repressive complex 1 and de-repression of its target genes known to block mesenchymal differentiation. Our findings are mirrored in human undifferentiated sarcomas in which novel K36M/I mutations in H3.1 are identified.

Single-cell RNA-seq reveals that glioblastoma recapitulates a normal neurodevelopmental hierarchy

Cancer stem cells are critical for cancer initiation, development, and treatment resistance. Our understanding of these processes, and how they relate to glioblastoma heterogeneity, is limited. To overcome these limitations, we performed single-cell RNA sequencing on 53586 adult glioblastoma cells and 22637 normal human fetal brain cells, and compared the lineage hierarchy of the developing human brain to the transcriptome of cancer cells. We find a conserved neural tri-lineage cancer hierarchy centered around glial progenitor-like cells. We also find that this progenitor population contains the majority of the cancer's cycling cells, and, using RNA velocity, is often the originator of the other cell types. Finally, we show that this hierarchal map can be used to identify therapeutic targets specific to progenitor cancer stem cells. Our analyses show that normal brain development reconciles glioblastoma development, suggests a possible origin for glioblastoma hierarchy, and helps to identify cancer stem cell-specific targets.

Germline and somatic FGFR1 abnormalities in dysembryoplastic neuroepithelial tumors

Dysembryoplastic neuroepithelial tumor (DNET) is a benign brain tumor associated with intractable drug-resistant epilepsy. In order to identify underlying genetic alterations and molecular mechanisms, we examined three family members affected by multinodular DNETs as well as 100 sporadic tumors from 96 patients, which had been referred to us as DNETs. We performed whole-exome sequencing on 46 tumors and targeted sequencing for hotspot FGFR1 mutations and BRAF p.V600E was used on the remaining samples. FISH, copy number variation assays and Sanger sequencing were used to validate the findings. By whole-exome sequencing of the familial cases, we identified a novel germline FGFR1 mutation, p.R661P. Somatic activating FGFR1 mutations (p.N546K or p.K656E) were observed in the tumor samples and further evidence for functional relevance was obtained by in silico modeling. The FGFR1 p.K656E mutation was confirmed to be in cis with the germline p.R661P variant. In 43 sporadic cases, in which the diagnosis of DNET could be confirmed on central blinded neuropathology review, FGFR1 alterations were also frequent and mainly comprised intragenic tyrosine kinase FGFR1 duplication and multiple mutants in cis (25/43; 58.1 %) while BRAF p.V600E alterations were absent (0/43). In contrast, in 53 cases, in which the diagnosis of DNET was not confirmed, FGFR1 alterations were less common (10/53; 19 %; p < 0.0001) and hotspot BRAF p.V600E (12/53; 22.6 %) (p < 0.001) prevailed. We observed overexpression of phospho-ERK in FGFR1 p.R661P and p.N546K mutant expressing HEK293 cells as well as FGFR1 mutated tumor samples, supporting enhanced MAP kinase pathway activation under these conditions. In conclusion, constitutional and somatic FGFR1 alterations and MAP kinase pathway activation are key events in the pathogenesis of DNET. These findings point the way towards existing targeted therapies.

Recessive osteogenesis imperfecta caused by missense mutations in SPARC

Secreted protein, acidic, cysteine-rich (SPARC) is a glycoprotein that binds to collagen type I and other proteins in the extracellular matrix. Using whole-exome sequencing to identify the molecular defect in two unrelated girls with severe bone fragility and a clinical diagnosis of osteogenesis imperfecta type IV, we identified two homozygous variants in SPARC (GenBank: NM_003118.3; c.497G>A [p.Arg166His] in individual 1; c.787G>A [p.Glu263Lys] in individual 2). Published modeling and site-directed mutagenesis studies had previously shown that the residues substituted by these mutations form an intramolecular salt bridge in SPARC and are essential for the binding of SPARC to collagen type I. The amount of SPARC secreted by skin fibroblasts was reduced in individual 1 but appeared normal in individual 2. The migration of collagen type I alpha chains produced by these fibroblasts was mildly delayed on SDS-PAGE gel, suggesting some overmodification of collagen during triple helical formation. Pulse-chase experiments showed that collagen type I secretion was mildly delayed in skin fibroblasts from both individuals. Analysis of an iliac bone sample from individual 2 showed that trabecular bone was hypermineralized on the material level. In conclusion, these observations show that homozygous mutations in SPARC can give rise to severe bone fragility in humans.

QTL for several metabolic traits map to loci controlling growth and body composition in an F2 intercross between high- and low-growth chicken lines

Quantitative trait loci (QTL) for metabolic and body composition traits were mapped at 7 and 9 wk, respectively, in an F2 intercross between high-growth and low-growth chicken lines. These lines also diverged for abdominal fat percentage (AFP) and plasma insulin-like growth factor-I (IGF-I), insulin, and glucose levels. Genotypings were performed with 129 microsatellite markers covering 21 chromosomes. A total of 21 QTL with genomewide level of significance were detected by single-trait analyses for body weight (BW), breast muscle weight (BMW) and percentage (BMP), AF weight (AFW) and percentage (AFP), shank length (ShL) and diameter (ShD), fasting plasma glucose level (Gluc), and body temperature (Tb). Other suggestive QTL were identified for these parameters and for plasma IGF-I and nonesterified fatty acid levels. QTL controlling adiposity and Gluc were colocalized on GGA3 and GGA5 and QTL for BW, ShL and ShD, adiposity, and Tb on GGA4. Multitrait analyses revealed two QTL controlling Gluc and AFP on GGA5 and Gluc and Tb on GGA26. Significant effects of the reciprocal cross were observed on BW, ShD, BMW, and Gluc, which may result from mtDNA and/or maternal effects. Most QTL regions for Gluc and adiposity harbor genes for which alleles have been associated with increased susceptibility to diabetes and/or obesity in humans. Identification of genes responsible for these metabolic QTL will increase our understanding of the constitutive “hyperglycemia” found in chickens. Furthermore, a comparative approach could provide new information on the genetic causes of diabetes and obesity in humans.

APRDX1 mutant allele causes a MMACHC secondary epimutation in cblC patients

To date, epimutations reported in man have been somatic and erased in germlines. Here, we identify a cause of the autosomal recessive cblC class of inborn errors of vitamin B₁₂ metabolism that we name “epi-cblC”. The subjects are compound heterozygotes for a genetic mutation and for a promoter epimutation, detected in blood, fibroblasts, and sperm, at the MMACHC locus; 5-azacytidine restores the expression of MMACHC in fibroblasts. MMACHC is flanked by CCDC163P and PRDX1, which are in the opposite orientation. The epimutation is present in three generations and results from PRDX1 mutations that force antisense transcription of MMACHC thereby possibly generating a H3K36me3 mark. The silencing of PRDX1 transcription leads to partial hypomethylation of the epiallele and restores the expression of MMACHC. This example of epi-cblC demonstrates the need to search for compound epigenetic-genetic heterozygosity in patients with typical disease manifestation and genetic heterozygosity in disease-causing genes located in other gene trios.

Inborn errors of vitamin B₁₂ metabolism of the cblC class are caused by mutations in the MMACHC gene. Here, Guéant et al. report epi-cblC, a class of cblC in which patients are compound heterozygous for a genetic mutation and a secondary epimutation at the MMACHC locus.

Whole-exome sequencing as a diagnostic tool: current challenges and future opportunities

Whole-exome sequencing (WES) represents a significant breakthrough in the field of human genetics. This technology has largely contributed to the identification of new disease-causing genes and is now entering clinical laboratories. WES represents a powerful tool for diagnosis and could reduce the 'diagnostic odyssey' for many patients. In this review, we present a technical overview of WES analysis, variants annotation and interpretation in a clinical setting. We evaluate the usefulness of clinical WES in different clinical indications, such as rare diseases, cancer and complex diseases. Finally, we discuss the efficacy of WES as a diagnostic tool and the impact on patient management.

DGCR8 microprocessor defect characterizes familial multinodular goiter with schwannomatosis

BACKGROUNDDICER1 is the only miRNA biogenesis component associated with an inherited tumor syndrome, featuring multinodular goiter (MNG) and rare pediatric-onset lesions. Other susceptibility genes for familial forms of MNG likely exist.METHODSWhole-exome sequencing of a kindred with early-onset MNG and schwannomatosis was followed by investigation of germline pathogenic variants that fully segregated with the disease. Genome-wide analyses were performed on 13 tissue samples from familial and nonfamilial DGCR8-E518K-positive tumors, including MNG, schwannomas, papillary thyroid cancers (PTCs), and Wilms tumors. miRNA profiles of 4 tissue types were compared, and sequencing of miRNA, pre-miRNA, and mRNA was performed in a subset of 9 schwannomas, 4 of which harbor DGCR8-E518K.RESULTSWe identified c.1552G>A;p.E518K in DGCR8, a microprocessor component located in 22q, in the kindred. The variant identified is a somatic hotspot in Wilms tumors and has been identified in 2 PTCs. Copy number loss of chromosome 22q, leading to loss of heterozygosity at the DGCR8 locus, was found in all 13 samples harboring c.1552G>A;p.E518K. miRNA profiling of PTCs, MNG, schwannomas, and Wilms tumors revealed a common profile among E518K hemizygous tumors. In vitro cleavage demonstrated improper processing of pre-miRNA by DGCR8-E518K. MicroRNA and RNA profiling show that this variant disrupts precursor microRNA production, impacting populations of canonical microRNAs and mirtrons.CONCLUSIONWe identified DGCR8 as the cause of an unreported autosomal dominant mendelian tumor susceptibility syndrome: familial multinodular goiter with schwannomatosis.FUNDINGCanadian Institutes of Health Research, Compute Canada, Alex's Lemonade Stand Foundation, the Mia Neri Foundation for Childhood Cancer, Cassa di Sovvenzioni e Risparmio fra il Personale della Banca d'Italia, and the KinderKrebsInitiative Buchholz/Holm-Seppensen.

Molecular analyses reveal close similarities between small cell carcinoma of the ovary, hypercalcemic type and atypical teratoid/rhabdoid tumor

Small cell carcinoma of the ovary, hypercalcemic type (SCCOHT) is the most common undifferentiated ovarian malignancy diagnosed in women under age 40. We and others recently determined that germline and/or somatic deleterious mutations in SMARCA4 characterize SCCOHT. Alterations in this gene, or the related SWI/SNF chromatin remodeling gene SMARCB1, have been previously reported in atypical teratoid/rhabdoid tumors (ATRTs) and malignant rhabdoid tumors (MRTs). To further describe the somatic landscape of SCCOHT, we performed whole exome sequencing on 14 tumors and their matched normal tissues and compared their genomic alterations with those in ATRT and ovarian high grade serous carcinoma (HGSC). We confirmed that SMARCA4 is the only recurrently mutated gene in SCCOHT, and show that recurrent allelic imbalance is observed exclusively on chromosome 19p, where SMARCA4 resides. By comparing genomic alterations between SCCOHT, ATRT and HGSC, we demonstrate that SCCOHTs, like ATRTs, have a remarkably simple genome and harbor significantly fewer somatic protein-coding mutations and chromosomal alterations than HGSC. Furthermore, a comparison of global DNA methylation profiles of 45 SCCOHTs, 65 ATRTs, and 92 HGSCs demonstrates a strong epigenetic correlation between SCCOHT and ATRT. Our results further confirm that the genomic and epigenomic signatures of SCCOHT are more similar to those of ATRT than HGSC, supporting our previous hypothesis that SCCOHT is a rhabdoid tumor and should be renamed MRT of the ovary. Furthermore, we conclude that SMARCA4 inactivation is the main cause of SCCOHT, and that new distinct therapeutic approaches should be developed to specifically target this devastating tumor.

Loss-of-Function Mutation in APC2 Causes Sotos Syndrome Features

Sotos syndrome, characterized by intellectual disability and characteristic facial features, is caused by haploinsufficiency in the NSD1 gene. We conducted an etiological study on two siblings with Sotos features without mutations in NSD1 and detected a homozygous frameshift mutation in the APC2 gene by whole-exome sequencing, which resulted in the loss of function of cytoskeletal regulation in neurons. Apc2-deficient (Apc2^-/-) mice exhibited impaired learning and memory abilities along with an abnormal head shape. Endogenous Apc2 expression was downregulated by the knockdown of Nsd1, indicating that APC2 is a downstream effector of NSD1 in neurons. Nsd1 knockdown in embryonic mouse brains impaired the migration and laminar positioning of cortical neurons, as observed in Apc2^-/- mice, and this defect was rescued by the forced expression of Apc2. Thus, APC2 is a crucial target of NSD1, which provides an explanation for the intellectual disability associated with Sotos syndrome.

Functionally Null RAD51D Missense Mutation Associates Strongly with Ovarian Carcinoma

RAD51D is a key player in DNA repair by homologous recombination (HR), and RAD51D truncating variant carriers have an increased risk for ovarian cancer. However, the contribution of nontruncating RAD51D variants to cancer predisposition remains uncertain. Using deep sequencing and case-control genotyping studies, we show that in French Canadians, the missense RAD51D variant c.620C>T;p.S207L is highly prevalent and is associated with a significantly increased risk for ovarian high-grade serous carcinoma (HGSC; 3.8% cases vs. 0.2% controls). The frequency of the p.S207L variant did not significantly differ from that of controls in breast, endometrial, pancreas, or colorectal adenocarcinomas. Functionally, we show that this mutation impairs HR by disrupting the RAD51D-XRCC2 interaction and confers PARP inhibitor sensitivity. These results highlight the importance of a functional RAD51D-XRCC2 interaction to promote HR and prevent the development of HGSC. This study identifies c.620C>T;p.S207L as the first bona fide pathogenic RAD51D missense cancer susceptibility allele and supports the use of targeted PARP-inhibitor therapies in ovarian cancer patients carrying deleterious missense RAD51D variants. Cancer Res; 77(16); 4517-29. ©2017 AACR.

Invasive growth associated with Cold-Inducible RNA-Binding Protein expression drives recurrence of surgically resected brain metastases

BACKGROUND

Sixty percent of surgically resected brain metastases (BrM) recur within 1 year. These recurrences have long been thought to result from the dispersion of cancer cells during surgery. We tested the alternative hypothesis that invasion of cancer cells into the adjacent brain plays a significant role in local recurrence and shortened overall survival.

METHODS

We determined the invasion pattern of 164 surgically resected BrM and correlated with local recurrence and overall survival. We performed single-cell RNA sequencing (scRNAseq) of >15,000 cells from BrM and adjacent brain tissue. Validation of targets was performed with a novel cohort of BrM patient-derived xenografts (PDX) and patient tissues.

RESULTS

We demonstrate that invasion of metastatic cancer cells into the adjacent brain is associated with local recurrence and shortened overall survival. scRNAseq of paired tumor and adjacent brain samples confirmed the existence of invasive cancer cells in the tumor-adjacent brain. Analysis of these cells identified Cold-Inducible RNA-Binding Protein (CIRBP) overexpression in invasive cancer cells compared to cancer cells located within the metastases. Applying PDX models that recapitulate the invasion pattern observed in patients, we show that CIRBP is overexpressed in highly invasive BrM and is required for efficient invasive growth in the brain.

CONCLUSIONS

These data demonstrate peritumoral invasion as a driver of treatment failure in BrM that is functionally mediated by CIRBP. These findings improve our understanding of the biology underlying post-operative treatment failure and lay the groundwork for rational clinical trial development based upon invasion pattern in surgically resected brain metastases.

Transcriptome profiling of granulosa and theca cells during dominant follicle development in the horse

Several aspects of equine ovarian physiology are unique among domestic species. Moreover, follicular growth patterns are very similar between horses and humans. This study aimed to characterize, for the first time, global gene expression profiles associated with growth and preovulatory (PO) maturation of equine dominant follicles. Granulosa cells (GCs) and theca interna cells (TCs) were harvested from follicles (n = 5) at different stages of an ovulatory wave in mares corresponding to early dominance (ED; diameter ≥22 mm), late dominance (LD; ≥33 mm) and PO stage (34 h after administration of crude equine gonadotropins at LD stage), and separately analyzed on a horse gene expression microarray, followed by validation using quantitative PCR and immunoblotting/immunohistochemistry. Numbers of differentially expressed transcripts (DETs; ≥2-fold; P < 0.05) during the ED-LD and LD-PO transitions were 546 and 2419 in GCs and 5 and 582 in TCs. The most prominent change in GCs was the down-regulation of transcripts associated with cell division during both ED-LD and LD-PO. In addition, DET sets during LD-PO in GCs were enriched for genes involved in cell communication/adhesion, antioxidation/detoxification, immunity/inflammation, and cholesterol biosynthesis. In contrast, the largest change in TCs during the LD-PO transition was an up-regulation of genes involved in immune activation, with other DET sets mapping to GPCR/cAMP signaling, lipid/amino acid metabolism, and cell proliferation/survival and differentiation. In conclusion, distinct expression profiles were identified between growing and PO follicles and, particularly, between GCs and TCs within each stage. Several DETs were identified that have not been associated with follicle development in other species.

Glioblastoma scRNAseq Shows Treatment-induced, Immune-dependent Rise In Mesenchymal Cancer Cells, and Structural Variants in Distal Neural Stem Cells

BACKGROUND

Glioblastoma is a treatment-resistant brain cancer. Its hierarchical cellular nature and its tumour microenvironment (TME) before, during, and after treatments remain unresolved.

METHODS

Here, we used single-cell RNA-sequencing to analyze new and recurrent glioblastoma, and the nearby subventricular zone (SVZ).

RESULTS

We found four glioblastoma neural lineages are present in new and recurrent glioblastoma with an enrichment of the cancer mesenchymal lineage, immune cells, and reactive astrocytes in early recurrences. Cancer lineages were hierarchically organized around cycling oligodendrocytic and astrocytic progenitors that are transcriptomically similar but distinct to SVZ neural stem cells (NSCs). Furthermore, NSCs from the SVZ of glioblastoma patients harbored glioblastoma chromosomal anomalies. Lastly, mesenchymal cancer cells and TME reactive astrocytes shared similar gene signatures which were induced by radiotherapy in a myeloid-dependent fashion in vivo.

CONCLUSION

These data reveal the dynamic, immune-dependent nature of glioblastoma's response to treatments and identify distant NSCs as likely cells of origin.

Molecular characterization of DICER1-mutated pituitary blastoma

Pituitary blastoma (PitB) has recently been identified as a rare and potentially lethal pediatric intracranial tumor. All cases that have been studied molecularly possess at least one DICER1 pathogenic variant. Here, we characterized nine pituitary samples, including three fresh frozen PitBs, three normal fetal pituitary glands and three normal postnatal pituitary glands using small-RNA-Seq, RNA-Seq, methylation profiling, whole genome sequencing and Nanostring® miRNA analyses; an extended series of 21 pituitary samples was used for validation purposes. These analyses demonstrated that DICER1 RNase IIIb hotspot mutations in PitBs induced improper processing of miRNA precursors, resulting in aberrant 5p-derived miRNA products and a skewed distribution of miRNAs favoring mature 3p over 5p miRNAs. This led to dysregulation of hundreds of 5p and 3p miRNAs and concomitant dysregulation of numerous mRNA targets. Gene expression analysis revealed PRAME as the most significantly upregulated gene (500-fold increase). PRAME is a member of the Retinoic Acid Receptor (RAR) signaling pathway and in PitBs, the RAR, WNT and NOTCH pathways are dysregulated. Cancer Hallmarks analysis showed that PI3K pathway is activated in the tumors. Whole genome sequencing demonstrated a quiet genome with very few somatic alterations. The comparison of methylation profiles to publicly available data from ~ 3000 other central nervous system tumors revealed that PitBs have a distinct methylation profile compared to all other tumors, including pituitary adenomas. In conclusion, this comprehensive characterization of DICER1-related PitB revealed key molecular underpinnings of PitB and identified pathways that could potentially be exploited in the treatment of this tumor.

NHLRC2 variants identified in patients with fibrosis, neurodegeneration, and cerebral angiomatosis (FINCA): characterisation of a novel cerebropulmonary disease

A novel multi-organ disease that is fatal in early childhood was identified in three patients from two non-consanguineous families. These children were born asymptomatic but at the age of 2 months they manifested progressive multi-organ symptoms resembling no previously known disease. The main clinical features included progressive cerebropulmonary symptoms, malabsorption, progressive growth failure, recurrent infections, chronic haemolytic anaemia and transient liver dysfunction. In the affected children, neuropathology revealed increased angiomatosis-like leptomeningeal, cortical and superficial white matter vascularisation and congestion, vacuolar degeneration and myelin loss in white matter, as well as neuronal degeneration. Interstitial fibrosis and previously undescribed granuloma-like lesions were observed in the lungs. Hepatomegaly, steatosis and collagen accumulation were detected in the liver. A whole-exome sequencing of the two unrelated families with the affected children revealed the transmission of two heterozygous variants in the NHL repeat-containing protein 2 (NHLRC2); an amino acid substitution p.Asp148Tyr and a frameshift 2-bp deletion p.Arg201GlyfsTer6. NHLRC2 is highly conserved and expressed in multiple organs and its function is unknown. It contains a thioredoxin-like domain; however, an insulin turbidity assay on human recombinant NHLRC2 showed no thioredoxin activity. In patient-derived fibroblasts, NHLRC2 levels were low, and only p.Asp148Tyr was expressed. Therefore, the allele with the frameshift deletion is likely non-functional. Development of the Nhlrc2 null mouse strain stalled before the morula stage. Morpholino knockdown of nhlrc2 in zebrafish embryos affected the integrity of cells in the midbrain region. This is the first description of a fatal, early-onset disease; we have named it FINCA disease based on the combination of pathological features that include fibrosis, neurodegeneration, and cerebral angiomatosis.

ExomeAI: detection of recurrent allelic imbalance in tumors using whole-exome sequencing data

SUMMARY

Whole-exome sequencing (WES) has extensively been used in cancer genome studies; however, the use of WES data in the study of loss of heterozygosity or more generally allelic imbalance (AI) has so far been very limited, which highlights the need for user-friendly and flexible software that can handle low-quality datasets. We have developed a statistical approach, ExomeAI, for the detection of recurrent AI events using WES datasets, specifically where matched normal samples are not available.

AVAILABILITY

ExomeAI is a web-based application, publicly available at: http://genomequebec.mcgill.ca/exomeai.

CONTACT

JavadNadaf@gmail.com or somayyeh.fahiminiya@mcgill.ca

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Effect of the prior distribution of SNP effects on the estimation of total breeding value

Background

Five main methods, commonly applied in genomic selection, were used to estimate the GEBV on the 15^th QTLMAS workshop dataset: GBLUP, LASSO, Bayes A and two Bayes B type of methods (BBn and BBt). GBLUP is a mixed model approach where GEBV are obtained using a relationship matrix calculated from the SNP genotypes. The remaining methods are regression-based approaches where the SNP effects are first estimated and, then GEBV are calculated given the individuals' genotypes.

Methods

The differences between the regression-based methods are in their prior distributions for the SNP effects. The prior distribution for LASSO is a Laplace distribution, for Bayes A is a scaled Student-t distribution, and the Bayes B type methods have a Spike and Slab prior where only a proportion (π) of SNP has an effect, following a given distribution. In this study, two different distributions were considered for the Bayes B type methods: (i) normal and (ii) scaled Student-t. They are referred here as the BBn and BBt methods, respectively. These prior distributions are defined by one or more parameters controlling their scale/rate (λ), shape (df) or proportion of SNP with effect (π). LASSO requires one (λ); two for Bayes A (λ, df) and Bayes Bn (λ, π); and three for Bayes Bt (λ, df, π). In this study, all parameters were estimated from the data. An extra scenario for Bayes A and BBt was included where df was not estimated but fixed to 4 (suffixed _4df). The implementation of GBLUP was done using ASREML, the heritability was also estimated from the data. All other methods were implemented using a MCMC approach.

Results

All Bayes A and B methods showed accuracy (correlation between True and Estimated BV) as high as 0.94 except for BA_4df (r = 0.91). Compared to the traditional BLUP using pedigree information, these methods improved the accuracy between 50 and 55%. GBLUP and LASSO were less accurate (0.81 and 0.85 respectively) and the improvements were 34 and 40% compared to BLUP.

Conclusions

Results of all methods were consistent and the accuracies for GEBV ranged between 0.81 and 0.94. When all parameters were estimated the results were similar for the Bayes A and Bayes B methods. Results showed that Bayes A was more sensitive to the changes in the shape parameter, and the parameter changes led to change in the accuracy of GEBV. However BBt was more robust to the change in this parameter. This may be explained by the fact that BBt estimates one extra parameter and it can buffer against a non-proper shape parameter.

A functionally impaired missense variant identified in French Canadian families implicates FANCI as a candidate ovarian cancer-predisposing gene

Background

Familial ovarian cancer (OC) cases not harbouring pathogenic variants in either of the BRCA1 and BRCA2 OC-predisposing genes, which function in homologous recombination (HR) of DNA, could involve pathogenic variants in other DNA repair pathway genes.

Methods

Whole exome sequencing was used to identify rare variants in HR genes in a BRCA1 and BRCA2 pathogenic variant negative OC family of French Canadian (FC) ancestry, a population exhibiting genetic drift. OC cases and cancer-free individuals from FC and non-FC populations were investigated for carrier frequency of FANCI c.1813C>T; p.L605F, the top-ranking candidate. Gene and protein expression were investigated in cancer cell lines and tissue microarrays, respectively.

Results

In FC subjects, c.1813C>T was more common in familial (7.1%, 3/42) than sporadic (1.6%, 7/439) OC cases (P = 0.048). Carriers were detected in 2.5% (74/2950) of cancer-free females though female/male carriers were more likely to have a first-degree relative with OC (121/5249, 2.3%; Spearman correlation = 0.037; P = 0.011), suggesting a role in risk. Many of the cancer-free females had host factors known to reduce risk to OC which could influence cancer risk in this population. There was an increased carrier frequency of FANCI c.1813C>T in BRCA1 and BRCA2 pathogenic variant negative OC families, when including the discovery family, compared to cancer-free females (3/23, 13%; OR = 5.8; 95%CI = 1.7–19; P = 0.005). In non-FC subjects, 10 candidate FANCI variants were identified in 4.1% (21/516) of Australian OC cases negative for pathogenic variants in BRCA1 and BRCA2, including 10 carriers of FANCI c.1813C>T. Candidate variants were significantly more common in familial OC than in sporadic OC (P = 0.04). Localization of FANCD2, part of the FANCI-FANCD2 (ID2) binding complex in the Fanconi anaemia (FA) pathway, to sites of induced DNA damage was severely impeded in cells expressing the p.L605F isoform. This isoform was expressed at a reduced level, destabilized by DNA damaging agent treatment in both HeLa and OC cell lines, and exhibited sensitivity to cisplatin but not to a poly (ADP-ribose) polymerase inhibitor. By tissue microarray analyses, FANCI protein was consistently expressed in fallopian tube epithelial cells and only expressed at low-to-moderate levels in 88% (83/94) of OC samples.

Conclusions

This is the first study to describe candidate OC variants in FANCI, a member of the ID2 complex of the FA DNA repair pathway. Our data suggest that pathogenic FANCI variants may modify OC risk in cancer families.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00998-5.

Publisher Correction: A PRDX1 mutant allele causes a MMACHC secondary epimutation in cblC patients

The original version of this Article contained an error in the title, which was incorrectly given as 'APRDX1 mutant allele causes a MMACHC secondary epimutation in cblC patients'. This has now been corrected in both the PDF and HTML versions of the Article to read 'A PRDX1 mutant allele causes a MMACHC secondary epimutation in cblC patients'.

The Genetic Causes of Nonsyndromic Congenital Retinal Detachment: A Genetic and Phenotypic Study of Pakistani Families

Purpose

To evaluate consanguineous pedigrees from Pakistan with a clinical diagnosis of nonsyndromic congenital retinal nonattachment (NCRNA) and identify genes responsible for the disease as currently only one NCRNA gene is known (atonal basic helix-loop-helix transcription factor 7: ATOH7).

Methods

We implemented a three-step genotyping platform: single nucleotide polymorphism genotyping to identify loss of heterozygosity regions in patients, Retinal Information Network panel screening for mutations in currently known retinal genes. Negative patients were then subjected to whole exome sequencing.

Results

We evaluated 21 consanguineous NCRNA pedigrees and identified the causal mutations in known retinal genes in 13 out of our 21 families. We found mutations in ATOH7 in three families. Surprisingly, we then found mutations in familial exudative vitreoretinopathy (FEVR) genes; low-density lipoprotein receptor-related protein 5 mutations (six families), tetraspanin 12 mutations (two families), and NDP mutations (two families). Thus, 62% of the patients were successfully genotyped in our study with seven novel and six previously reported mutations in known retinal genes.

Conclusions

Although the clinical diagnosis of all children was NCRNA with severe congenital fibrotic retinal detachments, the molecular diagnosis determined that the disease process was in fact a very severe form of FEVR in 10 families. Because severe congenital retinal detachment has not been previously associated with all the FEVR genes, we have thus expanded the phenotypic spectrum of FEVR, a highly variable retinal detachment phenotype that has clinical overlap with NCRNA. We identified seven novel mutations. We also established for the first time genetic overlap between the Iranian and Pakistani populations. We identified eight NCRNA families that do not harbor mutations in any known retinal genes, suggesting novel causal genes in these families.

PRAME protein expression in DICER1-related tumours

DICER1 syndrome is an autosomal dominant tumour predisposition syndrome usually affecting persons under 30 years of age. Many of the associated benign and malignant lesions occur almost exclusively in DICER1 syndrome. One such tumour, pituitary blastoma (pitB), overexpresses PRAME 500x above control levels. PRAME (PReferentially expressed Antigen in MElanoma) is expressed in malignancies that are not DICER1-related (e.g. melanoma). To address whether PRAME expression is part of the DICER1 phenotype, or simply a feature of pitB, a series of 75 DICER1-mutated specimens and 33 non-mutated specimens was surveyed using immunohistochemistry for PRAME, together with EZH2, which complexes with PRAME. In DICER1-mutated specimens, positive staining for PRAME was only seen in malignant tumours; 7 of 11 histological types and 34/62 individual tumours were positive, while non-tumourous lesions were always negative. Pleuropulmonary blastoma (PPB) showed a continuum in staining, with type I lesions being PRAME negative (n = 7) but all type II and type III lesions PRAME positive (n = 7). Similarly, cystic nephroma (CN) was negative (n = 8), with anaplastic sarcoma of the kidney being positive (n = 2). However, one atypical CN with mesenchymal cell proliferation was PRAME-positive. Embryonal rhabdomyosarcoma (RMS) with DICER1 pathogenic variants (PVs) was positive for PRAME (5/6), but the same tumour type without DICER1 PVs was also positive (9/15). Staining for EZH2 corresponded to that seen with PRAME, validating the latter. This study leads us to conclude that (1) PRAME expression occurs in two-thirds of DICER1-related malignancies; (2) PRAME may be a marker for the progression that certain DICER1-related lesions are thought to undergo, such as PPB and CN; and (3) PRAME expression in some tumours, such as RMS, appears to be an intrinsic feature of the tumour, rather than specifically related to DICER1 PVs. Therapy directed against PRAME may offer novel treatment options in patients with the DICER1 syndrome.

Exome sequencing study of partial agenesis of the corpus callosum in men with developmental delay, epilepsy, and microcephaly

Background

This study reports the genetic features of four Caucasian males from the Saguenay‒Lac‐St‐Jean region affected by partial agenesis of the corpus callosum (ACC) with hypotonia, epilepsy, developmental delay, microcephaly, hypoplasia, and autistic behavior.

Methods

We performed whole exome sequencing (WES) to identify new genes involved in this pathological phenotype. The regions of interest were subsequently sequenced for family members.

Results

Single‐nucleotide variations (SNVs) and insertions or deletions were detected in genes potentially implicated in brain defects observed in these patients. One patient did not have mutations in genes related to ACC, but carried a de novo pathogenic mutation in Mucolipin‐1 (MCOLN1) and was diagnosed with mucolipidosis type IV. Among the other probands, missense SNVs were observed in DCLK2 (Doublecortin Like Kinase 2), HERC2 (HECT And RLD Domain Containing E3 Ubiquitin Protein Ligase 2), and KCNH3 (Potassium channel, voltage‐gated, subfamily H, member 3). One patient also carried a non‐frameshift insertion in CACNA1A (Cav2.1(P/Q‐type) calcium channels).

Conclusion

Although no common genetic defect was observed in this study, we provide evidence for new avenues of investigation for ACC, such as molecular pathways involving HERC2, CACNA1A, KCNH3, and more importantly DCLK2. We also allowed to diagnose an individual with mucolipidosis type IV.