# Combined mutation and rearrangement screening by quantitative PCR high-resolution melting: is it relevant for hereditary recurrent Fever genes?

N. Pallares-Ruiz, L. Philibert, B. Dumont, Aude Fabre, L. Cuisset, E. Cointin, C. Rittore, S. Soler, I. Touitou

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

## ► To cite this version:

N. Pallares-Ruiz, L. Philibert, B. Dumont, Aude Fabre, L. Cuisset, et al.. Combined mutation and rearrangement screening by quantitative PCR high-resolution melting: is it relevant for hereditary recurrent Fever genes?. PLoS ONE, 2010, 5 (11), pp.e14096. 10.1371/journal.pone.0014096. hal-00575068

HAL Id: hal-00575068

https://hal.science/hal-00575068v1

Submitted on 26 May 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Combined Mutation And Rearrangement Screening by Quantitative PCR High-Resolution Melting: Is It Relevant for Hereditary Recurrent Fever Genes?

Nathalie Pallares-Ruiz¹, Laurent Philibert¹, Bruno Dumont¹, Aurélie Fabre¹, Laurence Cuisset², Elodie Cointin¹, Cécile Rittore¹, Stéphan Soler¹, Isabelle Touitou¹,³*

1 Unité Médicale des Maladies Auto-Inflammatoires, Laboratoire de Génétique, Hôpital A de Villeneuve, Montpellier, France, 2 Assistance Publique Hôpitaux de Paris, Hôpital Cochin et Institut Cochin, Université Paris Descartes, INSERM, Paris, France, 3 Université Montpellier 1, Montpellier, France

## Abstract

The recent identification of genes implicated in hereditary recurrent fevers has allowed their specific diagnosis. So far however, only punctual mutations have been identified and a significant number of patients remain with no genetic confirmation of their disease after routine molecular approaches such as sequencing. The possible involvement of sequence rearrangements in these patients has only been examined in familial Mediterranean fever and was found to be unlikely. To assess the existence of larger genetic alterations in 3 other concerned genes, *MVK* (Mevalonate kinase), *NLRP3* (Nod like receptor family, pyrin domain containing 3) and *TNFRSF1A* (TNF receptor superfamily 1A), we adapted the qPCR-HRM method to study possible intragenic deletions and duplications. This single-tube approach, combining both qualitative (mutations) and quantitative (rearrangement) screening, has proven effective in Lynch syndrome diagnosis. Using this approach, we studied 113 unselected (prospective group) and 88 selected (retrospective group) patients and identified no intragenic rearrangements in the 3 genes. Only qualitative alterations were found with a sensitivity similar to that obtained using classical molecular techniques for screening punctual mutations. Our results support that deleterious copy number alterations in *MVK*, *NLRP3* and *TNFRSF1A* are rare or absent from the mutational spectrum of hereditary recurrent fevers, and demonstrate that a routine combined method such as qPCR-HRM provides no further help in genetic diagnosis. However, quantitative approaches such as qPCR or SQF-PCR did prove to be quick and effective and could still be useful after non contributory punctual mutation screening in the presence of clinically evocative signs.

**Citation:** Pallares-Ruiz N, Philibert L, Dumont B, Fabre A, Cuisset L, et al. (2010) Combined Mutation And Rearrangement Screening by Quantitative PCR High-Resolution Melting: Is It Relevant for Hereditary Recurrent Fever Genes? PLoS ONE 5(11): e14096. doi:10.1371/journal.pone.0014096

**Editor:** Jörn Coers, Duke University Medical Center, United States of America

**Received** May 18, 2010; **Accepted** October 27, 2010; **Published** November 23, 2010

**Copyright:** © 2010 Pallares-Ruiz 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.

**Funding:** This work was supported by the French Ministry of Health, the University Hospital of Montpellier, and the Coordination Theme 1 (Health) of the European Community's FP7, Grant Agreement Number HEALTH-F2-2008-200923. 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.

* E-mail: i-touitou@chu-montpellier.fr

## Introduction

