Interrogating eleven fast-evolving genes for signatures of recent positive selection in worldwide human populations

Different signatures of natural selection persist over varying time scales in our genome, revealing possible episodes of adaptative evolution during human history. Here, we identify genes showing signatures of ancestral positive selection in the human lineage and investigate whether some of those genes have been evolving adaptatively in extant human populations. Specifically, we compared more than 11,000 human genes with their orthologs inchimpanzee, mouse, rat and dog and applied a branch-site likelihood method to test for positive selection on the human lineage. Among the significant cases, a robust set of 11 genes were then further explored for signatures of recent positive selection using SNP data. We genotyped 223 SNPs in 39 worldwide populations from the HGDP Diversity panel and supplemented this information with available genotypes for up to 4,814 SNPs distributed along 2 Mb centered on each gene. After exploring the allele frequency spectrum, population differentiation and the maintainance of long unbroken haplotypes, we found signals of recent adaptative phenomena in only one of the 11 candidate gene regions. However, the signal ofrecent selection in this region may come from a different, neighbouring gene (CD5) ratherthan from the candidate gene itself (VPS37C). For this set of positively-selected genes in thehuman lineage, we find no indication that these genes maintained their rapid evolutionarypace among human populations. Based on these data, it therefore appears that adaptation forhuman-specific and for population-specific traits may have involved different genes.

Functional diversity of bacterial genes associated with aromatic hydrocarbon degradation in anthropogenic dark earth of Amazonia

The objective of this work was to evaluate the catabolic gene diversity for the bacterial degradation of aromatic hydrocarbons in anthropogenic dark earth of Amazonia (ADE) and their biochar (BC). Functional diversity analyses in ADE soils can provide information on how adaptive microorganisms may influence the fertility of soils and what is their involvement in biogeochemical cycles. For this, clone libraries containing the gene encoding for the alpha subunit of aromatic ring-hydroxylating dioxygenases (α-ARHD bacterial gene) were constructed, totaling 800 clones. These libraries were prepared from samples of an ADE soil under two different land uses, located at the Caldeirão Experimental Station - secondary forest (SF) and agriculture (AG) -, and the biochar (SF_BC and AG_BC, respectively). Heterogeneity estimates indicated greater diversity in BC libraries; and Venn diagrams showed more unique operational protein clusters (OPC) in the SF_BC library than the ADE soil, which indicates that specific metabolic processes may occur in biochar. Phylogenetic analysis showed unidentified dioxygenases in ADE soils. Libraries containing functional gene encoding for the alpha subunit of the aromatic ring-hydroxylating dioxygenases (ARHD) gene from biochar show higher diversity indices than those of ADE under secondary forest and agriculture.

dKDM5/LID regulates H3K4me3 dynamics at the transcription-start site (TSS) of actively transcribed developmental genes

H3K4me3 is a histone modification that accumulates at the transcription-start site (TSS) of active genes and is known to be important for transcription activation. The way in which H3K4me3 is regulated at TSS and the actual molecular basis of its contribution to transcription remain largely unanswered. To address these questions, we have analyzed the contribution of dKDM5/LID, the main H3K4me3 demethylase in Drosophila, to the regulation of the pattern of H3K4me3. ChIP-seq results show that, at developmental genes, dKDM5/LID localizes at TSS and regulates H3K4me3. dKDM5/LID target genes are highly transcribed and enriched in active RNApol II and H3K36me3, suggesting a positive contribution to transcription. Expression-profiling show that, though weakly, dKDM5/LID target genes are significantly downregulated upon dKDM5/LID depletion. Furthermore, dKDM5/LID depletion results in decreased RNApol II occupancy, particularly by the promoter-proximal Pol lloser5 form. Our results also show that ASH2, an evolutionarily conserved factor that locates at TSS and is required for H3K4me3, binds and positively regulates dKDM5/LID target genes. However, dKDM5/LID and ASH2 do not bind simultaneously and recognize different chromatin states, enriched in H3K4me3 and not, respectively. These results indicate that, at developmental genes, dKDM5/LID and ASH2 coordinately regulate H3K4me3 at TSS and that this dynamic regulation contributes to transcription.

Mutations in the EXT1 and EXT2 genes in Spanish patients with multiple osteochondromas

