Blood Pressure Regulation Evolved from Basic Homeostatic Components

biomedicines

Article Blood Pressure Regulation Evolved from Basic Homeostatic Components

Alon Botzer 1, Yoram Finkelstein 2 and Ron Unger 1,*

1 The Mina & Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan 5290002, Israel; [email protected] 2 Neurology and Toxicology Service and Unit, Shaare Zedek Medical Center, Jerusalem 9103102, Israel; yoram.ﬁ[email protected] * Correspondence: [email protected]

Abstract: Blood pressure (BP) is determined by several physiological factors that are regulated by a range of complex neural, endocrine, and paracrine mechanisms. This study examined a collection of 198 human genes related to BP regulation, in the biological processes and functional prisms, as well as gene expression in organs and tissues. This was made in conjunction with an orthology analysis performed in 19 target organisms along the phylogenetic tree. We have demonstrated that transport and signaling, as well as homeostasis in general, are the most prevalent biological processes associated with BP gene orthologs across the examined species. We showed that these genes and their orthologs are expressed primarily in the kidney and adrenals of complex organisms (e.g., high order vertebrates) and in the nervous system of low complexity organisms (e.g., flies, nematodes). Furthermore, we have determined that basic functions such as ion transport are ancient and appear in all organisms, while more complex regulatory functions, such as control of extracellular volume emerged in high order organisms. Thus, we conclude that the complex system of BP regulation evolved from simpler components that were utilized to maintain specific homeostatic functions that Citation: Botzer, A.; Finkelstein, Y.; play key roles in existence and survival of organisms. Unger, R. Blood Pressure Regulation Evolved from Basic Homeostatic Keywords: blood pressure; homeostasis; evolution; gene orthologs Components. Biomedicines 2021, 9, 469. https://doi.org/10.3390/ biomedicines9050469 1. Introduction Academic Editor: Juan Sahuquillo The multifactorial nature of BP regulation involves many genes with a widespread distribution across numerous cellular subsystems, posing significant challenges in the Received: 7 March 2021 effort to decipher its complex mechanisms [1]. Here we implement a comparative approach Accepted: 22 April 2021 Published: 25 April 2021 to review the phylogenetic history of the BP system via an analysis of BP associated genes and their orthologs in 19 organisms, enabling us to gain an evolutionary perspective of the

