The fine-scale architecture of structural variants in 17 mouse genomes.

BACKGROUND: Accurate catalogs of structural variants (SVs) in mammalian genomes are necessary to elucidate the potential mechanisms that drive SV formation and to assess their functional impact. Next generation sequencing methods for SV detection are an advance on array-based methods, but are almost exclusively limited to four basic types: deletions, insertions, inversions and copy number gains. RESULTS: By visual inspection of 100 Mbp of genome to which next generation sequence data from 17 inbred mouse strains had been aligned, we identify and interpret 21 paired-end mapping patterns, which we validate by PCR. These paired-end mapping patterns reveal a greater diversity and complexity in SVs than previously recognized. In addition, Sanger-based sequence analysis of 4,176 breakpoints at 261 SV sites reveal additional complexity at approximately a quarter of structural variants analyzed. We find micro-deletions and micro-insertions at SV breakpoints, ranging from 1 to 107 bp, and SNPs that extend breakpoint micro-homology and may catalyze SV formation. CONCLUSIONS: An integrative approach using experimental analyses to train computational SV calling is essential for the accurate resolution of the architecture of SVs. We find considerable complexity in SV formation; about a quarter of SVs in the mouse are composed of a complex mixture of deletion, insertion, inversion and copy number gain. Computational methods can be adapted to identify most paired-end mapping patterns.

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.

The Eukaryotic Promoter Database EPD: the impact of in silico primer extension.

The Eukaryotic Promoter Database (EPD) is an annotated non-redundant collection of eukaryotic POL II promoters, experimentally defined by a transcription start site (TSS). There may be multiple promoter entries for a single gene. The underlying experimental evidence comes from journal articles and, starting from release 73, from 5' ESTs of full-length cDNA clones used for so-called in silico primer extension. Access to promoter sequences is provided by pointers to TSS positions in nucleotide sequence entries. The annotation part of an EPD entry includes a description of the type and source of the initiation site mapping data, links to other biological databases and bibliographic references. EPD is structured in a way that facilitates dynamic extraction of biologically meaningful promoter subsets for comparative sequence analysis. Web-based interfaces have been developed that enable the user to view EPD entries in different formats, to select and extract promoter sequences according to a variety of criteria and to navigate to related databases exploiting different cross-references. Tools for analysing sequence motifs around TSSs defined in EPD are provided by the signal search analysis server. EPD can be accessed at http://www.epd. isb-sib.ch.

Vaccination with a Melan-A peptide selects an oligoclonal T cell population with increased functional avidity and tumor reactivity.

Both the underlying molecular mechanisms and the kinetics of TCR repertoire selection following vaccination against tumor Ags in humans have remained largely unexplored. To gain insight into these questions, we performed a functional and structural longitudinal analysis of the TCR of circulating CD8(+) T cells specific for the HLA-A2-restricted immunodominant epitope from the melanocyte differentiation Ag Melan-A in a melanoma patient who developed a vigorous and sustained Ag-specific T cell response following vaccination with the corresponding synthetic peptide. We observed an increase in functional avidity of Ag recognition and in tumor reactivity in the postimmune Melan-A-specific populations as compared with the preimmune blood sample. Improved Ag recognition correlated with an increase in the t(1/2) of peptide/MHC interaction with the TCR as assessed by kinetic analysis of A2/Melan-A peptide multimer staining decay. Ex vivo analysis of the clonal composition of Melan-A-specific CD8(+) T cells at different time points during vaccination revealed that the response was the result of asynchronous expansion of several distinct T cell clones. Some of these T cell clones were also identified at a metastatic tumor site. Collectively, these data show that tumor peptide-driven immune stimulation leads to the selection of high-avidity T cell clones of increased tumor reactivity that independently evolve within oligoclonal populations.

Soil fungal communities of grasslands are environmentally structured at a regional scale in the Alps.

Studying patterns of species distributions along elevation gradients is frequently used to identify the primary factors that determine the distribution, diversity and assembly of species. However, despite their crucial role in ecosystem functioning, our understanding of the distribution of below-ground fungi is still limited, calling for more comprehensive studies of fungal biogeography along environmental gradients at various scales (from regional to global). Here, we investigated the richness of taxa of soil fungi and their phylogenetic diversity across a wide range of grassland types along a 2800 m elevation gradient at a large number of sites (213), stratified across a region of the Western Swiss Alps (700 km(2)). We used 454 pyrosequencing to obtain fungal sequences that were clustered into operational taxonomic units (OTUs). The OTU diversity-area relationship revealed uneven distribution of fungal taxa across the study area (i.e. not all taxa are everywhere) and fine-scale spatial clustering. Fungal richness and phylogenetic diversity were found to be higher in lower temperatures and higher moisture conditions. Climatic and soil characteristics as well as plant community composition were related to OTU alpha, beta and phylogenetic diversity, with distinct fungal lineages suggesting distinct ecological tolerances. Soil fungi, thus, show lineage-specific biogeographic patterns, even at a regional scale, and follow environmental determinism, mediated by interactions with plants.

