# Oligogenic Inheritance Underlying Incomplete Penetrance of PROKR2 Mutations in Hypogonadotropic Hypogonadism

Rahma Mkaouar, Lamia Cherif Ben Abdallah, Chokri Naouali, Saida Lahbib, Zinet Turki, Sahar Elouej, Yosra Bouyacoub, Maali Somai, Ken Mcelreavey, Anu Bashamboo, et al.

► To cite this version:

Rahma Mkaouar, Lamia Cherif Ben Abdallah, Chokri Naouali, Saida Lahbib, Zinet Turki, et al. Oligogenic Inheritance Underlying Incomplete Penetrance of PROKR2 Mutations in Hypogonadotropic Hypogonadism. Frontiers in Genetics, 2021, 12, pp.665174. 10.3389/fgene.2021.665174. pasteur-03521812

HAL Id: pasteur-03521812

https://pasteur.hal.science/pasteur-03521812

Submitted on 11 Jan 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
# Oligogenic Inheritance Underlying Incomplete Penetrance of PROKR2 Mutations in Hypogonadotropic Hypogonadism

Rahma Mkaouar¹,²†, Lamia Cherif Ben Abdallah¹†, Chokri Naouali¹, Saida Lahbib¹, Zinet Turki³, Sahar Elouej¹, Yosra Bouyacoub¹, Maali Somai³, Kenneth Mcelreavey⁴, Anu Bashamboo⁴, Sonia Abdelhak¹ and Olfa Messaoud¹*

OPEN ACCESS

**Edited by:**
Luigi Maione,
Université Paris-Saclay, France

**Reviewed by:**
Georgios E. Papadakis,
Centre Hospitalier Universitaire
Vaudois (CHUV), Switzerland
Antonio Agostino Sinisi,
University of Campania Luigi Vanvitelli,
Italy

*Correspondence:
Olfa Messaoud
olfa.messaoud@pasteur.tn;
olfaezeddine@yahoo.fr

† These authors have contributed
equally to this work and share first
authorship

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

**Received:** 07 February 2021
**Accepted:** 05 July 2021
**Published:** 03 September 2021

**Citation:**
Mkaouar R, Abdallah LCB, Naouali C, Lahbib S, Turki Z, Elouej S, Bouyacoub Y, Somai M, Mcelreavey K, Bashamboo A, Abdelhak S and Messaoud O (2021) Oligogenic Inheritance Underlying Incomplete Penetrance of PROKR2 Mutations in Hypogonadotropic Hypogonadism. *Front. Genet.* 12:665174.
doi: 10.3389/fgene.2021.665174

¹ Laboratoire de Génomique Biomédicale et Oncogénétique, Institut Pasteur de Tunis, Tunis, Tunisia, ² Faculté des Sciences Mathématiques, Physiques et Naturelles de Tunis, Université de Tunis El Manar, Tunis, Tunisia, ³ Département d'Endocrinologie et de Technologie Alimentaire, Institut de Nutrition, Tunis, Tunisia, ⁴ Génétique du Développement Humain, Institut Pasteur, Paris, France

The role of the prokineticin 2 pathway in human reproduction, olfactory bulb morphogenesis, and gonadotropin-releasing hormone secretion is well established. Recent studies have highlighted the implication of di/oligogenic inheritance in this disorder. In the present study, we aimed to identify the genetic mechanisms that could explain incomplete penetrance in hypogonadotropic hypogonadism (HH). This study involved two unrelated Tunisian patients with HH, which was triggered by identifying a homozygous p.(Pro290Ser) mutation in the *PROKR2* gene in a girl (HH1) with Kallmann syndrome (KS). The functional effect of this variant has previously been well demonstrated. Unexpectedly, her unaffected father (HH1P) and brother (HH1F) also carried this genetic variation at a homozygous state. In the second family, we identified a heterozygous p.(Lys205del) mutation in *PROKR2*, both in a male patient with normosmic idiopathic IHH (HH12) and his asymptomatic mother. Whole-exome sequencing in the three HH1 family members allowed the identification of additional variants in the prioritized genes. We then carried out digenic combination predictions using the oligogenic resource for variant analysis (ORVAL) software. For HH1, we found the highest number of disease-causing variant pairs. Notably, a *CCDC141* variant (c.2803C > T) was involved in 18 pathogenic digenic combinations. The *CCDC141* variant acts in an autosomal recessive inheritance mode, based on the digenic effect prediction data. For the second patient (HH12), prediction by ORVAL allowed the identification of an interesting pathogenic digenic combination between *DUSP6* and *SEMA7A* genes, predicted as “dual molecular diagnosis.” The *SEMA7A* variant p.(Glu436Lys) is novel and predicted as a VUS by Varsome. Sanger validation revealed the absence of this variant in the healthy mother. We hypothesize that disease expression in HH12 could be induced by the digenic transmission of the *SEMA7A* and *DUSP6* variants or a monogenic inheritance involving only the *SEMA7A*
VUS if further functional assays allow its reclassification into pathogenic. Our findings confirm that homozygous loss-of-function genetic variations are insufficient to cause KS, and that oligogenism is most likely the main transmission mode involved in Congenital Hypogonadotropic Hypogonadism.

**Keywords:** hypogonadotropic hypogonadism, Kallmann syndrome, oligogenism, pathogenic combination, digenic score

## INTRODUCTION

