# Factors Predisposing the Response to Lumacaftor/Ivacaftor in People with Cystic Fibrosis

Julie Mésinèle, Manon Ruffin, Loïc Guillot, Pierre-Yves Boëlle, Harriet Corvol

► To cite this version:

Julie Mésinèle, Manon Ruffin, Loïc Guillot, Pierre-Yves Boëlle, Harriet Corvol. Factors Predisposing the Response to Lumacaftor/Ivacaftor in People with Cystic Fibrosis. Journal of Personalized Medicine, 2022, 12 (2), pp.252. 10.3390/jpm12020252. hal-03829149

HAL Id: hal-03829149

https://hal.science/hal-03829149

Submitted on 25 Jan 2023

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
<figure><img src="image_2.png" /></figure>

Article

# Factors Predisposing the Response to Lumacaftor/Ivacaftor in People with Cystic Fibrosis

Julie Mésinèle<fnref n="1" />,<fnref n="2" /> ID, Manon Ruffin<fnref n="1" /> ID, Loïc Guillot<fnref n="1" />,*, ID, Pierre-Yves Boëlle<fnref n="2" />,† ID, Harriet Corvol<fnref n="3" />,*,† ID
and on behalf of the French CF Modifier Gene Study Investigators ‡

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

check for
updates

**Citation:** Mésinèle, J.; Ruffin, M.; Guillot, L.; Boëlle, P.-Y.; Corvol, H.; on behalf of the French CF Modifier Gene Study Investigators. Factors Predisposing the Response to Lumacaftor/Ivacaftor in People with Cystic Fibrosis. *J. Pers. Med.* **2022**, *12*, 252. https://doi.org/10.3390/jpm12020252

Academic Editors: Miquéias Lopes-Pacheco, Neeraj Sharma and Shafagh Waters

Received: 16 December 2021

Accepted: 7 February 2022

Published: 10 February 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Lumacaftor/ivacaftor (LUMA-IVA) therapy is prescribed to people with cystic fibrosis (pwCF) homozygous for the Phe508del-CFTR variant to restore CFTR protein function. There is, however, large inter-individual variability in treatment response. Here, we seek to identify clinical and/or genetic factors that may modulate the response to this CFTR modulator therapy. A total of 765 pwCF older than 12 years under LUMA-IVA therapy and with lung function and nutritional measurements available before and after treatment initiation were included. Response to treatment was determined by the change in lung function and nutritional status, from baseline and over the first two years after initiation, and it was assessed by weighted generalized estimating equation models. Gains in lung function and nutritional status were observed after 6 months of treatment (on average 2.11 ± 7.81% for percent predicted FEV₁ and 0.44 ± 0.77 kg/m<fnref n="2" /> for BMI) and sustained over the 2 years. We observed that the more severe patients gained the most in lung function and nutritional status. While females started with a nutritional status more impaired than males, they had a larger response and regained BMI Z-score values similar to men after 2 years of treatment. We observed no association between variants in solute carrier (SLC) genes and the respiratory function response to LUMA-IVA, but the SLC6A14 rs12839137 variant was associated with the nutritional response. Further investigations, including other genomic regions, will be needed to fully explore the inter-individual variability of the response to LUMA-IVA.

**Keywords:** CFTR modulator therapy; long-term effect; lung function; nutritional status; modifiers genes; solute carrier family genes; cystic fibrosis

## 1. Introduction

Cystic fibrosis (CF) is a rare, autosomal recessive, life-shortening genetic disease that affects more than 90,000 people worldwide [1]. It is caused by variants in the gene encoding the cystic fibrosis transmembrane conductance regulator (CFTR), a chloride channel expressed in epithelial cells throughout the body. Until one decade ago, only symptomatic treatments were available for people with CF (pwCF). In the last decade, considerable efforts have led to the development of therapies that target the CFTR protein named CFTR modulators. Since 2012, pwCF carrying some *CFTR* gating variants can be treated with ivacaftor, a potentiator therapy which increases the probability of CFTR-channel opening.

Following this, lumacaftor, a CFTR corrector that improves the processing and trafficking of the Phe508del-CFTR protein, was combined with ivacaftor to treat pwCF ho-

<footnote n="1">¹ Centre de Recherche Saint-Antoine (CRSA), Inserm, Sorbonne Université, 75012 Paris, France; julie.mesinele@inserm.fr (J.M.); manon.ruffin@inserm.fr (M.R.)</footnote>

<footnote n="2">² Hôpital Saint-Antoine, AP-HP, Institut Pierre Louis d'Epidémiologie et de Santé Publique (iPLESP), Inserm, Sorbonne Université, 75012 Paris, France; pierre-yves.boelle@sorbonne-universite.fr</footnote>

<footnote n="3">³ Service de Pneumologie Pédiatrique, Hôpital Trousseau, AP-HP, Centre de Recherche Saint-Antoine (CRSA), Inserm, Sorbonne Université, 75012 Paris, France</footnote>

<footnote>* Correspondence: loic.guillot@inserm.fr (L.G.); harriet.corvol@aphp.fr (H.C.)</footnote>

<footnote>† These authors contributed equally to this work.</footnote>

<footnote>‡ Membership of the French CF Modifier Gene Study Investigators is provided in the Acknowledgments.</footnote>
mozygous for the Phe508del-CFTR variant. A phase 2 trial of this association, lumacaftor-ivacaftor (LUMA-IVA), demonstrated improvements of the lung function, measured by an increase of the percent-predicted forced expiratory volume in one second (ppFEV₁), and of the nutritional status, evaluated by the body mass index (BMI) [2]. Marketing authorizations were granted in 2015.

Since then, clinical benefits of LUMA-IVA have been questioned, particularly because of its high cost [3,4]. In addition, phase 3 trials and real-life studies over the first 2 years of treatment demonstrated highly heterogeneous responses in lung function and nutritional status [5–11]. These studies also highlighted the inter-individual variability of the airway response and the limited tolerance to treatment, with high discontinuation rates ranging from 17.2% to 28.9% [5,7,9]. Predisposing factors for interruption and for the response variability were shown to be baseline lung function, age at treatment initiation and gender [5,7–11].

In a previous study, we found that lung response variability to ivacaftor was associated with variants in the *Solute Carrier Family 26 Member 9 (SLC26A9)* gene in a cohort of French pwCF, confirming results observed in Canadians [12,13]. However, these results were controverted in North American pwCF [14]. Nevertheless, this gene was also shown to be involved in lung function variability [13,15] and to meconium ileus susceptibility in pwCF [16]. Furthermore, other variants in the *SLC* gene family (*SLC6A14* and *SLC9A3*) have been shown to modify several CF phenotypes, suggesting a pleiotropic effect of this gene family of genes [15–18].

