The role of gene duplication in the origin and evolution of new biological functions

Abstract : Gene duplication is an essential source of material for the origin of genetic novelty and the evolution of lineage- or species-specific phenotypic traits. The reverse transcription of source gene mRNA followed by the genomic insertion of the resulting cDNA - retroposition - has provided the human genome with a significant number of gene copies during the last ~63 million years (MYA) of primate evolution. We estimated that at least 1 new functional gene (retrogene) per MYA emerged by retroposition in the primate lineage leading to humans. Using a combination of comparative sequencing and evolutionary simulations, we obtained strong evidence of functionality for 7 primate specific retrogenes. Most of these genes are specifically expressed in testis suggesting that retroposition has contributed with genetic raw material necessary for the evolution ofmale-specific functions in primates. We characterized CDC14Bretro (identified in the previous survey) that originated from the retroposition of a cell cycle gene - CDC14B - in the common ancestor of humans and apes. We demonstrate that CDC14Bretro experienced a period of intense positive selection in the African ape ancestor. By virtue of the amino acid substitutions that occurred during this period CDC 14Bretro adapted to a new subcellular compartment in African apes. Further analyses indicate that this subcellular shift reflects the evolution of anew functional role of CDC 14Bretro. Prompted by this result, we used yeast (Saccharomyces cerevisiae) to investigate on a global scale the extent of functional diversification of duplicate genes through the subcellular adaptation of their encoded proteins. We found that duplicate proteins frequently evolved new cellular localization patterns, either by partitioning of ancestral localizations ("sublocalization"), or more frequently by relocalization to previously unoccupied compartments ("neolocalization"). Interestingly, proteins involved in processes with a wider subcellular distribution more frequently evolved new localization patterns suggesting that subcellular localization changes are dependent on progenitor gene functions. Relocated proteins adapted to their new subcellular environments and evolved new functional roles through changes of their physio-chemical properties, expression levels, and interaction partners. Our work suggests an important role of subcellular adaptation for the emergence of new gene functions.

Frequent clonal loss of heterozygosity but scarcity of microsatellite instability at chromosomal breakpoint cluster regions in adult leukemias.

Microsatellites are important highly polymorphic genetic markers dispersed in the human genome. Using a panel of 22 (CA)n repeat microsatellite markers mapped to recurrent breakpoint cluster regions specifically involved in leukemia, we investigated 114 adult leukemias (25 acute lymphocytic leukemia [ALL], 32 acute myeloid leukemia [AML], 36 chronic lymphocytic leukemia [CLL], and 21 chronic myeloid leukemia [CML] in chronic phase) for somatic mutations at these loci. In each patient, DNA from fresh leukemia samples was analyzed alongside normal constitutive DNA from buccal epithelium. We detected loss of heterozygosity (LOH) in 81 of 114 patients (ALL 16/25, AML 25/32, CLL 30/36, CML 10/21). Deletions were most often seen in ALL at 11q23 and 19p13; in AML at 8q22 and 11q23; in CLL at 13q14.3, 11q13, and 11q23; and in CML at 3q26. Only six deletions were reported in 74 karyotypes analyzed, whereas in these same cases, 91 LOH events were detected by microsatellites. Of 26 leukemias with a normal karyotype, 16 nevertheless showed at least one LOH by microsatellite analysis. Replication errors were found in 10 of 114 patients (8.8%). Thus, microsatellite instability is rare in leukemia in contrast to many solid tumors. Our findings suggest that in adult leukemia, LOH may be an important genetic event in addition to typical chromosomal translocations. LOH may point to the existence of tumor suppressor genes involved in leukemogenesis to a degree that has hitherto been underestimated.

Characterization of NALP2 and NALP3 IL-1β processing inflammasomes

