# Case Report: Identification of Novel Variants in ERCC4 and DDB2 Genes in Two Tunisian Patients With Atypical Xeroderma Pigmentosum Phenotype

Imen Nabouli, Asma Chikhaoui, Houcemeddine Othman, Sahar Elouej, Meriem Jones, Arnaud Lagarde, Meriem Ben Rekaya, Olfa Messaoud, Mohamed Zghal, Valerie Delague, et al.

► To cite this version:

Imen Nabouli, Asma Chikhaoui, Houcemeddine Othman, Sahar Elouej, Meriem Jones, et al. Case Report: Identification of Novel Variants in ERCC4 and DDB2 Genes in Two Tunisian Patients With Atypical Xeroderma Pigmentosum Phenotype. Frontiers in Genetics, 2021, 12, pp.650639. 10.3389/fgene.2021.650639. pasteur-03550460

HAL Id: pasteur-03550460

https://riip.hal.science/pasteur-03550460

Submitted on 1 Feb 2022

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
# Case Report: Identification of Novel Variants in *ERCC4* and *DDB2* Genes in Two Tunisian Patients With Atypical Xeroderma Pigmentosum Phenotype

OPEN ACCESS

**Edited by:**

Emiliano González Vioque,
University Clinical Hospital of
Santiago, Spain

**Reviewed by:**

Liang-Liang Fan,
Central South University, China
Wei-Ling Tsou,
Wayne State University, United States

*Correspondence:

Houda Yacoub-Youssef
houda.yacoub@pasteur.utm.tn
orcid.org/0000-0001-5351-4775

**Specialty section:**

This article was submitted to
Genetics of Common and Rare
Diseases,
a section of the journal
Frontiers in Genetics

**Received:** 07 January 2021

**Accepted:** 29 April 2021

**Published:** 31 May 2021

**Citation:**

Nabouli I, Chikhaoui A, Othman H,
Elouej S, Jones M, Lagarde A,
Rekaya MB, Messaoud O, Zghal M,
Delague V, Levy N, De
Sandre-Giovannoli A, Abdelhak S and
Yacoub-Youssef H (2021) Case
Report: Identification of Novel Variants
in ERCC4 and DDB2 Genes in Two
Tunisian Patients With Atypical
Xeroderma Pigmentosum Phenotype.
Front. Genet. 12:650639.
doi: 10.3389/fgene.2021.650639

Imen Nabouli<fnref n="1" />, Asma Chikhaoui<fnref n="1" />, Houcemeddine Othman², Sahar Elouej<fnref n="1" />,³,
Meriem Jones<fnref n="1" />,⁴, Arnaud Lagarde³, Meriem Ben Rekaya<fnref n="1" />, Olfa Messaoud<fnref n="1" />,
Mohamed Zghal⁴, Valerie Delague³, Nicolas Levy³,⁵, Annachiara De Sandre-Giovannoli³,⁶,
Sonia Abdelhak<fnref n="1" /> and Houda Yacoub-Youssef<fnref n="1" />*

Xeroderma Pigmentosum (XP) is a rare genetic disorder affecting the nucleotide excision repair system (NER). It is characterized by an extreme sensitivity to sunlight that induces cutaneous disorders such as severe sunburn, freckling and cancers. In Tunisia, six complementation groups have been already identified. However, the genetic etiology remains unknown for several patients. In this study, we investigated clinical characteristics and genetic defects in two families with atypical phenotypes originating from the central region in Tunisia. Clinical investigation revealed mild cutaneous features in two patients who develop multiple skin cancers at later ages, with no neurological disorders. Targeted gene sequencing revealed that they carried novel variants. A homozygous variation in the *ERCC4* gene c.1762G>T, p.V588F, detected in patient XP21. As for patient XP134, he carried two homozygous mutations in the *DDB2* gene c.613T>C, p.C205R and c.618C>A, p.S206R. Structural modeling of the protein predicted the identified *ERCC4* variant to mildly affect protein stability without affecting its functional domains. As for the case of DDB2 double mutant, the second variation seems to cause a mild effect on the protein structure unlike the first variation which does not seem to have an effect on it. This study contributes to further characterize the mutation spectrum of XP in Tunisian families. Targeted gene sequencing accelerated the identification of rare unexpected genetic defects for diagnostic testing and genetic counseling.

**Keywords:** xeroderma pigmentosum, NER defects, skin cancer, ERCC4/XPF, DDB2 gene

<footnote n="1">¹ Laboratoire de Génomique Biomédicale et Oncogénétique, Institut Pasteur de Tunis, LR16IPT05, Université Tunis ElManar, Tunis, Tunisia, ² Faculty of Health Sciences, Sydney Brenner Institute for Molecular Bioscience, University of the Witwatersrand, Johannesburg, South Africa, ³ Aix Marseille Univ, INSERM, MMG, U1251, Marseille, France, ⁴ Service de dermatologie, Hôpital Charles Nicolle, Tunis, Tunisia, ⁵ Departement of Medical Genetics, Assistance Publique Hôpitaux de Marseille, La Timone Children's Hospital, Marseille, France, ⁶ Biological Resource Center (CRB-TAC), Assistance Publique Hôpitaux de Marseille, La Timone Children's Hospital, Marseille, France</footnote>
# INTRODUCTION

