Diversity in the Glucose Transporter-4 Gene (SLC2A4) in Humans Reflects the Action of Natural Selection along the Old-World Primates Evolution

Background

Glucose is an important source of energy for living organisms. In vertebrates it is ingested with the diet and transported into the cells by conserved mechanisms and molecules, such as the trans-membrane Glucose Transporters (GLUTs). Members of this family have tissue specific expression, biochemical properties and physiologic functions that together regulate glucose levels and distribution. GLUT4 –coded by SLC2A4 (17p13) is an insulin-sensitive transporter with a critical role in glucose homeostasis and diabetes pathogenesis, preferentially expressed in the adipose tissue, heart muscle and skeletal muscle. We tested the hypothesis that natural selection acted on SLC2A4.

Methodology/Principal Findings

We re-sequenced SLC2A4 and genotyped 104 SNPs along a ∼1 Mb region flanking this gene in 102 ethnically diverse individuals. Across the studied populations (African, European, Asian and Latin-American), all the eight common SNPs are concentrated in the N-terminal region upstream of exon 7 (∼3700 bp), while the C-terminal region downstream of intron 6 (∼2600 bp) harbors only 6 singletons, a pattern that is not compatible with neutrality for this part of the gene. Tests of neutrality based on comparative genomics suggest that: (1) episodes of natural selection (likely a selective sweep) predating the coalescent of human lineages, within the last 25 million years, account for the observed reduced diversity downstream of intron 6 and, (2) the target of natural selection may not be in the SLC2A4 coding sequence.

Conclusions

We propose that the contrast in the pattern of genetic variation between the N-terminal and C-terminal regions are signatures of the action of natural selection and thus follow-up studies should investigate the functional importance of differnet regions of the SLC2A4 gene.     

Phylogenesis and Biological Characterization of a New Glucose Transporter in the Chicken (Gallus gallus), GLUT12

In mammals, insulin-sensitive GLUTs, including GLUT4, are recruited to the plasma membrane of adipose and muscle tissues in response to insulin. The GLUT4 gene is absent from the chicken genome, and no functional insulin-sensitive GLUTs have been characterized in chicken tissues to date. A nucleotide sequence is predicted to encode a chicken GLUT12 ortholog and, interestingly, GLUT12 has been described to act as an insulin-sensitive GLUT in mammals. It encodes a 596 amino acid protein exhibiting 71% identity with human GLUT12. First, we present the results of a phylogenetic study showing the stability of this gene during evolution of vertebrates. Second, tissue distribution of chicken SLC2A12 mRNA was characterized by RT-PCR. It was predominantly expressed in skeletal muscle and heart. Protein distribution was analysed by Western blotting using an anti-human GLUT12 antibody directed against a highly conserved region (87% of identity). An immuno-reactive band of the expected size (75kDa) was detected in the same tissues. Third a physiological characterization was performed: SLC2A12 mRNA levels were significantly lowered in fed chickens subjected to insulin immuno-neutralization. Finally, recruitment of immuno-reactive GLUT12 to the muscle plasma membrane was increased following 1h of intraperitoneal insulin administration (compared to a control fasted state). Thus insulin administration elicited membrane GLUT12 recruitment. In conclusion, these results suggest that the facilitative glucose transporter protein GLUT12 could act in chicken muscle as an insulin-sensitive transporter that is qualitatively similar to GLUT4 in mammals.     

Rabbit SLC15A1, SLC7A1 and SLC1A1 genes are affected by site of digestion,stage of development and dietary protein content

Peptide transporter 1 (SLC15A1, PepT1), excitatory amino acid transporter 3 (SLC1A1, EAAT3) and cationic amino acid transporter 1 (SLC7A1, CAT1) were identified as genes responsible for the transport of small peptides and amino acids. The tissue expression pattern of rabbit (SLC15A1, SLC7A1 and SLC1A1) across the digestive tract remains unclear. The present study investigated SLC15A1, SLC7A1 and SLC1A1 gene expression patterns across the digestive tract at different stages of development and in response to dietary protein levels. Real time-PCR results indicated that SLC15A1, SLC7A1 and SLC1A1 genes throughout the rabbits’ entire development and were expressed in all tested rabbit digestive sites, including the stomach, duodenum, jejunum, ileum, colon and cecum. Furthermore, SLC7A1 and SLC1A1 mRNA expression occurred in a tissue-specific and time-associated manner, suggesting the distinct transport ability of amino acids in different tissues and at different developmental stages. The most highly expressed levels of all three genes were in the duodenum, ileum and jejunum in all developmental stages. All increased after lactation. With increased dietary protein levels, SLC7A1 mRNA levels in small intestine and SLC1A1 mRNA levels in duodenum and ileum exhibited a significant decreasing trend. Moreover, rabbits fed a normal level of protein had the highest levels of SLC15A1 mRNA in the duodenum and jejunum (P<0.05). In conclusion, gene mRNA differed across sites and with development suggesting time and sites related differences in peptide and amino acid absorption in rabbits. The effects of dietary protein on expression of the three genes were also site specific.     

