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.

Mutation of yeast Ku genes disrupts the subnuclear organization of telomeres.

The mammalian Ku70 and Ku86 proteins form a heterodimer that binds to the ends of double-stranded DNA in vitro and is required for repair of radiation-induced strand breaks and V(D)J recombination [1,2]. Deletion of the Saccharomyces cerevisiae genes HDF1 and HDF2--encoding yKu70p and yKu80p, respectively--enhances radiation sensitivity in a rad52 background [3,4]. In addition to repair defects, the length of the TG-rich repeat on yeast telomere ends shortens dramatically [5,6]. We have shown previously that in yeast interphase nuclei, telomeres are clustered in a limited number of foci near the nuclear periphery [7], but the elements that mediate this localization remained unknown. We report here that deletion of the genes encoding yKu70p or its partner yKu80p altered the positioning of telomeric DNA in the yeast nucleus. These are the first mutants shown to affect the subnuclear localization of telomeres. Strains deficient for either yKu70p or yKu80p lost telomeric silencing, although they maintained repression at the silent mating-type loci. In addition, the telomere-associated silencing factors Sir3p and Sir4p and the TG-repeat-binding protein Rap1p lost their punctate pattern of staining and became dispersed throughout the nucleoplasm. Our results implicate the yeast Ku proteins directly in aspects of telomere organization, which in turn affects the repression of telomere-proximal genes.

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.

Organization of lin genes and IS6100 among different strains of hexachlorocyclohexane-degrading Sphingomonas paucimobilis: evidence for horizontal gene transfer.

The organization of lin genes and IS6100 was studied in three strains of Sphingomonas paucimobilis (B90A, Sp+, and UT26) which degraded hexachlorocyclohexane (HCH) isomers but which had been isolated at different geographical locations. DNA-DNA hybridization data revealed that most of the lin genes in these strains were associated with IS6100, an insertion sequence classified in the IS6 family and initially found in Mycobacterium fortuitum. Eleven, six, and five copies of IS6100 were detected in B90A, Sp+, and UT26, respectively. IS6100 elements in B90A were sequenced from five, one, and one regions of the genomes of B90A, Sp+, and UT26, respectively, and were found to be identical. DNA-DNA hybridization and DNA sequencing of cosmid clones also revealed that S. paucimobilis B90A contains three and two copies of linX and linA, respectively, compared to only one copy of these genes in strains Sp+ and UT26. Although the copy number and the sequence of the remaining genes of the HCH degradative pathway (linB, linC, linD, and linE) were nearly the same in all strains, there were striking differences in the organization of the linA genes as a result of replacement of portions of DNA sequences by IS6100, which gave them a strange mosaic configuration. Spontaneous deletion of linD and linE from B90A and of linA from Sp+ occurred and was associated either with deletion of a copy of IS6100 or changes in IS6100 profiles. The evidence gathered in this study, coupled with the observation that the G+C contents of the linA genes are lower than that of the remaining DNA sequence of S. paucimobilis, strongly suggests that all these strains acquired the linA gene through horizontal gene transfer mediated by IS6100. The association of IS6100 with the rest of the lin genes further suggests that IS6100 played a role in shaping the current lin gene organization.

Revised definitions of invasive fungal disease from the European Organization for Research and Treatment of Cancer/Invasive Fungal Infections Cooperative Group and the National Institute of Allergy and Infectious Diseases Mycoses Study Group (EORTC/MSG) Consensus Group.