Comparative transcriptomics of rice reveals an ancient pattern of response to microbial colonization.

Glomalean fungi induce and colonize symbiotic tissue called arbuscular mycorrhiza on the roots of most land plants. Other fungi also colonize plants but cause disease not symbiosis. Whole-transcriptome analysis using a custom-designed Affymetrix Gene-Chip and confirmation with real-time RT-PCR revealed 224 genes affected during arbuscular mycorrhizal symbiosis. We compared these transcription profiles with those from rice roots that were colonized by pathogens (Magnaporthe grisea and Fusarium moniliforme). Over 40% of genes showed differential regulation caused by both the symbiotic and at least one of the pathogenic interactions. A set of genes was similarly expressed in all three associations, revealing a conserved response to fungal colonization. The responses that were shared between pathogen and symbiont infection may play a role in compatibility. Likewise, the responses that are different may cause disease. Some of the genes that respond to mycorrhizal colonization may be involved in the uptake of phosphate. Indeed, phosphate addition mimicked the effect of mycorrhiza on 8% of the tested genes. We found that 34% of the mycorrhiza-associated rice genes were also associated with mycorrhiza in dicots, revealing a conserved pattern of response between the two angiosperm classes.

The ENCODE (ENCyclopedia Of DNA Elements) Project.

The ENCyclopedia Of DNA Elements (ENCODE) Project aims to identify all functional elements in the human genome sequence. The pilot phase of the Project is focused on a specified 30 megabases (approximately 1%) of the human genome sequence and is organized as an international consortium of computational and laboratory-based scientists working to develop and apply high-throughput approaches for detecting all sequence elements that confer biological function. The results of this pilot phase will guide future efforts to analyze the entire human genome.

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.

An exogenous mouse mammary tumor virus with properties of Mls-1a (Mtv-7).

The classical minor lymphocyte stimulating (Mls) antigens, which induce a strong primary T cell response in vitro, are closely linked to endogenous copies of mouse mammary tumor viruses (MMTV). Expression of Mls genes leads to clonal deletion of T cell subsets expressing specific T cell receptor (TCR) V beta chains. We describe the isolation and characterization of a new exogenous (infectious) MMTV with biological properties similar to the Mls antigen Mls-1a. In vivo administration of either Mls-1a-expressing B cells or the infectious MMTV (SW) led to an increase of T cells expressing V beta 6 followed by their deletion. Surprisingly, different kinetics of deletion were observed with the exogenous virus depending upon the route of infection. Infection through the mucosa led to a slow deletion of V beta 6+ T cells, whereas deletion was rapid after subcutaneous infection. Sequence analysis of the open reading frames in the 3' long terminal repeat of both this exogenous MMTV (SW) and of Mtv-7 (which is closely linked to Mls-1a) revealed striking similarities, particularly in the COOH terminus, which has been implicated in TCR V beta recognition. The identification of an infectious MMTV with the properties of a strong Mls antigen provides a new, powerful tool to study immunity and tolerance in vivo.

STAT3 links IL-22 signaling in intestinal epithelial cells to mucosal wound healing.

Signal transducer and activator of transcription (STAT) 3 is a pleiotropic transcription factor with important functions in cytokine signaling in a variety of tissues. However, the role of STAT3 in the intestinal epithelium is not well understood. We demonstrate that development of colonic inflammation is associated with the induction of STAT3 activity in intestinal epithelial cells (IECs). Studies in genetically engineered mice showed that epithelial STAT3 activation in dextran sodium sulfate colitis is dependent on interleukin (IL)-22 rather than IL-6. IL-22 was secreted by colonic CD11c(+) cells in response to Toll-like receptor stimulation. Conditional knockout mice with an IEC-specific deletion of STAT3 activity were highly susceptible to experimental colitis, indicating that epithelial STAT3 regulates gut homeostasis. STAT3(IEC-KO) mice, upon induction of colitis, showed a striking defect of epithelial restitution. Gene chip analysis indicated that STAT3 regulates the cellular stress response, apoptosis, and pathways associated with wound healing in IECs. Consistently, both IL-22 and epithelial STAT3 were found to be important in wound-healing experiments in vivo. In summary, our data suggest that intestinal epithelial STAT3 activation regulates immune homeostasis in the gut by promoting IL-22-dependent mucosal wound healing.