A novel compound heterozygous variant Of SLC12A3 Gene in A Pedigree with Gitelman Syndrome Co-Existent with Thyroid Dysfunction

Objective: Gitelman syndrome (GS) is an autosomal recessive disorder characterized by salt wasting and hypokalemia resulting from mutations in the SLC12A3 (solute carrier family 12 member 3) gene, which encodes the thiazide-sensitive sodium-chloride cotransporter. To date, more than 488 mutations of the SLC12A3 gene have been discovered in patients with GS. In this study, we reported a GS pedigree complicated by thyroid diseases or thyroid dysfunction.Methods: Sanger sequencing and next-generation sequencing analysis were performed to determine the SLC12A3 gene mutations in a GS pedigree including the 16-year old male patient with GS and his family members within 3 generations. Chemiluminescence immunoassays were used to detect thyroid hormone and antibody concentrations.Results: Genetic analysis of the SLC12A3 gene identified 2 mutations in the 16-year old male patient with GS concomitant with Graves disease (GD) and his younger sister accompanied by abnormal thyroid function. Additionally, one mutation site (c.1456G>A) in SLC12A3 gene was found in his father, paternal uncle and elder female cousin, who were complicated by subclinical hypothyroidism or autoantibody against thyroid. The other mutation site (c.2102_2107 delACAAGA) in SLC12A3 gene, a novel mutated variant of SLC12A3 gene, was carried by his mother and maternal grandfather.Conclusion: Two mutation sites were documented in the pedigree with GS, and one has not been reported before. Moreover, we found a mutation at nucleotide c.1456 G>A in the SLC12A3 gene that may affect thyroid function. However, further studies are needed to explore the underlying molecular mechanisms.Abbreviations: FT3 = free triiodothyronine; FT4 = free tetraiodothyronine; GD = Graves disease; GS = Gitelman syndrome; SLC12A3 = solute carrier family 12 member 3; TGAb = thyroglobulin antibody; TPOAb = thyroid peroxidase antibody; TSH = thyroid-stimulating hormone; TT3 = total triiodothyronine; TT4 = total tetraiodothyronine     

Genome-Wide Association Analysis Identifies a Mutation in the Thiamine Transporter 2 (SLC19A3) Gene Associated with Alaskan Husky Encephalopathy

Alaskan Husky Encephalopathy (AHE) has been previously proposed as a mitochondrial encephalopathy based on neuropathological similarities with human Leigh Syndrome (LS). We studied 11 Alaskan Husky dogs with AHE, but found no abnormalities in respiratory chain enzyme activities in muscle and liver, or mutations in mitochondrial or nuclear genes that cause LS in people. A genome wide association study was performed using eight of the affected dogs and 20 related but unaffected control AHs using the Illumina canine HD array. SLC19A3 was identified as a positional candidate gene. This gene controls the uptake of thiamine in the CNS via expression of the thiamine transporter protein THTR2. Dogs have two copies of this gene located within the candidate interval (SLC19A3.2 – 43.36–43.38 Mb and SLC19A3.1 – 43.411–43.419 Mb) on chromosome 25. Expression analysis in a normal dog revealed that one of the paralogs, SLC19A3.1, was expressed in the brain and spinal cord while the other was not. Subsequent exon sequencing of SLC19A3.1 revealed a 4bp insertion and SNP in the second exon that is predicted to result in a functional protein truncation of 279 amino acids (c.624 insTTGC, c.625 C>A). All dogs with AHE were homozygous for this mutation, 15/41 healthy AH control dogs were heterozygous carriers while 26/41 normal healthy AH dogs were wild type. Furthermore, this mutation was not detected in another 187 dogs of different breeds. These results suggest that this mutation in SLC19A3.1, encoding a thiamine transporter protein, plays a critical role in the pathogenesis of AHE.     

