Sex-stratified genome-wide association studies including 270,000 individuals show sexual dimorphism in genetic loci for anthropometric traits.

Given the anthropometric differences between men and women and previous evidence of sex-difference in genetic effects, we conducted a genome-wide search for sexually dimorphic associations with height, weight, body mass index, waist circumference, hip circumference, and waist-to-hip-ratio (133,723 individuals) and took forward 348 SNPs into follow-up (additional 137,052 individuals) in a total of 94 studies. Seven loci displayed significant sex-difference (FDR<5%), including four previously established (near GRB14/COBLL1, LYPLAL1/SLC30A10, VEGFA, ADAMTS9) and three novel anthropometric trait loci (near MAP3K1, HSD17B4, PPARG), all of which were genome-wide significant in women (P<5×10(-8)), but not in men. Sex-differences were apparent only for waist phenotypes, not for height, weight, BMI, or hip circumference. Moreover, we found no evidence for genetic effects with opposite directions in men versus women. The PPARG locus is of specific interest due to its role in diabetes genetics and therapy. Our results demonstrate the value of sex-specific GWAS to unravel the sexually dimorphic genetic underpinning of complex traits.

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.

Exhaustion of tumor-specific CD8⁺ T cells in metastases from melanoma patients.

In chronic viral infections, CD8⁺ T cells become functionally deficient and display multiple molecular alterations. In contrast, only little is known of self- and tumor-specific CD8⁺ T cells from mice and humans. Here we determined molecular profiles of tumor-specific CD8⁺ T cells from melanoma patients. In peripheral blood from patients vaccinated with CpG and the melanoma antigen Melan-A/MART-1 peptide, we found functional effector T cell populations, with only small but nevertheless significant differences in T cells specific for persistent herpesviruses (EBV and CMV). In contrast, Melan-A/MART-1-specific T cells isolated from metastases from patients with melanoma expressed a large variety of genes associated with T cell exhaustion. The identified exhaustion profile revealed extended molecular alterations. Our data demonstrate a remarkable coexistence of effector cells in circulation and exhausted cells in the tumor environment. Functional T cell impairment is mediated by inhibitory receptors and further molecular pathways, which represent potential targets for cancer therapy.

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.

Genome-wide meta-analysis for serum calcium identifies significantly associated SNPs near the calcium-sensing receptor (CASR) gene.

Calcium has a pivotal role in biological functions, and serum calcium levels have been associated with numerous disorders of bone and mineral metabolism, as well as with cardiovascular mortality. Here we report results from a genome-wide association study of serum calcium, integrating data from four independent cohorts including a total of 12,865 individuals of European and Indian Asian descent. Our meta-analysis shows that serum calcium is associated with SNPs in or near the calcium-sensing receptor (CASR) gene on 3q13. The top hit with a p-value of 6.3 x 10(-37) is rs1801725, a missense variant, explaining 1.26% of the variance in serum calcium. This SNP had the strongest association in individuals of European descent, while for individuals of Indian Asian descent the top hit was rs17251221 (p = 1.1 x 10(-21)), a SNP in strong linkage disequilibrium with rs1801725. The strongest locus in CASR was shown to replicate in an independent Icelandic cohort of 4,126 individuals (p = 1.02 x 10(-4)). This genome-wide meta-analysis shows that common CASR variants modulate serum calcium levels in the adult general population, which confirms previous results in some candidate gene studies of the CASR locus. This study highlights the key role of CASR in calcium regulation.

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.

Characterization of alternatively spliced products and tissue-specific isoforms of USP28 and USP25

Background: The ubiquitin-dependent protein degradation pathway is essential for the proteolysis of intracellular proteins and peptides. Deubiquitinating enzymes constitute a complex protein family involved in a multitude of cellular processes. The ubiquitin-specific proteases (UBP) are a group of enzymes whose predicted function is to reverse the ubiquitinating reaction by removing ubiquitin from a large variety of substrates. We have lately reported the characterization of human USP25, a specific-ubiquitin protease gene at 21q11.2, with a specific pattern of expression in murine fetal brains and adult testis. Results: Database homology searches at the DNA and protein levels and cDNA library screenings led to the identification of a new UBP member in the human genome, named USP28, at 11q23. This novel gene showed preferential expression in heart and muscle. Moreover, cDNA, expressed sequence tag and RT-PCR analyses provided evidence for alternatively spliced products and tissue-specific isoforms. Concerning function, USP25 overexpression in Down syndrome fetal brains was shown by real-time PCR. Conclusions: On the basis of the genomic and protein sequence as well as the functional data, USP28 and USP25 establish a new subfamily of deubiquitinating enzymes. Both genes have alternatively spliced exons that could generate protein isoforms with distinct tissue-specific activity. The overexpression of USP25 in Down syndrome fetal brains supports the gene-dosage effects suggested for other UBP members related to aneuploidy syndromes.

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.

