Transcriptional differentiation of Trypanosoma brucei during in vitro acquisition of resistance to acoziborole.

Subspecies of the protozoan parasite Trypanosoma brucei are the causative agents of Human African Trypanosomiasis (HAT), a debilitating neglected tropical disease prevalent across sub-Saharan Africa. HAT case numbers have steadily decreased since the start of the century, and sustainable elimination of one form of the disease is in sight. However, key to this is the development of novel drugs to combat the disease. Acoziborole is a recently developed benzoxaborole, currently in advanced clinical trials, for treatment of stage 1 and stage 2 HAT. Importantly, acoziborole is orally bioavailable, and curative with one dose. Recent studies have made significant progress in determining the molecular mode of action of acoziborole. However, less is known about the potential mechanisms leading to acoziborole resistance in trypanosomes. In this study, an in vitro-derived acoziborole-resistant cell line was generated and characterised. The AcoR line exhibited significant cross-resistance with the methyltransferase inhibitor sinefungin as well as hypersensitisation to known trypanocides. Interestingly, transcriptomics analysis of AcoR cells indicated the parasites had obtained a procyclic- or stumpy-like transcriptome profile, with upregulation of procyclin surface proteins as well as differential regulation of key metabolic genes known to be expressed in a life cycle-specific manner, even in the absence of major morphological changes. However, no changes were observed in transcripts encoding CPSF3, the recently identified protein target of acoziborole. The results suggest that generation of resistance to this novel compound in vitro can be accompanied by transcriptomic switches resembling a procyclic- or stumpy-type phenotype.

Defining the sediment prokaryotic communities of the Indian River Lagoon, FL, USA, an Estuary of National Significance.

The Indian River Lagoon, located on the east coast of Florida, USA, is an Estuary of National Significance and an important economic and ecological resource. The Indian River Lagoon faces several environmental pressures, including freshwater discharges through the St. Lucie Estuary; accumulation of anoxic, fine-grained, organic-rich sediment; and metal contamination from agriculture and marinas. Although the Indian River Lagoon has been well-studied, little is known about its microbial communities; thus, a two-year 16S amplicon sequencing study was conducted to assess the spatiotemporal changes of the sediment bacterial and archaeal groups. In general, the Indian River Lagoon exhibited a prokaryotic community that was consistent with other estuarine studies. Statistically different communities were found between the Indian River Lagoon and St. Lucie Estuary due to changes in porewater salinity causing microbes that require salts for growth to be higher in the Indian River Lagoon. The St. Lucie Estuary exhibited more obvious prokaryotic seasonality, such as a higher relative abundance of Betaproteobacteriales in wet season and a higher relative abundance of Flavobacteriales in dry season samples. Distance-based linear models revealed these communities were more affected by changes in total organic matter and copper than changes in temperature. Anaerobic prokaryotes, such as Campylobacterales, were more associated with high total organic matter and copper samples while aerobic prokaryotes, such as Nitrosopumilales, were more associated with low total organic matter and copper samples. This initial study fills the knowledge gap on the Indian River Lagoon bacterial and archaeal communities and serves as important data for future studies to compare to determine possible future changes due to human impacts or environmental changes.

Genomic instability at the locus of sterol C24-methyltransferase promotes amphotericin B resistance in Leishmania parasites.

Amphotericin B is an increasingly important tool in efforts to reduce the global disease burden posed by Leishmania parasites. With few other chemotherapeutic options available for the treatment of leishmaniasis, the potential for emergent resistance to this drug is a considerable threat. Here we characterised four novel amphotericin B-resistant Leishmania mexicana lines. All lines exhibited altered sterol biosynthesis, and hypersensitivity to pentamidine. Whole genome sequencing demonstrated resistance-associated mutation of the sterol biosynthesis gene sterol C5-desaturase in one line. However, in three out of four lines, RNA-seq revealed loss of expression of sterol C24-methyltransferase (SMT) responsible for drug resistance and altered sterol biosynthesis. Additional loss of the miltefosine transporter was associated with one of those lines. SMT is encoded by two tandem gene copies, which we found to have very different expression levels. In all cases, reduced overall expression was associated with loss of the 3' untranslated region of the dominant gene copy, resulting from structural variations at this locus. Local regions of sequence homology, between the gene copies themselves, and also due to the presence of SIDER1 retrotransposon elements that promote multi-gene amplification, correlate to these structural variations. Moreover, in at least one case loss of SMT expression was not associated with loss of virulence in primary macrophages or in vivo. Whilst such repeat sequence-mediated instability is known in Leishmania genomes, its presence associated with resistance to a major antileishmanial drug, with no evidence of associated fitness costs, is a significant concern.

