Loss of single HLA class I allospecificities in melanoma cells due to selective genomic abbreviations.

Expression of human leucocyte antigen (HLA) Class I molecules is essential for the recognition of malignant melanoma (MM) cells by CD8(+) T lymphocytes. A complete or partial loss of HLA Class I molecules is a potent strategy for MM cells to escape from immunosurveillance. In 2 out of 55 melanoma cell cultures we identified a complete phenotypic loss of HLA allospecificities. Both patients have been treated unsuccessfully with HLA-A2 peptides. To identify the reasons underlying the loss of single HLA-A allospecificities, we searched for genomic alterations at the locus for HLA Class I alpha-chain on chromosome 6 in melanoma cell cultures established from 2 selected patients with MM in advanced stage. This deficiency was associated with alterations of HLA-A2 gene sequences as determined by polymerase chain reaction-sequence specific primers (PCR-SSP). Karyotyping revealed a chromosomal loss in Patient 1, whereas melanoma cell cultures established from Patient 2 displayed 2 copies of chromosome 6. Loss of heterozygosity (LOH) using markers located around position 6p21 was detected in both cases. By applying group-specific primer-mixes spanning the 5'-flanking region of the HLA-A2 gene locus the relevant region was amplified by PCR and subsequent sequencing allowed alignment with the known HLA Class I reference sequences. Functional assays using HLA-A2-restricted cytotoxic T-cell clones were performed in HLA-A2 deficient MM cultures and revealed a drastically reduced susceptibility to CTL lysis in HLA-A2 negative cells. We could document the occurrence of selective HLA-A2 deficiencies in cultured advanced-stage melanoma metastases and identify their molecular causes as genomic alterations within the HLA-A gene locus.

DIP-STR: Highly Sensitive Markers for the Analysis of Unbalanced Genomic Mixtures.

Samples containing highly unbalanced DNA mixtures from two individuals commonly occur both in forensic mixed stains and in peripheral blood DNA microchimerism induced by pregnancy or following organ transplant. Because of PCR amplification bias, the genetic identification of a DNA that contributes trace amounts to a mixed sample represents a tremendous challenge. This means that standard genetic markers, namely microsatellites, also referred as short tandem repeats (STR), and single-nucleotide polymorphism (SNP) have limited power in addressing common questions of forensic and medical genetics. To address this issue, we developed a molecular marker, named DIP-STR that relies on pairing deletion-insertion polymorphisms (DIP) with STR. This novel analytical approach allows for the unambiguous genotyping of a minor component in the presence of a major component, where DIP-STR genotypes of the minor were successfully procured at ratios up to 1:1,000. The compound nature of this marker generates a high level of polymorphism that is suitable for identity testing. Here, we demonstrate the power of the DIP-STR approach on an initial set of nine markers surveyed in a Swiss population. Finally, we discuss the limitations and potential applications of our new system including preliminary tests on clinical samples and estimates of their performance on simulated DNA mixtures.

Evolution of a cellular immune response in Drosophila: a phenotypic and genomic comparative analysis.

Understanding the genomic basis of evolutionary adaptation requires insight into the molecular basis underlying phenotypic variation. However, even changes in molecular pathways associated with extreme variation, gains and losses of specific phenotypes, remain largely uncharacterized. Here, we investigate the large interspecific differences in the ability to survive infection by parasitoids across 11 Drosophila species and identify genomic changes associated with gains and losses of parasitoid resistance. We show that a cellular immune defense, encapsulation, and the production of a specialized blood cell, lamellocytes, are restricted to a sublineage of Drosophila, but that encapsulation is absent in one species of this sublineage, Drosophila sechellia. Our comparative analyses of hemopoiesis pathway genes and of genes differentially expressed during the encapsulation response revealed that hemopoiesis-associated genes are highly conserved and present in all species independently of their resistance. In contrast, 11 genes that are differentially expressed during the response to parasitoids are novel genes, specific to the Drosophila sublineage capable of lamellocyte-mediated encapsulation. These novel genes, which are predominantly expressed in hemocytes, arose via duplications, whereby five of them also showed signatures of positive selection, as expected if they were recruited for new functions. Three of these novel genes further showed large-scale and presumably loss-of-function sequence changes in D. sechellia, consistent with the loss of resistance in this species. In combination, these convergent lines of evidence suggest that co-option of duplicated genes in existing pathways and subsequent neofunctionalization are likely to have contributed to the evolution of the lamellocyte-mediated encapsulation in Drosophila.

Post-genomic update on a classical candidate gene for coronary artery disease: ESR1.

