Complex signatures of locus-specific selective pressures and gene conversion on Human Growth Hormone/Chorionic Somatomammotropin genes

Growth hormone proteoformics atlas created to promote predictive, preventive, and personalized approach in overall management of pituitary neuroendocrine tumors

Human growth hormone (GH) is the indispensable hormone for the maintenance of normal physiological functions of the human body, including the growth, development, metabolism, and even immunoregulation. The GH is synthesized, secreted, and stored by somatotroph cells in adenohypophysis. Abnormal GH is associated with various GH-related diseases, such as acromegaly, dwarfism, diabetes, and cancer. Currently, some studies found there are dozens or even hundreds of GH proteoforms in tissue and serum as well as a series of GH-binding protein (GHBP) proteoforms and GH receptor (GHR) proteoforms were also identified. The structure-function relationship of protein hormone proteoforms is significantly important to reveal their overall physiological and pathophysiological mechanisms. We propose the use of proteoformics to study the relationship between every GH proteoform and different physiological/pathophysiological states to clarify the pathogenic mechanism of GH-related disease such as pituitary neuroendocrine tumor and conduct precise molecular classification to promote predictive preventive personalized medicine (PPPM / 3P medicine). This article reviews GH proteoformics in GH-related disease such as pituitary neuroendocrine tumor, which has the potential role to provide novel insight into pathogenic mechanism, discover novel therapeutic targets, identify effective GH proteoform biomarker for patient stratification, predictive diagnosis, and prognostic assessment, improve therapy method, and further accelerate the development of 3P medicine.

Network-based analysis of key regulatory genes implicated in Type 2 Diabetes Mellitus and Recurrent Miscarriages in Turner Syndrome

The information on the genotype-phenotype relationship in Turner Syndrome (TS) is inadequate because very few specific candidate genes are linked to its clinical features. We used the microarray data of TS to identify the key regulatory genes implicated with TS through a network approach. The causative factors of two common co-morbidities, Type 2 Diabetes Mellitus (T2DM) and Recurrent Miscarriages (RM), in the Turner population, are expected to be different from that of the general population. Through microarray analysis, we identified nine signature genes of T2DM and three signature genes of RM in TS. The power-law distribution analysis showed that the TS network carries scale-free hierarchical fractal attributes. Through local-community-paradigm (LCP) estimation we find that a strong LCP is also maintained which means that networks are dynamic and heterogeneous. We identified nine key regulators which serve as the backbone of the TS network. Furthermore, we recognized eight interologs functional in seven different organisms from lower to higher levels. Overall, these results offer few key regulators and essential genes that we envisage have potential as therapeutic targets for the TS in the future and the animal models studied here may prove useful in the validation of such targets.

Human growth hormone proteoform pattern changes in pituitary adenomas: Potential biomarkers for 3P medical approaches

Relevance

Human growth hormone (hGH) is synthesized, stored, and secreted by somatotroph cells in the pituitary gland, and promotes human growth and metabolism. Compared to a normal pituitary, a GH-secreting pituitary adenoma can secrete excessive GH to cause pathological changes in body tissues. GH proteoform changes would be associated with GH-related disease pathogenesis.

Purpose

This study aimed to elucidate changes in GH proteoforms between GH-secreting pituitary adenomas and control pituitaries for the predictive diagnostics, targeted prevention, and personalization of medical services.

Methods

The isoelectric point (pI) and relative molecular mass (Mr) are two basic features of a proteoform that can be used to effectively array and detect proteoforms with two-dimensional gel electrophoresis (2DGE) and 2DGE-based western blot. GH proteoforms were characterized with liquid chromatography (LC) and mass spectrometry (MS). Phosphoproteomics, ubiquitinomics, acetylomics, and bioinformatics were used to analyze post-translational modifications (PTMs) of GH proteoforms in GH-secreting pituitary adenoma tissues and control pituitaries.

Results

Sixty-six 2D gel spots were found to contain hGH, including 46 spots (46 GH proteoforms) in GH-secreting pituitary adenomas and 35 spots (35 GH proteoforms) in control pituitaries. Further, 35 GH proteoforms in control pituitary tissues were matched with 35 of 46 GH proteoforms in GH-secreting pituitary adenoma tissues; and 11 GH proteoforms were presented in only GH-secreting pituitary adenoma tissues but not in control pituitary tissues. The matched 35 GH proteoforms showed quantitative changes in GH-secreting pituitary adenomas compared to the controls. The quantitative levels of those 46 GH proteoforms in GH-secreting pituitary adenomas were significantly different from those 35 GH proteoforms in control pituitaries. Meanwhile, different types of PTMs were identified among those GH proteoforms. Phosphoproteomics identified phosphorylation at residues Ser77, Ser132, Ser134, Thr174, and Ser176 in hGH. Ubiquitinomics identified ubiquitination at residue Lys96 in hGH. Acetylomics identified acetylation at reside Lys171 in hGH. Deamination was identified at residue Asn178 in hGH.

