Rare genomic structural variants in complex disease: lessons from the replication of associations with obesity.

The limited ability of common variants to account for the genetic contribution to complex disease has prompted searches for rare variants of large effect, to partly explain the 'missing heritability'. Analyses of genome-wide genotyping data have identified genomic structural variants (GSVs) as a source of such rare causal variants. Recent studies have reported multiple GSV loci associated with risk of obesity. We attempted to replicate these associations by similar analysis of two familial-obesity case-control cohorts and a population cohort, and detected GSVs at 11 out of 18 loci, at frequencies similar to those previously reported. Based on their reported frequencies and effect sizes (OR≥25), we had sufficient statistical power to detect the large majority (80%) of genuine associations at these loci. However, only one obesity association was replicated. Deletion of a 220 kb region on chromosome 16p11.2 has a carrier population frequency of 2×10(-4) (95% confidence interval [9.6×10(-5)-3.1×10(-4)]); accounts overall for 0.5% [0.19%-0.82%] of severe childhood obesity cases (P = 3.8×10(-10); odds ratio = 25.0 [9.9-60.6]); and results in a mean body mass index (BMI) increase of 5.8 kg.m(-2) [1.8-10.3] in adults from the general population. We also attempted replication using BMI as a quantitative trait in our population cohort; associations with BMI at or near nominal significance were detected at two further loci near KIF2B and within FOXP2, but these did not survive correction for multiple testing. These findings emphasise several issues of importance when conducting rare GSV association, including the need for careful cohort selection and replication strategy, accurate GSV identification, and appropriate correction for multiple testing and/or control of false discovery rate. Moreover, they highlight the potential difficulty in replicating rare CNV associations across different populations. Nevertheless, we show that such studies are potentially valuable for the identification of variants making an appreciable contribution to complex disease.

Sympatric speciation in palms on an oceanic island.

The origin of species diversity has challenged biologists for over two centuries. Allopatric speciation, the divergence of species resulting from geographical isolation, is well documented. However, sympatric speciation, divergence without geographical isolation, is highly controversial. Claims of sympatric speciation must demonstrate species sympatry, sister relationships, reproductive isolation, and that an earlier allopatric phase is highly unlikely. Here we provide clear support for sympatric speciation in a case study of two species of palm (Arecaceae) on an oceanic island. A large dated phylogenetic tree shows that the two species of Howea, endemic to the remote Lord Howe Island, are sister taxa and diverged from each other well after the island was formed 6.9 million years ago. During fieldwork, we found a substantial disjunction in flowering time that is correlated with soil preference. In addition, a genome scan indicates that few genetic loci are more divergent between the two species than expected under neutrality, a finding consistent with models of sympatric speciation involving disruptive/divergent selection. This case study of sympatric speciation in plants provides an opportunity for refining theoretical models on the origin of species, and new impetus for exploring putative plant and animal examples on oceanic islands.

The tectonically controlled emplacement of a vertically sheeted gabbro-pyroxenite intrusion: Feeder-zone of an ocean-island volcano (Fuerteventura, Canary Islands)