The hereditary recurrent fever syndromes are a group of auto-inflammatory disorders characterized by self-limited episodes of fever accompanied by inflammation of the serosal membranes, without apparent infectious etiology [1]. Included in this group are familial Mediterranean fever (FMF) and mevalonate kinase deficiencies (MKD; MIM#610377 for mevalonic aciduria and 260920 for hyper-IgD syndrome), all of which are recessively inherited. Dominantly inherited syndromes have also been described, such as TNF receptor associated periodic syndrome (TRAPS; MIM# 142680) and cryopyrin associated periodic syndromes (CAPS; MIM#120100 for familial cold autoinflammatory syndrome 1, 191900 for Muckle-Wells syndrome and 607115 for CINCA syndrome) [2].

The vast majority (FMF 97%; CAPS 96%; TRAPS 95%; MKD 81%) of the mutations implicated in the hereditary recurrent fever syndromes are punctual substitutions, and have been recorded in the dedicated online database INFEVERS (http://fmf.igh.cnrs.fr/ISSAID/infevers/). However, most laboratories use molecular techniques such as sequence analysis or exon screening (DGGE,

Denaturing Gradient Gel Electrophoresis, dHLPC, denaturing high performance chromatography, HRM, High Resolution Melting analysis) that only allow for the identification of micro-rearrangements. On the other hand, a significant number of apparently mutation-free patients show evocative symptoms of hereditary recurrent fevers. The possible involvement of sequence rearrangements in these diseases has only been examined in FMF and was found to be unlikely [3]. Several methodological approaches have been have been described to search for such gene alterations: SQF-PCR, semi-quantitative fluorescent PCR [4]; MLPA®, multiplex ligation-dependent probe amplification [5]; MAPH, multiplex amplification probe hybridization [5]; qPCR, quantitative PCR [6]. More recently, innovative approaches have been developed that allow the simultaneous, rapid and accurate detection of both micro and macro-rearrangements in genes responsible for hereditary cancers. This strategy, called qPCR-HRM, is based on a real time PCR assay (gene dosage quantitative analysis) combined with HRM curves (qualitative analysis of punctual mutations) in a single tube [7].

To search for sequence rearrangements in our patients, we adapted the qPCR-HRM method for the study of intragenic

<figure><img src="image_3.png" /></figure>
deletions and duplications in the *MVK* (Mevalonate kinase, MIM#251170; NM_000431.1), *NLRP3* (Nod like receptor family, pyrin domain containing 3, MIM#606416; NM_004895.3) and *TNFRSF1A* (TNF receptor superfamily 1A, MIM#191190; NM_001065.2) genes, responsible for MKD, CAPS and TRAPS respectively. We used SQF-PCR as a confirmatory approach.

# Materials and Methods

## Study group

The study includes a total of 201 patients referred for the genetic diagnosis of hereditary recurrent fevers in our reference laboratory specialized in auto-inflammatory disorders (Montpellier, France). Informed consent was obtained from all participating subjects.

The study design consisted in the delineation of two groups of patients (summarized in Figure 1). The prospective study included patients unselectively recruited over a period of 6 months in 2009 (N = 113). The retrospective study included 88 selected patients previously genotyped using our routine diagnosis strategy (screening by sequencing and/or DGGE of hot spot mutations, i.e. *MVK*, all exons; *NLRP3*: exon 3). The selection criteria for MKD (N = 12) was the detection of only a single mutation and for CAPS (N = 76), it was evocative signs, i.e. recurrent fever, urticaria, arthralgia and sometimes neurosensorial deficiency. No TRAPS patients were retrospectively selected since the clinical spectrum for this disease is too unspecific.

## DNA extraction

Human genomic DNA was prepared from whole-blood samples using the QIAmp DNA blood kit (Qiagen, Courtaboeuf, France)

according to the manufacturer's standard protocol. The quality and concentration of the DNA was assessed using a Nanodrop® spectrophotometer (Coleman Technologies, Orlando, FL). DNA working solutions were prepared at a concentration of 5ng/ul for qPCR-HRM and qPCR, and 50ng/ul for SQF-PCR.

## qPCR-HRM and qPCR

**Primer design.** Common primers were used for qPCR-HRM and qPCR for the screening of exons containing mutation hot spots. Specific primers were designed for the other exons (retrospective study). Since effective primer design is an important component of HRM analysis, all amplicons were chosen to contain a single melting domain as assessed using the Poland program (http://www.biophys.uni-duesseldorf.de/local/POLAND/poland.html [8]). The PrimerBlast software (NCBI, http://www.ncbi.nlm.nih.gov/tools/primer-blast/, date of access 03/20/2010) was used to design primers with the same melting temperature (60°C), and to result in amplicons ranging between 100 and 400 bp in size. All primer sequences used in this work are listed in Table 1. The primer pairs were validated if they matched the following conditions: PCR efficiency range between 80-100% as determined using a DNA concentration gradient (40ng-20ng-10ng-5ng-2.5ng), Cycle threshold (Ct) lower than 30 cycles, and all pathological genotypes detected in the validation set.

## Assay conditions

qPCR-HRM and qPCR were performed in a single run using a LightCycler 480 instrument (Roche diagnostics, Meylan, France). Reaction mixtures contained 10 ng genomic DNA, 0.25µM of each primer, 3mM MgCl2 and 1x LightCycler 480 High

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Flow chart for the screening of mutations in the genes *MVK*, *NLRP3* and *TNFRSF1A*. In our patient cohort, we performed a simultaneous search for punctual mutations and gene rearrangements using qPCR-HRM in the prospective group and qPCR alone in the retrospective group (as these patients had already been assessed for qualitative alterations), then SQF-PCR in both groups as a confirmatory quantitative approach.
doi:10.1371/journal.pone.0014096.g001</figcaption></figure>

<figure><img src="image_5.png" /></figure>
<table><caption>**Table 1.** Primers used in this study.</caption><thead><tr><th></th><th>Gene/Exon</th><th>Forward Primer</th><th>Reverse Primer</th><th>Amplicon size (bp)</th><th>Multiplex set</th></tr></thead><tbody><tr><td rowspan="15">qPCR-HRM<sup>1</sup></td><td>B2M/2</td><td>cactgaaaagatgagtatgcc</td><td>Aacattccctgacaatccc</td><td>231</td><td></td></tr><tr><td>DMD/16</td><td>tctatgcaaatgagcaatacacgc</td><td>ggtatcactaacctgtgctgtactc</td><td>290</td><td></td></tr><tr><td>MVK/2</td><td>tcttcctgctggttctgaca</td><td>gcatcaggagaaggggta</td><td>203</td><td></td></tr><tr><td>MVK/3</td><td>tcacctcaggcttattgct</td><td>tgggttctcctccttgcact</td><td>210</td><td></td></tr><tr><td>MVK/4</td><td>ccctctcacccacttgttt</td><td>aggctgaatctggactccttc</td><td>226</td><td></td></tr><tr><td>MVK/5</td><td>cgggagagtacgtttcac</td><td>gacactggccaggtaaggac</td><td>221</td><td></td></tr><tr><td>MVK/6</td><td>ccacctcactgccacag</td><td>ctcttgggcacctaccattg</td><td>221</td><td></td></tr><tr><td>MVK/7</td><td>cctctctccaagtagcacag</td><td>ccctgcacacatcaatgc</td><td>170</td><td></td></tr><tr><td>MVK/8</td><td>tgagttcagtgtggacctgc</td><td>taattgtgtcctggccttcc</td><td>181</td><td></td></tr><tr><td>MVK/9<sup>2</sup></td><td>gaacctcctccctccac</td><td>ttctgagcacagccagattg</td><td>295</td><td></td></tr><tr><td>MVK/10</td><td>aagtggaacagatggaacct</td><td>ccaatgaggaagcaagaacc</td><td>346</td><td></td></tr><tr><td>MVK/11</td><td>gtcaagggtgacctgccttc</td><td>gcctctccagcagtgtcag</td><td>300</td><td></td></tr><tr><td>NLRP3/3a</td><td>gatgtgtgtatactttccccctaa</td><td>aaacctgtcttggtagagtgtcc</td><td>399</td><td></td></tr><tr><td>NLRP3/3c</td><td>agcctcatcagaaagaagctg</td><td>tgaggtcggaactcctcaaac</td><td>469</td><td></td></tr><tr><td>NLRP3/3d</td><td>tctttggctgcagatggaat</td><td>ttgctgagagatcttcaactta</td><td>399</td><td></td></tr><tr><td>NLRP3/3e</td><td>tttcctctttggcctggtaa</td><td>agaagaagctggcgaggaag</td><td>477</td><td></td></tr><tr><td>NLRP3/3c<sup>3</sup></td><td>gcctctgctgcatcacgac</td><td>gctcttGCCactctccatct</td><td>264</td><td></td></tr><tr><td>TNFRSF1A/2</td><td>ttgatggtgtctcctctatctga</td><td>aagaagcagcacccccagac</td><td>240</td><td></td></tr><tr><td>TNFRSF1A/3</td><td>gggctccttccttgtgttct</td><td>cacatagacaggcacccaca</td><td>247</td><td></td></tr><tr><td>TNFRSF1A/4</td><td>agaaatgggtcaggtggagat</td><td>gccagagaggatttgtggtgt</td><td>280</td><td></td></tr><tr><td rowspan="10">qPCR</td><td>MVK/9<sup>2</sup></td><td>accagccgttcctctttttt</td><td>ttctgagcacagccagattg</td><td>230</td><td></td></tr><tr><td>NLRP3/2</td><td>ggggtcctcctctcatgc</td><td>ctagaagcaccaccccccagtc</td><td>242</td><td></td></tr><tr><td>NLRP3/4</td><td>tcgaggctgattttcttttctg</td><td>gcactcacacagatcacatgc</td><td>248</td><td></td></tr><tr><td>NLRP3/5</td><td>tctgatgctttctgcctctgt</td><td>cgctggcagaacttccttag</td><td>244</td><td></td></tr><tr><td>NLRP3/6</td><td>tgactgacattctgccatctct</td><td>tctggtaaagacacccatgaaga</td><td>227</td><td></td></tr><tr><td>NLRP3/7</td><td>ggaacagctgggtactgagg</td><td>cttgtggccactctgccta</td><td>295</td><td></td></tr><tr><td>NLRP3/8</td><td>aatcaccccccctttttgcag</td><td>agagccatcctggatttttga</td><td>274</td><td></td></tr><tr><td>NLRP3/9a</td><td>agtgcaacccaggctttcta</td><td>tcacagagctgtggtcttgg</td><td>240</td><td></td></tr><tr><td>SQF-PCR<sup>4</sup></td><td>GFAP/3</td><td>gaggaaaggattgatggcca</td><td>gaggaggagatccggttctt</td><td>249</td><td>Control (A to H)</td></tr><tr><td></td><td>DMD/28</td><td>tgcattttgaattaccctgctaca</td><td>agtaccaaatagaagcaaatccaaag</td><td>361</td><td>Control (A to H)</td></tr><tr><td></td><td>MVK/2</td><td>tattatgatgggctgaactagg</td><td>tgcctcaggtgtcctttta</td><td>292</td><td>A</td></tr><tr><td></td><td>MVK/3</td><td>cttcttagcacgtggtgcct</td><td>ctctctgtaggctctt</td><td>286</td><td>A</td></tr><tr><td></td><td>MVK/4</td><td>gtcgattttctgttgtctgttgtt</td><td>agagcatgtgcattctccag</td><td>337</td><td>B</td></tr><tr><td></td><td>MVK/5</td><td>ctggaccagatgcttggagt</td><td>aagccacgtccctgtcctg</td><td>336</td><td>A</td></tr><tr><td></td><td>MVK/6</td><td>gagtggaacttgttcttcttgagc</td><td>agactcttgggcacctacca</td><td>328</td><td>B</td></tr><tr><td></td><td>MVK/7</td><td>ttcctgaatggggcaaaat</td><td>ctgcctcctatggtacttccc</td><td>266</td><td>C</td></tr><tr><td></td><td>MVK/8</td><td>gtatcaggtggggcggcttcc</td><td>gaggagacctcgaaaatcc</td><td>279</td><td>A</td></tr><tr><td></td><td>MVK/9</td><td>ggctgtgtgaacacctcctc</td><td>cacagccagattgcagagcca</td><td>314</td><td>B</td></tr><tr><td></td><td>MVK/10</td><td>tctccagccaacaactgtca</td><td>tcaagggaattctccaggtg</td><td>348</td><td>C</td></tr><tr><td></td><td>MVK/11</td><td>ctgggcttttgccctgaat</td><td>aataatccagaaaggggcatc</td><td>304</td><td>C</td></tr><tr><td></td><td>TNFRSF1A/1</td><td>accagccgtgatctctatg</td><td>cactcttccctttgtccctg</td><td>224</td><td>H</td></tr><tr><td></td><td>TNFRSF1A/2</td><td>aggacttgagccagggaagt</td><td>tttccttggggacacacact</td><td>198</td><td>D</td></tr><tr><td></td><td>TNFRSF1A/3</td><td>ctggctgttgtccctagcat</td><td>cacccacacaccactcaaga</td><td>256</td><td>E</td></tr><tr><td></td><td>TNFRSF1A/4</td><td>agaaatgggtcaggtggagat</td><td>ttggttgtcagccacacaga</td><td>268</td><td>D</td></tr><tr><td></td><td>TNFRSF1A/5</td><td>acaaccaacttcctctctggc</td><td>atctgttGCCagctaatgg</td><td>277</td><td>I</td></tr><tr><td></td><td>TNFRSF1A/6</td><td>caccagtgccgtctcttctt</td><td>atagatggatgggtgggatg</td><td>186</td><td>I</td></tr><tr><td></td><td>TNFRSF1A/7</td><td>aacacctgcttgtgtgcag</td><td>acctcttGCCagagtccc</td><td>268</td><td>H</td></tr></tbody></table>

<figure><img src="image_6.png" /></figure>
<table><caption>**Table 1. Cont.**</caption><thead><tr><th>Gene/Exon</th><th>Forward Primer</th><th>Reverse Primer</th><th>Amplicon size (bp)</th><th>Multiplex set</th></tr></thead><tbody><tr><td>TNFRSF1A/8-9</td><td>gggaaatcgacacctgaaaa</td><td>aagctccccctgaaagagag</td><td>233</td><td>H</td></tr><tr><td>TNFRSF1A/9-10</td><td>cccttcagaagtggaggac</td><td>gatcgatctcgtggtcgct</td><td>305</td><td>I</td></tr><tr><td>NLRP3/1</td><td>ccagagccttcagtttggag</td><td>aggagtgtgtcctgagccat</td><td>269</td><td>G</td></tr><tr><td>NLRP3/2</td><td>ccactgtgatatgccaggaa</td><td>gcattccaaagagcaggaac</td><td>208</td><td>F</td></tr><tr><td>NLRP3/3</td><td>gcatctcaggtggatgtgtg</td><td>caggctcagaatgctcatca</td><td>307</td><td>G</td></tr><tr><td>NLRP3/4</td><td>cctcacttccagtttttgcc</td><td>gcactcacacagatcacatgc</td><td>189</td><td>G</td></tr><tr><td>NLRP3/5</td><td>tctgtgtgtggactgaagc</td><td>ctttccccacgacaaacact</td><td>228</td><td>F</td></tr><tr><td>NLRP3/6</td><td>tcagtattgagcaccagcca</td><td>tgaccaaaagtaaccccccatc</td><td>213</td><td>G</td></tr><tr><td>NLRP3/7</td><td>gagtcaaaagcagctgcaca</td><td>ccaccatgtgttctcattgc</td><td>283</td><td>F</td></tr><tr><td>NLRP3/8</td><td>gggatggttaaggggacatt</td><td>caggcccaacctaatcttga</td><td>339</td><td>G</td></tr><tr><td>NLRP3/9a</td><td>tgcaacccaggctttctatt</td><td>tgctgtcattgtcctggtgt</td><td>348</td><td>F</td></tr></tbody></table>

Resolution Melting Master Mix (Roche Diagnostics, Meylan, France) containing Taq polymerase, nucleotides and the dye Resolight in a total volume of 10µl. We used the same cycling conditions for all reactions: an activation step at 95°C for 10 min followed by 50 cycles of 95°C for 15s, a touchdown from 65 to 55°C for 15s (0.5°C/cycle) and 72°C for 20s. Following amplification, products were heated to 95°C for 1 min and cooled to 40°C for 1 min to favor heteroduplex formation. For the melting step, the temperature was raised from 75°C to 97°C with a ramp rate of 0.02°C/s and 25 acquisitions/°C. All reactions were performed at least twice. Each experiment included male or female wild-type DNA (for quantitative calibration and as a qualitative negative control), and positive control DNA with a known point mutation (as a qualitative positive control). All controls are summarized in Table 2.

## Qualitative analysis for punctual mutation screening

Upon completion of the run of the qPCR-HRM (approximately 2 hours), HRM curve analysis was performed using the Light-Cycler 480 Software version 1.5 supplied with the device. The melting curves were normalized and temperature shifted to allow the direct comparison of samples. Difference plots were generated by selecting a negative control, converting the melting profile to a horizontal line and normalizing the melting profiles of the other samples against the test sample. Significant differences in fluorescence from the horizontal baseline were indicative of mutations. Differences were judged as significant if the replicates fell outside the range of variation seen in the wild type samples.

## Quantitative analysis for sequence rearrangement screening

Standard curves were generated for all exons investigated (Figure 1) and for two reference genes (exon 16 of the X-linked gene *DMD*, and exon 2 of *B2M*) to derive sample and calibrator (negative control) relative DNA amounts. The ratio (R) representing the relative copy number of the target exon in the patient was automatically calculated by the LightCycler 480 data analysis software (LC480; Roche) according to the formula:

$$R = \frac{[\text{target exon concentration/reference exon concentration}]_{\text{patient}}}{[\text{target exon concentration/reference exon concentration}]_{\text{control}}}$$

The gene *DMD* was analyzed in each experiment as a control of correct gene dosage: the ratio *DMD/B2M* should be equal to two for female samples, and one for male samples.

## SQF-PCR

**Primer design.** The method of SQF-PCR is based on the comparison of fluorescent profiles of multiplex PCR fragments obtained from different samples (patients and control)<fnref n="4" />. Primers labeled with 6FAM fluorochrome (Table 1) for exons 2–10 of *MVK*, exons 1–10 of *TNFRSFA1*, exons 1–9 of *NLRP3*, and 2 internal control genes (exon 28 of *DMD* and exon 3 of *GFAP*) were designed using the PrimerBlast software (http://www.ncbi.nlm.nih.gov/tools/primer-blast/). We chose allele pairs with the same annealing temperature to perform multiplex PCR resulting in fragments of sizes differing by at least 20bp (186–361 bp). We checked that all exons were adequately co-amplified.

## Assay conditions

The PCR reactions were performed using the QIAGEN Multiplex kit (Qiagen, Courtaboeuf, France) in a total volume of 20 µl containing 100 ng of genomic DNA and 5 to 20 pMol of each primer. The amplification was stopped at the exponential phase (22 cycles). The PCR cycling conditions used were those suggested in the manufacturer's instructions. The PCR products were then separated on an ABI3130xL DNA Sequencer (Applied Biosystems, Foster city, CA, USA).

## Semi-quantitative analysis

Electrophoregrams were analyzed using the GeneMapper software v4.0 (Applied Biosystems, Foster City, CA, USA). In each experiment, the genes *GFAP* and *DMD* were analyzed to determine the relative gene dosage for *MVK*, *NLRP3* and *TNFRSFA1*. To calculate the relative copy number (R) of each target region, each peak value including height and area from the patient samples was compared with those from the normal control

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

<footnote n="1">¹All primers (except HRM-qPCR MVK/9 primer) showed a PCR efficiency range between 80 and 100% (IC+/−0.2%).</footnote>

<footnote n="2">²qPCR and HRM of MVK exon 9 could not be performed in the same run due to low PCR efficiency. We used therefore two different primer pairs.</footnote>

<footnote n="3">³p.T348M allele specific primer.</footnote>

<footnote n="4">⁴Forward primers were 6FAM labeled. Nine distinct PCR sets (Multiplex A to I) were designed, each yielding a pattern of 3 to 6 peaks including the 2 control fragments (DMD and GFAP).</footnote>

<footnote>doi:10.1371/journal.pone.0014096.t001</footnote>
<table><caption>**Table 2.** Controls used in this work.</caption><thead><tr><th colspan="4">qPCR-HRM</th><th colspan="2">SQF-PCR</th></tr><tr><th colspan="3">HRM (Qualitative)</th><th>qPCR (Quantitative)</th><th colspan="2">(Semi-Quantitative)</th></tr></thead><tbody><tr><th colspan="6">Validation Set</th></tr><tr><td>Negative Controls</td><td>N = 48</td><td>Wild-type DNAs</td><td>N = 10</td><td>Asymptomatic individual</td><td>N = 10</td><td>Asymptomatic individual</td></tr><tr><td>Positive Controls</td><td>N = 49</td><td>DNAs with known punctual mutations:</td><td>N = 3</td><td>DNAs with or mimicking a CNV†:</td><td>N = 1</td><td>DNA with a deletion</td></tr><tr><td></td><td>MVK (N = 19)</td><td>p.S52N<sup>(1) (2)</sup>,c.371+8C>T<sup>(1) (2)</sup>,<br/>p.S135L<sup>(1)</sup>, p.S135S<sup>(1)</sup>,<br/>p.D170D<sup>(1) (2)</sup>,c.632-18A>G<sup>(1) (2)</sup>,<br/>c.769-7_769-<br/>6dupT<sup>(1) *</sup>,c.885+24G>A<sup>(1) (2)</sup>,<br/>p.G309S<sup>(1)</sup>, p.S329N<sup>(1)</sup>, p.V377I<sup>(1) (2)</sup>,<br/>p.R388X<sup>(1)</sup>, p.E296GfsX14<sup>(1)</sup></td><td>MVK (N = 1)</td><td>Chr12q24.1del</td><td>MVK<br/>(N = 1)</td><td>Chr12q24.1<br/>del</td></tr><tr><td></td><td>NLRP3 (N = 23)</td><td>c.278-45T>C<sup>(1)</sup>, p.Y141Y<sup>(1) *</sup>,<br/>p.R168Q<sup>(3)</sup>, p.V198M<sup>(1)</sup>, p.T219T<sup>(1)</sup>,<br/>p.A242A<sup>(1) (2) *</sup>, p.P340P<sup>(1) (3)</sup>,<br/>p.L344L<sup>(1) (3)</sup>, p.A439P<sup>(1)</sup>, p.H463H<sup>(1)</sup>,<br/>p.T587I<sup>(1)</sup>, p.Q703K<sup>(1) (2)</sup>,<br/>p.Y859C<sup>(1)</sup>,c.3005+25C>T<sup>(1)</sup>,<br/>p.T348M<sup>(1)</sup>, p.L411L<sup>(3)</sup>, p.S434S<sup>(1) (2)</sup>,<br/>p.T436P<sup>(3)</sup></td><td>NLRP3 (N = 1)</td><td>p.T348M artificial<br/>deletion (allele<br/>specific primer)</td><td></td><td></td></tr><tr><td></td><td>TNFRSF1A (N = 7)</td><td>p.C30R<sup>(1)</sup>, p.C30F<sup>(1)</sup>, p.C73W<sup>(1)</sup>,<br/>p.H69fs<sup>(1)</sup>, p.E54E<sup>(1)</sup>, p.P46L<sup>(1)</sup>,<br/>p.R92Q<sup>(1)</sup></td><td>ITGB5 (N = 1)</td><td>Chr3q21.2dup</td><td></td><td></td></tr><tr><th colspan="6">Study group</th></tr><tr><td>Negative Controls</td><td>N = 1</td><td>Wild-type DNA</td><td>N = 1</td><td>Asymptomatic individual</td><td>N = 1</td><td>Asymptomatic individual</td></tr><tr><td>Positive Controls</td><td>N = 1</td><td>DNA with a known variation per exon tested where available</td><td>N = 1</td><td>DMD (Internal control)</td><td>N = 1</td><td>DMD (Internal control)</td></tr></tbody></table>

after normalization with the reference gene *GFAP*. The data (area and height of peak) were exported to Microsoft Excel and inserted into the equation:

$$R = \frac{\text{[target exon height or peak area]/reference exon height or peak area}]_{\text{patient}}}{\text{[target exon height or peak area]/reference exon height or peak area}]_{\text{control}}$$

The *DMD* gene was used as a control of correct gene dosage: the ratio (R) *DMD/B2M* should be equal to two for female samples and one for male samples.

## Results

We first implemented qPCR-HRM as a valid method for screening hereditary recurrent fever genes in our laboratory. Then in our patient cohort, we conducted a simultaneous search for punctual mutations and gene rearrangements using qPCR-HRM in the prospective group and qPCR alone in the retrospective group (as these patients had already been assessed for qualitative alterations), then SQF-PCR in both groups as a confirmatory quantitative approach (Figure 1).

## Validation of qPCR-HRM

To implement a strategy that could be used on a routine basis, we set the method up to focus on mutation hot spots (*MVK*, exons 2–11; *NLRP3*, exon 3; *TNFRSF1A*, exons 2–4). The qPCR-HRM

method could be performed for all regions except *MVK* exon 9 for which the Ct could not be obtained below 30 cycles. For this exon, we therefore conducted separate HRM and qPCR. The validation process included a number of DNA controls to assess the successful detection of a variety of gene alterations.

To validate the mutation analysis (HRM), we used 49 DNA samples carrying known mutations and 48 wild-type controls (Table 2). The sensitivity of the technique in our conditions (93.8%) was comparable to other mutation screening methods [9]. Indeed, all pathological mutations (homozygous or heterozygous) were detected, and only 3 polymorphisms (c.769-7_769-6dupT in *MVK* and p.A242A, p.Y141Y in *NLRP3*) went undetected (not shown).

To validate the gene dosage (qPCR), we used DNA samples from 3 individuals (Table 2): 1) One patient with no recurrent fever in whom was found a large deletion (2,13Mb, Chr12:107813096..109948089, version: NCBI36/hg18) including the *MVK* gene in an unrelated study (micro-array screening for mental retardation); 2) One patient with a p.T348M (c.1043C>T) mutation in *NLRP3* in whom we performed allele specific PCR using a wild-type primer to mimic a deletion; 3) One patient with 3 copies of the gene *ITGB5* (154.9kb, Chr3:125946672..126101558, version NCBI36/hg18), since we had no available sample with a duplication in our genes of interest. As normal controls, 10 DNA samples from asymptomatic individuals were screened. Expected copy numbers were obtained (Figure 2).