Conclusion

These findings provide the first hGH proteoform pattern changes in GH-secreting pituitary adenoma tissues compared to control pituitary tissues, and the status of partial PTMs in hGH proteoforms. Those data provide in-depth insights into biological roles of hGH in GH-related diseases, and identify hGH proteoform pattern biomarkers for treatment of a GH-secreting pituitary adenoma in the context of 3P medicine -predictive diagnostics, targeted prevention, and personalization of medical services.

Supplementary information

The online version contains supplementary material available at 10.1007/s13167-021-00232-7.

Extensive load of somatic CNVs in the human placenta

Placenta is a temporary, but indispensable organ in mammalian pregnancy. From its basic nature, it exhibits highly invasive tumour-like properties facilitating effective implantation through trophoblast cell proliferation and migration, and a critical role in pregnancy success. We hypothesized that similarly to cancer, somatic genomic rearrangements are promoted in the support of placental function. Here we present the first profiling of copy number variations (CNVs) in human placental genomes, showing an extensive load of somatic CNVs, especially duplications and suggesting that this phenomenon may be critical for normal gestation. Placental somatic CNVs were significantly enriched in genes involved in cell adhesion, immunity, embryonic development and cell cycle. Overrepresentation of imprinted genes in somatic duplications suggests that amplified gene copies may represent an alternative mechanism to support parent-of-origin specific gene expression. Placentas from pregnancy complications exhibited significantly altered CNV profile compared to normal gestations, indicative to the clinical implications of the study.

Coherent somatic mutation in autoimmune disease

BACKGROUND

Many aspects of autoimmune disease are not well understood, including the specificities of autoimmune targets, and patterns of co-morbidity and cross-heritability across diseases. Prior work has provided evidence that somatic mutation caused by gene conversion and deletion at segmentally duplicated loci is relevant to several diseases. Simple tandem repeat (STR) sequence is highly mutable, both somatically and in the germ-line, and somatic STR mutations are observed under inflammation.

RESULTS

Protein-coding genes spanning STRs having markers of mutability, including germ-line variability, high total length, repeat count and/or repeat similarity, are evaluated in the context of autoimmunity. For the initiation of autoimmune disease, antigens whose autoantibodies are the first observed in a disease, termed primary autoantigens, are informative. Three primary autoantigens, thyroid peroxidase (TPO), phogrin (PTPRN2) and filaggrin (FLG), include STRs that are among the eleven longest STRs spanned by protein-coding genes. This association of primary autoantigens with long STR sequence is highly significant (p<3.0x10(-7)). Long STRs occur within twenty genes that are associated with sixteen common autoimmune diseases and atherosclerosis. The repeat within the TTC34 gene is an outlier in terms of length and a link with systemic lupus erythematosus is proposed.

CONCLUSIONS

The results support the hypothesis that many autoimmune diseases are triggered by immune responses to proteins whose DNA sequence mutates somatically in a coherent, consistent fashion. Other autoimmune diseases may be caused by coherent somatic mutations in immune cells. The coherent somatic mutation hypothesis has the potential to be a comprehensive explanation for the initiation of many autoimmune diseases.

Population-specific differences in gene conversion patterns between human SUZ12 and SUZ12P are indicative of the dynamic nature of interparalog gene conversion

Nonallelic homologous gene conversion (NAHGC) resulting from interparalog recombination without crossover represents an important influence on the evolution of duplicated sequences in the human genome. In 17q11.2, different paralogous sequences mediate large NF1 deletions by nonallelic homologous recombination with crossover (NAHR). Among these paralogs are SUZ12 and its pseudogene SUZ12P which harbour the breakpoints of type-2 (1.2-Mb) NF1 deletions. Such deletions are caused predominantly by mitotic NAHR since somatic mosaicism with normal cells is evident in most patients. Investigating whether SUZ12 and SUZ12P have also been involved in NAHGC, we observed gene conversion tracts between these paralogs in both Africans (AFR) and Europeans (EUR). Since germline type-2 NF1 deletions resulting from meiotic NAHR are very rare, the vast majority of the gene conversion tracts in SUZ12 and SUZ12P are likely to have resulted from mitotic recombination during premeiotic cell divisions of germ cells. A higher number of gene conversion tracts were noted within SUZ12 and SUZ12P in AFR as compared to EUR. Further, the distinctive signature of NAHGC (a high number of SNPs per paralog and a high number of shared SNPs between paralogs), a characteristic of many actively recombining paralogs, was observed in both SUZ12 and SUZ12P but only in AFR and not in EUR. A novel polymorphic 2.3-kb deletion in SUZ12P was identified which exhibited a high allele frequency in EUR. We postulate that this interparalog structural difference, together with low allelic recombination rates, could have caused a reduction in NAHGC between SUZ12 and SUZ12P during human evolution.

