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

A bacterial artificial chromosome library for Biomphalaria glabrata, intermediate snail host of Schistosoma mansoni

To provide a novel resource for analysis of the genome of Biomphalaria glabrata, members of the international Biomphalaria glabrata Genome Initiative (biology.unm.edu/biomphalaria-genome.html), working with the Arizona Genomics Institute (AGI) and supported by the National Human Genome Research Institute (NHGRI), produced a high quality bacterial artificial chromosome (BAC) library. The BB02 strain B. glabrata, a field isolate (Belo Horizonte, Minas Gerais, Brasil) that is susceptible to several strains of Schistosoma mansoni, was selfed for two generations to reduce haplotype diversity in the offspring. High molecular weight DNA was isolated from ovotestes of 40 snails, partially digested with HindIII, and ligated into pAGIBAC1 vector. The resulting B. glabrata BAC library (BG_BBa) consists of 61824 clones (136.3 kb average insert size) and provides 9.05 × coverage of the 931 Mb genome. Probing with single/low copy number genes from B. glabrata and fingerprinting of selected BAC clones indicated that the BAC library sufficiently represents the gene complement. BAC end sequence data (514 reads, 299860 nt) indicated that the genome of B. glabrata contains ~ 63% AT, and disclosed several novel genes, transposable elements, and groups of high frequency sequence elements. This BG_BBa BAC library, available from AGI at cost to the research community, gains in relevance because BB02 strain B. glabrata is targeted whole genome sequencing by NHGRI.

Methylation alterations are not a major cause of PTTG1 misregulation.

BACKGROUND On its physiological cellular context, PTTG1 controls sister chromatid segregation during mitosis. Within its crosstalk to the cellular arrest machinery, relies a checkpoint of integrity for which gained the over name of securin. PTTG1 was found to promote malignant transformation in 3T3 fibroblasts, and further found to be overexpressed in different tumor types. More recently, PTTG1 has been also related to different processes such as DNA repair and found to trans-activate different cellular pathways involving c-myc, bax or p53, among others. PTTG1 over-expression has been correlated to a worse prognosis in thyroid, lung, colorectal cancer patients, and it can not be excluded that this effect may also occur in other tumor types. Despite the clinical relevance and the increasing molecular characterization of PTTG1, the reason for its up-regulation remains unclear. METHOD We analysed PTTG1 differential expression in PC-3, DU-145 and LNCaP tumor cell lines, cultured in the presence of the methyl-transferase inhibitor 5-Aza-2'-deoxycytidine. We also tested whether the CpG island mapping PTTG1 proximal promoter evidenced a differential methylation pattern in differentiated thyroid cancer biopsies concordant to their PTTG1 immunohistochemistry status. Finally, we performed whole-genome LOH studies using Affymetix 50 K microarray technology and FRET analysis to search for allelic imbalances comprising the PTTG1 locus. CONCLUSION Our data suggest that neither methylation alterations nor LOH are involved in PTTG1 over-expression. These data, together with those previously reported, point towards a post-transcriptional level of misregulation associated to PTTG1 over-expression.

Regulatory elements involved in the post-transcriptional control of stage-specific gene expression in Trypanosoma cruzi: a review

Trypanosoma cruzi, a protozoan parasite that causes Chagas disease, exhibits unique mechanisms for gene expression such as constitutive polycistronic transcription of protein-coding genes, RNA editing and trans-splicing. In the absence of mechanism controlling transcription initiation, organized subsets of T. cruzi genes must be post-transcriptionally co-regulated in response to extracellular signals. The mechanisms that regulate stage-specific gene expression in this parasite have become much clearer through sequencing its whole genome as well as performing various proteomic and microarray analyses, which have demonstrated that at least half of the T. cruzi genes are differentially regulated during its life cycle. In this review, we attempt to highlight the recent advances in characterising cis and trans-acting elements in the T. cruzi genome that are involved in its post-transcriptional regulatory machinery.

ALDH1A3 loss of function causes bilateral anophthalmia/microphthalmia and hypoplasia of the optic nerve and optic chiasm.