The Miocene PX1 gabbro-pyroxenite pluton, Fuerteventura, Canary Islands, is a 3.5 x 5.5 km shallow-level intrusion (0.15-0.2 GPa and 1100-1120 degrees C), interpreted as the feeder-zone to an ocean-island volcano. It displays a vertical magmatic banding expressed in five 50 to 100 metre-wide NNE-SSW trending alkaline gabbro sequences alternating with pyroxenites. This emplacement geometry was controlled by brittle to ductile shear zones, generated by a regional E-W extensional tectonic setting that affected Fuerteventura during the Miocene. At a smaller scale, the PX1 gabbro and pyroxenite bands consist of metre-thick differentiation units, which suggest emplacement by periodic injection of magma pulses as vertical dykes that amalgamated, similarly to a sub-volcanic sheeted dyke complex. Individual dykes underwent internal differentiation following a solidification front parallel to the dyke edges. This solidification front may have been favoured by a significant lateral/horizontal thermal gradient, expressed by the vertical banding in the gabbros, the fractionation asymmetry within individual dykes and the migmatisation of the wall rocks. Pyroxenitic layers result from the fractionation and accumulation of clinopyroxene +/- olivine +/- plagioclase crystals from a mildly alkaline basaltic liquid. They are interpreted as truncated differentiation sequences, from which residual melts were extracted at various stages of their chemical evolution by subsequent dyke intrusions, either next to or within the crystallising unit. Compaction and squeezing of the crystal mush is ascribed to the incoming and inflating magma pulses. The expelled interstitial liquid was likely collected and erupted along with the magma flowing through the newly injected dykes. Clinopyroxene mineral orientation - as evidenced by EBSD and micro X-ray tomography investigations - displays a marked pure-shear component, supporting the interpretation of the role of compaction in the generation of the pyroxenites. Conversely, gabbro sequences underwent minor melt extraction and are believed to represent crystallised coalesced magma batches emplaced at lower rates at the end of eruptive cycles. Clinopyroxene orientations in gabbros record a simple shear component suggesting syn-magmatic deformation parallel to observed NNE-SSW trending shear zones induced by the regional tensional stress field. This emplacement model implies a crystallisation time of 1 to 5 years for individual dykes, consistent with PX1 emplacement over less than 0.5 My. A minimum amount of approximately 150 km(3) of magma is needed to generate the pluton, part of it having been erupted through the Central Volcanic Centre of Fuerteventura. If the regional extensional tectonic regime controls the PX1 feeder-zone initiation and overall geometry, rates and volumes of magma depend on other, source-related factors. High injection rates are likely to induce intrusion growth rates larger than could be accommodated by the regional extension. In this case, dyke intrusion by propagation of a weak tip, combined with the inability of magma to circulate through previously emplaced and crystallised dykes could result in an increase of non-lithostatic pressure on previously emplaced mushy dyke walls; thus generating strong pure-shear compaction within the pluton feeder-zone and interstitial melt expulsion. These compaction-dominated processes are recorded by the cumulitic pyroxenite bands. (C) 2010 Elsevier B.V. All rights reserved.

Soft tissue sarcomas with complex genomic profiles.

Soft tissue sarcomas (STS) with complex genomic profiles (50% of all STS) are predominantly composed of spindle cell/pleomorphic sarcomas, including leiomyosarcoma, myxofibrosarcoma, pleomorphic liposarcoma, pleomorphic rhabdomyosarcoma, malignant peripheral nerve sheath tumor, angiosarcoma, extraskeletal osteosarcoma, and spindle cell/pleomorphic unclassified sarcoma (previously called spindle cell/pleomorphic malignant fibrous histiocytoma). These neoplasms show, characteristically, gains and losses of numerous chromosomes or chromosome regions, as well as amplifications. Many of them share recurrent aberrations (e.g., gain of 5p13-p15) that seem to play a significant role in tumor progression and/or metastatic dissemination. In this paper, we review the cytogenetic, molecular genetic, and clinicopathologic characteristics of the most common STS displaying complex genomic profiles. Features of diagnostic or prognostic relevance will be discussed when needed.

Genomic, Epigenomic and Transcriptomic Molecular Characterization of Splenic Marginal Zone Lymphoma Reveals Recurrent Abnormalities in miR-182

