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.

Draft genome of the globally widespread and invasive Argentine ant (Linepithema humile).

Ants are some of the most abundant and familiar animals on Earth, and they play vital roles in most terrestrial ecosystems. Although all ants are eusocial, and display a variety of complex and fascinating behaviors, few genomic resources exist for them. Here, we report the draft genome sequence of a particularly widespread and well-studied species, the invasive Argentine ant (Linepithema humile), which was accomplished using a combination of 454 (Roche) and Illumina sequencing and community-based funding rather than federal grant support. Manual annotation of >1,000 genes from a variety of different gene families and functional classes reveals unique features of the Argentine ant's biology, as well as similarities to Apis mellifera and Nasonia vitripennis. Distinctive features of the Argentine ant genome include remarkable expansions of gustatory (116 genes) and odorant receptors (367 genes), an abundance of cytochrome P450 genes (>110), lineage-specific expansions of yellow/major royal jelly proteins and desaturases, and complete CpG DNA methylation and RNAi toolkits. The Argentine ant genome contains fewer immune genes than Drosophila and Tribolium, which may reflect the prominent role played by behavioral and chemical suppression of pathogens. Analysis of the ratio of observed to expected CpG nucleotides for genes in the reproductive development and apoptosis pathways suggests higher levels of methylation than in the genome overall. The resources provided by this genome sequence will offer an abundance of tools for researchers seeking to illuminate the fascinating biology of this emerging model organism.

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

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.

Novel microsatellite loci for Sebaea aurea (Gentianaceae) and cross-amplification in related species.

Premise of the study: Microsatellite loci were developed in Sebaea aurea (Gentianaceae) to investigate the functional role of diplostigmaty (i.e., the presence of additional stigmas along the style). Methods and Results: One hundred seventy-four and 180 microsatellite loci were isolated through 454 shotgun sequencing of genomic and microsatellite-enriched DNA libraries, respectively. Sixteen polymorphic microsatellite loci were characterized, and 12 of them were selected to genotype individuals from two populations. Microsatellite amplification was conducted in two multiplex groups, each containing six microsatellite loci. Cross-species amplification was tested in seven other species of Sebaea. The 12 novel microsatellite loci amplified only in the two most closely related species to S. aurea (i.e., S. ambigua and S. minutiflora) and were also polymorphic in these two species. Conclusions: These results demonstrate the usefulness of this set of newly developed microsatellite loci to investigate the mating system and population genetic structure in S. aurea and related species.

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.

Positional scanning-synthetic peptide library-based analysis of self- and pathogen-derived peptide cross-reactivity with tumor-reactive Melan-A-specific CTL.

Synthetic combinatorial peptide libraries in positional scanning format (PS-SCL) have recently emerged as a useful tool for the analysis of T cell recognition. This includes identification of potentially cross-reactive sequences of self or pathogen origin that could be relevant for the understanding of TCR repertoire selection and maintenance, as well as of the cross-reactive potential of Ag-specific immune responses. In this study, we have analyzed the recognition of sequences retrieved by using a biometric analysis of the data generated by screening a PS-SCL with a tumor-reactive CTL clone specific for an immunodominant peptide from the melanocyte differentiation and tumor-associated Ag Melan-A. We found that 39% of the retrieved peptides were recognized by the CTL clone used for PS-SCL screening. The proportion of peptides recognized was higher among those with both high predicted affinity for the HLA-A2 molecule and high predicted stimulatory score. Interestingly, up to 94% of the retrieved peptides were cross-recognized by other Melan-A-specific CTL. Cross-recognition was at least partially focused, as some peptides were cross-recognized by the majority of CTL. Importantly, stimulation of PBMC from melanoma patients with the most frequently recognized peptides elicited the expansion of heterogeneous CD8(+) T cell populations, one fraction of which cross-recognized Melan-A. Together, these results underline the high predictive value of PS-SCL for the identification of sequences cross-recognized by Ag-specific T cells.

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.

Genetics of cadmium hyperaccumulation variance in Thlaspi caerulescens natural populations from Switzerland

