Tandem chimeric transcription as a means to increase protein diversity in the human genome

The “one-gene, one-protein” rule, coined by Beadle and Tatum, has been fundamental to molecular biology. The rule implies that the genetic complexity of an organism depends essentially on its gene number. The discovery, however, that alternative gene splicing and transcription are widespread phenomena dramatically altered our understanding of the genetic complexity of higher eukaryotic organisms; in these, a limited number of genes may potentially encode a much larger number of proteins. Here we investigate yet another phenomenon that may contribute to generate additional protein diversity. Indeed, by relying on both computational and experimental analysis, we estimate that at least 4%–5% of the tandem gene pairs in the human genome can be eventually transcribed into a single RNA sequence encoding a putative chimeric protein. While the functional significance of most of these chimeric transcripts remains to be determined, we provide strong evidence that this phenomenon does not correspond to mere technical artifacts and that it is a common mechanism with the potential of generating hundreds of additional proteins in the human genome.

Gene Promoter Evolution Targets the Center of the Human Protein Interaction Network

Assessing the contribution of promoters and coding sequences to gene evolution is an important step toward discovering the major genetic determinants of human evolution. Many specific examples have revealed the evolutionary importance of cis-regulatory regions. However, the relative contribution of regulatory and coding regions to the evolutionary process and whether systemic factors differentially influence their evolution remains unclear. To address these questions, we carried out an analysis at the genome scale to identify signatures of positive selection in human proximal promoters. Next, we examined whether genes with positively selected promoters (Prom+ genes) show systemic differences with respect to a set of genes with positively selected protein-coding regions (Cod+ genes). We found that the number of genes in each set was not significantly different (8.1% and 8.5%, respectively). Furthermore, a functional analysis showed that, in both cases, positive selection affects almost all biological processes and only a few genes of each group are located in enriched categories, indicating that promoters and coding regions are not evolutionarily specialized with respect to gene function. On the other hand, we show that the topology of the human protein network has a different influence on the molecular evolution of proximal promoters and coding regions. Notably, Prom+ genes have an unexpectedly high centrality when compared with a reference distribution (P = 0.008, for Eigenvalue centrality). Moreover, the frequency of Prom+ genes increases from the periphery to the center of the protein network (P = 0.02, for the logistic regression coefficient). This means that gene centrality does not constrain the evolution of proximal promoters, unlike the case with coding regions, and further indicates that the evolution of proximal promoters is more efficient in the center of the protein network than in the periphery. These results show that proximal promoters have had a systemic contribution to human evolution by increasing the participation of central genes in the evolutionary process.

Hepatocyte nuclear factor-4 alpha regulates the human apolipoprotein AV gene: Identification of a novel response element and involvement in the control by peroxisome proliferator-activated receptor-gamma coactivator-1 alpha, AMP-activated protein kinase, and mitogen-activated protein kinase pathway

The recently discovered apolipoprotein AV (apoAV) gene has been reported to be a key player in modulating plasma triglyceride levels. Here we identify the hepatocyte nuclear factor-4 (HNF-4 ) as a novel regulator of human apoAV gene. Inhibition of HNF-4 expression by small interfering RNA resulted in down-regulation of apoAV. Deletion, mutagenesis, and binding assays revealed that HNF-4 directly regulates human apoAV promoter through DR1 [a direct repeat separated by one nucleotide (nt)], and via a novel element for HNF-4 consisting of an inverted repeat separated by 8 nt (IR8). In addition, we show that the coactivator peroxisome proliferator-activated receptor- coactivator-1 was capable of stimulating the HNF-4 -dependent transactivation of apoAV promoter. Furthermore, analyses in human hepatic cells demonstrated that AMP-activated protein kinase (AMPK) and the MAPK signaling pathway regulate human apoAV expression and suggested that this regulation may be mediated, at least in part, by changes in HNF-4 . Intriguingly, EMSAs and mice with a liver-specific disruption of the HNF-4 gene revealed a species-distinct regulation of apoAV by HNF-4 , which resembles that of a subset of HNF-4 target genes. Taken together, our data provide new insights into the binding properties and the modulation of HNF-4 and underscore the role of HNF-4 in regulating triglyceride metabolism.