**Prospective study: qPCR-HRM has little relevance for routine genetic diagnosis of hereditary recurrent fevers.** To

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

<footnote>(1) heterozogote, (2) homozygote, (3) compound heterozygote.</footnote>

<footnote>*These polymorphisms were not detected by HRM.</footnote>

<footnote>†Copy number variation.</footnote>

<footnote>doi:10.1371/journal.pone.0014096.t002</footnote>
<figure><img src="image_9.png" /><figcaption>**Figure 2.** Results of *MVK*, *NLRP3* and *TNFRSF1A* gene dosage analysis by qPCR-HRM (prospective study) and qPCR (retrospective study). We screened by qPCR-HRM (i) exons 2–8 and 10–11 of the *MVK* gene in 10 controls, 1 DNA sample with a 12q23.2 deletion, and 49 patients, (ii) the exon 3 of the *NLRP3* gene in 9 controls, 10 patients, and the artificial positive p.T348M patient and (iii) exons 2–4 of the *TNFRSF1A* gene in 13 controls and 54 patients. We screened by qPCR (i) exon 9 of the *MVK* gene in 10 controls, 1 DNA sample with a 12q23.2q24.11 deletion, and 12 patients and (ii) exons 2–9 of the *NLRP3* gene in 9 controls and 38 patients.</figcaption></figure>