Splenic marginal zone lymphoma (SMZL) is a low grade B-cell non-Hodgkin's lymphoma. The molecular pathology of this entity remains poorly understood. To characterise this lymphoma at the molecular level, we performed an integrated analysis of 1) genome wide genetic copy number alterations 2) gene expression profiles and 3) epigenetic DNA methylation profiles.We have previously shown that SMZL is characterised by recurrent alterations of chromosomes 7q, 6q, 3q, 9q and 18; however, gene resolution oligonucleotide array comparative genomic hybridisation did not reveal evidence of cryptic amplification or deletion in these regions. The most frequently lost 7q32 region contains a cluster of miRNAs. qRT-PCR revealed that three of these (miR-182/96/183) show underexpression in SMZL, and miR-182 is somatically mutated in >20% of cases of SMZL, as well as in >20% of cases of follicular lymphoma, and between 5-15% of cases of chronic lymphocytic leukaemia, MALT-lymphoma and hairy cell leukaemia. We conclude that miR-182 is a strong candidate novel tumour suppressor miRNA in lymphoma.The overall gene expression signature of SMZL was found to be strongly distinct fromthose of other lymphomas. Functional analysis of gene expression data revealed SMZL to be characterised by abnormalities in B-cell receptor signalling (especially through the CD19/21-PI3K/AKT pathway) and apoptotic pathways. In addition, genes involved in the response to viral infection appeared upregulated. SMZL shows a unique epigenetic profile, but analysis of differentially methylated genes showed few with methylation related transcriptional deregulation, suggesting that DNA methylation abnormalities are not a critical component of the SMZL malignant phenotype.

Systematic identification of genomic markers of drug sensitivity in cancer cells.

Clinical responses to anticancer therapies are often restricted to a subset of patients. In some cases, mutated cancer genes are potent biomarkers for responses to targeted agents. Here, to uncover new biomarkers of sensitivity and resistance to cancer therapeutics, we screened a panel of several hundred cancer cell lines--which represent much of the tissue-type and genetic diversity of human cancers--with 130 drugs under clinical and preclinical investigation. In aggregate, we found that mutated cancer genes were associated with cellular response to most currently available cancer drugs. Classic oncogene addiction paradigms were modified by additional tissue-specific or expression biomarkers, and some frequently mutated genes were associated with sensitivity to a broad range of therapeutic agents. Unexpected relationships were revealed, including the marked sensitivity of Ewing's sarcoma cells harbouring the EWS (also known as EWSR1)-FLI1 gene translocation to poly(ADP-ribose) polymerase (PARP) inhibitors. By linking drug activity to the functional complexity of cancer genomes, systematic pharmacogenomic profiling in cancer cell lines provides a powerful biomarker discovery platform to guide rational cancer therapeutic strategies.

Characterization of a Trypanosoma cruzi genomic fragment complementary to several species-specific mRNAs but different from the spliced leader sequence

We have isolated a clone of Trypanosoma cruzi genimic DNA, lambda 3b2-5, which contains sequences that are reiterated in the genome. Northtern blot analysis showed that clone 3b2-5 hybridizes to 1,200-5,000 bases different mRNA species. The number of mRNAs species hybridized to clone 3b2-5 exceeds its coding capacity showing that this clone carries sequences that are common to several mRNAs species and conserved in the poly A(+) RNA. These sequences are not homologous to the T. cruzi spliced leader sequence, since clone 3b2-5 hybridize to a synthetic 20 nucleotice complementary to the spliced leader sequence. Clone 3b2-5 does not hybridize to DNA and RNA from several genera of Trypanosomatidae and other Trypanosoma species indicating that it carries T. cruzi species-specific sequences.

The newly defined "Greenstone Unit'' of the Aiguilles Rouges massif (western Alps): remnant of an Early Palaeozoic oceanic island-arc ?

In the southwestern part of the Aiguilles Rouges massif (pre-Alpine basement of the Helvetic realm, western Alps), a metavolcanic sequence, newly defined as the ``Greenstone Unit'',is exposed in two NS trending belts of several 100 metres in thickness. It consists of epidote amphibolites, partly epidote and/or calcic amphibole-bearing greenschists, and small amounts of alkali feldspar-bearing greenschists, which underwent low- to medium-grade metamorphism during Visean oblique collision. Metamorphic calcic amphiboles and epidotes show strong chemical zoning, whereas metamorphic plagioclase is exclusively albitic in composition (An 1-3). The SiO2 content of the subalkaline tholeiitic to calc-alkaline suite ranges continuously from 44 wt% to 73 wt%,but andesitic rocks predominate. The majority of samples have chemical compositions close to recent subduction-related lavas; some are even restricted to recent oceanic arcs (extremely low Ta and Nb contents, high La/Nb and Th/Ta ratios). But several basaltic to basalto-andesitic samples resemble continental tholeiites (low Th/Ta, La/Nb ratio). As it is very probable that both lava types are to some extent contemporaneous, it is proposed that the Greenstone Unit represents a former oceanic volcanic are which temporarily underwent extension during which emplacement of continental tholeiite-like rocks occurred. The cause of the extension remains ambiguous. Considering palaeotectonic significance and age of other metavolcanic units in the Aiguilles Rouges massif, the Greenstone Unit most likely formed in the Early Palaeozoic.

