Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 27 November 2020
Sec. Human and Medical Genomics
This article is part of the Research Topic Recent Advances in the Prediction of Functional Effects of Non-Coding Variants and their Clinical Applications View all 12 articles

A Variant in the Nicotinic Acetylcholine Receptor Alpha 3 Subunit Gene Is Associated With Hypertension Risks in Hypogonadic Patients

\r\nTao Wu&#x;Tao Wu1†Yujia Wang,&#x;Yujia Wang2,3†Wei Shi&#x;Wei Shi2†Bi-Qi ZhangBi-Qi Zhang1John RaelsonJohn Raelson2Yu-Mei YaoYu-Mei Yao4Huan-Dong WuHuan-Dong Wu1Zao-Xian XuZao-Xian Xu1Francois-Christophe Marois-BlanchetFrancois-Christophe Marois-Blanchet2Jonathan LedouxJonathan Ledoux5Rikard BlunckRikard Blunck6Jian-Zhong Sheng*Jian-Zhong Sheng7*Shen-Jiang Hu*Shen-Jiang Hu1*Hongyu Luo*Hongyu Luo2*Jiangping Wu,*Jiangping Wu2,8*
  • 1Institute of Cardiology, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, China
  • 2Research Centre, Centre Hospitalier de l’Université de Montréal (CHUM), Montreal, QC, Canada
  • 3Children’s Hospital, Zhejiang University School of Medicine, Hangzhou, China
  • 4Department of Cardiology, The Third Affiliated Hospital of Zhejiang Chinese Medical University, Hangzhou, China
  • 5Montreal Heart Institute, Université de Montréal, Montreal, QC, Canada
  • 6Department of Physics, University of Montreal, Montreal, QC, Canada
  • 7Department of Pathology and Physiopathology, College of Medicine, Zhejiang University, Hangzhou, China
  • 8Nephrology Service, Centre Hospitalier de l’Université de Montréal (CHUM), Montreal, QC, Canada

Ephb6 gene knockout causes hypertension in castrated mice. EPHB6 controls catecholamine secretion by adrenal gland chromaffin cells (AGCCs) in a testosterone-dependent way. Nicotinic acetylcholine receptor (nAChR) is a ligand-gated Ca2+/Na+ channel, and its opening is the first signaling event leading to catecholamine secretion by AGCCs. There is a possibility that nAChR might be involved in EPHB6 signaling, and thus sequence variants of its subunit genes are associated with hypertension risks. CHRNA3 is the major subunit of nAChR used in human and mouse AGCCs. We conducted a human genetic study to assess the association of CHRNA3 variants with hypertension risks in hypogonadic males. The study cohort included 1,500 hypogonadic Chinese males with (750 patients) or without (750 patients) hypertension. The result revealed that SNV rs3743076 in the fourth intron of CHRNA3 was significantly associated with hypertension risks in the hypogonadic males. We further showed that EPHB6 physically interacted with CHRNA3 in AGCCs, providing a molecular basis for nAChR being in the EPHB6 signaling pathway.

Introduction

Erythropoietin-producing hepatocellular kinases (EPHs) comprise the largest family of receptor tyrosine kinases (Eph Nomenclature Committee, 1997). Their ligands are called ephrins (EFNs), which are also cell surface molecules. When EPH and EFN on neighboring cells interact, both molecules are capable of transducing signals into their respective cells. Interactions between EPHs and EFNs are promiscuous. One EPH can interact with multiple EFNs and vice versa. In general, EPHA members bind preferentially to EFNA members, and EPHB members, to EFNB members (Pasquale, 2008).

Erythropoietin-producing hepatocellular kinases and EFNs have functions in many biological systems and processes, such as the nervous (Flanagan and Vanderhaeghen, 1998; Wilkinson, 2000), the vascular (Kuijper et al., 2007), the digestive (Batlle et al., 2002), the immune (Wu and Luo, 2005), and the endocrine systems (Konstantinova et al., 2007). Our lab is the first to report the function of EPHs and EFNs in the immune system (Luo et al., 2001, 2002). We have recently discovered the previously unknown functions of EPHs and EFNs in regulating blood pressure (Luo et al., 2012; Wu et al., 2012; Wang et al., 2015). Our human genetics studies have provided corroborating evidence that five single-nucleotide variants (SNVs) in the EFNB2 (Wang et al., 2016) and two SNVs in the EFNB3 (Tremblay et al., 2017) are significantly associated with human hypertension in a sex-specific manner.

The deletion of EPHB6 results in increased vascular smooth muscle cell (VSMC) contractility but does not affect blood pressure in male Ephb6 knockout (KO) mice unless the mice are castrated: Ephb6 KO plus castration leads to increased blood pressure (Luo et al., 2012). We have shown that, in addition to its effects on VSMCs, EPHB6 also regulates catecholamine biosynthesis (Shi et al., 2019) and secretion by adrenal gland chromaffin cells (AGCCs) (Luo et al., 2012). Male Ephb6 KO mice have reduced 24-hour urine catecholamine levels, counteracting the outcome of increased VSMC contractility, and the sum of these two effects results in normal blood pressure. The castration of male Ephb6 KO mice leads to a return to normal levels of catecholamine secretion (Luo et al., 2012). This change, concomitantly with enhanced VSMC contractility, results in blood pressure elevation in the castrated Ephb6 KO mice. These results suggest that EPHB6 and male sex hormones act in concert to regulate catecholamine secretion and blood pressure.

Our further investigation revealed that AGCCs from male Ephb6 KO mice presented reduced acetylcholine (ACh)-triggered Ca2+ influx, depending on the non-genomic effects of testosterone (Wang et al., 2018). We demonstrated that the decreased Ca2+ influx resulted from enhanced large-conductance calcium-activated potassium (BK) currents in these cells (Wang et al., 2018).

In addition to the abnormality of BK channels in the Ephb6-KO AGCCs, EPHB6 might regulate other ion channels such as the ACh receptor (AChR), whose opening is the first event in triggering catecholamine release in AGCC. There are two types of AChRs, nicotinic AChR (nAChR) and muscarinic AChR, based on their agonists (nicotine and muscarine). nAChR is a ligand-gated ion channel mediating fast responses of ACh (Itier and Bertrand, 2001), while muscarinic AChR (mAChR) is a G-protein-coupled metabotropic receptor with a slower signaling pace (Eglen, 2006).

nAChR is a pentamer with five homomeric or heteromeric subunits, depending on species and cell types. The subunits consist of α, β, δ. ε and γ types. The α type has ten members (α1–10), β type five members (β1–5). In the mouse AGCCs, α1, α3, α4, α7, α9, β1, β2, and β4 are detectable, with α3 and β4 being most abundant (Mousavi and Nordberg, 2006; Wu et al., 2010), similar to that in rat and bovine AGCCs (Campos-Caro et al., 1997; Di Angelantonio et al.,2003). According to subunit-specific blocker/agonists and the measurement of mRNA, α3b4 nAChR is also present and functional in human AGCCs (Mousavi et al., 2001; Perez-Alvarez et al., 2012).

Multiple proteins can interact with nAChR and positively or negatively modulate its function (Jones et al., 2010). EPHB6 might be one such regulator, and nAChR might lie in the EPHB6 signaling pathway explaining its role in blood pressure control. If so, variants in nAChR subunits should be associated with hypertension risks in a testosterone-dependent way, similar to those of EPHB6. In this study, we conducted a human genetic study to assess the association of variants in CHRNA3, the major subunit of nAChR in AGCCs, with hypertension risks in hypogonadic, hypertensive patients. We also investigated the possible physical interaction of EPHB6 with nAChR.

Materials and Methods

Patient Population