$$R = \frac{[\text{target exon concentration/reference exon concentration}]_{\text{patient}}}{[\text{target exon concentration/reference exon concentration}]_{\text{control}}}$$

The red bars indicate the range limits of the control ratios, between 0.75 and 1.28.

doi:10.1371/journal.pone.0014096.g002

assess the usefulness of this approach for routine genetic diagnosis, we applied it to our prospective study group consisting of 113 unselected patients referred for genetic diagnosis. All amplicons were sequenced in parallel. We confirmed diagnosis in 6 patients, all of whom bore already known mutations: three associated with MKD, p.[V377I]+[V377I], p.[V377I]+[G226S], and p.[T237S]+[T237S]; one with CAPS, p.T348M; and two with TRAPS, p.C70R and p.R92Q. Representative HRM difference plots are shown in Figure 3. The qualitative step (melting curves) detected all the mutations except the T237S homozygous genotype which was only detected after equimolar mixing with a wild type sample using the HRM method alone (not shown). The quantitative step detected no gene dosage abnormalities (Figure 2). A second method (SQF-PCR) was performed on a subset of these patients and confirmed the absence of copy number variants in this cohort (Figure 4).

**Retrospective study: sequence rearrangements are likely absent in hereditary recurrent fever genes.** To confirm this absence of copy number alterations in our genes of interest and enhance the probability of finding new mutations, we extended our series to a panel of selected patients with recurrent fever (N = 88, Figure 1) with either one *MVK* punctual mutation for MKD, or urticaria and/or neurosensorial symptoms for CAPS. Indeed, FMF had previously been addressed³ and TRAPS does not show a clear cut phenotype. In this retrospective series, punctual mutations in all coding exons and intron boundaries had