Résumé II y a cinq ans, la découverte d'un nouveau domaine, le PYD domaine, lié aux domaines de la mort, a permis la description de la nouvelle famille des NALP protéines. L'analyse structurelle de cette famille de protéines révéla la présence de deux autres domaines, impliqués dans l'oligomerisation, NACHT, et la détection des ligands, Leucine rich repeats ou LRR. Cette architecture protéique est homologue à celle qui est décrite pour les NODs, les Tol1 récepteurs et tes protéines de résistance chez les plantes. Cette homologie suggère une possible implication des NALPs dans la régulation de l'immunité innée. Premièrement, nous avons décrit les composants minimaux qui permettent à l'inflammasomeNALP3 d'activer la caspase pro-inflammatoire, caspase-1. En comparaison à NALP1, NALP3 ne contient pas de FIIND domaine, ni de CARD domaine en C-terminus et n'interagit pas avec caspase-5. Nous avons découvert une protéine très homologue au C-terminus de NALP1, Cardinal, qui se lie au NACHT domaine de NALP2 et NALP3 par l'intermédiaire de son FIIND domaine. Cardinal possède la capacité d'interagir avec caspase-l, mais seul ASC semble être nécessaire à la maturation de la prointerleukine-1β suite à la stimulation de NALP3. Deuxièmement, notre étude s'est concentrée sur la nature du stimulus capable d'induire la formation et l'activation de l'inflammasome-NALP3. Nous avons démontré que l'ajout de muramyl dipeptide (MDP), produit à partir de la digestion enzymatique de peptidoglycaris bactériens, induit à la fois l'expression de la proIL-1β par la voie NOD2 et sa maturation en IL-1β active par la voie NALP3. Bien que le MDP active l'inflammasome-NALP3, il est incapable d'induire la sécrétion de l'IL-1β mature dans la lignée cellulaire THP1, comparé aux monocytes primaires humains. Cette différence pourrait être liée à l'absence, dans les THP1, de la protéine Filamin, qui est proposée d'interagir avec Cardinal. L'implication de NALP3 dans la maturation de l'IL-lb est confirmée suite à la découverte de mutations sur le gène CIAS1/NALP3/cryopyrin associées à trois maladies auto-inflammatoires : le syndrome de Muckle-Wells (MWS), l'urticaire familial au froid (FCU) et le syndrome CINCA/NOMID. Une élévation constitutive de la maturation et de la sécrétion de la proIL-1β en absence de stimulation MDP est détectée dans les macrophages des patients Muckle-Wells. En conclusion, nos études ont démontré que l'inflammasome-NALP3 doit être finement régulé pour éviter une activité incontrôlée qui représente la base moléculaire des symptômes associés aux syndromes auto-inflammatoires liés à NALP3. Summary Five years ago, the description of the NALP family originated from the discovery of a new death-domain fold family, the PYD domain. NALP contains aprotein-protein interaction domain (PYD), an oligomerization domain (NACHT) and a ligand-sensing domain, leucine rich repeats or LRR. This protein architecture shares similarity with receptors involved in immunity, such as NODS, Toll receptors (TLRs) and related plant resistance proteins, and points to an important role of NALPs in defense mechanisms. We first described the minimal complex involved in the pro-inflammatory Interleukin-1beta (IL-1β) cytokine maturation, called the inflammasome, which contains NALP3. In contrast to NALP1, NALP3, like other members of the NALP family, is devoid of C-terminal FIIND and CARD domains and does not interact with the pro-inflammatory caspase-5. Interestingly, a homolog of the C-terminal portion of NALP1 was found in the human genome and was named Cardinal. We found that NALP2 and NALP3 interact with the CARD-containing proteins Cardinal. Cardinal is able to bind to caspase-1 but is not required for IL-1β maturation through NALP3 activation, as demonstrated for the adaptor ASC. Secondly, our study focused on the stimuli involved in the activation of the NALP3 inflammasome. MDP was shown to induce the expression of proIL1β through NOD2 and then the maturation into active IL-1β by activation of the NALP3 inflammasome. However, in the monocytic THP1 cell line, secretion of IL-1β upon MDP stimulation seems to be independent of the inflammasome activation compared to human primary monocytes. This difference might be linked to a Cardinal-interacting protein, filamin. Until now, the role of Cardinal and filamin is still unknown and remains to be elucidated. Finally, mutations in the NALP3/cryopyrin/CIAS1 gene are associated with three autoinflammatory diseases: Muckle-Wells syndrome, familial cold autoinflammatory syndrome, and CINCA. Constitutive, elevated IL-1β maturation and secretion, even in the absence of MDP stimulation, was observed in macrophages from Muckle-Wells patients and confirmed a key role for the NALP3 inflammasome in innate immunity In conclusion, our studies describes the formation of the NALP3 inflammasome and suggests that this complex has to be tightly regulated to avoid an increased deregulated inflammasome activity that is the molecular basis for the symptoms associated with NALP3-dependent autoinflammatory disorders.