Effect of training-sample size and classification difficulty on the accuracy of genomic predictors.

Introduction: As part of the MicroArray Quality Control (MAQC)-II project, this analysis examines how the choice of univariate feature-selection methods and classification algorithms may influence the performance of genomic predictors under varying degrees of prediction difficulty represented by three clinically relevant endpoints. Methods: We used gene-expression data from 230 breast cancers (grouped into training and independent validation sets), and we examined 40 predictors (five univariate feature-selection methods combined with eight different classifiers) for each of the three endpoints. Their classification performance was estimated on the training set by using two different resampling methods and compared with the accuracy observed in the independent validation set. Results: A ranking of the three classification problems was obtained, and the performance of 120 models was estimated and assessed on an independent validation set. The bootstrapping estimates were closer to the validation performance than were the cross-validation estimates. The required sample size for each endpoint was estimated, and both gene-level and pathway-level analyses were performed on the obtained models. Conclusions: We showed that genomic predictor accuracy is determined largely by an interplay between sample size and classification difficulty. Variations on univariate feature-selection methods and choice of classification algorithm have only a modest impact on predictor performance, and several statistically equally good predictors can be developed for any given classification problem.

Trypanosoma rangeli (Tejera, 1920) isolated from a sylvatic rodent (Echimys dasythrix) in Santa Catarina island, Santa Catarina state: first report of this trypanosome in southern Brazil

A trypanosome strain isolated from a sylvatic rodent (Echimys dasythrix) from Santa Catarina Island (Santa Catarina State, Brazil) was characterized by the following methods: experimental transmission and development in invertebrate hosts, morphometry, cross protection, complement sensitivity, lectin agglutination and isoenzyme profiles. Comparasions were made with standard Trypanosoma cruzi and T. rangeli strains. All methods except isoenzyne analysis led to the identification of the isolate as T. rangeli. The isoenzyme differences found could be explained on the basis of polymorphism. Therefore this is the first report of T. rangeli in southern Brazil, increasing the geographical distribution of this parasite.

Genomic determinants of the efficiency of internal ribosomal entry sites of viral and cellular origin

