- Research article
- Open Access
The genomic profiling and MAMLD1 expression in human and canines with Cushing’s disease
BMC Endocrine Disorders volume 21, Article number: 185 (2021)
Cushing’s disease (CD) is defined as hypercortisolemia caused by adrenocorticotropic hormone (ACTH)-secreting pituitary adenomas (corticotroph PA) that afflicts humans and dogs. In order to map common aberrant genomic features of CD between humans and dogs, we performed genomic sequencing and immunostaining on corticotroph PA.
For inclusion, humans and dog were diagnosed with CD. Whole exome sequencing (WES) was conducted on 6 human corticotroph PA. Transcriptome RNA-Seq was performed on 6 human and 7 dog corticotroph PA. Immunohistochemistry (IHC) was complete on 31 human corticotroph PA. Corticotroph PA were compared with normal tissue and between species analysis were also performed.
Eight genes (MAMLD1, MNX1, RASEF, TBX19, BIRC5, TK1, GLDC, FAM131B) were significantly (P < 0.05) overexpressed across human and canine corticotroph PA. IHC revealed MAMLD1 to be positively (3+) expressed in the nucleus of ACTH-secreting tumor cells of human corticotroph PA (22/31, 70.9%), but absent in healthy human pituitary glands.
In this small exploratory cohort, we provide the first preliminary insights into profiling the genomic characterizations of human and dog corticotroph PA with respect to MAMLD1 overexpression, a finding of potential direct impact to CD microadenoma diagnosis. Our study also offers a rationale for potential use of the canine model in development of precision therapeutics.
Corticotrope pituitary adenomas (PA) are benign tumors of the anterior pituitary gland that secrete excessive amounts of adrenocorticotropin hormone (ACTH), causing Cushing’s disease (CD) . The precursor polyprotein to ACTH, proopiomelanocortin (POMC), is synthesized within residential anterior pituitary corticotropic cells . When ACTH is secreted, zona fasciculata cells within the adrenal cortex release cortisol, which exerts diverse physiological systemic effects [3, 4]. Mortality is increased up to 5 times in CD compared with the general population, but significantly improves after achievement of normocortisolemia by treatment [5, 6].
CD treatment options are diverse. Surgery is the first choice treatment and removal of the tumor results in remission in 75–90% of patients with a 35% biochemical recurrence [4, 7,8,9,10]. The remainder have either small tumors not identified at surgery or nonresectable tumors that invade the cavernous sinus . Other options include radiotherapy and bilateral adrenalectomy, however these procedures have risks of permanent cortisol and aldosterone deficiency along with potentiation of aggressive tumor growth and inevasible behavior (Corticotroph Tumor Progression) . Finally, while medical treatment options have expanded to include somatostatin receptor ligands, cortisol synthesis inhibitors, glucocorticoid receptor blockers, and dopamine agonists, long-term safety data continues to be lacking [12,13,14,15]. Taken together, many of the available compounds target the adrenal gland or suppress ACTH secretion, but not the source of the secreting PA, therefore there is a strong need for further research on this topic.
Because CD is a rare tumor (incidence of 4:1,000,000 [16, 17]), research has yielded a limited number of dysfunctional genes for diagnostic or treatment purposes. While genes (POMC, t-box transcription factor 19 (TBX19), ubiquitin carboxyl-terminal hydrolase 8 (USP8), ubiquitin carboxyl-terminal hydrolase 48 (USP48) and cyclin-dependent kinase 5 and abelson murine leukemia viral oncogene homolog 1 enzyme substrate 1 (CABLES1)) associated to PA have been identified [2, 18,19,20,21,22,23,24,25,26,27], preclinical studies suggest that only USP8 PA mutation and TBX19 expression may respond to epidermal growth factor receptor (EGFR) inhibition and cyclin-dependent kinase (CDK) inhibitors, respectively [18,19,20,21]. With limited results, a broader catalogue of PA mutations are needed to explore a wider variety of treatment diagnosis options.
Dogs may be viable models for targeted interventions in CD. Dogs share an 84% genetic overlap with humans, and they naturally develop Cushing’s disease from the intermediate pituitary lobe with an incidence of 1500:1,000,000 and recurrence after hypophysectomy of 27%, analogous to humans [28,29,30,31,32]. A canine disease model has advantages over transgenic mouse models, particularly since dogs share similar environmental exposures and stressors to humans . However, while TBX19 expression were demonstrated in dog PA, USP8 mutations were not identified as frequently mutant . This prompted us to investigate a more complete catalogue of genomic aberrations that may lead to dysfunctional protein expression observed in CD and could potentially lead to the development of targeted pharmaceutical interventions.
Using both human and dog corticotroph PA, our objective was to describe the common aberrant genomic features of CD between humans and dogs. We profiled genetic expression using whole transcriptome sequencing of both species. Whole exome sequencing (WES) was also performed on the human corticotroph PA and an identified target was confirmed by immunohistochemistry when compared to normal human pituitary tissue. Across experimentation, corticotroph PA were compared with normal tissue, and when possible between species analysis were also performed. Furthermore, although genetic studies have been previously conducted in humans with CD [15, 18, 19, 24, 34,35,36,37], our results provide a broad overview of the genetic profiles between human and dog corticotroph PA.
For eligibility, human and dogs were required to have a CD diagnosis. Genomic sequencing and immunostaining were performed on corticotroph PA. Across experimentation, corticotroph PA were compared with normal tissue, and when possible between species analysis were also performed.
Thirty-seven CD African American patients (average 43 +/− 12 years old) with intermediate anterior sized PA (avg lateral 1.16 +/− 0.44 cm x AP 0.88 +/− 0.29 cm x cranial 0.90 +/− 0.49 cm) were recruited through the Department of Neurosurgery at Emory University Hospital. In all cases, the diagnosis of Cushing’s disease entailed the following steps: 1. two or more abnormal screening tests (including lack of suppression to low-dose dexamethasone, higher than normal urinary free cortisol levels and/or higher than normal late-night salivary cortisol levels), 2. ACTH-dependent hypercortisolism (i.e. high or high-normal ACTH level) and 3. localization tests supporting autonomous pituitary ACTH secretion (cortisol suppression by more than 50% or below 5 mcg/dL after to high-dose dexamethasone and/or ACTH stimulation after CRH test). In addition, depending on pituitary MRI findings, selected patients underwent petrosal sinus sampling which indicated a baseline central-to-periphery ACTH ratio > 2 and/or post-CRH ratio > 3). For all cases, pathology reports showed immunopositive staining for ACTH, and negative immunostaining for tumor protein p53 (TP53), follicle-stimulating hormone (FSH), thyroid stimulating hormone (TSH), luteinizing hormone (LH), growth hormone (GH), and E3 ubiquitin protein-ligase (MIB-1) proliferation index < 3%. During microsurgical transsphenoidal hypophysectomy, the PA tissue was removed with no complications.
In total, 39 human specimens were utilized (37 corticotroph PA, 2 normals). For WES and RNA-seq, frozen specimens were available from 6 African American patients (1 male, 5 females) and normal human pituitary specimens from 1 male (ND01199–09) and 1 female (ND01218–06) were purchased from the Coriell Biorepository (Camden, NJ). For IHC, paraffin specimens were available from 31 African American patients with 1 normal anterior pituitary tissue.”
CD Dogs from canine client owners were recruited through the Department of Clinical Sciences of Companion Animals at the Utrecht University, the Netherlands. Suspicion of hypercortisolism was raised based on the characteristic clinical signs including polyphagia, polydipsia/polyuria, skin atrophy, thin hair coat, calcinosis cutis, truncal obesity (pot belly), depression, and exercise intolerance. Serum chemistry showed typical changes: alkaline phosphatase, alanine aminotransferase and total cholesterol were commonly increased in canine Cushing’s disease.
Preliminary diagnosis of pituitary-dependent hypercortisolism was based on increased urinary corticoid-to-creatinine ratios (UCCRs) in the first two morning urinary samples collected at home and more than 50% suppression of the UCCR in the third urine sample in the oral high dose (0.1 mg/kg) dexamethasone suppression test . CD was further confirmed by measurement of elevated plasma ACTH concentration, visualization of symmetrically enlarged adrenal glands by ultrasonography, and visualization of pituitary gland enlargement with computed tomography (CT) or MRI. During microsurgical transsphenoidal hypophysectomy in client-owned dogs, the PA tissue was removed with no complications. Immediately after collection, specimens of pars distalis PA tissue were fixed in 4% neutral buffered formaldehyde, embedded in paraffin, and consecutive sections were used for histology and immunohistochemistry for ACTH, α-melanocyte-stimulating hormone (MSH), and GH . Representative adenoma tissue samples were snap-frozen and stored in liquid nitrogen until analysis after histology confirmed a basophilic adenoma with ACTH immunostaining.
Pituitary specimens from 6 healthy beagle dogs and 7 client-owned dogs with CD that underwent hypophysectomy were sequenced. The cohort of client-owned dog patients included different breeds (i.e., Jack Russell Terrier, Scottish Shepherd, English Springer Spaniel, Mixed breed, French Bulldog, Chesapeake Bay Retriever, and American Staffordshire Terrier), 3 female (of which two castrated) and 4 male (of which one castrated) dogs, with a median age of 8.2 years (range 5.6–10.7 years). All owners consented to the hypophysectomy as treatment for CD. As control tissue, anterior lobes of normal pituitary glands were obtained from 6 healthy Beagle dogs euthanized in unrelated experiments (all female; age range 1.6–1.8 years).
Nucleic acid extraction and sequencing
DNA and RNA were extracted from 6 intermediate sized frozen human PA tissues using E.Z.N.A. Ⓡ kits (Omega Bio-tek, Norcross, GA). RNA was extracted from the canine healthy and tumor pituitary samples using the miRCURY cell and plant RNA isolation kit according to the manufacturer’s instruction including an on-column DNAse digestion (Qiagen, RNAse free DNAse, 79,254). Specimen quantity and quality were assessed using Qubit and Agilent Bioanalyzer.
Human exome sequence analysis
Human genomic DNA was extracted from fresh-frozen PA and normal anterior pituitary as described above. Whole human exome libraries were prepared using the Agilent SureSelect Human All Exon (version 5) per manufacturer’s protocols. Libraries were paired-end sequenced at 50x coverage using and Illumina HiSeq 2000. FASTQ files were aligned to the human hg19 reference genome using BWA 0.7.5.a and duplicate reads were removed with Picard tools (Version 1.1.1) [32, 38]. De-duplicated aligned (BAM) files were used for copy number estimation using Control-FREEC with recommended settings for exome sequencing . Normalized FREEC output values were segmented with DNACopy . Copy number assessments for 6 PA samples were matched to one normal. Human PA specimen 1, had a matched normal sample and mutations were called using Varscan2 in somatic mode under standard parameters . Mutations in the other 5 samples were called using Varscan2 without a matched normal. Human PA specimen 1 was also analyzed with Mutect (version 1.1.4) using standard parameters following probabilistic indel realignment and base recalibration using GATK (version 3.3.0) . GISTIC-like analysis using cghMCR were used to identify potential regions of focal gain and loss . Identified single nucleotide variants (SNV) were compared to somatic variants in the COSMIC database (version 81) and gnomAD database [38, 44]. All predicted mutations were annotated using ANNOVAR .
RNA sequence analysis
Total RNA was processed using the Illumina TruSeq RNA kit and paired-end sequenced (2 × 75) at 100 M reads per specimen on a HiSeq 2000 instrument. FASTQ files were aligned to the hg19 human or CanFam3 reference genomes using Tophat 2.0.1 using standard parameters . RefSeq, CanFam3 and Ensembl transcripts were quantified using Cufflinks (version 2.0.1) . Human RefSeq transcripts were also quantified using HTSeq (version 0.6.1) . Mutations were called in the canine and human samples using Varscan2 as described above. Identification of differentially expressed genes (> 2 fold change) between tumor and normal tissue by simple hierarchical clustering of normalized fragments per kilobase of transcript per million mapped reads (FPKM) values were assessed separately for human and dog Cushing samples and dichotomized into over or under gene expression. With the feasible utilization of IHC for conformational protein upregulation, only shared overexpressed genes were analyzed.
Immunostaining was performed on 5 μm formalin fixed paraffin embedded (FFPE) sections of 31 corticotroph PA samples and 1 normal anterior pituitary tissue. Sections were deparaffinized in xylene and tissue was hydrated by a descending ethanol sequence. After rehydration, sections were incubated with 3% H2O2 to inactivate endogenous peroxidases and blocked with 1% BSA for 10 min. Due to limited tissue quantity, a random choice generator (SOCR) selected four of the eight cross-species significantly overexpressed genes (stained - TBX19, RASEF, MAMLD1, MNX1; non-stained - BIRC5, TK1, GLDC, FAM131B) as well as POMC, USP8, and ubiquitin carboxyl-terminal hydrolase 48 (USP48) were selected from reference studies [2, 18, 19, 49]. In total, seven genes were selected for staining. Antibodies against TBX19 (Sigma Aldrich, HPA072686, rabbit 1:2000), MAMLD1 (Sigma Aldrich, HPA003923, rabbit 1:250), MNX1 (Sigma Aldrich, ABN174, rabbit 1:500), RASEF (Sigma Aldrich, WH0158158M1, mouse 1:200), USP8 (Sigma Aldrich, HPA004869, rabbit 1:100), USP48 (Sigma Aldrich, HPA030046, rabbit 1:100), and POMC (Abcam, AB210605, rabbit 1:8000) were used. The antibodies were then incubated at 4 °C overnight and visualized with diaminobenzidine. After staining, 2 tissues were fully exhausted. Image analysis was performed under a light microscope at × 400 magnification. Immunopositivity was graded by two independent reviewers (AW, SGN) with a 0 to 3+ scoring system:
(0) no immunostaining.
(1) mild positive: weak immunostaining of less than 30% of tumor cells.
(2) positive: complete membranous staining and either uniform or weak in at least 50% of tumor cells.
(3+) strong positive: uniform intense nuclear and/or cytoplasmic staining in at least 80% of tumor cells.
Grading disagreements were reexamined and reassessed collectively by the same reviewers.
Differentially expressed transcripts were converted to normalized fragments per kilobase million (FPKM) values using tools described above [47, 48]. Average log transformed FPKM values for both humans and dog transcripts are expressed with standard deviations. Comparisons of each subset of genes between ACTH secreting PA and normal pituitary tissues across either the human or dog groups were considered statistically significant when P < 0.05 in a one-tail Student T-test (i.e. overexpression tail). Immunohistochemistry results were expressed as proportions, mean, standard deviations (SD) and median. Scoring agreements were expressed as a weighted Cohen’s Kappa coefficient. With no censored or missing data, a sensitivity analysis were performed between stained and non-stained genes. Stata/IC (v15.0; StataCorp LP, College Station, TX, USA) statistical software was used for all analyses.
Whole exome sequencing
To identify recurrent somatic mutations and copy number abnormalities in corticotroph PA, whole-exome sequencing was performed in human PA.
No specimens with mutations of USP8 were identified. Two human specimens with USP48 (NM_032236) SNV were found. Human PA specimen 1 had a USP48 c.1245 G > A, p.M415I (COSM904151) SNV and specimen 6 had a USP48 c.1243 A > G, p.M415V (COSM25000) SNV.
Two human PA (3 and 4) had guanine nucleotide-binding protein G(q) subunit alpha (GNAQ) c.286A > T, p.T96S substitutions.
Copy number analysis using the whole exome sequencing data was also performed to ascertain recurrent gains of whole chromosomes or chromosome arms (Fig. 1). Human PA specimens 2 and 5 had similar patterns with gain of 1q, 5, 7, 8q, 9, 12, 13 and 14. Specimen 1 shared gain of 7, 12 and 14 with specimens 2 and 5. Loss of 19 was shared by specimens 1, 3 and 6. No focal gains or losses enriched in PA in non-benign copy number variations were observed.
Gene expression between Cushing’s tumor and normal tissue was assessed in human and dogs.
Unsupervised hierarchical clustering showed the normal pituitary specimens clustered separately from the PA specimens for both species (Fig. 2). Heterogeneity of gene expression was found in individual PA specimens in both species, and a subset of genes that were over or under expressed relative to normal pituitary in both species organisms clustered (Fig. 2).
Differential gene expression showed 517 overexpressed and 1670 underexpressed genes in human PA, and 127 overexpressed and 953 underexpressed genes in dog PA relative to normal pituitary tissues (Fig. 3). Intersection of the human and dog gene lists showed 14 overexpressed genes and 239 under expressed genes shared between species. The 14 highly expressed genes were: POMC, arginine vasopressin receptor 1B (AVPR1B), baculoviral inhibitor of apoptosis repeat-containing 5 (BIRC5), centrosomal protein of 55 kDa (CEP55), family with sequence similarity 131 member B (FAM131B), frizzled-9 (FZD9), glycine decarboxylase (GLDC), mastermind like domain containing 1 (MAMLD1), motor neuron and pancreas homeobox 1 (MNX1), prostaglandin E2 receptor 4 (PTGER4), RAS and EF-hand domain-containing protein (RASEF), t-box transcription factor (TBX19), thymidine kinase 1 (TK1) and vasoactive intestinal peptide receptor 2 (VIPR2).
Gene level summaries were performed for the 14 highly expressed genes, comparing PA to normal pituitary (Fig. 4) (SP Fig. 1, 2 and 3). Human samples revealed significant (P < 0.05) results for 10 genes: TBX19, BIRC5, AVPR1B, RASEF, TK1, GLDC, FZD9, FAM131B, MAMLD1 and MNX1. In dogs, 10 genes: POMC, TBX19, BIRC5, RASEF, TK1, GLDC, CEP55, FAM131B, MAMLD1 and MNX1 yielded significant results (P < 0.05). Eight genes: MAMLD1, MNX1, RASEF, TBX19, BIRC5, TK1, GLDC, and FAM131B were significant (P < 0.05) in PA from both species.
To examine the association between significantly expressed genes and protein translation, immunohistochemistry was performed on four of the eight significant genes in human PA (Fig. 5).
Consistently strong (3+) staining was observed for POMC (26/31, 2.77+/− 0.56), MAMLD1 (22/31, 2.58+/− 0.72), USP8 (22/31, 2.52+/− 0.85), TBX19 (16/31, 2.23+/− 0.94), USP48 (16/31, 2.32+/− 0.83) and MNX1 (12/31, 2.10+/− 0.83) (SP Table. 1). RASEF (1.55+/− 0.93) displayed weak staining (< 3). Nuclear staining was observed for TBX19, MAMLD1, RASEF, and USP48. POMC staining was cytoplasmic. MNX1 and USP8 staining was nuclear and cytoplasmic. Specimens with high (3+) MAMLD1 expression shared strong (3+) staining in POMC (21/22, 95.5%) and USP8 (20/22, 90.9%). USP48 (13/22 59.1%), TBX19 (12/22, 54.5%), MNX1 (11/22, 50.0%), and RASEF (5/22, 22.7%) showed limited high coexpressions among high MAMLD1 specimens (SP Fig. 1). While weak (1+) TBX19 and POMC staining was observed in normal anterior pituitary, no staining for MAMDL1, MNX1, RASEF, USP8 and USP48 was found (SP Table. 1).
The strength of agreement between reviewers was almost perfect (k = 0.87). Sensitivity analysis between stained and non-stained genes outcomes were not statistically significant (P > 0.05).
Our WES provides the first broad overview of RNA-seq gene expression profiles of PA in dog and man with CD. Mutations in USP8 as the main driver and Guanine nucleotide-binding protein G(s) subunit alpha isoforms short (GNAS) have been implicated in PA tumorigenesis [18, 19, 34].
While USP8 even by manual inspection of aligned DNA and RNA sequencing reads were not found, two separate non-synonymous changes were found at codon 415 of USP48 (p.M415I and p.M415V). Both have been reported as somatic variants in recent Cushing’s PA and in the COSMIC database in cancers arising from the urinary tract and prostate . These variants were unlikely to be either common or rare population variants as they were absent from the gnomAD database. In these variants, methionine 415 is a conserved residue falling within the Peptidase C19 ubiquitinyl hydrolase domain of USP48 and bioinformatic predictions using Polyphen 2 suggested that both changes may affected protein function.
Two specimens had Guanine nucleotide-binding protein G(q) subunit alpha q polypeptide (GNAQ) p.T96S SNVs instead of GNAS mutations. T96S has been reported in 20 COSMIC samples (COSM404628) - predominantly large intestine and skin - and may represent a third hotspot within the protein additional to the well-characterized Q209 and R183 positions. There have been more recent studies that have shown that no significantly recurrent mutations are identified by WES in PA [27, 50]. Combined with our results, these discrepancy illustrates that PA, although rare, may also be a heterogeneous disease from the genetic perspective. Larger cohorts are needed to determine if mutations or copy number abnormalities can be used to stratify and target therapies for PA.
Along with TBX19, our preliminary findings show MAMLD1 was found to be strongly associated to PA. While a relationship between gonadotroph pituitary adenomas and MAMLD1 has been elucidated , little is known about the function of MAMLD1 - to the best of our knowledge this is the first report link to Cushing’s disease. The Mastermind Like (MAML) family of genes vary significantly in function. MAMLD1 has a transcriptional co-activator role in regulating NOTCH signaling [39,40,41, 52]. Recently, another MAML family gene, Mastermind Like Transcriptional Coactivator 3 (MAML3) has been found to be associated with hereditary paragangliomas (PGL), a type of neuroendocrine tumor that is usually derived from extra-adrenal chromaffin cells as opposed to the pituitary . However, despite being in the same family, MAML3 functions as a fusion protein, whereas MAMLD1 functions as an inhibitor protein. Evidence suggests MAMLD1 may bind the recombining binding protein suppressor of hairless (RBJP) repressor protein and subsequently inhibit hairy and enhancer of split gene (HES) expression . During neuronal development, HES family members play a role in corticotrophic proliferation and differentiation , which HES1 knockouts has lead to an underdeveloped hypermorphic pituitary . Our study underscores the role MAML proteins seem to play in regulating neuroendocrine function and encourages further functional investigations.
There are several limitations to this study. First, due to the rare quantity of corticotrophic PA, only four of the eight significantly overexpressed cross-species genes in RNA-seq were able to have been randomly selected and stained for IHC from archival specimens. Although a sensitivity analysis showed no difference, continued IHC studies should be performed for non-stained genes, for MAMLD1 replication in dog PA specimens, and for exploring the value of the commonly underexpressed genes. Second, while USP8 mutations where found in 40% of human PA (Caucasian and Asian populations), similar to other studies, neither of our cohorts (2.4% chance) found USP8 mutations or statistically significant DGE of POMC (due to wider expression variance) despite staining immuno-positive for well-established CD genes (POMC, TBX19) from intermediated sized PA female samples (African American population) [18,19,20,21]. Further exploration behind different populations and PA models may be of interest to fully illustrate the genetic variations of PA. Third, appropriate tissue was not available to confirm protein expression in the dog specimens and contrary to humans canine CD rise predominately from intermediate pituitary lobe as opposed to the anterior. However, 31 archival human PA provided strong evidence that co-overexpression of MAMLD1 and TBX19 occurred in approximately 70.9% of PA. Finally, although transcriptome analysis (Affymetrix arrays) on human PA has been complete [57, 58] and no functional MAMLD1 assays were complete, our current investigation was neither a discovery cohort for direct comparison nor a functional assessment study, respectively. Rather, our cohort aimed to broadly profile the somatic aberrations of PA between dog and man.
In spite of these genetic heterogeneities, differentially expressed genes were identified that were common between human and dog PA. Our study not only highlights and contributes to the growing complex understandings of genetic mechanism of CD, but also underscores the potential clinical implications of shared mechanisms between dog and man. Since MAMLD1 appears to stain exclusively for corticotrophic tumor cells as compared to TBX19 (may stain both normal and tumor corticotrophic cells), MAMLD1 staining may be helpful in diagnosing CD, particularly in microadenomas with aberrant in symptomology and potentially recurrent when tissue quantity is limited. This may improve CD clinical outcomes and help better characterize different hormone secreting pituitary tumors. Our results also suggests that, because MAMLD1 expression appears absent in normal pituitary, targeted inhibition of MAMLD1 may be a potential potent targeted strategy for inhibiting the growth of corticotroph PA that merits further functional exploration . Furthermore, dog PA may provide suitable veterinary cohorts to test novel diagnostic and therapeutic approaches that may reduce CD burden in both dog and man.
We highlight the first preliminary insights into profiling the genomic characterizations of human and dog corticotroph PA with respect to MAMLD1 overexpression, a finding of potential direct impact to CD microadenoma diagnosis. Our study also offers a rationale for potential use of the canine model in development of precision therapeutics.
Availability of data and materials
We are currently analyzing the data from a different perspective and planning a related study. Therefore, the data and material are not shared in the current state. However, the datasets used and/or analyzed during the current study are available from the corresponding author on reasonable written request. After the conclusion of the final study, all sequencing data will be made available in NCBI SRA database.
Mastermind Like Domain Containing 1
Whole Exome Sequencing
- RNA-Seq :
Whole transcriptome shotgun sequencing
Lamberts SW, de Lange SA, Stefanko SZ. Adrenocorticotropin-secreting pituitary adenomas originate from the anterior or the intermediate lobe in Cushing’s disease: differences in the regulation of hormone secretion. J Clin Endocrinol Metab. 1982;54(2):286–91 https://doi.org/10.1210/jcem-54-2-286.
Slominski A, Wortsman J, Luger T, Paus R, Solomon S. Corticotropin releasing hormone and proopiomelanocortin involvement in the cutaneous response to stress. Physiol Rev. 2000;80(3):979–1020 https://doi.org/10.1152/physrev.2000.80.3.979.
Aguilera G. Regulation of pituitary ACTH secretion during chronic stress. Front Neuroendocrinol. 1994;15(4):321–50. https://doi.org/10.1006/frne.1994.1013.
Valassi E, Biller BMK, Swearingen B, Pecori Giraldi F, Losa M, Mortini P, et al. Delayed remission after transsphenoidal surgery in patients with Cushing’s disease. J Clin Endocrinol Metab. 2010;95(2):601–10 https://doi.org/10.1210/jc.2009-1672.
Clayton RN, Raskauskiene D, Reulen RC, Jones PW. Mortality and morbidity in Cushing’s disease over 50 years in stoke-on-Trent, UK: audit and meta-analysis of literature. J Clin Endocrinol Metab. 2011;96(3):632–42 https://doi.org/10.1210/jc.2010-1942.
Valassi E, Tabarin A, Brue T, Feelders RA, Reincke M, Netea-Maier R, et al. High mortality within 90 days of diagnosis in patients with Cushing’s syndrome: results from the ERCUSYN registry. Eur J Endocrinol. 2019;181(5):461–72 https://doi.org/10.1530/EJE-19-0464.
Nieman LK, Biller BMK, Findling JW, Murad MH, Newell-Price J, Savage MO, et al. Treatment of Cushing’s syndrome: an Endocrine Society clinical practice guideline. J Clin Endocrinol Metab. 2015;100(8):2807–31 https://doi.org/10.1210/jc.2015-1818.
Mampalam TJ. Transsphenoidal microsurgery for Cushing disease: a report of 216 cases. Ann Intern Med. 1988;109(6):487–93 https://doi.org/10.7326/0003-4819-109-6-487.
Trainer PJ, Lawrie HS, Verhelst J, Howlett TA, Lowe DG, Grossman AB, et al. Transsphenoidal resection in Cushing’s disease: undetectable serum cortisol as the definition of successful treatment. Clin Endocrinol. 1993;38(1):73–8 https://doi.org/10.1111/j.1365-2265.1993.tb00975.x.
Patil CG, Prevedello DM, Lad SP, Vance ML, Thorner MO, Katznelson L, et al. Late recurrences of Cushing’s disease after initial successful transsphenoidal surgery. J Clin Endocrinol Metab. 2008;93(2):358–62 https://doi.org/10.1210/jc.2007-2013.
Guiot G, Derome P. Surgical problems of pituitary adenomas. In: Krayenbühl H, Brihaye J, Loew F, et al., editors. Advances and technical standards in neurosurgery. Vienna: Springer Vienna; 1976. p. 3–33. https://doi.org/10.1007/978-3-7091-7080-9_1.
Fleseriu M. Medical treatment of Cushing disease: new targets, new hope. Endocrinol Metab Clin N Am. 2015;44(1):51–70 https://doi.org/10.1016/j.ecl.2014.10.006.
Feelders RA, de Bruin C, Pereira AM, Romijn JA, Netea-Maier RT, Hermus AR, et al. Pasireotide alone or with cabergoline and ketoconazole in Cushing’s disease. N Engl J Med. 2010;362(19):1846–8 https://doi.org/10.1056/NEJMc1000094.
Hamrahian AH, Yuen KCJ, Hoffman AR, AACE Neuroendocrine And Pituitary Scientific Committee (2014) AACE/ACE disease state clinical review: medical Management of Cushing Disease. Endocr Pract 20:746–757. https://doi.org/10.4158/EP14147.RA.
Bilodeau S, Vallette-Kasic S, Gauthier Y, Figarella-Branger D, Brue T, Berthelet F, et al. Role of Brg1 and HDAC2 in GR trans-repression of the pituitary POMC gene and misexpression in Cushing disease. Genes Dev. 2006;20(20):2871–86 https://doi.org/10.1101/gad.1444606.
Molitch ME. Diagnosis and treatment of pituitary adenomas: a review. JAMA. 2017;317(5):516–24. https://doi.org/10.1001/jama.2016.19699.
Steffensen C, Bak AM, Rubeck KZ, Jørgensen JOL. Epidemiology of Cushing’s syndrome. Neuroendocrinology. 2010;92(Suppl 1):1–5 https://doi.org/10.1159/000314297.
Ma Z-Y, Song Z-J, Chen J-H, Wang YF, Li SQ, Zhou LF, et al. Recurrent gain-of-function USP8 mutations in Cushing’s disease. Cell Res. 2015;25(3):306–17 https://doi.org/10.1038/cr.2015.20.
Reincke M, Sbiera S, Hayakawa A, Theodoropoulou M, Osswald A, Beuschlein F, et al. Mutations in the deubiquitinase gene USP8 cause Cushing’s disease. Nat Genet. 2015;47(1):31–8 https://doi.org/10.1038/ng.3166.
Vallette-Kasic S, Figarella-Branger D, Grino M, Pulichino AM, Dufour H, Grisoli F̧, et al. Differential regulation of proopiomelanocortin and pituitary-restricted transcription factor (TPIT), a new marker of normal and adenomatous human corticotrophs. J Clin Endocrinol Metab. 2003;88(7):3050–6 https://doi.org/10.1210/jc.2002-021934.
Liu N-A, Araki T, Cuevas-Ramos D, Hong J, Ben-Shlomo A, Tone Y, et al. Cyclin E-mediated human proopiomelanocortin regulation as a therapeutic target for Cushing disease. J Clin Endocrinol Metab. 2015;100(7):2557–64 https://doi.org/10.1210/jc.2015-1606.
Sbiera S, Perez-Rivas LG, Taranets L, Weigand I, Flitsch J, Graf E, et al. Driver mutations in USP8 wild-type Cushing’s disease. Neuro-oncology. 2019;21(10):1273–83 https://doi.org/10.1093/neuonc/noz109.
Perez-Rivas LG, Theodoropoulou M, Ferraù F, Nusser C, Kawaguchi K, Stratakis CA, et al. The gene of the ubiquitin-specific protease 8 is frequently mutated in adenomas causing Cushing’s disease. J Clin Endocrinol Metab. 2015;100(7):E997–1004 https://doi.org/10.1210/jc.2015-1453.
Hernández-Ramírez LC, Gam R, Valdés N, Lodish MB, Pankratz N, Balsalobre A, et al. Loss-of-function mutations in the CABLES1 gene are a novel cause of Cushing’s disease. Endocr Relat Cancer. 2017;24(8):379–92 https://doi.org/10.1530/ERC-17-0131.
Neou M, Villa C, Armignacco R, Jouinot A, Raffin-Sanson ML, Septier A, et al. Pangenomic classification of pituitary neuroendocrine tumors. Cancer Cell. 2020;37(1):123–34 e5. https://doi.org/10.1016/j.ccell.2019.11.002.
Salomon MP, Wang X, Marzese DM, Hsu SC, Nelson N, Zhang X, et al. The Epigenomic landscape of pituitary adenomas reveals specific alterations and differentiates among acromegaly, Cushing’s disease and endocrine-inactive subtypes. Clin Cancer Res. 2018;24(17):4126–36 https://doi.org/10.1158/1078-0432.CCR-17-2206.
Bi WL, Horowitz P, Greenwald NF, Abedalthagafi M, Agarwalla PK, Gibson WJ, et al. Landscape of genomic alterations in pituitary adenomas. Clin Cancer Res. 2017;23(7):1841–51 https://doi.org/10.1158/1078-0432.CCR-16-0790.
de Bruin C, Meij BP, Kooistra HS, Hanson JM, Lamberts SWJ, Hofland LJ. Cushing’s disease in dogs and humans. Horm Res. 2009;71(Suppl 1):140–3 https://doi.org/10.1159/000178058.
Karlsson EK, Lindblad-Toh K. Leader of the pack: gene mapping in dogs and other model organisms. Nat Rev Genet. 2008;9(9):713–25. https://doi.org/10.1038/nrg2382.
Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45. https://doi.org/10.1101/gr.092759.109.
Lindblad-Toh K, Wade CM, Mikkelsen TS, et al. Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature. 2005;438(7069):803–19. https://doi.org/10.1038/nature04338.
van Rijn SJ, Galac S, Tryfonidou MA, et al. The influence of pituitary size on outcome after Transsphenoidal Hypophysectomy in a large cohort of dogs with pituitary-dependent Hypercortisolism. J Vet Intern Med. 2016;30(4):989–95. https://doi.org/10.1111/jvim.14367.
Hanson JM, Mol JA, Leegwater PA, Bilodeau S, Drouin J, Meij BP. Expression and mutation analysis of Tpit in the canine pituitary gland and corticotroph adenomas. Domest Anim Endocrinol 2008;34(3):217–222. doi: https://doi.org/10.1016/j.domaniend.2007.03.002. Epub 2007 May 21. PMID: 17544240.
Fragoso MCBV, Domenice S, Latronico AC, Martin RM, Pereira MAA, Zerbini MCN, et al. Cushing’s syndrome secondary to adrenocorticotropin-independent macronodular adrenocortical hyperplasia due to activating mutations of GNAS1 gene. J Clin Endocrinol Metab. 2003;88(5):2147–51. https://doi.org/10.1210/jc.2002-021362.
Stratakis CA, Tichomirowa MA, Boikos S, Azevedo MF, Lodish M, Martari M, et al. The role of germline AIP, MEN1, PRKAR1A, CDKN1B and CDKN2C mutations in causing pituitary adenomas in a large cohort of children, adolescents, and patients with genetic syndromes. Clin Genet. 2010;78(5):457–63. https://doi.org/10.1111/j.1399-0004.2010.01406.x.
Albani A, Theodoropoulou M, Reincke M. Genetics of Cushing’s disease. Clin Endocrinol. 2018;88(1):3–12. https://doi.org/10.1111/cen.13457.
Hernández-Ramírez LC, Stratakis CA. Genetics of Cushing’s syndrome. Endocrinol Metab Clin N Am. 2018;47(2):275–97. https://doi.org/10.1016/j.ecl.2018.02.007.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. 2020). https://doi.org/10.1038/s41586-020-2308-7;581(7809):434–43.
Boeva V, Popova T, Bleakley K, Chiche P, Cappo J, Schleiermacher G, et al. Control-FREEC: a tool for assessing copy number and allelic content using next-generation sequencing data. Bioinformatics. 2012;28(3):423–5. https://doi.org/10.1093/bioinformatics/btr670.
Olshen AB, Bengtsson H, Neuvial P, Spellman PT, Olshen RA, Seshan VE. Parent-specific copy number in paired tumor-normal studies using circular binary segmentation. Bioinformatics. 2011;27(15):2038–46. https://doi.org/10.1093/bioinformatics/btr329.
Koboldt DC, Zhang Q, Larson DE, Shen D, McLellan MD, Lin L, et al. VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 2012;22(3):568–76. https://doi.org/10.1101/gr.129684.111.
Cibulskis K, Lawrence MS, Carter SL, Sivachenko A, Jaffe D, Sougnez C, et al. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat Biotechnol. 2013;31(3):213–9. https://doi.org/10.1038/nbt.2514.
Zhang J, Feng B. cghMCR: find chromosome regions showing common gains/losses. R package version 1.48.0 2020.
John G Tate, Sally Bamford, Harry C Jubb, Zbyslaw Sondka, David M Beare, Nidhi Bindal, Harry Boutselakis, Charlotte G Cole, Celestino Creatore, Elisabeth Dawson, Peter Fish, Bhavana Harsha, Charlie Hathaway, Steve C Jupe, Chai Yin Kok, Kate Noble, Laura Ponting, Christopher C Ramshaw, Claire E Rye, Helen E Speedy, Ray Stefancsik, Sam L Thompson, Shicai Wang, Sari Ward, Peter J Campbell, Simon A Forbes, COSMIC: the Catalogue Of Somatic Mutations In Cancer, Nucleic Acids Research, Volume 47, Issue D1, 08 January 2019, Pages D941–D947, https://doi.org/10.1093/nar/gky1015
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164. https://doi.org/10.1093/nar/gkq603.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11. https://doi.org/10.1093/bioinformatics/btp120.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012;7(3):562–78. https://doi.org/10.1038/nprot.2012.016.
Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9. https://doi.org/10.1093/bioinformatics/btu638.
Chen J, Jian X, Deng S, Ma Z, Shou X, Shen Y, et al. Identification of recurrent USP48 and BRAF mutations in Cushing’s disease. Nat Commun. 2018;9(1):3171. https://doi.org/10.1038/s41467-018-05275-5.
Bi WL, Greenwald NF, Ramkissoon SH, Abedalthagafi M, Coy SM, Ligon KL, et al. Clinical identification of oncogenic drivers and copy-number alterations in pituitary tumors. Endocrinology. 2017;158(7):2284–91. https://doi.org/10.1210/en.2016-1967.
Qi J, Ni W. Attenuation of MAMLD1 expression suppresses the growth and migratory properties of Gonadotroph pituitary adenomas. Pathol Oncol Res. 2020). https://doi.org/10.1007/s12253-019-00615-2”;26(2):937–46.
Nam Y, Sliz P, Song L, Aster JC, Blacklow SC. Structural basis for cooperativity in recruitment of MAML coactivators to notch transcription complexes. Cell. 2006;124(5):973–83. https://doi.org/10.1016/j.cell.2005.12.037.
Fishbein L, Leshchiner I, Walter V, Danilova L, Robertson AG, Johnson AR, et al. Comprehensive molecular characterization of Pheochromocytoma and Paraganglioma. Cancer Cell. 2017;31(2):181–93. https://doi.org/10.1016/j.ccell.2017.01.001.
Fukami M, Wada Y, Okada M, Kato F, Katsumata N, Baba T, Morohashi KI, Laporte J, Kitagawa M, Ogata T (2008) Mastermind-like domain-containing 1 (MAMLD1 or CXorf6) transactivates the Hes3 promoter, augments testosterone production, and contains the SF1 target sequence. J Biol Chem 283:5525–5532. M703289200, 9. https://doi.org/10.1074/jbc.
Monahan P, Rybak S, Raetzman LT. The notch target gene HES1 regulates cell cycle inhibitor expression in the developing pituitary. Endocrinology. 2009;150(9):4386–94. https://doi.org/10.1210/en.2009-0206.
Goldberg LB, Aujla PK, Raetzman LT. Persistent expression of activated notch inhibits corticotrope and melanotrope differentiation and results in dysfunction of the HPA axis. Dev Biol. 2011;358(1):23–32. https://doi.org/10.1016/j.ydbio.2011.07.004.
Evans C-O, Young AN, Brown MR, Brat DJ, Parks JS, Neish AS, et al. Novel patterns of gene expression in pituitary adenomas identified by complementary deoxyribonucleic acid microarrays and quantitative reverse transcription-polymerase chain reaction. None. 2001;86(7):3097–107. https://doi.org/10.1210/jcem.86.7.7616.
Morris DG, Musat M, Czirják S, et al. Differential gene expression in pituitary adenomas by oligonucleotide array analysis. Eur J Endocrinol. 2005;153(1):143–51. https://doi.org/10.1530/eje.1.01937.
We appreciate the support of Emory University, Winship – Pathology Corelab services, specifically Dianne Alexis, MPH, who supported section histology and staining.
Al Lerner Chair Award, Department of Neurosurgery, Emory University School of Medicine [AW, NMO], supported experimentation expenses.
Ethics approval and consent to participate
This study was submitted and approved by the review and ethics committee of each institution.
This study protocol for the acquisition of human samples was approved by the Emory University Institutional Review Board (IRB00045827). Written and verbal informed consent was obtained from the participating patients preoperatively, and all human participant research were performed in accordance with relevant guidelines and regulations.
The study protocol for the acquisition of all dog samples was approved by the Ethics Committee on Animal Experimentation of Utrecht University, the Netherlands, in accordance to the 3R (replacement, reduction, and refinement of animal experimentation)-policy (DEC2012.III.02.017). Written and verbal informed consent was obtained from all client owners, and canine treatment adhered to the legal and institutional animal welfare guidelines. All live vertebrae experiments were performed in accordance with relevant guidelines and regulations, and released to the client owners.
Consent for publication
All authors had access to the study data and reviewed and approved the final manuscript. All authors report no conflict of interests.
All authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Post-agreement Immunohistochemistry Chemistry Scoring. Individualized immunohistochemistry scoring (0 to 3+) of Cushing’s pituitary adenoma human specimens from 7 antibodies (k = 0.87). Fig. S1. Difference in Magnitude of RNA-seq T-test Results. Normalized averaged FPKM values and standard errors are shown for genes suspected to be overexpressed by RNA-seq in human corticotroph pituitary adenomas (HC PA, n = 6) and dog corticotroph pituitary adenomas (DC PA, n = 7). Fig. S2. GeneCards® Summary Characteristics of Highly Expressed Genes in Humans and Dogs with Cushing’s disease. Fig. S3. Theoretical Schematic Representation of the Highly Expressed Genes in Humans and Dogs with Cushing’s disease.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Wang, A., Neill, S.G., Newman, S. et al. The genomic profiling and MAMLD1 expression in human and canines with Cushing’s disease. BMC Endocr Disord 21, 185 (2021). https://doi.org/10.1186/s12902-021-00845-z
- Cushing’s disease
- Pituitary adenomas
- Dog and man sequencing