A Leishmania infantum genetic marker associated with miltefosine treatment failure for visceral leishmaniasis.

BACKGROUND: Miltefosine has been used successfully to treat visceral leishmaniasis (VL) in India, but it was unsuccessful for VL in a clinical trial in Brazil. METHODS: To identify molecular markers that predict VL treatment failure whole genome sequencing of 26 L. infantum isolates, from cured and relapsed patients allowed a GWAS analysis of SNPs, gene and chromosome copy number variations. FINDINGS: A strong association was identified (p = 0·0005) between the presence of a genetically stable L. infantumMiltefosine Sensitivity Locus (MSL), and a positive response to miltefosine treatment. The risk of treatment failure increased 9·4-fold (95% CI 2·11-53·54) when an isolate did not have the MSL. The complete absence of the MSL predicted miltefosine failure with 0·92 (95% CI 0·65-0·996) sensitivity and 0·78 (95% CI 0·52-0·92) specificity. A genotyping survey of L. infantum (n = 157) showed that the frequency of MSL varies in a cline from 95% in North East Brazil to <5% in the South East. The MSL was found in the genomes of all L. infantum and L. donovani sequenced isolates from the Old World (n = 671), where miltefosine can have a cure rate higher than 93%. INTERPRETATION: Knowledge on the presence or absence of the MSL in L. infantum will allow stratification of patients prior to treatment, helping to establish better therapeutic strategies for VL treatment. FUND: CNPq, FAPES, GCRF MRC and Wellcome Trust.

Sterol 14α-demethylase mutation leads to amphotericin B resistance in Leishmania mexicana.

Amphotericin B has emerged as the therapy of choice for use against the leishmaniases. Administration of the drug in its liposomal formulation as a single injection is being promoted in a campaign to bring the leishmaniases under control. Understanding the risks and mechanisms of resistance is therefore of great importance. Here we select amphotericin B-resistant Leishmania mexicana parasites with relative ease. Metabolomic analysis demonstrated that ergosterol, the sterol known to bind the drug, is prevalent in wild-type cells, but diminished in the resistant line, where alternative sterols become prevalent. This indicates that the resistance phenotype is related to loss of drug binding. Comparing sequences of the parasites' genomes revealed a plethora of single nucleotide polymorphisms that distinguish wild-type and resistant cells, but only one of these was found to be homozygous and associated with a gene encoding an enzyme in the sterol biosynthetic pathway, sterol 14α-demethylase (CYP51). The mutation, N176I, is found outside of the enzyme's active site, consistent with the fact that the resistant line continues to produce the enzyme's product. Expression of wild-type sterol 14α-demethylase in the resistant cells caused reversion to drug sensitivity and a restoration of ergosterol synthesis, showing that the mutation is indeed responsible for resistance. The amphotericin B resistant parasites become hypersensitive to pentamidine and also agents that induce oxidative stress. This work reveals the power of combining polyomics approaches, to discover the mechanism underlying drug resistance as well as offering novel insights into the selection of resistance to amphotericin B itself.

Comparative genomics of Tunisian Leishmania major isolates causing human cutaneous leishmaniasis with contrasting clinical severity.

