B-cell immune deficiency in twin sisters expands the phenotype of MOPDI.

Microcephalic osteodysplastic primordial dwarfism type I (MOPDI) is a very rare and severe autosomal recessive disorder characterized by marked intrauterine growth retardation, skeletal dysplasia, microcephaly and brain malformations. MOPDI is caused by biallelic mutations in RNU4ATAC, a non-coding gene involved in U12-type splicing of 1% of the introns in the genome, which are recognized by their specific splicing consensus sequences. Here, we describe a unique observation of immunodeficiency in twin sisters with mild MOPDI, who harbor a novel n.108_126del mutation, encompassing part of the U4atac snRNA 3' stem-loop and Sm protein binding site, and the previously reported n.111G>A mutation. Interestingly, both twin sisters show mild B-cell anomalies, including low naive B-cell counts and increased memory B-cell and plasmablasts counts, suggesting partial and transitory blockage of B-cell maturation and/or excessive activation of naive B-cells. Hence, the localization of a mutation in stem II of U4atac snRNA, as observed in another RNU4ATAC-opathy with immunodeficiency, that is, Roifman syndrome (RFMN), is not required for the occurrence of an immune deficiency. Finally, we emphasize the importance of considering immunodeficiency in MOPDI management to reduce the risk of serious infectious episodes.

PHGDH-related microcephalic dwarfism in two fetuses: Expanding the phenotypical spectrum of L-serine biosynthesis defect.

Defects in L-serine biosynthesis are a group of autosomal recessive diseases resulting in a wide phenotypic spectrum ranging from viable to lethal presentations and caused by variants in the three genes encoding the L-serine biosynthesis enzymes, PHGDH, PSAT1, and PSPH. Neu-Laxova syndrome (NLS) is the fetal form of this group, characterized by multiple congenital anomalies including severe intrauterine growth retardation, cutaneous lesions extending from ichthyosis to severe restrictive dermopathy with ectropion and eclabion, edema, microcephaly, central nervous system abnormalities, and flexion contractures. Here we report on two unrelated fetuses with an attenuated phenotype of NLS, that initially evoked Taybi-Linder syndrome. They carry biallelic pathogenic variants in the PHGDH gene. These observations expand the phenotypic continuum of L-serine biosynthesis defects, and illustrate the phenotypic overlap between NLS and microcephalic primordial dwarfism.

Deficiency of the minor spliceosome component U4atac snRNA secondarily results in ciliary defects in human and zebrafish.

In the human genome, about 750 genes contain one intron excised by the minor spliceosome. This spliceosome comprises its own set of snRNAs, among which U4atac. Its noncoding gene, RNU4ATAC, has been found mutated in Taybi-Linder (TALS/microcephalic osteodysplastic primordial dwarfism type 1), Roifman (RFMN), and Lowry-Wood (LWS) syndromes. These rare developmental disorders, whose physiopathological mechanisms remain unsolved, associate ante- and post-natal growth retardation, microcephaly, skeletal dysplasia, intellectual disability, retinal dystrophy, and immunodeficiency. Here, we report bi-allelic RNU4ATAC mutations in five patients presenting with traits suggestive of the Joubert syndrome (JBTS), a well-characterized ciliopathy. These patients also present with traits typical of TALS/RFMN/LWS, thus widening the clinical spectrum of RNU4ATAC-associated disorders and indicating ciliary dysfunction as a mechanism downstream of minor splicing defects. Intriguingly, all five patients carry the n.16G>A mutation, in the Stem II domain, either at the homozygous or compound heterozygous state. A gene ontology term enrichment analysis on minor intron-containing genes reveals that the cilium assembly process is over-represented, with no less than 86 cilium-related genes containing at least one minor intron, among which there are 23 ciliopathy-related genes. The link between RNU4ATAC mutations and ciliopathy traits is supported by alterations of primary cilium function in TALS and JBTS-like patient fibroblasts, as well as by u4atac zebrafish model, which exhibits ciliopathy-related phenotypes and ciliary defects. These phenotypes could be rescued by WT but not by pathogenic variants-carrying human U4atac. Altogether, our data indicate that alteration of cilium biogenesis is part of the physiopathological mechanisms of TALS/RFMN/LWS, secondarily to defects of minor intron splicing.

