# A Parent-of-Origin Effect Impacts the Phenotype in Low Penetrance Retinoblastoma Families Segregating the c.1981C>T/p.Arg661Trp Mutation of RB1

Philippine Eloy, Catherine Dehainault, Meriem Sefta, Isabelle Aerts, François

Doz, Nathalie Cassoux, Livia Lumbroso Le Rouic, Dominique

Stoppa-Lyonnet, François Radvanyi, Gaël a. Millot, et al.

► To cite this version:

Philippine Eloy, Catherine Dehainault, Meriem Sefta, Isabelle Aerts, François Doz, et al. A Parent-of-Origin Effect Impacts the Phenotype in Low Penetrance Retinoblastoma Families Segregating the c.1981C>T/p.Arg661Trp Mutation of RB1. PLoS Genetics, 2016, 12 (2), pp.e1005888. 10.1371/journal.pgen.1005888. hal-01286009

HAL Id: hal-01286009

https://hal.sorbonne-universite.fr/hal-01286009

Submitted on 10 Mar 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# A Parent-of-Origin Effect Impacts the Phenotype in Low Penetrance Retinoblastoma Families Segregating the c.1981C>T/p.Arg661Trp Mutation of RB1

Philippine Eloy¹, Catherine Dehainault¹, Meriem Sefta², Isabelle Aerts³, François Doz³,⁴, Nathalie Cassoux⁴,⁵, Livia Lumbroso le Rouic⁵, Dominique Stoppa-Lyonnet¹,⁴,⁶, François Radvanyi², Gaël A. Millot⁷,⁸,⁹, Marion Gauthier-Villars¹, Claude Houdayer¹,⁴,⁶*

1 Service de Génétique, Institut Curie, Paris, France, 2 CNRS UMR144, centre de recherche de l'Institut Curie, Paris, France, 3 Département d'oncologie pédiatrique, adolescents jeunes adultes, Institut Curie, Paris, France, 4 Université Paris Descartes, Sorbonne Paris Cité, Paris, France, 5 Département d'oncologie chirurgicale, service d'Ophtalmologie, Institut Curie, Paris, France, 6 INSERM U830, centre de recherche de l'Institut Curie, Paris, France, 7 Institut Curie, PSL Research University, Paris, France, 8 CNRS UMR 3244, Paris, France, 9 Sorbonne Universités, UPMC Univ Paris 06, Paris, France

* claude.houdayer@curie.fr

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Eloy P, Dehainault C, Sefta M, Aerts I, Doz F, Cassoux N, et al. (2016) A Parent-of-Origin Effect Impacts the Phenotype in Low Penetrance Retinoblastoma Families Segregating the c.1981C>T/ p.Arg661Trp Mutation of *RB1*. PLoS Genet 12(2): e1005888. doi:10.1371/journal.pgen.1005888

**Editor:** Sharon E. Plon, Baylor College of Medicine, UNITED STATES

**Received:** July 28, 2015

**Accepted:** January 30, 2016

**Published:** February 29, 2016

**Copyright:** © 2016 Eloy et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by RETINOSTOP and Programme Incitatif et Coopératif Institut Curie. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

## Abstract

Retinoblastoma (Rb), the most common pediatric intraocular neoplasm, results from inactivation of both alleles of the *RB1* tumor suppressor gene. The second allele is most commonly lost, as demonstrated by loss of heterozygosity studies. *RB1* germline carriers usually develop bilateral tumors, but some Rb families display low penetrance and variable expressivity. In order to decipher the underlying mechanisms, 23 unrelated low penetrance pedigrees segregating the common c.1981C>T/p.Arg661Trp mutation and other low penetrance mutations were studied. In families segregating the c.1981C>T mutation, we demonstrated, for the first time, a correlation between the gender of the transmitting carrier and penetrance, as evidenced by Fisher's exact test: the probability of being unaffected is 90.3% and 32.5% when the mutation is inherited from the mother and the father, respectively (p-value = 7.10⁻⁷). Interestingly, a similar correlation was observed in families segregating other low penetrance alleles. Consequently, we investigated the putative involvement of an imprinted, modifier gene in low penetrance Rb. We first ruled out a *MED4*-driven mechanism by *MED4* methylation and expression analyses. We then focused on the differentially methylated CpG85 island located in intron 2 of *RB1* and showing parent-of-origin-specific DNA methylation. This differential methylation promotes expression of the maternal c.1981C>T allele. We propose that the maternally inherited c.1981C>T/p.Arg661Trp allele retains sufficient tumor suppressor activity to prevent retinoblastoma development. In contrast, when the mutation is paternally transmitted, the low residual activity would mimic a null mutation and subsequently lead to retinoblastoma. This implies that the c.1981C>T mutation is not deleterious *per se* but needs to be destabilized in order to reach pRb haploinsufficiency and initiate tumorigenesis. We suggest that this
phenomenon might be a general mechanism to explain phenotypic differences in low penetrance Rb families.

## Author Summary

Complex genotype-phenotype correlations lead to clinically and emotionally difficult situations. Improved understanding of these correlations is of utmost importance in medical genetics. Low penetrance retinoblastoma families segregating the c.1981C>T / p. Arg661Trp mutation are a good model as germline carriers develop bilateral, unilateral retinoblastoma, benign retinomas or remain unaffected. The c.1981C>T mutation results in a mutant pRb protein that is partially inactivated which may explain the reduced severity observed. However it is still unclear why this would be the case in some family members and not in others. We have demonstrated a parent-of-origin effect in c.1981C>T / p. Arg661Trp pedigrees and have concluded that overexpressed maternally inherited p. Arg661Trp alleles retain sufficient tumor suppressor activity to prevent Rb development. This might be a general phenomenon driving low penetrance retinoblastoma. Our findings shed light on genotype-phenotype correlations in low penetrance retinoblastoma and are of special relevance for genetic counselling.

## Introduction