The details of the patient population were previously described (Wu et al., 2018) but are presented here again for the convenience of readers. A total of 4,480 male patients ≥40 years old from the Cardiology Ward, Endocrinology Ward, and Physical Examination Center of First Affiliated Hospital, College of Medicine, Zhejiang University in Hangzhou, China, were recruited for this study. They were tested for total plasma testosterone levels. Those with hypogonadism [plasma total testosterone levels <346 ng/dL, the International Society’s cut-off level recommended for the Study of the Aging Male (Traish et al., 2011)] were retained. Among those hypogonadic individuals, 982 were diagnosed with primary hypertension and were considered as cases. The primary hypertension phenotype was defined as having recorded systolic pressure >140 mm Hg or diastolic pressure >90 mmHg or having been actively treated for hypertension, excluding known conditions or medication use that could cause BP to increase. All other medical conditions were allowed as long as they were not likely to cause hypertension. All types of medication in the last 3 months were allowed except those that are known to affect testosterone levels, such as testosterone replacement therapy. Seven hundred and eighty-eight normotensive hypogonadic patients among the cohort were retained as controls. Seven hundred and fifty cases and 750 controls were selected with an attempt to match their ages as closely as possible. The systolic and diastolic BP and the values of potential covariate parameters (i.e., age, plasma testosterone level, heart rates, body mass index (BMI), serum uric acid levels, and smoking status), which could be implicated hypertension risks, of the case and control groups are shown in Supplementary Tables 1, 2.

Blood Sample Collection and Plasma Total Testosterone Measurements

Venous blood samples were drawn from all subjects after an overnight fast of at least 8 h. Five mL of blood was collected into vacuum tubes with the anticoagulant EDTA-K+ and centrifuged at the collection site within 1 h. Cell pellets were frozen until DNA extraction.

Plasma total testosterone levels were measured with Siemens Immulite 2000 Total Testosterone Kits on Siemens Immulite 2000 Immunoassay Analyzer according to the manufacturer’s protocols.

DNA Extraction and Purification

