Germline or somatic GPR101 duplication leads to X-linked acrogigantism: a clinico-pathological and genetic study
- Donato Iacovazzo1,
- Richard Caswell2,
- Benjamin Bunce2,
- Sian Jose3,
- Bo Yuan4,
- Laura C. Hernández-Ramírez1, 5,
- Sonal Kapur1,
- Francisca Caimari1,
- Jane Evanson6,
- Francesco Ferraù1,
- Mary N. Dang1,
- Plamena Gabrovska1,
- Sarah J. Larkin7,
- Olaf Ansorge7,
- Celia Rodd8,
- Mary L. Vance9,
- Claudia Ramírez-Renteria10,
- Moisés Mercado10,
- Anthony P. Goldstone11,
- Michael Buchfelder12,
- Christine P. Burren13,
- Alper Gurlek14,
- Pinaki Dutta15,
- Catherine S. Choong16,
- Timothy Cheetham17,
- Giampaolo Trivellin5,
- Constantine A. Stratakis5,
- Maria-Beatriz Lopes18,
- Ashley B. Grossman19,
- Jacqueline Trouillas20,
- James R. Lupski4, 21, 22, 23,
- Sian Ellard2,
- Julian R. Sampson†3,
- Federico Roncaroli†24 and
- Márta Korbonits†1Email author
© The Author(s). 2016
Received: 17 March 2016
Accepted: 19 May 2016
Published: 1 June 2016
Non-syndromic pituitary gigantism can result from AIP mutations or the recently identified Xq26.3 microduplication causing X-linked acrogigantism (XLAG). Within Xq26.3, GPR101 is believed to be the causative gene, and the c.924G > C (p.E308D) variant in this orphan G protein-coupled receptor has been suggested to play a role in the pathogenesis of acromegaly.
We studied 153 patients (58 females and 95 males) with pituitary gigantism. AIP mutation-negative cases were screened for GPR101 duplication through copy number variation droplet digital PCR and high-density aCGH. The genetic, clinical and histopathological features of XLAG patients were studied in detail. 395 peripheral blood and 193 pituitary tumor DNA samples from acromegaly patients were tested for GPR101 variants.
We identified 12 patients (10 females and 2 males; 7.8 %) with XLAG. In one subject, the duplicated region only contained GPR101, but not the other three genes found to be duplicated in the previously reported patients, defining a new smallest region of overlap of duplications. While females presented with germline mutations, the two male patients harbored the mutation in a mosaic state. Nine patients had pituitary adenomas, while three had hyperplasia. The comparison of the features of XLAG, AIP-positive and GPR101&AIP-negative patients revealed significant differences in sex distribution, age at onset, height, prolactin co-secretion and histological features. The pathological features of XLAG-related adenomas were remarkably similar. These tumors had a sinusoidal and lobular architecture. Sparsely and densely granulated somatotrophs were admixed with lactotrophs; follicle-like structures and calcifications were commonly observed. Patients with sporadic or familial acromegaly did not have an increased prevalence of the c.924G > C (p.E308D) GPR101 variant compared to public databases.
In conclusion, XLAG can result from germline or somatic duplication of GPR101. Duplication of GPR101 alone is sufficient for the development of XLAG, implicating it as the causative gene within the Xq26.3 region. The pathological features of XLAG-associated pituitary adenomas are typical and, together with the clinical phenotype, should prompt genetic testing.
KeywordsXLAG Gigantism GPR101 CNV mutation Pituitary
X-linked acrogigantism (XLAG) is a recently identified cause of early-onset pituitary gigantism [1, 2]. The condition often manifests during the first year of life, occurs more frequently in females usually as sporadic disease with only two families reported so far . XLAG patients develop pituitary hyperplasia or mixed somatotroph/lactotroph adenomas both resulting in significant growth hormone (GH) excess . All previously published patients harbor Xq26.3 microduplications encompassing a region of approximately 500Kb [1, 2]. This region contains the locus of the GPR101 gene, encoding an orphan G protein-coupled receptor (GPCR) that is significantly overexpressed in the pituitary samples of XLAG patients. In addition, the c.924G > C (p.E308D) GPR101 missense variant was identified in 4.4 % of a series of patients with sporadic acromegaly. This variant was suggested to represent a disease-associated mutation, as it increases cell proliferation and GH release in vitro .
The genetic background of non-syndromic gigantism also includes inactivating germline mutations in the aryl hydrocarbon receptor-interacting protein (AIP) gene. These mutations are found in 29–50 % of gigantism cases [3–5], causing a low-penetrance disease typically manifesting clinically in the second decade of life, either sporadically or in the setting of Familial Isolated Pituitary Adenoma (FIPA).
In this study, we aimed to evaluate the prevalence of Xq26.3 microduplication, a copy number variation (CNV) gain including GPR101, in a large cohort of 153 patients with non-syndromic pituitary gigantism, who were all screened for AIP mutations. We provide a detailed clinical and histopathological description of this condition and compare the clinical characteristics of gigantism patients with GPR101, AIP and without GPR101&AIP mutations. We also assessed the prevalence of GPR101 variants in a large series of patients with acromegaly.
Materials and methods
The study cohort consisted of 153 patients (58 females and 95 males) diagnosed with pituitary gigantism. Patients were recruited via the International FIPA consortium network (http://www.fipapatients.org/fipaconsortium). Pituitary gigantism was defined as GH excess associated with accelerated growth velocity (> + 2 SDS, standard deviation score) or abnormally tall stature (> + 3 SDS above the normal mean height, i.e. Z-score > +3, or > +2 SDS above the mid-parental height) [3, 6]. Detailed clinical and biochemical data were obtained from the referring physicians. SDS for height and BMI were calculated based on the Centers for Disease Control (CDC) growth charts and, when available, by country-specific growth charts. Other conditions predisposing to pituitary gigantism, including McCune-Albright syndrome, Multiple Endocrine Neoplasia type 1 (MEN1) and Carney complex, were excluded based on the clinical data and, when appropriate, by genetic testing. Two patients from our cohort  were previously reported to carry an Xq26.3 microduplication [1, 2]. Additionally, the prevalence of GPR101 variants was studied in a cohort of 579 acromegaly patients (sporadic and familial) recruited through the FIPA consortium. The study was approved by the local Ethical Committee and informed consent was obtained from all individual participants included in the study.
The AIP gene was tested in all patients through Sanger sequencing and, in most patients, by MLPA dosage as well, as previously described [7, 8]. Patients who were not found to carry an AIP mutation were screened for GPR101 duplication through CNV droplet digital PCR (CNV ddPCR) for GPR101 in leukocyte- or saliva-derived DNA and, in a subset of patients, in DNA isolated from other sources (pituitary, palatine tonsil, skin and buccal cells). When available, DNA samples from both parents of identified XLAG patients were also tested for GPR101 duplication using the CNV ddPCR. Positive results were confirmed by means of standard and high-density array comparative genomic hybridization (HD-aCGH) and breakpoint junction analysis. Further methodological details are shown in the Additional file 1.
Three hundred and ninety-five leukocyte- and 193 pituitary tumor-derived DNA samples from patients with acromegaly (total number of patients = 579) were tested for the c.924G > C (p.E308D) GPR101 variant through Sanger sequencing. Leukocyte-derived DNA samples were also tested for the previously reported c.1098C > A (p.D366E) variant [9, 10]. Sanger sequencing of the whole coding region of GPR101 was also performed in a subset of DNA samples isolated from 42 randomly selected somatotroph adenomas. Primer sequences are available upon request.
Pituitary adenoma tissue from six XLAG patients was available for review and for further studies. Details of the methods for immunohistochemistry, double immunofluorescence and electron microscopy are reported in the Additional file 1. Mitoses were counted in at least 30 fields at the magnification of x40 using a Nikon Eclipse E600 microscope (Nikon UK, Kingston upon Thames, UK) equipped with a Plan Fluor x40/0.75 objective. Quantification of immunoreactions was performed on images taken at the magnification of x20 with a Leica DM5500 microscope (Leica, Wetzlar, Germany) coupled with a Leica DFC295 camera (Leica). Scoring was performed independently by two co-authors (FR and DI). Ki-67 labelling index was evaluated in at least five representative areas counting at least 1000 cells; the value was given as the percentage of positive nuclei. Somatostatin receptor (SSTR) expression was scored as previously described , taking into account the localization and extent of the staining .
Parametric data are presented as mean ± standard deviation (SD), and nonparametric data as median [interquartile range, IQR]. Normal distribution was assessed using the Shapiro-Wilk test. Data were analyzed through univariate tests (Chi-square, ANOVA and Kruskal-Wallis tests with Bonferroni and Dunn’s post-hoc tests, as appropriate) using the software Prism v5 (GraphPad Software Inc, La Jolla, CA, USA). Significance was set for P values <0.05.
Genetics and family history
Of the 153 patients, 63 had AIP mutations (AIPpos, 41.2 %) (27 were sporadic and 36 familial) and 78 patients (GPR101&AIPneg, 51 %) were negative for both AIP mutations and GPR101 duplication CNV. Twelve patients (10 females, 2 males) with GPR101 duplication were identified, accounting for 7.8 % of the whole cohort and 17.2 % of females. The duplication appeared to be germline in the ten female patients, while the two male patients harbored the mutation in a mosaic state. In one of these patients (case IX), the duplication was identified by means of the CNV ddPCR in pituitary-, skin- and palatine tonsil-derived but not in leukocyte-, saliva- or buccal cell-derived DNA . In the other case (case VIII), the duplication was found in the pituitary tissue, while saliva-derived DNA tested negative. Analysis of leukocyte-derived DNA in this patient using the CNV ddPCR showed intermediate results between normal and duplicated dosage, suggesting the presence of a heterogeneous blood cell population as a result of somatic mosaicism (Additional file 1: Figure S1).
Clinical and biochemical features
XLAG patients were generally born at full term and had normal birth length and weight. One patient (case IX) was born large for gestational age (+2.4 SDS). All patients presented with accelerated growth starting as early as seven months of age with a median age at onset of 1.9 years [1.1–2.4]. The median age at diagnosis of gigantism was 4.4 years [2.7–6.7]. Median delay between onset of the disease and diagnosis was 2.6 years [1–3.7]. The median height SDS was +5.4 [4–6.3]. BMI SDS was increased in four out of seven patients with available data, with a median of +2.2 [1.2–3]. Details of symptoms at presentation are shown in the Additional file 1.
Basal GH levels were increased in all patients. IGF-1 levels were increased at 2.9xULN (upper limit of normal) [2–3.9]. Oral glucose tolerance test showed unsuppressed GH (available for 8/12 patients) with a mean change in GH levels of −14.5 % ± 25.5; two patients presented a paradoxical increase of GH levels in response to the glucose load. Patient IV showed a paradoxical rise of GH after thyrotropin-releasing hormone (TRH), and an increase of both GH and prolactin (PRL) levels after the administration of growth hormone-releasing hormone (GHRH) . PRL was elevated in ten patients with a median of 6.9xULN [3.3–11.3]. Circulating GHRH was measured in the three patients with pituitary hyperplasia and was within the reference range. None of the other pituitary axes was affected at diagnosis in any of the patients. No differences in clinical and biochemical parameters were evident between patients carrying germline or somatic GPR101 duplication.
Tumor size and extension
Nine of the 12 XLAG patients (75 %) had macroadenomas. Median maximum tumor size was 18 mm [14–25.5]. All the adenomas showed suprasellar extension, and three lesions extended into the cavernous sinus. MRI of three patients (patient II, IV and IX) (25 % of our XLAG cases) showed diffuse enlargement of the gland suggestive of pituitary hyperplasia rather than a distinct adenoma. Representative MRI images are shown in the Additional file 1: Figure S2.
Histological and immunohistochemical features
Scattered adenoma cells showed expression of the common α-subunit (Fig. 4c), while the immunoreactions for ACTH and FSH, LH and TSH β-subunits were negative. Double immunofluorescence highlighted two distinct populations of GH- and PRL-positive cells and no co-localization of the two hormones (Fig. 4d). The transcription factor PIT-1 was expressed in over 90 % of tumor cells (Fig. 5c). The Ki-67 labelling index was lower than 3 % in all cases (Fig. 5d). Weak nuclear p53 expression was noted in all XLAG adenomas, and the degree of expression was similar to sporadic somatotroph adenomas tested as controls. SSTR2a and SSTR5 expression was variable and both cytoplasmic and membranous (scores 1 and 2), although none of the adenomas reached score 3 indicating circumferential membranous staining in more than 50 % of tumor cells. AIP expression was moderate to strong while GHRH immunostaining was negative in all cases. Histological examination of the previously published cases of hyperplasia (case IV and IX) revealed similar features, including enlarged acini with preservation of the reticulin network, substantial increase in somatotroph and lactotroph cells, and scattered cells expressing ACTH and FSH, LH and TSH β-subunits [13, 14].
On electron microscopy, the main cell population in XLAG adenomas was represented by DG cells with well-developed Golgi complexes and rough endoplasmic reticulum containing granules ranging between 250-600 nm in size resembling DG somatotroph cells. Most of the SG cells resembled lactotrophs, showing peripheral granules measuring between 150 and 300 nm and abundant endoplasmic reticulum, while scattered SG somatotroph cells containing fibrous bodies and granules between 200 and 450 nm were also identified (Figure S4). In the hyperplasia cases, a non-tumoral adenohypophysis was seen, dominated by cells with ultrastructural features of somatotrophs and lactotrophs. Many cells containing large pleomorphic secretory granules and exocytosis were identified, showing co-localization of GH and PRL within the same cell and the same granules, as assessed by immunoelectron microscopy in case IV , in keeping with mammosomatotroph cells.
Treatment and outcomes
The clinical features, treatment and outcomes of XLAG patients are detailed in the Additional file 1: Table S2. The median number of treatments per patient was 3.5 [2–4.7]. Five patients received medical therapy as first-line treatment, either with somatostatin analogues (SSAs) combined with dopamine agonists (DAs) (four patients) or DAs alone (one patient), with overall poor results. The GH receptor antagonist pegvisomant was effective in normalizing IGF-1 levels and growth velocity in three of the four patients in whom it was used. Overall, eight out of the 11 patients with surgical or radiotherapeutical intervention developed partial or complete hypopituitarism. Median duration of follow-up was 10.6 years [4.4–22.9]. At the last follow up, overall disease control had been achieved in 11/12 patients, including remission (n = 4), disease control (normal age-adjusted IGF-1 in patients on pegvisomant, n = 3) and partial control (normal age-adjusted IGF-1 and random GH >1 ng/ml in patients on SSAs and/or DAs, n = 4). No other tumors were found to date; the oldest patient with GPR101 duplication in our series is now 50 years old and has not developed any other disease manifestation.
Comparison of clinical characteristics of XLAG, AIPpos and GPR101&AIPneg patients
Sex distribution, clinical and biochemical parameters, and MRI findings observed in XLAG patients were compared to those of AIPpos and GPR101&AIPneg patients (Additional file 1: Table S3). Ten of the 12 XLAG patients were females, as opposed to 31.7 % of AIPpos and 35.9 % of GPR101&AIPneg cases (P < 0.001 for both comparisons). The median age at onset and diagnosis in XLAG patients was significantly earlier compared to AIPpos and GPR101&AIPneg patients (P < 0.001 for both comparisons). The height SDS was also significantly higher in XLAG patients compared to GPR101&AIPneg (P < 0.01) and AIPpos cases (P < 0.05).
No significant difference was found in IGF-1 levels at diagnosis. Hyperprolactinemia was significantly more prevalent in XLAG (83.3 %), compared to AIPpos (23.5 %, P < 0.001) and GPR101&AIPneg patients (32.3 %, P < 0.01). Median maximum tumor diameter and proportion of macroadenomas was not significantly different in the three groups, although no giant adenomas were found in XLAG patients, while they represented 25–30 % of the cases in the other two groups (P < 0.05 XLAG vs GPR101&AIPneg). Pituitary hyperplasia was significantly more common in XLAG patients (25 %), compared to AIPpos (2.6 %, P < 0.05) and GPR101&AIPneg cases (none, P < 0.01). No difference was found in the rate of suprasellar extension and cavernous sinus invasion. Pituitary apoplexy was not observed in XLAG patients, while it occurred in 14.5 % of AIPpos and 2.9 % of GPR101&AIPneg cases. There was no significant difference in the rate of hypopituitarism and number of treatments.
GPR101 variants in acromegaly patients
We sequenced DNA samples from 579 acromegaly patients (395 leukocyte- and 193 pituitary tumor-derived) and identified four patients (0.69 %) harboring the germline GPR101 c.924G > C (p.E308D) variant. None of these patients had a family history of pituitary adenoma. The allele frequency of the c.924G > C (p.E308D) variant in our series (0.45 %) was similar to that reported in the Exome Aggregation Consortium Database (ExAC) (0.37 %, P = 0.69). The c.1098C > A (p.D366E) variant was not identified in any of the 395 leukocyte-derived DNA samples. The full coding region of the GPR101 gene was sequenced in 42 unselected sporadic somatotroph adenomas. Two common single nucleotide variants were identified: c.370G > T (p.V124L) (rs1190736, minor allele frequency ExAC database 38 %) and c.1127 T > C (p.L376P) (rs5931046, minor allele frequency 17.2 %). No rare or novel variants were found.
In this large international cohort of patients with non-syndromic pituitary gigantism, we identified germline or somatic GPR101 duplication causing XLAG in 8 % of patients. This frequency, consistent with previously published findings , supports the need for testing for GPR101 duplication in patients with early-onset pituitary gigantism.
Our results provide significant and clinically relevant mechanistic insights on XLAG. Firstly, we provide proof that duplication of GPR101 is sufficient to cause the XLAG phenotype. HD-aCGH in one of the patients from our series has allowed the definition of a new smallest region of overlap in XLAG patients to an area encompassing GPR101 only, but not the other three genes (CD40LG, ARHGEF6 and RBMX) previously identified in the duplicated region [1, 2, 15, 16]. Furthermore, our results indicate that, as such small duplications cannot be found on a standard CGH array, alternative diagnostic methods should be employed, such as CNV ddPCR for GPR101 or HD-aCGH.
Secondly, as recently described [13, 16], XLAG can result from a duplication occurring during an early postzygotic stage, resulting in somatic mosaicism. The only two male subjects in our series were mosaic for the GPR101 duplication. In one of these patients, the CNV ddPCR did not show the presence of the duplication in leukocyte-derived DNA , indicating that, in the presence of a clinical phenotype suggestive of XLAG, analysis of pituitary- or other tissues-derived DNA samples should be considered in order to confirm the diagnosis.
Thirdly, we have described in detail the histopathological features of XLAG adenomas. Most XLAG patients develop mixed somatotroph/lactotroph adenomas that show a characteristic sinusodoidal and lobular architecture and contain both DG and SG somatotroph cells. Microcalcifications and follicle-like structures are commonly observed. Mitotic activity was negligible in five of our cases and the Ki-67 labelling index never exceeded 3 %.
Adenomas observed in patients with XLAG differ from other GH-secreting adenomas. The classification of somatotroph adenomas is based on morphological, ultrastructural and immunophenotypical features, and includes the common DG and SG somatotroph adenomas, mixed somatotroph/lactotroph adenomas , as well as the less frequent mammosomatotroph , silent type III  and acidophilic stem cell adenomas . Neoplastic cells of DG somatotroph adenomas typically show acidophilic cytoplasm with intense and widespread immunoreactivity for GH. At electron microscopy, cells resemble normal somatotrophs. Other anterior pituitary hormones are often expressed, including PRL, the common α-subunit and β-subunits of FSH, LH and TSH in various combinations. Prolactin immunoreactivity is by far the commonest in DG adenomas and, in fact, about 50 % of acromegaly patients present with signs and symptoms of hyperprolactinemia . SG somatotroph adenomas are composed of smaller cells, have chromophobic cytoplasm with faint and focal immunoreactivity for GH, reflecting the paucity of secretory granules seen at electron microscopy. SG adenomas contain distinctive paranuclear cytokeratin-positive inclusions known as fibrous bodies. Immunoreactivity for other pituitary hormones is considerably less common than DG adenomas and, when present, their expression is weak. In about 25 % of somatotroph adenomas, tumor cells show intermediate features between DG and SG cells , a feature we did not observe in XLAG. In mixed somatotroph/lactotroph adenomas, the degree of hormone expression and distribution of GH and PRL-positive cells vary among cases. Mixed adenomas contain either DG or SG somatotrophs admixed with lactotroph cells [17, 22]; the coexistence of DG and SG somatotroph cells observed in XLAG-related adenomas has never been formally described to the best of our knowledge.
The distinction between nodular hyperplasia and adenoma is certainly challenging, given the lobular, acinar and sinusoidal/cordonal architecture of XLAG-related adenomas. Several features are, however, in keeping with a neoplastic rather than a hyperplastic process. Firstly, reticulin fibers are almost exclusively limited to perivascular spaces, which appeared thickened, distorted and stretched. Secondly, we are not aware of any description of somatotroph hyperplasia containing SG somatotroph cells [23, 24]. In addition, proliferative activity identified with Ki-67, and the presence of a few mitoses, are further features indicative of an adenoma rather than hyperplasia.
It is important to consider that GH excess was caused by mammosomatroph hyperplasia in two of our XLAG patients [13, 14], while one hyperplasia patient had pure GH excess (case II). Interestingly, none of the adenomas examined in this study demonstrated a transition from hyperplasia to adenoma or showed any tumor cells with GH and PRL co-localization. Sampling may explain the lack of hyperplasia, although we cannot exclude that mammosomatotroph hyperplasia and mixed somatotroph/lactotroph adenomas represent distinct pathological manifestations of XLAG.
The mechanisms underlying the molecular pathogenesis of XLAG are unknown. GPR101 encodes an orphan GPCR . In mice, GPR101 is expressed throughout the central nervous system, with higher levels observed in the hypothalamus [1, 25]. The finding of increased circulating GHRH levels in some patients with XLAG [1, 2, 26] and in other cases with a phenotype highly suggestive of XLAG [27, 28], suggests that upregulation of hypothalamic GHRH might play a pathogenic role. Notably, adenomatous transformation to a SG somatotroph adenoma was reported in a case of pituitary hyperplasia and ectopic acromegaly due to a pituitary metastasis from a GHRH-secreting neuroendocrine tumor . Mammosomatotroph cells account for a sizeable proportion of the human fetal anterior pituitary and are considered precursors of GH and PRL cells in the adult anterior pituitary [30, 31]. These bi-hormonal cells are frequently found in gigantism , while they are rare in adult acromegaly patients . The expansion of mammosomatotrophs observed in XLAG patients with pituitary hyperplasia could potentially result from prenatal exposure to increased GHRH levels. Interestingly, mice transgenic for GHRH develop mammosomatotroph hyperplasia , and occurrence of somatotroph adenomas has been described in older animals . Moreover, the pituitary cells from one of our patients were found to normally respond to octreotide and bromocriptine treatment in vitro, despite only a marginal biochemical response in vivo , further supporting the role for a stimulatory factor hampering responsiveness to medical treatment. Duplication of GPR101 might potentially affect the GH axis both at the pituitary and the hypothalamic level. In fact, GPR101 is coupled to the stimulatory G protein  and, in GH3 mammosomatotroph cells, overexpression of GPR101 has been found to increase cAMP levels , which represents a key factor involved in the regulation of GH secretion and cell proliferation in response to GHRH [35, 36].
Patients with XLAG have a unique clinical phenotype, whether it is a result of a germline or somatic GPR101 duplication. The disease presents very early; excessive growth started before four years of age and was the presenting feature in all patients. We confirmed that XLAG is more common in females, occurs earlier and is more frequently associated with hyperprolactinemia, compared with gigantism due to AIP mutations. Pituitary adenomas in AIP mutation carriers are frequently large and invasive, invariably SG and often characterized by a high proliferation index [3, 5, 37]. Considering the marginal overlap we observed between XLAG and AIPpos patients in regards to the age at disease onset, in the rare case of young-onset gigantism due to an AIP mutation , neuroimaging and pathological features can help in the differential diagnosis. Interestingly, two of our youngest AIPpos cases (aged 4 and 6) presented with pituitary apoplexy [8, 38], which was not observed among XLAG patients. Moreover, although described in the setting of AIP-related pituitary disease , hyperplasia is an unusual finding in AIPpos cases, occurring only in one patient in our series .
XLAG patients frequently require multi-modal treatment. First generation SSAs are usually ineffective, while dopamine agonists effectively control PRL excess in the cases where appropriate doses are used, although they don’t seem to have a significant effect on GH and IGF-1 levels. A combination of surgery and radiation was necessary in most patients, and, in some of them, allowed the control of GH excess with further use of SSAs. Pegvisomant was an effective treatment both (i) in patients whose disease was not controlled despite multi-modal treatment, and (ii) in patients with pituitary hyperplasia, where extensive pituitary surgery was not considered due to the risk of hypopituitarism.
Finally, we investigated the prevalence of GPR101 variants in a large series of acromegaly patients. The frequency of the previously reported c.924G > C (p.E308D) variant in our cohort was found to be similar to that reported in the ExAC database. No other rare or novel coding variants were identified, either at the germline or somatic level, suggesting that GPR101 variants do not occur frequently and might not play a significant role in the pathogenesis of acromegaly. These results are in line with recently published studies [10, 40].
XLAG accounts for a significant proportion of patients with non-syndromic gigantism, and results from either germline or somatic duplications involving GPR101, which we have proved to be the causative gene within the Xq26.3 region. Our data suggest that specific testing for GPR101 duplication and, in some cases, testing of alternative tissue-derived DNA samples, should be employed for the genetic diagnosis of XLAG. Most XLAG patients develop pituitary adenomas showing remarkably similar histopathological features, including sinusoidal and lobular architecture, the presence of both DG and SG somatotrophs admixed with lactotroph cells, follicle-like structures and calcifications. These features, together with the clinical phenotype, should raise the suspicion of XLAG. Further studies are needed to untangle the molecular mechanisms involved in the pathogenesis of this condition.
XLAG, X-linked acrogigantism; GH, growth hormone; GPCR, G protein-coupled receptor; AIP, aryl hydrocarbon receptor-interacting protein; FIPA, Familial Isolated Pituitary Adenoma; CNV, copy number variation; MEN1, Multiple Endocrine Neoplasia type 1; ddPCR, droplet digital PCR; aCGH, array comparative genomic hybridization; HD-aCGH, high-density array comparative genomic hybridization; PRL, prolactin; TRH, thyrotropin-relasing hormone; GHRH, growth hormone-releasing hormone; DG, densely granulated; SG, sparsely granulated; SSTR2a, somatostatin receptor 2a; SSTR5, somatostatin receptor 5; SSA, somatostatin analogue; DA, dopamine agonist.
We gratefully acknowledge patients, their family members and all healthcare professionals involved in the study (http://www.fipapatients.org/fipaconsortium). We acknowledge Dr Susan Van Noorden for the help in retrieving one of our patients’ data. We also acknowledge Prof Kalman Kovacs, Toronto, Canada for advice on the pathological findings and Prof Michael Besser, London, UK, for the help in the management of one of the XLAG patients. Prof John Jane, Charlottesville, Virginia, US, is gratefully acknowledged for retrieving the MRI images of one patient.
Grant support from the Medical Research Council of the UK (MRC), Wellcome Trust, National Institute of Health Research (NIHR), Barts and The London Charity, Royal Society, Fundación Alfonso Martín Escudero, DST-UKIERI, Wales Gene Park and Pfizer Ltd is gratefully acknowledged. This work was also supported in part by the Intramural Research Program of the Eunice Kennedy Shriver National Institute of Child Health & Human Development, Bethesda, MD (to GT and CAS), the US National Institute of Neurological Disorders and Stroke (http://www.ninds.nih.gov/, R01NS058529 to JRL), by the US National Human Genome Research Institute (NHGRI)/National Heart Lung and Blood Institute (NHLBI) (grant number U54HG006542 to the Baylor-Hopkins Center for Mendelian Genomics, BHCMG) and by the National Institute for Health Research (NIHR) Oxford Biomedical Research Centre based at Oxford University Hospitals NHS Trust and University of Oxford. We also acknowledge the Oxford Brain Bank, supported by the Medical Research Council (MRC) and Brains for Dementia Research (BDR).
The trial of Genetics of Endocrine Tumors-Familial Isolated Pituitary Adenoma-FIPA is registered on clinicaltrials.gov with the identifier of NCT00461188.
The funding bodies had no role in the study design, collection, analysis and interpretation of the data or in the preparation of the manuscript.
DI collected and entered the clinical data in the database, analyzed the data, performed most of the stainings, scored them, contributed to the sequencing for GPR101, performed the statistical analyses and prepared the manuscript. RC and BB set up the CNV ddPCR, undertook the testing and analyzed the data. BB also performed the sequencing for GPR101 for most samples. SJ performed the standard aCGH and analyzed the data. BY performed the HD-aCGH and breakpoint recombination analyses. LCHR and FC collected and entered the clinical and genetic data in the database. JE reviewed imaging investigations. MND and PG collected clinical data, managed the database and the ethical permissions, contacted the collaborators and extracted DNA samples. SK performed histological stainings and sequencing for GPR101 variants in somatotroph adenomas. FF, SJL, ABG and OA contributed to the sequencing for GPR101. CR, MLV, CRR, MM, ABG, MB, CPB, AG, PD, CSC, TC and MBL provided clinical data and samples from XLAG patients. MBL also reviewed pathology and EM features of one of the cases. GT and CAS assisted in the setup of the HD-aCGH. JT undertook the stainings for GHRH and reviewed the pathology. JRL participated in the study setup. SE participated in the study design and her lab performed the CNV ddPCR. JRS participated in the setup of the study and his lab performed the standard aCGH. FR scored the stainings, reviewed the pathology and prepared the manuscript. MK designed and coordinated the study, recruited the patients, cared for some of the patients, collected and entered the clinical and genetic data in the database and prepared the manuscript. All authors reviewed the manuscript. All authors read and approved the final manuscript.
ABG has received speaker fees from Novartis and Ipsen. JRL holds stock ownership in 23andMe and Lasergen, is a paid consultant for Regeneron Pharmaceuticals, and is a co-inventor of multiple United States and European patents related to molecular diagnostics. The Department of Molecular and Human Genetics at Baylor College of Medicine derives revenue from molecular genetic testing offered in the Baylor Miraca Genetics Laboratories (BMGL; http://bmgl.com/). JRL is a member of the Scientific Advisory Board of BMGL. SE is a Wellcome Trust Senior Investigator. MK is a member of the Scientific Advisory Board of Pfizer and received grant support from Pfizer and Novartis. The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Trivellin G, Daly AF, Faucz FR, Yuan B, Rostomyan L, Larco DO, Schernthaner-Reiter MH, Szarek E, Leal LF, Caberg JH, Castermans E, Villa C, Dimopoulos A, Chittiboina P, Xekouki P, Shah N, Metzger D, Lysy PA, Ferrante E, Strebkova N, Mazerkina N, Zatelli MC, Lodish M, Horvath A, de Alexandre RB, Manning AD, Levy I, Keil MF, Sierra Mde L, Palmeira L, Coppieters W, Georges M, Naves LA, Jamar M, Bours V, Wu TJ, Choong CS, Bertherat J, Chanson P, Kamenicky P, Farrell WE, Barlier A, Quezado M, Bjelobaba I, Stojilkovic SS, Wess J, Costanzi S, Liu P, Lupski JR, Beckers A, Stratakis CA. Gigantism and acromegaly due to Xq26 microduplications and GPR101 mutation. N Engl J Med. 2014;371:2363–74. doi:10.1056/NEJMoa1408028.View ArticlePubMedPubMed CentralGoogle Scholar
- Beckers A, Lodish MB, Trivellin G, Rostomyan L, Lee M, Faucz FR, Yuan B, Choong CS, Caberg JH, Verrua E, Naves LA, Cheetham TD, Young J, Lysy PA, Petrossians P, Cotterill A, Shah NS, Metzger D, Castermans E, Ambrosio MR, Villa C, Strebkova N, Mazerkina N, Gaillard S, Barra GB, Casulari LA, Neggers SJ, Salvatori R, Jaffrain-Rea ML, Zacharin M, Santamaria BL, Zacharieva S, Lim EM, Mantovani G, Zatelli MC, Collins MT, Bonneville JF, Quezado M, Chittiboina P, Oldfield EH, Bours V, Liu P, W WdH, Pellegata N, Lupski JR, Daly AF, Stratakis CA. X-linked acrogigantism syndrome: clinical profile and therapeutic responses. Endocr Relat Cancer. 2015;22:353–67. doi:10.1530/ERC-15-0038.View ArticlePubMedPubMed CentralGoogle Scholar
- Hérnandez-Ramírez LC, Gabrovska P, Denes J, Stals K, Trivellin G, Tilley D, Ferraù F, Evanson J, Ellard S, Grossman AB, Roncaroli F, Gadelha MR, Korbonits M, International FC. Landscape of familial isolated and young-onset pituitary adenomas: prospective diagnosis in AIP mutation carriers. J Clin Endocrinol Metab. 2015;100:E1242–54. doi:10.1210/jc.2015-1869.View ArticlePubMedPubMed CentralGoogle Scholar
- Cuny T, Pertuit M, Sahnoun-Fathallah M, Daly A, Occhi G, Odou MF, Tabarin A, Nunes ML, Delemer B, Rohmer V, Desailloud R, Kerlan V, Chabre O, Sadoul JL, Cogne M, Caron P, Cortet-Rudelli C, Lienhardt A, Raingeard I, Guedj AM, Brue T, Beckers A, Weryha G, Enjalbert A, Barlier A. Genetic analysis in young patients with sporadic pituitary macroadenomas: besides AIP don't forget MEN1 genetic analysis. Eur J Endocrinol. 2013;168:533–41. doi:10.1530/EJE-12-0763.View ArticlePubMedGoogle Scholar
- Rostomyan L, Daly AF, Petrossians P, Nachev E, Lila AR, Lecoq AL, Lecumberri B, Trivellin G, Salvatori R, Moraitis AG, Holdaway I, Kranenburg-van Klaveren DJ, Chiara Zatelli M, Palacios N, Nozieres C, Zacharin M, Ebeling T, Ojaniemi M, Rozhinskaya L, Verrua E, Jaffrain-Rea ML, Filipponi S, Gusakova D, Pronin V, Bertherat J, Belaya Z, Ilovayskaya I, Sahnoun-Fathallah M, Sievers C, Stalla GK, Castermans E, Caberg JH, Sorkina E, Auriemma RS, Mittal S, Kareva M, Lysy PA, Emy P, De Menis E, Choong CS, Mantovani G, Bours V, De Herder W, Brue T, Barlier A, Neggers SJ, Zacharieva S, Chanson P, Shah NS, Stratakis CA, Naves LA, Beckers A. Clinical and genetic characterization of pituitary gigantism: an international collaborative study in 208 patients. Endocr Relat Cancer. 2015;22:745–57. doi:10.1530/ERC-15-0320.View ArticlePubMedGoogle Scholar
- Melmed S, Kleinberg D. Pituitary masses and tumors, Williams textbook of endocrinology. 12th ed. Philadelphia: Elsevier Saunders; 2011.Google Scholar
- Leontiou CA, Gueorguiev M, van der Spuy J, Quinton R, Lolli F, Hassan S, Chahal HS, Igreja SC, Jordan S, Rowe J, Stolbrink M, Christian HC, Wray J, Bishop-Bailey D, Berney DM, Wass JA, Popovic V, Ribeiro-Oliveira A, Jr., Gadelha MR, Monson JP, Akker SA, Davis JR, Clayton RN, Yoshimoto K, Iwata T, Matsuno A, Eguchi K, Musat M, Flanagan D, Peters G, Bolger GB, Chapple JP, Frohman LA, Grossman AB, Korbonits M. The role of the aryl hydrocarbon receptor-interacting protein gene in familial and sporadic pituitary adenomas. J Clin Endocrinol Metab. 2008;93:2390–401. doi:10.1210/jc.2007-2611.View ArticlePubMedGoogle Scholar
- Igreja S, Chahal HS, King P, Bolger GB, Srirangalingam U, Guasti L, Chapple JP, Trivellin G, Gueorguiev M, Guegan K, Stals K, Khoo B, Kumar AV, Ellard S, Grossman AB, Korbonits M, International FC. Characterization of aryl hydrocarbon receptor interacting protein (AIP) mutations in familial isolated pituitary adenoma families. Hum Mutat. 2010;31:950–60. doi:10.1002/humu.21292.View ArticlePubMedPubMed CentralGoogle Scholar
- Kamenicky P, Bouligand J, Chanson P. Gigantism, acromegaly, and GPR101 mutations. N Engl J Med. 2015;372:1264. doi:10.1056/NEJMc1500340#SA1.View ArticlePubMedGoogle Scholar
- Lecoq AL, Bouligand J, Hage M, Cazabat L, Salenave S, Linglart A, Young J, Guiochon-Mantel A, Chanson P, Kamenicky P. Very low frequency of germline GPR101 genetic variation and no biallelic defects with AIP in a large cohort of patients with sporadic pituitary adenomas. Eur J Endocrinol. 2016. doi:10.1530/EJE-15-1044.PubMedGoogle Scholar
- Volante M, Brizzi MP, Faggiano A, La Rosa S, Rapa I, Ferrero A, Mansueto G, Righi L, Garancini S, Capella C, De Rosa G, Dogliotti L, Colao A, Papotti M. Somatostatin receptor type 2A immunohistochemistry in neuroendocrine tumors: a proposal of scoring system correlated with somatostatin receptor scintigraphy. Mod Pathol. 2007;20:1172–82. doi:10.1038/modpathol.3800954.View ArticlePubMedGoogle Scholar
- Iacovazzo D, Carlsen E, Lugli F, Chiloiro S, Piacentini S, Bianchi A, Giampietro A, Mormando M, Clear AJ, Doglietto F, Anile C, Maira G, Lauriola L, Rindi G, Roncaroli F, Pontecorvi A, Korbonits M, De Marinis L. Factors predicting pasireotide responsiveness in somatotroph pituitary adenomas resistant to first-generation somatostatin analogues: an immunohistochemical study. Eur J Endocrinol. 2016;174:241–50. doi:10.1530/EJE-15-0832.View ArticlePubMedGoogle Scholar
- Rodd C, Millette M, Iacovazzo D, Stiles C, Barry S, Evanson J, Albrecht S, Caswell R, Bunce B, Jose S, Trouillas J, Roncaroli F, Sampson J, Ellard S, Korbonits M. Somatic GPR101 duplication causing X-linked acrogigantism (X-LAG) - diagnosis and management. J Clin Endocrinol Metab. 2016. doi:10.1210/jc.2015-4366.PubMedPubMed CentralGoogle Scholar
- Moran A, Asa SL, Kovacs K, Horvath E, Singer W, Sagman U, Reubi JC, Wilson CB, Larson R, Pescovitz OH. Gigantism due to pituitary mammosomatotroph hyperplasia. N Engl J Med. 1990;323:322–7. doi:10.1056/NEJM199008023230507.View ArticlePubMedGoogle Scholar
- Naves LA, Daly AF, Dias LA, Yuan B, Zakir JC, Barra GB, Palmeira L, Villa C, Trivellin G, Junior AJ, Neto FF, Liu P, Pellegata NS, Stratakis CA, Lupski JR, Beckers A. Aggressive tumor growth and clinical evolution in a patient with X-linked acro-gigantism syndrome. Endocrine. 2015;51:236–44. doi:10.1007/s12020-015-0804-6.View ArticlePubMedGoogle Scholar
- Daly AF, Yuan B, Fina F, Caberg JH, Trivellin G, Rostomyan L, de Herder WW, Naves LA, Metzger D, Cuny T, Rabl W, Shah N, Jaffrain-Rea ML, Zatelli MC, Faucz FR, Castermans E, Nanni-Metellus I, Lodish M, Muhammad A, Palmeira L, Potorac I, Mantovani G, Neggers SJ, Klein M, Barlier A, Liu P, Ouafik L, Bours V, Lupski JR, Stratakis CA, Beckers A. Somatic mosaicism underlies X-linked acrogigantism syndrome in sporadic male subjects. Endocr Relat Cancer. 2016;23:221–33. doi:10.1530/ERC-16-0082.View ArticlePubMedGoogle Scholar
- Horvath E, Kovacs K. Pathology of acromegaly. Neuroendocrinology. 2006;83:161–5. doi:10.1159/000095524.View ArticlePubMedGoogle Scholar
- Erickson D, Scheithauer B, Atkinson J, Horvath E, Kovacs K, Lloyd RV, Young WF, Jr. Silent subtype 3 pituitary adenoma: a clinicopathologic analysis of the Mayo Clinic experience. Clin Endocrinol. 2009;71:92–9. doi:10.1111/j.1365-2265.2008.03514.x.View ArticleGoogle Scholar
- Horvath E, Kovacs K, Singer W, Smyth HS, Killinger DW, Erzin C, Weiss MH. Acidophil stem cell adenoma of the human pituitary: clinicopathologic analysis of 15 cases. Cancer. 1981;47:761–71.View ArticlePubMedGoogle Scholar
- Lopes MB. Growth hormone-secreting adenomas: pathology and cell biology. Neurosurg Focus. 2010;29, E2. doi:10.3171/2010.7.FOCUS10169.View ArticlePubMedGoogle Scholar
- Obari A, Sano T, Ohyama K, Kudo E, Qian ZR, Yoneda A, Rayhan N, Mustafizur Rahman M, Yamada S. Clinicopathological features of growth hormone-producing pituitary adenomas: difference among various types defined by cytokeratin distribution pattern including a transitional form. Endocr Pathol. 2008;19:82–91. doi:10.1007/s12022-008-9029-z.View ArticlePubMedGoogle Scholar
- Corenblum B, Sirek AM, Horvath E, Kovacs K, Ezrin C. Human mixed somatotrophic and lactotrophic pituitary adenomas. J Clin Endocrinol Metab. 1976;42:857–63. doi:10.1210/jcem-42-5-857.View ArticlePubMedGoogle Scholar
- Al-Gahtany M, Horvath E, Kovacs K. Pituitary hyperplasia. Hormones. 2003;2:149–58.View ArticlePubMedGoogle Scholar
- Horvath E, Kovacs K, Scheithauer BW. Pituitary hyperplasia. Pituitary. 1999;1:169–79.View ArticlePubMedGoogle Scholar
- Bates B, Zhang L, Nawoschik S, Kodangattil S, Tseng E, Kopsco D, Kramer A, Shan Q, Taylor N, Johnson J, Sun Y, Chen HM, Blatcher M, Paulsen JE, Pausch MH. Characterization of Gpr101 expression and G-protein coupling selectivity. Brain Res. 2006;1087:1–14. doi:10.1016/j.brainres.2006.02.123.View ArticlePubMedGoogle Scholar
- Daly AF, Lysy P, Defilles C, Rostomyan L, Mohamed A, Caberg JH, Raverot V, Castermans E, Marbaix E, Maiter D, Brunelle C, Trivellin G, Stratakis CA, Bours V, Raftopoulos C, Beauloye V, Barlier A, Beckers A. Growth hormone releasing hormone excess and blockade in X-LAG syndrome. Endocr Relat Cancer. 2015;23:161–70. doi:10.1530/ERC-15-0478.View ArticlePubMedGoogle Scholar
- Zimmerman D, Young Jr WF, Ebersold MJ, Scheithauer BW, Kovacs K, Horvath E, Whitaker MD, Eberhardt NL, Downs TR, Frohman LA. Congenital gigantism due to growth hormone-releasing hormone excess and pituitary hyperplasia with adenomatous transformation. J Clin Endocrinol Metab. 1993;76:216–22. doi:10.1210/jcem.76.1.8421089.PubMedGoogle Scholar
- Dubuis JM, Deal CL, Drews RT, Goodyer CG, Lagace G, Asa SL, Van Vliet G, Collu R. Mammosomatotroph adenoma causing gigantism in an 8-year old boy: a possible pathogenetic mechanism. Clin Endocrinol. 1995;42:539–49.View ArticleGoogle Scholar
- Nasr C, Mason A, Mayberg M, Staugaitis SM, Asa SL. Acromegaly and somatotroph hyperplasia with adenomatous transformation due to pituitary metastasis of a growth hormone-releasing hormone-secreting pulmonary endocrine carcinoma. J Clin Endocrinol Metab. 2006;91:4776–80. doi:10.1210/jc.2006-0610.View ArticlePubMedGoogle Scholar
- Mulchahey JJ, Jaffe RB. Detection of a potential progenitor cell in the human fetal pituitary that secretes both growth hormone and prolactin. J Clin Endocrinol Metab. 1988;66:24–32. doi:10.1210/jcem-66-1-24.View ArticlePubMedGoogle Scholar
- Asa SL, Kovacs K, Horvath E, Losinski NE, Laszlo FA, Domokos I, Halliday WC. Human fetal adenohypophysis. Electron microscopic and ultrastructural immunocytochemical analysis. Neuroendocrinology. 1988;48:423–31.View ArticlePubMedGoogle Scholar
- Scheithauer BW, Kovacs KT, Stefaneanu L, Horvath E, Kane LA, Young Jr WF, Lloyd RV, Randall RV, Davis DH. The pituitary in gigantism. Endocr Pathol. 1995;6:173–87.View ArticlePubMedGoogle Scholar
- Stefaneanu L, Kovacs K, Horvath E, Asa SL, Losinski NE, Billestrup N, Price J, Vale W. Adenohypophysial changes in mice transgenic for human growth hormone-releasing factor: a histological, immunocytochemical, and electron microscopic investigation. Endocrinology. 1989;125:2710–8. doi:10.1210/endo-125-5-2710.View ArticlePubMedGoogle Scholar
- Asa SL, Kovacs K, Stefaneanu L, Horvath E, Billestrup N, Gonzalez-Manchon C, Vale W. Pituitary adenomas in mice transgenic for growth hormone-releasing hormone. Endocrinology. 1992;131:2083–9. doi:10.1210/endo.131.5.1425411.PubMedGoogle Scholar
- Struthers RS, Vale WW, Arias C, Sawchenko PE, Montminy MR. Somatotroph hypoplasia and dwarfism in transgenic mice expressing a non-phosphorylatable CREB mutant. Nature. 1991;350:622–4. doi:10.1038/350622a0.View ArticlePubMedGoogle Scholar
- Billestrup N, Swanson LW, Vale W. Growth hormone-releasing factor stimulates proliferation of somatotrophs in vitro. Proc Natl Acad Sci U S A. 1986;83:6854–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Daly AF, Tichomirowa MA, Petrossians P, Heliovaara E, Jaffrain-Rea ML, Barlier A, Naves LA, Ebeling T, Karhu A, Raappana A, Cazabat L, De Menis E, Montanana CF, Raverot G, Weil RJ, Sane T, Maiter D, Neggers S, Yaneva M, Tabarin A, Verrua E, Eloranta E, Murat A, Vierimaa O, Salmela PI, Emy P, Toledo RA, Sabate MI, Villa C, Popelier M, Salvatori R, Jennings J, Longas AF, Labarta Aizpun JI, Georgitsi M, Paschke R, Ronchi C, Valimaki M, Saloranta C, De Herder W, Cozzi R, Guitelman M, Magri F, Lagonigro MS, Halaby G, Corman V, Hagelstein MT, Vanbellinghen JF, Barra GB, Gimenez-Roqueplo AP, Cameron FJ, Borson-Chazot F, Holdaway I, Toledo SP, Stalla GK, Spada A, Zacharieva S, Bertherat J, Brue T, Bours V, Chanson P, Aaltonen LA, Beckers A. Clinical characteristics and therapeutic responses in patients with germ-line AIP mutations and pituitary adenomas: an international collaborative study. J Clin Endocrinol Metab. 2010;95:E373–83. doi:10.1210/jc.2009-2556.View ArticlePubMedGoogle Scholar
- Korbonits M, Dutta P, Reddy KS, Bhansali A, Gupta P, Rai A, Radotra BD, Mukherjee KK. Exome sequencing reveals double hit by AIP gene mutation and copy loss of chromosome 11 but negative X-LAG in a pituitary adenoma of a 4 year old child with gigantism treated with multimodal therapy, Abstract presented at the Endocrine Society's 98th Annual Meeting and Expo. 2016.Google Scholar
- Villa C, Lagonigro MS, Magri F, Koziak M, Jaffrain-Rea ML, Brauner R, Bouligand J, Junier MP, Di Rocco F, Sainte-Rose C, Beckers A, Roux FX, Daly AF, Chiovato L. Hyperplasia-adenoma sequence in pituitary tumorigenesis related to aryl hydrocarbon receptor interacting protein gene mutation. Endocr Relat Cancer. 2011;18:347–56. doi:10.1530/ERC-11-0059.View ArticlePubMedGoogle Scholar
- Ferraù F, Romeo PD, Puglisi S, Ragonese M, Torre ML, Scaroni C, Occhi G, De Menis E, Arnaldi G, Trimarchi F, Cannavo S. Analysis of GPR101 and AIP genes mutations in acromegaly: a multicentric study. Endocrine. 2016. doi:10.1007/s12020-016-0862-4.Google Scholar