Structural and functional diversity of connexin genes in the mouse and human genome

Gap junctions are clustered channels between contacting cells through which direct intercellular communication via diffusion of ions and metabolites can occur. Two hemichannels, each built up of six connexin protein subunits in the plasma membrane of adjacent cells, can dock to each other to form conduits between cells. We have recently screened mouse and human genomic data bases and have found 19 connexin (Cx) genes in the mouse genome and 20 connexin genes in the human genome. One mouse connexin gene and two human connexin genes do not appear to have orthologs in the other genome. With three exceptions, the characterized connexin genes comprise two exons whereby the complete reading frame is located on the second exon. Targeted ablation of eleven mouse connexin genes revealed basic insights into the functional diversity of the connexin gene family. In addition, the phenotypes of human genetic disorders caused by mutated connexin genes further complement our understanding of connexin functions in the human organism. In this review we compare currently identified connexin genes in both the mouse and human genome and discuss the functions of gap junctions deduced from targeted mouse mutants and human genetic disorders.

Connexin genes in the mouse and human genome

Gap junctions serve for direct intercellular communication by docking of two hemichannels in adjacent cells thereby forming conduits between the cytoplasmic compartments of adjacent cells. Connexin genes code for subunit proteins of gap junction channels and are members of large gene families in mammals. So far, 17 connexin (Cx) genes have been described and characterized in the murine genome. For most of them, orthologues in the human genome have been found (see White and Paul 1999; Manthey et al. 1999; Teubner et al. 2001; Söhl et al. 2001). We have recently performed searches for connexin genes in murine and human gene libraries available at EMBL/Heidelberg, NCBI and the Celera company that have increased the number of identified connexins to 19 in mouse and 20 in humans. For one mouse connexin gene and two human connexin genes we did not find orthologues in the other genome. Here we present a short overview on distinct connexin genes which we found in the mouse and human genome and which may include all members of this gene family, if no further connexin gene will be discovered in the remaining non-sequenced parts (about 1-5%) of the genomes.

A highly conserved interspecies V(H) in the human genome.

Idiotype conservation between human and mouse antibodies has been observed in association with various infectious and autoimmune diseases. We have isolated a human anti-idiotypic antibody to a mouse monoclonal anti-IgE antibody (BSW17) suggesting a conserved interspecies idiotype associated with an anti-IgE response. To find the homologue of BSW17 in the human genome we applied the guided selection strategy. Combining V(H) of BSW17 with a human V(L) repertoire resulted in three light chains. The three V(L) chains were then combined with a human V(H) repertoire resulting in three clones specific for human IgE. Surprisingly, one clone, Hu41, had the same epitope specificity and functional in vitro activity as BSW17 and V(H) complementarity-determining regions identical with that of BSW17. Real-time PCR analysis confirmed the presence of the Hu41 V(H) sequence in the human genome. These data document the first example of the isolation of a human antibody where high sequence similarity to the original murine V(H) sequence is associated with common antigen and epitope specificity.

Global diversity, population stratification and selection of human copy number variation

In order to explore the diversity and selective signatures of duplication and deletion human copy number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single nucleotide variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load.

Epigenetic regulation of somatic angiotensin-converting enzyme by DNA methylation and histone acetylation

Somatic angiotensin-converting enzyme (sACE) is crucial in cardiovascular homeostasis and displays a tissue-specific profile. Epigenetic patterns modulate genes expression and their alterations were implied in pathologies including hypertension. However, the influence of DNA methylation and chromatin condensation state on the expression of sACE is unknown. We examined whether such epigenetic mechanisms could participate in the control of sACE expression in vitro and in vivo. We identified two CpG islands in the human ace-1 gene 3 kb proximal promoter region. Their methylation abolished the luciferase activity of ace-1 promoter/reporter constructs transfected into human liver (HepG2), colon (HT29), microvascular endothelial (HMEC-1) and lung (SUT) cell lines (p < 0.001). Bisulphite sequencing revealed a cell-type specific basal methylation pattern of the ace-1 gene -1,466/+25 region. As assessed by RT-qPCR, inhibition of DNA methylation by 5-aza-2'-deoxycytidine and/or of histone deacetylation by trichostatin A highly stimulated sACE mRNA expression cell-type specifically (p < 0.001 vs. vehicle treated cells). In the rat, in vivo 5-aza-cytidine injections demethylated the ace-1 promoter and increased sACE mRNA expression in the lungs and liver (p = 0.05), but not in the kidney. In conclusion, the expression level of somatic ACE is modulated by CpG-methylation and histone deacetylases inhibition. The basal methylation pattern of the promoter of the ace-1 gene is cell-type specific and correlates to sACE transcription. DNMT inhibition is associated with altered methylation of the ace-1 promoter and a cell-type and tissue-specific increase of sACE mRNA levels. This study indicates a strong influence of epigenetic mechanisms on sACE expression.

Evolutionary forces shaping genomic islands of population differentiation in humans