BACKGROUND: Invasive fungal diseases are important causes of morbidity and mortality. Clarity and uniformity in defining these infections are important factors in improving the quality of clinical studies. A standard set of definitions strengthens the consistency and reproducibility of such studies. METHODS: After the introduction of the original European Organization for Research and Treatment of Cancer/Invasive Fungal Infections Cooperative Group and the National Institute of Allergy and Infectious Diseases Mycoses Study Group (EORTC/MSG) Consensus Group definitions, advances in diagnostic technology and the recognition of areas in need of improvement led to a revision of this document. The revision process started with a meeting of participants in 2003, to decide on the process and to draft the proposal. This was followed by several rounds of consultation until a final draft was approved in 2005. This was made available for 6 months to allow public comment, and then the manuscript was prepared and approved. RESULTS: The revised definitions retain the original classifications of "proven," "probable," and "possible" invasive fungal disease, but the definition of "probable" has been expanded, whereas the scope of the category "possible" has been diminished. The category of proven invasive fungal disease can apply to any patient, regardless of whether the patient is immunocompromised, whereas the probable and possible categories are proposed for immunocompromised patients only. CONCLUSIONS: These revised definitions of invasive fungal disease are intended to advance clinical and epidemiological research and may serve as a useful model for defining other infections in high-risk patients.

Complex organization of CTF/NF-I, C/EBP, and HNF3 binding sites within the promoter of the liver-specific vitellogenin gene.

Vitellogenin genes are expressed specifically in the liver of female oviparous vertebrates under the strict control of estrogen. To explain this tissue-specific expression, we performed a detailed analysis of the Xenopus laevis vitellogenin gene B1 promoter by DNase I footprinting and gel mobility-shift assays. We characterized five binding sites for the ubiquitous factor CTF/NF-I. Two of these sites are close to the TATA-box, whereas the others are located on both sides of the estrogen responsive unit formed by two imperfect estrogen response elements. Moreover two liver-enriched factors, C/EBP and HNF3, were found to interact with multiple closely spaced proximal promoter elements in the first 100 base pairs upstream of the TATA-box. To confirm the physiological significance of this in vitro analysis, in vivo DNase I footprinting experiments were carried out using the ligation-mediated polymerase chain reaction technique. The various cis-elements characterized in vitro as binding sites for known transcription factors and more particularly for liver-enriched transcription factors are efficiently recognized in vivo as well, suggesting that they play an important role in the control of the liver-specific vitellogenin gene B1 expression.

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.

Unique genome replication mechanism of the archaeal virus AFV1.

The exceptional genomic content and genome organization of the Acidianus filamentous virus 1 (AFV1) that infects the hyperthermophilic archaeon Acidianus hospitalis suggest that this virus might exploit an unusual mechanism of genome replication. An analysis of replicative intermediates of the viral genome by two-dimensional (2D) agarose gel electrophoresis revealed that viral genome replication starts by the formation of a D-loop and proceeds via strand displacement replication. Characterization of replicative intermediates using dark-field electron microscopy, in combination with the 2D agarose gel electrophoresis data, suggests that recombination plays a key role in the termination of AFV1 genome replication through the formation of terminal loops. A terminal protein was found to be attached to the ends of the viral genome. The results allow us to postulate a model of genome replication that relies on recombination events for initiation and termination.

Novel markers of the human follicle-associated epithelium identified by genomic profiling and microdissection.

BACKGROUND & AIMS: Regulation of gene expression in the follicle-associated epithelium (FAE) over Peyer's patches is largely unknown. CCL20, a chemokine that recruits immature dendritic cells, is one of the few FAE-specific markers described so far. Lymphotoxin beta (LTalpha1beta2) expressed on the membrane of immune cells triggers CCL20 expression in enterocytes. In this study, we measured expression profiles of LTalpha1beta2-treated intestinal epithelial cells and selected CCL20 -coregulated genes to identify new FAE markers. METHODS: Genomic profiles of T84 and Caco-2 cell lines treated with either LTalpha1beta2, flagellin, or tumor necrosis factor alpha were measured using the Affymetrix GeneChip U133A. Clustering analysis was used to select CCL20 -coregulated genes, and laser dissection microscopy and real-time polymerase chain reaction on human biopsy specimens was used to assess the expression of the selected markers. RESULTS: Applying a 2-way analysis of variance, we identified regulated genes upon the different treatments. A subset of genes involved in inflammation and related to the nuclear factor kappaB pathway was coregulated with CCL20 . Among these genes, the antiapoptotic factor TNFAIP3 was highly expressed in the FAE. CCL23 , which was not coregulated in vitro with CCL20 , was also specifically expressed in the FAE. CONCLUSIONS: We have identified 2 novel human FAE specifically expressed genes. Most of the CCL20 -coregulated genes did not show FAE-specific expression, suggesting that other signaling pathways are critical to modulate FAE-specific gene expression.