The GENCODE pseudogene resource.

BACKGROUND: Pseudogenes have long been considered as nonfunctional genomic sequences. However, recent evidence suggests that many of them might have some form of biological activity, and the possibility of functionality has increased interest in their accurate annotation and integration with functional genomics data. RESULTS: As part of the GENCODE annotation of the human genome, we present the first genome-wide pseudogene assignment for protein-coding genes, based on both large-scale manual annotation and in silico pipelines. A key aspect of this coupled approach is that it allows us to identify pseudogenes in an unbiased fashion as well as untangle complex events through manual evaluation. We integrate the pseudogene annotations with the extensive ENCODE functional genomics information. In particular, we determine the expression level, transcription-factor and RNA polymerase II binding, and chromatin marks associated with each pseudogene. Based on their distribution, we develop simple statistical models for each type of activity, which we validate with large-scale RT-PCR-Seq experiments. Finally, we compare our pseudogenes with conservation and variation data from primate alignments and the 1000 Genomes project, producing lists of pseudogenes potentially under selection. CONCLUSIONS: At one extreme, some pseudogenes possess conventional characteristics of functionality; these may represent genes that have recently died. On the other hand, we find interesting patterns of partial activity, which may suggest that dead genes are being resurrected as functioning non-coding RNAs. The activity data of each pseudogene are stored in an associated resource, psiDR, which will be useful for the initial identification of potentially functional pseudogenes.

Evolution of the P-type II ATPase gene family in the fungi and presence of structural genomic changes among isolates of Glomus intraradices.

BACKGROUND: The P-type II ATPase gene family encodes proteins with an important role in adaptation of the cell to variation in external K+, Ca2+ and Na2+ concentrations. The presence of P-type II gene subfamilies that are specific for certain kingdoms has been reported but was sometimes contradicted by discovery of previously unknown homologous sequences in newly sequenced genomes. Members of this gene family have been sampled in all of the fungal phyla except the arbuscular mycorrhizal fungi (AMF; phylum Glomeromycota), which are known to play a key-role in terrestrial ecosystems and to be genetically highly variable within populations. Here we used highly degenerate primers on AMF genomic DNA to increase the sampling of fungal P-Type II ATPases and to test previous predictions about their evolution. In parallel, homologous sequences of the P-type II ATPases have been used to determine the nature and amount of polymorphism that is present at these loci among isolates of Glomus intraradices harvested from the same field. RESULTS: In this study, four P-type II ATPase sub-families have been isolated from three AMF species. We show that, contrary to previous predictions, P-type IIC ATPases are present in all basal fungal taxa. Additionally, P-Type IIE ATPases should no longer be considered as exclusive to the Ascomycota and the Basidiomycota, since we also demonstrate their presence in the Zygomycota. Finally, a comparison of homologous sequences encoding P-type IID ATPases showed unexpectedly that indel mutations among coding regions, as well as specific gene duplications occur among AMF individuals within the same field. CONCLUSION: On the basis of these results we suggest that the diversification of P-Type IIC and E ATPases followed the diversification of the extant fungal phyla with independent events of gene gains and losses. Consistent with recent findings on the human genome, but at a much smaller geographic scale, we provided evidence that structural genomic changes, such as exonic indel mutations and gene duplications are less rare than previously thought and that these also occur within fungal populations.

MGMT methylation status: the advent of stratified therapy in glioblastoma?

Glioblastomas are the most malignant gliomas with median survival times of only 15 months despite modern therapies. All standard treatments are palliative. Pathogenetic factors are diverse, hence, stratified treatment plans are warranted considering the molecular heterogeneity among these tumors. However, most patients are treated with "one fits all" standard therapies, many of them with minor response and major toxicities. The integration of clinical and molecular information, now becoming available using new tools such as gene arrays, proteomics, and molecular imaging, will take us to an era where more targeted and effective treatments may be implemented. A first step towards the design of such therapies is the identification of relevant molecular mechanisms driving the aggressive biological behavior of glioblastoma. The accumulation of diverse aberrations in regulatory processes enables tumor cells to bypass the effects of most classical therapies available. Molecular alterations underlying such mechanisms comprise aberrations on the genetic level, such as point mutations of distinct genes, or amplifications and deletions, while others result from epigenetic modifications such as aberrant methylation of CpG islands in the regulatory sequence of genes. Epigenetic silencing of the MGMT gene encoding a DNA repair enzyme was recently found to be of predictive value in a randomized clinical trial for newly diagnosed glioblastoma testing the addition of the alkylating agent temozolomide to standard radiotherapy. Determination of the methylation status of the MGMT promoter may become the first molecular diagnostic tool to identify patients most likely to respond that will allow individually tailored therapy in glioblastoma. To date, the test for the MGMT-methylation status is the only tool available that may direct the choice for alkylating agents in glioblastoma patients, but many others may hopefully become part of an arsenal to stratify patients to respective targeted therapies within the next years.