A fish-specific transposable element shapes the repertoire of p53 target genes in zebrafish.

Transposable elements, as major components of most eukaryotic organisms' genomes, define their structural organization and plasticity. They supply host genomes with functional elements, for example, binding sites of the pleiotropic master transcription factor p53 were identified in LINE1, Alu and LTR repeats in the human genome. Similarly, in this report we reveal the role of zebrafish (Danio rerio) EnSpmN6_DR non-autonomous DNA transposon in shaping the repertoire of the p53 target genes. The multiple copies of EnSpmN6_DR and their embedded p53 responsive elements drive in several instances p53-dependent transcriptional modulation of the adjacent gene, whose human orthologs were frequently previously annotated as p53 targets. These transposons define predominantly a set of target genes whose human orthologs contribute to neuronal morphogenesis, axonogenesis, synaptic transmission and the regulation of programmed cell death. Consistent with these biological functions the orthologs of the EnSpmN6_DR-colonized loci are enriched for genes expressed in the amygdala, the hippocampus and the brain cortex. Our data pinpoint a remarkable example of convergent evolution: the exaptation of lineage-specific transposons to shape p53-regulated neuronal morphogenesis-related pathways in both a hominid and a teleost fish.

A user's guide to the encyclopedia of DNA elements (ENCODE).

The mission of the Encyclopedia of DNA Elements (ENCODE) Project is to enable the scientific and medical communities to interpret the human genome sequence and apply it to understand human biology and improve health. The ENCODE Consortium is integrating multiple technologies and approaches in a collective effort to discover and define the functional elements encoded in the human genome, including genes, transcripts, and transcriptional regulatory regions, together with their attendant chromatin states and DNA methylation patterns. In the process, standards to ensure high-quality data have been implemented, and novel algorithms have been developed to facilitate analysis. Data and derived results are made available through a freely accessible database. Here we provide an overview of the project and the resources it is generating and illustrate the application of ENCODE data to interpret the human genome.

Identification of novel functional TBP-binding sites and general factor repertoires.

Our current knowledge of the general factor requirement in transcription by the three mammalian RNA polymerases is based on a small number of model promoters. Here, we present a comprehensive chromatin immunoprecipitation (ChIP)-on-chip analysis for 28 transcription factors on a large set of known and novel TATA-binding protein (TBP)-binding sites experimentally identified via ChIP cloning. A large fraction of identified TBP-binding sites is located in introns or lacks a gene/mRNA annotation and is found to direct transcription. Integrated analysis of the ChIP-on-chip data and functional studies revealed that TAF12 hitherto regarded as RNA polymerase II (RNAP II)-specific was found to be also involved in RNAP I transcription. Distinct profiles for general transcription factors and TAF-containing complexes were uncovered for RNAP II promoters located in CpG and non-CpG islands suggesting distinct transcription initiation pathways. Our study broadens the spectrum of general transcription factor function and uncovers a plethora of novel, functional TBP-binding sites in the human genome.

Prominent use of distal 5' transcription start sites and discovery of a large number of additional exons in ENCODE regions.