Developmental programming of growth: genetic variant in GH2 gene encoding placental growth hormone contributes to adult height determination

INTRODUCTION

Given the physiological role of placental growth hormone (PGH) during intrauterine development and growth, genetic variation in the coding Growth hormone 2 (GH2) gene may modulate developmental programming of adult stature. Two major GH2 variants were described worldwide, determined by single polymorphism (rs2006123; c.171 + 50C > A). We sought to study whether GH2 variants may contribute to adult anthropometric measurements.

METHODS

Genotyping of GH2 SNP rs2006123 by RFLP, testing its genetic association with adult height and Body Mass Index (BMI) by linear regression analysis, and combining the results of three individual study samples in meta-analysis.

STUDY SAMPLES

HYPEST (Estonia), n = 1464 (506 men/958 women), CADCZ (Czech), n = 871 (518/353); UFA (Bashkortostan), n = 954 (655/299); meta-analysis, n = 3289 (1679/1610).

RESULTS

Meta-analysis across HYPEST, CADCZ and UFA samples (n = 3289) resulted in significant association of GH2 rs2006123 with height (recessive model: AA-homozygote effect: beta (SE) = 1.26 (0.46), P = 5.90 × 10⁻³; additive model: A-allele effect: beta (SE) = 0.45 (0.18), P = 1.40 × 10⁻²). Among men (n = 1679), the association of the A-allele with taller stature remained significant after multiple-testing correction (additive effect: beta = 0.86 (0.28), P = 1.83 × 10⁻³). No association was detected with BMI. Notably, rs2006123 was in strong LD (r² ≥ 0.87) with SNPs significantly associated with height (rs2665838, rs7209435, rs11658329) and mapped near GH2 in three independent meta-analyses of GWA studies.

CONCLUSIONS

This is the first study demonstrating a link between a placental gene variant and programming of growth potential in adulthood. The detected association between PGH encoding GH2 and adult height promotes further research on the role of placental genes in prenatal programming of human metabolism.

Differential placental expression profile of human Growth Hormone/Chorionic Somatomammotropin genes in pregnancies with pre-eclampsia and gestational diabetes mellitus

The human GH/CSH cluster consisting of one pituitary-expressed (GH1) and four placenta-expressed loci has been implicated in maternal metabolic adaptation to pregnancy, regulation of intrauterine and postnatal growth. We investigated how the mRNA expression profile of placental GH2, CSH1 and CSH2 genes and their alternative transcripts correlates with maternal pre-eclampsia (PE) and/or gestational diabetes mellitus (GD). The expression of studied genes in PE placentas (n=17) compared to controls (n=17) exhibited a trend for reduced transcript levels. The alternative transcripts retaining intron 4, GH2-2 and CSH1-2 showed significantly reduced expression in PE cases without growth restriction (P=0.007, P=0.008, respectively). In maternal GD (n=23), a tendency of differential expression was detected only for the GH2 gene and in pregnancies with large-for-gestational-age newborns. Our results, together with those reported by others, are consistent with a pleiotropic effect of placental hGH/CSH genes at the maternal-fetal interface relating to the regulation of fetal growth and the risk of affected maternal metabolism.

Chimeric genes as a source of rapid evolution in Drosophila melanogaster