already been investigated using classical one-step routine analysis (HRM and/or sequencing). No genetic rearrangements were identified (Figure 2). We also performed a second SQF-PCR method in this retrospective panel which confirmed the absence of intragenic duplications or deletions (Figure 4).

Discussion

The existence of a significant number of patients with convincing symptoms though not carrying point mutations suggests that other genes and/or mutational mechanisms, such as gene dosage alteration, could be involved in the etiology of hereditary recurrent fever syndromes. The aim of this study was therefore to explore the possible existence of copy number variants in hereditary recurrent fever genes and estimate the usefulness of qPCR-HRM in routine diagnosis. Duplications and deletions of the genes involved in MKD, a recessive enzymopathy [10], and TRAPS, thought to be a protein misfolding disorder could theoretically lead to these two diseases [11]. CAPS is due to gain-of-function mutations in NLRP3, a major protein of the inflammasome platform that activates pro-caspase-1 into its active form caspase 1, which in turn promotes the processing of the potent pro-inflammatory cytokine interleukin-1β [12]. Classical gain-of-function rearrangements (e.g. duplications) have never been recorded in CAPS. However, a nonsense mutation (p.Arg554X) was identified in *NLRP3* putatively resulting in a

<figure><img src="image_10.png" /></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 3. Pathological genotypes detected by qPCR-HRM in the prospective study.** Qualitative changes were assessed using the gene scanning module of the LightCycler 480 Software version 1.5. Wild type sequences were used to define baselines. Difference plots revealed disease associated variants.
doi:10.1371/journal.pone.0014096.g003</figcaption></figure>