Zoonotic cutaneous leishmaniasis caused by Leishmania (L.) major parasites affects urban and suburban areas in the center and south of Tunisia where the disease is endemo-epidemic. Several cases were reported in human patients for which infection due to L. major induced lesions with a broad range of severity. However, very little is known about the mechanisms underlying this diversity. Our hypothesis is that parasite genomic variability could, in addition to the host immunological background, contribute to the intra-species clinical variability observed in patients and explain the lesion size differences observed in the experimental model. Based on several epidemiological, in vivo and in vitro experiments, we focused on two clinical isolates showing contrasted severity in patients and BALB/c experimental mice model. We used DNA-seq as a high-throughput technology to facilitate the identification of genetic variants with discriminating potential between both isolates. Our results demonstrate that various levels of heterogeneity could be found between both L. major isolates in terms of chromosome or gene copy number variation (CNV), and that the intra-species divergence could surprisingly be related to single nucleotide polymorphisms (SNPs) and Insertion/Deletion (InDels) events. Interestingly, we particularly focused here on genes affected by both types of variants and correlated them with the observed gene CNV. Whether these differences are sufficient to explain the severity in patients is obviously still open to debate, but we do believe that additional layers of -omic information is needed to complement the genomic screen in order to draw a more complete map of severity determinants.

Genome-wide mapping reveals single-origin chromosome replication in Leishmania, a eukaryotic microbe.

BACKGROUND: DNA replication initiates on defined genome sites, termed origins. Origin usage appears to follow common rules in the eukaryotic organisms examined to date: all chromosomes are replicated from multiple origins, which display variations in firing efficiency and are selected from a larger pool of potential origins. To ask if these features of DNA replication are true of all eukaryotes, we describe genome-wide origin mapping in the parasite Leishmania. RESULTS: Origin mapping in Leishmania suggests a striking divergence in origin usage relative to characterized eukaryotes, since each chromosome appears to be replicated from a single origin. By comparing two species of Leishmania, we find evidence that such origin singularity is maintained in the face of chromosome fusion or fission events during evolution. Mapping Leishmania origins suggests that all origins fire with equal efficiency, and that the genomic sites occupied by origins differ from related non-origins sites. Finally, we provide evidence that origin location in Leishmania displays striking conservation with Trypanosoma brucei, despite the latter parasite replicating its chromosomes from multiple, variable strength origins. CONCLUSIONS: The demonstration of chromosome replication for a single origin in Leishmania, a microbial eukaryote, has implications for the evolution of origin multiplicity and associated controls, and may explain the pervasive aneuploidy that characterizes Leishmania chromosome architecture.

The Comparative Genomics and Phylogenomics of Leishmania amazonensis Parasite.

Leishmaniasis is an infectious disease caused by Leishmania species. Leishmania amazonensis is a New World Leishmania species belonging to the Mexicana complex, which is able to cause all types of leishmaniasis infections. The L. amazonensis reference strain MHOM/BR/1973/M2269 was sequenced identifying 8,802 codifying sequences (CDS), most of them of hypothetical function. Comparative analysis using six Leishmania species showed a core set of 7,016 orthologs. L. amazonensis and Leishmania mexicana share the largest number of distinct orthologs, while Leishmania braziliensis presented the largest number of inparalogs. Additionally, phylogenomic analysis confirmed the taxonomic position for L. amazonensis within the "Mexicana complex", reinforcing understanding of the split of New and Old World Leishmania. Potential non-homologous isofunctional enzymes (NISE) were identified between L. amazonensis and Homo sapiens that could provide new drug targets for development.

A cascade of DNA-binding proteins for sexual commitment and development in Plasmodium.

Commitment to and completion of sexual development are essential for malaria parasites (protists of the genus Plasmodium) to be transmitted through mosquitoes. The molecular mechanism(s) responsible for commitment have been hitherto unknown. Here we show that PbAP2-G, a conserved member of the apicomplexan AP2 (ApiAP2) family of DNA-binding proteins, is essential for the commitment of asexually replicating forms to sexual development in Plasmodium berghei, a malaria parasite of rodents. PbAP2-G was identified from mutations in its encoding gene, PBANKA_143750, which account for the loss of sexual development frequently observed in parasites transmitted artificially by blood passage. Systematic gene deletion of conserved ApiAP2 genes in Plasmodium confirmed the role of PbAP2-G and revealed a second ApiAP2 member (PBANKA_103430, here termed PbAP2-G2) that significantly modulates but does not abolish gametocytogenesis, indicating that a cascade of ApiAP2 proteins are involved in commitment to the production and maturation of gametocytes. The data suggest a mechanism of commitment to gametocytogenesis in Plasmodium consistent with a positive feedback loop involving PbAP2-G that could be exploited to prevent the transmission of this pernicious parasite.