Mito-nuclear genetic comparison in a Wolbachia infected weevil: insights on reproductive mode, infection age and evolutionary forces shaping genetic variation

Background

The composition and expression of vertebrate gene families is shaped by species specific gene loss in combination with a number of gene and genome duplication events (R1, R2 in all vertebrates, R3 in teleosts) and depends on the ecological and evolutionary context. In this study we analyzed the evolutionary history of the solute carrier 1 (SLC1) gene family. These genes are supposed to be under strong selective pressure (purifying selection) due to their important role in the timely removal of glutamate at the synapse.

Results

In a genomic survey where we manually annotated and analyzing sequences from more than 300 SLC1 genes (from more than 40 vertebrate species), we found evidence for an interesting evolutionary history of this gene family. While human and mouse genomes contain 7 SLC1 genes, in prototheria, sauropsida, and amphibia genomes up to 9 and in actinopterygii up to 13 SLC1 genes are present. While some of the additional slc1 genes in ray-finned fishes originated from R3, the increased number of SLC1 genes in prototheria, sauropsida, and amphibia genomes originates from specific genes retained in these lineages. Phylogenetic comparison and microsynteny analyses of the SLC1 genes indicate, that theria genomes evidently lost several SLC1 genes still present in the other lineage. The genes lost in theria group into two new subfamilies of the slc1 gene family which we named slc1a8/eaat6 and slc1a9/eaat7.

Conclusions

The phylogeny of the SLC1/EAAT gene family demonstrates how multiple genome reorganization and duplication events can influence the number of active genes. Inactivation and preservation of specific SLC1 genes led to the complete loss of two subfamilies in extant theria, while other vertebrates have retained at least one member of two newly identified SLC1 subfamilies.     

SLC2A9 is a high-capacity urate transporter in humans

Background

Serum uric acid levels in humans are influenced by diet, cellular breakdown, and renal elimination, and correlate with blood pressure, metabolic syndrome, diabetes, gout, and cardiovascular disease. Recent genome-wide association scans have found common genetic variants of SLC2A9 to be associated with increased serum urate level and gout. The SLC2A9 gene encodes a facilitative glucose transporter, and it has two splice variants that are highly expressed in the proximal nephron, a key site for urate handling in the kidney. We investigated whether SLC2A9 is a functional urate transporter that contributes to the longstanding association between urate and blood pressure in man.

Methods and Findings

We expressed both SLC2A9 splice variants in Xenopus laevis oocytes and found both isoforms mediate rapid urate fluxes at concentration ranges similar to physiological serum levels (200–500 μM). Because SLC2A9 is a known facilitative glucose transporter, we also tested whether glucose or fructose influenced urate transport. We found that urate is transported by SLC2A9 at rates 45- to 60-fold faster than glucose, and demonstrated that SLC2A9-mediated urate transport is facilitated by glucose and, to a lesser extent, fructose. In addition, transport is inhibited by the uricosuric benzbromarone in a dose-dependent manner (K _i = 27 μM). Furthermore, we found urate uptake was at least 2-fold greater in human embryonic kidney (HEK) cells overexpressing SLC2A9 splice variants than nontransfected kidney cells. To confirm that our findings were due to SLC2A9, and not another urate transporter, we showed that urate transport was diminished by SLC2A9-targeted siRNA in a second mammalian cell line. In a cohort of men we showed that genetic variants of SLC2A9 are associated with reduced urinary urate clearance, which fits with common variation at SLC2A9 leading to increased serum urate. We found no evidence of association with hypertension (odds ratio 0.98, 95% confidence interval [CI] 0.9 to 1.05, p > 0.33) by meta-analysis of an SLC2A9 variant in six case–control studies including 11,897 participants. In a separate meta-analysis of four population studies including 11,629 participants we found no association of SLC2A9 with systolic (effect size −0.12 mm Hg, 95% CI −0.68 to 0.43, p = 0.664) or diastolic blood pressure (effect size −0.03 mm Hg, 95% CI −0.39 to 0.31, p = 0.82).

Conclusions

This study provides evidence that SLC2A9 splice variants act as high-capacity urate transporters and is one of the first functional characterisations of findings from genome-wide association scans. We did not find an association of the SLC2A9 gene with blood pressure in this study. Our findings suggest potential pathogenic mechanisms that could offer a new drug target for gout.     