Multiple osteochondromas is an autosomal dominant skeletal disorder characterized by the formation of multiple cartilage-capped tumours. Two causal genes have been identified, EXT1 and EXT2, which account for 65% and 30% of cases, respectively. We have undertaken a mutation analysis of the EXT1 and EXT2 genes in 39 unrelated Spanish patients, most of them with moderate phenotype, and looked for genotype-phenotype correlations. We found the mutant allele in 37 patients, 29 in EXT1 and 8 in EXT2. Five of the EXT1 mutations were deletions identified by MLPA. Two cases of mosaicism were documented. We detected a lower number of exostoses in patients with missense mutation versus other kinds of mutations. In conclusion, we found a mutation in EXT1 or in EXT2 in 95% of the Spanish patients. Eighteen of the mutations were novel.

Mutations in the EXT1 and EXT2 genes in Spanish patients with multiple osteochondromas

Multiple osteochondromas is an autosomal dominant skeletal disorder characterized by the formation of multiple cartilage-capped tumours. Two causal genes have been identified, EXT1 and EXT2, which account for 65% and 30% of cases, respectively. We have undertaken a mutation analysis of the EXT1 and EXT2 genes in 39 unrelated Spanish patients, most of them with moderate phenotype, and looked for genotype-phenotype correlations. We found the mutant allele in 37 patients, 29 in EXT1 and 8 in EXT2. Five of the EXT1 mutations were deletions identified by MLPA. Two cases of mosaicism were documented. We detected a lower number of exostoses in patients with missense mutation versus other kinds of mutations. In conclusion, we found a mutation in EXT1 or in EXT2 in 95% of the Spanish patients. Eighteen of the mutations were novel.

Transcriptional regulatory patterns of the myelin basic protein and malic enzyme genes by the thyroid hormone receptors alpha1 and beta1.

While there is evidence that the two ubiquitously expressed thyroid hormone (T3) receptors, TRalpha1 and TRbeta1, have distinct functional specificities, the mechanism by which they discriminate potential target genes remains largely unexplained. In this study, we demonstrate that the thyroid hormone response elements (TRE) from the malic enzyme and myelin basic protein genes (METRE and MBPTRE) respectively, are not functionally equivalent. The METRE, which is a direct repeat motif with a 4-base pair gap between the two half-site hexamers binds thyroid hormone receptor as a heterodimer with 9-cis-retinoic acid receptor (RXR) and mediates a high T3-dependent activation in response to TRalpha1 or TRbeta1 in NIH3T3 cells. In contrast, the MBPTRE, which consists of an inverted palindrome formed by two hexamers spaced by 6 base pairs, confers an efficient transactivation by TRbeta1 but a poor transactivation by TRalpha1. While both receptors form heterodimers with RXR on MBPTRE, the poor transactivation by TRalpha1 correlates also with its ability to bind efficiently as a monomer. This monomer, which is only observed with TRalpha1 bound to MBPTRE, interacts neither with N-CoR nor with SRC-1, explaining its functional inefficacy. However, in Xenopus oocytes, in which RXR proteins are not detectable, the transactivation mediated by TRalpha1 and TRbeta1 is equivalent and independent of a RXR supply, raising the question of the identity of the thyroid hormone receptor partner in these cells. Thus, in mammalian cells, the binding characteristics of TRalpha1 to MBPTRE (i.e. high monomer binding efficiency and low transactivation activity) might explain the particular pattern of T3 responsiveness of MBP gene expression during central nervous system development.

Identification of cancer/testis-antigen genes by massively parallel signature sequencing.

Massively parallel signature sequencing (MPSS) generates millions of short sequence tags corresponding to transcripts from a single RNA preparation. Most MPSS tags can be unambiguously assigned to genes, thereby generating a comprehensive expression profile of the tissue of origin. From the comparison of MPSS data from 32 normal human tissues, we identified 1,056 genes that are predominantly expressed in the testis. Further evaluation by using MPSS tags from cancer cell lines and EST data from a wide variety of tumors identified 202 of these genes as candidates for encoding cancer/testis (CT) antigens. Of these genes, the expression in normal tissues was assessed by RT-PCR in a subset of 166 intron-containing genes, and those with confirmed testis-predominant expression were further evaluated for their expression in 21 cancer cell lines. Thus, 20 CT or CT-like genes were identified, with several exhibiting expression in five or more of the cancer cell lines examined. One of these genes is a member of a CT gene family that we designated as CT45. The CT45 family comprises six highly similar (>98% cDNA identity) genes that are clustered in tandem within a 125-kb region on Xq26.3. CT45 was found to be frequently expressed in both cancer cell lines and lung cancer specimens. Thus, MPSS analysis has resulted in a significant extension of our knowledge of CT antigens, leading to the discovery of a distinctive X-linked CT-antigen gene family.