A new locus for congenital cataract, microcornea, microphthalmia, and atypical iris coloboma maps to chromosome 2.

OBJECTIVE: To report a novel phenotype of autosomal dominant atypical congenital cataract associated with variable expression of microcornea, microphthalmia, and iris coloboma linked to chromosome 2. Molecular analysis of this phenotype may improve our understanding of anterior segment development. DESIGN: Observational case study, genome linkage analysis, and gene mutation screening. PARTICIPANTS: Three families, 1 Egyptian and 2 Belgians, with a total of 31 affected were studied. METHODS: Twenty-one affected subjects and 9 first-degree relatives underwent complete ophthalmic examination. In the Egyptian family, exclusion of PAX6, CRYAA, and MAF genes was demonstrated by haplotype analysis using microsatellite markers on chromosomes 11, 16, and 21. Genome-wide linkage analysis was then performed using 385 microsatellite markers on this family. In the 2 Belgian families, the PAX6 gene was screened for mutations by direct sequencing of all exons. MAIN OUTCOME MEASURES: Phenotype description, genome-wide linkage of the phenotype, linkage to the PAX6, CRYAA, and MAF genes, and mutation detection in the PAX6 gene. RESULTS: Affected members of the 3 families had bilateral congenital cataracts inherited in an autosomal dominant pattern. A novel form of hexagonal nuclear cataract with cortical riders was expressed. Among affected subjects with available data, 95% had microcornea, 39% had microphthalmia, and 38% had iris coloboma. Seventy-five percent of the colobomata were atypical, showing a nasal superior location in 56%. A positive lod score of 4.86 was obtained at theta = 0 for D2S2309 on chromosome 2, a 4.9-Mb common haplotype flanked by D2S2309 and D2S2358 was obtained in the Egyptian family, and linkage to the PAX6, CRYAA, or MAF gene was excluded. In the 2 Belgian families, sequencing of the junctions and all coding exons of PAX6 did not reveal any molecular change. CONCLUSIONS: We describe a novel phenotype that includes the combination of a novel form of congenital hexagonal cataract, with variably expressed microcornea, microphthalmia, and atypical iris coloboma, not caused by PAX6 and mapping to chromosome 2. FINANCIAL DISCLOSURE(S): The authors have no proprietary or commercial interest in any materials discussed in this article.

Networking of differentially expressed genes in human cancer cells resistant to methotrexate

BACKGROUND: The need for an integrated view of data obtained from high-throughput technologies gave rise to network analyses. These are especially useful to rationalize how external perturbations propagate through the expression of genes. To address this issue in the case of drug resistance, we constructed biological association networks of genes differentially expressed in cell lines resistant to methotrexate (MTX). METHODS: Seven cell lines representative of different types of cancer, including colon cancer (HT29 and Caco2), breast cancer (MCF-7 and MDA-MB-468), pancreatic cancer (MIA PaCa-2), erythroblastic leukemia (K562) and osteosarcoma (Saos-2), were used. The differential expression pattern between sensitive and MTX-resistant cells was determined by whole human genome microarrays and analyzed with the GeneSpring GX software package. Genes deregulated in common between the different cancer cell lines served to generate biological association networks using the Pathway Architect software. RESULTS: Dikkopf homolog-1 (DKK1) is a highly interconnected node in the network generated with genes in common between the two colon cancer cell lines, and functional validations of this target using small interfering RNAs (siRNAs) showed a chemosensitization toward MTX. Members of the UDP-glucuronosyltransferase 1A (UGT1A) family formed a network of genes differentially expressed in the two breast cancer cell lines. siRNA treatment against UGT1A also showed an increase in MTX sensitivity. Eukaryotic translation elongation factor 1 alpha 1 (EEF1A1) was overexpressed among the pancreatic cancer, leukemia and osteosarcoma cell lines, and siRNA treatment against EEF1A1 produced a chemosensitization toward MTX. CONCLUSIONS: Biological association networks identified DKK1, UGT1As and EEF1A1 as important gene nodes in MTX-resistance. Treatments using siRNA technology against these three genes showed chemosensitization toward MTX.