Sample DNA was extracted using DNeasy Blood and Tissue Kit (Cat. 69506, QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. Purified DNA quantity and quality were assessed by Qubit®2.0 Fluorometer (Q32866, Invitrogen, Carlsbad, CA, United States) and 1% agarose gel electrophoresis. Samples with DNA quantity ≥ 2 μg and optical density (OD) 260 nm/280 nm = 1.8–2.0 were submitted to SNV assay.

Candidate Tag Single Nucleotide Variants (SNVs)

Tag SNVs were chosen from the chromosomal regions at 15q25.1 between positions 78,837,801 and 78,959,037 (Build 37/hg 19), containing the CHRNA3 gene and 50 kb on either side of the gene, using the Tagger Program (De Bakker et al., 2005) with Han Chinese in Beijing (CHB) linkage-disequilibrium (LD) data. Tag SNVs were chosen with minimum r2 > 0.80 and minor allele frequency >0.05. Candidate tag SNVs were then submitted for analysis by Illumina Software for compatibility with the GoldenGate multiplexing process. Alternative Tag SNVs were chosen for those Tag SNVs determined to be incompatible with the multiplexing technology, and the new Tag SNVs were then re-submitted for a new Illumina software analysis to determine the compatibility of the new set of SNVs including previously compatible SNVs and the new alternative Tag SNVs. This process continued iteratively until all Tag SNVs were found to be compatible for use by the multiplexing technology. Ultimately, 14 tag SNVs were chosen for this region. The Bonferroni-corrected critical p-values (Pcrit = 0.05/14 = 0.0037) were calculated for SNVs with the hypertension association test, assuming tag SNVs represent independent statistical tests, when performing association analysis across the region.

Tag SNV Genotyping

The tag SNVs were genotyped by the Shanghai Biotechnology Corporation using the Illumina GoldenGate genotyping platform according to the manufacturer’s instructions. Those SNVs with a call rate of less than 90% were filtered out and not analyzed.

Association Analysis

The genotyped SNVs were tested for Hardy-Weinberg equilibrium and were analyzed for association with hypertensive versus normotensive status using the PLINK program (Purcell et al., 2007) and a logistic regression model both with and without covariates, the choice of which is described in the Results section. Each SNV was assigned a reference allele and an alternative allele, the SNV of the reference allele being given a numeric code of 1, and the alternative allele, 0. The sum of the numeric codes of the whole reference allele in each individual was used as the genotype value. The sums of the genotype values of all cases and controls were calculated, respectively, and were entered into the logistic regression equation as the genotype terms for each group. This approach is analogous to an additive genetic model.

Reverse Transcription-Quantitative Polymerase Chain Reaction (RT-qPCR)

mRNA levels of Chrna3 and Chrnb4 in the mouse adrenal glands were measured by RT-qPCR. Total RNA from adrenal glands of male Ephb6 KO and WT mice was extracted with TRIzol® (Invitrogen, Burlington, ON, Canada) and reverse-transcribed with iScriptTM cDNA Synthesis Kit (Bio-Rad Laboratories (Canada) Ltd., Mississauga, ON, Canada). The sequences of qPCR primers used were listed in Table 1. The qPCR condition was as follows: 2 min at 50°C, 2 min at 95°C, followed by 40 cycles of 10 s at 94°C, 20 s at 58°C, and 20 s at 72°C. β-actin mRNA levels were used as internal controls. Data from 20 to 30 cycles of amplification were used. Samples were assayed in duplicate. Data were expressed as signal ratios of target RNA/β-actin mRNA.

TABLE 1
www.frontiersin.org

Table 1. RT-qPCR primer sequences.

Primary AGCC Culture

The adrenal glands from 8- to 10-week old mice were isolated, and the fat and the cortex were removed from these glands. Papain (P4762, Sigma-Aldrich, Oakville, ON, Canada) was activated with 5 mM L-cysteine. Adrenal gland medullae were digested with papain in Hank’s buffer (2 medullae/100 ml Hank’s buffer containing 4 units of activated papain) at 37°C for 25 min. They were washed twice with Hank’s buffer and then triturated by pipetting in 300 ml Hank’s buffer until they became feather-like. Cells were pelleted at 3,700 g for 3 min and resuspended in DMEM containing 15% FCS for culture.

Immunoprecipitation and Immunoblotting

HEK293 cells were co-transfected with plasmid pCEP4-HA-mEphb6 (expressing the HA-tagged mouse EPHB6 intracellular domain (aa 2,268 to 3,537)) and plasmid MR226689 (OriGene, Rockville, MD, United States; expressing Myc-tagged mouse CHRNA3), or control empty pCMV6-Entry Vector (PS100001; OriGene). After 48 h, the cells were lysed using radioimmunoprecipitation assay buffer (RIPA), which contained a cocktail of protease inhibitors (Roche Applied Science, Meylan, France). The lysates were pre-cleared with 30 μl protein G-agarose beads (17061801; GE Healthcare Life Sciences, Mississauga, ON, Canada) and then precipitated with mouse anti-Myc mAb (sc-40; Santa Cruz Biotechnology, Dallas, United States) or normal mouse IgG (sc-2025; Santa Cruz Biotechnology) plus protein G-agarose beads at 4°C with gentle rotation overnight. The protein-bead complexes were mixed with SDS-PAGE-loading buffer and boiled for 5 min. Samples were resolved in 10% SDS-PAGE and transferred to PVDF membranes, which were blotted with rabbit anti-HA mAb (3724; Cell Signaling Technology, Danvers, MA, United States), mouse anti-Myc Ab (sc-40; Santa Cruz Biotechnology), and rabbit anti-b-actin Ab (4967; Cell Signaling Technology). Blots were washed and incubated with fluorophore-conjugated secondary Abs (IRDye®680RD-conjugated donkey anti-mouse IgG Ab, P/N 925-68072, LI-COR Biosciences, Lincoln, NE, United States; or IRDye®800CW-conjugated donkey anti-rabbit IgG Ab, P/N 925-32213, LI-COR Biosciences). All the Abs were used at the manufacturer’s recommended dilutions. Signals were visualized by the Odyssey infrared imaging system (LI-COR Biosciences).

Immunofluorescence Microscopy

Adrenal gland chromaffin cells were cultured in 6-well plates with cover glass placed at the bottom of the wells. After one day, the cells were washed once with PBS and fixed with paraformaldehyde (4%) for 20 min. The cells were blocked with 10% FBS in PBS for 20 min and then incubated with goat anti-mouse EPHB6 Ab (AF611; R&D Systems, Minneapolis, MN, United States) and rabbit anti-mouse CHRNA3 Ab (ABN281; Millipore, Etobicoke, ON, Canada) overnight at 4°C. Cells were then reacted with Alexa Fluor-488-conjugated donkey anti-goat IgG Ab (A-11055; Thermo Fisher Scientific, Waltham, MA, United States) for EPHB6 and with rhodamine-conjugated goat anti-rabbit IgG Ab (31670; Thermo Fisher Scientific) for CHRNA3 at room temperature for 2 h. The cells were embedded with ProLong® Gold anti-fade reagent (Invitrogen). Cell staining was examined with a Zeiss microscope.

Fluorescence Resonance Energy Transfer Microscopy (FRET)

Adrenal gland chromaffin cells were cultured in 6-well plates with cover glass placed at the bottom of the wells. After one day, the cells were washed once with PBS and fixed with paraformaldehyde (4%) for 20 min. The cells were blocked with 10% FBS in PBS for 20 min and then incubated with goat anti-mouse EPHB6 Ab (AF611; R&D Systems, Minneapolis, MN, United States) and rabbit anti-mouse CHRNA3 Ab (ABN281; Millipore, Etobicoke, ON, Canada) overnight at 4°C. The cells were then reacted with Alexa Fluor-488-conjugated donkey anti-goat IgG Ab (A-11055; Thermo Fisher Scientific, Waltham, MA, United States) for EPHB6 and with rhodamine-conjugated goat anti-rabbit IgG Ab (31670; Thermo Fisher Scientific) for CHRNA3 at room temperature for 2 h. The cells were embedded with ProLong® Gold anti-fade reagent (Invitrogen). The FRET signal was examined under a Leica TCS SP5 laser-scanning confocal microscope (Leica Microsystems Inc., Concord, ON, Canada). Rhodamine was the donor fluorophore, and Alexa Fluor-488, the acceptor fluorophore.

Fluorescence resonance energy transfer microscopy was measured with acceptor photobleaching (AB) by FRET AB Wizard software (Leica Microsystems Inc.). All necessary controls for AB, such as cells with relevant single fluorescence staining, were performed to satisfy background deductions in calculating FRET efficiency, as required by the software.

Acceptor photobleaching FRET efficiency was calculated by fluorescence intensity of the donor before (Dpre) and after (Dpost) acceptor-selective photobleaching, according to the following formula:

AB FRET efficiency = (Dpost − Dpre) / Dpost

Results

We used 750 hypogonadic men with hypertension as cases and 750 age-matched normotensive hypogonadic men as controls to assess the association of CHRNA3 variants with hypertension. The mean measured blood pressure (systolic as well as diastolic) was significantly higher in the cases (125.42 + 0.54/76.45 + 0.36 mmHg, mean + SE) than the controls (116.6 + 0.43/72.06 + 0.32 mmHg) (p = 2.2 × 10–16) for both systolic and diastolic pressures (Supplementary Table 1), but the case blood pressure did not reach the hypertension diagnostic criteria, i.e., systolic pressure >140 mmHg and/or diastolic pressure >90 mmHg. This was due to the fact that the majority of the cases (93.7%) were previously diagnosed and were under anti-hypertension medication, which controlled their measured blood pressure within the normal range.

We performed a covariate analysis to test the effects of non-genetic factors in the determination of case and control status (hypertension) in our sample (Supplementary Table 2). We conducted simple logistic regression analysis on several phenotype candidate covariates that we suspected might differ between cases and controls (variable 1 for cases, 0 for controls). These were age, plasma testosterone levels, heart rates, body mass index (BMI), serum uric acid levels, and smoking status. Each regression was done for a single phenotype covariate value only without any genotype factor included.

Age, heart rate, BMI, and serum uric acid levels were significantly different between cases and controls in the simple regression analyses and were retained for further analysis. In addition, the testosterone level was nearly significant at p = 0.05, so it is also included since the assumption of hypogonadism in the whole cohort is central to this analysis.

Realizing that some of these covariates could be correlated among themselves, we attempted to consider their combined effect on case and control status. In order to do so, we performed a multiple logistic regression with all of the retained candidate covariates but without any genotype factor. All of the retained covariates from the simple regression analyses, when combined in a multiple regression model without a genotype factor, remained significantly different between cases and controls, again with the exception of testosterone levels, which were nearly significant (p = 0.0506). (Supplementary Table 3). Therefore, we retained age, heart rate, BMI, serum uric acid levels, and blood testosterone levels in the genetic regression models subsequently performed for each SNV separately. All covariate factors were included in each of the genetic multiple regression analysis. There may have been slight differences in covariate measures for each SNV due only to differences in missing values, but this effect will have been slight as the data genotype completeness was controlled and variation kept within strict limits. Essentially, the covariate effects remained constant for all SNV analyses, and the significance of each genotype term was tested once the consistent covariate effects were removed.

The genetic association was tested for all SNVs using logistic regression models, both with and without covariates. SNV rs3743706 (chromosome 15, position 78,909,227; build 37/hg 19), which is intronic to CHRNA3 lying 138 bp 5′ of exon 5, was significant at p-values below the critical multiple testing p-value (p-crit) of 0.0037 for models both with and without covariates (p = 0.00190 and p = 0.00143, respectively). Table 2 presents the genotyped tag SNVs along with their positions, reference alleles, minor alleles, minor allele frequencies, odds ratios with respect to the reference alleles, standard errors of those odds ratios, and the p-values for the significance of the odds ratios for association in the models with or without covariates. The SNV with a p-value below p-crit along with its associated parameters is in bold. Figure 1 shows LocusZoom plots illustrating association p-values for all the 14 SNVs calculated without or with covariates, with rs3743706 showing the highest −log10 p-value. The plots also illustrate the location of this significant SNV in the CHRNA3 and the recombination rate in the analyzed region. Analysis with covariates did not reduce the p-value of rs3743706, compared to the analysis without covariates, suggesting that this SNV’s association with hypertension is independent of the phenotypes used in the logistic multiple regression model.

TABLE 2
www.frontiersin.org

Table 2. Logistic regression association tests for SNVs in the CHRNA3 gene region with and without covariates.

FIGURE 1
www.frontiersin.org

Figure 1. LocusZoom plots for an association of all SNVs assayed across the CHRNA3 gene region for regression models with and without covariates. Left Y-axis: –Log10 (p-value) of the SNV association with hypertension. Right Y-axis: recombination rate (cM/Mb). X-axis: the position of SNVs and the CHRNA3 gene on chromosome 15.

The results of the human genetic study were compatible with our hypothesis that AChR in AGCCs lies within the EPHB6 signaling pathway. If so, EPHB6 might also regulate AChR expression or interact with AChR directly to regulate its function. We first investigated whether Ephb6 KO resulted in altered expression of the major nAChR subunits α3 and β4 in AGCCs. As shown in Figure 2, the mRNA expression of these two major subunits in KO AGCCs from male mice remained unchanged, compared to their wild type (WT) counterparts. This excludes the possibility that EPHB6 regulates AChR function by modulating the expression of its major components α3 and β4 at the gene level and favors the possibility that the putative regulation might occur through EPHB6 and AChR interaction at the protein level.

FIGURE 2
www.frontiersin.org

Figure 2. Similar Chrna3 and Chrnb4 mRNA expression levels in AGCCs from WT and EPHB6 KO mice. Mouse Chrna3 and Chrnb4 mRNA levels were measured by RT-qPCR, using total RNA from adrenal glands of male Ephb6 KO and WT mice. β-actin mRNA levels were employed as internal controls. Each sample was analyzed in duplicate. Three independent experiments were conducted, and the pooled data of these experiments were expressed as signal ratios (means ± SE) of target R.N.A./β-actin mRNA. There was no significant difference between WT and Ephb6 KO AGCCs in their Chrna3 and Chrnb4 mRNA levels (two-way paired Student’s t-test).

We then examined the possible interaction between EPHB6 and AChR. HEK293 cells were transfected with the HA-tagged mouse EPHB6 intra-cellular domain and Myc-tagged mouse CHRNA3. As shown in Figure 3, EPHB6-HA and CHRNA3-Myc both presented in the lysates of HEK293 cells that were transfected with EPHB6-HA or CHRNA3-Myc, but not in cells transfected with the control plasmid. EPHB6-HA was found in the anti-Myc (for CHRNA3-Myc) precipitates (first lane) of cells transfected with both EPHB6-HA and CHRNA3-Myc expressing plasmids, but not in the control precipitates from cells transfected with Myc-vector (second lane). CHRNA3 (the two major bands indicated by arrows) only appeared in lysates and precipitates from cells transfected with CHRNA3-Myc plasmid (lanes 1 and 5). There were some minor non-specific bands observed in anti-My Ab blotting. The abundance of the non-specific proteins was different in the lysates and immunoprecipitates, so the patterns of the non-specific protein bands were different in lanes 1 and 5. Nevertheless, the two CHRNA3 bands were present in both the lysates and immunoprecipitates. Non-specific Ab did not precipitate CHRNA3-Myc or EPHB6-HA (middle two lanes). Beta-actin immunoblotting was used to show similar input of lysates. For immunoprecipitates, the similar intensity of b-actin bands were present in all the lanes, including those controls of normal IgG and empty vectors. The b-actin in the immunoprecipitates was due to non-specific carryover by the beads, but such non-specific carryover indicated that the absence of EPHB6 signals in the control lanes was not due to lack of or unequal protein inputs. This result indicates that EPHB6 physically interacts with AChR via its major subunit CHRNA3. Since the EPHB6 intracellular domain-expressing plasmid was used in this experiment, the results also suggest that the EPHB6 intracellular domain is responsible and sufficient to bind to CHRNA3.

FIGURE 3
www.frontiersin.org

Figure 3. Mouse EPHB6 physically interacts with CHRNA3 according to immunoprecipitation. HEK293 cells were co-transfected with plasmids expressing HA-tagged EPHB6 intracellular domain (aa 2268 to 3537) (EPHB6-HA) and Myc-tagged CHRNA3 (CHRNA-MYC), or an empty control vector (MYC-vector), as indicated. The cell lysates were precipitated with anti-Myc mAb, or normal mouse IgG, as shown. The precipitates were resolved in 10% SDS-PAGE and immunoblotted with anti-Myc (first row) and anti-HA (second row) Abs followed by fluorophore-conjugated secondary Abs. The membranes were also blotted with anti-b-actin Ab (row 3) to show equal loading. Signals were visualized by the Odyssey infrared imaging system (LI-COR Biosciences). The experiments were conducted three times, and representative results are presented. Arrows point to CHRNA3 and EPHB6 bands.

To assess the EPHB6 and CHRNA3 interaction in a more physiological setting, we examined EPHB6 and CHRNA3 expression in mouse AGCCs using confocal microscopy and found that these two molecules were colocalized (Figure 4).

FIGURE 4
www.frontiersin.org

Figure 4. EPHB6 co-localized with CHRNA3 in mouse AGCCs according to confocal microscopy. AGCCs from male WT mice were fixed and stained with anti-EPHB6 and anti-CHRNA3 Abs followed by Alexa Fluor-488-conjugated donkey anti-goat IgG Ab and rhodamine-conjugated goat anti-rabbit IgG Ab, respectively. Three independent experiments were conducted, and micrographs from a representative experiment are shown. The rectangles in the left panels were enlarged at the right. EPHB6 was in pseudo green, and CHRNA3, pseudo red.

The possible physical interaction between EPHB6 and CHRNA3 in WT mouse AGCCs was further examined using a different approach: acceptor photobleaching (AB) FRET. During AB, fluorescence from the acceptor EPHB6 was selectively depleted without influencing donor CHRNA3 signals. The donor and acceptor fluorescence intensity before and after photobleaching in the bleached region (rectangles in Figure 5A) and the unbleached region of AGCCs were quantified. Data from a representative cell are listed in the inset table (Figure 5B). After AB, the signal intensity of donor CHRNA3 was significantly increased in the bleached but not in the unbleached area, compared to that observed before photobleaching. This indicates that donor fluorophore energy could not be transferred to the acceptor due to the acceptor fluorophore bleaching, demonstrating that the donor and acceptor existed in proximity. The acceptor fluorophore intensity was reduced in the photobleached but not unbleached area, again as expected. FRET efficiency from CHRNA3 to EPHB6 of areas inside and outside the bleached regions for more than 7–10 randomly selected cells from 3 independent experiments is shown in Figure 5C. The results show that CHRNA3 and EPHB6 in mouse AGCCs are in very close proximity of <10 nm, which is conventionally considered as a direct association.

FIGURE 5
www.frontiersin.org

Figure 5. Mouse EPHB6 interacts with CHRNA3 in WT AGCCs according to FRET. FRET was measured with acceptor photobleaching (AB). All the necessary controls for AB, such as cells with relevant single fluorescence staining, were performed to satisfy background deductions in calculating FRET efficiency. AB FRET efficiency was calculated by fluorescence intensity of the donor before (Dpre) and after (Dpost) acceptor-selective photobleaching, according to the following formula: AB FRET efficiency = (Dpost – Dpre) / Dpost. (A) Merged image of an EPHB6-(pseudo-green) and CHRNA3-(pseudo-red) stained AGCC from a WT male mouse. The rectangle indicates the bleached area. (B) Donor and acceptor signal measurements and AB FRET efficiency of a representative AGCC from a WT male mouse. (C) AB FRET efficiency of bleached and unbleached areas. Data of 7–10 AGCCs from 3 independent experiments were pooled and presented with means ± SE. The differences between the AB FRET efficiency values of bleached areas and unbleached areas were highly significant (p < 0.01; Two-way Student’s t-test).

Discussion

We previously demonstrated that EPHB6 and testosterone jointly regulated adrenal gland catecholamine secretion and blood pressure in mice (Luo et al., 2012). The major event triggering catecholamine secretion in AGCCs from the male EPHB6 KO, Ca2+ influx, was compromised, and this was accompanied by decreased catecholamine release (Wang et al., 2018). We found that one of the causes for the compromised Ca2+ influx was the earlier closure of the BK channel in AGCCs. This premature closure led to quicker AGCC repolarization, hence, to earlier termination of the Ca2+ influx and reduced total Ca2+ influx (Wang et al., 2018). We wondered whether EPHB6 also affected the Ca2+ influx by interacting with other ion channels. nAChR is a ligand-gated ion channel. The opening of nAChR after ligand Ach binding is the first event in AGCC depolarization, which results in the subsequent large Ca2+ influx followed by catecholamine secretion. Thus, nAChR in AGCCs is possibly involved in EPHB6 signaling pathways related to blood pressure regulation. This prompted us to conduct a human genetic study to assess the association of variants of a major nAChR subunit CHRNA3 with hypertension risks in hypogonadic male hypertensive patients. These patients had augmented chances of having detrimental variants in the EPHB6 signaling pathways (including AChR subunits) due to their hypogonadic condition because hypertension was only revealed in castrated Ephb6 KO mice (Luo et al., 2012). Indeed, SNV rs3743706 within CHRNA3 was found to be significantly associated with hypertension risks in this cohort.

The significant SNV, rs3743706, is located in the fourth intron of CHRNA3. This SNV might and might not be a functional variant. This SNV itself or nearby functionally active variants in linkage disequilibrium (LD) with it in this intron could be part of an intronic enhancer/repressor regulating CHRNA3 expression levels. Such intronic regulators are not uncommon. For example, HCN4 has a potent enhancer in its first intron controlling its expression (Kuratomi et al., 2009). The genomic sequence from 100-bp upstream to 100-bp downstream of this SNV was surveyed for the presence of possible regulatory elements using software available in the UCSC Genome Browser (EPDnew, EPDnewNC, ENCODE Candidate Cis-Regulatory Elements (cCREs), GeneHancer, CpG islands, H3K27Ac Mark, ORegAnno, NCBI RefSeq Functional Elements, Double Elite, Hi-C and Micro-C). No conserved enhancer/repressor motifs or transcription factor-binding sites were found in this region (Supplementary Figure 1). We cannot exclude the existence of cryptic and non-canonical regulator elements for which rs3743706 is a part.

Alternatively, the functional variants could be in LD with this SNV and located in the nearby exons or introns. If so, sequence variants in CHRNA3 protein or enhancer/repressor of its gene might be responsible for the alteration of its function or expression, respectively.

CHRNA3 SNVs have not been found to be significantly associated with hypertension risks in many genome-wide association studies (Wellcome Trust Case Control Consortium, 2007; Adeyemo et al., 2009; Levy et al., 2009; Org et al., 2009; Hiura et al., 2010; Padmanabhan et al., 2010; International Consortium for Blood Pressure Genome-Wide Association Studies et al., 2011; Slavin et al., 2011; Guo et al., 2012). However, testosterone levels are not stratified in these studies. Since the majority of the males are testosterone sufficient in the general population, CHRNA3 variants would not be associated with blood pressure in the majority of individuals in these studies, while the smaller subpopulation of patients with both hypogonadism and CHRNA3 variants would be obscured statistically and so could not reach significance in a GWAS even if they were hypertensive, unless samples were stratified according to testosterone levels.

Our genetic study was performed using an ethnic Han cohort in China. Additional studies using a different Han cohort and a cohort of different genetic backgrounds will be needed to confirm and generalize our conclusion.

There are different views regarding whether catecholamines released from the adrenal glands play a critical role in primary hypertension. The adrenal glands produce all three major species of catecholamines: epinephrine, norepinephrine, and dopamine. However, while plasma epinephrine is mainly derived from the adrenal gland, plasma norepinephrine predominantly comes from sympathetic nerve endings. Dopamine levels in plasma are very low. Some early studies showed that epinephrine levels are associated with hypertension in humans and animal models (Axelrod, 1976; Buhler et al., 1982; Dominiak and Grobecker, 1982; Goldstein, 1983; Borkowski and Quinn, 1984; Jablonskis and Howe, 1994). In recent decades, chronic stress (Sparrenberger et al., 2009; Liu et al., 2017) and systemic low-level sterile inflammation (De Miguel et al., 2015) are positively correlated to primary hypertension. Such conditions are normally associated with elevation of plasma norepinephrine, implicating the involvement of the nervous system, although in many cases, augmentation of epinephrine levels is also observed (; Rivier et al., 1989; Pende et al., 1990; Corssmit et al., 1996; Kannan et al., 1996; Esler et al., 2008a,b; Byrne et al., 2018). In a mouse model with genetic manipulation, TRPM4 deletion resulted in increased catecholamine secretion by AGCCs and hypertension (Mathar et al., 2010). These data collectively suggest that the catecholamine release from the adrenal glands contributes to primary hypertension.

The role of testosterone in blood pressure is controversial. Testosterone increases blood pressure in many animal studies (Masubuchi et al., 1982; Chen and Meng, 1991; Reckelhoff et al., 1998). However, hypogonadism in humans also increases the risk of hypertension (Hughes et al., 1989; Phillips et al., 1993; Jaffe et al., 1996; Liu et al., 2003; Zitzmann, 2009; Garcia-Cruz et al., 2012). In several human studies, testosterone replacement therapy increases cardiovascular disease risks (Miner et al., 2014; Anaissie et al., 2017), but the effect sizes are too small to be convincing. More importantly, such studies cannot exclude the possibility that testosterone can reduce blood pressure for a subpopulation of males with disease-associated variants in molecules of the EPHB6 signaling pathway (including EPHB6 and CHRNA3) but that such a subpopulation, being a minor one, is obscured by others in statistical analysis in studies without stratification based on testosterone levels.

Our study with stratification of testosterone levels has enhanced the association of variants of molecules in the EPHB6 signaling pathway with hypertension risks and revealed the association of a CHRNA3 variant with hypertension risks. Our results further suggest that for the subpopulation carrying the rs3743706 SNV, testosterone replacement therapy could be a personalized treatment for the cause of their hypertension. Although testosterone treatment might raise the concern of increased prostate cancer risks, this could be mitigated by routine examination of the patients for elevated serum prostate-specific antigen (PSA) levels. Moreover, a meta-analysis of testosterone replacement therapy outcome has shown no evidence of increased prostate cancer risks and no significant increase in PSA levels (Kang and Li, 2015).

AChR function, which is affected by its assembly, cycling, and clustering in addition to its opening, is known to be modulated by its associated molecules (Christianson and Green, 2004; Dobbins et al., 2008; Jones et al., 2010; Baier et al., 2011). We used three different experimental approaches (i.e., immunoprecipitation, colocalization, and FRET) to show that EPHB6 physically interacted with CHRNA3 in AGCCs. This provides a molecular basis for EPHB6 to function as a modulator of nAChR. AChR opening is the very first triggering event leading to a subsequent large Ca2+ influx, which is required for proper adrenaline secretion. Therefore, the physical interaction between EPHB6 and a dominant AChR subunit CHRNA3 implies the possibility that EPHB6 modulates AChR opening and hence the subsequent Ca2+ influx that controls adrenaline secretion.

To demonstrate that such interaction has functional relevance, we need to prove experimentally that KO AGCCs have reduced initial Ca2+/N+ influx through nAChR upon acetylcholine stimulation. We spent significant efforts to demonstrate this but were not able to succeed due to the technical difficulties of patching primary AGCCs. We also need to test in human patients with adequate testosterone levels but with loss-of-function variants of EPHB6 and/or its downstream signaling molecules whether less AGCC-derived catecholamines are produced according to 24-hour urine tests.

We previously demonstrated that EPHB6 KO in castrated mice led to hypertension (Luo et al., 2012). This prior knowledge plus our current findings from the human genetic study and mouse AGCC experiments allow us to propose the following model for the implication of the regulation of AGCC function by CHRNA3 and EPHB6 in human hypertension pathogenesis.

EPHB6 and molecules in its signaling pathways (e.g., CHRNA3) normally have a positive effect on catecholamine secretion by AGCCs. The results of our current study suggest that such a positive effect could be exerted by two different mechanisms. (1) A proper expression or function of molecules (e.g., CHRNA3) in the EPHB6 signaling pathway; or (2) an appropriate direct interaction at the protein level between EPHB6 and its signaling molecules (e.g., the physical interaction between EPHB6 and CHRNA3). Compromises in either of these two mechanisms will reduce catecholamine secretion by AGCCs, and such reduction has a protective effect against hypertension. The significant variant found in the CHRNA3 gene is related to mechanism 1. The physical interaction between EPHB6 and CHRNA3 in AGCCs implicates mechanism 2. Such a protective effect depends on the presence of sufficient testosterone, and hence only occurs in testosterone-sufficient males (Wang et al., 2018). At the same time, EPHB6 has other targets (e.g., vascular smooth muscle cells), and EPHB6 KO results in increased blood vessel contractility (Luo et al., 2012). Although we can only speculate at this time, the deleterious CHRNA3 variant might also have other targets leading to increased blood pressure. As a result, the protective effect due to the reduced catecholamine secretion is neutralized by the above-described blood increasing effects. Hence, males with deleterious variants in EPHB6 or its signaling molecules (e.g., CHRNA3) have normal BP. However, when their testosterone level is decreased due to either old ages or some pathological conditions, the protective effect due to reduced catecholamine secretion is lost. This is evidenced in our previous report that the castrated EPHB6 KO mice presented hypertension (Luo et al., 2012). Consistent with this hypothesis, we demonstrated that hypogonadic patients with deleterious CHRNA3 variants had increased hypertension risks according to our human genetic study.

Molecules in the EPHB6 signaling pathway (including EPHB6, EFNBs, AChR, BK, etc.) are also expressed in catecholamine-producing neuronal cells in the nervous system. Whether their sequence variants alone or in concert with sex hormones compromise catecholamine release and consequently modify the pathophysiology of diseases caused by dysfunctional catecholamine secretion (e.g., Parkinson’s disease) is an area worth investigating.

Data Availability Statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee of the First Affiliated Hospital of Zhejiang University (No. 2013-145). The patients/participants provided their written informed consent to participate in this study. The animal study was reviewed and approved by Animal Protection Committee of CHUM Research Center. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author Contributions

JW, HL, S-JH, and J-ZS generated the concept and initiated this project. TW, YW, WS, B-QZ, Y-MY, H-DW, Z-XX, JW, J-ZS, S-JH, JL, and RB conducted the experiments. JR and FM-B performed the genetic data analysis. YW and WS analyzed the animal data. TW, JR, YW, WS, HL, and JW drafted the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Fonds de recherche du Québec – Santé, the Natural Sciences and Engineering Research Council of Canada (RGPIN-2017-04790), the Canadian Institutes of Health Research (MOP 123389), and the J.-Louis Lévesque Foundation to JW. This study was also made possible by a group grant from the National Sciences Foundation of China (#81361120264) to J-ZS, S-JH, TW, and JW, and a grant from the National Natural Science Foundation of China (#81700363) to YW.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

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

Abbreviations

ACh, acetylcholine; AChR, acetylcholine receptors; AGCCs, adrenal gland chromaffin cells; BK channel, large-conductance calcium-activated potassium channel; BMI, body mass index; C.H.B., Han Chinese in Beijing; CHRNA3, nicotinic acetylcholine receptor alpha 3 subunit; EFNs, ephrins; EPHs, erythropoietin-producing hepatocellular kinases; KO, knockout; LD, linkage-disequilibrium; mAChR, muscarinic AChR; PSA, prostate-specific antigen; SNVs, single-nucleotide variants; VSMC, vascular smooth muscle cell; WT, wild type.

References

Adeyemo, A., Gerry, N., Chen, G., Herbert, A., Doumatey, A., Huang, H., et al. (2009). A genome-wide association study of hypertension and blood pressure in African Americans. PLoS Genet. 5:e1000564. doi: 10.1155/2018/7259704

PubMed Abstract | CrossRef Full Text | Google Scholar

Anaissie, J., Roberts, N. H., Wang, P., and Yafi, F. A. (2017). Testosterone replacement therapy and components of the metabolic syndrome. Sex Med. Rev. 5, 200–210. doi: 10.1016/j.sxmr.2017.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Axelrod, J. (1976). Catecholamines and hypertension. Clin. Sci. Mol. Med. Suppl. 3, 415s–421s.

Google Scholar

Baier, C. J., Fantini, J., and Barrantes, F. J. (2011). Disclosure of cholesterol recognition motifs in transmembrane domains of the human nicotinic acetylcholine receptor. Sci. Rep. 1:69.

Google Scholar

Batlle, E., Henderson, J. T., Beghtel, H., Van Den Born, M. M., Sancho, E., Huls, G., et al. (2002). Beta-catenin and TCF mediate cell positioning in the intestinal epithelium by controlling the expression of EphB/ephrinB. Cell 111, 251–263. doi: 10.1016/s0092-8674(02)01015-2

CrossRef Full Text | Google Scholar

Borkowski, K. R., and Quinn, P. (1984). Adrenaline and the development of genetic hypertension. J. Hypertens. Suppl. 2, S81–S83.

Google Scholar

Buhler, F. R., Amann, F. W., Bolli, P., Hulthen, L., Kiowski, W., Landmann, R., et al. (1982). Elevated adrenaline and increased alpha-adrenoceptor-mediated vasoconstriction in essential hypertension. J. Cardiovasc. Pharmacol. 4(Suppl. 1), S134–S138.

Google Scholar

Byrne, C. J., Khurana, S., Kumar, A., and Tai, T. C. (2018). Inflammatory signaling in hypertension: regulation of adrenal catecholamine biosynthesis. Front. Endocrinol. 9:343. doi: 10.3389/fendo.2018.00343

PubMed Abstract | CrossRef Full Text | Google Scholar

Campos-Caro, A., Smillie, F. I., Dominguez Del Toro, E., Rovira, J. C., Vicente-Agullo, F., Chapuli, J., et al. (1997). Neuronal nicotinic acetylcholine receptors on bovine chromaffin cells: cloning, expression, and genomic organization of receptor subunits. J. Neurochem. 68, 488–497. doi: 10.1046/j.1471-4159.1997.68020488.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y. F., and Meng, Q. C. (1991). Sexual dimorphism of blood pressure in spontaneously hypertensive rats is androgen dependent. Life Sci. 48, 85–96. doi: 10.1016/0024-3205(91)90428-e

CrossRef Full Text | Google Scholar

Christianson, J. C., and Green, W. N. (2004). Regulation of nicotinic receptor expression by the ubiquitin-proteasome system. EMBO J. 23, 4156–4165. doi: 10.1038/sj.emboj.7600436

PubMed Abstract | CrossRef Full Text | Google Scholar

Corssmit, E. P., Heijligenberg, R., Endert, E., Ackermans, M. T., Sauerwein, H. P., and Romijn, J. A. (1996). Endocrine and metabolic effects of interferon-alpha in humans. J. Clin. Endocrinol. Metab. 81, 3265–3269. doi: 10.1210/jc.81.9.3265

CrossRef Full Text | Google Scholar

De Bakker, P. I., Yelensky, R., Pe’er, I., Gabriel, S. B., Daly, M. J., and Altshuler, D. (2005). Efficiency and power in genetic association studies. Nat. Genet. 37, 1217–1223. doi: 10.1038/ng1669

PubMed Abstract | CrossRef Full Text | Google Scholar

De Miguel, C., Rudemiller, N. P., Abais, J. M., and Mattson, D. L. (2015). Inflammation and hypertension: new understandings and potential therapeutic targets. Curr. Hypertens. Rep. 17:507.

Google Scholar

Di Angelantonio, S., Matteoni, C., Fabbretti, E., and Nistri, A. (2003). Molecular biology and electrophysiology of neuronal nicotinic receptors of rat chromaffin cells. Eur. J. Neurosci. 17, 2313–2322. doi: 10.1046/j.1460-9568.2003.02669.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dobbins, G. C., Luo, S., Yang, Z., Xiong, W. C., and Mei, L. (2008). alpha-Actinin interacts with rapsyn in agrin-stimulated AChR clustering. Mol Brain 1:18. doi: 10.1186/1756-6606-1-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Dominiak, P., and Grobecker, H. (1982). Elevated plasma catecholamines in young hypertensive and hyperkinetic patients: effect of pindolol. Br. J. Clin. Pharmacol. 13, 381S–390S.

Google Scholar

Eglen, R. M. (2006). Muscarinic receptor subtypes in neuronal and non-neuronal cholinergic function. Auton. Autacoid. Pharmacol. 26, 219–233. doi: 10.1111/j.1474-8673.2006.00368.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Eph Nomenclature Committee (1997). Unified nomenclature for Eph family receptors and their ligands, the ephrins. Cell 90, 403–404. doi: 10.1016/s0092-8674(00)80500-0

CrossRef Full Text | Google Scholar

Esler, M., Eikelis, N., Schlaich, M., Lambert, G., Alvarenga, M., Dawood, T., et al. (2008a). Chronic mental stress is a cause of essential hypertension: presence of biological markers of stress. Clin. Exp. Pharmacol. Physiol. 35, 498–502. doi: 10.1111/j.1440-1681.2008.04904.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Esler, M., Schwarz, R., and Alvarenga, M. (2008b). Mental stress is a cause of cardiovascular diseases: from scepticism to certainty. Stress Health 24, 175–180. doi: 10.1002/smi.1198

CrossRef Full Text | Google Scholar

Flanagan, J. G., and Vanderhaeghen, P. (1998). The ephrins and Eph receptors in neural development. Annu. Rev. Neurosci. 21, 309–345. doi: 10.1146/annurev.neuro.21.1.309

PubMed Abstract | CrossRef Full Text | Google Scholar

Garcia-Cruz, E., Piqueras, M., Huguet, J., Perez-Marquez, M., Gosalbez, D., Peri, L., et al. (2012). Hypertension, dyslipidemia and overweight are related to lower testosterone levels in a cohort of men undergoing prostate biopsy. Int. J. Impot. Res. 24, 110–113. doi: 10.1038/ijir.2011.55

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldstein, D. S. (1983). Plasma catecholamines and essential hypertension. An analytical review. Hypertension 5, 86–99. doi: 10.1161/01.hyp.5.1.86

CrossRef Full Text | Google Scholar

Guo, Y., Tomlinson, B., Chu, T., Fang, Y. J., Gui, H., Tang, C. S., et al. (2012). A genome-wide linkage and association scan reveals novel loci for hypertension and blood pressure traits. PLoS One 7:e31489. doi: 10.1371/journal.pone.0031489

PubMed Abstract | CrossRef Full Text | Google Scholar

Hiura, Y., Tabara, Y., Kokubo, Y., Okamura, T., Miki, T., Tomoike, H., et al. (2010). A genome-wide association study of hypertension-related phenotypes in a Japanese population. Circ. J. 74, 2353–2359. doi: 10.1253/circj.cj-10-0353

PubMed Abstract | CrossRef Full Text | Google Scholar

Hughes, G. S., Mathur, R. S., and Margolius, H. S. (1989). Sex steroid hormones are altered in essential hypertension. J. Hypertens. 7, 181–187.

Google Scholar

International Consortium for Blood Pressure Genome-Wide Association Studies, Ehret, G. B., Munroe, P. B., Rice, K. M., Bochud, M., Johnson, A. D., et al. (2011). Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk. Nature 478, 103–109. doi: 10.1038/nature10405

PubMed Abstract | CrossRef Full Text | Google Scholar

Itier, V., and Bertrand, D. (2001). Neuronal nicotinic receptors: from protein structure to function. FEBS Lett. 504, 118–125. doi: 10.1016/s0014-5793(01)02702-8

CrossRef Full Text | Google Scholar

Jablonskis, L. T., and Howe, P. R. (1994). Elevated plasma adrenaline in spontaneously hypertensive rats. Blood Press. 3, 106–111. doi: 10.3109/08037059409101529

PubMed Abstract | CrossRef Full Text | Google Scholar

Jaffe, A., Chen, Y., Kisch, E. S., Fischel, B., Alon, M., and Stern, N. (1996). Erectile dysfunction in hypertensive subjects. Assessment of potential determinants. Hypertension 28, 859–862. doi: 10.1161/01.hyp.28.5.859

CrossRef Full Text | Google Scholar

Jones, A. K., Buckingham, S. D., and Sattelle, D. B. (2010). Proteins interacting with nicotinic acetylcholine receptors: expanding functional and therapeutic horizons. Trends Pharmacol. Sci. 31, 455–462. doi: 10.1016/j.tips.2010.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, D. Y., and Li, H. J. (2015). The effect of testosterone replacement therapy on prostate-specific antigen (PSA) levels in men being treated for hypogonadism: a systematic review and meta-analysis. Medicine 94:e410. doi: 10.1097/md.0000000000000410

PubMed Abstract | CrossRef Full Text | Google Scholar

Kannan, H., Tanaka, Y., Kunitake, T., Ueta, Y., Hayashida, Y., and Yamashita, H. (1996). Activation of sympathetic outflow by recombinant human interleukin-1 beta in conscious rats. Am. J. Physiol. 270, R479–R485.

Google Scholar

Konstantinova, I., Nikolova, G., Ohara-Imaizumi, M., Meda, P., Kucera, T., Zarbalis, K., et al. (2007). EphA-Ephrin-A-mediated beta cell communication regulates insulin secretion from pancreatic islets. Cell 129, 359–370. doi: 10.1016/j.cell.2007.02.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuijper, S., Turner, C. J., and Adams, R. H. (2007). Regulation of angiogenesis by Eph-ephrin interactions. Trends Cardiovasc. Med. 17, 145–151. doi: 10.1016/j.tcm.2007.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuratomi, S., Ohmori, Y., Ito, M., Shimazaki, K., Muramatsu, S., Mizukami, H., et al. (2009). The cardiac pacemaker-specific channel Hcn4 is a direct transcriptional target of MEF2. Cardiovasc. Res. 83, 682–687. doi: 10.1093/cvr/cvp171

PubMed Abstract | CrossRef Full Text | Google Scholar

Levy, D., Ehret, G. B., Rice, K., Verwoert, G. C., Launer, L. J., Dehghan, A., et al. (2009). Genome-wide association study of blood pressure and hypertension. Nat. Genet. 41, 677–687.

Google Scholar

Liu, M. Y., Li, N., Li, W. A., and Khan, H. (2017). Association between psychosocial stress and hypertension: a systematic review and meta-analysis. Neurol. Res. 39, 573–580. doi: 10.1080/01616412.2017.1317904

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P. Y., Death, A. K., and Handelsman, D. J. (2003). Androgens and cardiovascular disease. Endocr. Rev. 24, 313–340.

Google Scholar

Luo, H., Wan, X., Wu, Y., and Wu, J. (2001). Cross-linking of EphB6 resulting in signal transduction and apoptosis in Jurkat cells. J. Immunol. 167, 1362–1370. doi: 10.4049/jimmunol.167.3.1362

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, H., Wu, Z., Tremblay, J., Thorin, E., Peng, J., Lavoie, J. L., et al. (2012). Receptor tyrosine kinase Ephb6 regulates vascular smooth muscle contractility and modulates blood pressure in concert with sex hormones. J. Biol. Chem. 287, 6819–6829. doi: 10.1074/jbc.m111.293365

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, H., Yu, G., Wu, Y., and Wu, J. (2002). EphB6 crosslinking results in costimulation of T cells. J. Clin. Invest. 110, 1141–1150. doi: 10.1172/jci0215883

CrossRef Full Text | Google Scholar

Masubuchi, Y., Kumai, T., Uematsu, A., Komoriyama, K., and Hirai, M. (1982). Gonadectomy-induced reduction of blood pressure in adult spontaneously hypertensive rats. Acta Endocrinol. 101, 154–160. doi: 10.1530/acta.0.1010154

PubMed Abstract | CrossRef Full Text | Google Scholar

Mathar, I., Vennekens, R., Meissner, M., Kees, F., Van Der Mieren, G., Camacho Londono, J. E., et al. (2010). Increased catecholamine secretion contributes to hypertension in TRPM4-deficient mice. J. Clin. Invest. 120, 3267–3279. doi: 10.1172/jci41348

PubMed Abstract | CrossRef Full Text | Google Scholar

Miner, M., Barkin, J., and Rosenberg, M. T. (2014). Testosterone deficiency: myth, facts, and controversy. Can. J. Urol. 21(Suppl. 2), 39–54.

Google Scholar

Mousavi, M., Hellstrom-Lindahl, E., Guan, Z. Z., Bednar, I., and Nordberg, A. (2001). Expression of nicotinic acetylcholine receptors in human and rat adrenal medulla. Life Sci. 70, 577–590. doi: 10.1016/s0024-3205(01)01427-8

CrossRef Full Text | Google Scholar

Mousavi, M., and Nordberg, A. (2006). Expression of the alpha7, alpha4 and alpha3 nicotinic receptor subtype in the brain and adrenal medulla of transgenic mice carrying genes coding for human AChE and beta-amyloid. Int. J. Dev. Neurosci. 24, 269–273. doi: 10.1016/j.ijdevneu.2006.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Org, E., Eyheramendy, S., Juhanson, P., Gieger, C., Lichtner, P., Klopp, N., et al. (2009). Genome-wide scan identifies CDH13 as a novel susceptibility locus contributing to blood pressure determination in two European populations. Hum. Mol. Genet. 18, 2288–2296. doi: 10.1093/hmg/ddp135

PubMed Abstract | CrossRef Full Text | Google Scholar

Padmanabhan, S., Melander, O., Johnson, T., Di Blasio, A. M., Lee, W. K., Gentilini, D., et al. (2010). Genome-wide association study of blood pressure extremes identifies variant near UMOD associated with hypertension. PLoS Genet 6:e1001177. doi: 10.1371/journal.pgen.1001177

PubMed Abstract | CrossRef Full Text | Google Scholar

Pasquale, E. B. (2008). Eph-ephrin bidirectional signaling in physiology and disease. Cell 133, 38–52. doi: 10.1016/j.cell.2008.03.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Pende, A., Musso, N. R., Vergassola, C., Puppo, F., Ioverno, A., Criscuolo, D., et al. (1990). Neuroendocrine effects of interferon alpha 2-a in healthy human subjects. J. Biol. Regul. Homeost. Agents 4, 67–72.

Google Scholar

Perez-Alvarez, A., Hernandez-Vivanco, A., Mcintosh, J. M., and Albillos, A. (2012). Native alpha6beta4 nicotinic receptors control exocytosis in human chromaffin cells of the adrenal gland. FASEB J. 26, 346–354. doi: 10.1096/fj.11-190223

PubMed Abstract | CrossRef Full Text | Google Scholar

Phillips, G. B., Jing, T. Y., Resnick, L. M., Barbagallo, M., Laragh, J. H., and Sealey, J. E. (1993). Sex hormones and hemostatic risk factors for coronary heart disease in men with hypertension. J. Hypertens. 11, 699–702. doi: 10.1097/00004872-199307000-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575. doi: 10.1086/519795

PubMed Abstract | CrossRef Full Text | Google Scholar

Reckelhoff, J. F., Zhang, H., and Granger, J. P. (1998). Testosterone exacerbates hypertension and reduces pressure-natriuresis in male spontaneously hypertensive rats. Hypertension 31, 435–439. doi: 10.1161/01.hyp.31.1.435

CrossRef Full Text | Google Scholar

Rivier, C., Vale, W., and Brown, M. (1989). In the rat, interleukin-1 alpha and -beta stimulate adrenocorticotropin and catecholamine release. Endocrinology 125, 3096–3102. doi: 10.1210/endo-125-6-3096

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, W., Wang, Y., Peng, J., Qi, S., Vitale, N., Kaneda, N., et al. (2019). EPHB6 controls catecholamine biosynthesis by up-regulating tyrosine hydroxylase transcription in adrenal gland chromaffin cells. J. Biol. Chem. 294, 6871–6887. doi: 10.1074/jbc.ra118.005767

PubMed Abstract | CrossRef Full Text | Google Scholar

Slavin, T. P., Feng, T., Schnell, A., Zhu, X., and Elston, R. C. (2011). Two-marker association tests yield new disease associations for coronary artery disease and hypertension. Hum. Genet. 130, 725–733. doi: 10.1007/s00439-011-1009-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Sparrenberger, F., Cichelero, F. T., Ascoli, A. M., Fonseca, F. P., Weiss, G., Berwanger, O., et al. (2009). Does psychosocial stress cause hypertension? A systematic review of observational studies. J. Hum. Hypertens. 23, 12–19. doi: 10.1038/jhh.2008.74

PubMed Abstract | CrossRef Full Text | Google Scholar

Traish, A. M., Miner, M. M., Morgentaler, A., and Zitzmann, M. (2011). Testosterone deficiency. Am. J. Med. 124, 578–587.

Google Scholar

Tremblay, J., Wang, Y., Raelson, J., Marois-Blanchet, F. C., Wu, Z., Luo, H., et al. (2017). Evidence from single nucleotide polymorphism analyses of ADVANCE study demonstrates EFNB3 as a hypertension risk gene. Sci. Rep. 7:44114.

Google Scholar

Wang, Y., Hamet, P., Thorin, E., Tremblay, J., Raelson, J., Wu, Z., et al. (2016). Reduced blood pressure after smooth muscle EFNB2 deletion and the potential association of EFNB2 mutation with human hypertension risk. Eur. J. Hum. Genet. 24, 1817–1825.

Google Scholar

Wang, Y., Shi, W., Blanchette, A., Peng, J., Qi, S., Luo, H., et al. (2018). EPHB6 and testosterone in concert regulate epinephrine release by adrenal gland chromaffin cells. Sci. Rep. 8:842.

Google Scholar

Wang, Y., Thorin, E., Luo, H., Tremblay, J., Lavoie, J. L., Wu, Z., et al. (2015). EPHB4 protein expression in vascular smooth muscle cells regulates their contractility, and EPHB4 deletion leads to hypotension in mice. J. Biol. Chem. 290, 14235–14244. doi: 10.1074/jbc.m114.621615

PubMed Abstract | CrossRef Full Text | Google Scholar

Wellcome Trust Case Control Consortium (2007). Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 447, 661–678. doi: 10.1038/nature05911

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilkinson, D. G. (2000). Eph receptors and ephrins: regulators of guidance and assembly. Int. Rev. Cytol. 196, 177–244. doi: 10.1016/s0074-7696(00)96005-4

CrossRef Full Text | Google Scholar

Wu, J., and Luo, H. (2005). Recent advances on T-cell regulation by receptor tyrosine kinases. Curr. Opin. Hematol. 12, 292–297. doi: 10.1097/01.moh.0000166497.26397.9f

CrossRef Full Text | Google Scholar

Wu, P. C., Fann, M. J., and Kao, L. S. (2010). Characterization of Ca2+ signaling pathways in mouse adrenal medullary chromaffin cells. J. Neurochem. 112, 1210–1222. doi: 10.1111/j.1471-4159.2009.06533.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, T., Zhang, B. Q., Raelson, J., Yao, Y. M., Wu, H. D., Xu, Z. X., et al. (2018). Analysis of the association of EPHB6, EFNB1 and EFNB3 variants with hypertension risks in males with hypogonadism. Sci. Rep. 8:14497.

Google Scholar

Wu, Z., Luo, H., Thorin, E., Tremblay, J., Peng, J., Lavoie, J. L., et al. (2012). Possible role of Efnb1 protein, a ligand of Eph receptor tyrosine kinases, in modulating blood pressure. J. Biol. Chem. 287, 15557–15569. doi: 10.1074/jbc.m112.340869

PubMed Abstract | CrossRef Full Text | Google Scholar

Zitzmann, M. (2009). Testosterone deficiency, insulin resistance and the metabolic syndrome. Nat. Rev. Endocrinol. 5, 673–681. doi: 10.1038/nrendo.2009.212

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: CHRNA3, acetylcholine receptor, adrenal gland chromaffin cells, single nucleotide variant, hypogonadism, hypertension

Citation: Wu T, Wang Y, Shi W, Zhang B-Q, Raelson J, Yao Y-M, Wu H-D, Xu Z-X, Marois-Blanchet F-C, Ledoux J, Blunck R, Sheng J-Z, Hu S-J, Luo H and Wu J (2020) A Variant in the Nicotinic Acetylcholine Receptor Alpha 3 Subunit Gene Is Associated With Hypertension Risks in Hypogonadic Patients. Front. Genet. 11:539862. doi: 10.3389/fgene.2020.539862

Received: 02 March 2020; Accepted: 10 November 2020;
Published: 27 November 2020.

Edited by:

Jingyun Yang, Rush University Medical Center, United States

Reviewed by:

Chuntao Zhao, Cincinnati Children’s Hospital Medical Center, United States
Paulina Carmona-Mora, University of California, Davis, United States

Copyright © 2020 Wu, Wang, Shi, Zhang, Raelson, Yao, Wu, Xu, Marois-Blanchet, Ledoux, Blunck, Sheng, Hu, Luo and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jian-Zhong Sheng, shengjz@zju.edu.cn; Shen-Jiang Hu, s0hu0001@hotmail.com; Hongyu Luo, hongyu.luo@umontreal.ca; Jiangping Wu, jianping.wu@umontreal.ca

These authors have contributed equally to this work

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