Epigenetic regulatory elements associate with specific histone modifications to prevent silencing of telomeric genes.

In eukaryotic cells, transgene expression levels may be limited by an unfavourable chromatin structure at the integration site. Epigenetic regulators are DNA sequences which may protect transgenes from such position effect. We evaluated different epigenetic regulators for their ability to protect transgene expression at telomeres, which are commonly associated to low or inconsistent expression because of their repressive chromatin environment. Although to variable extents, matrix attachment regions (MARs), ubiquitous chromatin opening element (UCOE) and the chicken cHS4 insulator acted as barrier elements, protecting a telomeric-distal transgene from silencing. MARs also increased the probability of silent gene reactivation in time-course experiments. Additionally, all MARs improved the level of expression in non-silenced cells, unlike other elements. MARs were associated to histone marks usually linked to actively expressed genes, especially acetylation of histone H3 and H4, suggesting that they may prevent the spread of silencing chromatin by imposing acetylation marks on nearby nucleosomes. Alternatively, an UCOE was found to act by preventing deposition of repressive chromatin marks. We conclude that epigenetic DNA elements used to enhance and stabilize transgene expression all have specific epigenetic signature that might be at the basis of their mode of action.

Revisit of multiple epiphyseal dysplasia: ethnic difference in genotypes and comparison of radiographic features linked to the COMP and MATN3 genes.

Multiple epiphyseal dysplasia (MED) is a genetically heterogeneous group of diseases characterized by variable degrees of epiphyseal abnormality primarily involving the hip and knee joints. The purpose of this study was to investigate the frequency of mutations in individuals with a clinical and radiographic diagnosis of MED and to test the hypothesis that characteristic radiological findings may be helpful in predicting the gene responsible. The radiographs of 74 Korean patients were evaluated by a panel of skeletal dysplasia experts. Six genes known to be associated with MED (COMP, MATN3, COL9A1, COL9A2, COL9A3, and DTDST) were screened by sequencing. Mutations were found in 55 of the 63 patients (87%). MATN3 mutations were found in 30 patients (55%), followed by COMP mutations in 23 (41%), and COL9A2 and DTDST mutations in one patient (2%) each. Comparisons of radiographic findings in patients with COMP and MATN3 mutations showed that albeit marked abnormalities in hip and knee joints were observed in both groups, the degree of involvement and the morphology of dysplastic epiphyses differed markedly. The contour of the pelvic acetabulum, the presence of metaphyseal vertical striations, and/or the brachydactyly of the hand were also found to be highly correlated with the genotypes. The study confirms that MATN3 and COMP are the genes most frequently responsible for MED and that subtle radiographic signs may give precious indications on which gene(s) should be prioritized for mutational screening in a given individual.

Expressão do mRNA de genes mitocondriais e desempenho produtivo de codornas alimentadas com glicerol

O objetivo deste trabalho foi avaliar o efeito de dietas com glicerol no desempenho produtivo de codornas japonesas de corte e na expressão do mRNA de genes mitocondriais da proteína adenina nucleotídeo translocase (ANT) e da proteína desacopladora (UCP), envolvidas no metabolismo energético e na resposta ao estresse oxidativo. As codornas foram alimentadas com dietas contendo 0, 8 e 12% de glicerol, em substituição parcial ao milho. Aos 28 dias de idade, o RNA total foi extraído de amostras do músculo do peito e a síntese do cDNA foi feita por meio de qRT‑PCR com iniciadores específicos para genes da ANT e UCP, obtidos de Gallus gallus. A conversão alimentar e o consumo de ração foram avaliados para as três dietas testadas. A adição de 8% de glicerol não afetou o desempenho dos animais. No entanto, a adição de 12% aumentou o consumo de ração e piorou a conversão alimentar. O ganho de peso não foi afetado pela inclusão de glicerol na dieta. No grupo alimentado com 8% de glicerol, a expressão da UCP aumentou, mas a da ANT não variou, em comparação ao controle. A expressão da UCP foi menor e a da ANT foi maior no grupo alimentado com 12% de glicerol. A inclusão de 8% de glicerol na dieta não afeta o desempenho de codornas de corte, embora aumente a expressão da UCP. A inclusão de 12% de glicerol piora o desempenho e aumenta a expressão da ANT.