Consanguinity plays a relevant role in the emergence of rare genetic diseases especially in North Africa and Middle East countries (Romdhane et al., 2012). Xeroderma pigmentosum (XP) is a rare autosomal recessive disorder characterized by predisposition to cutaneous malignancies. Its prevalence is estimated to 1/250.000 in Europe (Kleijer et al., 2008). However, it is more frequent in North Africa and especially in Tunisia (1/10.000) (Zghal et al., 2005). XP displays a wide clinical and genetic heterogeneity. It results from mutations in eight different genes: XPA through XPG, which encode Nucleotide Excision Repair (NER) genes, and XPV that encodes the Translesion Synthesis (TLS) DNA polymerase eta (Kleijer et al., 2008).

XP prognosis was widely improved due to increased knowledge about the frequent forms of the disease in Tunisia as XP-A, XP-C and XP-V (Ben Rekaya et al., 2009, 2014; Messaoud et al., 2010), Recently, thanks to the use of novel technologies such as high throughput sequencing, extremely rare forms are increasingly being identified as XP-D, XP-E, XP-G (Ben Rekaya et al., 2017; Chikhaoui et al., 2019).

The identification of disease-causing mutations in XP patients has an important impact on patient care, as it provides early dermatological follow up to prevent skin cancers and provides genetic counseling for affected family members. However, genetic diagnosis is difficult to carry out for a disease that shows such genetic heterogeneity. In addition, Sanger sequencing for known recurrent mutations that are common in specific geographic area does not always succeed in finding the genetic causes of atypical forms of the disease.

In this study, we describe clinical features and report novel genetic defects using targeted gene sequencing in two patients with particular forms of XP. One associated to XP-F, phenotype which we identify for the first time in the Tunisian population and the other consists of a double homozygous mutation in the *DDB2* gene associated to XP-E form.

# Case Presentation

## Clinical Presentation

### Patient XP21

Patient XP21 is a 56-year-old woman, second child of healthy first-degree consanguineous parents originating from North-West of Tunisia (**Figure 1A**). When she was 6 years old, she developed achromic maculas on sun-exposed area of her skin. She consulted a doctor at 26, when she developed her fist basal cell carcinoma (BCC). XP21 had a normal psychomotor and mental development with no neurological manifestations. In her last consultation at age 56, when this study was conducted, she had developed several actinic keratosis. Her skin was dry and presented hyperpigmented areas. It is important to note that XP21 did not suffer from acute sunburn reaction nor had any photosensitivity manifestations. In addition, she had a short stature, but no microcephaly or neurological disorders suggestive of Cockayne syndrome were observed (**Table 1**). Only the reported patient and her healthy sister underwent clinical examination in a medical center. As mentioned by the patient, other family members with similar phenotypes

came to our attention during the genetic inquiry about familial history, who presented mainly achromic maculas on their sun exposed area and that's in 2 brothers and one sister (**Figure 1A**). Unfortunately, they were out of reach as they lived far from medical centers.

### Patient XP134

Patient XP134 was a 15-year old girl born from a consanguineous marriage (**Figure 1B**), with parents originating from the central region of Tunisia. The patient did not show any neurological abnormalities. As she was not protected from UV radiations, she presented an important number of actinic keratosis on sun-exposed areas and numerous squamous cell carcinomas and basal cell carcinomas on the face especially around the nose. She presented few hyperpigmented and hypopigmented patterns macules on her skin (**Table 1**). We tried to reach XP134's sister who was described as having similar clinical manifestations as pigmentation problems on sun exposed area, but she was out of reach. The proband's parents as well.

# Genetic Investigations

## Sanger Sequencing for Recurrent Mutations

Genetic pre-screening for recurrent and founder mutations observed in Tunisian XP patients was conducted in both patients XP21 and XP134 using Sanger Sequencing. Firstly, the patients were screened for already known mutations that are recurrent in Tunisian population as: p.Arg228* in XPA gene, p.Val548AlafsX25 in XPC gene and the deletion g.36847_40771del3925 in POLH gene. Secondly, we screened the mutations associated with other mild forms of XP with neurological manifestations as: p.Arg683Gln in the *ERCC2* gene (XP-D phenotype), P.Leu778Pro in *ERCC5* gene (XP-G phenotype) and p.Lys381Argfs*2 mutation in *DDB2* gene (XP-E phenotype). None of these mutations was found.

Therefore, samples underwent targeted gene sequencing to cover other genes of the DNA repair pathways.

## Targeted Gene Sequencing Results

Through targeted gene sequencing of 87 genes involved in NER pathway, we identified in patient XP21a novel homozygous missense variation in *ERCC4*, within exon 8, (NM_005236.2 c.1762G>T p.Val588Phe) (**Figure 2A**). Bioinformatic prediction tools (SIFT, POLYPHEN, mutation taster...) indicated that this mutation was pathogenic (**Table 2**). Moreover, through the prediction algorithm of human splicing finder tools, this variation was suggested to create an exonic splicing silencer (ESS) probably altering the splicing process.