Chimeric genes form through the combination of portions of existing coding sequences to create a new open reading frame. These new genes can create novel protein structures that are likely to serve as a strong source of novelty upon which selection can act. We have identified 14 chimeric genes that formed through DNA-level mutations in Drosophila melanogaster, and we investigate expression profiles, domain structures, and population genetics for each of these genes to examine their potential to effect adaptive evolution. We find that chimeric gene formation commonly produces mid-domain breaks and unites portions of wholly unrelated peptides, creating novel protein structures that are entirely distinct from other constructs in the genome. These new genes are often involved in selective sweeps. We further find a disparity between chimeric genes that have recently formed and swept to fixation versus chimeric genes that have been preserved over long periods of time, suggesting that preservation and adaptation are distinct processes. Finally, we demonstrate that chimeric gene formation can produce qualitative expression changes that are difficult to mimic through duplicate gene formation, and that extremely young chimeric genes (d(S) < 0.03) are more likely to be associated with selective sweeps than duplicate genes of the same age. Hence, chimeric genes can serve as an exceptional source of genetic novelty that can have a profound influence on adaptive evolution in D. melanogaster.

Gene conversions in the growth hormone gene family of primates: stronger homogenizing effects in the Hominidae lineage

In humans, the growth hormone/chorionic somatomammotropin gene family is composed of five highly similar genes. We characterized the gene conversions that occurred between the growth hormone genes of 11 primate species. We detected 48 conversions using GENECONV and others were only detected using phylogenetic analyses. Gene conversions were detected in all species analyzed, their average size (±standard deviation) is 197.8±230.4 nucleotides, the size of the conversions is correlated with sequence similarity and converted regions are significantly more GC-rich than non-converted regions. Gene conversions have a stronger homogenizing effect in Hominidae genes than in other primate species. They are also less frequent in conserved gene regions and towards functionally important genes. This suggests that the high degree of sequence similarity observed between the growth hormone genes of primate species is a consequence of frequent gene conversions in gene regions which are under little selective constraints.

Differential expression profile of growth hormone/chorionic somatomammotropin genes in placenta of small- and large-for-gestational-age newborns

CONTEXT

The human growth hormone/chorionic somatomammotropin (hGH/CSH) locus at 17q22-24, consisting of one pituitary-expressed postnatal (GH1) and four placenta-expressed genes (GH2, CSH1, CSH2, and CSHL1), is implicated in regulation of postnatal and intrauterine growth. A positive correlation has been reported between the offspring's birth weight and serum placental GH (coded by GH2) and placental lactogen (coded by CSH1, CSH2) levels in pregnant women.

OBJECTIVE

The objective of the study was the investigation of the hypothesis that the mRNA expression profile of placental hGH/CSH genes contributes to the determination of birth weight.

DESIGN AND SUBJECTS

We developed a sensitive, fluorescent-labeled semiquantitative RT-PCR assay coupled with gene-specific restriction analysis, capable of distinguishing alternative splice-products of individual placental hGH/CSH genes and quantification of their relative expression levels. The detailed profile of alternative transcripts of GH2, CSH1, CSH2, and CSHL1 genes in placenta from uncomplicated term pregnancies of the REPROMETA sample collection was addressed in association with the birth weight of newborns, grouped as appropriate for gestational age (AGA; n = 23), small for gestational age (SGA; n = 15), and large for gestational age (LGA; n = 34).

RESULTS

The majority of pregnancies with SGA newborn showed down-regulation of the entire hGH/CSH cluster in placenta, whereas in the case of LGA, the expression of CSH1-1, CSH2-1, and CSHL1-4 mRNA transcripts in placenta was significantly increased compared with AGA newborns (P < 0.0001, P = 0.009, P = 0.002, respectively).

CONCLUSION

The expression profile of placental hGH/CSH genes in placenta is altered in pregnancies accompanied by SGA and LGA compared with AGA newborns, and thus, it may directly affect the circulating fetal and maternal placental GH and placental lactogen levels.

Growth hormone (GH1) gene variation and the growth hormone receptor (GHR) exon 3 deletion polymorphism in a West-African population

Among Europeans, functionally significant GH1 gene variants occur not only in individuals with idiopathic growth hormone (GH) deficiency and/or short stature but also fairly frequently in the general population. To assess the generality of these findings, 163 individuals from Benin, West Africa were screened for mutations and polymorphisms in their GH1 genes. A total of 37 different sequence variants were identified in the GH1 gene region, 24 of which occurred with a frequency of >1%. Although four of these variants were novel missense substitutions (Ala13Val, Arg19His, Phe25Tyr and Ser95Arg), none of these had any measurable effect on either GH function or secretion in vitro. Some 37 different GH1 promoter haplotypes were identified, 23 of which are as yet unreported in Europeans. The mean in vitro expression level of the GH1 promoter haplotypes observed in the African population was significantly higher than that previously measured in Britons (p<0.001). A gene conversion in the GH1 promoter, previously reported in a single individual of British origin, was found to occur at polymorphic frequency (5%) in the West-African population and was associated with a 1.7-fold increase in promoter activity relative to the wild-type. The d3 allele of the GHR exon 3 deletion polymorphism, known to be associated with increased GH responsiveness, was also found to occur at an elevated frequency in these individuals from Benin. We speculate that both elevated GH1 gene expression and increased GHR-mediated GH responsiveness may constitute adaptive responses to the effects of scarce food supply in this West-African population since increased circulating GH appears to form part of a physiological response to nutritional deprivation.