Retinoblastoma (Rb) is the most common pediatric intraocular neoplasm and occurs in 1 of every 15,000 births. It results from the biallelic inactivation of the *RB1* tumor suppressor gene, located on 13q14 [1]. *RB1* encodes the nuclear phosphoprotein pRB, which plays a prominent role during the G1/S phase transition [2].

In tumors, both *RB1* alleles can be inactivated via diverse mechanisms including point mutations, large rearrangements, promoter hypermethylation and, most frequently, loss of the second allele demonstrated by loss of heterozygosity studies. In non-hereditary retinoblastoma, both *RB1* mutations are somatic and occur in the same retinal cell that develops into a tumor. In contrast, in hereditary retinoblastoma, germline mutation of one allele is associated with predisposition to Rb, while the second mutation on the other allele is somatic, usually acquired during early childhood. Non-hereditary retinoblastomas are usually unilateral (one eye affected) with a median age at diagnosis of 2 years, whereas hereditary cases are usually bilateral (both eyes affected) with a median age at diagnosis of 1 year and an increased risk for second tumors.

Familial hereditary Rb is defined as two or more carriers of an *RB1* germline gene mutation in a family and represents 10% of all retinoblastomas. Predisposition to Rb is transmitted as an autosomal dominant trait with almost complete penetrance (over 90%) and germline carriers usually develop bilateral or multifocal tumors. However, some Rb families display low penetrance (unaffected carriers) and variable expressivity (carriers develop bilateral, unilateral Rb or even benign retinomas). Deciphering the mechanisms underlying low penetrance (LP) Rb is of utmost importance, as it will impact the clinical management of these families and furthers our understanding of Rb as a disease.

The well-known c.1981C>T / p.Arg661Trp low penetrance mutation in exon 20 of *RB1* results in a mutant pRb that is partially inactivated [3,4] which may explain the reduced severity observed. However, why this would be the case in some family members and not in others
remains unclear. Based on the collection of large families, we have demonstrated, that in the context of c.1981C>T/p.Arg661Trp low penetrance, a parent-of-origin effect impacts on Rb phenotype. When the mutation is inherited from the paternal side, offspring are retinoblastoma-prone. In contrast, when the mutation is inherited from the maternal side, offspring mostly remain unaffected. Based on these observations, the involvement of a putative modifier, imprinted gene should be considered. Two alternative hypotheses were tested. Firstly, we searched for a possible involvement of the *MED4* gene, located in the flanking centromeric region of *RB1*, as we have recently demonstrated that *MED4* expression is required for Rb development [5]. We postulated maternal imprinting for *MED4*, which results in monoallelic expression from the paternal allele. As a result, when the *RB1* c.1981C>T/p.Arg661Trp mutation is inherited from the mother, loss of the contralateral paternal allele in the tumor would switch off *MED4* expression and prevent retinoblastoma development in the context of a low penetrance mutation. Secondly, we focused on a differentially methylated CpG island showing parent-of-origin-specific DNA methylation at the *RB1* gene and located in *RB1* intron 2 (called CpG85 hereafter) [6,7]. Differential methylation of CpG85 skews *RB1* expression in favor of the maternal allele [6]. Our results on a series of germline, tumor DNAs and RNAs did not support any involvement of *MED4* in the low penetrance phenotype, but confirmed the differentially methylated status of *RB1* CpG85. It was therefore concluded that overexpressed maternally inherited p.Arg661Trp alleles retain sufficient tumor suppressor activity to prevent Rb development. On the other hand, when the mutation is paternally transmitted, the low residual activity would mimic a null mutation, leading to haploinsufficiency and Rb development.

# Results

## Description of the families

We reviewed the records of 49 pedigrees from Institut Curie with a family history of Rb. Thirty-four of these families segregated high penetrance mutations and 15 families segregated low penetrance mutations. Eight low penetrance families derived from the literature were also found by PubMed search and were added to the study (Table 1). All first generation carriers were excluded to avoid any bias in DER calculation (Disease Eye Ratio, see “Patients and Methods” section).

**High penetrance families.** Unilateral Rb and bilateral Rb were identified in 7 (3 males and 4 females) and 71 patients (34 males and 37 females), respectively. Mean DER was 1.94.

**p.Arg661Trp low penetrance families.** Six families were derived from Institut Curie and 4 families were derived from the literature, corresponding to 85 germline carriers: 49 males and 36 females. Unilateral Rb was identified in 21 patients (12 males and 9 females), bilateral Rb was identified in 12 patients (9 males and 3 females) and retinoma was identified in 2 patients. The remaining 50 carriers were unaffected. Mean DER was 0.60. One pedigree segregating the c.1981C>T/p.Arg661Trp mutation is shown (Fig 1).

**Other low penetrance families.** We found 9 families from Institut Curie and 6 families from the literature segregating a non-c.1981C>T/p.Arg661Trp low penetrance *RB1* mutation (i.e. c.-193T>G, c.10A>T, c.19del, c.43_65dup, c.45_79dup, c.607+1G>T, c.862-10T>C, c.1331A>G, c.1422-2A>G, c.1960G>A, c.1960G>C), corresponding to 127 carriers of a germline mutation. Unilateral Rb and bilateral Rb were identified in 34 and 15 patients, respectively, and 68 carriers were unaffected. Mean DER was 0.53

## Statistical analysis