Mutations in the non-coding RNU4ATAC gene affect the homeostasis and function of the Integrator complex.

Various genetic diseases associated with microcephaly and developmental defects are due to pathogenic variants in the U4atac small nuclear RNA (snRNA), a component of the minor spliceosome essential for the removal of U12-type introns from eukaryotic mRNAs. While it has been shown that a few RNU4ATAC mutations result in impaired binding of essential protein components, the molecular defects of the vast majority of variants are still unknown. Here, we used lymphoblastoid cells derived from RNU4ATAC compound heterozygous (g.108_126del;g.111G>A) twin patients with MOPD1 phenotypes to analyze the molecular consequences of the mutations on small nuclear ribonucleoproteins (snRNPs) formation and on splicing. We found that the U4atac108_126del mutant is unstable and that the U4atac111G>A mutant as well as the minor di- and tri-snRNPs are present at reduced levels. Our results also reveal the existence of 3'-extended snRNA transcripts in patients' cells. Moreover, we show that the mutant cells have alterations in splicing of INTS7 and INTS10 minor introns, contain lower levels of the INTS7 and INTS10 proteins and display changes in the assembly of Integrator subunits. Altogether, our results show that compound heterozygous g.108_126del;g.111G>A mutations induce splicing defects and affect the homeostasis and function of the Integrator complex.

Diagnostic chest X-rays and breast cancer risk among women with a hereditary predisposition to breast cancer unexplained by a BRCA1 or BRCA2 mutation.

BACKGROUND: Diagnostic ionizing radiation is a risk factor for breast cancer (BC). BC risk increases with increased dose to the chest and decreases with increased age at exposure, with possible effect modification related to familial or genetic predisposition. While chest X-rays increase the BC risk of BRCA1/2 mutation carriers compared to non-carriers, little is known for women with a hereditary predisposition to BC but who tested negative for a BRCA1 or BRCA2 (BRCA1/2) mutation. METHODS: We evaluated the effect of chest X-rays from diagnostic medical procedures in a dataset composed of 1552 BC cases identified through French family cancer clinics and 1363 unrelated controls. Participants reported their history of X-ray exposures in a detailed questionnaire and were tested for 113 DNA repair genes. Logistic regression and multinomial logistic regression models were used to assess the association with BC. RESULTS: Chest X-ray exposure doubled BC risk. A 3% increased BC risk per additional exposure was observed. Being 20 years old or younger at first exposure or being exposed before first full-term pregnancy did not seem to modify this risk. Birth after 1960 or carrying a rare likely deleterious coding variant in a DNA repair gene other than BRCA1/2 modified the effect of chest X-ray exposure. CONCLUSION: Ever/never chest X-ray exposure increases BC risk 2-fold regardless of age at first exposure and, by up to 5-fold when carrying 3 or more rare variants in a DNA repair gene. Further studies are needed to evaluate other DNA repair genes or variants to identify those which could modify radiation sensitivity. Identification of subpopulations that are more or less susceptible to ionizing radiation is important and potentially clinically relevant.

Altered regulation of BRCA1 exon 11 splicing is associated with breast cancer risk in carriers of BRCA1 pathogenic variants.

Germline pathogenic variants in BRCA1 confer a high risk of developing breast and ovarian cancer. The BRCA1 exon 11 (formally exon 10) is one of the largest exons and codes for the nuclear localization signals of the corresponding gene product. This exon can be partially or entirely skipped during pre-mRNA splicing, leading to three major in-frame isoforms that are detectable in most cell types and tissue, and in normal and cancer settings. However, it is unclear whether the splicing imbalance of this exon is associated with cancer risk. Here we identify a common genetic variant in intron 10, rs5820483 (NC_000017.11:g.43095106_43095108dup), which is associated with exon 11 isoform expression and alternative splicing, and with the risk of breast cancer, but not ovarian cancer, in BRCA1 pathogenic variant carriers. The identification of this genetic effect was confirmed by analogous observations in mouse cells and tissue in which a loxP sequence was inserted in the syntenic intronic region. The prediction that the rs5820483 minor allele variant would create a binding site for the splicing silencer hnRNP A1 was confirmed by pull-down assays. Our data suggest that perturbation of BRCA1 exon 11 splicing modifies the breast cancer risk conferred by pathogenic variants of this gene.