Summary : Internal ribosome entry sites (IRES) are used by viruses as a strategy to bypass inhibition of cap-dependent translation that commonly results from viral infection. IRES are also used in eukaryotic cells to control mRNA translation under conditions of cellular stress (apoptosis, heat shock) or during the G2 phase of the cell cycle when general protein synthesis is inhibited. Variation in cellular expression levels has been shown to be inherited. Expression is controlled, among others, by transcriptional factors and by the efficiency of cap-mediated translation and ribosome activity. We aimed at identifying genomic determinants of variability in IRES-mediated translation of two representative IRES [Encephalomyocarditis virus (EMCV) and X-linked Inhibitor-of-Apoptosis (XIAP) IRES]. We used bicistronic lentiviral constructions expressing two fluorescent reporter transgenes. Lentiviruses were used to transduce seven different laboratory cell lines and B lymphoblastoid cell lines from the Centre d'Etude du Polymorphisme Humain (CEPH; 15 pedigrees; n=209); representing an in vitro approach to family structure allowing genome scan analyses. The relative expression of the two markers was assessed by FACS. IRES efficiency varies according to cellular background, but also varies, for a same cell type, among individuals. The control of IRES activity presents an inherited component (h2) of 0.47 and 0.36 for EMCV and XIAP IRES, respectively. A genome scan identified a suggestive Quantitative Trait Loci (LOD 2.35) involved in the control of XIAP IRES activity. Résumé : Les sites internes d'entrée des ribosomes (IRES = internal ribosome entry sites) sont utilisés par les virus comme une stratégie afin d'outrepasser l'inhibition de traduction qui résulte communément d'une infection virale. Les IRES sont également utilisés par les cellules eucaryotes pour contrôler la traduction de l'ARN messager dans des conditions de stress cellulaire (apoptose, choc thermique) ou durant la phase G2 du cycle cellulaire, situations durant lesquelles la synthèse générale des protéines est inhibée. La variation des niveaux d'expression cellulaire de transcription est un caractère héréditaire. L'expression des gènes est contrôlée entre autre par les facteurs de transcription et par l'efficacité de la traduction initiée par la coiffe ainsi que par l'activité des ribosomes. Durant cette étude nous avons eu pour but d'identifier les déterminants génomiques responsables de la variabilité de la traduction contrôlée par l'IRES. Ceci a été effectué en étudiant deux IRES représentatifs : l'IRES du virus de l'encéphalomyocardite (EMCV) et l'IRES de l'inhibiteur de l'apoptose XIAP (X-linked Inhibitor-of-Apoptosis). Nous avons utilisés des lentivirus délivrant un transgène bicistronique codant pour deux gènes rapporteurs fluorescents. Ces lentivirus ont été utilisés pour transduire sept différentes lignées cellulaires de laboratoire et des lignées cellulaires lymphoblastoïdes B du Centre d'Etude du Polymorphisme Humain (CEPH; 15 pedigrees; n=209) qui représentent une approche in vitro de la structure familiale et qui permettent des analyses par balayage du génome. L'expression relative des deux marqueurs fluorescents a été analysée par FACS. Nos résultats montrent que l'efficacité des IRES varie en fonction du type de cellules. Il varie aussi, pour le même type de cellules, selon les individus. Le contrôle de l'activité de l'IRES est un caractère héritable (héritabilité h2) de 0.47 et 0.36 pour les IRES de EMCV et XIAP respectivement. Le balayage du génome a permis l'identification d'un locus à effets quantitatifs [QTL Quantitative Trait Loci (LOD 2.35)] impliqué dans le contôle de l'activité de l'IRES de XIAP.

A robust genomic signature for the detection of colorectal cancer patients with microsatellite instability phenotype and high mutation frequency.

Microsatellite instability (MSI) occurs in 10-20% of colorectal tumours and is associated with good prognosis. Here we describe the development and validation of a genomic signature that identifies colorectal cancer patients with MSI caused by DNA mismatch repair deficiency with high accuracy. Microsatellite status for 276 stage II and III colorectal tumours has been determined. Full-genome expression data was used to identify genes that correlate with MSI status. A subset of these samples (n = 73) had sequencing data for 615 genes available. An MSI gene signature of 64 genes was developed and validated in two independent validation sets: the first consisting of frozen samples from 132 stage II patients; and the second consisting of FFPE samples from the PETACC-3 trial (n = 625). The 64-gene MSI signature identified MSI patients in the first validation set with a sensitivity of 90.3% and an overall accuracy of 84.8%, with an AUC of 0.942 (95% CI, 0.888-0.975). In the second validation, the signature also showed excellent performance, with a sensitivity 94.3% and an overall accuracy of 90.6%, with an AUC of 0.965 (95% CI, 0.943-0.988). Besides correct identification of MSI patients, the gene signature identified a group of MSI-like patients that were MSS by standard assessment but MSI by signature assessment. The MSI-signature could be linked to a deficient MMR phenotype, as both MSI and MSI-like patients showed a high mutation frequency (8.2% and 6.4% of 615 genes assayed, respectively) as compared to patients classified as MSS (1.6% mutation frequency). The MSI signature showed prognostic power in stage II patients (n = 215) with a hazard ratio of 0.252 (p = 0.0145). Patients with an MSI-like phenotype had also an improved survival when compared to MSS patients. The MSI signature was translated to a diagnostic microarray and technically and clinically validated in FFPE and frozen samples.