The parental origin of the c.1981C>T/p.Arg661Trp mutant allele was documented in 71 of the 85 carriers. In this series of 71 carriers, 31 and 40 received the mutant allele from their mother
<table><caption>**Table 1. Description of low penetrance families.** DER: disease-eye ratio (see text for details). Nomenclature follows HGVS rules using the reference sequence NM_000321.2. Previously published families are indicated. Pedigrees F6, F7, F16, F17, F20-22 were from our series and have been published in part (see text for details). Families F14 and F15 were removed from statistical analysis (see text for details).</caption><thead><tr><th>Family</th><th>Mutation description</th><th>Expected consequence</th><th>Number affected</th><th>Total number of carriers</th><th>DER</th><th>Comments</th></tr></thead><tbody><tr><td>F1[8]</td><td>c.1981C&gt;T</td><td>p.Arg661Trp</td><td>5</td><td>6</td><td>1</td><td>4 Unilateral Rb<br>1 Bilateral Rb</td></tr><tr><td>F2[8]</td><td>c.1981C&gt;T</td><td>p.Arg661Trp</td><td>3</td><td>5</td><td>0.6</td><td>3 Unilateral Rb</td></tr><tr><td>F3[9]</td><td>c.1981C&gt;T</td><td>p.Arg661Trp</td><td>7</td><td>18</td><td>0.56</td><td>3 Bilateral Rb<br>4 Unilateral Rb</td></tr><tr><td>F4[10]</td><td>c.1981C&gt;T</td><td>p.Arg661Trp</td><td>6</td><td>10</td><td>0.8</td><td>2 Bilateral Rb<br>2 Unilateral Rb<br>2 retinomas</td></tr><tr><td>F5</td><td>c.1981C&gt;T</td><td>p.Arg661Trp</td><td>2</td><td>11</td><td>0.18</td><td>2 Unilateral Rb</td></tr><tr><td>F6[11]</td><td>c.1981C&gt;T</td><td>p.Arg661Trp</td><td>4</td><td>8</td><td>1</td><td>4 Bilateral Rb</td></tr><tr><td>F7[11]</td><td>c.1981C&gt;T</td><td>p.Arg661Trp</td><td>5</td><td>18</td><td>0.33</td><td>4 Unilateral Rb<br>1 Bilateral Rb</td></tr><tr><td>F8</td><td>c.1981C&gt;T</td><td>p.Arg661Trp</td><td>1</td><td>2</td><td>0.5</td><td>1 Unilateral Rb</td></tr><tr><td>F9</td><td>c.1981C&gt;T</td><td>p.Arg661Trp</td><td>1</td><td>3</td><td>0.33</td><td>1 Unilateral Rb</td></tr><tr><td>F10</td><td>c.1981C&gt;T</td><td>p.Arg661Trp</td><td>1</td><td>3</td><td>0.67</td><td>1 Bilateral Rb</td></tr><tr><td>F11[12]</td><td>c.1960G&gt;C</td><td>p.Val654Leu</td><td>7</td><td>16</td><td>0.44</td><td>7 Unilateral Rb</td></tr><tr><td>F12</td><td>c.1960G&gt;A</td><td>p.Val654Met</td><td>1</td><td>4</td><td>0.5</td><td>1 Bilateral Rb</td></tr><tr><td>F13</td><td>c.10A&gt;T</td><td>p.Lys4*</td><td>1</td><td>3</td><td>0.33</td><td>1 Unilateral Rb</td></tr><tr><td>F14[13]</td><td>c.607+1G&gt;T</td><td>Exon 6 skipped</td><td>13</td><td>25</td><td>0.84</td><td>5 Unilateral Rb<br>8 Bilateral Rb</td></tr><tr><td>F15[13]</td><td>c.607+1G&gt;T</td><td>Exon 6 skipped</td><td>3</td><td>10</td><td>0.4</td><td>2 Unilateral Rb<br>1 Bilateral Rb</td></tr><tr><td>F16[11]</td><td>c.607+1G&gt;T</td><td>Exon 6 skipped</td><td>2</td><td>5</td><td>0.4</td><td>2 Unilateral Rb</td></tr><tr><td>F17[11]</td><td>c.607+1G&gt;T</td><td>Exon 6 skipped</td><td>3</td><td>5</td><td>1</td><td>1 Unilateral Rb<br>2 Bilateral Rb</td></tr><tr><td>F18</td><td>c.1696-2A&gt;G</td><td></td><td>2</td><td>4</td><td>0.5</td><td>2 Unilateral Rb</td></tr><tr><td>F19[14]</td><td>c.1331A&gt;G</td><td>Exon 13 skipped</td><td>2</td><td>8</td><td>0.25</td><td>2 Unilateral Rb</td></tr><tr><td>F20[11]</td><td>c.45_79dup</td><td>p.Pro27Leufs*50</td><td>1</td><td>6</td><td>0.17</td><td>1 Unilateral Rb</td></tr><tr><td>F21[11]</td><td>c.1422-2A&gt;G</td><td>Exon 16 skipped</td><td>3</td><td>4</td><td>1</td><td>2 Unilateral Rb<br>1 Bilateral Rb</td></tr><tr><td>F22[11]</td><td>c.-193T&gt;G</td><td>Promoter</td><td>2</td><td>3</td><td>1</td><td>1 Unilateral Rb<br>1 Bilateral Rb</td></tr><tr><td>F23</td><td>c.19del</td><td>p.Arg7Glufs*58</td><td>1</td><td>5</td><td>0.2</td><td>1 Unilateral Rb</td></tr><tr><td>F24[15]</td><td>c.43_65dup</td><td>p.Pro23Leufs*50</td><td>4</td><td>10</td><td>0.4</td><td>3 Unilateral Rb<br>1 retinoma</td></tr><tr><td>F25[16]</td><td>c.862-10T&gt;C</td><td>Exon 9 skipped</td><td>4</td><td>9</td><td>0.55</td><td>3 Unilateral Rb<br>1 Bilateral Rb</td></tr></tbody></table>

doi:10.1371/journal.pgen.1005888.t001