Mutations in Glucose Transporter 9 Gene SLC2A9 Cause Renal Hypouricemia

Renal hypouricemia is an inherited disorder characterized by impaired renal urate (uric acid) reabsorption and subsequent low serum urate levels, with severe complications such as exercise-induced acute renal failure and nephrolithiasis. We previously identified SLC22A12, also known as URAT1, as a causative gene of renal hypouricemia. However, hypouricemic patients without URAT1 mutations, as well as genome-wide association studies between urate and SLC2A9 (also called GLUT9), imply that GLUT9 could be another causative gene of renal hypouricemia. With a large human database, we identified two loss-of-function heterozygous mutations in GLUT9, which occur in the highly conserved “sugar transport proteins signatures 1/2.” Both mutations result in loss of positive charges, one of which is reported to be an important membrane topology determinant. The oocyte expression study revealed that both GLUT9 isoforms showed high urate transport activities, whereas the mutated GLUT9 isoforms markedly reduced them. Our findings, together with previous reports on GLUT9 localization, suggest that these GLUT9 mutations cause renal hypouricemia by their decreased urate reabsorption on both sides of the renal proximal tubules. These findings also enable us to propose a physiological model of the renal urate reabsorption in which GLUT9 regulates serum urate levels in humans and can be a promising therapeutic target for gout and related cardiovascular diseases.     

Complex Analysis of Urate Transporters SLC2A9, SLC22A12 and Functional Characterization of Non-Synonymous Allelic Variants of GLUT9 in the Czech Population: No Evidence of Effect on Hyperuricemia and Gout

Objective

Using European descent Czech populations, we performed a study of SLC2A9 and SLC22A12 genes previously identified as being associated with serum uric acid concentrations and gout. This is the first study of the impact of non-synonymous allelic variants on the function of GLUT9 except for patients suffering from renal hypouricemia type 2.

Methods

The cohort consisted of 250 individuals (150 controls, 54 nonspecific hyperuricemics and 46 primary gout and/or hyperuricemia subjects). We analyzed 13 exons of SLC2A9 (GLUT9 variant 1 and GLUT9 variant 2) and 10 exons of SLC22A12 by PCR amplification and sequenced directly. Allelic variants were prepared and their urate uptake and subcellular localization were studied by Xenopus oocytes expression system. The functional studies were analyzed using the non-parametric Wilcoxon and Kruskall-Wallis tests; the association study used the Fisher exact test and linear regression approach.

Results

We identified a total of 52 sequence variants (12 unpublished). Eight non-synonymous allelic variants were found only in SLC2A9: rs6820230, rs2276961, rs144196049, rs112404957, rs73225891, rs16890979, rs3733591 and rs2280205. None of these variants showed any significant difference in the expression of GLUT9 and in urate transport. In the association study, eight variants showed a possible association with hyperuricemia. However, seven of these were in introns and the one exon located variant, rs7932775, did not show a statistically significant association with serum uric acid concentration.

Conclusion

Our results did not confirm any effect of SLC22A12 and SLC2A9 variants on serum uric acid concentration. Our complex approach using association analysis together with functional and immunohistochemical characterization of non-synonymous allelic variants did not show any influence on expression, subcellular localization and urate uptake of GLUT9.     

Mutations in the SLC3A1 Transporter Gene in Cystinuria

Cystinuria is an autosomal recessive disease characterized by the development of kidney stones. Guided by the identification of the SLC3A1 amino acid–transport gene on chromosome 2, we recently established genetic linkage of cystinuria to chromosome 2p in 17 families, without evidence for locus heterogeneity. Other authors have independently identified missense mutations in SLC3A1 in cystinuria patients. In this report we describe four additional cystinuria-associated mutations in this gene: a frameshift, a deletion, a transversion inducing a critical amino acid change, and a nonsense mutation. The latter stop codon was found in all of eight Ashkenazi Jewish carrier chromosomes examined. This report brings the number of disease-associated mutations in this gene to 10. We also assess the frequency of these mutations in our 17 cystinuria families.     

Evolutionary Genomics and Adaptive Evolution of the Hedgehog Gene Family (Shh,Ihh and Dhh) in Vertebrates