This report presents systematic empirical annotation of transcript products from 399 annotated protein-coding loci across the 1% of the human genome targeted by the Encyclopedia of DNA elements (ENCODE) pilot project using a combination of 5' rapid amplification of cDNA ends (RACE) and high-density resolution tiling arrays. We identified previously unannotated and often tissue- or cell-line-specific transcribed fragments (RACEfrags), both 5' distal to the annotated 5' terminus and internal to the annotated gene bounds for the vast majority (81.5%) of the tested genes. Half of the distal RACEfrags span large segments of genomic sequences away from the main portion of the coding transcript and often overlap with the upstream-annotated gene(s). Notably, at least 20% of the resultant novel transcripts have changes in their open reading frames (ORFs), most of them fusing ORFs of adjacent transcripts. A significant fraction of distal RACEfrags show expression levels comparable to those of known exons of the same locus, suggesting that they are not part of very minority splice forms. These results have significant implications concerning (1) our current understanding of the architecture of protein-coding genes; (2) our views on locations of regulatory regions in the genome; and (3) the interpretation of sequence polymorphisms mapping to regions hitherto considered to be "noncoding," ultimately relating to the identification of disease-related sequence alterations.

Genetics of normal and pathological sleep in humans.

The complexity of sleep-wake regulation, in addition to the many environmental influences, includes genetic predisposing factors, which begin to be discovered. Most of the current progress in the study of sleep genetics comes from animal models (dogs, mice, and drosophila). Multiple approaches using both animal models and different genetic techniques are needed to follow the segregation and ultimately to identify 'sleep genes' and molecular bases of sleep disorders. Recent progress in molecular genetics and the development of detailed human genome map have already led to the identification of genetic factors in several complex disorders. Only a few genes are known for which a mutation causes a sleep disorder. However, single gene disorders are rare and most common disorders are complex in terms of their genetic susceptibility, environmental factors, gene-gene, and gene-environment interactions. We review here the current progress in the genetics of normal and pathological sleep and suggest a few future perspectives.

eVOC: a controlled vocabulary for unifying gene expression data.

Expression data contribute significantly to the biological value of the sequenced human genome, providing extensive information about gene structure and the pattern of gene expression. ESTs, together with SAGE libraries and microarray experiment information, provide a broad and rich view of the transcriptome. However, it is difficult to perform large-scale expression mining of the data generated by these diverse experimental approaches. Not only is the data stored in disparate locations, but there is frequent ambiguity in the meaning of terms used to describe the source of the material used in the experiment. Untangling semantic differences between the data provided by different resources is therefore largely reliant on the domain knowledge of a human expert. We present here eVOC, a system which associates labelled target cDNAs for microarray experiments, or cDNA libraries and their associated transcripts with controlled terms in a set of hierarchical vocabularies. eVOC consists of four orthogonal controlled vocabularies suitable for describing the domains of human gene expression data including Anatomical System, Cell Type, Pathology and Developmental Stage. We have curated and annotated 7016 cDNA libraries represented in dbEST, as well as 104 SAGE libraries,with expression information,and provide this as an integrated, public resource that allows the linking of transcripts and libraries with expression terms. Both the vocabularies and the vocabulary-annotated libraries can be retrieved from http://www.sanbi.ac.za/evoc/. Several groups are involved in developing this resource with the aim of unifying transcript expression information.

Genetic basis and variable phenotypic expression of Kallmann syndrome: towards a unifying theory.

Idiopathic hypogonadotropic hypogonadism (IHH) is defined by absent or incomplete puberty and characterised biochemically by low levels of sex steroids, with low or inappropriately normal gonadotropin hormones. IHH is frequently accompanied by non-reproductive abnormalities, most commonly anosmia, which is present in 50-60% of cases and defines Kallmann syndrome. The understanding of IHH has undergone rapid evolution, both in respect of genetics and breadth of phenotype. Once considered in monogenic Mendelian terms, it is now more coherently understood as a complex genetic condition. Oligogenic and complex genetic-environmental interactions have now been identified, with physiological and environmental factors interacting in genetically susceptible individuals to alter the clinical course and phenotype. These potentially link IHH to ancient evolutionary pressures on the ancestral human genome.

Defining the genomic signature of the parous breast.