Here, we analysed the evolution of pwCF treated with LUMA-IVA and included the "French CF Gene Modifier Study" to identify clinical and/or genetic factors that may change the response to this CFTR modulator therapy.

## 2. Materials and Methods

### 2.1. *Study Design and Patients*

The "French CF modifier gene study" is a nationwide observational study of pwCF in France with prospective data collection from participants. Here, we aimed to evaluate the lung and nutritional response to LUMA-IVA over the first two years of treatment. As of 31 January 2021, 4975 pwCF attending one of the 47 French CF centres were included in this national cohort (corresponding to ~70% of all French pwCF) [19].

Since 2001, neonatal CF screening has been underway in France and specialized care for pwCF has been delivered in expert CF centres according to national recommendations [20,21]. Standardized longitudinal data have been prospectively collected for each pwCF, making the current analysis possible. Data collected in electronic or paper medical records over time in the participating CF centres were input in a national database.

All pwCF older than 12 years, homozygous for the Phe508del-CFTR variant and treated by LUMA-IVA out of clinical trials were eligible for inclusion, provided they had at least a 6-month follow-up after treatment initiation and had not received a lung or liver transplant (n = 878) (Figure 1). We excluded pwCF without FEV₁ and/or BMI values at baseline (n = 22) and/or over the first 6 months after treatment initiation due to missing data (n = 42) or early treatment interruption (n = 49). Finally, 765 pwCF were included in the analyses (Figure 1).

### 2.2. *Lung and Nutritional Response to Combined Lumacaftor/IVacaftor Therapy*

According to CF care recommendations, lung function testing including spirometry and anthropometric measurements are performed at every visit in pwCF [21,22]. To assess the lung function, measurements of FEV₁ were expressed as percent-predicted values using the Global Lung Function Initiative equations [23] or transformed to the Kulich Normalized Mortality Adjusted CF-specific lung phenotype (SaKnorm Z-value) [24,25]. SaKnorm is a quantitative phenotype that allows the direct comparison of lung phenotypes between pwCF and accounts for differential survival. To evaluate the nutritional response, body mass index (BMI) measurements were Z-score transformed according to WHO Child Growth Standards [26]. For each pwCF, SaKnorm Z-value and BMI Z-score baseline values
were computed as the measurements' average over the 6 months prior to the treatment initiation.

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

FEV₁: Forced expiratory volume in one second

BMI: Body mass index

<caption>**Figure 1.** Flowchart of patient selection. Flowchart describing the selection of people with cystic fibrosis included in the analysis of the lung and nutritional response to combined lumacaftor/ivacaftor therapy.</caption>

We classified the baseline values according to terciles as follows: (i) SaKnorm Z-values from -1.620 to 0.168 (1st tercile), from 0.169 to 0.791 (2nd tercile) and from 0.792 to 2.520 (3rd tercile) for lung response analysis, (ii) BMI Z-score from -4.330 to -1.080 (1st tercile), from -1.079 to -0.307 (2nd tercile) and from -0.306 to 2.890 (3rd tercile) for nutritional response analysis.

Gender, age at treatment initiation, meconium ileus, CF-related diabetes (CFRD), CF-related liver disease (CFLD) and presence of Pseudomonas aeruginosa chronic colonization (*Pa*-CC) were analysed (see Supplementary Information for CFRD, CFLD and *Pa*-CC definitions). Age at treatment initiation was categorized in 3 classes: "≤20" years old, "20–30" years old, and ">30" years old. Changes in FEV₁ and BMI courses over time according to age at initiation was assessed relative to the "≤20" years old level, indicated as the "reference" level.

## 2.3. Modifiers of Response to Combined Lumacaftor/IVacaftor Therapy

We analysed single nucleotide polymorphisms (SNPs) located in or near SLC genes that had already been shown to be associated with lung response variability to ivacaftor [12,13] and/or several CF phenotypes [15–18]. Using Kompetitive Allele Specific PCR (LGC Group, Teddington, UK), we genotyped 9 SNPs of the following genes: SLC26A9 (rs7512462, rs1874361, rs4077468, rs4077469, rs7419153, rs12047830), SLC9A3 (rs57221529) and SLC6A14 (rs3788766 and rs12839137).

## 2.4. Statistical Analysis

Descriptive statistics were reported as mean ± standard deviation (SD) or ± standard deviation of mean (SDM) or percentages and 95% confidence intervals (95% CIs), as appropriate. The lung function and nutritional response over the first 2 years of LUMA-
IVA therapy were analysed by generalized estimating equation mixed models adapted to longitudinal data. Measurements of SaKnorm Z-value and BMI Z-score were averaged over successive 6-month periods, starting 6 months before treatment initiation (baseline value) and extending over the 2 years following treatment initiation. To account for missing data in patients' follow-ups, we used generalized estimating equations with weights according to the probability of missingness [27,28]. The missingness model included 4 variables: time since treatment initiation, age at treatment initiation, change in lung function in the previous semester and the baseline lung function value. These criteria included factors associated with treatment discontinuation [7].