5' Region Large Genomic Rearrangements in the BRCA1 Gene in French Families: Identification of a Tandem Triplication and Nine Distinct Deletions with Five Recurrent Breakpoints.

BACKGROUND: Large genomic rearrangements (LGR) in BRCA1 consisting of deletions/duplications of one or several exons have been found throughout the gene with a large proportion occurring in the 5' region from the promoter to exon 2. The aim of this study was to better characterize those LGR in French high-risk breast/ovarian cancer families. METHODS: DNA from 20 families with one apparent duplication and nine deletions was analyzed with a dedicated comparative genomic hybridization (CGH) array, high-resolution BRCA1 Genomic Morse Codes analysis and Sanger sequencing. RESULTS: The apparent duplication was in fact a tandem triplication of exons 1 and 2 and part of intron 2 of BRCA1, fully characterized here for the first time. We calculated a causality score with the multifactorial model from data obtained from six families, classifying this variant as benign. Among the nine deletions detected in this region, eight have never been identified. The breakpoints fell in six recurrent regions and could confirm some specific conformation of the chromatin. CONCLUSIONS: Taken together, our results firmly establish that the BRCA1 5' region is a frequent site of different LGRs and highlight the importance of the segmental duplication and Alu sequences, particularly the very high homologous region, in the mechanism of a recombination event. This also confirmed that those events are not systematically deleterious.

Gene- and pathway-level analyses of iCOGS variants highlight novel signaling pathways underlying familial breast cancer susceptibility.

Single-nucleotide polymorphisms (SNPs) in over 180 loci have been associated with breast cancer (BC) through genome-wide association studies involving mostly unselected population-based case-control series. Some of them modify BC risk of women carrying a BRCA1 or BRCA2 (BRCA1/2) mutation and may also explain BC risk variability in BC-prone families with no BRCA1/2 mutation. Here, we assessed the contribution of SNPs of the iCOGS array in GENESIS consisting of BC cases with no BRCA1/2 mutation and a sister with BC, and population controls. Genotyping data were available for 1281 index cases, 731 sisters with BC, 457 unaffected sisters and 1272 controls. In addition to the standard SNP-level analysis using index cases and controls, we performed pedigree-based association tests to capture transmission information in the sibships. We also performed gene- and pathway-level analyses to maximize the power to detect associations with lower-frequency SNPs or those with modest effect sizes. While SNP-level analyses identified 18 loci, gene-level analyses identified 112 genes. Furthermore, 31 Kyoto Encyclopedia of Genes and Genomes and 7 Atlas of Cancer Signaling Network pathways were highlighted (false discovery rate of 5%). Using results from the "index case-control" analysis, we built pathway-derived polygenic risk scores (PRS) and assessed their performance in the population-based CECILE study and in a data set composed of GENESIS-affected sisters and CECILE controls. Although these PRS had poor predictive value in the general population, they performed better than a PRS built using our SNP-level findings, and we found that the joint effect of family history and PRS needs to be considered in risk prediction models.

A de novo frameshift pathogenic variant in TBR1 identified in autism without intellectual disability.

BACKGROUND: In order to be able to provide accurate genetic counseling to patients with Autism Spectrum Disorder (ASD), it is crucial to identify correlations between heterogeneous phenotypes and genetic alterations. Among the hundreds of de novo pathogenic variants reported in ASD, single-nucleotide variations and small insertions/deletions were reported in TBR1. This gene encodes a transcription factor that plays a key role in brain development. Pathogenic variants in TBR1 are often associated with severe forms of ASD, including intellectual disability and language impairment. METHODS: Adults diagnosed with ASD but without intellectual disability (diagnosis of Asperger syndrome, according to the DSM-IV) took part in a genetic consultation encompassing metabolic assessments, a molecular karyotype and the screening of a panel of 268 genes involved in intellectual disability, ASD and epilepsy. In addition, the patient reported here went through a neuropsychological assessment, structural magnetic resonance imaging and magnetic resonance spectroscopy measurements. RESULTS: Here, we report the case of a young adult male who presents with a typical form of ASD. Importantly, this patient presents with no intellectual disability or language impairment, despite a de novo heterozygous frameshift pathogenic variant in TBR1, leading to an early premature termination codon (c.26del, p.(Pro9Leufs*12)). CONCLUSION: Based on this case report, we discuss the role of TBR1 in general brain development, language development, intellectual disability and other symptoms of ASD. Providing a detailed clinical description of the individuals with such pathogenic variants should help to understand the genotype-phenotype relationships in ASD.