ABSTRACT: BACKGROUND: It is accepted that a woman's lifetime risk of developing breast cancer after menopause is reduced by early full term pregnancy and multiparity. This phenomenon is thought to be associated with the development and differentiation of the breast during pregnancy. METHODS: In order to understand the underlying molecular mechanisms of pregnancy induced breast cancer protection, we profiled and compared the transcriptomes of normal breast tissue biopsies from 71 parous (P) and 42 nulliparous (NP) healthy postmenopausal women using Affymetrix Human Genome U133 Plus 2.0 arrays. To validate the results, we performed real time PCR and immunohistochemistry. RESULTS: We identified 305 differentially expressed probesets (208 distinct genes). Of these, 267 probesets were up- and 38 down-regulated in parous breast samples; bioinformatics analysis using gene ontology enrichment revealed that up-regulated genes in the parous breast represented biological processes involving differentiation and development, anchoring of epithelial cells to the basement membrane, hemidesmosome and cell-substrate junction assembly, mRNA and RNA metabolic processes and RNA splicing machinery. The down-regulated genes represented biological processes that comprised cell proliferation, regulation of IGF-like growth factor receptor signaling, somatic stem cell maintenance, muscle cell differentiation and apoptosis. CONCLUSIONS: This study suggests that the differentiation of the breast imprints a genomic signature that is centered in the mRNA processing reactome. These findings indicate that pregnancy may induce a safeguard mechanism at post-transcriptional level that maintains the fidelity of the transcriptional process.

Ultra high throughput sequencing excludes MDH1 as candidate gene for RP28-linked retinitis pigmentosa.

PURPOSE: Mutations in IDH3B, an enzyme participating in the Krebs cycle, have recently been found to cause autosomal recessive retinitis pigmentosa (arRP). The MDH1 gene maps within the RP28 arRP linkage interval and encodes cytoplasmic malate dehydrogenase, an enzyme functionally related to IDH3B. As a proof of concept for candidate gene screening to be routinely performed by ultra high throughput sequencing (UHTs), we analyzed MDH1 in a patient from each of the two families described so far to show linkage between arRP and RP28. METHODS: With genomic long-range PCR, we amplified all introns and exons of the MDH1 gene (23.4 kb). PCR products were then sequenced by short-read UHTs with no further processing. Computer-based mapping of the reads and mutation detection were performed by three independent software packages. RESULTS: Despite the intrinsic complexity of human genome sequences, reads were easily mapped and analyzed, and all algorithms used provided the same results. The two patients were homozygous for all DNA variants identified in the region, which confirms previous linkage and homozygosity mapping results, but had different haplotypes, indicating genetic or allelic heterogeneity. None of the DNA changes detected could be associated with the disease. CONCLUSIONS: The MDH1 gene is not the cause of RP28-linked arRP. Our experimental strategy shows that long-range genomic PCR followed by UHTs provides an excellent system to perform a thorough screening of candidate genes for hereditary retinal degeneration.

Understanding mutational and selective processes that govern gene duplication in mammals

Abstract : Gene duplication is an essential source of material for the origin of genetic novelties. The reverse transcription of source gene mRNA followed by the genomic insertion of the resulting cDNA - retroposition - has provided the human genome with at least ~3600 detectable retrocopies. We find that ~30% of these retrocopies are transcribed, generally in testes. Their transcription often relies on preexisting regulatory elements (or open chromatin) close to their insertion site, which is illustrated by mRNA molecules containing retrocopies fused to their neighboring genes. Retrocopies appear to have been profoundly shaped by selection. Consistently, human retrocopies with an intact open reading (ORF) are more often transcribed than retropseudogenes, which leads to a minimal estimate of 120 functional retrogenes present in our genome. We also performed an analysis of Ka/Ks for human retrocopies. This analysis demonstrates that several intact retrocopies evolved under purifying selection and yields an estimated formation rate of ~1 retrogene per million year in the primate lineage. Using DNA sequencing and evolutionary simulations, we have identified 7 such primate-specific retrogenes that emerged on the lineage leading to humans In therian genomes, we found an excess of retrogenes with X-linked parents. Expression analyses support the idea that this "out of X" movement was driven by natural selection to produce autosomal functional counterparts for X-linked genes, which are silenced during male meiosis. Phylogenetic dating of this "out of X" movement suggests that our sex chromosomes arose about 180 MYA ago and are thus much younger than previously thought. Finally, we have also analyzed young gene duplications (and deletions) that arose by non allelic-homologous recombination and are not fixed in species. Using wild-caught and laboratory animals, we detected thousands of DNA segments that are polymorphic in copy number in mice. These copy number variants were found to profoundly alter the transcriptome of several mouse tissues. Strikingly, their influence on gene expression is not limited to the gene they contain but seems to extend to genes located up to 1.5 million bases away.