For the genetic association study, we applied additive SNP coding, and reference alleles were determined as those with the highest frequency in the European population (http://www.ensembl.org (accessed on 1 June 2021)). SNPs in the chromosome X (SLC6A14 rs3788766 and rs12839137) were also additively coded (0 or 1 or 2 for women and 0 or 2 for men). Fisher's exact test was used to test the conformance of the allele frequencies with the Hardy-Weinberg equilibrium.

All analyses were carried out using the R software (version 4.1.0, http://www.R-project.org (accessed on 1 June 2021)) using the package "wgeesel".

## 3. Results

### 3.1. Study Population

Demographics and baseline characteristics of the eligible pwCF on LUMA-IVA are summarized in Table 1. Distribution is similar across the analysed and the excluded groups by gender, origin, clinical characteristics, lung disease and nutritional severity. Compared to analysed pwCF, excluded pwCF were older (p-value < 0.001). Among the 765 pwCF analysed, 14.5% (111 pwCF) discontinued treatment.

<table><caption>**Table 1.** Demographics and baseline (prior to combined lumacaftor/ivacaftor therapy) characteristics of the 878 eligible people with cystic fibrosis.</caption><thead><tr><th></th><th>Patients<br>Analysed<br>n = 765</th><th>Patients<br>Excluded<br>n = 113</th><th>p-Value</th></tr></thead><tbody><tr><td>Male, % (n)</td><td>56.21% (430)</td><td>46.02% (52)</td><td>0.042</td></tr><tr><td>Caucasian origin, % (n)</td><td>97.51% (745)</td><td>97.35% (110)</td><td>0.915</td></tr><tr><td>Age at treatment initiation (years), mean ± SD</td><td>22.2 ± 9.0</td><td>25.2 ± 9.4</td><td>&lt;0.001</td></tr><tr><td>Age at treatment initiation (years), % (n)</td><td></td><td></td><td></td></tr><tr><td>≤20</td><td>50.07% (383)</td><td>31.86% (36)</td><td>0.001</td></tr><tr><td>20–30</td><td>28.37% (217)</td><td>38.05% (43)</td><td></td></tr><tr><td>&gt;30</td><td>21.57% (165)</td><td>30.09% (34)</td><td></td></tr><tr><td>Presence of meconium ileus % (n)</td><td>16.57% (118)</td><td>25.47% (27)</td><td>0.025</td></tr><tr><td>Presence of CFRD *, % (n)</td><td>25.93% (195)</td><td>31.13% (33)</td><td>0.256</td></tr><tr><td>Presence of CFLD *, % (n)</td><td>36.04% (275)</td><td>45.54% (51)</td><td>0.052</td></tr><tr><td>Presence of Pa-CC *, % (n)</td><td>38.56% (295)</td><td>44.25% (50)</td><td>0.248</td></tr><tr><td>Lung disease severity ‡ (SaKnorm Z-value),<br>mean ± SD</td><td>0.447 ± 0.702</td><td>0.457 ± 0.710</td><td>0.942</td></tr><tr><td>Nutritional severity † (BMI Z-score), mean ±<br>SD</td><td>-0.677 ± 0.932</td><td>-0.807 ± 1.060</td><td>0.072</td></tr></tbody></table>

### 3.2. Overall Changes in Lung Function and Nutritional Status with LUMA-IVA

Following the initiation of LUMA-IVA, the lung function and nutritional status improved with gains in SaKnorm Z-value of 0.106 ± 0.015 (p-value < 0.0001) and in BMI Z-score of 0.108 ± 0.017 (p-value < 0.0001). Improvements were sustained over the 2 years

<footnote>* Before treatment initiation; ‡ Over the past 3 years before treatment initiation, forced expiratory volume in one second (FEV₁) measurements were expressed as Kulich Normalized Mortality Adjusted CF-specific lung phenotype (SaKnorm Z-value) (1,2); † Over the past 3 years before treatment initiation, body mass index (BMI) measurements were Z-score transformed according to WHO Child Growth Standards (3). Among the 878 eligible patients, 1 had missing data for ethnicity, 60 for meconium ileus, 20 for CFRD and 3 for CFLD. Abbreviations: CFRD: cystic fibrosis-related diabetes, CFLD: cystic fibrosis liver disease, Pa-CC: *Pseudomonas aeruginosa* chronic colonization.</footnote>
of treatment with change slopes of 0.054 ± 0.010 (p-value < 0.0001) for the lung function and of 0.057 ± 0.015 (p-value = 0.0001) for the nutritional response. In the first semester, the average absolute gains from baseline were 2.11 ± 0.14% predicted for ppFEV₁ and 0.44 ± 0.01 kg/m² for BMI (Table S1).

Treatment response was similar according to age at initiation, presence of meconium ileus, CFRD, CFLD or *Pa*-CC status, but changed by pre-treatment values. Indeed, pwCF with the lowest SaKnorm and BMI baseline values had the greater gains in lung function and nutritional status (Table 2 and Figure 2).

<table><caption>**Table 2.** Lung function and nutritional response of combined lumacaftor/ivacaftor therapy, according to clinical and demographic characteristics in 765 people with cystic fibrosis.</caption><thead><tr><th></th><th>Change in SaKnorm(Z-Value) ‡ ± SD</th><th>p-Value</th><th>Change in BMI(Z-Score) ‡ ± SD</th><th>p-Value</th></tr></thead><tbody><tr><td>Female</td><td>-0.015 ± 0.023</td><td>0.5203</td><td>0.069 ± 0.034</td><td>0.0415</td></tr><tr><td>Age at initiation (Years)</td><td></td><td></td><td></td><td></td></tr><tr><td>≤20</td><td>Reference</td><td>Reference</td><td>Reference</td><td>Reference</td></tr><tr><td>20–30</td><td>0.011 ± 0.026</td><td>0.6748</td><td>0.006 ± 0.037</td><td>0.8705</td></tr><tr><td>>30</td><td>-0.068 ± 0.023</td><td>0.0029</td><td>-0.048 ± 0.038</td><td>0.2047</td></tr><tr><td>Presence of meconium ileus</td><td>-0.019 ± 0.031</td><td>0.5401</td><td>0.010 ± 0.048</td><td>0.8407</td></tr><tr><td>Presence of CFRD *</td><td>-0.003 ± 0.023</td><td>0.9066</td><td>-0.011 ± 0.039</td><td>0.7698</td></tr><tr><td>Presence of CFLD *</td><td>-0.004 ± 0.024</td><td>0.8839</td><td>0.003 ± 0.036</td><td>0.9414</td></tr><tr><td>Presence of Pa-CC *</td><td>0.018 ± 0.023</td><td>0.4380</td><td>-0.022 ± 0.034</td><td>0.5116</td></tr><tr><td>Baseline</td><td></td><td></td><td></td><td></td></tr><tr><td>3rd tercile</td><td>Reference</td><td>Reference</td><td>Reference</td><td>Reference</td></tr><tr><td>2nd tercile</td><td>0.105 ± 0.027</td><td>0.0001</td><td>0.080 ± 0.036</td><td>0.0273</td></tr><tr><td>1st tercile</td><td>0.145 ± 0.028</td><td>&lt;0.0001</td><td>0.237 ± 0.040</td><td>&lt;0.0001</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Change according to baseline terciles following initiation of combined lumacaftor/ivacaftor therapy with 95% CIs in (A) lung function (SaKnorm Z-value), and (B) nutritional status (BMI Z-score), in 765 people with cystic fibrosis (pwCF). For lung response analysis, the baseline SaKnorm (Z-value) terciles were from -1.620 to 0.168 (1st tercile in red), from 0.169 to 0.791 (2nd tercile in orange) and from 0.792 to 2.520 (3rd tercile in green). For nutritional response analysis, the baseline BMI (Z-score) terciles were from -4.330 to -1.080 (1st tercile in red), from -1.079 to -0.307 (2nd tercile in orange) and from -0.306 to 2.890 (3rd tercile in green). The black dotted lines show the model predicted change and the grey dotted lines the overall change in the 765 pwCF.</figcaption></figure>

<footnote>‡ Forced expiratory volume in one second (FEV₁) measurements were expressed as Kulich Normalized Mortality Adjusted CF-specific lung phenotype (SaKnorm Z-value) (1,2); † Body mass index (BMI) measurements were Z-score transformed according to WHO Child Growth Standards (3); * Before initiation of lumacaftor/ivacaftor combination therapy. Abbreviations: CFRD: cystic fibrosis-related diabetes, CFLD: cystic fibrosis liver disease, Pa-CC: *Pseudomonas aeruginosa* chronic colonization.</footnote>
At 6 months post-treatment, the pwCF in the baseline 1st tercile had an increased ppFEV₁ of 2.68 ± 0.25% predicted and of 0.49 ± 0.02 kg/m² for BMI (Table S1). Gender was not associated with lung function response, but females had a higher nutritional response (p-value = 0.04) (Table 2 and Figure 3A). At initiation of LUMA-IVA, the baseline nutritional status of women (−0.805 ± 0.992 of BMI Z-score) was more impaired than men (−0.642 ± 0.931 of BMI Z-score). After 2 years of LUMA-IVA, both men and women had a similar nutritional status (−0.511 ± 0.931 and −0.522 ± 0.922 for BMI Z-score for men and women, respectively) (Figure 3B).

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Nutritional status (BMI Z-score) over 2 years following initiation of combined lumacaftor/ivacaftor therapy according to gender with 95% CIs, (A) mean of absolute changes from baseline, (B) mean of BMI Z-score, in 765 people with cystic fibrosis (pwCF). Men are in blue, women in orange and the dotted lines show the overall change in the 765 pwCF.</figcaption></figure>

### 3.3. Genetic Analysis

Genotype distributions and results of weighted generalized estimating equation mixed models are shown in Table 3. There was no evidence of association between lung function response to LUMA-IVA and variants of SLC26A9, SLC9A3 or SLC6A14 genes.

<table><caption>**Table 3.** Lung function and nutritional response of combined lumacaftor/ivacaftor therapy, according to SLC26A9, SLC9A3 and SLC6A14 variants, in 765 people with cystic fibrosis.</caption><thead><tr><th></th><th>Position <sup>£</sup></th><th>Alleles *</th><th>MAF EUR</th><th>MAF Cohort</th><th>HWE **</th><th>Change in SaKnorm Z-Value ‡ ± SD</th><th>p-Value</th><th>Change in BMI Z-Score ‡ ± SD</th><th>p-Value</th></tr></thead><tbody><tr><td>SLC26A9</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>rs1874361</td><td>1:205939058</td><td>A/C</td><td>0.48</td><td>0.46</td><td>0.419</td><td>-0.006 ± 0.017</td><td>0.7338</td><td>0.009 ± 0.022</td><td>0.6919</td></tr><tr><td>rs4077468</td><td>1:205945629</td><td>A/G</td><td>0.41</td><td>0.41</td><td>0.406</td><td>-0.011 ± 0.016</td><td>0.4976</td><td>-0.019 ± 0.023</td><td>0.4178</td></tr><tr><td>rs4077469</td><td>1:205945757</td><td>C/T</td><td>0.41</td><td>0.41</td><td>0.496</td><td>-0.013 ± 0.016</td><td>0.4177</td><td>-0.023 ± 0.023</td><td>0.3350</td></tr><tr><td>rs7419153</td><td>1:205948181</td><td>G/A</td><td>0.38</td><td>0.41</td><td>0.033</td><td>0.011 ± 0.017</td><td>0.4894</td><td>0.015 ± 0.024</td><td>0.5456</td></tr><tr><td>rs7512462</td><td>1:205930467</td><td>T/C</td><td>0.41</td><td>0.41</td><td>0.451</td><td>-0.006 ± 0.017</td><td>0.7102</td><td>-0.023 ± 0.023</td><td>0.3166</td></tr><tr><td>rs12047830</td><td>1:205947571</td><td>G/A</td><td>0.49</td><td>0.47</td><td>0.501</td><td>-0.011 ± 0.017</td><td>0.5290</td><td>-0.017 ± 0.024</td><td>0.4936</td></tr><tr><td>SLC9A3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>rs57221529</td><td>5:586509509</td><td>A/G</td><td>0.21</td><td>0.20</td><td>0.735</td><td>0.029 ± 0.019</td><td>0.1256</td><td>0.022 ± 0.029</td><td>0.4564</td></tr><tr><td>SLC6A14</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>rs3788766</td><td>X:116435671435671</td><td>G/A</td><td>0.36</td><td>0.37</td><td>0.388</td><td>0.000 ± 0.014</td><td>0.9719</td><td>0.008 ± 0.019</td><td>0.6797</td></tr><tr><td>rs12839137</td><td>X:116434382</td><td>G/A</td><td>0.22</td><td>0.21</td><td>0.862</td><td>0.000 ± 0.017</td><td>0.9852</td><td>0.045 ± 0.020</td><td>0.0276</td></tr></tbody></table>

<footnote>£ Physical position according to Ensembl GRCh38 (www.ensembl.org (accessed on 1 June 2021)); * Minor allele in European Population are in bold; ** Hardy-Weinberg equilibrium (HWE) p-values were computed by Fisher's exact test and among women for SLC6A14; ‡ Forced expiratory volume in one second (FEV1) measurements were expressed as Kulich Normalized Mortality Adjusted CF-specific lung phenotype (SaKnorm Z-value) (1,2); † Body mass index (BMI) measurements were Z-score transformed according to WHO Child Growth Standards (1). MAF: minor allele frequency.</footnote>
For *SLC26A9* rs7512462 variant and for each additional C-allele, the model estimated a decreased gain in SaKnorm Z-value of -0.006 (±0.017) (p-value = 0.71) and in BMI Z-score of -0.023 (±0.023) (p-value = 0.32) (Table 3 and Figure 2A). In the case of TT, TC and CC genotypes, over the 6 months post-treatment, the predicted average ppFEV₁ absolute gains from baseline were 1.95 ± 0.25%, 2.13 ± 0.20% and 2.40 ± 0.32%, respectively (Table S1).

We observed that only the *SLC6A14* rs12839137 variant was associated with the nutritional response to LUMA-IVA (Table 3). In the first semester, the absolute BMI gains were 0.38 kg/m² (±0.02) and 0.63 kg/m² (±0.03) for the homozygous genotypes GG and AA, respectively, and of 0.47 kg/m² (±0.04) for the GA heterozygous genotypes (Table S1).

**4. Discussion**

LUMA-IVA was the first promising targeted therapy available for the management of pwCF homozygous for the Phe508del-CFTR variant, the most common CF genotype worldwide. Its use in real life has shown a heterogeneous and relatively modest improvement in lung function with limited tolerance. This study analysed clinical factors as well as *solute carrier* variants, located on *SLC26A9*, *SLC9A3* and *SLC6A14* genes, as predictors of lung function and nutritional response to LUMA-IVA.

At 6 months post-treatment, we observed gains of lung function and nutritional status very close to findings of clinical trials that led to the marketing of the therapy [2]. Improvement in lung function severity was sustained over the first 2 years of treatment. The positive slope (+0.054 ± 0.010 of SaKnorm Z-value by year) was reflected in a reduction of ppFEV₁ decline over the 2 years post-treatment (Table S2), again observed in clinical trials [6]. Among the 878 pwCF who started LUMA-IVA, 19.7% (173 pwCF) interrupted the treatment, similarly to a previously reported rate of discontinuation of 17.2% due to adverse effects [5].

Since 2016 in France, LUMA-IVA treatment can be prescribed to pwCF homozygous for the Phe508del-CFTR variant without baseline ppFEV₁ requirements. We observed that pre-treatment lung function and nutritional status were predictors of treatment response. Indeed, as the severity of pwCF decreased, so did the gain in lung function. The absolute mean change in ppFEV₁ after 6 months of treatment was +2.68% in the more severe patients (1st SaKnorm Z-value tercile), and +1.89% and +1.77% for patients in the 2nd and 3rd terciles of SaKnorm Z-value, respectively. This pattern was observed in the efficacy analysis, in which improvement in ppFEV₁ after 6 months of treatment was +3.30% (95% CI: 0.20 to 6.40) in pwCF with baseline ppFEV₁ < 40% and +2.80% (95% CI: 1.70 to 3.80) in pwCF with baseline ppFEV₁ ≥ 40% [11]. However, it was also reported that, after 6 months of LUMA-IVA, pwCF with baseline ppFEV₁ < 40% had a non-significant loss of lung function of -0.40% (-1.90 to 1.10) [8], while those with baseline ppFEV₁ ≥ 90% had an unchanged lung function of 0.15% [10]. The differences in lung function improvement in these reports and our findings may be due to airways' response assessment. We used ppFEV₁ transformed in SaKnorm Z-value to classify pwCF at baseline. This adjustment is an indicator of lung function severity and provides for the direct comparison of pwCF of various ages. We used this specific phenotype to account for age-dependent clinical variables (CFRD, CFLD and Pa-CC). In addition, it should be noted that previous studies included fewer than 50 pwCF, whereas the current study analysed 765 pwCF of the "French CF modifier gene study", which includes more than 85% of all French adult pwCF [29].

Apart from baseline lung function, no demographic or clinical factors were found to be associated with the respiratory response to treatment. Only gender was associated with the nutritional response, with a statistically significant greater change in BMI Z-score in women after LUMA-IVA initiation. This result is in part related to a difference in nutritional severity between men and women before treatment initiation. Indeed, we observed that women gained the same level of nutritional status as men after 2 years of treatment. All studies agree that there is a significant gain in BMI, ranging from 0.21 to 0.96 kg/m², irrespective of the LUMA-IVA dose, lung function severity or treatment duration [2,6–8,11].
However, none of these studies showed differences in gender for the nutritional status response. So far, female gender was only associated with a higher discontinuation rate [5].

This study, with a very large cohort of pwCF homozygous for the Phe508del-CFTR variant, did not show any evidence of association between SLC genes' family and respiratory response to LUMA-IVA therapy. The SLC26A9 gene was suggested to modulate the airway response to ivacaftor in pwCF who harbour at least one of the targeted CFTR gating variants [12–14]. Moreover, ex vivo experiments with nasal epithelial cells from pwCF showed that the SLC26A9 rs7512462 variant was associated with CFTR function in response to LUMA-IVA [30]. However, we observed no association of this variant with either the respiratory or the nutritional response to LUMA-IVA. Among the variants analysed, only the SLC6A14 rs12839137 variant was associated with the nutritional response. Interestingly, two studies in Finnish and French showed that several SLC6A14 variants were associated with obesity [31,32]. As the amino acid transporter SLC6A14 could modulate tryptophan availability for serotonin synthesis, authors hypothesized that these variants might affect perception of hunger and satiety. Miranda et al. also showed that other SLC6A14 variants were associated with food intake in children at 7–8 years of age [33]. Finally, an in vitro study highlighted that the obesity-associated rs2011162 variant, located in the 3'UTR of the SLC6A14 gene, reduces the expression of the transporter and that Slc6a14⁻/⁻ mice develop obesity, fatty liver and metabolic syndrome when fed with a high-fat diet [34]. These findings support the hypothesis of the involvement of SLC6A14 in the regulation of the nutritional status in pwCF.

Our study has several limitations. First, missing information led to the exclusion of pwCF from the analysis of FEV₁ and BMI change. According to French CF care recommendations, lung function and anthropometric measurements are measured every quarter in pwCF [21], and, following recommendations of the French CF Learning Society, each pwCF who initiates LUMA-IVA has systematic visits at treatment initiation and at 1, 3, 6 and 12 months after initiation (with clinical assessment and a pulmonary function test) [9]. Missing data could be due to the collection date from patients' paper and electronic medical records. However, we did not observe any clinical differences between included and excluded patients. Secondly, we observed a 14.5% treatment interruption among the 765 pwCF analysed, of which 10.81% (n = 12) could be attributed to a lack of observed benefits by the patient and/or physician. It was also shown that the risk of discontinuation increased with age and severity of lung function [7]. However, to account for this selection bias and shorter follow-up times, we used weighted statistical models to estimate the probability of missing data, as a function of age at treatment initiation, change in lung function at the prior semester and the baseline lung function value.

**5. Conclusions**

LUMA-IVA shows modest benefits in terms of lung function and nutritional status in pwCF. However, benefits in quality of life should not be overlooked, as well as the reduction of respiratory exacerbations, which have a major role in the CF management. We showed that variants located in the SLC gene family (SLC26A9, SLC9A3 and SLC6A14) were not predictors of the airway treatment response and further investigations, including of other genomic regions, are needed to better explain inter-individual variability. Nevertheless, we observed that the variant rs12839137 of the SLC6A14 gene was associated with the nutritional response, which will be important to further explore. With the development of the next-generation CFTR modulators, such as the elexacaftor/ivacaftor/tezacaftor tritherapy, identifying genetic modifiers involved in the response to treatment is an opportunity to achieve predictive and personalized medicine in CF.
**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jpm1202025/s1, Supplementary Information (phenotype definition, Analysis). Table S1: Mean of lung and nutritional response to combined lumacaftor/ivacaftor therapy within 6 months following initiation, in 765 people with cystic fibrosis. Table S2: Lung function decline in people with cystic fibrosis according to lumacaftor/ivacaftor status. Reference [35] is cited in the Supplementary Materials.

**Author Contributions:** Conceptualization, L.G., M.R., H.C. and P.-Y.B.; methodology, P.-Y.B. and J.M.; writing—original draft preparation, J.M., H.C. and P.-Y.B.; writing—review and editing, J.M., M.R., L.G., H.C. and P.-Y.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Fondation pour la Recherche Médicale (grant number "PMJ201810007620" to Julie Mésinèle), Vaincre la Mucoviscidose, Blanche pour Vaincre la Mucoviscidose, Agir et Informer Contre la Mucoviscidose, Institut National de la Santé et de la Recherche Médicale.

**Institutional Review Board Statement:** The study was conducted in accordance with the Declaration of Helsinki and approved by the French ethical committee (CPP n°2004/15), and the information collection was approved by the Commission Nationale de L'informatique et des Libertés (n°04.404).

**Informed Consent Statement:** Informed consent in writing was obtained from each patient and/or guardian.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ethical restrictions.

**Acknowledgments:** We would like to thank the people with CF, care providers and clinic coordinators at CF Centres, as well as the French CF Modifier Gene Study Investigators, namely: Michel ABELY, CHU de Reims, Reims; Nathalie ALLOU, CHU de la Réunion—Félix Guyon, la Réunion; Baptiste ARNOUAT, CHU de Vannes, Vannes; Carole BAILLY PICCINI, CHU de Nice, Nice; Chantal BELLEGUIC, Hôpital Pontchaillou, Rennes; Katia BESSACI, CHU de Reims, Reims; Tiphaine BI-HOUEE, Hôpital Mère-Enfant, Nantes; Yves BILLON, Hôpital de Brabois, Vandoeuvre les Nancy, Nancy; François Brémont, Hôpital des Enfants de Toulouse, Toulouse; Jacques BROUARD, CHU de Caen, Caen; Stéphanie BUI, CHU de Bordeaux, Bordeaux; Pierre-Régis BURGEL, APHP—Cochin, Paris; Boubou CAMARA, Hôpital de la Tronche, Grenoble; Karine CAMPBELL, CHU de Caen, Caen; Raphael CHIRON, Hôpital Arnaud de Villeneuve, Montpellier; Emmanuelle COIRIER-DUET, Hôpital André Mignot, Le Chesnay; Laure COSSON, Hôpital Clocheville, Tours; Marie-Laure DAL-PHIN, CHU de Besançon, Besançon; Isabelle DANNER BOUCHER, CHU de Nantes, Nantes; Valérie DAVID, Hôpital Mère-Enfant, Nantes; Eric DENEUVILLE, CHU de Rennes—Hôpital sud, Rennes; Philippe DOMBLIDES, CHU de Bordeaux, Bordeaux; Stéphane DOMINIQUE, Centre Hospitalier Universitaire Charles Nicolle, Rouen; Jean-Christophe DUBUS, APHM—Hôpital d'Enfants de la Timone, Marseille; Nadine DUFEU, APHM—Hôpital Nord, Marseille; Isabelle DURIEU, Hospices Civils de Lyon—Groupe hospitalier sud, Lyon; Sandra DURY, Hôpital Maison Blanche, Reims; Ralph EPAUD, CH Intercommunal de Créteil, Créteil; Annlyse FANTON, CHU de Dijon, Dijon; Michael FAYON, CHU de Bordeaux, Bordeaux; Pierre FOUCAUD, Hôpital André Mignot, Le Chesnay; Michèle GERARDIN, APHP—Robert Debré, Paris; Jean-Louis GINIES, CHU d'Angers, Angers; Dominique GRENET, Hôpital Foch, Suresnes; Marie-Christine HERAUD, CH Estaing, Clermont-Ferrand; Véronique HOUDOUIN, APHP—Robert Debré, Paris; Dominique HUBERT, APHP—Cochin, Paris; Frédéric HUET, CHU de Dijon, Dijon; Romain KESSLER, CHU de Strasbourg, Strasbourg; Alice LADAURADE, CHU De Besançon, Besançon; Jeannes LANGUEPIN, CH de Limoges, Limoges; Muriel LAURANS, CHU de Caen, Caen; Pascal Le ROUX, CH Le Havre, Montivilliers; Olivier Le ROUZIC, CHU de Lille, Lille; Sylvie LEROY, CHU de Nice, Nice; Cathie LLERENA, CHU de Grenoble, Grenoble; Julie MACEY, CHU de Bordeaux, Bordeaux; Julie MANKIKIAN, CHU de Tours, Tours; Christophe MARGUET, CHU de Rouen, Rouen; Laurent MELY, Hospices Civils de Lyon—Renée Sabran, Giens; Marie MITTAINE, Hôpital des Enfants de Toulouse, Toulouse; Anne MUNCK, APHP—Robert Debré, Paris; Marlène MURRIS-ESPIN, CHU de Toulouse, Toulouse; Raphaëlle NOVE JOSSERAND, Hospices Civils de Lyon—Groupe hospitalier sud, Lyon; Caroline PERISSON, Groupe Hospitalier Sud Réunion, La Réunion; Isabelle PIN, CHU de Grenoble, Grenoble; Sophie RAMEL, Centre de Perharidy, Roscoff; Cinthia RAMES, Hôpital Nord, Amiens; Gilles RAULT, Centre de Perharidy, Roscoff; Philippe REIX, Hospices Civils de Lyon—Hôpital femme mère enfant, Lyon; Natacha REMUS, CH Intercommunal de Créteil, Créteil; Martine REYNAUD-GAUBERT,
APHM—Hôpital Nord, Marseille; Bénédicte RICHAUD-THIRIEZ, CHU de Besançon, Besançon; Manuëla SCALBERT-DUJARDIN, Centre Hospitalier de Dunkerque, Dunkerque; Isabelle SERMET-GAUDELUS, APHP—Necker, Paris; Nathalie STREMLER, APHM—Hôpital d'Enfants de la Timone, Marseille; Véronique STORNI, Centre Hospitalier Bretagne Atlantique, Vannes; Aurélie TATOPOU-LOS, Hôpital d'Enfants, Vandoeuvre les Nancy; Guillaume THOUVENIN, APHP-Trousseau, Paris; Françoise TROUSSIER, CHU d'Angers, Angers; Philippe VIGNERON, CH Bretagne Sud, Lorient; Laurence WEISS, CHU de Strasbourg, Strasbourg; Nathalie WIZLA, CHU de Lille, Lille.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. Bell, S.C.; Mall, M.A.; Gutierrez, H.; Macek, M.; Madge, S.; Davies, J.C.; Burgel, P.-R.; Tullis, E.; Castaños, C.; Castellani, C.; et al. The future of cystic fibrosis care: A global perspective. *Lancet Respir. Med.* **2020**, *8*, 65–124; Corrigendum in **2019**, *7*, e40. [CrossRef]

2. Wainwright, C.E.; Elborn, J.S.; Ramsey, B.W.; Marigowda, G.; Huang, X.; Cipolli, M.; Colombo, C.; Davies, J.C.; De Boeck, K.; Flume, P.A.; et al. Lumacaftor–Ivacaftor in Patients with Cystic Fibrosis Homozygous for Phe508del CFTR. *N. Engl. J. Med.* **2015**, *373*, 220–231. [CrossRef]

3. Jones, A.M.; Barry, P.J. Lumacaftor/ivacaftor for patients homozygous for Phe508del-CFTR: Should we curb our enthusiasm? *Thorax* **2015**, *70*, 615–616. [CrossRef] [PubMed]

4. Elborn, J.S.; Ramsey, B.; Wainwright, C.; Boyle, M. Response to: 'Lumacaftor/ivacaftor for patients homozygous for Phe508del-CFTR: Should we curb our enthusiasm? *Thorax* **2015**, *71*, 185–186. [CrossRef] [PubMed]

5. Jennings, M.T.; Dezube, R.; Paranjape, S.; West, N.E.; Hong, G.; Braun, A.; Grant, J.; Merlo, C.A.; Lechtzin, N. An Observational Study of Outcomes and Tolerances in Patients with Cystic Fibrosis Initiated on Lumacaftor/Ivacaftor. *Ann. Am. Thorac. Soc.* **2017**, *14*, 1662–1666. [CrossRef] [PubMed]

6. Konstan, M.W.; McKone, E.F.; Moss, R.B.; Marigowda, G.; Tian, S.; Waltz, D.; Huang, X.; Lubarsky, B.; Rubin, J.; Millar, S.J.; et al. Assessment of safety and efficacy of long-term treatment with combination lumacaftor and ivacaftor therapy in patients with cystic fibrosis homozygous for the F508del-CFTR mutation (PROGRESS): A phase 3, extension study. *Lancet Respir. Med.* **2016**, *5*, 107–118. [CrossRef]

7. Burgel, P.R.; Munck, A.; Durieu, I.; Chiron, R.; Mely, L.; Prevotat, A.; Murris-Espin, M.; Porzio, M.; Abely, M.; Reix, P.; et al. Real-Life Safety and Effectiveness of Lumacaftor–Ivacaftor in Patients with Cystic Fibrosis. *Am. J. Respir. Crit. Care Med.* **2020**, *201*, 188–197. [CrossRef] [PubMed]

8. Taylor-Cousar, J.L.; Jain, M.; Barto, T.L.; Haddad, T.; Atkinson, J.; Tian, S.; Tang, R.; Marigowda, G.; Waltz, D.; Pilewski, J.; et al. Lumacaftor/ivacaftor in patients with cystic fibrosis and advanced lung disease homozygous for F508del-CFTR. *J. Cyst. Fibros.* **2018**, *17*, 228–235. [CrossRef]

9. Burgel, P.-R.; Durieu, I.; Chiron, R.; Mely, L.; Prevotat, A.; Murris-Espin, M.; Porzio, M.; Abely, M.; Reix, P.; Marguet, C.; et al. Clinical response to lumacaftor-ivacaftor in patients with cystic fibrosis according to baseline lung function. *J. Cyst. Fibros.* **2020**, *20*, 220–227. [CrossRef]

10. Aalbers, B.; Groot, K.D.W.-D.; Arets, H.; Hofland, R.; de Kiviet, A.; Ven, M.V.O.-V.D.; Kruijswijk, M.; Schotman, S.; Michel, S.; van der Ent, C.; et al. Clinical effect of lumacaftor/ivacaftor in F508del homozygous CF patients with FEV1 ≥ 90% predicted at baseline. *J. Cyst. Fibros.* **2020**, *19*, 654–658. [CrossRef]

11. Elborn, J.S.; Ramsey, B.W.; Boyle, M.P.; Konstan, M.W.; Huang, X.; Marigowda, G.; Waltz, D.; Wainwright, C. Efficacy and safety of lumacaftor/ivacaftor combination therapy in patients with cystic fibrosis homozygous for Phe508del CFTR by pulmonary function subgroup: A pooled analysis. *Lancet Respir. Med.* **2016**, *4*, 617–626. [CrossRef]

12. Corvol, H.; Mésinèle, J.; Douksieh, I.-H.; Strug, L.J.; Boëlle, P.-Y.; Guillot, L. SLC26A9 Gene Is Associated With Lung Function Response to Ivacaftor in Patients With Cystic Fibrosis. *Front. Pharmacol.* **2018**, *9*, 828. [CrossRef] [PubMed]

13. Strug, L.J.; Gonska, T.; He, G.; Keenan, K.; Ip, W.; Boëlle, P.-Y.; Lin, F.; Panjwani, N.; Gong, J.; Li, W.; et al. Cystic fibrosis gene modifier SLC26A9 modulates airway response to CFTR-directed therapeutics. *Hum. Mol. Genet.* **2016**, *25*, 4590–4600. [CrossRef] [PubMed]

14. Eastman, A.C.; Pace, R.G.; Dang, H.; Aksit, M.A.; Vecchio-Pagán, B.; Lam, A.-T.N.; O'Neal, W.K.; Blackman, S.M.; Knowles, M.R.; Cutting, G.R. SLC26A9 SNP rs7512462 is not associated with lung disease severity or lung function response to ivacaftor in cystic fibrosis patients with G551D-CFTR. *J. Cyst. Fibros.* **2021**, *20*, 851–856. [CrossRef] [PubMed]

15. Pereira, S.V.-N.; Ribeiro, J.D.; Bertuzzo, C.S.; Marson, F.A.L. Association of clinical severity of cystic fibrosis with variants in the SLC gene family (SLC6A14, SLC26A9, SLC11A1 and SLC9A3). *Gene* **2017**, *629*, 117–126. [CrossRef] [PubMed]

16. Sun, L.; Rommens, J.M.; Corvol, H.; Li, W.; Li, X.; Chiang, T.A.; Lin, F.; Dorfman, R.; Busson, P.-F.; Parekh, R.V.; et al. Multiple apical plasma membrane constituents are associated with susceptibility to meconium ileus in individuals with cystic fibrosis. *Nat. Genet.* **2012**, *44*, 562–569. [CrossRef]
17. Corvol, H.; Blackman, S.M.; Boëlle, P.-Y.; Gallins, P.J.; Pace, R.G.; Stonebraker, J.R.; Accurso, F.J.; Clement, A.; Collaco, J.M.; Dang, H.; et al. Genome-wide association meta-analysis identifies five modifier loci of lung disease severity in cystic fibrosis. *Nat. Commun.* **2015**, *6*, 8382. [CrossRef]

18. Li, W.; Soave, D.; Miller, M.R.; Keenan, K.; Lin, F.; Gong, J.; Chiang, T.; Stephenson, A.L.; Durie, P.; Rommens, J.; et al. Unraveling the complex genetic model for cystic fibrosis: Pleiotropic effects of modifier genes on early cystic fibrosis-related morbidities. *Qual. Life Res.* **2014**, *133*, 151–161. [CrossRef]

19. Dehillotte, C.; Lemonnier, L. Registre Français de la Mucoviscidose—Bilan des Données 2018, Vaincre la Mucoviscidose. 2020. Available online: https://www.vaincrelamuco.org/registredelamuco (accessed on 1 June 2021).

20. HAS. Conférence de Consensus—Prise en Charge du Patient Atteint de Mucoviscidose. 2002. Available online: https://www.has-sante.fr/upload/docs/application/pdf/Mucovisc_pneumo_infectio_long.pdf (accessed on 1 June 2021).

21. HAS. Mucoviscidose: Protocole National de Diagnostic et de Soins Pour une Maladie Rare. Guide—Affection de Longue Durée. 2017. Available online: https://www.has-sante.fr/jcms/c_2792719/fr/mucoviscidose (accessed on 1 June 2021).

22. Castellani, C.; Duff, A.J.; Bell, S.C.; Heijerman, H.G.; Munck, A.; Ratjen, F.; Sermet-Gaudelus, I.; Southern, K.W.; Barben, J.; Flume, P.A.; et al. ECFS best practice guidelines: The 2018 revision. *J. Cyst. Fibros.* **2018**, *17*, 153–178. [CrossRef]

23. Quanjer, P.H.; Stanojevic, S.; Cole, T.J.; Baur, X.; Hall, G.L.; Culver, B.H.; Enright, P.L.; Hankinson, J.L.; Ip, M.S.M.; Zheng, J.; et al. Multi-ethnic reference values for spirometry for the 3–95-yr age range: The global lung function 2012 equations. *Eur. Respir. J.* **2012**, *40*, 1324–1343. [CrossRef]

24. Kulich, M.; Rosenfeld, M.; Campbell, J.; Kronmal, R.; Gibson, R.L.; Goss, C.H.; Ramsey, B. Disease-specific Reference Equations for Lung Function in Patients with Cystic Fibrosis. *Am. J. Respir. Crit. Care Med.* **2005**, *172*, 885–891. [CrossRef] [PubMed]

25. Taylor, C.; Commander, C.; Collaco, J.M.; Strug, L.J.; Li, W.; Wright, F.A.; Webel, A.D.; Pace, R.G.; Stonebraker, J.R.; Naughton, K.; et al. A novel lung disease phenotype adjusted for mortality attrition for cystic fibrosis Genetic modifier studies. *Pediatr. Pulmonol.* **2011**, *46*, 857–869. [CrossRef] [PubMed]

26. De Onis, M.; Onyango, A.W.; Borghi, E.; Siyam, A.; Nishida, C.; Siekmann, J. Development of a WHO growth reference for school-aged children and adolescents. *Bull. World Health Organ.* **2007**, *85*, 660–667. [CrossRef] [PubMed]

27. Robins, J.M.; Rotnitzky, A.; Zhao, L.P. Analysis of Semiparametric Regression Models for Repeated Outcomes in the Presence of Missing Data. *J. Am. Stat. Assoc.* **1995**, *90*, 106–121. [CrossRef]

28. Preisser, J.S.; Lohman, K.K.; Rathouz, P.J. Performance of weighted estimating equations for longitudinal binary data with drop-outs missing at random. *Stat. Med.* **2002**, *21*, 3035–3054. [CrossRef]

29. Boëlle, P.; Debray, D.; Guillot, L.; Clement, A.; Corvol, H.; on behalf of the French CF Modifier Gene Study Investigators. Cystic Fibrosis Liver Disease: Outcomes and Risk Factors in a Large Cohort of French Patients. *Hepatology* **2018**, *69*, 1648–1656. [CrossRef]

30. Kmit, A.; Marson, F.; Pereira, S.V.-N.; Vinagre, A.M.; Leite, G.S.; Servidoni, M.F.; Ribeiro, J.D.; Ribeiro, A.F.; Bertuzzo, C.S.; Amaral, M.D. Extent of rescue of F508del-CFTR function by VX-809 and VX-770 in human nasal epithelial cells correlates with SNP rs7512462 in SLC26A9 gene in F508del/F508del Cystic Fibrosis patients. *Biochim. Biophys. Acta (BBA) Mol. Basis Dis.* **2019**, *1865*, 1323–1331. [CrossRef]

31. Durand, E.; Boutin, P.; Meyre, D.; Charles, M.A.; Clement, K.; Dina, C.; Froguel, P. Polymorphisms in the Amino Acid Transporter Solute Carrier Family 6 (Neurotransmitter Transporter) Member 14 Gene Contribute to Polygenic Obesity in French Caucasians. *Diabetes* **2004**, *53*, 2483–2486. [CrossRef]

32. Suviolahti, E.; Oksanen, L.J.; Öhman, M.; Cantor, R.M.; Ridderstråle, M.; Tuomi, T.; Kaprio, J.; Rissanen, A.; Mustajoki, P.; Jousilahti, P.; et al. The SLC6A14 gene shows evidence of association with obesity. *J. Clin. Investig.* **2003**, *112*, 1762–1772. [CrossRef]

33. Miranda, R.C.; Vetter, S.B.; Genro, J.P.; Campagnolo, P.D.; Mattevi, V.S.; Vitolo, M.R.; Almeida, S. SLC6A14 and 5-HTR2C polymorphisms are associated with food intake and nutritional status in children. *Clin. Biochem.* **2015**, *48*, 1277–1282. [CrossRef]

34. Sivaprakasam, S.; Sikder, M.O.; Ramalingam, L.; Kaur, G.; Dufour, J.M.; Moustaid-Moussa, N.; Wachtel, M.S.; Ganapathy, V. SLC6A14 deficiency is linked to obesity, fatty liver, and metabolic syndrome but only under conditions of a high-fat diet. *Biochim. Biophys. Acta (BBA) Mol. Basis Dis.* **2021**, *1867*, 166087. [CrossRef] [PubMed]

35. Debray, D.; Kelly, D.; Houwen, R.; Strandvik, B.; Colombo, C. Best practice guidance for the diagnosis and management of cystic fibrosis-associated liver disease. *J. Cyst. Fibros.* **2011**, *10* (Suppl. 2), S29–S36. [CrossRef]