Idiopathic hypogonadotropic hypogonadism (IHH) (MIM ID #146110) is defined as the absence of sex steroid synthesis associated with the lack of appropriate gonadotropin secretion. This leads to a variable degree of impuberism, often diagnosed during childhood or adolescence. When the embryonic migration of gonadotropin-releasing hormone (GnRH) neurons from the nasal placode to their final destination in the hypothalamus is disrupted, the resulting phenotype is Kallmann syndrome (KS), which is defined as the association of HH with hyposmia or anosmia (Sonne and Lopez-Ojeda, 2021). Nineteen genes are known to be involved in KS (ANOS1, FGF8, FGFR1, FGF17, CHD7, IL17RD, DUSP6, SPRY4, FLRT3, KLB, SEMA3A, SEMA3E, NSMF, HS6ST1, WDR11, SOX10, FEZF1, IGSF10, PROK2, and PROKR2) (Boehm et al., 2015; Stamou and Georgopoulos, 2018). Normosmic idiopathic hypogonadotropic hypogonadism (nIHH), which is not associated with anosmia, and results from the dysfunction of the GnRH neurons that successfully completed their embryonic migration to the hypothalamus. As of today, 46 genes have been associated with IHH (Topaloğlu, 2017; Cangiano et al., 2020). The molecular alterations identified in these genes account for 40% of all IHH cases (Gach et al., 2020), thus suggesting that half of the IHH causal genes remain unknown. Molecular alterations have been identified for several neuropeptides or their corresponding receptors, which are involved in the physiology of the gonadotropic axis: GNRHR/GNRH1, KISS1R/GPR54, TAC3/TACR3, and PROK2/PROKR2 (Brioude et al., 2010; Topaloglu and Kotan, 2010). In addition to reproductive dysfunction, nIHH/KS patients may also manifest a variety of other non-reproductive disorders, such as midline facial defects, dental agenesis, renal agenesis, hearing loss, or bimanual synkinesis (Young et al., 2019). IHH may be inherited in an X-linked recessive, autosomal dominant, or autosomal recessive modes of inheritance. IHH has been predominantly detected in sporadic cases (Quaynor et al., 2011; Gach et al., 2020).

The PROK2 gene (NG_008275.1) is located on chromosome 3p21.1 and comprises four exons. The PROKR2 gene (NG_008132.1) maps to chromosome 20p13 and contains seven transmembrane-domain receptors. Both genes, PROK2 and PROKR2, belong to the family of prokineticins, and a group of multifunctional secreted proteins first identified in 2000 by Li et al. (2001). Their involvement in KS was strongly suggested in 2006 when homozygous Prokr2 knockout mice were shown to have hypogonadotropic hypogonadism due to hypothalamic GnRH deficiency and agenesis or hypoplasia of the olfactory bulbs (Matsumoto et al., 2006). This was confirmed in the same

year by the discovery of genetic variations in patients with KS (Dodé et al., 2006). More than 27 genetic variations in PROKR2 (Dode and Rondard, 2013) and over 11 genetic variations in PROK2 have been reported in patients with nIHH or KS (Martin et al., 2011). However, the KS and nIHH genetics are complex and still not well understood (Dode and Hardelin, 2010). In addition, genetic variations in PROK2 and PROKR2 genes reported in KS and nIHH patients were found at heterozygous, homozygous, and compound heterozygous states (Dodé et al., 2006; Hardelin and Dode, 2008). Homozygous loss-of-function genetic variations in the PROK2 gene cause nIHH in mice and humans (Abreu et al., 2008; Schöneberg and Liebscher, 2021). Thus, an autosomal recessive mode of transmission was established (Pitteloud et al., 2007). Later, other studies on different ethnic populations reported a large number of heterozygous genetic variations in PROK2 and PROKR2 genes with considerable clinical and molecular heterogeneity among several patients having both KS and nIHH (Cole et al., 2008; Sarfati et al., 2010a). Heterozygous genetic variations inherited from clinically unaffected first-degree relatives in the PROKR2 gene have been reported in some KS/nIHH patients, which was attributed to digenic/oligogenic transmission rather than a dominant negative effect of monoallelic PROKR2 genetic variations (Dodé et al., 2006; Abreu et al., 2008; Cole et al., 2008; Monnier et al., 2009; Quaynor et al., 2011; Lewkowitz-Shpuntoff et al., 2012). The hypothesis of oligogenic inheritance postulates that disease expression is induced by more than one mutated IHH gene (Pitteloud et al., 2005; Maione et al., 2018). Indeed the PROKR2 gene has been involved in several digenic and trigenic associations such as PROK2/PROKR2, FGFR1/PROKR2, PROK/GNRHR, and PROKR2/CHD7/FEZF1 (Cole et al., 2008; Canto et al., 2009; Sarfati et al., 2010b; Pablo Méndez et al., 2015; Zhang et al., 2020). However, the expression of deleterious alleles is considerably variable if we compare the phenotypes of patients carrying identical variations (Pitteloud et al., 2007; Boehm et al., 2015). Several studies conducted on large cohorts have shown that oligogenic heredity accounts for 2.5–15% of Congenital Hypogonadotropic Hypogonadism (CHH) patients (Cassatella et al., 2018).

Here, we report a clinical and genetic investigation of KS and nIHH in two Tunisian families after excluding the involvement of monogenic inheritance of PROKR2 gene variants.

## PATIENTS AND METHODS

### Patients

The current molecular analysis was conducted for two unrelated HH families (HH1—a woman with KS and HH12—a man with
nIHH). The index cases were recruited at the Endocrinology Department of the National Institute of Nutrition in Tunis. Both patients belong to simplex families (only one family member is affected, with healthy siblings). For the two patients, GnRH deficiency diagnosis was established based on puberty state (absent or incomplete), hormonal tests (testosterone in HH12 and serum gonadotropins levels), response to GnRH, and anterior pituitary function, which was evaluated by measuring the basal levels of free T4, TSH, and prolactin as well as the cortisol peak levels after ACTH stimulation. Pituitary imaging was performed by magnetic resonance imaging (MRI) to exclude acquired causes of nIHH. Olfactory testing for anosmia or hyposmia was also assessed.

## Molecular Investigation

Genomic DNA was extracted from peripheral blood leucocytes by FlexiGene DNA extraction kit (Qiagen) according to the manufacturer's instructions. The coding region and intron-exon boundaries of *PROK2* and *PROKR2* genes as well as exon 13 of the *SEMA7A* gene were amplified as previously reported (Dodé et al., 2006). The PCR products were sequenced with the Big Dye terminator kit (Applied Biosystems, Foster City, CA, United States) using one of the PCR primers on an ABI prism 3,130 DNA Genetic Analyzer (Applied Biosystems) in accordance with the recommendations of the manufacturer.

## In-silico Analysis

To gain insight into the effect of the newly identified genetic variation p.(Lys205del), we carried out an *in silico* analysis. We first tested the impact of genetic variation on splicing by creating or abolishing a splice site. This analysis was performed using the Human Splicing Finder program (Desmet et al., 2009)<fnref n="1" />. We then carried out an *in silico* 3D structure prediction using Phre2 tool (Kelley et al., 2009)<fnref n="2" /> in order to predict the effect of genetic variation on protein folding.

## Exome Analysis

An in-house pipeline analysis was used to generate VCF files. Then, the annotation and the prioritization of variants were carried out using VarAFT tool, version 2.04³. A disease-causing gene approach was used to extract variants located in 46 CHH-related genes extracted from the Online Mendelian Inheritance in Men database and recent literature (Supplementary Table 1). Functionally relevant variants located in exonic genomic regions and splice sites were then selected from the list of variants contained in the VCF file. Variants located in the regulatory regions (UTR, promoters, enhancers, and miRNA binding sites, etc.) flanking the *PROKR2* gene were also screened. Regulatory regions were retrieved using the Genome Browser tool provided by the University of California Santa Cruz database<fnref n="4" /> as well

as the GeneCards Human Gene database<fnref n="5" />. The functional effects of the genetic sequence variants were evaluated by *in silico* prediction tools including SIFT<fnref n="6" />, PolyPhen<fnref n="7" />, FATHMM<fnref n="8" />, MutationTaster<fnref n="9" />, MutationAssessor 1.0¹⁰, PROVEAN v1.1¹<fnref n="1" />, and Varsome<fnref n="12" />. Furthermore, in order to pinpoint the variants that could potentially be associated with disease expression in the index cases (HH1 and HH12), each selected variant of interest was compared to the exome data of 92 unrelated Tunisian individuals stored in the local database. This approach allowed us to avoid bias especially that public databases, such as ESP6500, GnomAD, or 1000Genomes, do not contain relevant information regarding the genetic background of individuals from the Middle East, North Africa populations.

## Digenic Effect Prediction

The Oligogenic Resource for Variant Analysis (ORVAL) online software<fnref n="13" /> was used to predict the effect of variant combinations in disease-causing genes. This tool is based on variant annotation and effect prediction of two predictive methods named VarCOPP and Digenic Effect predictor. The evaluation scores of digenic combinations include the support score (SS) and the classification score (CS). SS informs about the percentage of algorithms integrated in VarCOPP and Digenic Effect predictor that support the pathogenicity of a given bi-locus combination. Its value ranges from 0 to 100. For candidate disease-causing combinations, SS should be equal or greater than 50%. The CS should correspond to the median probability of a digenic combination to be disease-causing. Its value ranges from 0 to 1. Variant pairs are considered pathogenic when CS is greater than 0.489. To estimate the probability that a predicted disease-causing combination is a true positive, VarCOPP encompasses 95 and 99% confidence intervals which are delimited by minimal CS and SS values. The 95% confidence interval requires CS and SS values greater than or equal to 0.55 and 75, respectively. A digenic combination falls inside the 99% confidence interval when CS is greater than 0.74 and SS is equal to 100.

The variant combinations were categorized into three classes. The first one, termed “true digenic” involves two variants in two different genes to induce disease expression. Otherwise, if an individual carries only one variant, he will be considered unaffected. The second class is referred as “composite class.” In this case, an individual carrying only one variant expresses the disease, whereas the second genetic locus is a genetic modifier that modifies the severity of the clinical presentation or the age at onset. The third class, called “dual molecular diagnosis,” requires the presence of variants in two different genes inducing two independent clinical entities (Papadimitriou et al., 2019).

<footnote n="1">¹http://www.umd.be/HSF/</footnote>

<footnote n="2">²http://www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi?id=index</footnote>

<footnote n="3">³http://varaft.eu/</footnote>

<footnote n="4">⁴https://genome.ucsc.edu/</footnote>

<footnote n="5">⁵https://www.genecards.org/</footnote>

<footnote n="6">⁶https://sift.bii.a-star.edu.sg/</footnote>

<footnote n="7">⁷genetics.bwh.harvard</footnote>

<footnote n="8">⁸http://fathmm.biocompute.org.uk</footnote>

<footnote n="9">⁹http://mutationtaster.org/</footnote>

<footnote n="10">¹⁰http://mutationassessor.org/r3/</footnote>

<footnote n="11">¹¹provean.jcvi.org/</footnote>

<footnote n="12">¹²https://varsome.com</footnote>

<footnote n="13">¹³https://orval.ibsquare.be/</footnote>
Only the pathogenic bi-locus combinations for each individual were extracted from the list of variant pairs. Considering that variants found in common between the three HH1 family members would not be at the origin of disease expression, the digenic combinations that were specific to the index case were selected. The specificity criteria were related to the presence or absence of a variant or to its zygosity state, i.e., variants that are homozygous in the index case and heterozygous in the two asymptomatic family members.

# RESULTS

## Clinical Investigation

HH1 is a 21-year-old woman born to first-cousin phenotypically normal parents (Figure 1) and is originating from Central Tunisia. The patient presented congenital anosmia associated with obesity and absence of spontaneous puberty. Her basal

and stimulated GnRH-gonadotropin levels were flat, and her MRI of the hypothalamic-pituitary region showed an aplasia of the olfactory bulbs (Table 1). The 17-beta estradiol levels were not measured for this patient. Her family presents a history of lung cancer, congenital deafness, and colon polyps. Except for the cousin of the index case who had delayed puberty, the other family members had normal pubertal development and olfactory tests. HH12 is a 28-year-old man born to unrelated phenotypically normal parents (Figure 1). The patient suffers from nHHH associated with hypertension and obesity. He first consulted for impuberism. A clinical examination showed that his height is 1.78 m and his weight is 113 kg. The genitalia examination showed a micropenis with a gonadal volume of 3 ml. His GnRH response showed a normal response of follicle-stimulating hormone and luteinizing hormone. The testosterone serum level was 1.9 nmol/l (10–41.5 nmol/l). He presents a normal sense of smell and a normal MRI.

<figure><img src="image_2.png" /><figcaption>TABLE 1 | Clinical and hormonal features of HH1 and HH12 patients.</figcaption></figure>

<table><thead><tr><th rowspan="2">Case</th><th rowspan="2">Age (year)</th><th rowspan="2">Sex</th><th rowspan="2">Diagnosis</th><th rowspan="2">Associated clinical features</th><th rowspan="2">BMI (kg/m²)</th><th colspan="2">LH (IU/L)</th><th colspan="2">FSH (IU/L)</th><th rowspan="2">T (nmol/L)</th></tr><tr><th>Basal</th><th>Peak</th><th>Basal</th><th>Peak</th></tr></thead><tbody><tr><td>HH1</td><td>21</td><td>F</td><td>KS</td><td>Overweight</td><td>25.3</td><td>&lt;0.11</td><td>0.5</td><td>0.11</td><td>1.5</td><td></td></tr><tr><td>HH12</td><td>28</td><td>M</td><td>nHH</td><td>Hypertension and obesity</td><td>35.7</td><td>0.7</td><td>14.3</td><td>0.69</td><td>3.4</td><td>1.9</td></tr></tbody></table>

LH, luteinizing hormone; FSH, follicle-stimulating hormone; BMI, body mass index; T, testosterone; IU/L, international unit per liter.
## Molecular Analysis of the PROKR2 and PROKR2 Genes

For the HH1 family, after having identified a missense genetic variation, c.868C > T; p.(Pro290Ser), in the *PROKR2* gene at a homozygous state in the proband, we sequenced this gene in the relatives: the mother and the sister carried the p.(Pro290Ser) genetic variation at a heterozygous state, but, intriguingly, the father and the brother were homozygous. To confirm these molecular findings, blood sampling, PCR, and sequencing reactions were performed twice. To exclude contamination and paternity problems, all members of this family were genotyped using the identifier kit (Applied Biosystems). The analysis of allele segregation was in favor of paternity inclusion.

For HH12, the molecular analysis allowed the identification of a novel variation, c.614_616del; p.(Lys205del), in the *PROKR2* gene at a heterozygous state. Bioinformatics analysis using the Uniprot database showed that the lysine amino acid at position 205 is buried amid the second extracellular loop and that this residue is highly conserved in Prokr2 from mouse, rat, chimp, dog, cow, *Xenopus tropicalis*, and zebrafish and in the human PROKR1 protein. The deleterious effect of this genetic variation could be explained by the creation of a new acceptor splice site as illustrated by the results of the Human Splicing Finder and Genetic Variation Taster tools (donor gained, score = 0.83mut/=0.55 C). The prediction of the p.(Lys205del) variation effect performed with Phyre2 showed a deleterious effect on the 3D protein structure. The comparison between normal and mutant 3D structures was consistent with the deleterious effect. Indeed the lysine amino acid at position 205 interacts with other residues in the EC2, EC3, the third transmembrane domain (TM3), and the TM; thus, its deletion could reduce the PROKR2 protein stability (Figure 2). The genetic investigation

of the mother of HH12 showed that she was a carrier of the same likely pathogenic mutation [p.(Lys205del)] at a heterozygous state.

## Exome Analysis and Digenic Effect Prediction for the HH1 Family

In total, 43, 39, and 33 genetic variants were identified in HH1, her brother, and her father, respectively. For HH1, seven rare variants were identified with a minor allele frequency (MAF) < 0.1 according to the 1000Genomes and GnomAD databases. Four rare variants were harbored by the brother and the father. Besides the p.(Pro290Ser) in the *PROKR2* gene, the remaining rare variants in each case were predicted as benign or as VUS by Varsome search engine as well as 10 *in silico* prediction programs (Supplementary Tables 2–4). For each set of variants identified in the three family members, we proceeded to the digenic effect prediction of variant pairs using the ORVAL software. A total of 940, 828, and 668 digenic combinations were, respectively, obtained in the index case HH1, the asymptomatic brother HH1F, and the father HH1P. In the proband, 62 variant pairs were predicted as pathogenic (Supplementary Table 5), whereas 44 (Supplementary Table 6) and 41 (Supplementary Table 7) disease-causing combinations were found in the brother and father of HH1, respectively (Figure 3). The confidence intervals supporting the pathogenicity of the digenic combinations ranged from 90 to 99%.

The three family members carried 19 variants in 10 prioritized genes: PROKR2, DUSP6, NTN1, ANOS1, KISS1R, PROP1, DCC, IGSF10, SRA1, and PLXNA1. They shared almost 30 pathogenic digenic combinations involving these variants. In the case of HH1, the PROKR2 variant [p.(Pro290Ser)] yielded the highest number of bi-locus combinations (31%) predicted as disease-causing, with 25 genes and a median pathogenicity score

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

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

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

**FIGURE 2** | Phyre2 prediction of the deleterious effect of p.(Lys205del) on the 3D protein structure. The protein backbone is shown in cartoon representation, with the extracellular domain (EC) in blue, the transmembrane domain (TM) in hot pink, and the intracellular domain in yellow. (A) Mutated protein prediction. (B) WT protein prediction. (C) Zoomed-in image of residue 205 (localized in EC3) shown as black sticks, and the residue interactions around 4 Å are shown as green spheres: 131N and 136T (EC 2); 193T, 194E, 195T, 204E, 206E, 207F, and 223Y (EC3); and 140Y (TM3) and 227I (TM5). The figure was produced with Pymol.
Number of disease-causing digenic combinations in the three HH1 family members

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

**FIGURE 3** | Schematic representation of disease-causing digenic combinations in the three HH1 family members. The histograms correspond to the number of variants in prioritized genes, the total number of pathogenic digenic combinations, and the number of specific pathogenic variant pairs identified in the patient HH1, her asymptomatic brother HH1F, and her father HH1P.

of 0.62. We also found 13 pathogenic digenic combinations (Figure 3) involving variants that were either only carried by HH1 or those that were homozygous in HH1 and heterozygous in the asymptomatic cases (Table 2). These specific variants included four exonic missense variants in *PCSK1* (NM_000439), *IL17RD* (NM_017563), and *CCDC141* (NM_173648) genes, which were homozygous in the proband and heterozygous in the unaffected family members. Another exonic variant in the *FLRT3* gene (NM_013281), as well as two splice site variants in *WDR11* (NM_018117) and *SMCHD1* (NM_015295) genes, was only present in the index case HH1 (Table 2). All these variants were frequent, with a MAF ranging from 0.2 to 0.6 according to the GnomAD database, except for the missense variant [c.2803C > T; p.(Arg935Trp)] in the *CCDC141* gene (MAF = 0.06) (Supplementary Table 2). This variant was absent in the in-house control database. Despite being predicted as benign by Varsome, this variant [p.(Arg935Trp)] was involved in the second most important pathogenic combinations (15%; 13 genes) in HH1, after the *PROKR2* gene variant [p.(Pro290Ser)]. On the other hand, when the *CCDC141* variant was at a heterozygous state, prediction by ORVAL yielded only two pathogenic digenic combinations with *PROKR2* and *DUSP6* variants in HH1F and HH1P cases (Figure 4). The variant pair *CCDC141* (c.2803C > T)-*PROKR2* (c.868C > T) was classified by ORVAL as true digenic. The contribution of the three missense variants in *IL17RD* and *PCSK1* genes, which were homozygous in the index case and heterozygous in the asymptomatic cases (Table 2), to the total number of pathogenic digenic combinations did not differ among the three family members

(Figure 4). This indicates their minor contribution to disease expression in HH1. The splice site variant (c.5476 + 10A > G) in the *SMCDH1* gene (NM_015295), only present in the index case (Table 2), was not involved in any disease-causing digenic combination. The two other genetic variations identified in *FLRT3* (NM_013281) and *WDR11* (NM_018117) genes, present only in HH1 case, represented only 2% of the overall disease-causing variant pairs.

We also evaluated the involvement of other variants located in the flanking regulatory regions of the *PROKR2* gene, which could be responsible for the incomplete penetrance in the family members of HH1. On the basis of predictions provided by the Encyclopedia of DNA Elements and the GeneHancer database, we focused on regulatory regions adjacent to the *PROKR2* gene, likely covered by whole-exome sequencing (WES). The closest promoter and enhancer regions are located in the exon 1-intron 1 junction of *PROKR2*. This region is rich with CpG islands and overlaps with the *AX746654* gene. Two transcription factors, namely, HDAC2 and CTCF, are known to bind to this region. The WES analysis did not reveal any variants in the regulatory region of *PROKR2* nor in the genes encoding its transcription factors.

## Exome Analysis and Digenic Effect Prediction for the HH12 Family

For HH12, 35 variants were found in 25 prioritized genes. Three variants in three genes were rare, including the *PROKR2* gene mutation [p.(Lys205del)], a novel heterozygous missense variant [c.1801G > A; p.(Glu436Lys)] in the *SEMA7A* gene
<table><caption>**TABLE 2** | Disease-causing digenic combinations involving specific variants found in the three HH1 family members.</caption><thead><tr><th>Gene_A</th><th>Gene_B</th><th>GeneA_Aleles<br>(genomic location)</th><th>GeneB_Aleles<br>(genomic location)</th><th>Classification<br>_score</th><th>Support_score</th><th>Predicted_class</th><th>Confiden<br>_zone</th><th>Digenic effect<br>class</th></tr></thead><tbody><tr><td>PROKR2</td><td>PCSK1</td><td>20:5282973:G:A</td><td>5:95728974:G:C</td><td>0.79</td><td>100</td><td>Disease-causing</td><td>95%-zone</td><td>Dual molecular<br>diagnosis</td></tr><tr><td>PROKR2</td><td>CCDC141</td><td>20:5282973:G:A</td><td>2:179721046:G:A</td><td>0.77</td><td>100</td><td>Disease-causing</td><td>95%-zone</td><td>True digenic</td></tr><tr><td>TUBB3</td><td>PROKR2</td><td>16:90002290:T:C</td><td>20:5282973:G:A</td><td>0.73</td><td>99.8</td><td>Disease-causing</td><td>95%-zone</td><td>True digenic</td></tr><tr><td>PROKR2</td><td>KISS1</td><td>20:5282973:G:A</td><td>1:204159787:G:C</td><td>0.71</td><td>98.8</td><td>Disease-causing</td><td>95%-zone</td><td>True digenic</td></tr><tr><td>PROKR2</td><td>KISS1</td><td>20:5282973:G:A</td><td>1:204159922:T:C</td><td>0.7</td><td>98.8</td><td>Disease-causing</td><td>95%-zone</td><td>True digenic</td></tr><tr><td>TACR3</td><td>PROKR2</td><td>4:104510766:G:A</td><td>20:5282973:G:A</td><td>0.62</td><td>87.2</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>CCDC141</td><td>KISS1</td><td>2:179721046:G:A</td><td>1:204159787:G:C</td><td>0.6</td><td>76.4</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>CCDC141</td><td>KISS1</td><td>2:179721046:G:A</td><td>1:204159922:T:C</td><td>0.59</td><td>75</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>TUBB3</td><td>CCDC141</td><td>16:90002290:T:C</td><td>2:179721046:G:A</td><td>0.57</td><td>69.39</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>PROKR2</td><td>IL17RD</td><td>20:5282973:G:A</td><td>3:57138419:G:A</td><td>0.56</td><td>65.6</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>WDR11</td><td>PROKR2</td><td>10:122618148:G:C</td><td>20:5282973:G:A</td><td>0.56</td><td>64.4</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>SPRY4</td><td>PROKR2</td><td>5:141693593:C:T</td><td>20:5282973:G:A</td><td>0.56</td><td>64.4</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>PROKR2</td><td>FLRT3</td><td>20:5282973:G:A</td><td>20:14306953:G:T</td><td>0.54</td><td>51.8</td><td>Disease-causing</td><td>90%-zone</td><td>Dual molecular<br>diagnosis</td></tr><tr><th colspan="9" style="text-align: center;">HH1F</th></tr><tr><td>PROKR2</td><td>PCSK1</td><td>20:5282973:G:A</td><td>5:95728974:G:C</td><td>0.79</td><td>100</td><td>Disease-causing</td><td>95%-zone</td><td>Dual molecular<br>diagnosis</td></tr><tr><td>PROKR2</td><td>CCDC141</td><td>20:5282973:G:A</td><td>2:179721046:G:A</td><td>0.74</td><td>99.6</td><td>Disease-causing</td><td>95%-zone</td><td>True digenic</td></tr><tr><td>SEMA3A</td><td>PROKR2</td><td>7:83606518:G:A</td><td>20:5282973:G:A</td><td>0.7</td><td>97.39</td><td>Disease-causing</td><td>95%-zone</td><td>True digenic</td></tr><tr><td>TACR3</td><td>PROKR2</td><td>4:104510766:G:A</td><td>20:5282973:G:A</td><td>0.62</td><td>87</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>PROKR2</td><td>IGSF10</td><td>20:5282973:G:A</td><td>3:151163138:G:A</td><td>0.57</td><td>68</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>PROKR2</td><td>IL17RD</td><td>20:5282973:G:A</td><td>3:57138419:G:A</td><td>0.56</td><td>64.8</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>PROKR2</td><td>IL17RD</td><td>20:5282973:G:A</td><td>3:57136585:C:T</td><td>0.55</td><td>60.19</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>PROKR2</td><td>SEMA3E</td><td>20:5282973:G:A</td><td>7:82996831:A:T</td><td>0.55</td><td>58.4</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>PROKR2</td><td>CCDC141</td><td>20:5282973:G:A</td><td>2:179753245:C:T</td><td>0.53</td><td>50.4</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><th colspan="9" style="text-align: center;">HH1P</th></tr><tr><td>PROKR2</td><td>PCSK1</td><td>20:5282973:G:A</td><td>5:95728974:G:C</td><td>0.79</td><td>100</td><td>Disease-causing</td><td>95%-zone</td><td>Dual molecular<br>diagnosis</td></tr><tr><td>PROKR2</td><td>HESX1</td><td>20:5282973:G:A</td><td>3:57232504:T:C</td><td>0.76</td><td>100</td><td>Disease-causing</td><td>95%-zone</td><td>True digenic</td></tr><tr><td>PROKR2</td><td>CCDC141</td><td>20:5282973:G:A</td><td>2:179721046:G:A</td><td>0.74</td><td>99.6</td><td>Disease-causing</td><td>95%-zone</td><td>True digenic</td></tr><tr><td>NTN1</td><td>PROKR2</td><td>17:9086293:G:C</td><td>20:5282973:G:A</td><td>0.65</td><td>91</td><td>Disease-causing</td><td>95%-zone</td><td>True digenic</td></tr><tr><td>TACR3</td><td>PROKR2</td><td>4:104510766:G:A</td><td>20:5282973:G:A</td><td>0.62</td><td>87</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>NR5A1</td><td>PROKR2</td><td>9:127262802:C:G</td><td>20:5282973:G:A</td><td>0.59</td><td>74.8</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>PROKR2</td><td>IL17RD</td><td>20:5282973:G:A</td><td>3:57138419:G:A</td><td>0.56</td><td>64.8</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>FEZF1</td><td>PROKR2</td><td>7:121942847:G:C</td><td>20:5282973:G:A</td><td>0.56</td><td>60.8</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr><tr><td>PROKR2</td><td>CCDC141</td><td>20:5282973:G:A</td><td>2:179753245:C:T</td><td>0.53</td><td>50.4</td><td>Disease-causing</td><td>90%-zone</td><td>True digenic</td></tr></tbody></table>

(NM_001146029), as well as a splice site variation in the PLXNA1 gene (NM_032242; MAF = 0.03 in GnomAD). The SEMA7A gene variant was predicted as a VUS according to Varsome, whereas the PLXNA1 gene variation was classified as benign. The SEMA7A variant [p.(Glu436Lys)] was absent in the 92 exomes of our local database (Supplementary Table 8).

The ORVAL prediction revealed five pathogenic variant pairs (confidence interval = 90–95%) involving DUSP6, ANOS1, DCC, PROP1, PLXNA1, and SEMA7A genes (Table 3 and Supplementary Table 9). On the other hand, no disease-causing digenic combinations included the PROKR2 gene variant

p.(Lys205del). The DUSP6 gene [c.340G > T; p.(Val114Leu)] was involved in all five disease-causing digenic combinations. Sanger sequencing showed that the SEMA7A variant [c.1759G > A; p.(Glu587Lys)] was only present in HH12 and absent in his asymptomatic mother (Figure 1). The variants located in the promoter region of PROKR2 were extracted, which revealed one common variant (c.–9 + 342A > G) in intron 1 with a MAF of 0.3 according to GnomAD. Varsome is the only prediction tool that annotated this variant as benign. We have also identified in the CTCF gene, predicted as a transcription factor of the PROKR2 gene, one variant (c.*29T > G) which is also frequent
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Contribution of CHH-related genes in the digenic inheritance in HH1 family members. Number of disease-causing digenic combinations per gene identified in the index case HH1 (A) in the symptomatic brother HH1F (B) and in the asymptomatic father HH1P (C).</figcaption></figure>

(MAF = 0.1) and predicted as benign by Varsome. Furthermore, the ORVAL prediction tool was employed to evaluate the digenic effect of both variants in combination with the p.(Lys205del) identified in *PROKR2* gene. However, no digenic combinations were identified.

DISCUSSION

Since 2006, many studies have reported *PROKR2/PROKR2* genetic variations mainly missense, nonsense, and frameshift genetic alterations (**Figure 5**; Sarfati et al., 2010b; Sbai et al., 2014). Most patients carry monoallelic genetic variations, but some patients may carry bi-allelic variants either in the same or in different genes, indicating a dominant, recessive, or oligogenic inheritance, respectively (Dodé et al., 2006; Pitteloud et al., 2007; Leroy et al., 2008; Sarfati et al., 2010b). The p.(Pro290Ser) variant in *PROKR2* gene has been described in CHH patients at heterozygous (Dodé et al., 2006) and homozygous (Sarfati et al., 2010a) states and is always associated with clinical presentations similar to those seen in HH1 (**Table 4**). The p.(Pro290Ser) mutation is located in the sixth transmembrane (TM) domain. TM domains are important for ligand binding and receptor trafficking (Tan et al., 2004). Functional assays predicted that this genetic variation impairs not only cell surface-targeting of the receptor but also its ability to activate

the Gq Protein (Monnier et al., 2009). In the present study, p.(Pro290Ser) in the *PROKR2* gene was identified both in the proband and her healthy father and brother, all at homozygous states. The identification of homozygous deleterious genetic variations among healthy relatives has been previously described for diverse genetic diseases (Bonyadi et al., 2009; Biegstraaten et al., 2011; Suk et al., 2011). The question that arises is how to explain that such a deleterious genetic variation could be seen in unaffected relatives. The most plausible explanations are (i) possible digenic or oligogenic interactions due to the co-occurrence of additional genetic variants (Sykiotis et al., 2010) and (ii) the involvement of epigenetic factors (Abreu et al., 2010; Dode and Hardelin, 2010; Noel and Kaiser, 2011; McCabe et al., 2013).

Several studies employing next-generation sequencing have revealed that more and more variations known to be causative of genetic diseases can also be found among controls, thus further demonstrating the complexity of heredity (Olson, 2008; Klassen et al., 2011; Nishiguchi and Rivolta, 2012). Hence, the correlation between oligogenicity and the presence of a severe phenotype is often unclear. Sykiotis et al. (2010) reported that the same CHH phenotype seen in a propositus, carrying digenic genetic variations in *FGFR1* and *PROKR2* genes, was also observed in another family member who harbored only the genetic variation in the *FGFR1* gene. This finding proves that carrying more or fewer genetic variations may not always correlate with the
<table><caption>**TABLE 3** | Disease-causing digenic combinations identified in HH12.</caption><thead><tr><th>Gene_A</th><th>Gene_B</th><th>GeneA_Alleles<br>(genomic location)</th><th>GeneB_Alleles<br>(genomic location)</th><th>Classification<br>_score</th><th>Support<br>_score</th><th>Predicted_class</th><th>Confidence_zone</th><th>Digenic_effect_class</th></tr></thead><tbody><tr><td>DUSP6</td><td>ANOS1</td><td>12:89745477:C:A</td><td>X:8504833:C:T</td><td>0.68</td><td>96.8</td><td>Disease-causing</td><td>95% zone</td><td>True digenic</td></tr><tr><td>PLXNA1</td><td>DUSP6</td><td>3:126751689:G:A</td><td>12:89745477:C:A</td><td>0.54</td><td>52.8</td><td>Disease-causing</td><td>90% zone</td><td>Monogenic + modifier</td></tr><tr><td>DCC</td><td>DUSP6</td><td>18:49867224:T:C</td><td>12:89745477:C:A</td><td>0.66</td><td>94.8</td><td>Disease-causing</td><td>95% zone</td><td>Dual molecular diagnosis</td></tr><tr><td>SEMA7A</td><td>DUSP6</td><td>15:74703165:C:T</td><td>12:89745477:C:A</td><td>0.59</td><td>74</td><td>Disease-causing</td><td>90% zone</td><td>Dual molecular diagnosis</td></tr><tr><td>PROP1</td><td>DUSP6</td><td>5:177422823:C:T</td><td>12:89745477:C:A</td><td>0.58</td><td>70.8</td><td>Disease-causing</td><td>90% zone</td><td>?</td></tr></tbody></table>

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

severity of disease expression, which further demonstrates the incomplete penetrance of many suspected variants.

Following exome sequencing in the family of HH1, we identified 42 additional variants in 25 candidate genes in the index case as well as 39 and 33 variants in her asymptomatic brother and father, respectively. The disease-causing digenic combination profiles between the HH1 family members were compared (Figure 4). The results showed that, in addition to the PROKR2 gene variant [p.(Pro290Ser)], a second variant in c.2803C > T in the CCDC141 gene was involved in the second highest number of pathogenic digenic combinations (15%), with 18 other variants in 13 genes. The CCDC141 variant was found at a homozygous state in the patient HH1 and at a heterozygous state in the asymptomatic cases. Our analysis indicated that the zygosity state of the c.2803C > T variant in the CCDC141 gene considerably influenced the rate of

pathogenic combinations. The same variant, when heterozygous in the asymptomatic cases, contributed to only 3% of the total number of pathogenic digenic combinations (Figure 4). In the same context, the implication of the IL17RD and PCSK1 variants in inducing the expression of CHH was ruled out in our propositus as they took part in the same number of digenic combinations with similar classification scores in all three family cases (Figure 4). The CCDC141 gene encodes for a coiled-coil domain-containing protein whose function is not clearly established. In mouse models, it has been shown that ccdc141 is expressed in GnRH neurons and olfactory fibers. The role of CCDC141 in the embryonic migration of GnRH neurons has been highlighted in human patients and mouse models (Hutchins et al., 2016). The prevalence of CCDC141 variants in HH patients was estimated to be 3.3% in the cohort studied by Turan et al. (2017). A nonsense variant p.(Arg724*)
<table><caption>**TABLE 4** | Clinical characteristics of patients having a monoallelic or a biallelic p.(Pro290Ser) mutation.</caption><thead><tr><th>References</th><th>Gender cases</th><th>Genotype</th><th>Familial cases</th><th>BMI</th><th>Sense of smell</th><th>Olfactory bulb MRI</th><th>Cryptorchidism/ microphallus</th><th>Spontaneous puberty</th><th>LH (IU/L) Basal-peak</th><th>FSH (IU/L) Basal-peak</th></tr></thead><tbody><tr><td rowspan="4">Sarfati et al., 2010a,b</td><td>M</td><td>p.(Pro290Ser)/ +</td><td>No</td><td>22.5</td><td>Anosmia</td><td>NA</td><td>No/no</td><td>No</td><td>0.1-0.1</td><td>0.2-1.2</td></tr><tr><td>M</td><td>p.(Pro290Ser)/ +</td><td>No</td><td>28</td><td>Hyposmia</td><td>NA</td><td>No/no</td><td>No</td><td>0.7-10.5</td><td>0.7-2.4</td></tr><tr><td>F</td><td>p.(Pro290Ser)/ +</td><td>Yes</td><td>NA</td><td>Hyposmia</td><td>NA</td><td></td><td>Partial breast development</td><td>4-NA</td><td>4-NA</td></tr><tr><td>F</td><td>p.(Pro290Ser)/p.(Pro290Ser)</td><td>No</td><td>24.3</td><td>Anosmia</td><td>Hypoplasia</td><td></td><td>No</td><td>0.1-0.5</td><td>0.2-1.5</td></tr><tr><td rowspan="2">Dodé et al., 2006</td><td>F</td><td>p.(Pro290Ser)/p.(Pro290Ser)</td><td>No</td><td>21.2</td><td>Anosmia</td><td>Aplasia</td><td></td><td>No</td><td>0.1-2.2</td><td>0.6-4.0</td></tr><tr><td>M</td><td>p.(Pro290Ser)/ +</td><td>No</td><td>28.3</td><td>Anosmia</td><td>Aplasia</td><td>Yes/no</td><td>No</td><td>0.2-0.5</td><td>0.5-2</td></tr><tr><td rowspan="2">Present study</td><td>M</td><td>p.(Pro290Ser)/ +</td><td>No</td><td>20.3</td><td>Hyposmia</td><td>—</td><td>Yes/yes</td><td>No</td><td>0.3-1.3</td><td>0.2-3.4</td></tr><tr><td>F</td><td>p.(Pro290Ser)/p.(Pro290Ser)</td><td>Yes</td><td>25.3</td><td>Anosmia</td><td>Hypoplasia</td><td></td><td>No</td><td>&gt;0.11-0.5</td><td>0.11-1.5</td></tr></tbody></table>

LH, luteinizing hormone; FSH, follicle-stimulating hormone; IU/L, international unit per liter; NA, not applicable.

was identified in 20 probands with KS along with a loss-of-function variant in the *FEZF1* gene (Hutchins et al., 2016). In 75% of the reported families, at least one additional likely pathogenic mutation in another causative gene was identified, hence underscoring the oligogenic involvement of *CCDC141* gene variants in CHH (Turan et al., 2017). The c.2803C > T variant identified in the HH1 family has not been previously associated with CHH. However, it has previously been annotated as damaging in genome-wide association studies (Marsman et al., 2014; van den Berg et al., 2017; Lin et al., 2018). In light of the *in silico* digenic effect prediction, we hypothesize that its implication in 18 disease-causing combinations, compared to only two pathogenic variant pairs in the healthy family members, may be at the origin of disease expression in the index case (HH1), while suggesting that the full penetrance of the *CCDC141* gene variant is associated to an autosomal recessive inheritance mode. It is also worth noting that the brain structure and function present differences between males and females. More than 2,000 genes show expression divergence between the two sexes at all developmental stages, especially during puberty (Shi et al., 2016). This may suggest that the variable expressivity in the HH1 family might be, in part, related to the different sex-specific gene expression profiles between our female patient and the asymptomatic male family members.

For the HH12 family, the healthy mother and her affected son were both heterozygous for the p.(Lys205del) variation in the *PROKR2* gene. The same situation was previously described among several patients (Dodé et al., 2006; Abreu et al., 2008; Cole et al., 2008; Sykiotis et al., 2010). The WES analysis revealed a total of 33 variants [33 variants in coding regions/splice sites (three rare); five in UTR] in 28 genes. The ORVAL prediction revealed five disease-causing digenic combinations involving *DUSP6*, *ANOS1*, *DCC*, *PLXNA1*, *PROP1*, and *SEMA7A* genes (Table 3). The absence of variant combinations involving the *PROKR2* gene variant p.(Lys205del) excludes its implication in digenic inheritance in the index case (HH12). Furthermore, since the variant is novel and has no functional evidence of pathogenicity, it is likely to be benign. Further molecular studies are needed to prove the deleterious character of the *PROKR2* Lys205del variant.

Except for the *SEMA7A* gene variant [p.(Glu436Lys)], mutations identified in *DUSP6*, *ANOS1*, *DCC*, *PLXNA1*, and *PROP1* genes were carried by HH1 family cases (HH1, HH1F, and HH1P) and involved in pathogenic digenic combinations with the *DUSP6* gene variant [p.(Val114Leu)]. Such findings bring into question their involvement in disease expression in HH12. The *SEMA7A* variant [p.(Glu436Lys)] was predicted as VUS by Varsome. Sanger validation revealed the absence of this mutation in the healthy mother. The *SEMA7A* and *DUSP6* genes were implicated in a digenic combination classified as “dual molecular diagnosis” by ORVAL. We hypothesize that the disease expression in HH12 may be explained by the digenic transmission of the *SEMA7A* and *DUSP6* variants or a monogenic inheritance involving only the *SEMA7A* VUS if further functional assays would allow its reclassification into pathogenic.

## CONCLUSION

The presence of homozygous *PROKR2* deleterious genetic variations in asymptomatic first-degree relatives and siblings is apparently insufficient to cause KS/nHH. This strongly suggests the incomplete penetrance trait of the disease. Further investigations are required to clarify the involvement of additional genetic and environmental factors in GnRH deficiency. Exome data analysis allowed us to explore the implication of oligogenism as a possible mechanism involved in the incomplete penetrance witnessed in the investigated families. Nonetheless, the molecular mechanisms that modulate the oligogenic interactions are far from being elucidated since the exact roles of some susceptibility genes in the physiological regulation of the GnRH nervous system are yet to be discovered.

## DATA AVAILABILITY STATEMENT

Processed data related to Sanger sequencing and digenic combinations are available in the article. Raw data related to WES are available from the corresponding author upon reasonable request. Note that in Tunisia, genetic data are considered as
personal private data, therefore the minimal dataset was submitted as **Supplementary Material**. The complete raw data may be made available upon request after obtaining IRB approval.

## ETHICS STATEMENT

Patients were interviewed by both a clinical investigator and a geneticist using a structured hypogonadotropic hypogonadism questionnaire. Written informed consent was obtained from both patients and controls in accordance with the Declaration of Helsinki and after being approved by the Institutional Review Board (Registration numbers IRB00005445 and FWA00010074).

## AUTHOR CONTRIBUTIONS

RM contributed to exome analysis, manuscript drafting, molecular investigation, quality control, and digenic score calculation. LCBA contributed to data analysis, manuscript drafting, recruitment of patients, molecular diagnosis, technical supervision, and quality control. SL contributed to functional *in silico* annotation and exome result interpretation. ZT and MS contributed to clinical diagnosis and characterization, and result evaluation. SE and CN contributed to whole exome data analysis and quality control of bioinformatic pipeline. YB contributed to molecular study, technical assistance, and 3D protein structure determination and interpretation. KM and AB contributed to data evaluation and quality control. SA contributed to study design and supervision, quality control, and manuscript revision. OM contributed to study design and supervision, exome data interpretation, quality control, and manuscript drafting and revision. All authors contributed to the article and approved the submitted version.

## FUNDING

This work has been funded by the Tunisian Ministry of Public Health, the Ministry of Higher Education

## REFERENCES

* Abreu, A. P., Kaiser, U. B., and Latronico, A. C. (2010). The role of prokineticins in the pathogenesis of hypogonadotropic hypogonadism. *Neuroendocrinology* 91, 283-290. doi: 10.1159/000308880

* Abreu, A. P., Trarbach, E. B., de Castro, M., Frade Costa, E. M., Versiani, B., Matias Baptista, M. T., et al. (2008). Loss-of-function mutations in the genes encoding prokineticin-2 or prokineticin receptor-2 cause autosomal recessive Kallmann syndrome. *J. Clin. Endocrinol. Metab.* 93, 4113-4118. doi: 10.1210/jc.2008-0958

* Biegstraaten, M., van Schaik, I. N., Aerts, J. M., Langeveld, M., Mannens, M. M., Bour, L. J., et al. (2011). A monozygotic twin pair with highly discordant Gaucher phenotypes. *Blood Cells Mol. Dis.* 46, 39-41. doi: 10.1016/j.bcmd.2010.10.007

* Boehm, U., Bouloux, P. M., Dattani, M. T., de Roux, N., Dodé, C., Dunkel, L., et al. (2015). European consensus statement on congenital hypogonadotropic hypogonadism—pathogenesis, diagnosis and treatment. *Nat. Rev. Endocrinol.* 11:547. doi: 10.1038/nrendo.2015.112

* Bonyadi, M., Esmaeili, M., Abhari, M., and Lotfi, A. (2009). Mutation analysis of familial GJB2-related deafness in iranian azeri turkish patients. *Genet. Test. Mol. Biomark.* 13, 689-692. doi: 10.1089/gtmb.2009.0026

and Scientific Research (LR20IPT05), and the E.C. grant agreement no. 295097 for FP7 project GM-NCD-Inco.

## ACKNOWLEDGMENTS

The authors would like to thank the patients for their cooperation throughout this study. The authors would also like to thank Catherine Dodé and Suhair Sunoqrot for their fruitful comments.

## SUPPLEMENTARY MATERIAL

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

* **Supplementary Table 1** | List of CHH-related genes prioritized in the exome analysis for HH1 and HH12 families.

* **Supplementary Table 2** | List of variants identified in prioritized genes in the index case HH1.

* **Supplementary Table 3** | List of variants identified in prioritized genes in the asymptomatic brother HH1F.

* **Supplementary Table 4** | List of variants identified in prioritized genes in the asymptomatic father HH1P.

* **Supplementary Table 5** | List of digenic combinations identified in the index case HH1.

* **Supplementary Table 6** | List of digenic combinations identified in the asymptomatic brother HH1F.

* **Supplementary Table 7** | List of digenic combinations identified in the asymptomatic father HH1P.

* **Supplementary Table 8** | List of variants identified in prioritized genes in the index case HH12.

* **Supplementary Table 9** | List of digenic combinations identified in the index case HH12.

* Brioude, F., Bouligand, J., Trabado, S., Francou, B., Salenave, S., Kamenicky, P., et al. (2010). Non-syndromic congenital hypogonadotropic hypogonadism: clinical presentation and genotype-phenotype relationships. *Eur. J. Endocrinol.* 162, 835-851. doi: 10.1530/eje-10-0083

* Cangiano, B., Quinton, R., and Bonomi, M. (2020). Genetics of congenital hypogonadotropic hypogonadism: peculiarities and phenotype of an oligogenic disease. *Hum. Genet.* 2, 1-35.

* Canto, P., Munguia, P., Söderlund, D., Castro, J., and Mendez, J. (2009). Genetic analysis in patients with Kallmann syndrome: coexistence of mutations in prokineticin receptor 2 and KAL1. *J. Androl.* 30, 41-45. doi: 10.2164/jandrol.108.005314

* Cassatella, D., Howard, S. R., Acierno, J. S., Xu, C., Papadakis, G. E., Santoni, F. A., et al. (2018). Congenital hypogonadotropic hypogonadism and constitutional delay of growth and puberty have distinct genetic architectures. *Eur. J. Endocrinol.* 178, 377-388.

* Cole, L. W., Sidis, Y., Zhang, C., Quinton, R., Plummer, L., Pignatelli, D., et al. (2008). Mutations in prokineticin 2 and prokineticin receptor 2 genes in human gonadotrophin-releasing hormone deficiency: molecular genetics and clinical spectrum. *J. Clin. Endocrinol. Metabol.* 93, 3551-3559. doi: 10.1210/jc.2007-2654
Desmet, F. O., Hamroun, D., Lalande, M., Collod-Béroud, G., Claustres, M., Béroud, C., et al. (2009). Human splicing finder: an online bioinformatics tool to predict splicing signals. *Nucleic Acids Res.* 37:e67 doi: 10.1093/nar/gkp215

Dode, C., and Hardelin, J. P. (2010). Clinical genetics of Kallmann syndrome. *Annal. Endocrinol.* 71, 149–157. doi: 10.1016/j.ando.2010.02.005

Dode, C., and Rondard, P. (2013). PROK2/PROKR2 signaling and Kallmann syndrome. *Front. Endocrinol.* 4:19.

Dodé, C., Teixeira, L., Levilliers, J., Fouveaut, C., Bouchard, P., Kottler, M.-L., et al. (2006). Kallmann syndrome: mutations in the genes encoding prokineticin-2 and prokineticin receptor-2. *PLoS Genet.* 2:e175. doi: 10.1371/journal.pgen.0020175

Gach, A., Pinkier, I., Salacińska, K., Szarras-Czapnik, M., Salachna, D., Kucińska, A., et al. (2020). Identification of gene variants in a cohort of hypogonadotropic hypogonadism: diagnostic utility of custom NGS panel and WES in unravelling genetic complexity of the disease. *Mol. Cell. Endocrinol.* 2:110968. doi: 10.1016/j.mce.2020.110968

Hardelin, J. P., and Dode, C. (2008). The complex genetics of Kallmann syndrome: KAL1, FGFR1, FGF8, PROKR2, PROK2, et al. *Sex. Dev.* 2, 181–193. doi: 10.1159/000152034

Hutchins, B. I., Kotan, L. D., Taylor-Burds, C., Ozkan, Y., Cheng, P. J., Gurbuz, F., et al. (2016). CCDC141 mutation identified in anosmic hypogonadotropic hypogonadism (Kallmann syndrome) alters GNRH neuronal migration. *Endocrinology* 157, 1956–1966. doi: 10.1210/en.2015-1846

Kelley, G. G., Chepurny, O. G., Schwede, F., Genieser, H. G., Leech, C. A., Roe, M. W., et al. (2009). Glucose-dependent potentiation of mouse islet insulin secretion by Epac activator 8-pCPT-2'-O-Me-cAMP-AM. *Islets* 1, 260–265. doi: 10.4161/isl.1.3.9645

Klassen, T., Davis, C., Goldman, A., Burgess, D., Chen, T., Wheeler, D., et al. (2011). Exome sequencing of ion channel genes reveals complex profiles confounding personal risk assessment in epilepsy. *Cell* 145, 1036–1048. doi: 10.1016/j.cell.2011.05.025

Leroy, C., Fouveaut, C., Leclercq, S., Jacquemont, S., Boullay, H. D., Lespinasse, J., et al. (2008). Biallelic mutations in the prokineticin-2 gene in two sporadic cases of Kallmann syndrome. *Eur. J. Hum. Genet.* EJHG 16, 865–868. doi: 10.1038/ehjg.2008.15

Lewkowitz-Shpuntoff, H. M., Hughes, V. A., Plummer, L., Au, M. G., Doty, R. L., Seminara, S. B., et al. (2012). Olfactory phenotypic spectrum in idiopathic hypogonadotropic hypogonadism: pathophysiological and genetic implications. *J. Clin. Endocrinol. Metab.* 97, E136–E144.

Li, M., Bullock, C. M., Knauer, D. J., Ehlert, F. J., and Zhou, Q. Y. (2001). Identification of two prokineticin cDNAs: recombinant proteins potently contract gastrointestinal smooth muscle. *Mol. Pharmacol.* 59, 692–698. doi: 10.1124/mol.59.4.692

Lin, H., van Setten, J., Smith, A. V., Bihlmeyer, N. A., Warren, H. R., Brody, J. A., et al. (2018). Common and rare coding genetic variation underlying the electrocardiographic PR interval. *Circulat. Genomic Precis. Med.* 11: e002037.

Maione, L., Dwyer, A. A., Francou, B., Guiochon-Mantel, A., Binart, N., Bouligand, J., et al. (2018). Genetic counseling for congenital hypogonadotropic hypogonadism and Kallmann syndrome: new challenges in the era of oligogenism and next generation sequencing. *Eur. J. Endocrinol.* 178, R55–R80.

Marsman, R. F., Tan, H. L., and Bezzina, C. R. (2014). Genetics of sudden cardiac death caused by ventricular arrhythmias. *Nat. Rev. Cardiol.* 11:96. doi: 10.1038/nrcardio.2013.186

Martin, C., Balasubramanian, R., Dwyer, A. A., Au, M. G., Sidis, Y., Kaiser, U. B., et al. (2011). The role of the prokineticin 2 pathway in human reproduction: evidence from the study of human and murine gene mutations. *Endocr. Rev.* 32, 225–246. doi: 10.1210/er.2010-0007

Matsumoto, S., Yamazaki, C., Masumoto, K.-h., Nagano, M., Naito, M., Soga, T., et al. (2006). Abnormal development of the olfactory bulb and reproductive system in mice lacking prokineticin receptor PKR2. *Proc. Natl. Acad. Sci. U.S.A.* 103, 4140–4145. doi: 10.1073/pnas.0508881103

McCabe, M. J., Gaston-Massuet, C., Gregory, L. C., Alatzoglou, K. S., Tziaferi, V., Sbai, O., et al. (2013). Variations in PROKR2, but not PROK2, are associated with hypopituitarism and septo-optic dysplasia. *J. Clin. Endocrinol. Metabol.* 98, E547–E557.

Monnier, C., Dode, C., Fabre, L., Teixeira, L., Labesse, G., Pin, J. P., et al. (2009). PROKR2 missense mutations associated with Kallmann syndrome impair receptor signalling activity. *Hum. Mol. Genet.* 18, 75–81. doi: 10.1093/hmg/ddn318

Nishiguchi, K. M., and Rivolta, C. (2012). Genes associated with retinitis pigmentosa and allied diseases are frequently mutated in the general population. *PLoS One* 7:e41902. doi: 10.1371/journal.pone.0041902

Noel, S. D., and Kaiser, U. B. (2011). G protein-coupled receptors involved in GNRH regulation: molecular insights from human disease. *Mol. Cell. Endocrinol.* 346, 91–101. doi: 10.1016/j.mce.2011.06.022

Olson, M. V. (2008). Human genetics: Dr Watson's base pairs. *Nature* 452, 819–820. doi: 10.1038/452819a

Pablo Méndez, J., Zenteno, J. C., Coronel, A., Soriano-Ursúa, M. A., Valencia-Villalvazo, E. Y., Soderlund, D., et al. (2015). Triallelic digenic mutation in the prokineticin 2 and GNRH receptor genes in two brothers with normosmic congenital hypogonadotropic hypogonadism. *Endocr. Res.* 40, 166–171. doi: 10.3109/07435800.2014.982327

Papadimitriou, S., Gazzo, A., Versbraegen, N., Nachtegael, C., Aerts, J., Moreau, Y., et al. (2019). Predicting disease-causing variant combinations. *Proc. Natl. Acad. Sci. U.S.A.* 116, 11878–11887.

Pitteloud, N., Acierno, J. S. Jr., Meysing, A. U., Dwyer, A. A., Hayes, F. J., Crowley, W. F. Jr., et al. (2005). Reversible Kallmann syndrome, delayed puberty, and isolated anosmia occurring in a single family with a mutation in the fibroblast growth factor receptor 1 gene. *J. Clin. Endocrinol. Metab.* 90, 1317–1322. doi: 10.1210/jc.2004-1361

Pitteloud, N., Zhang, C., Pignatelli, D., Li, J.-D., Raivio, T., Cole, L. W., et al. (2007). Loss-of-function mutation in the prokineticin 2 gene causes Kallmann syndrome and normosmic idiopathic hypogonadotropic hypogonadism. *Proc. Natl. Acad. Sci. U.S.A.* 104, 17447–17452. doi: 10.1073/pnas.0707173104

Quaynor, S. D., Kim, H. G., Cappello, E. M., Williams, T., Chorich, L. P., Bick, D. P., et al. (2011). The prevalence of digenic mutations in patients with normosmic hypogonadotropic hypogonadism and Kallmann syndrome. *Fertil. Steril.* 96, 1424–1430.e1426.

Sarfati, J., Anne, G.-M., Philippe, R., Isabelle, A., Alfons, G.-P., Slawomir, W., et al. (2010a). A comparative phenotypic study of Kallmann syndrome patients carrying monoallelic and biallelic mutations in the prokineticin 2 or prokineticin receptor 2 genes. *J. Clin. Endocrinol. Metab.* 95, 659–669. doi: 10.1210/jc.2009-0843

Sarfati, J., Dode, C., and Young, J. (2010b). Kallmann syndrome caused by mutations in the PROK2 and PROKR2 genes: pathophysiology and genotype-phenotype correlations. *Front. Horm. Res.* 39:121–132. doi: 10.1159/000312698

Sbai, O., Monnier, C., Dodé, C., Pin, J. P., Hardelin, J. P., and Rondard, P. (2014). Biased signaling through G-protein-coupled PROKR2 receptors harboring missense mutations. *FASEB J.* 28, 3734–3744. doi: 10.1096/fj.13-243402

Schöneberg, T., and Liebscher, I. (2021). Mutations in G protein-coupled receptors: mechanisms, pathophysiology and potential therapeutic approaches. *Pharmacol. Rev.* 73, 89–119. doi: 10.1124/pharmrev.120.000011

Shi, L., Zhang, Z., and Su, B. (2016). Sex biased gene expression profiling of human brains at major developmental stages. *Sci. Rep.* 6, 1–9.

Sonne, J., and Lopez-Ojeda, W. (2021). "Kallmann syndrome", in *StatPearls [Internet]*. Treasure Island, FL: StatPearls Publishing.

Stamou, M. I., and Georgopoulos, N. A. (2018). Kallmann syndrome: phenotype and genotype of hypogonadotropic hypogonadism. *Metabolism* 86, 124–134. doi: 10.1016/j.metabol.2017.10.012

Suk, E.-K., McEwen, G. K., Duitama, J., Nowick, K., Schulz, S., Palczewski, S., et al. (2011). A comprehensively molecular haplotype-resolved genome of a European individual. *Genome Res.* 21, 1672–1685. doi: 10.1101/gr.125047.111

Sykiotis, G. P., Plummer, L., Hughes, V. A., Au, M., Durrani, S., Nayak-Young, S., et al. (2010). Oligogenic basis of isolated gonadotropin-releasing hormone deficiency. *Proc. Natl. Acad. Sci. U.S.A.* 107, 15140–15144.

Tan, C. M., Brady, A. E., Nickols, H. H., Wang, Q., and Limbird, L. E. (2004). Membrane trafficking of G protein-coupled receptors. *Annu. Rev. Pharmacol. Toxicol.* 44, 559–609.

Topaloğlu, A. K. (2017). Update on the genetics of idiopathic hypogonadotropic hypogonadism. *J. Clin. Res. Pediatr. Endocrinol.* 9:113. doi: 10.4274/jcrpe.2017. s010
Topaloglu, A. K., and Kotan, L. D. (2010). Molecular causes of hypogonadotropic hypogonadism. *Curr. Opin. Obstetr. Gynecol.* 22, 264–270. doi: 10.1097/gco. 0b013e32833bb425

Turan, I., Hutchins, B. I., Hacihamdioglu, B., Kotan, L. D., Gurbuz, F., Ulubay, A., et al. (2017). CCDC141 mutations in idiopathic hypogonadotropic hypogonadism. *J. Clin. Endocrinol. Metabol.* 102, 1816–1825. doi: 10.1210/jc. 2016-3391

van den Berg, M. E., Warren, H. R., Cabrera, C. P., Verweij, N., Mifsud, B., Haessler, J., et al. (2017). Discovery of novel heart rate-associated loci using the Exome Chip. *Hum. Mol. Genet.* 26, 2346–2363.

Young, J., Xu, C., Papadakis, G. E., Acierno, J. S., Maione, L., Hietamäki, J., et al. (2019). Clinical management of congenital hypogonadotropic hypogonadism. *Endocr. Rev.* 40, 669–710.

Zhang, Q., He, H. -h., Janjua, M. U., Wang, F., Yang, Y. -b., Mo, Z. -h., et al. (2020). Identification of two novel mutations in three Chinese families with Kallmann syndrome using whole exome sequencing. *Andrologia* 52: e13594.

**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.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2021 Mkaouar, Abdallah, Naouali, Lahbib, Turki, Elouej, Bouyacoub, Somai, Mcelreavey, Bashamboo, Abdelhak and Messaoud. 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.