truncated protein lacking all leucine-rich repeats [13]. A partially
deleted protein could therefore result in a gain-of-function.

To address the issue of possible pathological copy number
variations in hereditary recurrent fevers, we took advantage of the
fact that the HRM mutation screening technique is already
routinely used in our laboratories, and implemented a recently
described strategy that combines point mutation and gene
rearrangement detection in a single assay (qPCR-HRM). This
approach proved useful for the genetic analysis of *MLH1* in Lynch
syndrome [7]. The main advantage of qPCR-HRM is the gain in

terms of cost and time. Another major benefit is the validation of
amplification of the two alleles through the quantitative measure,
which avoids the misinterpretation of primer mismatch of one
allele as a deletion, and validates homozygous genotypes. We
found the sensitivity of qPCR-HRM for qualitative analysis as
compared to sequencing (93.8%) to be similar to that of classical
qualitative methods such as DGGE or dHPLC [9]. The sensitivity
of qPCR-HRM and SQF-PCR for quantitative analysis was
100% for the positive controls. The known limit of qPCR-HRM
is the difficulty in designing primers efficient for both qualitative

<figure><img src="image_12.png" /></figure>
Result of MVK, NLRP3 and TNFRSF1A genes dosage by SQF-PCR

<figure><img src="image_13.png" /><figcaption>Figure 4. Results of MVK, NLRP3 and TNFRSF1A gene dosage analysis by SQF-PCR. We screened by SQF-PCR (i) all coding exons of the MVK gene in 9 controls, 1 DNA sample with a 12q23.2 deletion, and 28 patients (ii) all coding exons of the NLRP3 gene in 9 controls and 86 patients and (iii) all coding exons of the TNFRSF1A gene in 9 controls and 54 patients.</figcaption></figure>