Regarding the patient XP134, we found that he carried a double homozygous missense variation in *DDB2* gene in exon 5 (NM_000107.2: c.613T>C p.Cys205Arg) and c.618C>A p.Ser206Arg (rs759622121) (**Figure 2B**). Both variants were not described in any population database (Exac, GnomAD, 1000 Genomes...) and both are predicted to induce pathogenic effects using several prediction tools (SIFT, POLYPHEN, mutation taster...) (**Table 2**). Concerning the ACMG classification (Varsome), both variants were classed as of uncertain significance (class 3). Furthermore, this particular genetic profile was not
<figure><img src="image_2.png" /><figcaption>**FIGURE 1** | Pedigrees and clinical images describing two XP affected family members (XP21 and XP134) (A,B). Filled symbols represent affected individuals, open symbols represent unaffected individuals, gray symbols are for suspected individuals, the studied pro-band is indicated with an arrow.</figcaption></figure>

previously described in any XP patient. Molecular diagnosis was confirmed by Sanger sequencing for each variation.

## Molecular Modeling of the Mutant Protein Structures

### Modeling XP-F Mutation

The search for a template using BLAST-based approaches was not able to identify an appropriate template. Therefore, we used a fold recognition based approach for such an end. First, we predicted the repartition of the domains of XPF along the sequence. We found that the segment 374-623 matches a hit with a significant E-value of 7.60e-98 corresponding to a helicase insert domain superfamily in accordance with the functional properties of XPF. We then proceeded to model this domain using a threading approach. pGenTHREADER was able to identify several hits with high confidence values. We selected the PDB entry, 1WP9, which corresponds to the crystal structure of Pyrococcus furiosus Hef helicase domain (PubMed: 15642269) as a template to construct the model. Our initial attempt to predict the structure led to a

bad quality model consisting of severe distortions of the rigid secondary structure elements. In fact, the RMSD (Root Mean Square Distance) between the template and the model is of 6 Angstroms. When we inspected the model, we found that the problem is caused by the long exposed loops spanning between the alpha-helices of the protein. We, therefore, proceed by trimming the sequences of these loops (Supplemental Material) to keep at least 4 amino acids from the C and N termini ends. Such an approach might be capable of preserving the original properties of the exposed surface of the target protein. The Ramachandran plot confirmed a decent quality of the model with 90% of the dihedral angles are in the favored domain and 7% in the allowed region. Even with such a model, we could not be able to explain the effect of the mutation. Therefore, we inspected the possibility that the predicted structure of the helicase insert domain could establish interdomain contacts in the XPF protein. We proceeded first by predicting the residue-residue contact map using RaprorX (Wang et al., 2017). Indeed,
<table><caption>**TABLE 1** | Summary of the clinical examination of XP patients.</caption><thead><tr><th></th><th>XP21 (XP-F)</th><th>XP134-1 (XP-E)</th></tr></thead><tbody><tr><td>Age at last observation (year)</td><td>56</td><td>15</td></tr><tr><td>Age at onset of first symptoms (year)</td><td>6</td><td>5</td></tr><tr><td>Sex</td><td>F</td><td>F</td></tr><tr><td>Protection level</td><td>+</td><td>-</td></tr><tr><td>Geographical origin</td><td>Seliana (North West)</td><td>Kairouan (Sbikha)</td></tr><tr><td>Photophobia</td><td>-</td><td>±</td></tr><tr><td>Photosensitivity</td><td>-</td><td>-</td></tr><tr><td>Neurologic phenotype</td><td>-</td><td>-</td></tr><tr><td>Language delay</td><td>-</td><td>-</td></tr><tr><td>Sunburn</td><td>-</td><td>-</td></tr><tr><td>Irritation</td><td>±</td><td>±</td></tr><tr><td>Hyperpigmentation</td><td>+</td><td>+</td></tr><tr><td>Hypopigmented macula</td><td>-</td><td>+</td></tr><tr><td>Achromic macula</td><td>+</td><td>+</td></tr><tr><td>Dry skin</td><td>+</td><td>+</td></tr><tr><td>Cheilit</td><td>+</td><td>±</td></tr><tr><td>Age at consultation for skin cancer</td><td>26 years old for BCC</td><td>15 years old Multiple actinic keratosis, SCC and BCC</td></tr><tr><td>Actinic keratosis</td><td>+</td><td>++</td></tr><tr><td>Benign skin neoplasm, internal cancers</td><td>-</td><td>-</td></tr></tbody></table>

BCC, Basal Cell Carcinoma; SCC, Spino Cell Carcinoma.

the tool predicted a putative interaction of the helicase insert domain with residues 210-216 (EVVEIHV). According to the protein secondary structure prediction, this segment corresponds to a beta-sheet strand (Drozdetskiy et al., 2015). This region is also flanked by two alpha-helices with 6 and 24 residues. The equivalent aligned segment of the template structure corresponds to two helices of 5 and 24 amino acids flanking a beta stand segment interacting with the helicase insert domain.

Therefore, the structure of XPF was modeled in the presence of the 210-216 segment.

The position of V588 (R138 in the structure) corresponds to an alpha-helix (585-LTFVRQ LEIYR-596) rich in charged and polar amino acids. Residues of this helix are exposed to the solvent with only F589 and L591 making part of the hydrophilic core of the protein. We noticed that in the V588 in the wild type protein forms a hydrophobic cluster of residues with F607, A583, and L580. In the case of V588F mutation causes steric clashes with L580 residue due to its large bulky side chain compared to a valine amino acid in the wild type form (Figure 2C).