Etude des gènes codant les protéines de gouttelettes lipidiques chez les patients avec la stéatose hépatique[Study of genes encoding proteins of lipid droplets in patients with fatty liver]

Introduction: La prévalence de la «non-alcoholic fatty liver disease (NAFLD)» dans les pays industrialisés augment de manière exponentielle. La NAFLD se développe d'une simple stéatose hépatique jusqu'à l'hépatite, puis à la cirrhose. De plus, la stéatose hépatique est fréquemment accompagnée par une résistance à l'insuline, une des causes principales du diabète. Les lipides intermédiaires, tels que céramides et diacylglycérols, ont été décrits comme induisant la résistance à l'insuline. Cependant, nous avons démontré dans notre modèle de stéatose hépatique, que les souris présentant une invalidation de la protéine «microsomal triglyceride transfer protein» (Mtpp) au niveau hépatique, ne développent pas de résistance à l'insuline. Ceci suggère fortement l'existence d'autres mécanismes susceptibles d'induire la résistance à l'insuline. Résultats: Grâce à une analyse de Microarray, nous avons observé une augmentation de l'expression des gènes «cell-death inducing DFFA-like effector c (CIDEC)», «lipid storage droplet protein 5 (LSDP5)» et «Bernardinelli-Seip congenital lipodystrophy 2 homolog (Seipin)» dans le foie des souris Mttp. Ces gènes ont récemment été identifiés comme des protéines localisées autour des gouttelettes lipidiques. Nous avons également constaté que la souris Mttp développe plutôt une microstéatose (petites gouttelettes lipidiques) qu'une macrostéatose qui est normalement observée chez les patients avec NAFLD. Nous avons étudié l'expression des gènes associés aux gouttelettes lipidiques chez les patients obèses avec stéatose hépatique, avec ou sans résistance à l'insuline. Comparés aux sujets sains sans stéatose hépatique, les patients avec la stéatose ont une expression significativement plus élevée. De manière intéressante, les patients avec résistance à l'insuline ont une diminution de ces expressions. Conclusion : Ces données suggèrent que les gènes des gouttelettes lipidiques sont impliqués dans le développement de la stéatose hépatique chez l'homme et peut-être contribue à la mise en place de la résistance à l'insuline.

Mutational screening of splicing factor genes in cases with autosomal dominant retinitis pigmentosa.

PURPOSE: Mutations in genes encoding proteins from the tri-snRNP complex of the spliceosome account for more than 12% of cases of autosomal dominant retinitis pigmentosa (adRP). Although the exact mechanism by which splicing factor defects trigger photoreceptor death is not completely clear, their role in retinitis pigmentosa has been demonstrated by several genetic and functional studies. To test for possible novel associations between splicing factors and adRP, we screened four tri-snRNP splicing factor genes (EFTUD2, PRPF4, NHP2L1, and AAR2) as candidate disease genes. METHODS: We screened up to 303 patients with adRP from Europe and North America who did not carry known RP mutations. Exon-PCR and Sanger methods were used to sequence the NHP2L1 and AAR2 genes, while the sequences of EFTUD2 and PRPF4 were obtained by using long-range PCRs spanning coding and non-coding regions followed by next-generation sequencing. RESULTS: We detected novel missense changes in individual patients in the sequence of the genes PRPF4 and EFTUD2, but the role of these changes in relationship to disease could not be verified. In one other patient we identified a novel nucleotide substitution in the 5' untranslated region (UTR) of NHP2L1, which did not segregate with the disease in the family. CONCLUSIONS: The absence of clearly pathogenic mutations in the candidate genes screened in our cohort suggests that EFTUD2, PRPF4, NHP2L1, and AAR2 are either not involved in adRP or are associated with the disease in rare instances, at least as observed in this study in patients of European and North American origin.