The streamlined genome of Phytomonas spp. relative to human pathogenic kinetoplastids reveals a parasite tailored for plants.

Members of the family Trypanosomatidae infect many organisms, including animals, plants and humans. Plant-infecting trypanosomes are grouped under the single genus Phytomonas, failing to reflect the wide biological and pathological diversity of these protists. While some Phytomonas spp. multiply in the latex of plants, or in fruit or seeds without apparent pathogenicity, others colonize the phloem sap and afflict plants of substantial economic value, including the coffee tree, coconut and oil palms. Plant trypanosomes have not been studied extensively at the genome level, a major gap in understanding and controlling pathogenesis. We describe the genome sequences of two plant trypanosomatids, one pathogenic isolate from a Guianan coconut and one non-symptomatic isolate from Euphorbia collected in France. Although these parasites have extremely distinct pathogenic impacts, very few genes are unique to either, with the vast majority of genes shared by both isolates. Significantly, both Phytomonas spp. genomes consist essentially of single copy genes for the bulk of their metabolic enzymes, whereas other trypanosomatids e.g. Leishmania and Trypanosoma possess multiple paralogous genes or families. Indeed, comparison with other trypanosomatid genomes revealed a highly streamlined genome, encoding for a minimized metabolic system while conserving the major pathways, and with retention of a full complement of endomembrane organelles, but with no evidence for functional complexity. Identification of the metabolic genes of Phytomonas provides opportunities for establishing in vitro culturing of these fastidious parasites and new tools for the control of agricultural plant disease.

Regulators of Trypanosoma brucei cell cycle progression and differentiation identified using a kinome-wide RNAi screen.

The African trypanosome, Trypanosoma brucei, maintains an integral link between cell cycle regulation and differentiation during its intricate life cycle. Whilst extensive changes in phosphorylation have been documented between the mammalian bloodstream form and the insect procyclic form, relatively little is known about the parasite's protein kinases (PKs) involved in the control of cellular proliferation and differentiation. To address this, a T. brucei kinome-wide RNAi cell line library was generated, allowing independent inducible knockdown of each of the parasite's 190 predicted protein kinases. Screening of this library using a cell viability assay identified ≥42 PKs that are required for normal bloodstream form proliferation in culture. A secondary screen identified 24 PKs whose RNAi-mediated depletion resulted in a variety of cell cycle defects including in G1/S, kinetoplast replication/segregation, mitosis and cytokinesis, 15 of which are novel cell cycle regulators. A further screen identified for the first time two PKs, named repressor of differentiation kinase (RDK1 and RDK2), depletion of which promoted bloodstream to procyclic form differentiation. RDK1 is a membrane-associated STE11-like PK, whilst RDK2 is a NEK PK that is essential for parasite proliferation. RDK1 acts in conjunction with the PTP1/PIP39 phosphatase cascade to block uncontrolled bloodstream to procyclic form differentiation, whilst RDK2 is a PK whose depletion efficiently induces differentiation in the absence of known triggers. Thus, the RNAi kinome library provides a valuable asset for functional analysis of cell signalling pathways in African trypanosomes as well as drug target identification and validation.

Chromosome and gene copy number variation allow major structural change between species and strains of Leishmania.

Leishmania parasites cause a spectrum of clinical pathology in humans ranging from disfiguring cutaneous lesions to fatal visceral leishmaniasis. We have generated a reference genome for Leishmania mexicana and refined the reference genomes for Leishmania major, Leishmania infantum, and Leishmania braziliensis. This has allowed the identification of a remarkably low number of genes or paralog groups (2, 14, 19, and 67, respectively) unique to one species. These were found to be conserved in additional isolates of the same species. We have predicted allelic variation and find that in these isolates, L. major and L. infantum have a surprisingly low number of predicted heterozygous SNPs compared with L. braziliensis and L. mexicana. We used short read coverage to infer ploidy and gene copy numbers, identifying large copy number variations between species, with 200 tandem gene arrays in L. major and 132 in L. mexicana. Chromosome copy number also varied significantly between species, with nine supernumerary chromosomes in L. infantum, four in L. mexicana, two in L. braziliensis, and one in L. major. A significant bias against gene arrays on supernumerary chromosomes was shown to exist, indicating that duplication events occur more frequently on disomic chromosomes. Taken together, our data demonstrate that there is little variation in unique gene content across Leishmania species, but large-scale genetic heterogeneity can result through gene amplification on disomic chromosomes and variation in chromosome number. Increased gene copy number due to chromosome amplification may contribute to alterations in gene expression in response to environmental conditions in the host, providing a genetic basis for disease tropism.