## Modeling XP-E Mutation

The two mutations, C205R and S206R are located on a beta sheet that forms one of the blades of the WD40 beta-propeller domain. The mutations are close to the protein-DNA interface which are known to rely on electrostatic forces to stabilize the complex. In order to verify if the two mutations can disrupt

the electrostatic properties of the protein-DNA interface, we calculated the pKa for the ionizable side chains at positions 205 and 206 for both the wild type and the mutant forms at physiological pH of 7.4 using PROPKA 3.0 (10.1021/ct100578z). The pKa value of C205 side chain is predicted to be at 8.59. Therefore, this residue is probably deprotonated in the wild type form. Both residues R205 and R206 in the mutant form are positively charged with predicted pKa values of 10.41 and 8.63. To our surprise, the electrostatic potential at the protein DNA interface did not drastically change between the wild type and the mutant forms (Supplementary Data). The R205 ionizable group remains exposed to the central cavity of the WD40 beta-propeller domain. The same configuration is also demonstrated by the side chain of C205 residue in the wild type form. However, for R206, the side chain is bulky compared to S206 amino acid in the wild type form. Its relatively long side chain is unable to interact with the deepest most residues forming the interface of two bladders (Residues V225, V245 and T260) which is not the case for S206 in the wild type form (Figure 2D).

# MATERIALS AND METHODS

## Patients

This study was conducted according to the principles of the declaration of Helsinki and has obtained the ethical approval (IPT/LR05/ProjectPCI/22/2012/v2) from the institutional review board of Institut Pasteur de Tunis.

The two patients (XP21 and XP134) referred to dermatology department of Charles Nicolle hospital with clinical features suggestive of XP. Written informed consent was obtained for genetic analysis from patients or their parents.

## DNA Extraction, Library Construction, Sequencing, and Data Analysis

Genomic DNA was extracted from peripheral blood. Amplicon libraries were prepared from 1,000 ng of genomic DNA from patients. Custom design of DNA repair disorders' panel was performed using SureDesign (Agilent Technologies Inc.), with probes that cover the exons in 87 genes and 15 bp of the surrounding intronic sequences. Library preparation for NGS was done using the Agilent's HaloPlex<sup>HS</sup> workflow as a target enrichment method. Massively parallel sequencing was performed on an Ion Torrent PGM (Thermo Fisher Scientific).