Background Levels of differentiation among populations depend both on demographic and selective factors: genetic drift and local adaptation increase population differentiation, which is eroded by gene flow and balancing selection. We describe here the genomic distribution and the properties of genomic regions with unusually high and low levels of population differentiation in humans to assess the influence of selective and neutral processes on human genetic structure. Methods Individual SNPs of the Human Genome Diversity Panel (HGDP) showing significantly high or low levels of population differentiation were detected under a hierarchical-island model (HIM). A Hidden Markov Model allowed us to detect genomic regions or islands of high or low population differentiation. Results Under the HIM, only 1.5% of all SNPs are significant at the 1% level, but their genomic spatial distribution is significantly non-random. We find evidence that local adaptation shaped high-differentiation islands, as they are enriched for non-synonymous SNPs and overlap with previously identified candidate regions for positive selection. Moreover there is a negative relationship between the size of islands and recombination rate, which is stronger for islands overlapping with genes. Gene ontology analysis supports the role of diet as a major selective pressure in those highly differentiated islands. Low-differentiation islands are also enriched for non-synonymous SNPs, and contain an overly high proportion of genes belonging to the 'Oncogenesis' biological process. Conclusions Even though selection seems to be acting in shaping islands of high population differentiation, neutral demographic processes might have promoted the appearance of some genomic islands since i) as much as 20% of islands are in non-genic regions ii) these non-genic islands are on average two times shorter than genic islands, suggesting a more rapid erosion by recombination, and iii) most loci are strongly differentiated between Africans and non-Africans, a result consistent with known human demographic history.

Whole genome methylation array analysis reveals new aspects in Balkan endemic nephropathy etiology

Background Balkan endemic nephropathy (BEN) represents a chronic progressive interstitial nephritis in striking correlation with uroepithelial tumours of the upper urinary tract. The disease has endemic distribution in the Danube river regions in several Balkan countries. DNA methylation is a primary epigenetic modification that is involved in major processes such as cancer, genomic imprinting, gene silencing, etc. The significance of CpG island methylation status in normal development, cell differentiation and gene expression is widely recognized, although still stays poorly understood. Methods We performed whole genome DNA methylation array analysis on DNA pool samples from peripheral blood from 159 affected individuals and 170 healthy individuals. This technique allowed us to determine the methylation status of 27 627 CpG islands throughout the whole genome in healthy controls and BEN patients. Thus we obtained the methylation profile of BEN patients from Bulgarian and Serbian endemic regions. Results Using specifically developed software we compared the methylation profiles of BEN patients and corresponding controls and revealed the differently methylated regions. We then compared the DMRs between all patient-control pairs to determine common changes in the epigenetic profiles. SEC61G, IL17RA, HDAC11 proved to be differently methylated throughout all patient-control pairs. The CpG islands of all 3 genes were hypomethylated compared to controls. This suggests that dysregulation of these genes involved in immunological response could be a common mechanism in BEN pathogenesis in both endemic regions and in both genders. Conclusion Our data propose a new hypothesis that immunologic dysregulation has a place in BEN etiopathogenesis. Keywords: Epigenetics; Whole genome array analysis; Balkan endemic nephropathy

Comparative human-horse sequence analysis of the CYP3A subfamily gene cluster

Cytochrome P450 enzymes (CYP450s) represent a superfamily of haem-thiolate proteins. CYP450s are most abundant in the liver, a major site of drug metabolism, and play key roles in the metabolism of a variety of substrates, including drugs and environmental contaminants. Interaction of two or more different drugs with the same enzyme can account for adverse effects and failure of therapy. Human CYP3A4 metabolizes about 50% of all known drugs, but little is known about the orthologous CYP450s in horses. We report here the genomic organization of the equine CYP3A gene cluster as well as a comparative analysis with the human CYP3A gene cluster. The equine CYP450 genes of the 3A family are located on ECA 13 between 6.97-7.53 Mb, in a region syntenic to HSA 7 99.05-99.35 Mb. Seven potential, closely linked equine CYP3A genes were found, in contrast to only four genes in the human genome. RNA was isolated from an equine liver sample, and the approximately 1.5-kb coding sequence of six CYP3A genes could be amplified by RT-PCR. Sequencing of the RT-PCR products revealed numerous hitherto unknown single nucleotide polymorphisms (SNPs) in these six CYP3A genes, and one 6-bp deletion compared to the reference sequence (EquCab2.0). The presence of the variants was confirmed in a sample of genomic DNA from the same horse. In conclusion, orthologous genes for the CYP3A family exist in horses, but their number differs from those of the human CYP3A gene family. CYP450 genes of the same family show high homology within and between mammalian species, but can be highly polymorphic.

Distinct DNA methylation patterns in cirrhotic liver and hepatocellular carcinoma