Genetic factors underlying the risk of thalidomide-related neuropathy in patients with multiple myeloma.

PURPOSE: To indentify genetic variation that can modulate and predict the risk of developing thalidomide-related peripheral neuropathy (TrPN). PATIENTS AND METHODS: We analyzed DNA from 1,495 patients with multiple myeloma. Using a custom-built single nucleotide polymorphism (SNP) array, we tested the association of TrPN with 3,404 SNPs. The SNPs were selected in predicted functional regions within 964 genes spanning 67 molecular pathways thought to be involved in the pathogenesis, treatment response, and adverse effects associated with myeloma and its therapy. Patient cases and controls were derived from two large clinical trials that compared thalidomide with conventional-based treatment in myeloma patients (Medical Research Council Myeloma-IX and HOVON-50/GMMG-HD3). RESULTS: We report TrPN associations with SNPs-ABCA1 (rs363717), ICAM1 (rs1799969), PPARD (rs2076169), SERPINB2 (rs6103), and SLC12A6 (rs7164902)-where we show cross validation of the associations in both trials. To investigate whether TrPN SNP associations were related to exposure to thalidomide only or general drug-related peripheral neuropathy, we performed a second analysis on patients treated with vincristine. We report SNPs associated with vincristine neuropathy, with a seemingly distinct underlying genetic mechanism. CONCLUSION: Our results are consistent with the hypothesis that an individual's risk of developing a peripheral neuropathy after thalidomide treatment can be mediated by polymorphisms in genes governing repair mechanisms and inflammation in the peripheral nervous system. These findings will contribute to the development of future neuroprotective strategies with thalidomide therapy and the better use of this important compound.

A compendium of myeloma-associated chromosomal copy number abnormalities and their prognostic value.

To obtain a comprehensive genomic profile of presenting multiple myeloma cases we performed high-resolution single nucleotide polymorphism mapping array analysis in 114 samples alongside 258 samples analyzed by U133 Plus 2.0 expression array (Affymetrix). We examined DNA copy number alterations and loss of heterozygosity (LOH) to define the spectrum of minimally deleted regions in which relevant genes of interest can be found. The most frequent deletions are located at 1p (30%), 6q (33%), 8p (25%), 12p (15%), 13q (59%), 14q (39%), 16q (35%), 17p (7%), 20 (12%), and 22 (18%). In addition, copy number-neutral LOH, or uniparental disomy, was also prevalent on 1q (8%), 16q (9%), and X (20%), and was associated with regions of gain and loss. Based on fluorescence in situ hybridization and expression quartile analysis, genes of prognostic importance were found to be located at 1p (FAF1, CDKN2C), 1q (ANP32E), and 17p (TP53). In addition, we identified common homozygously deleted genes that have functions relevant to myeloma biology. Taken together, these analyses indicate that the crucial pathways in myeloma pathogenesis include the nuclear factor-κB pathway, apoptosis, cell-cycle regulation, Wnt signaling, and histone modifications. This study was registered at http://isrctn.org as ISRCTN68454111.

XBP1s levels are implicated in the biology and outcome of myeloma mediating different clinical outcomes to thalidomide-based treatments.

Immunoglobulin production by myeloma plasma cells depends on the unfolded protein response for protein production and folding. Recent studies have highlighted the importance of IRE1alpha and X box binding protein 1 (XBP1), key members of this pathway, in normal B-plasma cell development. We have determined the gene expression levels of IRE1alpha, XBP1, XBP1UNSPLICED (XBP1u), and XBP1SPLICED (XBP1s) in a series of patients with myeloma and correlated findings with clinical outcome. We show that IRE1alpha and XBP1 are highly expressed and that patients with low XBP1s/u ratios have a significantly better overall survival. XBP1s is an independent prognostic marker and can be used with beta2 microglobulin and t(4;14) to identify a group of patients with a poor outcome. Furthermore, we show the beneficial therapeutic effects of thalidomide in patients with low XBP1s/u ratios. This study highlights the importance of XBP1 in myeloma and its significance as an independent prognostic marker and as a predictor of thalidomide response.