$$R = \frac{[\text{target exon height or peak area}/\text{reference exon height or peak area}]_{\text{patient}}}{[\text{target exon height or peak area}/\text{reference exon height or peak area}]_{\text{control}}}$$

The red bars indicate the range limits of the control ratios, between 0.75 and 1.25.

doi:10.1371/journal.pone.0014096.g004

detection and quantitative dosage. Among the 21 amplicons we designed, only exon 9 of *MVK* could not be assessed by qPCR-HRM.

Having validated this new method being in our laboratories, we then used it to assess our patients in a prospective and retrospective study. We detected no gene dosage alteration in a total of 159 patients suspected of MKD, CAPS or TRAPS. A previous study in patients with FMF using MLPA also found no genetic rearrangements. A recent review demonstrated that copy number variants are more frequent in large genes [14]. The hereditary recurrent genes are in contrast rather small (*MVK*: 23434bp, *NLRP3*: 31054bp, *TNFRSF1A*: 13231bp). To investigate whether the chromosomal regions containing these genes could contain fragile sites, we consulted the Database of Genomics Variants (http://projects.tcag.ca/variation) and the Copy Number Variation project at the Children's Hospital of Philadelphia (http://cnv.chop.edu). We retrieved no copy variant within the regions of the *MVK* and *TNFRSF1A* genes. However, 3 variants have been reported in apparently healthy individuals in the 1q44 region of *NLRP3*, the largest of the 3 genes: 1) an insertion of 147bp in intron 3 of the *NLRP3* gene (position 245657612 on chromosome 1q, version NCBI36/hg18) in one Caucasian individual [15]; 2) a deletion spanning intron 4 to the 5' extremity of exon 5 (chr1:245663794..245664071) in one Asian individual [16]; and 3) complex rearrangements upstream of the *NLRP3* gene in 50 French individuals [17]. This may indicate the localization of

*NLRP3* in or close to a minor fragile site, although we identified no pathological rearrangements in our series.

In conclusion, our previous and present results support that deleterious intragenic rearrangements in *MEFV*, *MVK*, *NLRP3* and *TNFRSF1A* are rare or absent from the mutational spectrum of hereditary recurrent fevers, suggesting that other mutations in seemingly affected patients should be searched for either in other gene regions (e.g. promoter and introns) or in other candidate genes. The use of the combined qPCR-HRM as a first line approach for genetic diagnosis is not relevant for these diseases since only qualitative alterations were identified. However, rare pathological gene dosage alterations could still be found in future series and the method we have developed could be used punctually as a second approach since it proved effective and quick.

Acknowledgments

We thank the micro-array platform team for the positive samples with a deletion in Chr12q24.1 and duplication in Chr3q21.2. We are grateful to Dr S. Grandemange for her helpful reviewing of the manuscript, and to Angloscribe for editing the English language.

Author Contributions

Conceived and designed the experiments: NPR IP CR SS IT. Performed the experiments: NPR BD AF EC CR SS. Analyzed the data: NPR LP LC IT. Contributed reagents/materials/analysis tools: NPR LC. Wrote the paper: NPR LP BD IT. Reviewed the paper: LC.

<figure><img src="image_14.png" /></figure>
References

1. Simon A, van der Meer JW (2007) Pathogenesis of familial periodic fever syndromes or hereditary autoinflammatory syndromes. Am J Physiol Regul Integr Comp Physiol 292: R86-98.

2. Masters SL, Simon A, Aksentijevich I, Kastner DL (2009) Horror autoinflammatory: the molecular pathophysiology of autoinflammatory disease. Annu Rev Immunol 27: 621-668.

3. van Gijn ME, Soler S, de la Chapelle C, Mulder M, Ritorre C, et al. (2008) Search for copy number alterations in the MEFV gene using multiplex ligation probe amplification, experience from three diagnostic centres. Eur J Hum Genet 16: 1404-1406.

4. Giansily-Blaizot M, Thorel D, Khau Van Kien P, Behar C, Romey MC, et al. (2007) Characterisation of a large complex intragenic re-arrangement in the FVII gene (F7) avoiding misdiagnosis in inherited factor VII deficiency. Br J Haematol 138: 359-365.

5. Sellner LN, Taylor GR (2004) MLPA and MAPH: new techniques for detection of gene deletions. Hum Mutat 23: 413-419.

6. Laccone F, Junemann I, Whatley S, Morgan R, Butler R, et al. (2004) Large deletions of the MECP2 gene detected by gene dosage analysis in patients with Rett syndrome. Hum Mutat 23: 234-244.

7. Rouleau E, Lefol C, Bourdon V, Coulet F, Noguchi T, et al. (2009) Quantitative PCR high-resolution melting (qPCR-HRM) curve analysis, a new approach to simultaneously screen point mutations and large rearrangements: application to MLH1 germline mutations in Lynch syndrome. Hum Mutat 30: 867-875.

8. Steger G (1994) Thermal denaturation of double-stranded nucleic acids: prediction of temperatures critical for gradient gel electrophoresis and polymerase chain reaction. Nucleic Acids Res 22: 2760-2768.

9. Chou LS, Lyon E, Wittwer CT (2005) A comparison of high-resolution melting analysis with denaturing high-performance liquid chromatography for mutation

scanning: cystic fibrosis transmembrane conductance regulator gene as a model. Am J Clin Pathol 124: 330-338.

10. Cuisset L, Drenth JP, Simon A, Vincent MF, van der Velde Visser S, et al. (2001) Molecular analysis of MVK mutations and enzymatic activity in hyper-IgD and periodic fever syndrome. Eur J Hum Genet 9: 260-266.

11. McDermott MF, Aksentijevich I, Galon J, McDermott EM, Ogunkolade BW, et al. (1999) Germline mutations in the extracellular domains of the 55 kDa TNF receptor, TNFR1, define a family of dominantly inherited autoinflammatory syndromes. Cell 97: 133-144.

12. Hoffman HM, Mueller JL, Broide DH, Wanderer AA, Kolodner RD (2001) Mutation of a new gene encoding a putative pyrin-like protein causes familial cold autoinflammatory syndrome and Muckle-Wells syndrome. Nat Genet 29: 301-305.

13. Jeru I, Hayrapetyan H, Duquesnoy P, Sarkisian T, Amselem S (2006) PYPAF1 nonsense mutation in a patient with an unusual autoinflammatory syndrome: role of PYPAF1 in inflammation. Arthritis Rheum 54: 508-514.

14. Kalari KR, Casavant TL, Scheetz TE (2008) A knowledge-based approach to predict intragenic deletions or duplications. Bioinformatics 24: 1975-1979.

15. Levy S, Sutton G, Ng PC, Feuk L, Halpern AL, et al. (2007) The diploid genome sequence of an individual human. PLoS Biol 5: e254.

16. Wang J, Wang W, Li R, Li Y, Tian G, et al. (2008) The diploid genome sequence of an Asian individual. Nature 456: 60-65.

17. de Smith AJ, Tsalenko A, Sampas N, Scheffer A, Yamada NA, et al. (2007) Array CGH analysis of copy number variation identifies 1284 new genes variant in healthy white males: implications for association studies of complex diseases. Hum Mol Genet 16: 2783-2794.

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