Clinical interpretation of variants identified in RNU4ATAC, a non-coding spliceosomal gene.

Biallelic variants in RNU4ATAC, a non-coding gene transcribed into the minor spliceosome component U4atac snRNA, are responsible for three rare recessive developmental diseases, namely Taybi-Linder/MOPD1, Roifman and Lowry-Wood syndromes. Next-generation sequencing of clinically heterogeneous cohorts (children with either a suspected genetic disorder or a congenital microcephaly) recently identified mutations in this gene, illustrating how profoundly these technologies are modifying genetic testing and assessment. As RNU4ATAC has a single non-coding exon, the bioinformatic prediction algorithms assessing the effect of sequence variants on splicing or protein function are irrelevant, which makes variant interpretation challenging to molecular diagnostic laboratories. In order to facilitate and improve clinical diagnostic assessment and genetic counseling, we present i) an update of the previously reported RNU4ATAC mutations and an analysis of the genetic variations affecting this gene using the Genome Aggregation Database (gnomAD) resource; ii) the pathogenicity prediction performances of scores computed based on an RNA structure prediction tool and of those produced by the Combined Annotation Dependent Depletion tool for the 285 RNU4ATAC variants identified in patients or in large-scale sequencing projects; iii) a method, based on a cellular assay, that allows to measure the effect of RNU4ATAC variants on splicing efficiency of a minor (U12-type) reporter intron. Lastly, the concordance of bioinformatic predictions and cellular assay results was investigated.

New insights into minor splicing-a transcriptomic analysis of cells derived from TALS patients.

Minor intron splicing plays a central role in human embryonic development and survival. Indeed, biallelic mutations in RNU4ATAC, transcribed into the minor spliceosomal U4atac snRNA, are responsible for three rare autosomal recessive multimalformation disorders named Taybi-Linder (TALS/MOPD1), Roifman (RFMN), and Lowry-Wood (LWS) syndromes, which associate numerous overlapping signs of varying severity. Although RNA-seq experiments have been conducted on a few RFMN patient cells, none have been performed in TALS, and more generally no in-depth transcriptomic analysis of the ∼700 human genes containing a minor (U12-type) intron had been published as yet. We thus sequenced RNA from cells derived from five skin, three amniotic fluid, and one blood biosamples obtained from seven unrelated TALS cases and from age- and sex-matched controls. This allowed us to describe for the first time the mRNA expression and splicing profile of genes containing U12-type introns, in the context of a functional minor spliceosome. Concerning RNU4ATAC-mutated patients, we show that as expected, they display distinct U12-type intron splicing profiles compared to controls, but that rather unexpectedly mRNA expression levels are mostly unchanged. Furthermore, although U12-type intron missplicing concerns most of the expressed U12 genes, the level of U12-type intron retention is surprisingly low in fibroblasts and amniocytes, and much more pronounced in blood cells. Interestingly, we found several occurrences of introns that can be spliced using either U2, U12, or a combination of both types of splice site consensus sequences, with a shift towards splicing using preferentially U2 sites in TALS patients' cells compared to controls.

Novel diagnostic tool for prediction of variant spliceogenicity derived from a set of 395 combined in silico/in vitro studies: an international collaborative effort.