Abberrant DNA methylation is one of the hallmarks of cancerogenesis. Our study aims to delineate differential DNA methylation in cirrhosis and hepatic cancerogenesis. Patterns of methylation of 27,578 individual CpG loci in 12 hepatocellular carcinomas (HCCs), 15 cirrhotic controls and 12 normal liver samples were investigated using an array-based technology. A supervised principal component analysis (PCA) revealed 167 hypomethylated loci and 100 hypermethylated loci in cirrhosis and HCC as compared to normal controls. Thus, these loci show a "cirrhotic" methylation pattern that is maintained in HCC. In pairwise supervised PCAs between normal liver, cirrhosis and HCC, eight loci were significantly changed in all analyses differentiating the three groups (p < 0.0001). Of these, five loci showed highest methylation levels in HCC and lowest in control tissue (LOC55908, CELSR1, CRMP1, GNRH2, ALOX12 and ANGPTL7), whereas two loci showed the opposite direction of change (SPRR3 and TNFSF15). Genes hypermethylated between normal liver to cirrhosis, which maintain this methylation pattern during the development of HCC, are depleted for CpG islands, high CpG content promoters and polycomb repressive complex 2 (PRC2) targets in embryonic stem cells. In contrast, genes selectively hypermethylated in HCC as compared to nonmalignant samples showed an enrichment of CpG islands, high CpG content promoters and PRC2 target genes (p < 0.0001). Cirrhosis and HCC show distinct patterns of differential methylation with regards to promoter structure, PRC2 targets and CpG islands.

A human-horse comparative map based on equine BAC end sequences

In an effort to increase the density of sequence-based markers for the horse genome we generated 9473 BAC end sequences (BESs) from the CHORI-241 BAC library with an average read length of 677 bp. BLASTN searches with the BESs revealed 4036 meaningful hits (E human genome that provide useful markers for the human-horse comparative map. The 4036 BLASTN hits allowed the anchoring of 3079 BAC clones to the human genome, on average one corresponding equine BAC clone per megabase of human DNA. We used the BLASTN anchored BESs for an in silico prediction of the gene content and chromosome assignment of comparatively mapped equine BAC clones. As a first verification of our in silico mapping strategy we placed 19 equine BESs with matches to HSA6 onto the RH map. All markers were assigned to the predicted localizations on ECA10, ECA20, and ECA31, respectively.

Whole genome RNAi screens reveal a critical role of REV3 in coping with replication stress

REV3, the catalytic subunit of translesion polymerase zeta (polζ), is commonly associated with DNA damage bypass and repair. Despite sharing accessory subunits with replicative polymerase δ, very little is known about the role of polζ in DNA replication. We previously demonstrated that inhibition of REV3 expression induces persistent DNA damage and growth arrest in cancer cells. To reveal determinants of this sensitivity and obtain insights into the cellular function of REV3, we performed whole human genome RNAi library screens aimed at identification of synthetic lethal interactions with REV3 in A549 lung cancer cells. The top confirmed hit was RRM1, the large subunit of ribonucleotide reductase (RNR), a critical enzyme of de novo nucleotide synthesis. Treatment with the RNR-inhibitor hydroxyurea (HU) synergistically increased the fraction of REV3-deficient cells containing single stranded DNA (ssDNA) as indicated by an increase in replication protein A (RPA). However, this increase was not accompanied by accumulation of the DNA damage marker γH2AX suggesting a role of REV3 in counteracting HU-induced replication stress (RS). Consistent with a role of REV3 in DNA replication, increased RPA staining was confined to HU-treated S-phase cells. Additionally, we found genes related to RS to be significantly enriched among the top hits of the synthetic sickness/lethality (SSL) screen further corroborating the importance of REV3 for DNA replication under conditions of RS.

Expression of human cathepsin L or human cathepsin V in mouse thymus mediates positive selection of T helper cells in cathepsin L knock-out mice

A genetic deficiency of the cysteine protease cathepsin L (Ctsl) in mice results in impaired positive selection of conventional CD4+ T helper cells as a result of an incomplete processing of the MHC class II associated invariant chain or incomplete proteolytic generation of positively selecting peptide ligands. The human genome encodes, in contrast to the mouse genome, for two cathepsin L proteases, namely cathepsin L (CTSL) and cathepsin V (CTSV; alternatively cathepsin L2). In the human thymic cortex, CTSV is the predominately expressed protease as compared to CTSL or other cysteine cathepsins. In order to analyze the functions of CTSL and CTSV in the positive selection of CD4+ T cells we employed Ctsl knock-out mice crossed either with transgenic mice expressing CTSL under the control of its genuine human promoter or with transgenic mice expressing CTSV under the control of the keratin 14 (K14) promoter, which drives expression to the cortical epithelium. Both human proteases are expressed in the thymus of the transgenic mice, and independent expression of both CTSL and CTSV rescues the reduced frequency of CD4+ T cells in Ctsl-deficient mice. Moreover, the expression of the human cathepsins does not change the number of CD4+CD25+Foxp3+ regulatory T cells, but the normalization of the frequency of conventional CD4+ T cell in the transgenic mice results in a rebalancing of conventional T cells and regulatory T cells. We conclude that the functional differences of CTSL and CTSV in vivo are not mainly determined by their inherent biochemical properties, but rather by their tissue specific expression pattern.