The major active retinoid, all-trans retinoic acid, has long been recognized as critical for the development of several organs, including the eye. Mutations in STRA6, the gene encoding the cellular receptor for vitamin A, in patients with Matthew-Wood syndrome and anophthalmia/microphthalmia (A/M), have previously demonstrated the importance of retinol metabolism in human eye disease. We used homozygosity mapping combined with next-generation sequencing to interrogate patients with anophthalmia and microphthalmia for new causative genes. We used whole-exome and whole-genome sequencing to study a family with two affected brothers with bilateral A/M and a simplex case with bilateral anophthalmia and hypoplasia of the optic nerve and optic chiasm. Analysis of novel sequence variants revealed homozygosity for two nonsense mutations in ALDH1A3, c.568A>G, predicting p.Lys190*, in the familial cases, and c.1165A>T, predicting p.Lys389*, in the simplex case. Both mutations predict nonsense-mediated decay and complete loss of function. We performed antisense morpholino (MO) studies in Danio rerio to characterize the developmental effects of loss of Aldh1a3 function. MO-injected larvae showed a significant reduction in eye size, and aberrant axonal projections to the tectum were noted. We conclude that ALDH1A3 loss of function causes anophthalmia and aberrant eye development in humans and in animal model systems.

Next-generation sequencing of experimental mouse strains.

Since the turn of the century the complete genome sequence of just one mouse strain, C57BL/6J, has been available. Knowing the sequence of this strain has enabled large-scale forward genetic screens to be performed, the creation of an almost complete set of embryonic stem (ES) cell lines with targeted alleles for protein-coding genes, and the generation of a rich catalog of mouse genomic variation. However, many experiments that use other common laboratory mouse strains have been hindered by a lack of whole-genome sequence data for these strains. The last 5 years has witnessed a revolution in DNA sequencing technologies. Recently, these technologies have been used to expand the repertoire of fully sequenced mouse genomes. In this article we review the main findings of these studies and discuss how the sequence of mouse genomes is helping pave the way from sequence to phenotype. Finally, we discuss the prospects for using de novo assembly techniques to obtain high-quality assembled genome sequences of these laboratory mouse strains, and what advances in sequencing technologies may be required to achieve this goal.

Large-scale analysis of orthologs and paralogs under covarion-like and constant-but-different models of amino acid evolution.

Functional divergence between homologous proteins is expected to affect amino acid sequences in two main ways, which can be considered as proxies of biochemical divergence: a "covarion-like" pattern of correlated changes in evolutionary rates, and switches in conserved residues ("conserved but different"). Although these patterns have been used in case studies, a large-scale analysis is needed to estimate their frequency and distribution. We use a phylogenomic framework of animal genes to answer three questions: 1) What is the prevalence of such patterns? 2) Can we link such patterns at the amino acid level with selection inferred at the codon level? 3) Are patterns different between paralogs and orthologs? We find that covarion-like patterns are more frequently detected than "constant but different," but that only the latter are correlated with signal for positive selection. Finally, there is no obvious difference in patterns between orthologs and paralogs.

Functional diversification of duplicate genes through subcellular adaptation of encoded proteins.

BACKGROUND: Gene duplication is the primary source of new genes with novel or altered functions. It is known that duplicates may obtain these new functional roles by evolving divergent expression patterns and/or protein functions after the duplication event. Here, using yeast (Saccharomyces cerevisiae) as a model organism, we investigate a previously little considered mode for the functional diversification of duplicate genes: subcellular adaptation of encoded proteins. RESULTS: We show that for 24-37% of duplicate gene pairs derived from the S. cerevisiae whole-genome duplication event, the two members of the pair encode proteins that localize to distinct subcellular compartments. The propensity of yeast duplicate genes to evolve new localization patterns depends to a large extent on the biological function of their progenitor genes. Proteins involved in processes with a wider subcellular distribution (for example, catabolism) frequently evolved new protein localization patterns after duplication, whereas duplicate proteins limited to a smaller number of organelles (for example, highly expressed biosynthesis/housekeeping proteins with a slow rate of evolution) rarely relocate within the cell. Paralogous proteins evolved divergent localization patterns by partitioning of ancestral localizations ('sublocalization'), but probably more frequently by relocalization to new compartments ('neolocalization'). We show that such subcellular reprogramming may occur through selectively driven substitutions in protein targeting sequences. Notably, our data also reveal that relocated proteins functionally adapted to their new subcellular environments and evolved new functional roles through changes of their physico-chemical properties, expression levels, and interaction partners. CONCLUSION: We conclude that protein subcellular adaptation represents a common mechanism for the functional diversification of duplicate genes.