SUMMARY Heavy metal presence in the environment is a serious concern since some of them can be toxic to plants, animals and humans once accumulated along the food chain. Cadmium (Cd) is one of the most toxic heavy metal. It is naturally present in soils at various levels and its concentration can be increased by human activities. Several plants however have naturally developed strategies allowing them to grow on heavy metal enriched soils. One of them consists in the accumulation and sequestration of heavy metals in the above-ground biomass. Some plants present in addition an extreme strategy by which they accumulate a limited number of heavy metals in their shoots in amounts 100 times superior to those expected for a non-accumulating plant in the same conditions. Understanding the genetic basis of the hyperaccumulation trait - particularly for Cd - remains an important challenge which may lead to biotechnological applications in the soil phytoremediation. In this thesis, Thlaspi caerulescens J. & C. Presl (Brassicaceae) was used as a model plant to study the Cd hyperaccumulation trait, owing to its physiological and genetic characteristics. Twenty-four wild populations were sampled in different regions of Switzerland. They were characterized for environmental and soil parameters as well as intrinsic characteristics of plants (i.e. metal concentrations in shoots). They were as well genetically characterized by AFLPs, plastid DNA polymorphism and genes markers (CAPS and microsatellites) mainly developed in this thesis. Some of the investigated genes were putatively linked to the Cd hyperaccumulation trait. Since the study of the Cd hyperaccumulation in the field is important as it allows the identification of patterns of selection, the present work offered a methodology to define the Cd hyperaccumulation capacity of populations from different habitats permitting thus their comparison in the field. We showed that Cd, Zn, Fe and Cu accumulations were linked and that populations with higher Cd hyperaccumulation capacity had higher shoot and reproductive fitness. Using our genetic data, statistical methods (Beaumont & Nichols's procedure, partial Mantel tests) were applied to identify genomic signatures of natural selection related to the Cd hyperaccumulation capacity. A significant genetic difference between populations related to their Cd hyperaccumulation capacity was revealed based on somè specific markers (AFLP and candidate genes). Polymorphism at the gene encoding IRTl (Iron-transporter also participating to the transport of Zn) was suggested as explaining part of the variation in Cd hyperaccumulation capacity of populations supporting previous physiological investigations. RÉSUMÉ La présence de métaux lourds dans l'environnement est un phénomène préoccupant. En effet, certains métaux lourds - comme le cadmium (Cd) -sont toxiques pour les plantes, les animaux et enfin, accumulés le long de la chaîne alimentaire, pour les hommes. Le Cd est naturellement présent dans le sol et sa concentration peut être accrue par différentes activités humaines. Certaines plantes ont cependant développé des stratégies leur permettant de pousser sur des sols contaminés en métaux lourds. Parmi elles, certaines accumulent et séquestrent les métaux lourds dans leurs parties aériennes. D`autres présentent une stratégie encore plus extrême. Elles accumulent un nombre limité de métaux lourds en quantités 100 fois supérieures à celles attendues pour des espèces non-accumulatrices sous de mêmes conditions. La compréhension des bases génétiques de l'hyperaccumulation -particulièrement celle du Cd - représente un défi important avec des applications concrètes en biotechnologies, tout particulièrement dans le but appliqué de la phytoremediation des sols contaminés. Dans cette thèse, Thlaspi caerulescens J. & C. Presl (Brassicaceae) a été utilisé comme modèle pour l'étude de l'hyperaccumulation du Cd de par ses caractéristiques physiologiques et génétiques. Vingt-quatre populations naturelles ont été échantillonnées en Suisse et pour chacune d'elles les paramètres environnementaux, pédologique et les caractéristiques intrinsèques aux plantes (concentrations en métaux lourds) ont été déterminés. Les populations ont été caractérisées génétiquement par des AFLP, des marqueurs chloroplastiques et des marqueurs de gènes spécifiques, particulièrement ceux potentiellement liés à l'hyperaccumulation du Cd (CAPS et microsatellites). La plupart ont été développés au cours de cette thèse. L'étude de l'hyperaccumulation du Cd en conditions naturelles est importante car elle permet d'identifier la marque, éventuelle de sélection naturelle. Ce travail offre ainsi une méthodologie pour définir et comparer la capacité des populations à hyperaccumuler le Cd dans différents habitats. Nous avons montré que les accumulations du Cd, Zn, Fe et Cu sont liées et que les populations ayant une grande capacité d'hyperaccumuler le Cd ont également une meilleure fitness végétative et reproductive. Des méthodes statistiques (l'approche de Beaumont & Nichols, tests de Martel partiels) ont été utilisées sur les données génétiques pour identifier la signature génomique de la sélection naturelle liée à la capacité d'hyperaccumuler le Cd. Une différenciation génétique des populations liée à leur capacité d'hyperaccumuler le Cd a été mise en évidence sur certains marqueurs spécifiques. En accord avec les études physiologiques connues, le polymorphisme au gène codant IRT1 (un transporteur de Fe impliqué dans le transport du Zn) pourrait expliquer une partie de la variance de la capacité des populations à hyperaccumuler le Cd.

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.