Homozygous deletion mapping in myeloma samples identifies genes and an expression signature relevant to pathogenesis and outcome.

PURPOSE: Myeloma is a clonal malignancy of plasma cells. Poor-prognosis risk is currently identified by clinical and cytogenetic features. However, these indicators do not capture all prognostic information. Gene expression analysis can be used to identify poor-prognosis patients and this can be improved by combination with information about DNA-level changes. EXPERIMENTAL DESIGN: Using single nucleotide polymorphism-based gene mapping in combination with global gene expression analysis, we have identified homozygous deletions in genes and networks that are relevant to myeloma pathogenesis and outcome. RESULTS: We identified 170 genes with homozygous deletions and corresponding loss of expression. Deletion within the "cell death" network was overrepresented and cases with these deletions had impaired overall survival. From further analysis of these events, we have generated an expression-based signature associated with shorter survival in 258 patients and confirmed this signature in data from two independent groups totaling 800 patients. We defined a gene expression signature of 97 cell death genes that reflects prognosis and confirmed this in two independent data sets. CONCLUSIONS: We developed a simple 6-gene expression signature from the 97-gene signature that can be used to identify poor-prognosis myeloma in the clinical environment. This signature could form the basis of future trials aimed at improving the outcome of poor-prognosis myeloma.

MMSET deregulation affects cell cycle progression and adhesion regulons in t(4;14) myeloma plasma cells.

BACKGROUND: The recurrent immunoglobulin translocation, t(4;14)(p16;q32) occurs in 15% of multiple myeloma patients and is associated with poor prognosis, through an unknown mechanism. The t(4;14) up-regulates fibroblast growth factor receptor 3 (FGFR3) and multiple myeloma SET domain (MMSET) genes. The involvement of MMSET in the pathogenesis of t(4;14) multiple myeloma and the mechanism or genes deregulated by MMSET upregulation are still unclear. DESIGN AND METHODS: The expression of MMSET was analyzed using a novel antibody. The involvement of MMSET in t(4;14) myelomagenesis was assessed by small interfering RNA mediated knockdown combined with several biological assays. In addition, the differential gene expression of MMSET-induced knockdown was analyzed with expression microarrays. MMSET gene targets in primary patient material was analyzed by expression microarrays. RESULTS: We found that MMSET isoforms are expressed in multiple myeloma cell lines, being exclusively up-regulated in t(4;14)-positive cells. Suppression of MMSET expression affected cell proliferation by both decreasing cell viability and cell cycle progression of cells with the t(4;14) translocation. These findings were associated with reduced expression of genes involved in the regulation of cell cycle progression (e.g. CCND2, CCNG1, BRCA1, AURKA and CHEK1), apoptosis (CASP1, CASP4 and FOXO3A) and cell adhesion (ADAM9 and DSG2). Furthermore, we identified genes involved in the latter processes that were differentially expressed in t(4;14) multiple myeloma patient samples. CONCLUSIONS: In conclusion, dysregulation of MMSET affects the expression of several genes involved in the regulation of cell cycle progression, cell adhesion and survival.

Genome-wide co-expression analysis in multiple tissues.