Genetics of calcium homeostasis in humans: continuum between monogenic diseases and continuous phenotypes

Extracellular calcium participates in several key physiological functions, such as control of blood coagulation, bone calcification or muscle contraction. Calcium homeostasis in humans is regulated in part by genetic factors, as illustrated by rare monogenic diseases characterized by hypo or hypercalcaemia. Both serum calcium and urinary calcium excretion are heritable continuous traits in humans. Serum calcium levels are tightly regulated by two main hormonal systems, i.e. parathyroid hormone and vitamin D, which are themselves also influenced by genetic factors. Recent technological advances in molecular biology allow for the screening of the human genome at an unprecedented level of detail and using hypothesis-free approaches, such as genome-wide association studies (GWAS). GWAS identified novel loci for calcium-related phenotypes (i.e. serum calcium and 25-OH vitamin D) that shed new light on the biology of calcium in humans. The substantial overlap (i.e. CYP24A1, CASR, GATA3; CYP2R1) between genes involved in rare monogenic diseases and genes located within loci identified in GWAS suggests a genetic and phenotypic continuum between monogenic diseases of calcium homeostasis and slight disturbances of calcium homeostasis in the general population. Future studies using whole-exome and whole-genome sequencing will further advance our understanding of the genetic architecture of calcium homeostasis in humans. These findings will likely provide new insight into the complex mechanisms involved in calcium homeostasis and hopefully lead to novel preventive and therapeutic approaches. Keyword: calcium, monogenic, genome-wide association studies, genetics.

Functional analysis: evaluation of response intensities--tailoring ANOVA for lists of expression subsets.

Background: Microarray data is frequently used to characterize the expression profile of a whole genome and to compare the characteristics of that genome under several conditions. Geneset analysis methods have been described previously to analyze the expression values of several genes related by known biological criteria (metabolic pathway, pathology signature, co-regulation by a common factor, etc.) at the same time and the cost of these methods allows for the use of more values to help discover the underlying biological mechanisms. Results: As several methods assume different null hypotheses, we propose to reformulate the main question that biologists seek to answer. To determine which genesets are associated with expression values that differ between two experiments, we focused on three ad hoc criteria: expression levels, the direction of individual gene expression changes (up or down regulation), and correlations between genes. We introduce the FAERI methodology, tailored from a two-way ANOVA to examine these criteria. The significance of the results was evaluated according to the self-contained null hypothesis, using label sampling or by inferring the null distribution from normally distributed random data. Evaluations performed on simulated data revealed that FAERI outperforms currently available methods for each type of set tested. We then applied the FAERI method to analyze three real-world datasets on hypoxia response. FAERI was able to detect more genesets than other methodologies, and the genesets selected were coherent with current knowledge of cellular response to hypoxia. Moreover, the genesets selected by FAERI were confirmed when the analysis was repeated on two additional related datasets. Conclusions: The expression values of genesets are associated with several biological effects. The underlying mathematical structure of the genesets allows for analysis of data from several genes at the same time. Focusing on expression levels, the direction of the expression changes, and correlations, we showed that two-step data reduction allowed us to significantly improve the performance of geneset analysis using a modified two-way ANOVA procedure, and to detect genesets that current methods fail to detect.

Oviposition by pierid butterflies triggers defense responses in Arabidopsis.

Insect eggs represent a threat for the plant as hatching larvae rapidly start with their feeding activity. Using a whole-genome microarray, we studied the expression profile of Arabidopsis (Arabidopsis thaliana) leaves after oviposition by two pierid butterflies. For Pieris brassicae, the deposition of egg batches changed the expression of hundreds of genes over a period of 3 d after oviposition. The transcript signature was similar to that observed during a hypersensitive response or in lesion-mimic mutants, including the induction of defense and stress-related genes and the repression of genes involved in growth and photosynthesis. Deposition of single eggs by Pieris rapae caused a similar although much weaker transcriptional response. Analysis of the jasmonic acid and salicylic acid mutants coi1-1 and sid2-1 indicated that the response to egg deposition is mostly independent of these signaling pathways. Histochemical analyses showed that egg deposition is causing a localized cell death, accompanied by the accumulation of callose, and the production of reactive oxygen species. In addition, activation of the pathogenesis-related1::beta-glucuronidase reporter gene correlated precisely with the site of egg deposition and was also triggered by crude egg extract. This study provides molecular evidence for the detection of egg deposition by Arabidopsis plants and suggests that oviposition causes a localized response with strong similarity to a hypersensitive response.