Barrier-protective effects of activated protein C in human alveolar epithelial cells

Acute lung injury (ALI) is a clinical manifestation of respiratory failure, caused by lung inflammation and the disruption of the alveolar-capillary barrier. Preservation of the physical integrity of the alveolar epithelial monolayer is of critical importance to prevent alveolar edema. Barrier integrity depends largely on the balance between physical forces on cell-cell and cell-matrix contacts, and this balance might be affected by alterations in the coagulation cascade in patients with ALI. We aimed to study the effects of activated protein C (APC) on mechanical tension and barrier integrity in human alveolar epithelial cells (A549) exposed to thrombin. Cells were pretreated for 3 h with APC (50 mg/ml) or vehicle (control). Subsequently, thrombin (50 nM) or medium was added to the cell culture. APC significantly reduced thrombin-induced cell monolayer permeability, cell stiffening, and cell contraction, measured by electrical impedance, optical magnetic twisting cytometry, and traction microscopy, respectively, suggesting a barrier-protective response. The dynamics of the barrier integrity was also assessed by western blotting and immunofluorescence analysis of the tight junction ZO-1. Thrombin resulted in more elongated ZO-1 aggregates at cell-cell interface areas and induced an increase in ZO-1 membrane protein content. APC attenuated the length of these ZO-1 aggregates and reduced the ZO-1 membrane protein levels induced by thrombin. In conclusion, pretreatment with APC reduced the disruption of barrier integrity induced by thrombin, thus contributing to alveolar epithelial barrier protection.

Fusion of the human gene for the polyubiquitination co-effector UEV-1 with Kua, a newly identified gene

