Computational and statistical analysis of copy number variants in normal and cancer genomes

AbstractAlthough the genomes from any two human individuals are more than 99.99% identical at the sequence level, some structural variation can be observed. Differences between genomes include single nucleotide polymorphism (SNP), inversion and copy number changes (gain or loss of DNA). The latter can range from submicroscopic events (CNVs, at least 1kb in size) to complete chromosomal aneuploidies. Small copy number variations have often no (lethal) consequences to the cell, but a few were associated to disease susceptibility and phenotypic variations. Larger re-arrangements (i.e. complete chromosome gain) are frequently associated with more severe consequences on health such as genomic disorders and cancer. High-throughput technologies like DNA microarrays enable the detection of CNVs in a genome-wide fashion. Since the initial catalogue of CNVs in the human genome in 2006, there has been tremendous interest in CNVs both in the context of population and medical genetics. Understanding CNV patterns within and between human populations is essential to elucidate their possible contribution to disease. But genome analysis is a challenging task; the technology evolves rapidly creating needs for novel, efficient and robust analytical tools which need to be compared with existing ones. Also, while the link between CNV and disease has been established, the relative CNV contribution is not fully understood and the predisposition to disease from CNVs of the general population has not been yet investigated.During my PhD thesis, I worked on several aspects related to CNVs. As l will report in chapter 3, ! was interested in computational methods to detect CNVs from the general population. I had access to the CoLaus dataset, a population-based study with more than 6,000 participants from the Lausanne area. All these individuals were analysed on SNP arrays and extensive clinical information were available. My work explored existing CNV detection methods and I developed a variety of metrics to compare their performance. Since these methods were not producing entirely satisfactory results, I implemented my own method which outperformed two existing methods. I also devised strategies to combine CNVs from different individuals into CNV regions.I was also interested in the clinical impact of CNVs in common disease (chapter 4). Through an international collaboration led by the Centre Hospitalier Universitaire Vaudois (CHUV) and the Imperial College London I was involved as a main data analyst in the investigation of a rare deletion at chromosome 16p11 detected in obese patients. Specifically, we compared 8,456 obese patients and 11,856 individuals from the general population and we found that the deletion was accounting for 0.7% of the morbid obesity cases and was absent in healthy non- obese controls. This highlights the importance of rare variants with strong impact and provides new insights in the design of clinical studies to identify the missing heritability in common disease.Furthermore, I was interested in the detection of somatic copy number alterations (SCNA) and their consequences in cancer (chapter 5). This project was a collaboration initiated by the Ludwig Institute for Cancer Research and involved other groups from the Swiss Institute of Bioinformatics, the CHUV and Universities of Lausanne and Geneva. The focus of my work was to identify genes with altered expression levels within somatic copy number alterations (SCNA) in seven metastatic melanoma ceil lines, using CGH and SNP arrays, RNA-seq, and karyotyping. Very few SCNA genes were shared by even two melanoma samples making it difficult to draw any conclusions at the individual gene level. To overcome this limitation, I used a network-guided analysis to determine whether any pathways, defined by amplified or deleted genes, were common among the samples. Six of the melanoma samples were potentially altered in four pathways and five samples harboured copy-number and expression changes in components of six pathways. In total, this approach identified 28 pathways. Validation with two external, large melanoma datasets confirmed all but three of the detected pathways and demonstrated the utility of network-guided approaches for both large and small datasets analysis.RésuméBien que le génome de deux individus soit similaire à plus de 99.99%, des différences de structure peuvent être observées. Ces différences incluent les polymorphismes simples de nucléotides, les inversions et les changements en nombre de copies (gain ou perte d'ADN). Ces derniers varient de petits événements dits sous-microscopiques (moins de 1kb en taille), appelés CNVs (copy number variants) jusqu'à des événements plus large pouvant affecter des chromosomes entiers. Les petites variations sont généralement sans conséquence pour la cellule, toutefois certaines ont été impliquées dans la prédisposition à certaines maladies, et à des variations phénotypiques dans la population générale. Les réarrangements plus grands (par exemple, une copie additionnelle d'un chromosome appelée communément trisomie) ont des répercutions plus grave pour la santé, comme par exemple dans certains syndromes génomiques et dans le cancer. Les technologies à haut-débit telle les puces à ADN permettent la détection de CNVs à l'échelle du génome humain. La cartographie en 2006 des CNV du génome humain, a suscité un fort intérêt en génétique des populations et en génétique médicale. La détection de différences au sein et entre plusieurs populations est un élément clef pour élucider la contribution possible des CNVs dans les maladies. Toutefois l'analyse du génome reste une tâche difficile, la technologie évolue très rapidement créant de nouveaux besoins pour le développement d'outils, l'amélioration des précédents, et la comparaison des différentes méthodes. De plus, si le lien entre CNV et maladie a été établit, leur contribution précise n'est pas encore comprise. De même que les études sur la prédisposition aux maladies par des CNVs détectés dans la population générale n'ont pas encore été réalisées.Pendant mon doctorat, je me suis concentré sur trois axes principaux ayant attrait aux CNV. Dans le chapitre 3, je détaille mes travaux sur les méthodes d'analyses des puces à ADN. J'ai eu accès aux données du projet CoLaus, une étude de la population de Lausanne. Dans cette étude, le génome de plus de 6000 individus a été analysé avec des puces SNP et de nombreuses informations cliniques ont été récoltées. Pendant mes travaux, j'ai utilisé et comparé plusieurs méthodes de détection des CNVs. Les résultats n'étant pas complètement satisfaisant, j'ai implémenté ma propre méthode qui donne de meilleures performances que deux des trois autres méthodes utilisées. Je me suis aussi intéressé aux stratégies pour combiner les CNVs de différents individus en régions.Je me suis aussi intéressé à l'impact clinique des CNVs dans le cas des maladies génétiques communes (chapitre 4). Ce projet fut possible grâce à une étroite collaboration avec le Centre Hospitalier Universitaire Vaudois (CHUV) et l'Impérial College à Londres. Dans ce projet, j'ai été l'un des analystes principaux et j'ai travaillé sur l'impact clinique d'une délétion rare du chromosome 16p11 présente chez des patients atteints d'obésité. Dans cette collaboration multidisciplinaire, nous avons comparés 8'456 patients atteint d'obésité et 11 '856 individus de la population générale. Nous avons trouvés que la délétion était impliquée dans 0.7% des cas d'obésité morbide et était absente chez les contrôles sains (non-atteint d'obésité). Notre étude illustre l'importance des CNVs rares qui peuvent avoir un impact clinique très important. De plus, ceci permet d'envisager une alternative aux études d'associations pour améliorer notre compréhension de l'étiologie des maladies génétiques communes.Egalement, j'ai travaillé sur la détection d'altérations somatiques en nombres de copies (SCNA) et de leurs conséquences pour le cancer (chapitre 5). Ce projet fut une collaboration initiée par l'Institut Ludwig de Recherche contre le Cancer et impliquant l'Institut Suisse de Bioinformatique, le CHUV et les Universités de Lausanne et Genève. Je me suis concentré sur l'identification de gènes affectés par des SCNAs et avec une sur- ou sous-expression dans des lignées cellulaires dérivées de mélanomes métastatiques. Les données utilisées ont été générées par des puces ADN (CGH et SNP) et du séquençage à haut débit du transcriptome. Mes recherches ont montrées que peu de gènes sont récurrents entre les mélanomes, ce qui rend difficile l'interprétation des résultats. Pour contourner ces limitations, j'ai utilisé une analyse de réseaux pour définir si des réseaux de signalisations enrichis en gènes amplifiés ou perdus, étaient communs aux différents échantillons. En fait, parmi les 28 réseaux détectés, quatre réseaux sont potentiellement dérégulés chez six mélanomes, et six réseaux supplémentaires sont affectés chez cinq mélanomes. La validation de ces résultats avec deux larges jeux de données publiques, a confirmée tous ces réseaux sauf trois. Ceci démontre l'utilité de cette approche pour l'analyse de petits et de larges jeux de données.Résumé grand publicL'avènement de la biologie moléculaire, en particulier ces dix dernières années, a révolutionné la recherche en génétique médicale. Grâce à la disponibilité du génome humain de référence dès 2001, de nouvelles technologies telles que les puces à ADN sont apparues et ont permis d'étudier le génome dans son ensemble avec une résolution dite sous-microscopique jusque-là impossible par les techniques traditionnelles de cytogénétique. Un des exemples les plus importants est l'étude des variations structurales du génome, en particulier l'étude du nombre de copies des gènes. Il était établi dès 1959 avec l'identification de la trisomie 21 par le professeur Jérôme Lejeune que le gain d'un chromosome supplémentaire était à l'origine de syndrome génétique avec des répercussions graves pour la santé du patient. Ces observations ont également été réalisées en oncologie sur les cellules cancéreuses qui accumulent fréquemment des aberrations en nombre de copies (telles que la perte ou le gain d'un ou plusieurs chromosomes). Dès 2004, plusieurs groupes de recherches ont répertorié des changements en nombre de copies dans des individus provenant de la population générale (c'est-à-dire sans symptômes cliniques visibles). En 2006, le Dr. Richard Redon a établi la première carte de variation en nombre de copies dans la population générale. Ces découvertes ont démontrées que les variations dans le génome était fréquentes et que la plupart d'entre elles étaient bénignes, c'est-à-dire sans conséquence clinique pour la santé de l'individu. Ceci a suscité un très grand intérêt pour comprendre les variations naturelles entre individus mais aussi pour mieux appréhender la prédisposition génétique à certaines maladies.Lors de ma thèse, j'ai développé de nouveaux outils informatiques pour l'analyse de puces à ADN dans le but de cartographier ces variations à l'échelle génomique. J'ai utilisé ces outils pour établir les variations dans la population suisse et je me suis consacré par la suite à l'étude de facteurs pouvant expliquer la prédisposition aux maladies telles que l'obésité. Cette étude en collaboration avec le Centre Hospitalier Universitaire Vaudois a permis l'identification d'une délétion sur le chromosome 16 expliquant 0.7% des cas d'obésité morbide. Cette étude a plusieurs répercussions. Tout d'abord elle permet d'effectuer le diagnostique chez les enfants à naître afin de déterminer leur prédisposition à l'obésité. Ensuite ce locus implique une vingtaine de gènes. Ceci permet de formuler de nouvelles hypothèses de travail et d'orienter la recherche afin d'améliorer notre compréhension de la maladie et l'espoir de découvrir un nouveau traitement Enfin notre étude fournit une alternative aux études d'association génétique qui n'ont eu jusqu'à présent qu'un succès mitigé.Dans la dernière partie de ma thèse, je me suis intéressé à l'analyse des aberrations en nombre de copies dans le cancer. Mon choix s'est porté sur l'étude de mélanomes, impliqués dans le cancer de la peau. Le mélanome est une tumeur très agressive, elle est responsable de 80% des décès des cancers de la peau et est souvent résistante aux traitements utilisés en oncologie (chimiothérapie, radiothérapie). Dans le cadre d'une collaboration entre l'Institut Ludwig de Recherche contre le Cancer, l'Institut Suisse de Bioinformatique, le CHUV et les universités de Lausanne et Genève, nous avons séquencés l'exome (les gènes) et le transcriptome (l'expression des gènes) de sept mélanomes métastatiques, effectués des analyses du nombre de copies par des puces à ADN et des caryotypes. Mes travaux ont permis le développement de nouvelles méthodes d'analyses adaptées au cancer, d'établir la liste des réseaux de signalisation cellulaire affectés de façon récurrente chez le mélanome et d'identifier deux cibles thérapeutiques potentielles jusqu'alors ignorées dans les cancers de la peau.

Formation of a direct arterial blood supply to the anterior pituitary gland following complete or partial interruption of the hypophyseal portal vessels

If regions of the anterior pituitary gland received systemic blood via a direct arterial blood supply these regions would escape hypothalamic regulation and thus be a sequela in endocrine disorders. Since, in the untreated rat, all of the blood supply to the anterior pituitary gland is via the hypophyseal portal vessels, we hypothesized that partial interruption of the portal vessels could provoke the establishment of a direct arterial blood supply (arteriogenesis). We utilized the injection of polystyrene microspheres (15 or 9 micron diameter) into the left ventricle of the heart to test this hypothesis. Microspheres are trapped in the first capillary plexus they reach since they are too large to traverse the capillaries. No microspheres reached the anterior pituitary gland of control rats, a finding consistent with the fact that the anterior pituitary gland receives all of its blood supply via the hypophyseal portal blood vessels. Microspheres were observed in the primary portal capillary plexus in the infundibulum (median eminence), infundibular stalk (pituitary stalk), and infundibular process (pars nervosa), the first capillary plexus which they reached. A lesion of the medial basal hypothalamus (MBH) which destroyed the long portal vessels did not result in arteriogenesis since few, if any, microspheres were observed in the anterior pituitary gland. We confirmed, using vascular casts, that these lesions resulted in the long-term destruction of the primary portal capillaries in the infundibulum and of the long portal vessels. In MBH-lesioned animals it appears that all of the blood supply of the anterior pituitary gland is via short portal vessels arising from the infundibular stem and process.(ABSTRACT TRUNCATED AT 250 WORDS)

Six of 12 Relapsed or Refractory Indolent Lymphoma Patients Treated 10 Years Ago with 131I-Tositumomab Remain in Complete Remission.

The purpose of our study was to update the safety and efficacy results of radioimmunotherapy in relapsed or resistant indolent or transformed non-Hodgkin lymphoma. Methods: More than 9 y ago, we treated 12 indolent and 4 transformed, relapsed or refractory lymphoma patients with a single administration of nonmyeloablative therapy with tositumomab and I-131-tositumomab. The 16 patients had a mean of 3.1 (range, 1-6) previous chemotherapy and antibody treatments. Results: Six of 12 relapsed indolent lymphoma patients remain disease-free a mean of 9.8 y (range, 8.6-10.7 y) after radioimmunotherapy. Three of 4 transformed lymphoma patients progressed after radioimmunotherapy, and 1 patient had a partial response of 10 mo. Conclusion: Optimal patient benefit might be obtained in indolent lymphoma when administering radioimmunotherapy up-front in combination with chemotherapy and rituximab treatment. However, these results show that radioimmunotherapy alone achieved long-lasting remissions in 6 of 12 (50%) indolent lymphoma patients in relapse after 1 or multiple chemotherapies.

Tandem chimerism as a means to increase protein complexity 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.

Oviposition-induced changes in Arabidopsis genome expression: anticipating your enemy?

Plants have evolved exquisite ways to detect their enemies and are able to induce defenses responses tailored to their specific aggressors. Insect eggs deposited on a leaf represent a future threat as larvae hatching from the egg will ultimately feed on the plant. Although direct and indirect defenses towards oviposition have been documented, our knowledge of the molecular changes triggered by egg deposition is limited. Using a whole-genome microarray, we recently analyzed the expression profile of Arabidopsis thaliana leaves after oviposition by two pierid butterflies. Eggs laid by the large white Pieris brassicae modified the expression of hundreds of genes. The transcript signature included defense and stress-related genes that were also induced in plants experiencing localized cell death. Further analyses revealed that cellular changes associated with a hypersensitive response occur at the site of egg deposition and that they are triggered by egg-derived elicitors. Our study brings molecular evidence for previous observations of oviposition-induced necrosis in other plant species and might illustrate a direct defense of the plant against the egg. In this addendum, we discuss the relevance of the oviposition-induced gene expression changes and the possibility that plants use eggs as cues to anticipate their enemies.

Biochemical characterization of FANCD2 and FANCM in the Fanconi Anaemia pathway

Abstract : The maintenance of genome stability is a challenge for all living organisms. DNA is regularly subjected to chemical alterations by both endogenous and exogenous DNA damaging agents. If left unrepaired, these lesions will create mutations or lead to chromosomal instability. DNA crosslinking agents probably bring about the most toxic lesions. By linking covalently the two strands of DNA, crosslinking agents will impede essential cellular processes such as replication and transcription. Cells from Fanconi anaemia patients are extremely sensitive to these agents. Fanconi anaemia (FA) is a rare chromosomal instability disorder that leads to developmental defects, pancytopenia and cancer susceptibility. FA is a genetically heterogeneous disease with thirteen complementation groups identified. Proteins encoded by the FA genes work together in the FA pathway. Eight of these proteins form the FA core complex (FANC-A, B, C,E, F, G, L and -M), whose integrity is required to monoubiquitinate FANCD2 and FANCI in response to DNA damage. The hypersensitivity of FA cells to crosslinking agents, which perturb the progression of replication forks, has led to the hypothesis that FA proteins play a crucial role in the response to replication stress. However, at the molecular level, the functions of the FA pathway remain largely unknown. Our efforts were first focused on the characterization of FANCD2, "the key effector of the FA pathway". Using different substrates, we found that in vitro, purified hFANCD2 preferentially binds single strand DNA and double strand DNA extremities. Concomitantly, FANCM was identified as a new component of the FA core complex. Moreover FANCM was shown to have specific branch migration activities and probably a role as a "landing platform" on DNA for the other components of the core complex. By using FANCM mutants carrying deletions within the internal domain, we investigated the role of FANCM as a DNA anchor protein for the core complex. We observed that indeed, a specific part of the internal domain of FANCM interacts with components of the core complex. Finally, in collaboration with Weidong Wang's lab we characterized two new components of the FA pathway: FAAP10 and FAAP16. As a heterodimer these two proteins show affinity for dsDNA, and anneal complementary oligonucleotides in vitro. Moreover these proteins can associate with FANCM via a part of its internal domain. We find that FANCM, FAAP 10 and FAAP 16 can co-exist on the branch point of replication and recombination intermediates, and that FAAP10 and FAAP16 stimulate replication fork reversal by FANCM. These results suggest that FANCM may function as a landing platform for the core complex. After loading on DNA, the core complex can activate FANCD2 through monoubiquitination leading to its recruitment to the site of damage. Since ssDNA and double strand breaks are intermediates that are generated as a consequence of collapsed replication forks, FANCD2 by binding to ds DNA ends and ssDNA could protect such structures from the recombination repair machinery and prevent unscheduled recombination events. Alternatively, FANCD2 could avoid nucleases from gaining access to collapsed forks, preserving the DNA in state that can be used as a starting point for resumption of DNA synthesis. The overall comprehension of the FA pathway is far from been complete. Our results unravel new aspects of Fanconi Anaemia, which hopefully in the near future will address keys questions leading to a better understanding of the fascinating Fanconi Anaemia. Résumé : Le maintien de l'intégrité du génome est fondamentale chez tous les organismes vivants. L'ADN est constamment altéré par des composés aussi bien endogènes qu'exogènes. Si ces altérations ne sont pas réparées, elles peuvent conduire à l'apparition de mutations, ainsi qu'à une instabilité génomique accrue. Les lésions les plus sévères qui peuvent survenir sur l'ADN, sont les pontages inter caténaires. Des agents pontants en liant de façon covalente les deux brins d'ADN, vont empêcher le déroulement normal de processus cellulaires essentiels tels que la réplication ou la transcription. La compréhension des mécanismes permettant à la cellule de tolérer et réparer ces lésions est primordiale, notamment dans le cas des patients atteints de l'anémie de Fanconi qui présentent une très grande sensibilité à ces composés pontants. L'anémie de Fanconi est une maladie génétique rare appartenant à un groupe de pathologies associées à une grande instabilité chromosomique. Les patients atteints de l'anémie de Fanconi présentent des malformations du squelette, une pancytopénie et une forte propension à la survenue de cancer. L'anémie de Fanconi est génétiquement très hétérogène. À ce jour, 13 gènes codant pour 13 protéines FANC différentes ont été identifiés. Huit de ces protéines fonctionnent ensemble au sein d'un complexe (nommé le complexe FANC) ayant pour but de monoubiquitiner FANCD2 et FANCI en réponse à la formation de lésions sur l'ADN. L'extrême sensibilité des cellules de patients atteints de l'anémie de Fanconi à ces agents pontant l'ADN suggère l'implication des protéines FANC dans la réponse cellulaire suite à une stress réplicatif. Cependant, le rôle moléculaire exact de ces protéines demeure encore inconnu. Après purification, nous avons observé que FANCD2 était capable de lier l'ADN simple brin, ainsi que les extrémités d'ADN in vitro. Dans le même temps, FANCM fut identifié comme appartenant au complexe FANC. FANCM est décrit comme une translocase capable de promouvoir le déplacement de point de jonction dans des structures d'ADN spécifiques in vitro. De plus, en se liant à l'ADN, FANCM peut agir comme une plateforme pour les autres protéines FANC, leur permettant ainsi d'être adressées à l'ADN. En créant des protéines FANCM recombinantes ayant des délétions dans le domaine interne, nous avons pu observer que certaines protéines du complexe FANC se fixent à des sites spécifiques sur le domaine interne de FANCM. Enfin, au travers d'une collaboration, nous avons été amenés à caractériser deux nouvelles protéines appartenant au complexe FANC : FAAP 10 et FAAP16. Elles s'associent à FANCM par l'intermédiaire du domaine interne, et forment ainsi un hétérotrimére. La présence de FAAP10 et FAAP16 n'affecte pas la liaison de FANCM à l'ADN, mais semble potentialiser son activité de régression in vitro. FANCM semble donc fonctionner comme une plateforme pour les autres composants du complexe FANC. Ces derniers, une fois liés à l'ADN permettent la monoubiquitination de FANCD2 et son recrutement au site lésé de l'ADN. FANCD2 en se liant de façon préférentielle à l'ADN simple brin et aux extrémités d'ADN qui sont générés lors de l'arrêt et du démantèlement d'une fourche de réplication, pourrait protéger ces même fourches de réplication arrêtées, d'évènements de recombinaison aléatoires. Nos résultats apportent de nouveaux éléments concernant les mécanismes moléculaires de l'anémie de Fanconi. Enfin, l'étude de l'anémie de Fanconi permet aussi de mieux comprendre les mécanismes mis en place par la cellule pour tolérer des lésions survenant lors de la réplication.

Global DNA hypomethylation coupled to repressive chromatin domain formation and gene silencing in breast cancer.

While genetic mutation is a hallmark of cancer, many cancers also acquire epigenetic alterations during tumorigenesis including aberrant DNA hypermethylation of tumor suppressors, as well as changes in chromatin modifications as caused by genetic mutations of the chromatin-modifying machinery. However, the extent of epigenetic alterations in cancer cells has not been fully characterized. Here, we describe complete methylome maps at single nucleotide resolution of a low-passage breast cancer cell line and primary human mammary epithelial cells. We find widespread DNA hypomethylation in the cancer cell, primarily at partially methylated domains (PMDs) in normal breast cells. Unexpectedly, genes within these regions are largely silenced in cancer cells. The loss of DNA methylation in these regions is accompanied by formation of repressive chromatin, with a significant fraction displaying allelic DNA methylation where one allele is DNA methylated while the other allele is occupied by histone modifications H3K9me3 or H3K27me3. Our results show a mutually exclusive relationship between DNA methylation and H3K9me3 or H3K27me3. These results suggest that global DNA hypomethylation in breast cancer is tightly linked to the formation of repressive chromatin domains and gene silencing, thus identifying a potential epigenetic pathway for gene regulation in cancer cells.

A genome-wide screen for interactions reveals a new locus on 4p15 modifying the effect of waist-to-hip ratio on total cholesterol.

Recent genome-wide association (GWA) studies described 95 loci controlling serum lipid levels. These common variants explain ∼25% of the heritability of the phenotypes. To date, no unbiased screen for gene-environment interactions for circulating lipids has been reported. We screened for variants that modify the relationship between known epidemiological risk factors and circulating lipid levels in a meta-analysis of genome-wide association (GWA) data from 18 population-based cohorts with European ancestry (maximum N = 32,225). We collected 8 further cohorts (N = 17,102) for replication, and rs6448771 on 4p15 demonstrated genome-wide significant interaction with waist-to-hip-ratio (WHR) on total cholesterol (TC) with a combined P-value of 4.79×10(-9). There were two potential candidate genes in the region, PCDH7 and CCKAR, with differential expression levels for rs6448771 genotypes in adipose tissue. The effect of WHR on TC was strongest for individuals carrying two copies of G allele, for whom a one standard deviation (sd) difference in WHR corresponds to 0.19 sd difference in TC concentration, while for A allele homozygous the difference was 0.12 sd. Our findings may open up possibilities for targeted intervention strategies for people characterized by specific genomic profiles. However, more refined measures of both body-fat distribution and metabolic measures are needed to understand how their joint dynamics are modified by the newly found locus.

The complementary strand of the human T-cell leukemia virus type 1 RNA genome encodes a bZIP transcription factor that down-regulates viral transcription.

The RNA genome of the human T-cell leukemia virus type 1 (HTLV-1) codes for proteins involved in infectivity, replication, and transformation. We report in this study the characterization of a novel viral protein encoded by the complementary strand of the HTLV-1 RNA genome. This protein, designated HBZ (for HTLV-1 bZIP factor), contains a N-terminal transcriptional activation domain and a leucine zipper motif in its C terminus. We show here that HBZ is able to interact with the bZIP transcription factor CREB-2 (also called ATF-4), known to activate the HTLV-1 transcription by recruiting the viral trans-activator Tax on the Tax-responsive elements (TxREs). However, we demonstrate that the HBZ/CREB-2 heterodimers are no more able to bind to the TxRE and cyclic AMP response element sites. Taking these findings together, the functional inactivation of CREB-2 by HBZ is suggested to contribute to regulation of the HTLV-1 transcription. Moreover, the characterization of a minus-strand gene protein encoded by HTLV-1 has never been reported until now.

A genome-wide association study of anorexia nervosa.

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome-wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2907 cases with AN from 14 countries (15 sites) and 14 860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery data sets. Seventy-six (72 independent) single nucleotide polymorphisms were taken forward for in silico (two data sets) or de novo (13 data sets) replication genotyping in 2677 independent AN cases and 8629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication data sets comprised 5551 AN cases and 21 080 controls. AN subtype analyses (1606 AN restricting; 1445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01 × 10(-7)) in SOX2OT and rs17030795 (P=5.84 × 10(-6)) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76 × 10(-)(6)) between CUL3 and FAM124B and rs1886797 (P=8.05 × 10(-)(6)) near SPATA13. Comparing discovery with replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4 × 10(-6)), strongly suggesting that true findings exist but our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.

Genome-wide RNA polymerase II profiles and RNA accumulation reveal kinetics of transcription and associated epigenetic changes during diurnal cycles.

Interactions of cell-autonomous circadian oscillators with diurnal cycles govern the temporal compartmentalization of cell physiology in mammals. To understand the transcriptional and epigenetic basis of diurnal rhythms in mouse liver genome-wide, we generated temporal DNA occupancy profiles by RNA polymerase II (Pol II) as well as profiles of the histone modifications H3K4me3 and H3K36me3. We used these data to quantify the relationships of phases and amplitudes between different marks. We found that rhythmic Pol II recruitment at promoters rather than rhythmic transition from paused to productive elongation underlies diurnal gene transcription, a conclusion further supported by modeling. Moreover, Pol II occupancy preceded mRNA accumulation by 3 hours, consistent with mRNA half-lives. Both methylation marks showed that the epigenetic landscape is highly dynamic and globally remodeled during the 24-hour cycle. While promoters of transcribed genes had tri-methylated H3K4 even at their trough activity times, tri-methylation levels reached their peak, on average, 1 hour after Pol II. Meanwhile, rhythms in tri-methylation of H3K36 lagged transcription by 3 hours. Finally, modeling profiles of Pol II occupancy and mRNA accumulation identified three classes of genes: one showing rhythmicity both in transcriptional and mRNA accumulation, a second class with rhythmic transcription but flat mRNA levels, and a third with constant transcription but rhythmic mRNAs. The latter class emphasizes widespread temporally gated posttranscriptional regulation in the mouse liver.

The genome sequence of taurine cattle: a window to ruminant biology and evolution.

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

Genome-wide association study of recurrent major depressive disorder in two European case-control cohorts.

Major depressive disorder (MDD) is a highly prevalent disorder with substantial heritability. Heritability has been shown to be substantial and higher in the variant of MDD characterized by recurrent episodes of depression. Genetic studies have thus far failed to identify clear and consistent evidence of genetic risk factors for MDD. We conducted a genome-wide association study (GWAS) in two independent datasets. The first GWAS was performed on 1022 recurrent MDD patients and 1000 controls genotyped on the Illumina 550 platform. The second was conducted on 492 recurrent MDD patients and 1052 controls selected from a population-based collection, genotyped on the Affymetrix 5.0 platform. Neither GWAS identified any SNP that achieved GWAS significance. We obtained imputed genotypes at the Illumina loci for the individuals genotyped on the Affymetrix platform, and performed a meta-analysis of the two GWASs for this common set of approximately half a million SNPs. The meta-analysis did not yield genome-wide significant results either. The results from our study suggest that SNPs with substantial odds ratio are unlikely to exist for MDD, at least in our datasets and among the relatively common SNPs genotyped or tagged by the half-million-loci arrays. Meta-analysis of larger datasets is warranted to identify SNPs with smaller effects or with rarer allele frequencies that contribute to the risk of MDD.

Developmental constraints on vertebrate genome evolution.

Constraints in embryonic development are thought to bias the direction of evolution by making some changes less likely, and others more likely, depending on their consequences on ontogeny. Here, we characterize the constraints acting on genome evolution in vertebrates. We used gene expression data from two vertebrates: zebrafish, using a microarray experiment spanning 14 stages of development, and mouse, using EST counts for 26 stages of development. We show that, in both species, genes expressed early in development (1) have a more dramatic effect of knock-out or mutation and (2) are more likely to revert to single copy after whole genome duplication, relative to genes expressed late. This supports high constraints on early stages of vertebrate development, making them less open to innovations (gene gain or gene loss). Results are robust to different sources of data -- gene expression from microarrays, ESTs, or in situ hybridizations; and mutants from directed KO, transgenic insertions, point mutations, or morpholinos. We determine the pattern of these constraints, which differs from the model used to describe vertebrate morphological conservation ("hourglass" model). While morphological constraints reach a maximum at mid-development (the "phylotypic" stage), genomic constraints appear to decrease in a monotonous manner over developmental time.