Gene finding in the chicken genome.

BACKGROUND: Despite the continuous production of genome sequence for a number of organisms, reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularly true for genomes for which there is not a large collection of known gene sequences, such as the recently published chicken genome. We used the chicken sequence to test comparative and homology-based gene-finding methods followed by experimental validation as an effective genome annotation method. RESULTS: We performed experimental evaluation by RT-PCR of three different computational gene finders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram was computed and each component of it was evaluated. The results showed that de novo comparative methods can identify up to about 700 chicken genes with no previous evidence of expression, and can correctly extend about 40% of homology-based predictions at the 5' end. CONCLUSIONS: De novo comparative gene prediction followed by experimental verification is effective at enhancing the annotation of the newly sequenced genomes provided by standard homology-based methods.

Characterization of the rice PHO1 gene family reveals a key role for OsPHO1;2 in phosphate homeostasis and the evolution of a distinct clade in dicotyledons.

Phosphate homeostasis was studied in a monocotyledonous model plant through the characterization of the PHO1 gene family in rice (Oryza sativa). Bioinformatics and phylogenetic analysis showed that the rice genome has three PHO1 homologs, which cluster with the Arabidopsis (Arabidopsis thaliana) AtPHO1 and AtPHO1;H1, the only two genes known to be involved in root-to-shoot transfer of phosphate. In contrast to the Arabidopsis PHO1 gene family, all three rice PHO1 genes have a cis-natural antisense transcript located at the 5 ' end of the genes. Strand-specific quantitative reverse transcription-PCR analyses revealed distinct patterns of expression for sense and antisense transcripts for all three genes, both at the level of tissue expression and in response to nutrient stress. The most abundantly expressed gene was OsPHO1;2 in the roots, for both sense and antisense transcripts. However, while the OsPHO1;2 sense transcript was relatively stable under various nutrient deficiencies, the antisense transcript was highly induced by inorganic phosphate (Pi) deficiency. Characterization of Ospho1;1 and Ospho1;2 insertion mutants revealed that only Ospho1;2 mutants had defects in Pi homeostasis, namely strong reduction in Pi transfer from root to shoot, which was accompanied by low-shoot and high-root Pi. Our data identify OsPHO1;2 as playing a key role in the transfer of Pi from roots to shoots in rice, and indicate that this gene could be regulated by its cis-natural antisense transcripts. Furthermore, phylogenetic analysis of PHO1 homologs in monocotyledons and dicotyledons revealed the emergence of a distinct clade of PHO1 genes in dicotyledons, which include members having roles other than long-distance Pi transport.

Changes in the transcriptional profile of transporters in the intestine along the anterior-posterior and crypt-villus axes.

BACKGROUND: The purpose of this work was to characterize the expression of drug and nutrient carriers along the anterior-posterior and crypt-villus axes of the intestinal epithelium and to study the validity of utilizing whole gut tissue rather than purified epithelial cells to examine regional variations in gene expression. RESULTS: We have characterized the mRNA expression profiles of 76 % of all currently known transporters along the anterior-posterior axis of the gut. This is the first study to describe the expression profiles of the majority of all known transporters in the intestine. The expression profiles of transporters, as defined according to the Gene Ontology consortium, were measured in whole tissue of the murine duodenum, jejunum, ileum and colon using high-density microarrays. For nine transporters (Abca1, Abcc1, Abcc3, Abcg8, Slc10a2, Slc28a2, Slc2a1, Slc34a2 and Slc5a8), the mRNA profiles were further measured by RT-PCR in laser micro-dissected crypt and villus epithelial cells corresponding to the aforementioned intestinal regions. With respect to differentially regulated transporters, the colon had a distinct expression profile from small intestinal segments. The majority (59 % for p cutoff < or = 0.05) of transporter mRNA levels were constant across the intestinal sections studied. For the transporter subclass "carrier activity", which contains the majority of known carriers for biologically active compounds, a significant change (p < or = 0.05) along the anterior-posterior axis was observed. CONCLUSION: All nine transporters examined in laser-dissected material demonstrated good replication of the region-specific profiles revealed by microarray. Furthermore, we suggest that the distribution characteristics of Slc5a8 along the intestinal tract render it a suitable candidate carrier for monocarboxylate drugs in the posterior portion of the intestine. Our findings also predict that there is a significant difference in the absorption of carrier-mediated compounds in the different intestinal segments. The most pronounced differences can be expected between the adjoining segments ileum and colon, but the differences between the other adjoining segments are not negligible. Finally, for the examined genes, profiles measured in whole intestinal tissue extracts are representative of epithelial cell-only gene expression.