Variant interpretation is the key issue in molecular diagnosis. Spliceogenic variants exemplify this issue as each nucleotide variant can be deleterious via disruption or creation of splice site consensus sequences. Consequently, reliable in silico prediction of variant spliceogenicity would be a major improvement. Thanks to an international effort, a set of 395 variants studied at the mRNA level and occurring in 5' and 3' consensus regions (defined as the 11 and 14 bases surrounding the exon/intron junction, respectively) was collected for 11 different genes, including BRCA1, BRCA2, CFTR and RHD, and used to train and validate a new prediction protocol named Splicing Prediction in Consensus Elements (SPiCE). SPiCE combines in silico predictions from SpliceSiteFinder-like and MaxEntScan and uses logistic regression to define optimal decision thresholds. It revealed an unprecedented sensitivity and specificity of 99.5 and 95.2%, respectively, and the impact on splicing was correctly predicted for 98.8% of variants. We therefore propose SPiCE as the new tool for predicting variant spliceogenicity. It could be easily implemented in any diagnostic laboratory as a routine decision making tool to help geneticists to face the deluge of variants in the next-generation sequencing era. SPiCE is accessible at (https://sourceforge.net/projects/spicev2-1/).

Mutation screening of MIR146A/B and BRCA1/2 3'-UTRs in the GENESIS study.

Although a wide number of breast cancer susceptibility alleles associated with various levels of risk have been identified to date, about 50% of the heritability is still missing. Although the major BRCA1 and BRCA2 genes are being extensively screened for truncating and missense variants in breast and/or ovarian cancer families, potential regulatory variants affecting their expression remain largely unexplored. In an attempt to identify such variants, we focused our attention on gene regulation mediated by microRNAs (miRs). We screened two genes, MIR146A and MIR146B, producing miR-146a and miR-146b-5p, respectively, that regulate BRCA1, and the 3'- untranslated regions (3'-UTRs) of BRCA1 and BRCA2 in the GENESIS French national case/control study (BRCA1- and BRCA2-negative breast cancer cases with at least one sister with breast cancer and matched controls). We identified one rare variant in MIR146A, four in MIR146B, five in BRCA1 3'-UTR and one in BRCA2 3'-UTR in 716 index cases and 619 controls. Among these 11 rare variants, 7 were identified each in 1 index case. None of the three relevant MIR146A/MIR146B variants affected the pre-miR sequences. The potential causality of the four relevant BRCA1/BRCA2 3'-UTRs variants was evaluated with luciferase reporter assays and co-segregation studies, as well as with bioinformatics analyses to predict miRs-binding sites, RNA secondary structures and RNA accessibility. This is the first study to report the screening of miR genes and of BRCA2 3'-UTR in a large series of familial breast cancer cases. None of the variant identified in this study gave convincing evidence of potential pathogenicity.

BRCA2 Polymorphic Stop Codon K3326X and the Risk of Breast, Prostate, and Ovarian Cancers.

BACKGROUND: The K3326X variant in BRCA2 (BRCA2*c.9976A>T; p.Lys3326*; rs11571833) has been found to be associated with small increased risks of breast cancer. However, it is not clear to what extent linkage disequilibrium with fully pathogenic mutations might account for this association. There is scant information about the effect of K3326X in other hormone-related cancers. METHODS: Using weighted logistic regression, we analyzed data from the large iCOGS study including 76 637 cancer case patients and 83 796 control patients to estimate odds ratios (ORw) and 95% confidence intervals (CIs) for K3326X variant carriers in relation to breast, ovarian, and prostate cancer risks, with weights defined as probability of not having a pathogenic BRCA2 variant. Using Cox proportional hazards modeling, we also examined the associations of K3326X with breast and ovarian cancer risks among 7183 BRCA1 variant carriers. All statistical tests were two-sided. RESULTS: The K3326X variant was associated with breast (ORw = 1.28, 95% CI = 1.17 to 1.40, P = 5.9x10(-) (6)) and invasive ovarian cancer (ORw = 1.26, 95% CI = 1.10 to 1.43, P = 3.8x10(-3)). These associations were stronger for serous ovarian cancer and for estrogen receptor-negative breast cancer (ORw = 1.46, 95% CI = 1.2 to 1.70, P = 3.4x10(-5) and ORw = 1.50, 95% CI = 1.28 to 1.76, P = 4.1x10(-5), respectively). For BRCA1 mutation carriers, there was a statistically significant inverse association of the K3326X variant with risk of ovarian cancer (HR = 0.43, 95% CI = 0.22 to 0.84, P = .013) but no association with breast cancer. No association with prostate cancer was observed. CONCLUSIONS: Our study provides evidence that the K3326X variant is associated with risk of developing breast and ovarian cancers independent of other pathogenic variants in BRCA2. Further studies are needed to determine the biological mechanism of action responsible for these associations.