KRAB-zinc finger proteins and KAP1 can mediate long-range transcriptional repression through heterochromatin spreading.

Krüppel-associated box domain-zinc finger proteins (KRAB-ZFPs) are tetrapod-specific transcriptional repressors encoded in the hundreds by the human genome. In order to explore their as yet ill-defined impact on gene expression, we developed an ectopic repressor assay, allowing the study of KRAB-mediated transcriptional regulation at hundreds of different transcriptional units. By targeting a drug-controllable KRAB-containing repressor to gene-trapping lentiviral vectors, we demonstrate that KRAB and its corepressor KAP1 can silence promoters located several tens of kilobases (kb) away from their DNA binding sites, with an efficiency which is generally higher for promoters located within 15 kb or less. Silenced promoters exhibit a loss of histone H3-acetylation, an increase in H3 lysine 9 trimethylation (H3K9me3), and a drop in RNA Pol II recruitment, consistent with a block of transcriptional initiation following the establishment of silencing marks. Furthermore, we reveal that KRAB-mediated repression is established by the long-range spreading of H3K9me3 and heterochromatin protein 1 beta (HP1beta) between the repressor binding site and the promoter. We confirm the biological relevance of this phenomenon by documenting KAP1-dependent transcriptional repression at an endogenous KRAB-ZFP gene cluster, where KAP1 binds to the 3' end of genes and mediates propagation of H3K9me3 and HP1beta towards their 5' end. Together, our data support a model in which KRAB/KAP1 recruitment induces long-range repression through the spread of heterochromatin. This finding not only suggests auto-regulatory mechanisms in the control of KRAB-ZFP gene clusters, but also provides important cues for interpreting future genome-wide DNA binding data of KRAB-ZFPs and KAP1.

A genomewide scan for early-onset coronary artery disease in 438 families: the GENECARD Study.

A family history of coronary artery disease (CAD), especially when the disease occurs at a young age, is a potent risk factor for CAD. DNA collection in families in which two or more siblings are affected at an early age allows identification of genetic factors for CAD by linkage analysis. We performed a genomewide scan in 1,168 individuals from 438 families, including 493 affected sibling pairs with documented onset of CAD before 51 years of age in men and before 56 years of age in women. We prospectively defined three phenotypic subsets of families: (1) acute coronary syndrome in two or more siblings; (2) absence of type 2 diabetes in all affected siblings; and (3) atherogenic dyslipidemia in any one sibling. Genotypes were analyzed for 395 microsatellite markers. Regions were defined as providing evidence for linkage if they provided parametric two-point LOD scores >1.5, together with nonparametric multipoint LOD scores >1.0. Regions on chromosomes 3q13 (multipoint LOD = 3.3; empirical P value <.001) and 5q31 (multipoint LOD = 1.4; empirical P value <.081) met these criteria in the entire data set, and regions on chromosomes 1q25, 3q13, 7p14, and 19p13 met these criteria in one or more of the subsets. Two regions, 3q13 and 1q25, met the criteria for genomewide significance. We have identified a region on chromosome 3q13 that is linked to early-onset CAD, as well as additional regions of interest that will require further analysis. These data provide initial areas of the human genome where further investigation may reveal susceptibility genes for early-onset CAD.