Expression quantitative trait loci (eQTLs) represent genetic control points of gene expression, and can be categorized as cis- and trans-acting, reflecting local and distant regulation of gene expression respectively. Although there is evidence of co-regulation within clusters of trans-eQTLs, the extent of co-expression patterns and their relationship with the genotypes at eQTLs are not fully understood. We have mapped thousands of cis- and trans-eQTLs in four tissues (fat, kidney, adrenal and left ventricle) in a large panel of rat recombinant inbred (RI) strains. Here we investigate the genome-wide correlation structure in expression levels of eQTL transcripts and underlying genotypes to elucidate the nature of co-regulation within cis- and trans-eQTL datasets. Across the four tissues, we consistently found statistically significant correlations of cis-regulated gene expression to be rare (<0.9% of all pairs tested). Most (>80%) of the observed significant correlations of cis-regulated gene expression are explained by correlation of the underlying genotypes. In comparison, co-expression of trans-regulated gene expression is more common, with significant correlation ranging from 2.9%-14.9% of all pairs of trans-eQTL transcripts. We observed a total of 81 trans-eQTL clusters (hot-spots), defined as consisting of > or =10 eQTLs linked to a common region, with very high levels of correlation between trans-regulated transcripts (77.2-90.2%). Moreover, functional analysis of large trans-eQTL clusters (> or =30 eQTLs) revealed significant functional enrichment among genes comprising 80% of the large clusters. The results of this genome-wide co-expression study show the effects of the eQTL genotypes on the observed patterns of correlation, and suggest that functional relatedness between genes underlying trans-eQTLs is reflected in the degree of co-expression observed in trans-eQTL clusters. Our results demonstrate the power of an integrative, systematic approach to the analysis of a large gene expression dataset to uncover underlying structure, and inform future eQTL studies.

Deletions of CDKN2C in multiple myeloma: biological and clinical implications.

PURPOSE: Deletions of chromosome 1 have been described in 7% to 40% of cases of myeloma with inconsistent clinical consequences. CDKN2C at 1p32.3 has been identified in myeloma cell lines as the potential target of the deletion. We tested the clinical impact of 1p deletion and used high-resolution techniques to define the role of CDKN2C in primary patient material. EXPERIMENTAL DESIGN: We analyzed 515 cases of monoclonal gammopathy of undetermined significance (MGUS), smoldering multiple myeloma (SMM), and newly diagnosed multiple myeloma using fluorescence in situ hybridization (FISH) for deletions of CDKN2C. In 78 myeloma cases, we carried out Affymetrix single nucleotide polymorphism mapping and U133 Plus 2.0 expression arrays. In addition, we did mutation, methylation, and Western blotting analysis. RESULTS: By FISH we identified deletion of 1p32.3 (CDKN2C) in 3 of 66 MGUS (4.5%), 4 of 39 SMM (10.3%), and 55 of 369 multiple myeloma cases (15%). We examined the impact of copy number change at CDKN2C on overall survival (OS), and found that the cases with either hemizygous or homozygous deletion of CDKN2C had a worse OS compared with cases that were intact at this region (22 months versus 38 months; P = 0.003). Using gene mapping we identified three homozygous deletions at 1p32.3, containing CDKN2C, all of which lacked expression of CDKN2C. Cases with homozygous deletions of CDKN2C were the most proliferative myelomas, defined by an expression-based proliferation index, consistent with its biological function as a cyclin-dependent kinase inhibitor. CONCLUSIONS: Our results suggest that deletions of CDKN2C are important in the progression and clinical outcome of myeloma.

Genetic associations with thalidomide mediated venous thrombotic events in myeloma identified using targeted genotyping.

A venous thromboembolism (VTE) with the subsequent risk of pulmonary embolism is a major concern in the treatment of patients with multiple myeloma with thalidomide. The susceptibility to developing a VTE in response to thalidomide therapy is likely to be influenced by both genetic and environmental factors. To test genetic variation associated with treatment related VTE in patient peripheral blood DNA, we used a custom-built molecular inversion probe (MIP)-based single nucleotide polymorphism (SNP) chip containing 3404 SNPs. SNPs on the chip were selected in "functional regions" within 964 genes spanning 67 molecular pathways thought to be involved in the pathogenesis, treatment response, and side effects associated with myeloma therapy. Patients and controls were taken from 3 large clinical trials: Medical Research Council (MRC) Myeloma IX, Hovon-50, and Eastern Cooperative Oncology Group (ECOG) EA100, which compared conventional treatments with thalidomide in patients with myeloma. Our analysis showed that the set of SNPs associated with thalidomide-related VTE were enriched in genes and pathways important in drug transport/metabolism, DNA repair, and cytokine balance. The effects of the SNPs associated with thalidomide-related VTE may be functional at the level of the tumor cell, the tumor-related microenvironment, and the endothelium. The clinical trials described in this paper have been registered as follows: MRC Myeloma IX: ISRCTN68454111; Hovon-50: NCT00028886; and ECOG EA100: NCT00033332.