Mutation analysis of PALB2 gene in French breast cancer families.

Several population-based and family-based studies have demonstrated that germline mutations of the PALB2 gene (Partner and Localizer of BRCA2) are associated with an increased risk of breast cancer. Distinct mutation frequencies and spectrums have been described depending on the population studied. Here we describe the first complete PALB2 coding sequence screening in the French population. We screened the complete coding sequence and intron-exon boundaries of PALB2, using the EMMA technique, to assess the contribution of pathogenic mutations in a set of 835 familial breast cancer cases and 662 unrelated controls from the French national study GENESIS and the Paul Strauss Cancer Centre, all previously tested negative for BRCA1 and BRCA2 pathogenic mutations. Our analysis revealed the presence of four novel deleterious mutations: c.1186insT, c.1857delT and c.2850delC in three cases, c.3418dupT in one control. In addition, we identified two in-frame insertion/deletion, 19 missense substitutions (two of them predicted as pathogenic), 9 synonymous variants, 28 variants located in introns and 2 in UTRs, as well as frequent variants. Truncating PALB2 mutations were found in 0.36% of familial breast cancer cases, a frequency lower than the one detected in comparable studies in other populations (0.73-3.40%). This suggests a small but significant contribution of PALB2 mutations to the breast cancer susceptibility in the French population.

Targeted Sequencing of the Mitochondrial Genome of Women at High Risk of Breast Cancer without Detectable Mutations in BRCA1/2.

Breast Cancer is a complex multifactorial disease for which high-penetrance mutations have been identified. Approaches used to date have identified genomic features explaining about 50% of breast cancer heritability. A number of low- to medium penetrance alleles (per-allele odds ratio < 1.5 and 4.0, respectively) have been identified, suggesting that the remaining heritability is likely to be explained by the cumulative effect of such alleles and/or by rare high-penetrance alleles. Relatively few studies have specifically explored the mitochondrial genome for variants potentially implicated in breast cancer risk. For these reasons, we propose an exploration of the variability of the mitochondrial genome in individuals diagnosed with breast cancer, having a positive breast cancer family history but testing negative for BRCA1/2 pathogenic mutations. We sequenced the mitochondrial genome of 436 index breast cancer cases from the GENESIS study. As expected, no pathogenic genomic pattern common to the 436 women included in our study was observed. The mitochondrial genes MT-ATP6 and MT-CYB were observed to carry the highest number of variants in the study. The proteins encoded by these genes are involved in the structure of the mitochondrial respiration chain, and variants in these genes may impact reactive oxygen species production contributing to carcinogenesis. More functional and epidemiological studies are needed to further investigate to what extent variants identified may influence familial breast cancer risk.

Occurrence of a non deleterious gene conversion event in the BRCA1 gene.

The duplication in the primate lineage of a portion of the breast and ovarian cancer susceptibility gene BRCA1 has created a BRCA1 pseudogene 45 kb away. Non-allelic homologous recombination (NAHR) between BRCA1 and BRCA1P1 has generated recurrent deleterious germ-line 37-kb deletions encompassing the first two exons of BRCA1, accounting for several breast and ovarian cancer families in various populations. In principle, NAHR intermediates resolution could also lead through a non-crossover configuration to interlocus gene conversion (IGC), but none had been described as yet. Here, we report for the first time an IGC event identified in a breast and ovarian cancer family involving exactly the same segment as that involved in the 37-kb deletions. Close examination of the consequences of this IGC event showed that it does not impact BRCA1 expression. Detailed analysis of the regions of homology between BRCA1 and its pseudogene revealed the specificity of the segment where recombination systematically occurs.