Publisher’s Note: MDPI stays neutral development of the system. with regard to jurisdictional claims in Current evidence proposes that the blood vascular system initially emerged in an published maps and institutional afﬁl- ancestor of the tripoblasts over 600 million years ago, as a means to withstand the time- iations. distance constraints of diffusion [2]. This has advanced over the course of evolution and expanded to a more complex “machinery” to support the functional requirements of high order organisms. We suggest that the most fundamental organismal function of this machinery is to preserve the internal milieu, namely, to maintain a stable internal environment concerning Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. temperature, electrolytes and water concentrations, as external conditions perpetually This article is an open access article change. Internal environment, or “milieu intérieur” is a concept formulated by Claude distributed under the terms and Bernard, postulating that “the stability of the internal environment (the milieu intérieur) conditions of the Creative Commons is the condition for the free and independent life”. This is the fundamental principle of Attribution (CC BY) license (https:// homeostasis, a term coined later by Walter Bradford Cannon. creativecommons.org/licenses/by/ Claude Bernard also stated: “The constancy of the environment presupposes a perfec- 4.0/). tion of the organism such that external variations are at every instant compensated and

Biomedicines 2021, 9, 469. https://doi.org/10.3390/biomedicines9050469 https://www.mdpi.com/journal/biomedicines Biomedicines 2021, 9, 469 2 of 18

brought into balance. In consequence, far from being indifferent to the external world, the higher animal is on the contrary in a close and wise relation with it, so that its equilibrium results from a continuous and delicate compensation established as if the most sensitive of balances” [3]. This work has utilized bioinformatic tools and data mining techniques to investigate the BP system and elucidate the major homeostatic roles it played in the physiology of the corresponding organisms. In fact, we went back, tracing the evolution of genes that are involved in human BP regulation.

2. Materials and Methods 2.1. BP Genes, Orthologous Proteins and Target Organism Selection This study investigated a set of 198 genes associated with the BP regulation system. These genes were obtained and identified by means of bioinformatic analyses performed previously by Botzer et al. [4]. A set of orthologous proteins pertaining to these genes, from 19 organisms (see Figure1) were obtained from two different sources—STRING and InParanoid. STRING—Search Tool for the Retrieval of Interacting Genes/Proteins, (http://string- db.org/ (accessed on May 2017)) is a meta-resource that aggregates most of the available information on protein–protein associations and includes both direct physical interactions and protein interactions derived from literature. Since version 9.1, it contains pre-computed orthology relations imported from the eggNOG database [5]. eggnog—evolutionary geneal- ogy of genes: Non-supervised Orthologous Groups, (http://eggnogdb.embl.de (accessed on May 2017)) is a public resource that provides Orthologous Groups (OG’s) of proteins at different taxonomic levels. It provides pairwise orthology relationships within OG’s based on analysis of phylogenetic trees. Protein sequences from the selected organisms were extracted and used to compute an all-against-all pair-wise similarity matrix. The comparison uses Smith-Waterman alignments and computational adjustments of the scores, as in BLAST, to prevent spurious hits between low-complexity sequence regions [6]. In- Paranoid (http://InParanoid.sbc.su.se (accessed on February 2018))—gathers proteomes of completely sequenced eukaryotic species and calculates pairwise ortholog relationships among them. This tool implements a two-pass BLAST approach that makes use of high- precision compositional score matrix adjustment, but avoids the alignment truncation that sometimes follows [7]. Both tools demonstrated similar results, but we decided to pursue analyses with data retrieved from STRING, which included a more comprehensive and updated coverage. We then composed an organism-to-protein similarity map, consisting of similarity scores between the 198 human BP protein sequences and their ortholog proteins in the 19 target organisms, as obtained from STRING. This map was produced by calculation of pairwise sequence alignment, using EMBOSS Needle Global Alignment tool (https: //www.ebi.ac.uk/Tools/psa/emboss_needle/ (accessed on February 2018)). This tool creates an optimal global alignment of two sequences using the Needleman-Wunsch algorithm, returning a value representing the percent of similarity between sequences, enabling homology to be inferred and the evolutionary relationship between the sequences studied [8]. Target organism selection was made with the intention to create a span of a wide variety of complexities—starting with the multicellular organism T. adhaerens as the simplest, going through flies, nematodes, fish, amphibians and eventually complex organisms as primates and other mammals (see Figure1 for a complete list of organisms). We gave emphasis to well-researched model organisms, displaying sufficient and reliable orthology data, as well as extended annotations derived from established biological and genetic knowledge. Biomedicines 2021, 9, x FOR PEER REVIEW 3 of 19

Biomedicines 2021, 9, 469 3 of 18

Figure 1. Cont. Biomedicines 2021, 9, x FOR PEER REVIEW 4 of 19 Biomedicines 2021, 9, 469 4 of 18

CACNA2D1 99.9 67.4 94.1 96.2 95.3 79.3 83.8 94.2 77 84.5 81.8 75 28.3 45.7 43 44.5 42.4 38.9 41.1 KCNK2 100 57.5 97.7 98.8 99.5 96 99.5 95.5 95.1 93.9 89 84.6 39.6 35.9 35.8 17.1 39.2 40.5 31.4 NR3C2 99.7 67.7 91.4 95.3 96.2 94.8 97.8 94 93 88.8 76.1 63.2 20.5 26 16.1 25.3 21.1 17.2 5.8 LPL 100 88.3 96.2 95.6 94.8 96.2 95.6 90 88 81 83.3 71.8 18.8 46 32.4 39.6 40.4 0 42.4 PTGS2 99.7 98.5 93 93.2 95 94.2 94.7 91.9 92.5 90.9 86.8 85.2 57.9 62.7 6.4 17.1 33.8 15.4 31 PTGS1 94.7 93.7 92.9 94 94 94.8 95.3 88.1 79.5 84.2 84.4 79.1 58.1 55.1 25.5 32.2 33.1 14.2 30.4 NR3C1 99.9 97.4 47.9 92.3 95.9 96.3 95.6 89.7 88.9 84.9 75.8 63.4 24.4 29.5 23.3 29.1 25.5 22.7 19.8 ADRA1B 97.9 99 54 96.9 87.5 97.7 74.8 85.6 81.7 63 64.2 60.7 44.2 34.9 40.8 33.8 37.4 33 30.1 SCNN1G 90 93 92.9 91.9 92.2 91.4 93.6 86.6 87.1 80.4 74.5 33.4 38.7 32 39 33.9 28.9 34.7 35.1 SCNN1B 87.9 98.8 90.8 90.9 94.1 91.3 93.9 87.3 87.4 76 75 32.9 38.4 33 38.3 33.7 29.4 33.1 38.6 EDNRB 99.4 96.1 77.3 77.5 77.3 76.4 80.7 72 61.3 67.3 60.9 61 36.1 29.8 32.9 33.6 35.7 30.4 27.2 DRD5 99 97.1 87.2 87.2 91.2 89 80.9 81.5 47.6 77 73.7 64 45.6 34.6 43.7 39 42.3 36.7 35.1 F2R 99.8 98.4 81.2 87.9 88.6 87.4 91.6 82.4 43.6 71 70.2 65.4 35.5 33.7 37.1 29.9 33.8 32.8 38.6 WNK4 99.5 89 52.4 88.8 53.4 90.1 83.7 80.9 31.9 53.3 37.3 39 17.5 14.8 22.8 23.2 14.7 32.3 19.3 WNK1 99.3 97.7 77.6 77.3 92.3 86.9 84.3 83.5 59.9 72.5 62.2 53.1 9.6 8 36.7 35.5 10.3 35.9 10.1 PDX1 99.3 98.6 90.1 89.4 93.6 93 95.4 78.3 62.4 66.8 59.6 55.9 24.7 28.8 37.5 23.3 18.5 28.5 33.1

row III row SELE 99.5 97.7 75 81.4 86.1 70.9 88 42.3 72.1 45.4 48.1 43.1 36.8 13.4 32.3 22.5 23.5 23.3 16 SCNN1A 62.4 95.9 83.5 82.8 77 75.6 81.1 68.2 67.8 64.1 57.4 35.4 33.8 33.2 35.5 28.4 29.1 31.2 30.3 GIPR 60.2 85.7 83.4 83.4 88.8 89.7 90.6 72.8 52.3 53.8 57.3 59.4 45.2 45.1 37.5 37.5 37.4 39.8 24.9 HP 81.3 89.9 77.3 76.4 77.1 85.7 77.1 72.4 27.7 30.7 43.1 26.9 43 33.1 25.4 30.1 30.4 28.5 28.8 PTGER3 86.8 80.9 74 76.6 44.8 77.7 85.7 65.6 44.3 65.1 68.1 66.4 47.4 31.9 34 34 28.7 0 34.5 PTGIR 97.9 97.4 77.2 78.1 88.1 87.4 80.1 61.9 45.3 43.4 57.8 57.2 43.8 33.4 30.4 30.6 31 0 32.8 AGT 99 95.5 73.7 71.5 77.4 60 80.6 67.3 62.4 56.2 51.9 46.2 32.2 34.4 34.5 32.4 32.6 0 32.7 HSD11B1 100 97.9 7.4 87.7 71.2 91.8 91.9 43.8 70.7 72.8 69.5 60.9 56.2 52.8 49.5 44.1 40.1 43.5 40.9 CLCNKB 99.1 91.7 88.8 90 92.4 91.6 87.5 42.2 72.2 73.5 55.3 58.5 47.3 53.5 27.5 32.6 35.1 39.8 50.6 CLCNKA 97.2 86.9 89.4 90.8 93 91.7 86.6 42.2 71.8 73.1 54.9 49.8 47.2 52.7 27.5 32.8 35 40.3 51.6 LRP8 99.2 97.6 89.4 82.7 87.4 78.2 87.5 64.6 54.1 70.3 60.7 62.9 44.7 35.9 10.8 42.4 44.2 42.2 39.4 CMA1 98.8 98.4 80.2 77.3 86.3 80.5 66.4 71.4 64.8 56.7 51.6 54.3 47.2 48.7 10.4 44.5 45.1 39.8 47.7 CTSG 98.4 92.2 78.1 77.6 75 72.4 72.6 70 62.6 58 49.8 56.5 45.4 49.8 6 37 14.3 36.4 46.3 SERPINA4 99.5 98.8 73.3 59.4 81.2 70 63.1 60.3 66.6 63.9 61.2 54.6 47.4 46.3 47.3 45 44.2 42.1 40.5 CYP11B2 79.2 97.4 83.1 83 85.5 84.1 88.1 60.7 66.7 57.3 56.6 48.6 47.2 37.9 43.2 44.5 41.4 43.3 45.9 CYP11B1 82.1 95.8 81.5 81.2 84.5 82.1 86.3 60.5 66.3 56.4 58 48 48.9 38.4 42.3 44.6 42.2 42.8 40.6 Conservation HSD11B2 100 97.3 90.4 88.1 70.6 88.9 73.1 48.9 64 55.9 59.3 54.1 40.3 38.3 37.2 35.2 33.9 40.6 39.7 AVPR2 99.7 93 89.8 89.8 89.8 87.6 86.5 54.6 65.4 58 56.3 56.3 46.5 41.1 42.1 34.8 30.8 43.9 35 CA4 99.4 91.1 74.4 71.8 83.1 79.6 80.8 43.9 38.2 53.7 46.7 52.5 47.6 43.9 44.8 43.9 36.6 42.9 38.5 ADRA2C 54.7 69.6 93.7 93.5 90 58.3 57.7 71.7 37.4 71.5 69.7 65.8 43.4 48.3 39.3 38.4 41.6 45.2 24.7 ADRA1D 70.1 46.9 88.6 88.1 62 83.8 47.6 80.1 56.4 60.2 48.9 64 37.7 29.5 39.7 32.7 36.7 33.5 25.4 DRD4 70.2 92.6 77 77.3 47.5 48.9 45.7 53.9 43.4 58 52.7 62.9 37.2 34.9 36.5 37.3 37.2 43.9 27.3 KLK1 54.2 73.7 78.2 81.7 25.6 85.9 80.9 59.7 59.9 55.8 56.1 56.5 50.8 52.7 39 51.3 47.8 43.4 50.2 GOSR2 86.2 96.3 82.5 80.3 82.1 83.4 83 86 0 78.6 76.4 83.2 69.1 65.8 51.3 53.8 0 52.7 55.2 RENBP 99.5 95.1 89.1 89.6 90.5 86.1 90.6 81.3 63.8 0 74.8 70.8 67.4 60.5 66.5 0 0 0 0 PNMT 99.683.189.886.189.892.678.987.27867.370.753.149.5000040.30 MAS1 99.7 99.4 95.4 95.1 45.4 92.9 53.9 92.9 81 74.5 54 44.7 35.4 39.7 37.3 0 0 0 0 PTGER2 98.9 98 88.7 91.2 88.9 89.7 81.6 75.3 72.3 39.7 64.5 55 32.4 37.9 32.1 0 0 0 36.8 TBXA2R 99.3 90.2 68.8 68.3 77.1 75.9 46.3 70.1 65.1 43.7 58.3 44.4 43.4 30 29.8 0 0 0 0 IGF1 99 99 71.9 72.4 88.3 87.8 59.2 61.5 64.3 64.5 62.4 57.1 17.6 26.9 20.8 0 0 0 0 PLA2G1B 98.6 96.6 85.1 85.8 89.9 48 87.2 73.2 62 70.3 74.7 65.8 51.2 49.4 56 16.9 12.2 48.1 34.6 COMT 99.6 96.3 87.2 87.6 88.9 38 73.4 70.1 70.3 73.1 69.6 65.1 62.9 59.9 51.2 0 0 0 0 NPR3 99.4 99.1 93 93 93.7 95.6 96.5 83 40.3 83.4 74 65 17.3 40.2 15 18.6 19.3 21.2 17.7 VEGFC 99 96.9 91.2 91.7 94.8 92.2 83.1 73.9 48.8 83.8 73.1 68.1 32.4 0 0 25 0 0 0 CACNG1 10090.187.990.691.187.994.288.886.271.458.180.522.1000000 NPY 9910094.995.9999910084.590.791.888.777.30000000 AVP 96.3 60.2 84.5 85.1 91 91.6 68.3 80 73.1 68.3 67.9 61.1 39.7 0 33.2 0 0 0 0 APOA1 10032.575.780.592.587.387.666.745.866.361.347.60000000 PTHLH 84.183.3087.693.393.893.88784.280.960.3500000000 GNAS 99.5 83.1 31.8 72.4 80.6 33.4 33.7 31.9 30.9 34.1 34.5 34.2 29.2 28.9 30.9 30.6 31.1 28.9 30.7 SCNN1D 57.4 66.7 43.3 41 55 53.5 50.5 45.2 35.4 42.6 42.7 28 28.6 29.1 29.1 26.7 25.8 28.1 27.2 IGF2 78.8 23.3 70.8 69.5 81.6 67.1 65.8 24.8 28 53.7 54 51.9 25.8 28.2 17.3 0 0 0 0 GH2 55.357.144.544.544.841.645.542.8040.748.732.70000000 INS 99.183.687.385.530.181.888.20068.264.661.813.3000000

V SLC14A2 99.2 97.9 90.3 91.6 93.3 88.4 91.9 37.2 34.8 72.5 34.5 19.5 23.2 0 22.2 0 0 0 0 I LEPR 99.292.185.885.29189.590.859.464.464.159.341.322.4000000 w

o CALCA 95.792.9797984.470.882.343.860.462.24757.10000000 r EDN2 99.497.276.478.172.579.461.86641.656.940.148.90000000 ADM 10097.376.672.686.585.455.173.451.764.454.446.30000000 EDN1 98.194.876.177.57778.480.86862.255.453.250.70000000 GH1 95.489.976.67878.476.778.974.147.873.666.452.20000000 GCG 10095.19593.995.663.995.690.734.777.253.253.10000000 IAPP 97.896.673.175.389.979.187.676.452.148.90450000000 GAL 10075.680.682.376.885.48776.644.872.4069.10000000 VEGFB 100 98.6 60.1 90.3 45.5 87.4 54.8 45.3 46.2 39.1 35.9 40.6 18 0 0 0 34.2 0 0 APOE 95.97681.282.482.283.923.8625737.55050.90000000 VEGFA 21.693.285.284.771.360.5506438.444.829.534.229.9000000 NPPA 98.794.88584.386.188.289.578.37450000000000 LEP 10095.889.89172.190.49480.268.90000000000 NPPC 99.298.494.492.99693.774.683.127.804953.80000000 APOA2 999874.578.49285908000000000000 ICAM3 98.992.735.735.161.972.770.54852.4034.523.520000000 ICAM1 98.334.864.567.276.469.773.449.452.2035.424.60000000 KNG1 6295.769.167.750.880.28252.212.637.230.323.28.9000000 TRH 97.593.461.961.983.579.585.251.72444.136.144.50000000 EDN3 99.694.151.258.360.247.649.851.927.731.435.725.20000000 ICAM2 98.989.871.270.765.461.966.962.42021.421.832.90000000 BSND 90.99576.976.982.380.78222.816.80029.60000000 APOC4 95.389.871.769.368.572.474000000000000 APOC2 9991.176.276.251.980.281.2058.20000000000 APOC3 82.190.970.374.78173.3059.200000000000 NPPB 98.567.90061.767.267.2046.5045.700000000 FigureFigure 1. 1.HeatmapHeatmap displays displays thethe hierarchical hierarchical clustering clustering of similarity of similarity percentage percentage scores betweenscores between human BP human genes (verticalBP genes axis) (vertical and axis)their and orthologs their orthologs in 19 organisms in 19 (horizontal organisms axis). (horizontal Green color axis). denotes Green the morecolor conserved denotes genesthe more across conserved the 19 organisms, genes aligned across at the the 19 organisms,top of the aligned map (see at conservation the top of scalethe map arrow (see to theconservation right of map). scal Ae phylogenetic arrow to the tree right describing of map). the A developmental phylogenetic order tree describing is placed theabove developmental the map. order is placed above the map.

2.2. Creation of a Similarity Map and Hierarchical Clustering STRING—Search Tool for the Retrieval of Interacting Genes/Proteins, (http://string- db.org/ (accessedThe organism-to-protein on May 2017)) similarity is a meta-resource map enabled us that to implement aggregates a two-way most of hierarchi- the available cal clustering on the similarity matrix scores, also known as hierarchical cluster analysis. informationHierarchical on protein–protein clustering is an associations algorithm that and recursivelyincludes merges both direct objects physical based on interactions their andpair-wise protein distance.interactions Neighboring derived objectsfrom literature. are merged Since ﬁrst, while version objects 9.1, farthest it contains apart arepre-com- putedmerged orthology last. The relations ultimate imported result is a from set of the clusters, eggNOG where database each cluster [5]. iseggnog—evolution- distinct from ary eachgenealogy other cluster, of genes: and Non-supervised the objects within Orthologous each cluster are Groups, considerably (http://eggnogdb.embl.de similar to each (accessedother. on The May main 2017)) output is ofa public hierarchical resource clustering that provides is a dendrogram, Orthologous which Groups shows the(OG’s) of proteinshierarchical at different relationship taxonomic between levels. the clusters. It provi This calculation provided us with clusters of organisms on one axis, and of BP proteins in the other axis, suggesting many interesting des pairwise orthology relationships within OG’s based on analysis of phylogenetic conclusions relating to the evolutionary order by which the circulatory system has evolved. trees.Hierarchical Protein sequences clustering from was performed the selected using organisms the Broad Institutewere extracted Morpheus and tool used (https:// to compute an all-against-allsoftware.broadinstitute.org/morpheus/ pair-wise similarity (accessedmatrix. onThe May comparison 2018)), implementing uses Smith-Waterman Euclidian alignmentsdistance metricand computational and the Complete adjustments linkage method of onthe both scores, axes, as emphasizing in BLAST, the to maximum prevent spurious ofhits the between-clusterbetween low-complexity dissimilarities [9sequence]. regions [6]. InParanoid (http://InPara- noid.sbc.su.se (accessed on February 2018))—gathers proteomes of completely sequenced eukaryotic species and calculates pairwise ortholog relationships among them. This tool implements a two-pass BLAST approach that makes use of high-precision compositional score matrix adjustment, but avoids the alignment truncation that sometimes follows [7]. Both tools demonstrated similar results, but we decided to pursue analyses with data retrieved from STRING, which included a more comprehensive and updated coverage. Biomedicines 2021, 9, 469 5 of 18

2.3. Functional and Tissue Expression Enrichment Analyses For a detailed in-depth analysis, we chose ﬁve organisms that essentially reﬂect the entire span of complexities of the 19 organisms that appear in the heatmap—Homo sapiens, Mus musculus, Danio rerio, Drosophila melanogaster, and C. elegans. For each organism, we investigated the resulting clusters of orthologous proteins that were obtained from the hierarchical cluster analysis. These analyses consisted of functional GO term annotations (biological processes) as well as organ/tissue expression enrichment patterns. This was carried out by a set of tools, some of which specialize in human and mouse data, while others focus on other organisms. Table1 describes the list of the tools we utilized.

Table 1. Tissue enrichment and functional annotation tools utilized in this work.

Tool Link to URL Content Expression-based matching algorithm, provides gene https://ga.genecards.org/ (accessed GeneAnalytics annotations and enrichment in each specific tissue or organ in on February 2018) human or mouse. Provides a computational representation regarding the geneontology.org/ (accessed on GeneOntology functional enrichment of genes (proteins produced by genes) February 2018) from a variety of organisms from human to bacteria. Finds biological categories enriched in gene lists, functionality, https://biit.cs.ut.ee/gprofiler/gost g:Profiler conversions between gene identifiers and mappings to their (accessed on October 2018) orthologs. https: Defines tissue-specific genes using RNA-Seq data from the TissueEnrich //tissueenrich.gdcb.iastate.edu/ Human Protein Atlas, GTEx, and mouse ENCODE. (accessed on October 2018) https://bgee.org/ (accessed on A database used to retrieve and compare gene expression Bgee November 2018) patterns in multiple animal species. http: TSEA //genetics.wustl.edu/jdlab/tsea/ Tissue specific expression analysis for human and mouse. (accessed on October 2018) http://www.mousemine.org/ Provides queries for anatomy and gene ontology enrichments MouseMine mousemine/begin.do (accessed on as well as other analyses for mouse. November 2018) http: Searches for zebrafish related biological data including genes, ZebrafishMine //www.zebrafishmine.org/begin.do proteins, tissues, and more. (accessed on November 2018) http://zeogs.scienceontheweb.net/ ZEOGS index.html (accessed on November Predicts sites of preferential expression of zebrafish gene sets. 2018) Integrates many types of data for drosophila and other http://www.flymine.org/ (accessed FlyMine organisms. runs flexible queries for genes, proteins, ontology on November 2018) terms, and tissue enrichments. http://flybase.org/convert/id Curates and organizes a diverse array of genetic, molecular, FlyBase (accessed on December 2018) genomic, and developmental information in drosophila. https://www.wormbase.org/tools/ Provides information concerning the genetics, genomics, and WormBase enrichment/tea/tea.cgi (accessed on biology of C. elegans and related nematodes. December 2018) http://wormexp.zoologie.uni-kiel. Integrates all published expression data for C. elegans, runs WormExp de/wormexp/ (accessed on flexible queries for gene enrichment analyses on selected December 2018) gene sets.

GO term enrichment tools (primarily GeneOntology) are used to evaluate characteristics of sets of genes by comparing the frequency of GO terms in the sample gene set with the frequency of the same set of GO terms in a reference set, usually a whole genome. The tools apply the binomial test to identify over or under-represented terms in the sample gene set compared to the reference genome set. The default parameters also apply a Bonferroni correction for multiple comparisons. A similar principle is implemented in the tissue enrichment analyses tools utilized for the different organisms, as speciﬁed in Figure2. Biomedicines 2021, 9, 469 6 of 18

Since GO is roughly hierarchical, with “child” terms being more specialized than their “parent” terms, we were compelled to elect the terms in the 4–5 orders of the hierarchy in each cluster-species, in order to obtain broad deﬁnitions that enable comparison between the vast span of organisms in the analysis. As a control, the same GO term functional analysis was performed 30 times on randomly selected sets of 200 genes. This random Biomedicines 2021, 9, x FOR PEER REVIEW 8 of 19 analysis has not yielded any signiﬁcant enrichment results, emphasizing the validity of the results for the set of blood pressure genes.

Figure 2. Demonstrates a summary of functional analysis results overlaid on heatmap of Figure1 1.. InIn eacheach box,box, resultsresults pertain to the enriched biological processes for the respective genes genes-organisms-organisms of the cluster. For ease of interpretation, the colored text represents biological process families, reflecting their distribution throughout the map. For full functional the colored text represents biological process families, reﬂecting their distribution throughout the map. For full functional analysis results please refer to Supplementary Table S1. analysis results please refer to Supplementary Table S1. Tissue and organ expression enrichment analyses reveal that by and large genes in column3. Results I clusters are abundant in the kidney and somewhat in the adrenal gland, pan- creasThe, and heatmap liver. In column representing II, gene the expression results of is theevident two-way mainly clustered in the kidney similarity and nervous map is system,presented and in slightly Figure1 in. This the pancreas. map has yielded This expression three clusters pattern in thechanges organism significantly axis (columns for col- I umnto III) III, and where four gene clusters expression in the gene is clearly axis (rows noticeable I to IV). exclusively in nervous system tissues. Figure 3 demonstrates a summary of tissue expression for the various organisms investigated. For full expression analysis results please refer to Supplementary Table S2. Alto- gether, the functional and expression results of the analysis in zebrafish display the summation of results of analyses in the higher and lower order species. This reflects the unique spot of zebrafish in ontogeny that well supports the rationale to exploit zebrafish as a model organism in biological research [10]. Biomedicines 2021, 9, 469 7 of 18

The order by which the organism axis has aligned as a result of the clustering process is consistent with the natural known developmental order: (1) column I—high complexity mammals—from Pan troglodytes to Equus caballus, (2) column II—medium complexity— from O. anatinus to Danio rerio, (3) column III—low complexity—B. floridae to T. adhaerens. This order, presented in Figure1 as a phylogenetic tree, is consistent with the developmental order of the blood circulatory mechanism in respect to the affiliated genes; 4—chamber heart system in high-order mammals and avians, 3 and 2—chamber heart in fish and amphibians, diffusion and hemolymph-based systems in insects and simple multi-cellular organisms. The gene axis was aligned in four clusters, highlighting the more conserved genes across the 19 organisms at the top of the map (row I). The least conserved genes are shown at the bottom of the map (row IV), depicting a majority of zero value similarity scores, due to absence of orthologues genes in low complexity organisms. It is also amenable to discern that row I harbors genes associated with more basal functions, as ion transport, which is assumed to be ancient, hence appears in all organisms; while genes that are linked to a more complex regulatory functions, as control of extracellular volume or endothelial vasoconstriction appear in row IV, hence emerge in more recent organisms only. Functional gene ontology (GO) term annotation (a statement regarding the function of genes) analysis demonstrates that homeostasis is the most prevalent biological process across most heatmap clusters. Transport processes are significant in row I, while blood circulation and vasculature structure are prevalent in columns I and II. Signaling processes appear most in row IV of columns I and II, rows I—III of columns II and III. Response to stimulus is mainly displayed in column II and slightly in columns I and III. Figure2 demonstrates a summary of the functional analysis results for the various organisms investigated, displaying the enriched biological processes for each cluster. For full functional analysis results please refer to Supplementary Table S1. Tissue and organ expression enrichment analyses reveal that by and large genes in column I clusters are abundant in the kidney and somewhat in the adrenal gland, pancreas, and liver. In column II, gene expression is evident mainly in the kidney and nervous system, and slightly in the pancreas. This expression pattern changes significantly for column III, where gene expression is clearly noticeable exclusively in nervous system tissues. Figure3 demonstrates a summary of tissue expression for the various organisms investigated. For full expression analysis results please refer to Supplementary Table S2. Altogether, the functional and expression results of the analysis in zebrafish display the summation of results of analyses in the higher and lower order species. This reflects the unique spot of zebrafish in ontogeny that well supports the rationale to exploit zebrafish as a model organism in biological research [10]. Biomedicines 2021, 9, x FOR PEER REVIEW 9 of 19 Biomedicines 2021, 9, 469 8 of 18

Figure 3. Demonstrates a summary of expression analysis results overlaid on heatmap of Figure 11.. InIn eacheach box,box, resultsresults pertain to the enriched organ/tissue expression for the respective genes-organisms of the cluster. For ease of interpretation, pertain to the enriched organ/tissue expression for the respective genes-organisms of the cluster. For ease of interpretation, colored text represents tissue/organ types, reflecting their distribution throughout the map. For full functional analysis colored text represents tissue/organ types, reflecting their distribution throughout the map. For full functional analysis results please refer to Supplementary Table S2. results please refer to Supplementary Table S2. 4. Discussion 4. Discussion Regulation of BP is a complex systemic mechanism due to numerous physiological elements,Regulation involving: of BP pressure is a complex-volume systemic regulation, mechanism which is due tightly to numerous related to physiological pressure-na- triuresiselements, [11] involving:, rapid control pressure-volume of vessel resistance regulation, by the which central is nervous tightly relatedsystem to(CNS), pressure- spe- cificallynatriuresis modulated [11], rapid by controlboth sympathetic of vessel resistance and parasympathetic by the central nervous nervous system system—the (CNS), two branchesspecifically of the modulated autonomic by nervous both sympathetic system (ANS) and parasympathetic[12], neurotransmitters nervous (e.g. system—the, noradren- two branches of the autonomic nervous system (ANS) [12], neurotransmitters (e.g., no- aline (NA), adrenaline) and hormones (e.g., angiotensin) as well as the long-term activity radrenaline (NA), adrenaline) and hormones (e.g., angiotensin) as well as the long-term of the renin-angiotensin-aldosterone-system (RAAS). These mechanisms are efficient at activity of the renin-angiotensin-aldosterone-system (RAAS). These mechanisms are effi- maintaining BP within a normal physiological range at Systolic BP of <120 mmHg and cient at maintaining BP within a normal physiological range at Systolic BP of <120 mmHg Diastolic BP of <80 mmHg [13]. and Diastolic BP of <80 mmHg [13]. This study examined a collection of 198 human genes related to BP regulation in the This study examined a collection of 198 human genes related to BP regulation in biological processes and functional prisms, as well as gene expression in organ and tissue the biological processes and functional prisms, as well as gene expression in organ and perspective. In this context, it is unavoidable that the question of homeostasis will arise. tissue perspective. In this context, it is unavoidable that the question of homeostasis This observation, in conjunction with the orthology analysis we performed in 19 target will arise. This observation, in conjunction with the orthology analysis we performed in organisms along the phylogenetic tree, evokes numerous interesting conclusions and fur- 19 target organisms along the phylogenetic tree, evokes numerous interesting conclusions ther questions. and further questions. HomeostasisHomeostasis sustains sustains the the “ “internalinternal milieu” milieu” of of an an organism’s cells, tissues, tissues, organs organs,, andand whole whole body, body, within within limits limits that that are are compatible compatible with with survival. survival. This This internal internal milieu milieu reflects the composition of the primordial ocean, first and foremost its electrolyte and water Biomedicines 2021, 9, 469 9 of 18

reflects the composition of the primordial ocean, first and foremost its electrolyte and water balance [14,15]. Homeostasis in homeothermic land-dwelling organisms is profoundly developed since their fluid balance, blood pH, oxygen tension and particularly their body temperature must be maintained within tight boundaries amid all conditions of their life cycle and in all their habitats. Homeostasis will tend to stabilize BP, maintaining it at a steady state. The term stress defines any stimulus or succession of stimuli of such magnitude as to tend to disrupt the homeostasis of the organism. When mechanisms of adjustment fail or become incoordinate or disproportionate, the stress should be referred to as an insult. The response of the organism to stress depends largely on biochemical and physiological homeostatic mechanisms [16]. One of the most important mechanisms for maintaining homeostasis is the negative feedback system; if a physiological disturbance occurs, the body will counteract the disturbance via a negative feedback mechanism and attempt to return the body to its normal steady state. The more complex the organism, the more complex its adaptation to the ever-changing environmental conditions. Biological systems may undergo functional plasticity in the effort to adapt. This is the mainstay of evolutionary progression [17]. The results of this study may elucidate the diversity of the modulating response mechanisms during phylogenesis: from the response of primitive organisms to physical stimuli (such as osmotic or heat stressors) to the complex endocrine and neurobehavioral stress response to both physical, mental, and cognitive challenges in human.

4.1. The Function of the CNS in the Maintenance of Homeostasis The CNS is the organ that orchestrates the response to stressful stimuli through the autonomic, neuroendocrine, and immune systems, as well as through neurobehavioral responses such as fight-or-flight response. Adaptation to stress and to changing environmental conditions involves neural and humoral mediators: neurotransmitters and neuromodulators, hormones and cytokines of the immune system. The goal of this adaptation is to maintain homeostasis and promote survival of the organism [18]. The stress responses include primarily the activation of the mesolimbic system, the hypothalamo-pituitary-adrenal (HPA) axis and the ANS. The dopaminergic (DA-ergic) activity in the septo-hippocampus is a mainstay of the mesolimbic stress response, due to its close relationship with the hypothalamus (autonomic control), midbrain reticular formation (arousal) and multiple sensory pathways (attention). A DA-ergic pathway originates in the ventro-medial tegmentum of the midbrain and projects to the lateral septal nuclei with abundant indirect interconnections with the cholinergic nucleus of the medial septum, forming a final common pathway of the neural, psychic and endocrine stress response. Indeed, time honored studies corroborate the major importance of catecholamines in the ANS and of the glucocorticoids (GCs) in the adrenal cortex, as well as the CNS septo- hippocampal DA-ergic and cholinergic systems in the stress response. These complex mechanisms of adaptation are critical in the control of homeostasis [17,19]. Additionally, the mesolimbic pathway transports dopamine (DA) from the ventro-medial tegmentum to the nucleus accumbens and amygdala. The nucleus accumbens is found in the ventral medial portion of the striatum and is believed to play a role in reward, desire, and the placebo effect. The amygdala is a key component of the limbic system and is associated with emotion. Unabated DA-ergic stimulation (the absence of negative feedback) has been postulated to be associated with disorders such as binge eating or drug addiction [20,21]. The neuroendocrine system integrates the functions of two major control systems: the nervous system and the endocrine system. Thus, both internal and external fluctuations are monitored by the nervous sense organs, the neural signals of which are processed by the CNS and converted to endocrine outputs. Furthermore, associated with the neuroendocrine system are distinct circumventricu- lar fenestrations in the blood-brain barrier, mainly the chemoreceptor trigger zone (CTZ) in the area postrema of the medulla oblongata. The CTZ allows hormones and neurotrans- Biomedicines 2021, 9, 469 10 of 18

mitters to enter directly from the blood to the cerebrospinal fluid (CSF), hence enabling the feedback mechanisms with the neurosecretory system itself. The hypothalamus plays a principal role in the neuroendocrine control system, which together with its connections to the pituitary gland comprise the hypothalamic pituitary system—a neurosecretory system that produces releasing hormones (or releasing factors), emitting them to the pituitary gland. Important pituitary hormones are the adrenocorticotropic hormone (ACTH), which controls the adrenal cortex as well as growth hormone- inhibiting hormone (GHIH), which regulates the endocrine system and affects neurotrans- mission. Moreover, upper neurons of the ANS are located in the hypothalamus and control the adrenal medulla, an additional crucial endocrine gland. These roles of the CNS in the maintenance of homeostasis are dynamic, empowering the organism to function efficiently. As seen in the Results section and in Figure3, the adrenal gland and pancreas are strongly pronounced in human and mouse in rows II–IV, while the CNS and brain manifest rows I–III in fruitfly and roundworm. Interestingly, in the zebrafish gene enrichment in the CNS and the brain are prominent in all clusters, whereas the pancreas appears in rows III–IV, emphasizing the clear homology of these anatomic properties among the various species. This may suggest that organs that are enriched in expression of highly conserved genes across species reflect shared ancestry.

4.2. The Role of the Adrenal Glands The adrenal gland comprises two endocrine tissues that differ in function and em- bryonal origin: The adrenal cortex, a steroidogenic tissue that evolved from the coelomic epithelium; and the adrenal medulla, a catecholamine producing tissue composed of chromaffin cells. The medulla derived from the neuroectoderm, which migrated to the adrenal gland blastema. The secretion of the adrenal cortex hormones is controlled by ACTH in a closed loop feedback system. It produces numerous steroids secreted in widely varying amounts; these are classified as mineralocorticoids (MCs) and GCs including sex steroids. The main MC is aldosterone, which is the most potent steroid affecting active transport of sodium ions across membranes and thus is crucially important in maintaining electrolyte balance and as a result, in BP homeostasis. Cortisol is the principal GC in human and most mammalian primates and is quantitatively the main secretory product of the adrenal cortex. The medulla is basically a modified sympathetic ganglion that converts tyrosine into DA, NA, and adrenaline, which are secreted in postsynaptic response to direct neural inputs. These inputs originate in the CNS (cortical cerebral areas) and are transmitted to the adrenal medulla via cholinergic preganglionic sympathetic neurons [22]. As noted, the primordium of the adrenal cortex coincides with the mesonephric blastema. It develops within the kidney in fish and on the kidney in amphibians. The kidneys of both fish and amphibians are of mesonephric origin [23]. In fish, chromaffin and cortical adrenal cells are not placed in proximity, but rather intermingle. The interrenal cells, that correspond to the adrenal cortex in mammals, from which GCs and MCs are secreted, do not encapsulate the chromaffin cells, which are the functional equivalent of the adrenal medulla. NA is the main catecholamine stored by the chromaffin cells, albeit variable quantities of adrenaline can also be present. Amphibians show a closer relationship between cortical and chromaffin cells. The adrenals of amphibians lay on the ventral side of the kidneys, as the mesonephros—the analogue of the mammalian kidney, preserving the contact between the two organs during ontogenesis [24]. In mammals, where the evolution of the excretory system leads to the development of the metanephros, adrenal cortical cells completely surround the chromaffin cell mass— namely the medulla—and are grouped together to form the adrenal gland. Traces for this can be seen in our detailed expression analysis results (Supplementary Table S2), where row I–II are expressed in the adrenal glands of mammals, and in chromaffin cells and mesonephros of the zebrafish. Biomedicines 2021, 9, 469 11 of 18

4.3. Nervous Regulation of the Cardiovascular System The circulation is regulated partly by intrinsic and local mechanisms in each tissue and prominently by the nervous system, the action of which is extremely rapid and comprehensive. All the vascularized tissues in vertebrates are supplied with sympathetic nerve fibers. Sympathetic stimulation to the arterioles increases the resistance, thus changing the blood flow through the tissues. Sympathetic stimulation of the larger vessels decreases their blood volume. In a similar manner, cardiac output and BP are remarkably increased by increased sympathetic activity. Both heart rate and stroke volume (the two components of cardiac output) alongside vascular resistance are the main determinants of BP. The internal environment is monitored by sense organs and organelles: chemoreceptors sensitive to the partial oxygen pressure in the arterial blood, mechanoreceptors sensitive to blood pressure, and chemoreceptors within the central nervous system itself sensitive to hydrogen ion concentration or to various hormones. The input from the sense organs and organelles is transmitted to the CNS, in which it is processed and from which the appropriate outputs are sent to the effectors—muscles and glands. Above all, the vasomotor center (VMC)—in fact a network of neurons within the medulla oblongata, regulates BP and other homeostatic functions by a dual action: vasoconstriction and vasodilation. Tonic adrenergic discharge to the arterioles sustains arterial pressure, while fluctuations in this tonic discharge constitute the mechanism affecting the feedback regulation of BP carried out by the carotid sinus (baroreceptor) and the carotid body (chemoreceptor that modulates the cardiovascular and respiratory systems via sympathetic tone). The effects of this discharge are of substantial importance in the preparation of the organism to endure emergency.

4.4. Catecholamines and BP Homeostasis Dopamine (DA), noradrenaline (NA), and adrenaline (known as catecholamines) are physiologically active molecules, acting both as hormones and neurotransmitters, crucial for the sustainment of homeostasis via the ANS. This work considered DA receptor genes (DRs), (DR1-like subtype—DRD1, DRD5, and DR2-like—DRD2, DRD3, DRD4) and adrenergic receptor genes (ARs), (adrenoceptors type alpha—ADRA1A, ADRA1B, ADRA1D, ADRA2A, ADRA2B, ADRA2C, and beta—ADRB1, ADRB2, ADRB3) appearing in rows II–III, all which play key physiological roles in the nervous and cardiovascular systems, specifically in BP regulation. Almost all vasomotor nerves are adrenergic. The alpha-adrenoceptors are preeminent in innervation of vascular smooth muscles and also in the lower urinary tract. In the myocardium, beta1-adrenoceptors predominate and stimulate the rate and contractility. In both cases, neurotransmitter NA exerts its physiologic effects by binding to alpha ARs, whereas adrenaline reacts with both alpha and beta ARs, causing vasoconstriction and vasodilation, respectively. Generally, it is the alpha1-AR subtype, which is situated postsynaptically in smooth muscles, that causes vasoconstriction of blood vessels when stimulated. Sympathetic overactivity in hypertension results in an excess stimulation of postsynaptic alpha1 ARs [25]. DA, a major neurotransmitter in both CNS and PNS holds an essential function in the homeostatic control of BP, by regulation of vascular smooth muscle contraction, epithelial sodium transport, and reactive oxygen species (ROS) production. It is considered a major player in homeostatic regulation of extra-cellular fluid (ECF) volume and BP due to the vast effect it induces on renal hemodynamics as well as humoral agents such as catecholamines, aldosterone, renin, vasopressin and endothelin (e.g., DA inhibits NA release and acts as a vasodilator at normal concentrations by activation of D2 receptors; decreases aldosterone secretion via activation of D3 receptors, while significantly enhancing renin secretion by activation of D1 receptors). DA also adjusts sodium and fluid intake by means of activities within the gastroin- testinal tract and CNS, primarily by regulation of the cardiovascular control centers in Biomedicines 2021, 9, 469 12 of 18

the brain stem [26]. Each DA receptor subtype participates in BP regulation by specific mechanisms for the subtype; DRD2 and DRD5 act within the CNS and PNS. The mesolimbic DA system is implicated in diurnal profiles of the mean BP. A clear dip in the mean BP and heart rate occurs during the resting period (e.g., nocturnal dip in human) [27]. An increase in arterial BP can be seen during the transition from non- REM to REM sleep, displaying phasic surges throughout REM sleep that derive from the physiological phase of paradoxical sleep. Furthermore, the mesolimbic DA system is involved in the increases in REM-associated BP fluctuations [28], a homeostatic adaptation to the fluctuating states of the different sleep and arousal phases. Since the mesolimbic DA system surges are involved in REM-associated increases in BP, the homeostatic effect is the result of a greater ability to deal with emotional stress. The pineal gland translates light signals received by the retina to the rest of the body, for example through the synthesis of the hormone melatonin, which is produced and released at night and helps to regulate the body’s metabolic activity during sleep. NA is involved in regulating this synthesis and release of melatonin in the pineal gland. DA through ADRA1B-DRD4 and ADRB1-DRD4 receptor heteromers inhibits the effects of NE, resulting in a decrease in the production and release of melatonin [29]. Although there are two classes of DA receptors, DR1-like and DR2-like, the natriuretic effect of DA is primarily mediated by DR1-like receptors [30]. During sodium loading, DR2-like receptors may contribute to the natriuresis [31]. Of the three DR2-like receptors it is likely that it is the DRD3 receptor that interacts with the DRD1 receptor because it is the major DR2-like receptor expressed in the renal proximal tubule and the thick ascending limb of Henle [32]. The DRD1, DRD3, DRD4 receptors interact with the renin-angiotensin system (RAAS), affecting epithelial transport and control of secretion of multiple humoral agents and their receptors [33]. From the comparative perspective, our results show a high degree of conservation (>60%) for catecholamines in vertebrate species analyzed, especially in mammals, as well as in invertebrates such as C. elegans, Drosophila melanogaster, and Ciona intestinalis (~45%). The vertebrate DR1-like receptors are most closely related to the DOP1 group members of C. elegans and Drosophila melanogaster, both responsible for upregulation of intracellular cAMP in the presence of DA. The vertebrate DR2-like receptors share the most homology with the C. elegans CeDOP2 and the Drosophila melanogaster DD2R, sharing a similar response in the decrease of intracellular cAMP levels when treated with DA [34]. In Drosophila melanogaster, expression patterns of DOP1 (an ortholog of the vertebrate D1-type DA receptors) were observed in mushroom body neurons, subesophageal ganglion, and in unpaired abdominal ganglia [35], while DD2R is expressed in the larval and adult nervous systems, and in cell groups that include peptidergic neurons [36]. In C. elegans it was shown that CeDOP1 is expressed in mechanosensory neurons, motor neurons and the ventral nerve cord, as well as sensory support cells, which are glial-like cells that surround the sensilla (simple invertebrate sense organs that may take the form of a hair or bristle) [37], while CeDOP2 is expressed in DA-nergic neurons, suggesting it may act as an autoreceptor [38].

4.5. The Kidney’s Role in Maintaining Homeostasis Notably, genes in clusters of rows I–III are expressed in the kidney of zebraﬁsh, mouse, and human. Their functional analysis demonstrated that homeostasis and transport are the primary biological processes that characterize these clusters, while signaling and stimulus are secondary. The cells of the kidney contain many specialized ion channels and transporters, which act in concert to regulate volume and ionic concentration by absorption or secretion of ions into the urine [39]. The ion channels and transporters play essential roles in organelle to whole organism function. These roles range from regulation of cell volume, membrane excitability and pH, to control of systemic salt and water balance and behavior [40]. Our work speciﬁcally relates to several prominent gene families—solute carriers (primarily electroneutral potassium chloride cotransporters, glucose transporters, bicarbonate Biomedicines 2021, 9, 469 13 of 18

transporter proteins), calcium voltage-gated channels, sodium channels, epithelial chloride voltage-gated channels, potassium voltage-gated channels, DRs, and ARs, which are conserved across most species examined (similarity > 50%). It is especially important to underscore the largest gene family we analyzed in this work (21 members)—Solute Carriers (SLCs), which are responsible for the regulation of various types of substances over the cell membrane. They typically rely on an ion gradient over the cell membrane as a mechanism for transportation [41]. Interestingly, all SLCs that are found in human are also found in C. elegans and Drosophila melanogaster and even some in T. adhaerens (SLC8A2, SLC12A3, SLC12A4, SLC4A4, SLC6A2 >50% similarity), indicating that this superfamily is ancient and was present before the divergence of Bilateria (animals with bilateral symmetry). Another noteworthy gene family we analyzed is the voltage-gated ion channel superfamily, comprising calcium, potassium, and sodium ion channels. As can be seen from the heatmap in Figure1, KCNJ1, KCNJ11, KCNJ6 are well conserved (>45%) in all 19 species, including T. adhaerens, along with CACNA1C, CACNA1H, CACNA1A, CACNA2D1 (>40%) in all 19 species. Yu et al. suggested that this family has evolved from a bacterial ancestor channel and accumulated regulatory domains for ligand binding in the course of evolution. The resulting signaling mechanisms control most aspects of cell physiology, complex processes in the brain and movements of muscles, all of which allowed the development of complex multicellular organisms [42]. Multicellular organisms control their internal environment by altering either the electrolytes concentration (osmolality) or the extracellular volume, maintaining the acid- base balance. In mammals, the main site of this complex regulation is the kidney. The role of human kidney is maintaining the fragile equilibrium within the internal milieu, namely the ECF and the intra-cellular fluid (ICF). From the evolutionary point of view, the human kidney evolved from primitive mechanisms of the ionic and osmotic homeostasis in fish. The data presented in this work provide evidence that the genes that are implicated in these mechanisms appear in fish, but not in earlier stages of phylogeny. The ion composition of the ECF is essentially identical in all animal species (including fish, amphibians, reptiles, and mammals) where the kidneys are utilized to regulate homeostasis [15]. Since salt concentration of the ocean increased over time, saltwater fish began to contain hypo-osmolar internal milieu and constantly exhausted water to the hyperosmolar environment. To endure this challenge, saltwater fish consume considerable amounts of water [43]; hence the acquired capability of fish to actively preserve the internal milieu is remarkable. The zebrafish pronephros, which consists of two nephrons with glomeruli, contains two tubules that are analogous in many ways to the segments of the mammalian nephron [44]. The proximal straight segment of this nephric tubule displays a brush border- like columnar epithelial cells [45]. This segment plays a crucial role in reabsorption of electrolytes and small molecules through a glomerular filtration barrier [46]. This is the physiological principle that underlies the mammalian renal function. Kidney functions are achieved both independently and in conjunction with other organs, particularly the endocrine system via various endocrine hormones, including, among others: renin (REN), angiotensin II, aldosterone, antidiuretic hormone (ADH), and two atrial natriuretic peptides (NPPA, NPPB) which are evolutionarily conserved. These natriuretic peptide hormones are secreted by atrial myocytes in response to stretch, angiotensin II stimulation, endothelin, and sympathetic (beta adrenergic) stimulation [47]. Blocking ion channels is an important part of anti-hypertensive pharmacotherapy. This is exemplified by the commonly used drugs—nifedipine, amlodipine, diltiazem, and verapamil which block the voltage-gated calcium channels of vascular smooth muscle cells, thus lowering BP; and amiloride that blocks the epithelial sodium channels at the distal convoluted tubule and the collecting duct, thereby inhibiting sodium-potassium exchange, while lowering BP independent of aldosterone [48,49]. Biomedicines 2021, 9, 469 14 of 18

4.6. Sensing the Environment in Order to Maintain Homeostatic Balance An extensive network of sensory organs and organelles play an important role in maintaining homeostatic balance and monitoring the internal environment. These include: mechanoreceptors sensitive to BP, chemoreceptors that react to the partial oxygen pressure in the arterial blood and chemoreceptors within the CNS, which are sensitive to pH and to numerous hormones. Ion channels in the hypothalamus gauge the ECF osmolality and sodium concentration. Baroreceptors situated in numerous circulatory beds along with ion channels and various receptors within the distal tubule in the kidney nephron, regulate blood volume. The input from the sense organs and organelles is transmitted for processing to the CNS, from which appropriate outputs are forwarded to the effector tissues—muscles and glands. The volume effectors are the renin–angiotensin–aldosterone system (RAAS), alongside sympathetic nerve activation (SNA) as well as glomerular filtration rate and physical gradients along the nephron. Homeostatic control mechanisms maintain the balance between fluid gain and fluid loss. Body water homeostasis is mainly regulated through fluid ingestion. Sensory osmoreceptors, found primarily in the hypothalamus, detect changes in plasma osmolarity and contribute to fluid-balance regulation in the body. If the body is fluid deficient, increased plasma osmolarity is sensed by the osmoreceptors. When osmoreceptors detect high plasma osmolarity (a sign of low blood volume/dehydration), they send signals to the hypothalamus, which creates the sensation of thirst. This result also triggers an increase in ADH secretion, causing fluid retention by the kidneys and urine output to be reduced. Macula densa cells in the nephron’s ascending loop of Henle are another type of osmoreceptor. If the macula densa is stimulated by high osmolarity, the juxtaglomerular apparatus (JGA) releases renin into the bloodstream. The subsequent production of angiotensin II acts on the hypothalamus to cause thirst. Angiotensin II causes vasoconstriction and aldosterone release. Aldosterone increases expression of the nephron’s ATPase pumps resulting in increased water reabsorption through sodium cotransport [50].

4.7. RAAS Evolution The RAAS is the primary volume-regulating pathway in mammals, acting as the central controller of blood pressure homeostasis in humans. The renin-angiotensin-aldosterone system emerged approximately 450 million years ago, in the Paleozoic era, as marine organisms moved to the land and endured a strong selection pressure to conserve salt and maintain volume homeostasis, which were crucial for survival [51]. In this work we have examined the sequence homologs for 20 human genes associated with RAAS pathway in 19 model organisms. Interestingly, homologs for all genes were present in all species from Danio rerio (Zebrafish) up to Pan troglodytes (Chimpanzee), with lowest similarity score of 46% for angiotensinogen (AGT). Since our data do not provide evidence for the presence of AGT homologs prior to Danio rerio, and due to its central role in the RAAS, this could indicate that the pathway emerged approximately with the appearance of zebrafish. This observation, in conjunction with the fact that zebrafish is the taxon with the most primitive juxtaglomerular apparatus [52], agrees with the physiological evidence regarding its emergence in bony fishes [53]. Our work reveals several orthologs for ACE and ACE2 from B. floridae to T. adhaerens, representing primitive chordates, as well as insects and multicellular organisms, which lack most RAAS components and are characterized by an open circulatory system. No ortholog sequence was found in the C. elegans, a fact that challenges the advent and initial roles of the enzyme. Recent genomic data revealed the presence of orthologs in yet additional remote phyla, as placozoa and proteobacteria, suggesting that it emerged early in evolution. Thus, by and large, ACE exists from bacteria up to mammals, exhibiting conserved features such as structure, as well as molecular and biochemical characteristics. In higher order organisms, ACE exerts its effect mainly in the capillaries of the lungs, whereas the evolution of the respiratory system was a key step in the transition of the animal world from marine life to land habitation. Characteristics of Biomedicines 2021, 9, 469 15 of 18

mammalian ACE could therefore be the outcome of an extended course of evolutionary specialization of an ancient protease, the functions of which are not yet known [54]. Adrenal steroids, both MCs and GCs, influence almost all aspects of vertebrate development, as well as regulation of an extensive selection of physiological processes, including development, differentiation, reproduction, and homeostasis [24]. Their receptors belong to the nuclear receptor super-family of proteins, which are ligand-activated transcription factors, implicated in steroid activity in a myriad of physiological, morphological, and behavioral processes [55]. GRs and MRs encoded in human by the NR3C1 and NR3C2 genes respectively, are representatives of the two principal functional families of corticosteroid receptors in vertebrates. These have not been found in invertebrates, as confirmed in our work, where orthologs for NR3C1 and NR3C2 appear first only in Danio rerio with similarity >63%. The MR and the GR are known to have descended from the early corticosteroid receptor [56], which first appeared in jawless fish, at the base of the vertebrate lineage, while distinct MR first appeared in cartilaginous fish [57]. Interestingly, the ancient GR evolved through duplication, to develop into a GR due to two gene mutations [58]. It is hypothesized that this evolutionary process would have alleviated the difficulty to adapt to blood pressure elevation caused by increased gravitational effects of land habitation [59].

5. Conclusions To conclude this study, it appears that the genes that are involved in BP regulation are primarily implicated in mechanisms associated with the maintenance of homeostasis. Evidence for this can be seen in the conserved genes across the entire range of organisms investigated, performing functions that are notably related to homeostasis, e.g., ion transport, response to stimulus and more, as well as expression patterns in tissues and organs that support homeostatic processes in the respective species. This work is merely an attempt to cover the broad scope of the evolutionary dissection of the genetic machinery behind BP regulation; nonetheless it highlights the importance of the comparative approach that presupposes that the comprehension of our evolutionary past can highly enlighten our understanding of the genetic background of common disease.

Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/biomedicines9050469/s1, Table S1: GO term biological processes enrichment per cluster genes and organisms, Table S2: tissue/organ enrichment per cluster genes and organisms. Author Contributions: Conceptualization, A.B. and R.U.; methodology, A.B. and R.U.; formal analysis, A.B.; investigation, A.B.; resources, R.U.; data curation, A.B.; writing—original draft preparation, A.B.; writing—review and editing, Y.F. and R.U.; visualization, A.B.; supervision, Y.F. and R.U.; project administration, A.B. All authors have read and agreed to the published version of the manuscript. Funding: This research received no external funding. Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: Data supporting reported results of this study can be found in publicly archived datasets as specified in links throughout the manuscript. Acknowledgments: Special thanks to John Moult for constructive discussions and a critical perspective. The authors are also indebted to Gil Kotton, Tirza Doniger and Orit Adato for their valuable assistance throughout this work. Conflicts of Interest: The authors declare no conflict of interest. Biomedicines 2021, 9, 469 16 of 18

Abbreviations

CNS Central Nervous System ANS Autonomous Nervous System RAAS Renin Angiotensin Aldosterone System DA Dopamine NA Noradrenaline CTZ Chemoreceptor Trigger Zone ACTH Adrenocorticotropic Hormone GC Glucocorticoid GR Glucocorticoid Receptor AR Adrenergic Receptor ECF Extracellular Fluid SLC Solute Carrier GO Gene Ontology

References 1. Rossier, B.C.; Bochud, M.; Devuyst, O. The hypertension pandemic: An evolutionary perspective. Physiology 2017, 32, 112–125. [CrossRef][PubMed] 2. Monahan-Earley, R.; Dvorak, A.M.; Aird, W.C. Evolutionary origins of the blood vascular system and endothelium. J. Thromb. Haemost. 2013, 11, 46–66. [CrossRef][PubMed] 3. Bernard, C. Lectures on the Phenomena of Life Common to Animals and Plants; Charles C Thomas Publisher: Springfield, IL, USA, 1974; ISBN 0398028575. 4. Botzer, A.; Grossman, E.; Moult, J.; Unger, R. A system view and analysis of essential hypertension. J. Hypertens. 2018, 36, 1094–1103. [CrossRef][PubMed] 5. Szklarczyk, D.; Morris, J.H.; Cook, H.; Kuhn, M.; Wyder, S.; Simonovic, M.; Santos, A.; Doncheva, N.T.; Roth, A.; Bork, P.; et al. The STRING database in 2017: Quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017, 45, D362–D368. [CrossRef] 6. Huerta-Cepas, J.; Szklarczyk, D.; Forslund, K.; Cook, H.; Heller, D.; Walter, M.C.; Rattei, T.; Mende, D.R.; Sunagawa, S.; Kuhn, M.; et al. eggNOG 4.5: A hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences. Nucleic Acids Res. 2016, 44.[CrossRef][PubMed] 7. Sonnhammer, E.L.L.; Ostlund, G. InParanoid 8: Orthology analysis between 273 proteomes, mostly eukaryotic. Nucleic Acids Res. 2015, 43.[CrossRef] 8. Madeira, F.; Park, Y.M.; Lee, J.; Buso, N.; Gur, T.; Madhusoodanan, N.; Basutkar, P.; Tivey, A.R.N.; Potter, S.C.; Finn, R.D.; et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 2019, 47, W636–W641. [CrossRef] 9. Greenacre, M.J.; Primicerio, R. Multivariate Analysis of Ecological Data; Fundacion BBVA: Bilbao, Spain, 2013; ISBN 9788492937509 8492937505. 10. Lieschke, G.J.; Currie, P.D. Animal models of human disease: Zebrafish swim into view. Nat. Rev. Genet. 2007, 8, 353–367. [CrossRef] 11. Guyton, A.C.; Hall, J.E. Textbook of Medical Physiology, 9th ed.; W.B. Saunders Company: Philadelphia, PA, USA, 1996. 12. Walch, L.; Brink, C.; Norel, X. The muscarinic receptor subtypes in human blood vessels. Therapie 2001, 56, 223–226. 13. Whelton, P.K.; Carey, R.M.; Aronow, W.S.; Casey, D.E.; Collins, K.J.; Dennison Himmelfarb, C.; DePalma, S.M.; Gidding, S.; Jamerson, K.A.; Jones, D.W.; et al. 2017 ACC/AHA/AAPA/ABC/ACPM/AGS/APhA/ASH/ASPC/NMA/PCNA Guideline for the Prevention, Detection, Evaluation, and Management of High Blood Pressure in Adults: A Report of the American College of Cardiology/American Heart Association Task Force on Clinical Pr. Hypertension 2018, 71.[CrossRef] 14. Schulte, K.; Kunter, U.; Moeller, M.J. The evolution of blood pressure and the rise of mankind. Nephrol. Dial. Transplant. 2015, 30, 713–723. [CrossRef][PubMed] 15. Epstein, F.H. The sea within us. J. Exp. Zool. 1999, 284, 50–54. [CrossRef] 16. Nesse, R.M.; Bhatnagar, S.; Young, E.A. Evolutionary Origins and Functions of the Stress Response. In Encyclopedia of Stress; Academic Press: Cambridge, MA, USA, 2007; pp. 965–970. ISBN 9780123739476. 17. Finkelstein, Y.; Koffler, B.; Rabey, J.M.; Gilad, G.M. Dynamics of cholinergic synaptic mechanisms in rat hippocampus after stress. Brain Res. 1985, 343, 314–319. [CrossRef] 18. Dantzer, R. Neuroimmune interactions: From the brain to the immune system and vice versa. Physiol. Rev. 2018, 98, 477–504. [CrossRef][PubMed] 19. McEwen, B.S.; Nasca, C.; Gray, J.D. Stress effects on neuronal structure: Hippocampus, amygdala, and prefrontal cortex. Neuropsychopharmacology 2016, 41, 3–23. [CrossRef] 20. Blaess, S.; Stott, S.R.W.; Ang, S.-L. The generation of midbrain dopaminergic neurons. In Patterning and Cell Type Specification in the Developing CNS and PNS; Academic Press: Cambridge, MA, USA, 2020; pp. 369–398. Biomedicines 2021, 9, 469 17 of 18

21. Preedy, V.R. Neuropathology of Drug Addictions and Substance Misuse; Academic Press: Cambridge, MA, USA, 2016; Volume 2, ISBN 9780128003756. 22. Lumb, R.; Schwarz, Q. Sympathoadrenal neural crest cells: The known, unknown and forgotten? Dev. Growth Differ. 2015, 57, 146–157. [CrossRef] 23. Grassi Milano, E.; Basari, F.; Chimenti, C. Adrenocortical and adrenomedullary homologs in eight species of adult and developing teleosts: Morphology, histology, and immunohistochemistry. Gen. Comp. Endocrinol. 1997, 108, 483–496. [CrossRef] 24. Baker, M.E. Evolution of adrenal and sex steroid action in vertebrates: A ligand-based mechanism for complexity. BioEssays 2003, 25, 396–400. [CrossRef] 25. Michel, M.C.; Brodde, O.E.; Insel, P.A. Peripheral adrenergic receptors in hypertension. Hypertension 1990, 16, 107–120. [CrossRef] 26. Mannelli, M.; Pupilli, C.; Lanzillotti, R.; Ianni, L.; Serio, M. Catecholamines and blood pressure regulation. Horm. Res. 1990, 34, 156–160. [CrossRef] 27. Sakata, M.; Sei, H.; Toida, K.; Fujihara, H.; Urushihara, R.; Morita, Y. Mesolimbic dopaminergic system is involved in diurnal blood pressure regulation. Brain Res. 2002, 951, 150. [CrossRef] 28. Sei, H.; Morita, Y. Why does arterial blood pressure rise actively during REM sleep? J. Med. Investig. 1999, 46, 11–17. 29. González, S.; Moreno-Delgado, D.; Moreno, E.; Pérez-Capote, K.; Franco, R.; Mallol, J.; Cortés, A.; Casadó, V.; Lluís, C.; Ortiz, J.; et al. Circadian-related heteromerization of adrenergic and dopamine d4 receptors modulates melatonin synthesis and release in the pineal gland. PLoS Biol. 2012, 10.[CrossRef] 30. Hussain, T.; Lokhandwala, M.F. Renal dopamine receptor function in hypertension. Hypertension 1998, 32, 187–197. [CrossRef] [PubMed] 31. Jose, P.A.; Asico, L.D.; Eisner, G.M.; Pocchiari, F.; Semeraro, C.; Felder, R.A. Effects of costimulation of dopamine D1- and D2-like receptors on renal function. Am. J. Physiol. 1998, 275, R986–R994. [CrossRef][PubMed] 32. O’Connell, D.P.; Vaughan, C.J.; Aherne, A.M.; Botkin, S.J.; Wang, Z.-Q.; Felder, R.A.; Carey, R.M. Expression of the dopamine D3 receptor protein in the rat kidney. Hypertension 1998, 32, 886–895. [CrossRef][PubMed] 33. Jose, P.A.; Eisner, G.M.; Felder, R.A. Dopamine and the kidney: A role in hypertension? Curr. Opin. Nephrol. Hypertens. 2003, 12, 189–194. [CrossRef] 34. Mustard, J.A.; Beggs, K.T.; Mercer, A.R. Molecular biology of the invertebrate dopamine receptors. Arch. Insect Biochem. Physiol. 2005, 59, 103–117. [CrossRef] 35. Kim, Y.C.; Lee, H.G.; Seong, C.S.; Han, K.A. Expression of a D1 dopamine receptor dDA1/DmDOP1 in the central nervous system of Drosophila melanogaster. Gene Expr. Patterns 2003, 3, 237–245. [CrossRef] 36. Draper, I.; Kurshan, P.T.; McBride, E.; Jackson, F.R.; Kopin, A.S. Locomotor activity is regulated by D2-like receptors in Drosophila: An anatomic and functional analysis. Dev. Neurobiol. 2007, 67, 378–393. [CrossRef][PubMed] 37. Tsalik, E.L.; Niacaris, T.; Wenick, A.S.; Pau, K.; Avery, L.; Hobert, O. LIM homeobox gene-dependent expression of biogenic amine receptors in restricted regions of the C. elegans nervous system. Dev. Biol. 2003, 263, 81–102. [CrossRef] 38. Suo, S.; Sasagawa, N.; Ishiura, S. Identiﬁcation of a dopamine receptor from Caenorhabditis elegans. Neurosci. Lett. 2002, 319, 13–16. [CrossRef] 39. Kuo, I.Y.; Ehrlich, B.E. Ion channels in renal disease. Chem. Rev. 2012, 112, 6353–6372. [CrossRef][PubMed] 40. Strange, K. From genes to integrative physiology: Ion channel and transporter biology in Caenorhabditis elegans. Physiol. Rev. 2003, 83, 377–415. [CrossRef][PubMed] 41. Hoglund, P.J.; Nordstrom, K.J.V.; Schioth, H.B.; Fredriksson, R. The solute carrier families have a remarkably long evolutionary history with the majority of the human families present before divergence of bilaterian species. Mol. Biol. Evol. 2011, 28, 1531–1541. [CrossRef][PubMed] 42. Yu, F.H.; Catterall, W.A. The VGL-chanome: A protein superfamily specialized for electrical signaling and ionic homeostasis. Sci. STKE 2004, 2004.[CrossRef] 43. Hirano, T. Some factors regulating water intake by the eel, Anguilla japonica. J. Exp. Biol. 1974, 61, 737–747. 44. Wingert, R.A.; Davidson, A.J. The zebraﬁsh pronephros: A model to study nephron segmentation. Kidney Int. 2008, 73, 1120–1127. [CrossRef] 45. Alpern, R.; Moe, O.; Caplan, M. Seldin and Geibisch’s The Kidney; Academic Press: Cambridge, MA, USA, 2013; Volume 1, ISBN 9780123814623. 46. Carroll, T.J.; Wallingford, J.B. Induction, Development, and Physiology of the Pronephric Tubules. Kidney 2003, 19–50. [CrossRef] 47. Man, J.; Barnett, P.; Christoffels, V.M. Structure and function of the Nppa–Nppb cluster locus during heart development and disease. Cell. Mol. Life Sci. 2018, 75, 1435–1444. [CrossRef][PubMed] 48. Baker, E.H. Ion channels and the control of blood pressure. Br. J. Clin. Pharmacol. 2000, 49, 185–198. [CrossRef][PubMed] 49. Botzer, A.; Finkelstein, Y.; Grossman, E.; Moult, J.; Unger, R. Iatrogenic hypertension: A bioinformatic analysis. Pharmacogenom. J. 2018, 19, 1–34. [CrossRef][PubMed] 50. Water Balance—Boundless Anatomy and Physiology. Available online: https://courses.lumenlearning.com/boundless-ap/ chapter/water-balance/ (accessed on 20 April 2021). 51. Mansley, M.K.; Ivy, J.R.; Bailey, M.A. ISN Forefronts Symposium 2015: The evolution of hypertension—Old genes, new concepts. Kidney Int. Rep. 2016, 1, 197–203. [CrossRef][PubMed] 52. Sokabe, H.; Ogawa, M. Comparative studies of the Juxtaglomerular apparatus. Int. Rev. Cytol. 1974, 37, 271–327. [CrossRef] Biomedicines 2021, 9, 469 18 of 18

53. Nishimura, H.; Bailey, J.R. Intrarenal renin-angiotensin system in primitive vertebrates. Kidney Int. 1982, 12, S185–S192. 54. Rivière, G.; Michaud, A.; Deloffre, L.; Vandenbulcke, F.; Levoye, A.; Breton, C.; Corvol, P.; Salzet, M.; Vieau, D. Characterization of the ﬁrst non-insect invertebrate functional angiotensin-converting enzyme (ACE): Leech TtACE resembles the N-domain of mammalian ACE. Biochem. J. 2004, 382, 565–573. [CrossRef] 55. Bentley, P.J. Comparative Vertebrate Endocrinology; Cambridge University Press: Cambridge, MA, USA, 1998; ISBN 0521629985. 56. Thornton, J.W. Evolution of vertebrate steroid receptors from an ancestral estrogen receptor by ligand exploitation and serial genome expansions. Proc. Natl. Acad. Sci. USA 2001, 98, 5671–5676. [CrossRef] 57. Osório, J.; Rétaux, S. The lamprey in evolutionary studies. Dev. Genes Evol. 2008, 218, 221–235. [CrossRef] 58. Bridgham, J.T.; Carroll, S.M.; Thornton, J.W. Evolution of Hormone-Receptor Complexity by Molecular Exploitation. Science 2006, 312, 97–101. [CrossRef] 59. Kvetnansky, R.; Lu, X.; Ziegler, M.G. Stress-triggered changes in peripheral catecholaminergic systems. Adv. Pharmacol. 2013, 68, 359–397. [CrossRef][PubMed]