The data from the sequencing runs were analyzed using the Torrent Suite v4.0.2 analysis pipeline, and aligned using TMAPv.3. The output variant call format (VCF) file was then annotated using the in-house VarAft software version 2.5, which is available online (http://varaft.eu/index.php) (Desvignes et al., 2018). Sequence variants were compared with data in Exome Variant Server (http://evs.gs.washington.edu/EVS/), 1000 Genomes Project (http://www.1000genomes.org/), or GnomAD (https://gnomad.broadinstitute.org/), including Exome Aggregation Consortium database (ExAC), Cambridge, MA (http://exac.broadinstitute.org). A number of online tools were used to predict the functional impact and pathogenicity of the variants such as MutationTaster (http://www.mutationtaster.org/), PolyPhen (http://genetics.bwh.harvard.edu/pph2/),
<figure><img src="image_3.png" /><figcaption>**FIGURE 2** | Molecular findings: (A) Electropherogram showing the variation in *ERCC4* gene in exon 8 (c.1762 G>T p.V588F) at a homozygous state in the patient XP21. (B) Electropherogram showing double variation in *DDB2* gene (c.613 T>C p.C205R and c.618 C>A p.S206R) in exon 5 in XP134 patient. (C) Structures of the XPF protein and F588 XPF mutant suggesting a mild distortion in helix domain. (D) DNA DDB2 protein interaction for mutated DDB2 protein suggestion that R206 affect this process.</figcaption></figure>
<table><caption>**TABLE 2** | Prediction scores of the new detected variants in DDB2 and ERCC4 genes.</caption><thead><tr><th>Gene</th><th>Variant</th><th>Prediction tools</th><th>Score</th><th>Significance</th></tr></thead><tbody><tr><td rowspan="6">DDB2</td><td rowspan="6">c.613T>C; p.C205R</td><td>Mutation taster</td><td>0.99</td><td>Disease causing</td></tr><tr><td>Polyphen</td><td>1</td><td>Probably damaging</td></tr><tr><td>SIFT</td><td>2.77</td><td>damaging</td></tr><tr><td>M-CAP</td><td>0.066</td><td>Possibly damaging</td></tr><tr><td>Provean</td><td>-6.13</td><td>deleterious</td></tr><tr><td>MUpro</td><td>-1.40</td><td>Decrease stability</td></tr><tr><td rowspan="6">DDB2</td><td rowspan="6">c.618C&gt;A; p.S206R</td><td>Mutation taster</td><td>0.99</td><td>Disease causing</td></tr><tr><td>Polyphen</td><td>1</td><td>Probably damaging</td></tr><tr><td>SIFT</td><td>2.78</td><td>damaging</td></tr><tr><td>M-CAP</td><td>0.05</td><td>Possibly damaging</td></tr><tr><td>Provean</td><td>-2.9</td><td>deleterious</td></tr><tr><td>MUpro</td><td>-1.05</td><td>Decrease stability</td></tr><tr><td rowspan="6">ERCC4</td><td rowspan="6">c.1762G&gt;T; p.V588P</td><td>Mutation taster</td><td>0.99</td><td>Disease causing</td></tr><tr><td>Polyphen</td><td>1</td><td>Probably damaging</td></tr><tr><td>SIFT</td><td>2.8</td><td>damaging</td></tr><tr><td>M-CAP</td><td>0.177</td><td>Possibly damaging</td></tr><tr><td>Provean</td><td>-5.69</td><td>deleterious</td></tr><tr><td>MUpro</td><td>-1.64</td><td>Decrease stability</td></tr></tbody></table>

SIFT (http://sift.bii.a-star.edu.sg/), M-CAP (http://bejerano.stanford.edu/mcap/), Provean (http://provean.jcvi.org/index.php) and MU Pro for protein stability prediction (http://mupro.proteomics.ics.uci.edu/).

## Sanger Sequencing

The detected variants with a high prediction score of pathogenicity underwent targeted Sanger sequencing for validation using the following primers: **ERCC4F**: 5'GTAAGA TGTCTTCCCTTCGG 3'; **ERCC4R**: 5'CATAAGCAGCATCG TAACGG 3' and **DDB2F**: 5'ATGGAGCAGTCTGAATGTTC 3'; **DDB2R**: 5'CCACTCCTCTAGACAGG3') which cover both variants in DDB2 gene. Direct sequencing was performed using Big Dye terminator technology (ABI 3130), and sequences were analyzed using Bioedit packages. Genomic and protein sequence variants are described following the Human Genome Variations Society Guidelines available at https://varnomen.hgvs.org/. **ERCC4** and **DDB2** variants are, respectively, described relative to the transcript reference sequences NM_005236.2 and NM_ NM_000107.2. ACMG classification of the variants (Richards et al., 2015) was obtained using Varsome at https://varsome.com.

## Protein Modeling

The primary structure of XPF endonuclease was deposited in UiprotKB database under the accession Q92889. The sequence consists of 916 amino acids of which only the structure corresponding to segment 834-916 was solved experimentally. To solve the structure of the region containing the mutation we first searched for conserved domains of the entire sequence. The search for the appropriate template to construct the 3D structure of XPF was performed using the pGenTHREADER program. The sequences of the best hits

were then processed and aligned to the target sequence of XPF using 3D-expresso. The resulting alignment was then used to guide the model building process by MODELER (version 9.15). Twenty conformers were generated by the comparative modeling software of which we selected the best model according to the DOPE and the GA341 scores. Finally, the stereochemical quality of the model was assessed by constructing the Ramachandran plot to dress the repartition of the phi/psi angles in the favored, allowed, and non-tolerated dihedral domains.

The XP-E structure was retrieved from the co-crystal complex of PDB code 4E54 (PMID 22822215).

## DISCUSSION

XP is a rare autosomal recessive disorder, characterized by ubiquitous defects in the DNA repair system. This disorder is especially characterized by sunburn after brief sun exposure and pigmented macules on sun-exposed skin (Krieger and Berneburg, 2012). Targeted gene panels help to focus on a set of relevant candidate genes within the overlapping phenotypes with similar clinical manifestations, to identify the exact genetic etiology while providing cost and time advantages when compared to classical Sanger sequencing.

About 40% of XP patients do not show sunburn reactions (Bradford et al., 2011). As for XP-E, increased number of freckles in sun-exposed areas seems to be the unique symptom for clinical diagnosis of the disease. These maculae are usually ignored and only the development of primary malignant lesions seems to be the reason leading to clinical consultation. This often results in a late diagnosis of patients with mild XP forms.

We report in the region, the XP-F form which is among the rarest complementation groups of XP. To date, 29 XP-F patients were reported worldwide (Ferri et al., 2020). Most of them are Japanese (Tofuku et al., 2015). XP-F (Omim: #278760) is caused by mutations in **ERCC4** gene. Furthermore, several **ERCC4** mutations were associated to a combined XP/CS phenotype, to XFE progeroid syndrome, Cockayne syndrome, Cerebellar ataxia-dominant phenotype or to Fanconi Anemia disease complementation group Q (Kashiyama et al., 2013; Doi et al., 2018). Although the clinical features of patient XP21 indicated a phenotype of Xeroderma pigmentosum, a moderate form of Fanconi anemia type Q was suspected. This hypothesis was quickly dismissed as patient XP21, and despite her short stature, did not develop internal cancers or neurological disorders. Fertility was not explored as the patient was not married.

The protein XPF consists of 916 amino acids. It forms a heterodimeric complex with ERCC1 and plays the role of structure-specific endonuclease 5' incision during the activation of NER repair system upon UV induced DNA damage (Sijbers et al., 1996; Richards et al., 2015). Twenty-seven variations have been reported thus far in the **ERCC4** gene (Zhou et al., 2017). While most mutations in this gene are located in exon 8, suggesting that it may represent a hotspot mutation site, the variants in this area can lead to different disorders.
Indeed, patient XP21 harbored a novel G to T transversion at position c.1762. In the same exon, other variants as p.R589W give rise to combined XP/CS/FA, Fanconi anemia or mild XP forms. In our study, given the late onset of cutaneous cancer and the lack of associated neurological disorders, we suggest that patient XP21 was affected with a mild XP-F form.

The *ERCC4* endonuclease consists of an N-terminal DNA helicases domain (1-457), a domain containing the nuclease active site (656-813) and a C-terminal containing 2 tandem helix-hairpin-helix (HhH) domain (Manandhar et al., 2015). The previously reported *ERCC4* mutations are widely distributed over the entire length of the protein (Bogliolo et al., 2013; Kashiyama et al., 2013; Zhou et al., 2017; Mori et al., 2018). The novel reported variation in this study is located outside the nuclease active site.

The molecular modeling study suggests that the V588F mutation effect is probably caused by the steric clashes between F588 and L580. The latter amino acid belongs to the beta sheet which holds many of the core residues of the helicase insert domain. Of most importance, the 210-216 segment also belongs to this beta-sheet forming the outermost strand. The steric clash between F588 and L580 could cause a local distortion of the beta sheet. While this might not have a significant effect on the interaction of the helicase insert domain with the DNA, the fact that the 210-216 segment is implicated in the interaction network of the central beta sheet, suggests the cooperative functionality established between modeled helicase insert domain and other segments of the protein might be affected by the mutation.

Xeroderma pigmentosum type E (Omim: # 278740) caused by mutations in the *DDB2* gene. The DDB2 subunit is a component of a complex involved in ubiquitin-mediated proteolysis. Consisting of 427 residues, it plays a major role in DNA-damage recognition (Stoyanova et al., 2009). To date, more than 24 XP-E patients have been reported worldwide (Nichols et al., 1996; Itoh et al., 2000; Vahteristo et al., 2007; Oh et al., 2011; Karagün et al., 2020; Yang et al., 2020) and 17 variations were reported in *DDB2* gene (Yang et al., 2020) None of them is from North African countries except 3 Tunisian patients reported previously by our lab (Ben Rekaya et al., 2017). In this study, we report another novel variation in the Tunisian population emphasizing the complex genetic background of North African population.

It is a double homozygous patient harboring two novel mutations in the same exon of *DDB2* gene, responsible for XP-E, which was never described in XP patients. The emergence of rare autosomal genetic variants is frequent in people with common ethnic origin caused by consanguinity and endogamy. However, the appearance of double variations is an extremely rare event. These double mutations could explain the severe clinical phenotype in XP134 patient.

The DDB2 protein is composed from an N-terminal helix-loop helix segment (101 to 136) followed by a 7-bladed WD40 b propeller domain (residues 137 to 454) (Scrima et al., 2008). It is an ubiquitin ligase component of a multimeric complex involved in the degradation of DNA damage-response proteins

(Liu et al., 2009). The two C205R and S206R variants in XP134 are located on a beta sheet that forms one of the blades of the beta-propeller WD40 domain of the DDB2 protein. Mutations close to the protein-DNA interface are known to delay electrostatic forces to stabilize the complex (Scrima et al., 2008), which could explain the severity of the XP134 phenotype compared to previously described XP-E patients (Ben Rekaya et al., 2017).

Our protein modeling study also suggests that the effect of the C205R mutation is not probably responsible for the clinical phenotype since the mutation preserves both, the number of charges carried by the corresponding residue and the electrostatic potential at the protein-DNA interface. Instead, the S206R mutation might have a significant effect on the protein structure. This could be explained by the fact that the mutation causes the R205 positively charged side chain to be buried inside the hydrophobic core of the protein. Such configuration is highly unstable and could cause the local distortion of the protein at the mutation level, but also could have a more severe consequence on the interaction with the damaged DNA.

## CONCLUSION

In this report, we further characterize new XP complementation groups emerging in the Tunisian population and in the region. While the mutational findings have improved the understanding of XP phenotype, it is necessary to continue the molecular investigations for other underdiagnosed patients for whom the discovery of genetic etiology will improve the clinical follow up.

## DATA AVAILABILITY STATEMENT

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

## ETHICS STATEMENT

The studies involving human participants were reviewed and approved by Institut Pasteur Tunisia ethical committee. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin. Written informed consent was obtained from the minor(s)' legal guardian/next of kin for the publication of any potentially identifiable images or data included in this article.

## AUTHOR CONTRIBUTIONS

IN did experiments. AC did the analysis, interpretations of data, and drafted the manuscript. HO did *in silico* protein modeling. SE and AL did experiments, analysis, and interpretation of targeted gene sequencing data. MJ and MZ clinical investigation of patients and family members. MR and OM helped in patients pre-screening for known mutations. VD, AD, NL, and SA study
concept and design. AD, SA, and HY-Y critical revision of the manuscript. HY-Y supervision of the study. All authors contributed to the article and approved the submitted version.

## FUNDING

This work was supported by the Tunisian Ministry of Public Health, the Ministry of Higher Education and Scientific Research (LR16IPT05), the Projet collaborative interne (PCI_Melanoma, IPT), the RARE-MED project (A* MIDEX Initiative d'Excellence Aix Marseille Université) and the TRIM-RD project (Projet N° 19272) from AFM-Téléthon France.

## REFERENCES

Ben Rekaya, M., Laroussi, N., Messaoud, O., Jones, M., Jerbi, M., and Naouali, C. (2014). A founder large deletion mutation in Xeroderma pigmentosum-Variant form in Tunisia: implication for molecular diagnosis and therapy. *Biomed. Res. Int.* 2014:256245. doi: 10.1155/2014/256245

Ben Rekaya, M., Messaoud, O., Talmoudi, F., Nouira, S., Ouragini, H., Amouri, A., et al. (2009). High frequency of the V548A fs X572 XPC mutation in Tunisia: implication for molecular diagnosis. *J. Hum. Genet.* 54, 426–429. doi: 10.1038/jhg.2009.50

Ben Rekaya, M., Naouali, C., Messaoud, O., Jones, M., Bouyacoub, Y., Nagara, M., et al. (2017). Whole Exome Sequencing allows the identification of two novel groups of Xeroderma pigmentosum in Tunisia, XP-D and XP-E: impact on molecular diagnosis. *J. Dermatol. Sci.* 89, 172–180. doi: 10.1016/j.jdermsci.2017.10.015

Bogliolo, M., Schuster, B., Stoepker, C., Derkunt, B., Su, Y., Raams, A., et al. (2013). Mutations in ERCC4, encoding the DNA-repair endonuclease XPF, cause Fanconi anemia. *Am. J. Hum. Genet.* 92, 800–806. doi: 10.1016/j.ajhg.2013.04.002

Bradford, P. T., Goldstein, A. M., Tamura, D., Khan, S. G., Ueda, T., Boyle, J., et al. (2011). Cancer and neurologic degeneration in xeroderma pigmentosum: long term follow-up characterises the role of DNA repair. *J. Med. Genet.* 48, 168–176. doi: 10.1136/jmgb.2010.083022

Chikhaoui, A., Elouej, S., Nabouli, I., Jones, M., Lagarde, A., Ben Rekaya, M., et al. (2019). Identification of a ERCC5 c.2333T>C (L778P) variant in two tunisian siblings with mild xeroderma pigmentosum phenotype. *Front. Genet.* 10:111. doi: 10.3389/fgene.2019.00111

Desvignes, J. P., Bartoli, M., Delague, V., Krahn, M., Miltgen, M., Béroud, C., et al. (2018). VarAFT: a variant annotation and filtration system for human next generation sequencing data. *Nucleic Acids Res* 46, W545–w553. doi: 10.1093/nar/gky471

Doi, H., Koyano, S., Miyatake, S., Nakajima, S., Nakazawa, Y., Kunii, M., et al. (2018). Cerebellar ataxia-dominant phenotype in patients with ERCC4 mutations. *J. Hum. Genet.* 63, 417–423. doi: 10.1038/s10038-017-0408-5

Drozdetskiy, A., Cole, C., Procter, J., and Barton, G. J. (2015). JPred4: a protein secondary structure prediction server. *Nucleic Acids Res.* 43, W389–W394. doi: 10.1093/nar/gkv332

Ferri, D., Orioli, D., and Botta, E. (2020). Heterogeneity and overlaps in nucleotide excision repair disorders. *Clin. Genet.* 97, 12–24. doi: 10.1111/cge.13545

Itoh, T., Yamaizumi, M., Linn, S., and Ono, T. (2000). Reinvestigation of the classification of five cell strains of xeroderma pigmentosum group E with reclassification of three of them. *J. Investigative Dermatol.* 114, 1022–1029. doi: 10.1046/j.1523-1747.2000.00952.x

Karagün, E., Eroz, R., Gamsizkan, M., Baysak, S., Eyup, Y., and Ozcan, Y. (2020). Novel mutation identified in the DDB2 gene in patients with xeroderma pigmentosum group-E. *Int. J. Dermatol.* 59, 989–996. doi: 10.1111/ijd.14957

Kashiyama, K., Nakazawa, Y., Pilz, D. T., Guo, C., Shimada, M., Sasaki, K., et al. (2013). Malfunction of nuclease ERCC1-XPF results in diverse clinical manifestations and causes Cockayne syndrome, xeroderma pigmentosum, and Fanconi anemia. *Am. J. Hum. Genet.* 92, 807–819. doi: 10.1016/j.ajhg.2013.04.007

## ACKNOWLEDGMENTS

The authors would like to thank the patients and their families as well as the patients' support group Helping Xeroderma pigmentosum Children (http://www.xp-tunisie.org.tn/) for their collaboration.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.650639/full#supplementary-material

Kleijer, W. J., Laugel, V., Berneburg, M., Nardo, T., Fawcett, H., Gratchev, A., et al. (2008). Incidence of DNA repair deficiency disorders in western Europe: xeroderma pigmentosum, Cockayne syndrome and trichothiodystrophy. *DNA Repair.* 7, 744–750. doi: 10.1016/j.dnarep.2008.01.014

Krieger, L., and Berneburg, M. (2012). "Pigmentary lesions in patients with increased DNA damage due to defective DNA repair," in *Annales de Dermatologie et de Venereologie* (Tübingen: Elsevier), S130-S134. doi: 10.1016/S0151-9638(12)70124-2

Liu, L., Lee, S., Zhang, J., Peters, S. B., Hannah, J., Zhang, Y., et al. (2009). CUL4A abrogation augments DNA damage response and protection against skin carcinogenesis. *Mol. Cell* 34, 451–460. doi: 10.1016/j.molcel.2009.04.020

Manandhar, M., Boulware, K. S., and Wood, R. D. (2015). The ERCC1 and ERCC4 (XPF) genes and gene products. *Gene* 569, 153–161. doi: 10.1016/j.gene.2015.06.026

Messaoud, O., Ben Rekaya, M., Cherif, W., Talmoudi, F., Boussen, H., Mokhtar, I., et al. (2010). Genetic homogeneity of mutational spectrum of group-A xeroderma pigmentosum in Tunisian patients. *Int. J. Dermatol.* 49, 544–548. doi: 10.1111/j.1365-4632.2010.04421.x

Mori, T., Yousefzadeh, M. J., Faridounnia, M., Chong, J. X., Hisama, F. M., Hudgins, L., et al. (2018). ERCC4 variants identified in a cohort of patients with segmental progeroid syndromes. *Hum. Mutation* 39, 255–265. doi: 10.1002/humu.23367

Nichols, A. F., Ong, P., and Linn, S. (1996). Mutations specific to the xeroderma pigmentosum group E Ddb- phenotype. *J. Biol. Chem.* 271, 24317–24320. doi: 10.1074/jbc.271.40.24317

Oh, K.-S., Emmert, S., Tamura, D., Digiovanna, J. J., and Kraemer, K. H. (2011). Multiple skin cancers in adults with mutations in the XP-E (DDB2) DNA repair gene. *J. Investigative Dermatol.* 131:785. doi: 10.1038/jid.2010.352

Richards, S., Aziz, N., Bale, S., Bick, D., Das, S., Gastier-Foster, J., et al. (2015). Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. *Genet. Med.* 17, 405–424. doi: 10.1038/gim.2015.30

Romdhane, L., Kefi, R., Azaiez, H., Halim, N. B., Dellagi, K., and Abdelhak, S. (2012). Founder mutations in Tunisia: implications for diagnosis in North Africa and Middle East. *Orphanet J. Rare Dis.* 7:52. doi: 10.1186/1750-1172-7-52

Scrima, A., Konícková, R., Czyzewski, B. K., Kawasaki, Y., Jeffrey, P. D., Groisman, R., et al. (2008). Structural basis of UV DNA-damage recognition by the DDB1-DDB2 complex. *Cell* 135, 1213–1223. doi: 10.1016/j.cell.2008.10.045

Sijbers, A. M., De Laat, W. L., Ariza, R. R., Biggerstaff, M., Wei, Y.-F., Moggs, J. G., et al. (1996). Xeroderma pigmentosum group F caused by a defect in a structure-specific DNA repair endonuclease. *Cell* 86, 811–822. doi: 10.1016/S0092-8674(00)80155-5

Stoyanova, T., Roy, N., Kopanja, D., Raychaudhuri, P., and Bagchi, S. (2009). DDB2 (damaged-DNA binding protein 2) in nucleotide excision repair and DNA damage response. *Cell Cycle* 8, 4067–4071. doi: 10.4161/cc.8.24.10109
Tofuku, Y., Nobeyama, Y., Kamide, R., Moriwaki, S., and Nakagawa, H. (2015). Xeroderma pigmentosum complementation group F: report of a case and review of Japanese patients. *J. Dermatol.* 42, 897-899. doi: 10.1111/1346-8138.12942

Vahteristo, P., Kokko, A., Saksela, O., Aittomäki, K., and Aaltonen, L. A. (2007). Blood-derived gene-expression profiling in unravelling susceptibility to recessive disease. *J. Med. Genet.* 44, 718-720. doi: 10.1136/jmg.2007.051342

Wang, S., Sun, S., Li, Z., Zhang, R., and Xu, J. (2017). Accurate *de novo* prediction of protein contact map by ultra-deep learning model. *PLoS Comp. Biol.* 13:e1005324. doi: 10.1371/journal.pcbi.1005324

Yang, R., Kong, Q., Duan, Y., Li, W., and Sang, H. (2020). Identification of a novel DDB2 mutation in a Chinese Han family with Xeroderma pigmentosum group E:a case report and literature review. *BMC Med. Genet.* 21:67. doi: 10.1186/s12881-020-00997-0

Zghal, M., El-Fekih, N., Fazaa, B., Fredj, M., Zhioua, R., Mokhtar, I., et al. (2005). [Xeroderma pigmentosum. Cutaneous, ocular, and neurologic abnormalities in 49 Tunisian cases]. *La Tunisie Med.* 83, 760-763.

Zhou, E. Y., Wang, H., Lin, Z., Xu, G., Ma, Z., Zhao, J., et al. (2017). Clinical and molecular epidemiological study of xeroderma pigmentosum in China: a case series of 19 patients. *J. Dermatol.* 44, 71-75. doi: 10.1111/1346-8138.13576

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Nabouli, Chikhaoui, Othman, Elouej, Jones, Lagarde, Rekaya, Messaoud, Zghal, Delague, Levy, De Sandre-Giovannoli, Abdelhak and Yacoub-Youssef. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.