Stick insect genomes reveal natural selection's role in parallel speciation.

Natural selection can drive the repeated evolution of reproductive isolation, but the genomic basis of parallel speciation remains poorly understood. We analyzed whole-genome divergence between replicate pairs of stick insect populations that are adapted to different host plants and undergoing parallel speciation. We found thousands of modest-sized genomic regions of accentuated divergence between populations, most of which are unique to individual population pairs. We also detected parallel genomic divergence across population pairs involving an excess of coding genes with specific molecular functions. Regions of parallel genomic divergence in nature exhibited exceptional allele frequency changes between hosts in a field transplant experiment. The results advance understanding of biological diversification by providing convergent observational and experimental evidence for selection's role in driving repeatable genomic divergence.

Evaluation of the chicken transcriptome by SAGE of B cells and the DT40 cell line

Background: The understanding of whole genome sequences in higher eukaryotes depends to a large degree on the reliable definition of transcription units including exon/intron structures, translated open reading frames (ORFs) and flanking untranslated regions. The best currently available chicken transcript catalog is the Ensembl build based on the mappings of a relatively small number of full length cDNAs and ESTs to the genome as well as genome sequence derived in silico gene predictions.Results: We use Long Serial Analysis of Gene Expression (LongSAGE) in bursal lymphocytes and the DT40 cell line to verify the quality and completeness of the annotated transcripts. 53.6% of the more than 38,000 unique SAGE tags (unitags) match to full length bursal cDNAs, the Ensembl transcript build or the genome sequence. The majority of all matching unitags show single matches to the genome, but no matches to the genome derived Ensembl transcript build. Nevertheless, most of these tags map close to the 3' boundaries of annotated Ensembl transcripts.Conclusions: These results suggests that rather few genes are missing in the current Ensembl chicken transcript build, but that the 3' ends of many transcripts may not have been accurately predicted. The tags with no match in the transcript sequences can now be used to improve gene predictions, pinpoint the genomic location of entirely missed transcripts and optimize the accuracy of gene finder software.

Differential contribution of transcription factors to Arabidopsis thaliana defense against Spodoptera littoralis

In response to insect herbivory, Arabidopsis plants activate the synthesis of the phytohor- mone jasmonate-isoleucine, which binds to a complex consisting of the receptor COI1 and JAZ repressors. Upon proteasome-mediated JAZ degradation, basic helix-loop-helix tran- scription factors (TFs) MYC2, MYC3, and MYC4 become activated and this results in the expression of defense genes. Although the jasmonate (JA) pathway is known to be essen- tial for the massive transcriptional reprogramming that follows herbivory, there is however little information on other TFs that are required for defense against herbivores and whether they contribute significantly to JA-dependent defense gene expression. By transcriptome profiling, we identified 41TFs that were induced in response to herbivory by the generalist Spodoptera littoralis. Among them, nine genes, including WRKY18, WRKY40, ANAC019, ANAC055, ZAT10, ZAT12, AZF2, ERF13, and RRTF1, were found to play a significant role in resistance to S. littoralis herbivory. Compared to the triple mutant myc234 that is as sensitive as coi1-1 to herbivory, knockout lines of these nine TFs were only partially more sensitive to S. littoralis but, however, some displayed distinct gene expression changes at the whole-genome level. Data thus reveal that MYC2, MYC3, and MYC4 are master regu- lators of Arabidopsis resistance to a generalist herbivore and identify new genes involved in insect defense.

Extreme heterochiasmy and nascent sex chromosomes in European tree frogs.

We investigated sex-specific recombination rates in Hyla arborea, a species with nascent sex chromosomes and male heterogamety. Twenty microsatellites were clustered into six linkage groups, all showing suppressed or very low recombination in males. Seven markers were sex linked, none of them showing any sign of recombination in males (r=0.00 versus 0.43 on average in females). This opposes classical models of sex chromosome evolution, which envision an initially small differential segment that progressively expands as structural changes accumulate on the Y chromosome. For autosomes, maps were more than 14 times longer in females than in males, which seems the highest ratio documented so far in vertebrates. These results support the pleiotropic model of Haldane and Huxley, according to which recombination is reduced in the heterogametic sex by general modifiers that affect recombination on the whole genome.