Defining responses to therapy and study outcomes in clinical trials of invasive fungal diseases: Mycoses Study Group and European Organization for Research and Treatment of Cancer consensus criteria.

Invasive fungal diseases (IFDs) have become major causes of morbidity and mortality among highly immunocompromised patients. Authoritative consensus criteria to diagnose IFD have been useful in establishing eligibility criteria for antifungal trials. There is an important need for generation of consensus definitions of outcomes of IFD that will form a standard for evaluating treatment success and failure in clinical trials. Therefore, an expert international panel consisting of the Mycoses Study Group and the European Organization for Research and Treatment of Cancer was convened to propose guidelines for assessing treatment responses in clinical trials of IFDs and for defining study outcomes. Major fungal diseases that are discussed include invasive disease due to Candida species, Aspergillus species and other molds, Cryptococcus neoformans, Histoplasma capsulatum, and Coccidioides immitis. We also discuss potential pitfalls in assessing outcome, such as conflicting clinical, radiological, and/or mycological data and gaps in knowledge.

Genomic approaches to the study of HIV‐1 acquisition.

Host genome studies are increasingly available for the study of infectious disease susceptibility. Current technologies include large-scale genotyping, genome-wide screens such as transcriptome and silencing (silencing RNA) studies, and increasingly, the possibility to sequence complete genomes. These approaches are of interest for the study of individuals who remain uninfected despite documented exposure to human immunodeficiency virus type 1. The main limitation remains the ascertainment of exposure and establishing large cohorts of informative individuals. The pattern of enrichment for CCR5 Δ32 homozygosis should serve as the standard for assessing the extent to which a given cohort (of white subjects) includes a large proportion of exposed uninfected individuals.

Genometric analyses of the organization of circular chromosomes: a universal pressure determines the direction of ribosomal RNA genes transcription relative to chromosome replication.

Selective pressures related to gene function and chromosomal architecture are acting on genome sequences and can be revealed, for instance, by appropriate genometric methods. Cumulative nucleotide skew analyses, i.e., GC, TA, and ORF orientation skews, predict the location of the origin of DNA replication for 88 out of 100 completely sequenced bacterial chromosomes. These methods appear fully reliable for proteobacteria, Gram-positives, and spirochetes as well as for euryarchaeotes. Based on this genome architecture information, coorientation analyses reveal that in prokaryotes, ribosomal RNA (rRNA) genes encoding the small and large ribosomal subunits are all transcribed in the same direction as DNA replication; that is, they are located along the leading strand. This result offers a simple and reliable method for circumscribing the region containing the origin of the DNA replication and reveals a strong selective pressure acting on the orientation of rRNA genes similar to the weaker one acting on the orientation of ORFs. Rate of coorientation of transfer RNA (tRNA) genes with DNA replication appears to be taxon-specific. Analyzing nucleotide biases such as GC and TA skews of genes and plotting one against the other reveals a taxonomic clusterization of species. All ribosomal RNA genes are enriched in Gs and depleted in Cs, the only so far known exception being the rRNA genes of deuterostomian mitochondria. However, this exception can be explained by the fact that in the chromosome of the human mitochondrion, the model of the deuterostomian organelle genome, DNA replication, and rRNA transcription proceed in opposite directions. A general rule is deduced from prokaryotic and mitochondrial genomes: ribosomal RNA genes that are transcribed in the same direction as the DNA replication are enriched in Gs, and those transcribed in the opposite direction are depleted in Gs.