BACKGROUND: After age, sex is the most important risk factor for coronary artery disease (CAD). The mechanism through which women are protected from CAD is still largely unknown, but the observed sex difference suggests the involvement of the reproductive steroid hormone signaling system. Genetic association studies of the gene-encoding Estrogen Receptor α (ESR1) have shown conflicting results, although only a limited range of variation in the gene has been investigated. METHODS AND RESULTS: We exploited information made available by advanced new methods and resources in complex disease genetics to revisit the question of ESR1's role in risk of CAD. We performed a meta-analysis of 14 genome-wide association studies (CARDIoGRAM discovery analysis, N=≈87,000) to search for population-wide and sex-specific associations between CAD risk and common genetic variants throughout the coding, noncoding, and flanking regions of ESR1. In addition to samples from the MIGen (N=≈6000), WTCCC (N=≈7400), and Framingham (N=≈3700) studies, we extended this search to a larger number of common and uncommon variants by imputation into a panel of haplotypes constructed using data from the 1000 Genomes Project. Despite the widespread expression of ERα in vascular tissues, we found no evidence for involvement of common or low-frequency genetic variation throughout the ESR1 gene in modifying risk of CAD, either in the general population or as a function of sex. CONCLUSIONS: We suggest that future research on the genetic basis of sex-related differences in CAD risk should initially prioritize other genes in the reproductive steroid hormone biosynthesis system.

Host and invader impact of transfer of the clc genomic island into Pseudomonas aeruginosa PAO1

Genomic islands, large potentially mobile regions of bacterial chromosomes, are a major contributor to bacteria evolution. Here, we investigated the fitness cost and phenotypic differences between the bacterium Pseudomonas aeruginosa PAO1 and a derivative carrying one integrated copy of the clc element, a 103-kb genomic island [and integrative and conjugative element (ICE)] originating in Pseudomonas sp. strain B13 and a close relative of genomic islands found in clinical and environmental isolates of P. aeruginosa. By using a combination of whole genome transcriptome profiling, phenotypic arrays, competition experiments, and biofilm formation studies, only few differences became apparent, such as reduced biofilm growth and fourfold stationary phase repression of genes involved in acetoin metabolism in PAO1 containing the clc element. In contrast, PAO1 carrying the clc element acquired the capacity to grow on 3-chlorobenzoate and 2-aminophenol as sole carbon and energy substrates. No fitness loss >1% was detectable in competition experiments between PAO1 and PAO1 carrying the clc element. The genes from the clc element were not silent in PAO1, and excision was observed, although transfer of clc from PAO1 to other recipient bacteria was reduced by two orders of magnitude. Our results indicate that newly acquired mobile DNA not necessarily invoke an important fitness cost on their host. Absence of immediate detriment to the host may have contributed to the wide distribution of genomic islands like clc in bacterial genomes

Defining the role of common variation in the genomic and biological architecture of adult human height.

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Multiple adenosine 3',5'-cyclic [corrected] monophosphate response element DNA-binding proteins generated by gene diversification and alternative exon splicing.

The protein sequence deduced from the open reading frame of a human placental cDNA encoding a cAMP-responsive enhancer (CRE)-binding protein (CREB-327) has structural features characteristic of several other transcriptional transactivator proteins including jun, fos, C/EBP, myc, and CRE-BP1. Results of Southwestern analysis of nuclear extracts from several different cell lines show that there are multiple CRE-binding proteins, which vary in size in cell lines derived from different tissues and animal species. To examine the molecular diversity of CREB-327 and related proteins at the nucleic acid level, we used labeled cDNAs from human placenta that encode two different CRE-binding proteins (CREB-327 and CRE-BP1) to probe Northern and Southern blots. Both probes hybridized to multiple fragments on Southern blots of genomic DNA from various species. Alternatively, when a human placental c-jun probe was hybridized to the same blot, a single fragment was detected in most cases, consistent with the intronless nature of the human c-jun gene. The CREB-327 probe hybridized to multiple mRNAs, derived from human placenta, ranging in size from 2-9 kilobases. In contrast, the CRE-BP1 probe identified a single 4-kilobase mRNA. Sequence analyses of several overlapping human genomic cosmid clones containing CREB-327 sequences in conjunction with polymerase chain reaction indicates that the CREB-327/341 cDNAs are composed of at least eight or nine exons, and analyses of human placental cDNAs provide direct evidence for at least one alternatively spliced exon. Analyses of mouse/hamster-human hybridoma DNAs by Southern blotting and polymerase chain reaction localizes the CREB-327/341 gene to human chromosome 2. The results indicate that there is a dichotomy of CREB-like proteins, those that are related by overall structure and DNA-binding specificity as well as those that are related by close similarities of primary sequences.

Genomic organization, fine-mapping, and expression of the human islet-brain 1 (IB1)/c-Jun-amino-terminal kinase interacting protein-1 (JIP-1) gene.