The Hedgehog (Hh) gene family codes for a class of secreted proteins composed of two active domains that act as signalling molecules during embryo development, namely for the development of the nervous and skeletal systems and the formation of the testis cord. While only one Hh gene is found typically in invertebrate genomes, most vertebrates species have three (Sonic hedgehog – Shh; Indian hedgehog – Ihh; and Desert hedgehog – Dhh), each with different expression patterns and functions, which likely helped promote the increasing complexity of vertebrates and their successful diversification. In this study, we used comparative genomic and adaptive evolutionary analyses to characterize the evolution of the Hh genes in vertebrates following the two major whole genome duplication (WGD) events. To overcome the lack of Hh-coding sequences on avian publicly available databases, we used an extensive dataset of 45 avian and three non-avian reptilian genomes to show that birds have all three Hh paralogs. We find suggestions that following the WGD events, vertebrate Hh paralogous genes evolved independently within similar linkage groups and under different evolutionary rates, especially within the catalytic domain. The structural regions around the ion-binding site were identified to be under positive selection in the signaling domain. These findings contrast with those observed in invertebrates, where different lineages that experienced gene duplication retained similar selective constraints in the Hh orthologs. Our results provide new insights on the evolutionary history of the Hh gene family, the functional roles of these paralogs in vertebrate species, and on the location of mutational hotspots.     

SLC1A2 Variant Is Associated with Essential Tremor in Taiwanese Population

Essential tremor (ET), which is one of the most common movement disorders, may lead to severe interference in quality of life. The first genome-wide association study (GWAS) has identified an association of the LINGO1 variant (rs9652490) with ET in Americans and Europeans. Recently, a second GWAS that was performed in a European population has discovered a new variant (rs3794087) of the main glial glutamate transporter (SLC1A2) that increases the risk of ET with an odds ratio of about 1.4. SLC1A2 encodes for the major glial high-affinity glutamate reuptake transporter in the brain and is a potential ET susceptibility gene. Because replication in a different ethnic population is important for validating a finding, we conducted a case-control study to investigate the SLC1A2 variant in an Asian cohort with ET in Taiwan. A total of 542 subjects (273 ET patients and 269 controls) were included. The results showed that rs3794087 was associated with ET among the Taiwanese. The odds ratio was 1.37. Our results were similar to those of the second GWAS of ET in Europeans, and this confirms that SLC1A2 may be a good functional candidate gene for ET. A replication study in another independent population is of importance to validate this association.     

Influence of Familial Renal Glycosuria Due to Mutations in the SLC5A2 Gene on Changes in Glucose Tolerance over Time

Familial renal glycosuria is an inherited disorder resulting in glucose excretion in the urine despite normal blood glucose concentrations. It is most commonly due to mutations in the SLC5A2 gene coding for the glucose transporter SGLT2 in the proximal tubule. Several drugs have been introduced as means to lower glucose in patients with type 2 diabetes targeting SGLT2 resulting in renal glycosuria, but no studies have addressed the potential effects of decreased renal glucose reabsorption and chronic glycosuria on the prevention of glucose intolerance. Here we present data on a large pedigree with renal glycosuria due to two mutations (c.300-303+2del and p.A343V) in the SLC5A2 gene. The mutations, which in vitro affected glucose transport in a cell line model, and the ensuing glycosuria were not associated with better glycemic control during a follow-up period of more than 10 years. One individual, who was compound heterozygous for mutations in the SLC5A2 gene suffered from severe urogenital candida infections and postprandial hypoglycemia. In conclusion, in this family with familial glycosuria we did not find any evidence that chronic loss of glucose in the urine would protect from deterioration of the glucose tolerance over time.     

High SLC7A11 expression in normal skin of melanoma patients

BackgroundMelanoma is one of the highest metastatic cancers and its incidence is rapidly increasing. A great effort has been devoted to determine gene mutations and expression profiles in melanoma cells, but less attention has been given to the possible influence of melanin synthesis in melanocytes on melanomagenesis. SLC7A11 encodes the cystine/glutamate antiporter xCT and its expression increases the antioxidant capacity of cells by providing cysteine that may be used for glutathione (GSH) synthesis. Melanocytes, however, can also use cysteine for pheomelanin synthesis and pigmentation. Therefore, pheomelanin synthesis may lead to chronic oxidative stress. Possible consequences of this for melanomagenesis have never been explored.MethodsWe quantified the expression of SLC7A11 and other genes that are involved in the synthesis of pheomelanin but do not regulate the transport of cysteine from the extracellular medium to the cytosol (CTNS, MC1R, ASIP and SLC45A2) in non-tumorous skin of 45 patients of cutaneous melanoma and 50 healthy individuals. We controlled for the effects of Fitzpatrick skin type, age, gender, body mass, frequency of sun exposure and sunburns and number of melanocytic nevi, as well as for the intrinsic antioxidant capacity as given by the expression of the gene NFE2L2.ResultsThe expression of SLC7A11, but not of the other genes, was significantly higher in melanoma patients than in healthy individuals. This was independent of phenotypic factors and antioxidant capacity, thus supporting an effect of pheomelanin-induced oxidative stress on melanomagenesis.ConclusionOur findings indicate that SLC7A11 downregulation in normal epidermal melanocytes may represent a preventive treatment against melanoma.     