Chorionic gonadotropin beta-gene variants are associated with recurrent miscarriage in two European populations

CONTEXT

The incidence of recurrent miscarriage (RM) (>or=3 consecutive pregnancy losses) is estimated as 1-2% in fertile couples. Familial clustering of RM has suggested the contribution of a genetic component.

OBJECTIVE

A low level of human chorionic gonadotropin (HCG) in maternal serum during the first trimester of the pregnancy is a clinically accepted risk factor for miscarriage. We sought to study whether variation in chorionic gonadotropin beta-subunit genes (CGBs) expressed in placenta may contribute to the risk of RM.

DESIGN

Resequencing of CGB5 and CGB8, the two most actively transcribed loci of the four HCG beta-duplicate genes, was performed.

SETTING

A case-control study involving two sample sets, from Estonia (n = 194) and Finland (n = 185), was performed.

PATIENTS

RM patients (n = 184) and fertile controls (n = 195) participated in the study.

RESULTS

From 71 identified variants in CGB5 and CGB8, 48 polymorphisms were novel. Significant protective effect was associated with two single nucleotide polymorphisms located at identical positions in intron 2 in both CGB5 [P = 0.007; odds ratio (OR) = 0.53] and CGB8 (P = 0.042; OR = 0.15), and with four CGB5 promoter variants (P < 0.03; OR = 0.54-0.58). The carriers of minor alleles had a reduced risk of RM. The haplotype structure of the CGB8 promoter was consistent with balancing selection; a rare mutation in CGB8 initiator element was detected only among patients (n = 3). In addition, three rare nonsynonymous substitutions were identified among RM cases as possible variants increasing the risk of recurrent pregnancy loss.

CONCLUSION

The findings encourage studying the functional effect of the identified variants on CGB expression and HCG hormone activity to elucidate further the role of CGB variation in RM.

High divergence in primate-specific duplicated regions: human and chimpanzee chorionic gonadotropin beta genes

Background

Low nucleotide divergence between human and chimpanzee does not sufficiently explain the species-specific morphological, physiological and behavioral traits. As gene duplication is a major prerequisite for the emergence of new genes and novel biological processes, comparative studies of human and chimpanzee duplicated genes may assist in understanding the mechanisms behind primate evolution. We addressed the divergence between human and chimpanzee duplicated genomic regions by using Luteinizing Hormone Beta (LHB)/Chorionic Gonadotropin Beta (CGB) gene cluster as a model. The placental CGB genes that are essential for implantation have evolved from an ancestral pituitary LHB gene by duplications in the primate lineage.

Results

We shotgun sequenced and compared the human (45,165 bp) and chimpanzee (39,876 bp) LHB/CGB regions and hereby present evidence for structural variation resulting in discordant number of CGB genes (6 in human, 5 in chimpanzee). The scenario of species-specific parallel duplications was supported (i) as the most parsimonious solution requiring the least rearrangement events to explain the interspecies structural differences; (ii) by the phylogenetic trees constructed with fragments of intergenic regions; (iii) by the sequence similarity calculations. Across the orthologous regions of LHB/CGB cluster, substitutions and indels contributed approximately equally to the interspecies divergence and the distribution of nucleotide identity was correlated with the regional repeat content. Intraspecies gene conversion may have shaped the LHB/CGB gene cluster. The substitution divergence (1.8–2.59%) exceeded two-three fold the estimates for single-copy loci and the fraction of transversional mutations was increased compared to the unique sequences (43% versus ~30%). Despite the high sequence identity among LHB/CGB genes, there are signs of functional differentiation among the gene copies. Estimates for d_n/d_s rate ratio suggested a purifying selection on LHB and CGB8, and a positive evolution of CGB1.

Conclusion

If generalized, our data suggests that in addition to species-specific deletions and duplications, parallel duplication events may have contributed to genetic differences separating humans from their closest relatives. Compared to unique genomic segments, duplicated regions are characterized by high divergence promoted by intraspecies gene conversion and species-specific chromosomal rearrangements, including the alterations in gene copy number.