Islet-brain 1 (IB1), a regulator of the pancreatic beta-cell function in the rat, is homologous to JIP-1, a murine inhibitor of c-Jun amino-terminal kinase (JNK). Whether IB1 and JIP-1 are present in humans was not known. We report the sequence of the 2133-bp human IB1 cDNA, the expression, structure, and fine-mapping of the human IB1 gene, and the characterization of an IB1 pseudogene. Human IB1 is 94% identical to rat IB1. The tissue-specific expression of IB1 in human is similar to that observed in rodent. The IB1 gene contains 12 exons and maps to chromosome 11 (11p11.2-p12), a region that is deleted in DEFECT-11 syndrome. Apart from an IB1 pseudogene on chromosome 17 (17q21), no additional IB1-related gene was found in the human genome. Our data indicate that the sequence and expression pattern of IB1 are highly conserved between rodent and human and provide the necessary tools to investigate whether IB1 is involved in human diseases.

An assessment of donor-to-recipient transmission patterns of human cytomegalovirus by analysis of viral genomic variants.

BACKGROUND: We studied human cytomegalovirus (CMV) donor-to-recipient transmission patterns in organ transplantation by analyzing genomic variants on the basis of CMV glycoprotein B (gB) genotyping. METHODS: Organ transplant recipients were included in the study if they had CMV viremia, if they had received an organ from a CMV-seropositive donor, and if there was at least 1 other recipient of an organ from the same donor who developed CMV viremia. Genotypes (gB1-4) were determined by real-time polymerase chain reaction. RESULTS: Forty-seven recipients of organs from 21 donors developed CMV viremia. Twenty-three recipients had a pretransplant donor/recipient (D/R) CMV serostatus of D(+)/R(+), and 24 had a serostatus of D(+)/R(-). The prevalences of genotypes in recipients were as follows: for gB1, 51% (n = 24); for gB2, 19% (n = 9); for gB3, 9% (n = 4); for gB4, 0% (n = 0); and for mixed infection, 21% (n = 10). Recipients of an organ from a common donor had infection with CMV of the same gB genotype in 12 (57%) of 21 instances. Concordance between genotypes was higher among seronegative (i.e., D(+)/R(-)) recipients than among seropositive (D(+)/R(+)) recipients, although discordances resulting from the transmission of multiple strains were seen. In seropositive recipients, transmission of multiple strains from the donor could not be differentiated from reactivation of a recipient's own strains. CONCLUSION: Our analysis of strain concordance among recipients of organs from common donors showed that transmission of CMV has complex dynamic patterns. In seropositive recipients, transmission or reactivation of multiple CMV strains is possible.

Genomic and proteomic approaches towards an understanding of sleep.

The basic functions of sleep are still unclear, however, recent advances in genomics and proteomics have begun to contribute to our understanding of both normal and pathological sleep. In this review, we focus primarily on normal sleep and wake that have been studied in model organisms such as mice. Mice have been especially valuable since many different inbred strains exist that differ in sleep-related traits, and genes can be altered by either mutagenesis or targeted approaches. Advances in QTL (Quantitative Trait Loci) analysis have also helped to identify important sleep related genes, and several other QTLs have been mapped as a first step toward finding the genes that underlie basic sleep traits. In addition to more traditional genetic approaches, the abundance of different mRNAs across sleep and wake can now be studied and compared in different brain regions much more thoroughly using microarray methods. Progress at the protein level has been more difficult, but a few studies have begun to investigate changes in proteins during sleep and wake, and we present some of our own preliminary data in this area. A knowledge of which genes and proteins control or respond to changes in sleep will not only help answer fundamental questions, but may also suggest novel drug targets for improving multiple aspects of sleep and wake.

Intracellular excision and reintegration dynamics of the ICEclc genomic island of Pseudomonas knackmussii sp. strain B13.

Genomic islands are DNA elements acquired by horizontal gene transfer that are common to a large number of bacterial genomes, which can contribute specific adaptive functions, e.g. virulence, metabolic capacities or antibiotic resistances. Some genomic islands are still self-transferable and display an intricate life-style, reminiscent of both bacteriophages and conjugative plasmids. Here we studied the dynamical process of genomic island excision and intracellular reintegration using the integrative and conjugative element ICEclc from Pseudomonas knackmussii B13 as model. By using self-transfer of ICEclc from strain B13 to Pseudomonas putida and Cupriavidus necator as recipients, we show that ICEclc can target a number of different tRNA(Gly) genes in a bacterial genome, but only those which carry the GCC anticodon. Two conditional traps were designed for ICEclc based on the attR sequence, and we could show that ICEclc will insert with different frequencies in such traps producing brightly fluorescent cells. Starting from clonal primary transconjugants we demonstrate that ICEclc is excising and reintegrating at detectable frequencies, even in the absence of recipient. Recombination site analysis provided evidence to explain the characteristics of a larger number of genomic island insertions observed in a variety of strains, including Bordetella petri, Pseudomonas aeruginosa and Burkholderia.