UEV proteins are enzymatically inactive variants of the E2 ubiquitin-conjugating enzymes that regulate noncanonical elongation of ubiquitin chains. In Saccharomyces cerevisiae, UEV is part of the RAD6-mediated error-free DNA repair pathway. In mammalian cells, UEV proteins can modulate c-FOS transcription and the G2-M transition of the cell cycle. Here we show that the UEV genes from phylogenetically distant organisms present a remarkable conservation in their exon–intron structure. We also show that the human UEV1 gene is fused with the previously unknown gene Kua. In Caenorhabditis elegans and Drosophila melanogaster, Kua and UEV are in separated loci, and are expressed as independent transcripts and proteins. In humans, Kua and UEV1 are adjacent genes, expressed either as separate transcripts encoding independent Kua and UEV1 proteins, or as a hybrid Kua–UEV transcript, encoding a two-domain protein. Kua proteins represent a novel class of conserved proteins with juxtamembrane histidine-rich motifs. Experiments with epitope-tagged proteins show that UEV1A is a nuclear protein, whereas both Kua and Kua–UEV localize to cytoplasmic structures, indicating that the Kua domain determines the cytoplasmic localization of Kua–UEV. Therefore, the addition of a Kua domain to UEV in the fused Kua–UEV protein confers new biological properties to this regulator of variant polyubiquitination.[Kua cDNAs isolated by RT-PCR and described in this paper have been deposited in the GenBank data library under accession nos. AF1155120 (H. sapiens) and AF152361 (D. melanogaster). Genomic clones containing UEV genes: S. cerevisiae, YGL087c (accession no. Z72609); S. pombe, c338 (accession no. AL023781); P. falciparum, MAL3P2 (accession no. AL034558); A. thaliana, F26F24 (accession no. AC005292); C. elegans, F39B2 (accession no. Z92834); D. melanogaster, AC014908; and H. sapiens, 1185N5 (accession no. AL034423). Accession numbers for Kua cDNAs in GenBank dbEST: M. musculus, AA7853; T. cruzi, AI612534. Other Kua-containing sequences: A. thaliana genomic clones F10M23 (accession no. AL035440), F19K23 (accession no. AC000375), and T20K9 (accession no. AC004786).

Comparative gene prediction in human and mouse

The completion of the sequencing of the mouse genome promises to help predict human genes with greater accuracy. While current ab initio gene prediction programs are remarkably sensitive (i.e., they predict at least a fragment of most genes), their specificity is often low, predicting a large number of false-positive genes in the human genome. Sequence conservation at the protein level with the mouse genome can help eliminate some of those false positives. Here we describe SGP2, a gene prediction program that combines ab initio gene prediction with TBLASTX searches between two genome sequences to provide both sensitive and specific gene predictions. The accuracy of SGP2 when used to predict genes by comparing the human and mouse genomes is assessed on a number of data sets, including single-gene data sets, the highly curated human chromosome 22 predictions, and entire genome predictions from ENSEMBL. Results indicate that SGP2 outperforms purely ab initio gene prediction methods. Results also indicate that SGP2 works about as well with 3x shotgun data as it does with fully assembled genomes. SGP2 provides a high enough specificity that its predictions can be experimentally verified at a reasonable cost. SGP2 was used to generate a complete set of gene predictions on both the human and mouse by comparing the genomes of these two species. Our results suggest that another few thousand human and mouse genes currently not in ENSEMBL are worth verifying experimentally.

Signatures of selection in the human olfactory receptor OR5I1 gene

The human olfactory receptor repertoire is reduced in comparison to other mammalsand to other non-human primates. Nonetheless, this olfactory decline opens an opportunity forevolutionary innovation and improvement. In the present study, we focus on an olfactoryreceptor gene, OR5I1, which had previously been shown to present an excess of amino acidreplacement substitutions between humans and chimpanzees. We analyze the geneticvariation in OR5I1 in a large worldwide human panel and find an excess of derived allelessegregating at relatively high frequencies in all populations. Additional evidence for selectionincludes departures from neutrality in allele frequency spectra tests but no unusually extendedhaplotype structure. Moreover, molecular structural inference suggests that one of thenonsynonymous polymorphisms defining the presumably adaptive protein form of OR5I1may alter the functional binding properties of the olfactory receptor. These results arecompatible with positive selection having modeled the pattern of variation found in the OR5I1gene and with a relatively ancient, mild selective sweep predating the “Out of Africa”expansion of modern humans.

EGASP: the human ENCODE Genome Annotation Assessment Project

Background: We present the results of EGASP, a community experiment to assess the state-ofthe-art in genome annotation within the ENCODE regions, which span 1% of the human genomesequence. The experiment had two major goals: the assessment of the accuracy of computationalmethods to predict protein coding genes; and the overall assessment of the completeness of thecurrent human genome annotations as represented in the ENCODE regions. For thecomputational prediction assessment, eighteen groups contributed gene predictions. Weevaluated these submissions against each other based on a ‘reference set’ of annotationsgenerated as part of the GENCODE project. These annotations were not available to theprediction groups prior to the submission deadline, so that their predictions were blind and anexternal advisory committee could perform a fair assessment.Results: The best methods had at least one gene transcript correctly predicted for close to 70%of the annotated genes. Nevertheless, the multiple transcript accuracy, taking into accountalternative splicing, reached only approximately 40% to 50% accuracy. At the coding nucleotidelevel, the best programs reached an accuracy of 90% in both sensitivity and specificity. Programsrelying on mRNA and protein sequences were the most accurate in reproducing the manuallycurated annotations. Experimental validation shows that only a very small percentage (3.2%) of the selected 221 computationally predicted exons outside of the existing annotation could beverified.Conclusions: This is the first such experiment in human DNA, and we have followed thestandards established in a similar experiment, GASP1, in Drosophila melanogaster. We believe theresults presented here contribute to the value of ongoing large-scale annotation projects and shouldguide further experimental methods when being scaled up to the entire human genome sequence.

A short motif in Drosophila SECIS Binding Protein 2 provides differential binding affinity to SECIS RNA hairpins

Selenoproteins contain the amino acid selenocysteine which is encoded by a UGA Sec codon. Recoding UGA Sec requires a complex mechanism, comprising the cis-acting SECIS RNA hairpin in the 3′UTR of selenoprotein mRNAs, and trans-acting factors. Among these, the SECIS Binding Protein 2 (SBP2) is central to the mechanism. SBP2 has been so far functionally characterized only in rats and humans. In this work, we report the characterization of the Drosophila melanogaster SBP2 (dSBP2). Despite its shorter length, it retained the same selenoprotein synthesis-promoting capabilities as the mammalian counterpart. However, a major difference resides in the SECIS recognition pattern: while human SBP2 (hSBP2) binds the distinct form 1 and 2 SECIS RNAs with similar affinities, dSBP2 exhibits high affinity toward form 2 only. In addition, we report the identification of a K (lysine)-rich domain in all SBP2s, essential for SECIS and 60S ribosomal subunit binding, differing from the well-characterized L7Ae RNA-binding domain. Swapping only five amino acids between dSBP2 and hSBP2 in the K-rich domain conferred reversed SECIS-binding properties to the proteins, thus unveiling an important sequence for form 1 binding.

Identifying protein-coding genes in genomic sequences

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

Comparative gene finding in chicken indicates that we are closing in on the set of multi-exonic widely-expressed human genes

The recent availability of the chicken genome sequence poses the question of whether there are human protein-coding genes conserved in chicken that are currently not included in the human gene catalog. Here, we show, using comparative gene finding followed by experimental verification of exon pairs by RT–PCR, that the addition to the multi-exonic subset of this catalog could be as little as 0.2%, suggesting that we may be closing in on the human gene set. Our protocol, however, has two shortcomings: (i) the bioinformatic screening of the predicted genes, applied to filter out false positives, cannot handle intronless genes; and (ii) the experimental verification could fail to identify expression at a specific developmental time. This highlights the importance of developing methods that could provide a reliable estimate of the number of these two types of genes.

Design and evaluation of a panel os single-nucleotide polymorphisms in microRNA genomic regions for association studies in human disease

MicroRNAs (miRNA) are recognized posttranscriptional gene repressors involved in the control of almost every biological process. Allelic variants in these regions may be an important source of phenotypic diversity and contribute to disease susceptibility. We analyzed the genomic organization of 325 human miRNAs (release 7.1, miRBase) to construct a panel of 768 single-nucleotide polymorphisms (SNPs) covering approximately 1 Mb of genomic DNA, including 131 isolated miRNAs (40%) and 194 miRNAs arranged in 48 miRNA clusters, as well as their 5-kb flanking regions. Of these miRNAs, 37% were inside known protein-coding genes, which were significantly associated with biological functions regarding neurological, psychological or nutritional disorders. SNP coverage analysis revealed a lower SNP density in miRNAs compared with the average of the genome, with only 24 SNPs located in the 325 miRNAs studied. Further genotyping of 340 unrelated Spanish individuals showed that more than half of the SNPs in miRNAs were either rare or monomorphic, in agreement with the reported selective constraint on human miRNAs. A comparison of the minor allele frequencies between Spanish and HapMap population samples confirmed the applicability of this SNP panel to the study of complex disorders among the Spanish population, and revealed two miRNA regions, hsa-mir-26a-2 in the CTDSP2 gene and hsa-mir-128-1 in the R3HDM1 gene, showing geographical allelic frequency variation among the four HapMap populations, probably because of differences in natural selection. The designed miRNA SNP panel could help to identify still hidden links between miRNAs and human disease.

Inhibition of Mitogen-Activated Protein Kinase Erk1/2 Promotes Protein Degradation of ATP Binding Cassette Transporters A1 and G1 in CHO and in HuH7 cells.

Signal transduction modulates expression and activity of cholesterol transporters. We recently demonstrated that the Ras/mitogen-activated protein kinase (MAPK) signaling cascade regulates protein stability of Scavenger Receptor BI (SR-BI) through Proliferator Activator Receptor (PPARα) -dependent degradation pathways. In addition, MAPK (Mek/Erk 1/2) inhibition has been shown to influence liver X receptor (LXR) -inducible ATP Binding Cassette (ABC) transporter ABCA1 expression in macrophages. Here we investigated if Ras/MAPK signaling could alter expression and activity of ABCA1 and ABCG1 in steroidogenic and hepatic cell lines. We demonstrate that in Chinese Hamster Ovary (CHO) cells and human hepatic HuH7 cells, extracellular signal-regulated kinase 1/2 (Erk1/2) inhibition reduces PPARα-inducible ABCA1 protein levels, while ectopic expression of constitutively active H-Ras, K-Ras and MAPK/Erk kinase 1 (Mek1) increases ABCA1 protein expression, respectively. Furthermore, Mek1/2 inhibitors reduce ABCG1 protein levels in ABCG1 overexpressing CHO cells (CHO-ABCG1) and human embryonic kidney 293 (HEK293) cells treated with LXR agonist. This correlates with Mek1/2 inhibition reducing ABCG1 cell surface expression and decreasing cholesterol efflux onto High Density Lipoproteins (HDL). Real Time reverse transcriptase polymerase chain reaction (RT-PCR) and protein turnover studies reveal that Mek1/2 inhibitors do not target transcriptional regulation of ABCA1 and ABCG1, but promote ABCA1 and ABCG1 protein degradation in HuH7 and CHO cells, respectively. In line with published data from mouse macrophages, blocking Mek1/2 activity upregulates ABCA1 and ABCG1 protein levels in human THP1 macrophages, indicating opposite roles for the Ras/MAPK pathway in the regulation of ABC transporter activity in macrophages compared to steroidogenic and hepatic cell types. In summary, this study suggests that Ras/MAPK signaling modulates PPARα- and LXR-dependent protein degradation pathways in a cell-specific manner to regulate the expression levels of ABCA1 and ABCG1 transporters.

Insulin degradation by adipose tissue is increased in human obesity

White adipose tissue samples from obese and lean patients were used for the estimation ofinsulin protease and insulin:glutathione transhydrogenase using 1251-labeled insulin. There was no activity detected in the absence of reduced glutathione, which indicates that insulin is cleaved in human adipose "tissue through reduction of the disulfide bridge between the chains. O bese patients showed higher transhydrogenase activity (per U tissue protein wt, per U tissue wt, and in the total adipose tissue mass) than the lean group. There is a significant correlation between the activity per U tissue wt, and protein and total activity in the whole adipose tissue with respect to body mass index, with a higher activity in obese patients. The potential ofinsulin cleavage by adipose tissue in obese patients was a mean 5.6-fold higher than that in controla. The coexistence of high insulinemia and high cleavage capability implies that insulin secretion and turnover are increased in the o bese. Thus, white adipose tissue may be crucial in the control of energy availability through modulation ofinsulin cleavage.

Different effects of hyperlipidic diets in human lactation and adulthood: growth versus the development of obesity

After birth, the body shifts from glucose as primary energy substrate to milk-derived fats, with sugars from lactose taking a secondary place. At weaning, glucose recovers its primogeniture and dietary fat role decreases. In spite of human temporary adaptation to a high-fat (and sugars and protein) diet during lactation, the ability to thrive on this type of diet is lost irreversibly after weaning. We could not revert too the lactating period metabolic setting because of different proportions of brain/muscle metabolism in the total energy budget, lower thermogenesis needs and capabilities, and absence of significant growth in adults. A key reason for change was the limited availability of foods with high energy content at weaning and during the whole adult life of our ancestors, which physiological adaptations remain practically unchanged in our present-day bodies. Humans have evolved to survive with relatively poor diets interspersed by bouts of scarcity and abundance. Today diets in many societies are largely made up from choice foods, responding to our deeply ingrained desire for fats, protein, sugars, salt etc. Consequently our diets are not well adjusted to our physiological needs/adaptations but mainly to our tastes (another adaptation to periodic scarcity), and thus are rich in energy roughly comparable to milk. However, most adult humans cannot process the food ingested in excess because our cortical-derived craving overrides the mechanisms controlling appetite. This is produced not because we lack the biochemical mechanisms to use this energy, but because we are unprepared for excess, and wholly adapted to survive scarcity. The thrifty mechanisms compound the effects of excess nutrients and damage the control of energy metabolism, developing a pathologic state. As a consequence, an overflow of energy is generated and the disease of plenty develops.