Change in DNA Methylation Patterns of SLC6A4 Gene in the Gastric Mucosa in Functional Dyspepsia

Background

The neurochemical serotonin (5-HT) is an important signaling molecule in the gastrointestinal motor and sensory functions. A key regulator of 5-HT levels is the transmembrane serotonin transporter (5-HTT; SLC6A4) that governs the reuptake of 5-HT. Recent studies have indicated 5-HTT expression may be regulated by epigenetic mechanisms. We investigated DNA methylation status of SLC6A4 gene in the gastric mucosa from functional dyspepsia (FD) because of their potential role in dyspeptic symptoms.

Methods

Endoscopic gastric biopsies were obtained from 78 subjects with no upper abdominal symptoms and 79 patients with FD. Bisulfite Pyrosequencing was carried out to determine the methylation status of promoter CpG islands (PCGIs), promoter non-CpG islands (PNCGIs) and gene body non-CpG islands (NPNCGIs) in the SLC6A4 gene. Gene expression was examined by real-time PCR.

Results

In overall, methylation level of PCGIs was significantly lower in FD compared to control subjects (p = 0.04). On the other hand, methylation level of NPNCGIs was significantly higher in FD compared to control subjects (p = 0.03). Lower methylation level in PNCGIs was highlighted in the patients with PDS (p = 0.01), while higher methylation level in NPNCGIs was more prominent in the patients with EPS (p = 0.017). Methylation levels of PCGIs and PNCGIs were inversely correlated, while methylation levels of NPNCGIs was positively correlated with SLC6A4 mRNA levels in FD patients.

Conclusions

Our data suggest that change in DNA methylation pattern of SLC6A4 in the gastric mucosa may have a role for developing FD. A role of epigenetics for developing FD needs to be further evaluated.     

Gene Expression in Chicken Reveals Correlation with Structural Genomic Features and Conserved Patterns of Transcription in the Terrestrial Vertebrates

Background

The chicken is an important agricultural and avian-model species. A survey of gene expression in a range of different tissues will provide a benchmark for understanding expression levels under normal physiological conditions in birds. With expression data for birds being very scant, this benchmark is of particular interest for comparative expression analysis among various terrestrial vertebrates.

Methodology/Principal Findings

We carried out a gene expression survey in eight major chicken tissues using whole genome microarrays. A global picture of gene expression is presented for the eight tissues, and tissue specific as well as common gene expression were identified. A Gene Ontology (GO) term enrichment analysis showed that tissue-specific genes are enriched with GO terms reflecting the physiological functions of the specific tissue, and housekeeping genes are enriched with GO terms related to essential biological functions. Comparisons of structural genomic features between tissue-specific genes and housekeeping genes show that housekeeping genes are more compact. Specifically, coding sequence and particularly introns are shorter than genes that display more variation in expression between tissues, and in addition intergenic space was also shorter. Meanwhile, housekeeping genes are more likely to co-localize with other abundantly or highly expressed genes on the same chromosomal regions. Furthermore, comparisons of gene expression in a panel of five common tissues between birds, mammals and amphibians showed that the expression patterns across tissues are highly similar for orthologuous genes compared to random gene pairs within each pair-wise comparison, indicating a high degree of functional conservation in gene expression among terrestrial vertebrates.

Conclusions

The housekeeping genes identified in this study have shorter gene length, shorter coding sequence length, shorter introns, and shorter intergenic regions, there seems to be selection pressure on economy in genes with a wide tissue distribution, i.e. these genes are more compact. A comparative analysis showed that the expression patterns of orthologous genes are conserved in the terrestrial vertebrates during evolution.