Prospective study of insulin-like growth factor-I, insulin-like growth factor-binding protein 3, genetic variants in the IGF1 and IGFBP3 genes and risk of coronary artery disease.

Although experimental studies have suggested that insulin-like growth factor I (IGF-I) and its binding protein IGFBP-3 might have a role in the aetiology of coronary artery disease (CAD), the relevance of circulating IGFs and their binding proteins in the development of CAD in human populations is unclear. We conducted a nested case-control study, with a mean follow-up of six years, within the EPIC-Norfolk cohort to assess the association between circulating levels of IGF-I and IGFBP-3 and risk of CAD in up to 1,013 cases and 2,055 controls matched for age, sex and study enrolment date. After adjustment for cardiovascular risk factors, we found no association between circulating levels of IGF-I or IGFBP-3 and risk of CAD (odds ratio: 0.98 (95% Cl 0.90-1.06) per 1 SD increase in circulating IGF-I; odds ratio: 1.02 (95% Cl 0.94-1.12) for IGFBP-3). We examined associations between tagging single nucleotide polymorphisms (tSNPs) at the IGF1 and IGFBP3 loci and circulating IGF-I and IGFBP-3 levels in up to 1,133 cases and 2,223 controls and identified three tSNPs (rs1520220, rs3730204, rs2132571) that showed independent association with either circulating IGF-I or IGFBP-3 levels. In an assessment of 31 SNPs spanning the IGF1 or IGFBP3 loci, none were associated with risk of CAD in a meta-analysis that included EPIC-Norfolk and eight additional studies comprising up to 9,319 cases and 19,964 controls. Our results indicate that IGF-I and IGFBP-3 are unlikely to be importantly involved in the aetiology of CAD in human populations.

Immunopathogenesis of Asthma and Atopic Diseases – The specific Role of a selected Panel of Genes in human T helper Cell Differentiation

T helper cell (Th) functions are crucial for proper immune defence against various intra- and extracellular pathogens. According to the specific immune responses, Th cells can be classified into subtypes, Th1 and Th2 cells being the most frequently characterized classes. Th1 and Th2 cells interact with other immune cells by regulating their functions with specific cytokine production. IFN, IL-2 and TNF- are the cytokines predominantly produced by Th1 cells whereas Th2 cells produce Th2-type cytokines, such as IL-4, IL-5 and IL-13. Upon TCR activation and in the presence of polarizing cytokines, Th cells differentiate into effector subtypes from a common precursor cell. IFN and IL-12 are the predominant Th1 polarizing cytokines whereas IL-4 directs Th2 polarization. The cytokines mediate their effects through specific receptor signalling. The differentiation process is complex, involving various signalling molecules and routes, as well as functions of the specific transcription factors. The functions of the Th1/Th2 cells are tightly regulated; however, knowledge on human Th cell differentiation is, as yet, fairly poor. The susceptibility for many immune-mediated disorders often originates from disturbed Th cell responses. Thus, research is needed for defining the molecular mechanisms involved in the differentiation and balanced functions of the Th cells. Importantly, the new information obtained will be crucial for a better understanding of the pathogenesis of immune-mediated disorders, such as asthma or autoimmune diseases. In the first subproject of this thesis, the role of genetic polymorphisms in the human STAT6, GATA3 and STAT4 genes were investigated for asthma or atopy susceptibility in Finnish asthma families by association analysis. These genes code for key transcription factors regulating Th cell differentiation. The study resulted in the identification of a GATA3 haplotype that associated with asthma and related traits (high serum IgE level). In the second subproject, an optimized method for human primary T cell transfection and enrichment was established. The method can be utilized for functional studies for the selected genes of interest. The method was also utilized in the third subproject, which aimed at the identification of novel genes involved in early human Th cell polarization (0-48h) using genome-wide oligonucleotide arrays. As a result, numerous genes and ESTs with known or unknown functions were identified in the study. Using an shRNA knockdown approach, a panel of novel IL-4/STAT6 regulated genes were identified in the functional studies of the genes. Moreover, one of the genes, NDFIP2, with a previously uncharacterized role in the human Th differentiation, was observed to promote IFN production of the differentiated Th1 cells. Taken together, the results obtained have revealed potential new relevant candidate genes serving as a basis for further studies characterizing the detailed networks involved in the human Th cell differentiation as well as in the genetic susceptibility of Th-mediated immune disorders.