and father, respectively. Twenty-eight carriers who received the mutant allele from their mother remained unaffected (28/31, 90.3%), and only 3 developed Rb (3/31, 9.7%). In contrast, 13 carriers who received the mutant allele from their father remained unaffected (13/40, 32.5%) and 27 developed Rb (27/40, 67.5%). Consequently, inheriting the c.1981C>T/p.Arg661Trp mutation from the maternal side significantly prevented Rb development (p-value = 7.10⁻⁷, Fisher's exact test). In other words, the probability of being unaffected when the
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Family F7 segregating the *RB1* c.1981C>T/p.Arg661Trp mutation. Genotype is provided for tested members as m/n for heterozygous carriers and n/n for homozygous wild-type. OC indicates obligate carriers. Blackened symbols: bilateral Rb; half-blackened symbols: unilateral Rb; dotted symbols: unaffected carriers; dashed symbols: deceased.</figcaption></figure>

doi:10.1371/journal.pgen.1005888.g001

mutation is inherited from the maternal side is 90.3% versus only 32.5% when the mutation is
inherited from the paternal side.

We then looked for a similar disequilibrium in families segregating non-p.Arg661Trp low
penetrance mutant alleles (see Table 1). To avoid any bias, families F14 and F15 segregating
the c.607+1 G>T mutation were excluded from analysis since a parent of origin effect was pre-
viously described [13]. The parental origin of the mutant alleles was documented in 58 of the
82 remaining carriers. Seventeen carriers received the mutation from their mother and 41
received the mutation from their father. Thirteen carriers who received the mutant allele from
their mother were unaffected (13/17, 76.4%) and 4 developed Rb (4/17, 23.6%). Eighteen carri-
ers who received the mutation from their father were unaffected (18/41, 43.9%) and 23 devel-
oped Rb (23/41, 56.1%). Fisher's exact test demonstrated a disequilibrium between the gender
of the transmitting carrier parent and penetrance (p-value = 0.041). Lastly, families segregating
high penetrance mutations displayed no such correlation, as all 54 mutation carriers of known
parental origin developed retinoblastoma, regardless of the gender of the transmitting carrier.
As previously described, no preferential transmission of mutant or normal alleles from carrier
fathers or mothers was observed [17]. These results unambiguously demonstrate that, in the
context of low penetrance Rb, a parent-of-origin effect impacts on Rb phenotype.

RB1 CpG85 methylation analyses

**Blood samples.** To determine whether *RB1* CpG85 is differentially methylated in a par-
ent-of-origin-specific manner, we studied the methylation pattern of 9 CpG dinucleotides
within the CpG85 island using bisulfite treatment and pyrosequencing in DNAs extracted
from blood. For all non-deleted *RB1* samples (i.e. with 2 *RB1* alleles), the C-to-T ratio at the
CpG dinucleotides studied was close to 1:1, indicating 50% methylation at CpG85 (S1A Fig).
We then studied 6 Rb patients with a large *RB1* deletion of known parental origin. All 3
patients with loss of the maternal *RB1* allele showed absence of methylation at CpG85 (S1B
Fig). In contrast, CpG85 was fully methylated in all 3 patients with loss of the paternal *RB1*
allele (S1C Fig). These results confirmed that *RB1* is imprinted and that CpG85 is specifically
methylated on the maternal allele.

**Tumor samples.** To assess a putative imprinting defect at CpG85, the methylation level of 3 CpG islands within *RB1* was analysed in 2 normal retinas and 45 tumors. CpG106 is located in the promoter region, while CpG85 and CpG42 are located in intron 2. In normal retina DNAs, CpG106 was hypomethylated, while a high DNA methylation level was observed for CpG42 (Fig 2A). CpG85 displayed approximately 50% methylation (Fig 2A), in agreement with genomic imprinting at this locus. Interestingly, CpG85 was fully methylated in all but 3 tumor DNA samples (Fig 2B). These results strongly suggest, for the first time, loss of imprinting at CpG85 locus in retinoblastoma.

**RB1 allelic imbalance.** To assess specific expression imbalance of *RB1* according to the sex of the transmitting parent, a quantitative SNaPshot assay targeting the c.1981C>T/p. Arg661Trp mutation was used within a series of 20 carriers including the low penetrance family F5. To avoid any bias due to a putative exon skipping, exon 20 inclusion was first confirmed by a dedicated RNA study (S2 Fig). Next, using the mutated allele as a marker, we found allelic imbalance in favour of the maternal allele in all 20 patients, albeit to different extend (Table 2 and Fig 3). Surprisingly, in family F5, unaffected carriers 1, 4, 5, 6 and 8 showed higher expression of the mutant allele whereas affected probands 3 and 7 showed higher expression of the wild type allele. Allelic ratio was close to equilibrium for the unaffected carrier 2 and in favour of the mutant allele for the unaffected carrier 9. A similar expression pattern was found in the other c.1981C>T carriers as all affected individuals showed a lower expression of the mutant allele (Table 2). These results confirmed the higher expression of the maternally transmitted *RB1* allele but raised questions about genotype-phenotype correlation as higher expression of the mutant allele and lack of penetrance appeared to be linked.
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Methylation analyses of *RB1* CpG islands using methylation array. X axis represents the position on chromosome 13. Y axis represents overall methylation level. CpG106 localizing in *RB1* promoter is shown in green, CpG42 is shown in pink and CpG85 is shown in blue. For each sample, multiple CpGs are located within an island and each dot represents a single result. A: Normal retina. CpG85 showing approximately 50% of methylation. B: Tumor sample. CpG85 displaying a hypermethylated profile.</figcaption></figure>

doi:10.1371/journal.pgen.1005888.g002

**MED4 analyses.** To determine whether *MED4* CpG53 is differentially methylated in a parent-of-origin-specific manner, we used bisulfite treatment and pyrosequencing on 24 DNA samples extracted from blood. CpG53 from all 24 samples showed fully unmethylated pattern. Similarly, the 45 tumor DNAs and normal retina analyzed by bisulfite treatment and methylation array were found to be unmethylated at the CpG53 locus (S3 Fig). These results excluded a parent-of-origin regulation of *MED4* via *MED4* promoter differential methylation.

To assess whether *MED4* is imbalanced in a parent-of-origin manner, we looked for monoallelic expression using the *MED4* rs41284209 SNP as a marker. We identified 9 patients who were heterozygous carriers of the rs41284209 SNP. No significant allelic disequilibrium was detected by Sanger sequencing and SNaPhot analyses (S4 and S5 Figs). Overall these results showed that both parental alleles contributed equally to *MED4* expression and ruled out *MED4* parental imprinting.

Discussion

Deciphering the molecular basis of low penetrance retinoblastoma is of utmost importance for both researchers and clinicians, as it will shed light on retinoblastoma development, allow prognostic assessment in low penetrance families, and promote optimal genetic counseling and ophthalmological surveillance. In this study, we have identified, for the first time, a parent-of-origin effect in families segregating the c.1981C>T/p.Arg661Trp mutation. In these families, the probabilities of being unaffected for germline carriers were 90.3% and 32.5% when the mutation was inherited from the maternal and paternal side, respectively. Interestingly, a
<table><caption>**Table 2. Expression imbalance in 20 carriers of the c.1981C>T/p.Arg661Trp mutation.** Transmission in family F5 is detailed Fig 3. First degree relatives are indicated for the other families. See text for ratio calculation. (*) See Fig 3.</caption><thead><tr><th>Family</th><th>Patient</th><th>Carrier status</th><th>Parental origin of the c.1981C>T allele</th><th>Ratio c.1981C>T /WT</th></tr></thead><tbody><tr><td>F5</td><td>1</td><td>Unaffected</td><td>Maternal*</td><td>1.48</td></tr><tr><td>F5</td><td>2</td><td>Unaffected</td><td>Paternal*</td><td>0.95</td></tr><tr><td>F5</td><td>3</td><td>Unilateral</td><td>Paternal*</td><td>0.39</td></tr><tr><td>F5</td><td>4</td><td>Unaffected</td><td>Maternal*</td><td>1.85</td></tr><tr><td>F5</td><td>5</td><td>Unaffected</td><td>Maternal*</td><td>1.45</td></tr><tr><td>F5</td><td>6</td><td>Unaffected</td><td>Maternal*</td><td>1.86</td></tr><tr><td>F5</td><td>7</td><td>Unilateral</td><td>Paternal*</td><td>0.21</td></tr><tr><td>F5</td><td>8</td><td>Unaffected</td><td>Maternal*</td><td>1.25</td></tr><tr><td>F5</td><td>6</td><td>Unaffected</td><td>Paternal*</td><td>0.71</td></tr><tr><td>F6</td><td>1</td><td>Bilateral</td><td>Paternal</td><td>0.69</td></tr><tr><td>F6</td><td>2</td><td>Unaffected</td><td>Paternal</td><td>0.46</td></tr><tr><td>F6</td><td>3</td><td>Bilateral</td><td>Paternal (son of F6-2)</td><td>0.40</td></tr><tr><td>F7</td><td>1</td><td>Unaffected</td><td>Maternal</td><td>2.27</td></tr><tr><td>F7</td><td>2</td><td>Unilateral</td><td>Paternal</td><td>0.69</td></tr><tr><td>F7</td><td>3</td><td>Bilateral</td><td>Paternal</td><td>0.62</td></tr><tr><td>F8</td><td>1</td><td>Unilateral</td><td>First generation carrier</td><td>0.38</td></tr><tr><td>F8</td><td>2</td><td>Unaffected</td><td>Maternal (daughter of F8-1)</td><td>1.43</td></tr><tr><td>F9</td><td>1</td><td>Unaffected</td><td>First generation carrier</td><td>0.50</td></tr><tr><td>F9</td><td>2</td><td>Unilateral</td><td>Paternal (daughter of F9-1)</td><td>0.85</td></tr><tr><td>F9</td><td>3</td><td>Unaffected</td><td>Paternal (daughter of F9-1)</td><td>0.91</td></tr></tbody></table>

doi:10.1371/journal.pgen.1005888.t002

similar correlation was observed in families segregating other low penetrance alleles, albeit to a lesser extent: probabilities of being unaffected were 76.5% and 43.9% when the mutation was inherited from the maternal and paternal side, respectively. This finding echoes the maternal protective effect previously described in 2 families (F14 and F15 in this paper) in association with the c.607+1G>T low penetrance mutation [13]. Restoration of the maternal truncated transcript or mutation at an imprinted locus in *cis* were proposed to explain this observation. Our own results on a large number of pedigrees segregating a distinct low penetrance mutation rule out the first hypothesis, but support the second hypothesis.

We have recently shown that retinoblastoma *RB1* -/- cells cannot survive in the absence of *MED4*, both *in vitro* and in orthotopic xenograft models *in vivo*, therefore identifying *MED4* as a survival gene in retinoblastoma [5]. Consequently, we considered a *MED4*-driven general mechanism to explain low penetrance retinoblastoma. We postulated a parent-of-origin regulation of *MED4* that would be able to skew *MED4* expression in favor of the maternal allele. As a result, when the p.Arg661Trp mutation is inherited from the mother, loss of the contralateral paternal allele would dramatically decrease *MED4* expression and prevent retinoblastoma development in the context of a low penetrance mutation. However, methylation and expression studies both ruled out this mechanism to explain the parent-of-origin effect observed in p. Arg661Trp pedigrees.

A recent study demonstrated *RB1* imprinting by a differentially-methylated-region (DMR) at CpG85 in *RB1* intron 2. In humans, this DMR is methylated on the maternal allele and remains unmethylated on the paternal allele. Consequently, CpG85 acts as a weak promoter for an alternative, paternally expressed, *RB1* transcript (*RB1*-E2B) that competes with the main *RB1* transcript. This transcriptional interference skews *RB1* expression in favor of the maternal allele [6,18].
<figure><img src="image_6.png" /><figcaption>**Fig 3.** RB1 allelic imbalance in family F5. The normalized SNaPshot cDNA ratio between the mutant and the wild type alleles are indicated below each carrier individual with corresponding SNaPshot results. The c.1981C>T/p.Arg661Trp mutant allele “T” is indicated in green and the wild type allele “C” is indicated in blue. Dotted symbols: unaffected carriers; half-blackened symbols: unilateral Rb.</figcaption></figure>

doi:10.1371/journal.pgen.1005888.g003

In line with this previous report, our SNAPshot analyses targeting the c.1981C>T/p.Arg661Trp mutation demonstrated higher expression of the maternal *RB1* allele. Our results also demonstrated that, when this mutation is inherited from the maternal side, offspring mostly remain unaffected. Although counter-intuitive, this means that a high level of the c.1981C>T/p.Arg661Trp mutant allele would protect from retinoblastoma. A plausible explanation lies in the residual biochemical properties of p.Arg661Trp mutants, which lack E2F pocket protein-binding activity but retain E2F-independent tumor suppressor function and the wild-type ability to partially suppress colony growth of RB(-) cells and induce parameters of cell differentiation [19]. More broadly, an E2F-independent paradigm of tumor suppression is being developed for *RB1*[20]. Lastly, a study showed that certain LP alleles (p.Arg661Trp included) retain greater functional activity than expected, which is why additional cooperating events are needed to block this residual activity [21]. The competing *RB1*-E2B transcript that lowers *RB1* regular transcript on the paternal allele might constitute this additional event in low penetrance Rb families. Consequently, when the father transmits the mutation, the residual pRb activity is too low to prevent the development of Rb in the cell. The low residual activity
would mimic a null mutation, leading to genomic instability and Rb development. This also means that the c.1981C>T/p.Arg661Trp mutation is not deleterious *per se* but needs to be destabilized in order to reach pRb haploinsufficiency and initiate genomic instability and tumorigenesis [22,23]. Although our results on low penetrance families segregating other LP alleles reached borderline significance (p = 0.041), we propose this hypothesis as a general mechanism to explain disease occurrence in the context of low penetrance Rb.

Intriguingly, we have also reported, for the first time, a hypermethylated, deregulated *RB1* imprint in Rb tumors. Hypermethylation of CpG85 inhibits *RB1*-E2B transcription, therefore enhancing *RB1* main transcript expression. A plausible explanation would be that this loss of imprinting at the CpG85 locus might be used by tumor cells to attempt to increase the expression of pRB and thus restore its tumor suppressor activity.

Overall, we demonstrated that a parent-of-origin effect is involved in low penetrance Rb families segregating the c.1981C>T/p.Arg661Trp mutation of *RB1* and propose this phenomenon as a general mechanism to explain phenotypic differences in low penetrance Rb families.

# Materials and Methods

## Ethics statement

All patients have given written informed consent during genetic counselling sessions. The study was approved by the Groupe Thématique Transverse (GTT) "retinoblastome" of Institut Curie medical center (2013-2310).

## Patients

Institut Curie is the national referral center for retinoblastoma in France. Diagnosis of Rb is established on the basis of examinations by an ophthalmologist and histopathological criteria when treatment involves enucleation. All Rb patients are offered genetic counseling and *RB1* gene mutation analysis in constitutional and tumor DNA. When a germline mutation is found, molecular testing is extended to relatives. Individual written consent for genetic analysis was obtained from all participating patients or their legal guardians. The study was approved by our local ethic committee and retinoblastoma board.

In our series of 1,210 consecutively ascertained cases, we surveyed 49 pedigrees with a family history of Rb. Seven of the low penetrance families have been previously published in part [11]. We included family members for which the mutational status was ascertained by *RB1* analysis and obligate carriers when a DNA sample was not available. Relatives underwent routine fundus examination to look for the presence of retinomas (retinal scars). Since it has been described that retinoma develops after homozygous loss of *RB1* [24], individuals with retinoma were considered to be affected. Obligate carriers with normal fundus examination were considered to be non-penetrant or unaffected. Mutational mosaicism is known to explain the variable expressivity and penetrance in Rb patients [25]. Consequently, we excluded all first-generation carriers of a germline mutation displaying unilateral Rb or remaining unaffected, since retinal mosaicism could not be excluded in these patients. Clinical features included disease status (affected / unaffected) and diseased-eye ratio (DER). The DER is defined as the ratio of the sum of the eyes affected by tumors to the number of mutation carriers in a family. It provides a useful combination of penetrance and expressivity. Families with a DER ≥ 1.5 are considered to display complete penetrance. Families with a DER ≤ 1 are designated as LP [8].
## Statistical analysis

Fisher's exact test was performed using R statistical software v3.0.2 on i) 10 p.Arg661Trp families (6 from our series and 4 from the literature [8,9,10]), ii) 13 non-p.Arg661Trp low penetrance families (9 from our series and 4 from the literature [11,12,14,15,16]), iii) 34 high penetrance families from our series.

## RB1 CpG85 and MED4 CpG53 methylation analyses

**Blood samples and pyrosequencing.** Blood DNAs from 17 Rb patients (including 2 carriers of the p.Arg661Trp mutation and 6 carriers of a large deletion of known parental origin) and 2 controls underwent bisulfite conversion using the EZ DNA Methylation-Gold kit (Zymo research) according to the manufacturer's instructions. Pyrosequencing primers were designed to cover the whole CpG85 island using the PyroMark Assay Design Software v1.0.6 (Qiagen). Bisulfite conversion was assessed by PCR amplification using converted DNA specific primers and agarose gel electrophoresis. Samples were prepared with the PyroMark Q96 Vacuum Workstation and pyrosequencing was performed on a PyroMark Q96 (Qiagen) according to the manufacturer's instructions with subsequent analysis using the Analysis software package v2.5.7 (Qiagen).

**Tumor samples and CGH methylation analyses.** Tumor DNAs from 45 Rb patients (6 bilaterally and 39 unilaterally affected cases, respectively) and DNA samples from 2 normal retina were collected and hybridized on Infinium HumanMethylation 450 BeadChip arrays (Illumina, San Diego, CA). Prior to hybridization, DNA samples underwent bisulfite conversion using the EZ DNA Methylation Kit (Zymo Research). Four microliters of bisulfite-converted DNA were used for hybridization, following the Illumina Infinium HD Methylation protocol. Data were normalized using GenomeStudio (Illumina, Inc.) and R statistical software v3.0.2.

## Expression analyses

**RNA extraction and RT PCR.** Total RNAs were extracted from 200 µL frozen stabilized blood samples with the Nucleospin RNA Blood kit (Macherey-Nagel), according to the manufacturer's instructions. Nagel). RNA quality was controlled using the NanoDrop 1000 Spectrophotometer (Thermo Scientific). Reverse transcription (RT) was performed with random hexamers using the RNA PCR core kit GeneAmp (Applied Biosystems) according to the manufacturer's instructions.

**Allelic imbalance at the rs41284209 MED4 SNP.** To assess a putative MED4 expression imbalance, cDNAs from 9 patients heterozygous for the rs41284209 MED4 SNP were analyzed by Sanger sequencing and SNaPshot assay.

In order to avoid contamination by genomic DNA, we amplified a large cDNA fragment (3658 bp) containing the rs41284209 MED4 SNP and spanning exons 6 to the 3'UTR. Targeted sequencing was then performed using the BigDye Terminator Cycle Sequencing V1.1 Ready Reaction kit (Applied Biosystems) and following electrophoresis in an ABI 3500 Genetic Analyzer (Applied Biosystems). Sequence analyses were performed using Alamut version 2.4 (Interactive Biosoftware, Rouen, France) and FinchTV version 1.4.0 (Geospiza, Inc) softwares. The SNaPshot assay was performed as described below.

**Allelic imbalance at the c.1981C>T RB1 mutation.** Quantitative SNaPshot assay was performed using primers targeting the c.1981C>T RB1 mutation and the SNaPshot quantitative primer extension assay (Applied Biosystems), following a previously detailed protocol [26]. Briefly, to determine whether our assay was able to quantitatively measure allelic imbalance, tumor DNA homozygous for the p.Arg661Trp mutation and wild type DNA were mixed at the following ratios (100:0, 80:20, 60:40, 50:50, 40:60, 20:80 and 0:100), then SNaPshot was successfully tested. The peak ratios were measured between the two allelic versions that is, c.1981C
to c.1981T/p.Arg661Trp. cDNA ratios were then normalized with respect to the values obtained on genomic DNA (cDNA ratios/gDNA ratios) to correct from putative variations in dye incorporation induced by the nucleotide sequence (see S1 Table for an example). All experiments were performed in duplicate.

**Supporting Information**

**S1 Table. Allelic imbalance within a subset of patients from family F5.**

(DOC)

**S1 Fig.** RB1 CpG85 methylation analysis by pyrosequencing in blood samples. The sequence to analyzed is indicated at the top of each pyrogram; Y represents the 9 cytosine residues studied that were either methylated or unmethylated. Bisulfite treatment of DNA converts unmethylated cytosine residues to uracil, whereas 5-methylcytosine residues remain unchanged. Bisulfite-treated DNA sequences will then display a thymine or a cytosine at each CG dinucleotide depending on the methylation status of the cytosine. X axis represents the order of sequential dispensing of enzyme (E), substrate (S) and nucleotides [adenine (A), thymine (T), cytosine (C) and guanine (G)]. Y axis represents peak intensity, which is proportional to the number of dispensed nucleotides incorporated in the sequence. The CG dinucleotides analyzed are shaded on pyrograms. The percentage indicated in colored squares above the corresponding peaks represents the proportion of remaining cytosine residues at the corresponding CG dinucleotide, which in turn indicates the level of methylation at the CG site. The color of the squares above the corresponding peak reflects quality assessment. Yellow represents high quality and blue represents intermediate quality. A. Affected RB1 p.Arg661Trp carrier displaying approximately 50% CpG85 methylation. B. Patient with a large deletion of maternal RB1 allele showing no methylation at CpG85. C. Patient with a large deletion of paternal RB1 allele showing fully methylated CpG85.
(DOC)

**S2 Fig.** Exon 20 RNA analysis. Exon 20 contains the c.1981 C>T mutation. Exon 19/exon 20 junction is indicated at the top of the electropherograms. The c.1981C>T mutation is indicated by an arrow. Panel A: without puromycin. Panel B: with puromycin. Non sense mediated decay inhibition by puromycin didn't reveal any out of frame defect. Targeted RNA analysis showed exon 20 inclusion and absence of skipping.
(TIF)

**S3 Fig.** Methylation analyses of MED4 CpG53 using methylation array. X axis represents the position on chromosome 13. Y axis represents overall methylation level. SUCLA2 and NUDT15 are neighboring genes. MED4 CpG53 is represented in green. A: Normal retina. CpG53 is unmethylated. B: Tumor sample. CpG53 is unmethylated.
(TIF)

**S4 Fig.** MED4 expression analysis in lymphocytes analyzed by SNaPshot assay. The MED4 rs41284209 SNP (c.*783A>G) was used for allelic discrimination. Panel A, genomic results, panel B, cDNA results. No allelic disequilibrium was found.
(TIF)

**S5 Fig.** MED4 expression analysis in lymphocytes analyzed by Sanger sequencing. The MED4 rs41284209 SNP (c.*783A>G) was used for allelic discrimination. Electropherograms of 4 heterozygous carriers displayed no allelic disequilibrium on forward (A) and reverse (B) strands.
(TIF)
Acknowledgments

We thank the patients and their families for their cooperation and support, and Déborah Bourc'his for fruitful discussions.

Author Contributions

Conceived and designed the experiments: CD MGV CH. Performed the experiments: PE CD MS GAM. Analyzed the data: PE CD MS DSL MGV FR GAM CH. Contributed reagents/materials/analysis tools: IA FD NC LLIR DSL MGV FR GAM CH. Wrote the paper: PE CH. Contribution to general study design: PE MS IA FD NC LLIR DSL FR GAM. Critical revising and approval of the MS: CD MS IA FD NC LLIR DSL FR GAM MGV. Statistical analysis: GAM.

References

1. Friend SH, Bernards R, Rogelj S, Weinberg RA, Rapaport JM, et al. (1986) A human DNA segment with properties of the gene that predisposes to retinoblastoma and osteosarcoma. Nature 323: 643-646. PMID: 2877398

2. Chau BN, Wang JY (2003) Coordinated regulation of life and death by RB. Nat Rev Cancer 3: 130-138. PMID: 12563312

3. Otterson GA, Chen W, Coxon AB, Khleif SN, Kaye FJ (1997) Incomplete penetrance of familial retinoblastoma linked to germ-line mutations that result in partial loss of RB function. Proc Natl Acad Sci U S A 94: 12036-12040. PMID: 9342358

4. Otterson GA, Modi S, Nguyen K, Coxon AB, Kaye FJ (1999) Temperature-sensitive RB mutations linked to incomplete penetrance of familial retinoblastoma in 12 families. Am J Hum Genet 65: 1040-1046. PMID: 10486322

5. Dehainault C, Garancher A, Castera L, Cassoux N, Aerts I, et al. (2014) The survival gene MED4 explains low penetrance retinoblastoma in patients with large RB1 deletion. Hum Mol Genet 23: 5243-5250. doi: 10.1093/hmg/ddu245 PMID: 24858910

6. Kanber D, Berulava T, Ammerpohl O, Mitter D, Richter J, et al. (2009) The human retinoblastoma gene is imprinted. PLoS Genet 5: e1000790. doi: 10.1371/journal.pgen.1000790 PMID: 20041224

7. Buiting K, Kanber D, Horsthemke B, Lohmann D (2010) Imprinting of RB1 (the new kid on the block). Brief Funct Genomics 9: 347-353. doi: 10.1093/bfgp/elq014 PMID: 20551090

8. Lohmann DR, Brandt B, Hopping W, Passarge E, Horsthemke B (1994) Distinct RB1 gene mutations with low penetrance in hereditary retinoblastoma. Hum Genet 94: 349-354. PMID: 7927327

9. Abouzeid H, Munier FL, Thonney F, Schorderet DF (2007) Ten novel RB1 gene mutations in patients with retinoblastoma. Mol Vis 13: 1740-1745. PMID: 17960112

10. Onadim Z, Hogg A, Baird PN, Cowell JK (1992) Oncogenic point mutations in exon 20 of the RB1 gene in families showing incomplete penetrance and mild expression of the retinoblastoma phenotype. Proc Natl Acad Sci U S A 89: 6177-6181. PMID: 1352883

11. Taylor M, Dehainault C, Desjardins L, Doz F, Levy C, et al. (2007) Genotype-phenotype correlations in hereditary familial retinoblastoma. Hum Mutat 28: 284-293. PMID: 17096365

12. Hung CC, Lin SY, Lee CN, Chen CP, Lin SP, et al. (2011) Low penetrance of retinoblastoma for p. V654L mutation of the RB1 gene. BMC Med Genet 12: 76. doi: 10.1186/1471-2350-12-76 PMID: 21615945

13. Klutz M, Brockmann D, Lohmann DR (2002) A parent-of-origin effect in two families with retinoblastoma is associated with a distinct splice mutation in the RB1 gene. Am J Hum Genet 71: 174-179. PMID: 12016586

14. Scheffer H, Van Der Vlies P, Burton M, Verlind E, Moll AC, et al. (2000) Two novel germline mutations of the retinoblastoma gene (RB1) that show incomplete penetrance, one splice site and one missense. J Med Genet 37: E6. PMID: 10882758

15. Sanchez-Sanchez F, Ramirez-Castillejo C, Weekes DB, Beneyto M, Prieto F, et al. (2007) Attenuation of disease phenotype through alternative translation initiation in low-penetrance retinoblastoma. Hum Mutat 28: 159-167. PMID: 16988938

16. Lefevre SH, Chauveinc L, Stoppa-Lyonnet D, Michon J, Lumbroso L, et al. (2002) A T to C mutation in the polypyrimidine tract of the exon 9 splicing site of the RB1 gene responsible for low penetrance hereditary retinoblastoma. J Med Genet 39: E21. PMID: 12011162
17. Seminara SB, Dryja TP (1994) Unbiased transmission of mutant alleles at the human retinoblastoma locus. Hum Genet 93: 629-634. PMID: 8005586

18. Anwar SL, Krech T, Hasemeier B, Schipper E, Schweitzer N, et al. Deregulation of RB1 expression by loss of imprinting in human hepatocellular carcinoma. J Pathol 233: 392-401. doi: 10.1002/path.4376 PMID: 24838394

19. Fang W, Mori T, Cobrinik D (2002) Regulation of PML-dependent transcriptional repression by pRB and low penetrance pRB mutants. Oncogene 21: 5557-5565. PMID: 12165854

20. Dick FA, Rubin SM Molecular mechanisms underlying RB protein function. Nat Rev Mol Cell Biol 14: 297-306.

21. Park Y, Kubo A, Komiya T, Coxon A, Beebe K, et al. (2008) Low-penetrant RB allele in small-cell cancer shows geldanamycin instability and discordant expression with mutant ras. Cell Cycle 7: 2384-2391. PMID: 18677112

22. Zheng L, Flesken-Nikitin A, Chen PL, Lee WH (2002) Deficiency of Retinoblastoma gene in mouse embryonic stem cells leads to genetic instability. Cancer Res 62: 2498-2502. PMID: 11980640

23. Castera L, Sabbagh A, Dehainault C, Michaux D, Mansuet-Lupo A, et al. (2010) MDM2 as a modifier gene in retinoblastoma. J Natl Cancer Inst 102: 1805-1808. doi: 10.1093/jnci/djq416 PMID: 21051655

24. Dimaras H, Khetan V, Halliday W, Orlic M, Prigoda NL, et al. (2008) Loss of RB1 induces non-proliferative retinoma: increasing genomic instability correlates with progression to retinoblastoma. Hum Mol Genet 17: 1363-1372. doi: 10.1093/hmg/ddn024 PMID: 18211953

25. Sippel KC, Fraioli RE, Smith GD, Schalkoff ME, Sutherland J, et al. (1998) Frequency of somatic and germ-line mosaicism in retinoblastoma: implications for genetic counseling. Am J Hum Genet 62: 610-619. PMID: 9497263

26. Caux-Moncoutier V, Pages-Berhouet S, Michaux D, Asselain B, Castera L, et al. (2009) Impact of BRCA1 and BRCA2 variants on splicing: clues from an allelic imbalance study. Eur J Hum Genet 17: 1471-1480. doi: 10.1038/ejhg.2009.89 PMID: 19471317