Genomic acute myeloid leukemia-associated inv(16)(p13q22) breakpoints are tightly clustered.

The inv(16) and related t(16;16) are found in 10% of all cases with de novo acute myeloid leukemia. In these rearrangements the core binding factor beta (CBFB) gene on 16q22 is fused to the smooth muscle myosin heavy chain gene (MYH11) on 16p13. To gain insight into the mechanisms causing the inv(16) we have analysed 24 genomic CBFB-MYH11 breakpoints. All breakpoints in CBFB are located in a 15-Kb intron. More than 50% of the sequenced 6.2 Kb of this intron consists of human repetitive elements. Twenty-one of the 24 breakpoints in MYH11 are located in a 370-bp intron. The remaining three breakpoints in MYH11 are located more upstream. The localization of three breakpoints adjacent to a V(D)J recombinase signal sequence in MYH11 suggests a V(D)J recombinase-mediated rearrangement in these cases. V(D)J recombinase-associated characteristics (small nucleotide deletions and insertions of random nucleotides) were detected in six other cases. CBFB and MYH11 duplications were detected in four of six cases tested.

First insights into the transcriptome and development of new genomic tools of a widespread circum-Mediterranean tree species, Pinus halepensis Mill.

Aleppo pine (Pinus halepensis Mill.) is a relevant conifer species for studying adaptive responses to drought and fire regimes in the Mediterranean region. In this study, we performed Illumina next-generation sequencing of two phenotypically divergent Aleppo pine accessions with the aims of (i) characterizing the transcriptome through Illumina RNA-Seq on trees phenotypically divergent for adaptive traits linked to fire adaptation and drought, (ii) performing a functional annotation of the assembled transcriptome, (iii) identifying genes with accelerated evolutionary rates, (iv) studying the expression levels of the annotated genes and (v) developing gene-based markers for population genomic and association genetic studies. The assembled transcriptome consisted of 48,629 contigs and covered about 54.6 Mbp. The comparison of Aleppo pine transcripts to Picea sitchensis protein-coding sequences resulted in the detection of 34,014 SNPs across species, with a Ka /Ks average value of 0.216, suggesting that the majority of the assembled genes are under negative selection. Several genes were differentially expressed across the two pine accessions with contrasted phenotypes, including a glutathione-s-transferase, a cellulose synthase and a cobra-like protein. A large number of new markers (3334 amplifiable SSRs and 28,236 SNPs) have been identified which should facilitate future population genomics and association genetics in this species. A 384-SNP Oligo Pool Assay for genotyping with the Illumina VeraCode technology has been designed which showed an high overall SNP conversion rate (76.6%). Our results showed that Illumina next-generation sequencing is a valuable technology to obtain an extensive overview on whole transcriptomes of nonmodel species with large genomes.

Impact of genomic methylation on radiation sensitivity of colorectal carcinoma.

PURPOSE: To investigate the influence of demethylation with 5-aza-cytidine (AZA) on radiation sensitivity and to define the intrinsic radiation sensitivity of methylation deficient colorectal carcinoma cells. METHODS AND MATERIALS: Radiation sensitizing effects of AZA were investigated in four colorectal carcinoma cell lines (HCT116, SW480, L174 T, Co115), defining influence of AZA on proliferation, clonogenic survival, and cell cycling with or without ionizing radiation. The methylation status for cancer or DNA damage response-related genes silenced by promoter methylation was determined. The effect of deletion of the potential target genes (DNMT1, DNMT3b, and double mutants) on radiation sensitivity was analyzed. RESULTS: AZA showed radiation sensitizing properties at >or=1 micromol/l, a concentration that does not interfere with the cell cycle by itself, in all four tested cell lines with a sensitivity-enhancing ratio (SER) of 1.6 to 2.1 (confidence interval [CI] 0.9-3.3). AZA successfully demethylated promoters of p16 and hMLH1, genes associated with ionizing radiation response. Prolonged exposure to low-dose AZA resulted in sustained radiosensitivity if associated with persistent genomic hypomethylation after recovery from AZA. Compared with maternal HCT116 cells, DNMT3b-defcient deficient cells were more sensitive to radiation with a SER of 2.0 (CI 0.9-2.1; p = 0.03), and DNMT3b/DNMT1-/- double-deficient cells showed a SER of 1.6 (CI 0.5-2.7; p = 0.09). CONCLUSIONS: AZA-induced genomic hypomethylation results in enhanced radiation sensitivity in colorectal carcinoma. The mediators leading to sensitization remain unknown. Defining the specific factors associated with radiation sensitization after genomic demethylation may open the way to better targeting for the purpose of radiation sensitization.

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone.

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long-standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well-studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best-supported scenario was a long period of allopatric isolation during the first three-quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.