ALBERT

Description: Background: Despite their widespread use, the biological mechanisms underlying the efficacy of psychotropic drugs are still incompletely known; improved understanding of these is essential for development of novel more effective drugs and rational design of therapy. Given the large number of psychotropic drugs available and their differential pharmacological effects, it would be important to establish specific predictors of response to various classes of drugs. Results: To identify the molecular mechanisms that may initiate therapeutic effects, whole-genome expression profiling (using 324 Illumina Mouse WG-6 microarrays) of drug-induced alterations in the mouse brain was undertaken, with a focus on the time-course (1, 2, 4 and 8 h) of gene expression changes produced by eighteen major psychotropic drugs: antidepressants, antipsychotics, anxiolytics, psychostimulants and opioids. The resulting database is freely accessible at www.genes2mind.org. Bioinformatics approaches led to the identification of three main drug-responsive genomic networks and indicated neurobiological pathways that mediate the alterations in transcription. Each tested psychotropic drug was characterized by a unique gene network expression profile related to its neuropharmacological properties. Functional links that connect expression of the networks to the development of neuronal adaptations (MAPK signaling pathway), control of brain metabolism (adipocytokine pathway), and organization of cell projections (mTOR pathway) were found. Conclusions: The comparison of gene expression alterations between various drugs opened a new means to classify the different psychoactive compounds and to predict their cellular targets; this is well exemplified in the case of tianeptine, an antidepressant with unknown mechanisms of action. This work represents the first proof-of-concept study of a molecular classification of psychoactive drugs.

Description: Background: Loss of CpG dinucleotides in genomic DNA through methylation-induced mutation is characteristic of vertebrates and plants. However, these and other eukaryotic phyla show a range of other dinucleotide frequency biases with currently uncharacterized underlying mutational or selection mechanisms. We developed a parameterized Markov process to identify what neighbour context-dependent mutations best accounted for patterns of dinucleotide frequency biases in genomic and cytoplasmically expressed mRNA sequences of different vertebrates, other eukaryotic groups and RNA viruses that infect them. Results: Consistently, 11- to 14-fold greater frequencies of the methylation-association mutation of C to T upstream of G (depicted C[rightwards arrow]T,G) than other transitions best modelled dinucleotide frequencies in mammalian genomic DNA. However, further mutations such as G[rightwards arrow]T,T (5-fold greater than the default transversion rate) were required to account for the full spectrum of dinucleotide frequencies in mammalian sequence datasets. Consistent with modeling predictions for these two mutations, instability of both CpG and CpT dinucleotides was identified through SNP frequency analysis of human DNA sequences. Different sets of context-dependent mutations were modelled in other eukaryotes with non-methylated genomic DNA. In contrast to genomic DNA, best-fit models of dinucleotide frequencies in transcribed RNA sequences expressed in the cytoplasm from all organisms were dominated by mutations that eliminated UpA dinucleotides, observations consistent with cytoplasmically driven selection for mRNA stability. Surprisingly, mRNA sequences from organisms with methylated genomes showed evidence for additional selection against CpG through further context-dependent mutations (eg. C[rightwards arrow]A,G). Similar mutation or selection processes were identified among single-stranded mammalian RNA viruses; these potentially account for their previously described but unexplained under-representations of CpG and UpA dinucleotides. Conclusions: Methods we have developed identify mutational processes and selection pressures in organisms provide new insights into nucleotide compositional constraints and a wealth of biochemical and evolutionarily testable predictions for the future.

Description: Background: Temperature is one of key environmental parameters that affect the whole life of fishes and an increasing number of studies have been directed towards understanding the mechanisms of cold acclimation in fish. However, the adaptation of larvae to cold stress and the cold-specific transcriptional alterations in fish larvae remain largely unknown. In this study, we characterized the development of cold-tolerance in zebrafish larvae and investigated the transcriptional profiles under cold stress using RNA-seq. Results: Pre-exposure of 96hpf zebrafish larvae to cold stress (16[degree sign]C) for 24 h significantly increased their survival rates under severe cold stress (12[degree sign]C). RNA-seq generated 272 million raw reads from six sequencing libraries and about 92% of the processed reads were mapped to the reference genome of zebrafish. Differential expression analysis identified 1,431 up- and 399 down-regulated genes. Gene ontology enrichment analysis of cold-induced genes revealed that RNA splicing, ribosome biogenesis and protein catabolic process were the most highly overrepresented biological processes. Spliceosome, proteasome, eukaryotic ribosome biogenesis and RNA transport were the most highly enriched pathways for genes up-regulated by cold stress. Moreover, alternative splicing of 197 genes and promoter switching of 64 genes were found to be regulated by cold stress. A shorter isoform of stk16 that lacks 67 amino acids at the N-terminus was specifically generated by skipping the second exon in cold-treated larvae. Alternative promoter usage was detected for per3 gene under cold stress, which leading to a highly up-regulated transcript encoding a truncated protein lacking the C-terminal domains. Conclusions: These findings indicate that zebrafish larvae possess the ability to build cold-tolerance under mild low temperature and transcriptional and post-transcriptional regulations are extensively involved in this acclimation process.

Description: Background: Gene reshuffling, point mutations and horizontal gene transfer contribute to bacterial genome variation, but require the genome to rewire its transcriptional circuitry to ensure that inserted, mutated or reshuffled genes are transcribed at appropriate levels. The genomes of Epsilonproteobacteria display very low synteny, due to high levels of reshuffling and reorganisation of gene order, but still share a significant number of gene orthologs allowing comparison. Here we present the primary transcriptome of the pathogenic Epsilonproteobacterium Campylobacter jejuni, and have used this for comparative and predictive transcriptomics in the Epsilonproteobacteria. Results: Differential RNA-sequencing using 454 sequencing technology was used to determine the primary transcriptome of C. jejuni NCTC 11168, which consists of 992 transcription start sites (TSS), which included 29 putative non-coding and stable RNAs, 266 intragenic (internal) TSS, and 206 antisense TSS. Several previously unknown features were identified in the C. jejuni transcriptional landscape, like leaderless mRNAs and potential leader peptides upstream of amino acid biosynthesis genes. A cross-species comparison of the primary transcriptomes of C. jejuni and the related Epsilonproteobacterium Helicobacter pylori highlighted a lack of conservation of operon organisation, position of intragenic and antisense promoters or leaderless mRNAs. Predictive comparisons using 40 other Epsilonproteobacterial genomes suggests that this lack of conservation of transcriptional features is common to all Epsilonproteobacterial genomes, and is associated with the absence of genome synteny in this subdivision of the Proteobacteria. Conclusions: Both the genomes and transcriptomes of Epsilonproteobacteria are highly variable, both at the genome level by combining and division of multicistronic operons, but also on the gene level by generation or deletion of promoter sequences and 5[prime] untranslated regions. Regulatory features may have evolved after these species split from a common ancestor, with transcriptome rewiring compensating for changes introduced by genomic reshuffling and horizontal gene transfer.

Description: Background: External development and optical transparency of embryos make zebrafish exceptionally suitable for in vivo insertional mutagenesis using fluorescent proteins to visualize expression patterns of mutated genes. Recently developed Gene Breaking Transposon (GBT) vectors greatly improve the fidelity and mutagenicity of transposon-based gene trap vectors. Results: We constructed and tested a bipartite GBT vector with Gal4-VP16 as the primary gene trap reporter. Our vector also contains a UAS:eGFP cassette for direct detection of gene trap events by fluorescence. To confirm gene trap events, we generated a UAS:mRFP tester line. We screened 270 potential founders and established 41 gene trap lines. Three of our gene trap alleles display homozygous lethal phenotypes ranging from embryonic to late larval: nsftpl6, atp1a3atpl10 and flrtpl19. Our gene trap cassette is flanked by direct loxP sites, which enabled us to successfully revert nsftpl6, atp1a3atpl10 and flrtpl19 gene trap alleles by injection of Cre mRNA. The UAS:eGFP cassette is flanked by direct FRT sites. It can be readily removed by injection of Flp mRNA for use of our gene trap alleles with other tissue-specific GFP-marked lines. The Gal4-VP16 component of our vector provides two important advantages over other GBT vectors. The first is increased sensitivity, which enabled us to detect previously unnoticed expression of nsf in the pancreas. The second advantage is that all our gene trap lines, including integrations into non-essential genes, can be used as highly specific Gal4 drivers for expression of other transgenes under the control of Gal4 UAS. Conclusions: The Gal4-containing bipartite Gene Breaking Transposon vector presented here retains high specificity for integrations into genes, high mutagenicity and revertibility by Cre. These features, together with utility as highly specific Gal4 drivers, make gene trap mutants presented here especially useful to the research community.

Description: Background: Very little is known about manganese (Mn)-toxicity-responsive genes in citrus plants. Seedlings of 'Xuegan' (Citrus sinensis) and 'Sour pummelo' (Citrus grandis) were irrigated for 17 weeks with nutrient solution containing 2 muM (control) or 600 muM (Mn-toxicity) MnSO4. The objectives of this study were to understand the mechanisms of citrus Mn-tolerance and to identify differentially expressed genes, which might be involved in Mn-tolerance. Results: Under Mn-toxicity, the majority of Mn in seedlings was retained in the roots; C. sinensis seedlings accumulated more Mn in roots and less Mn in shoots (leaves) than C. grandis ones and Mn concentration was lower in Mn-toxicity C. sinensis leaves compared to Mn-toxicity C. grandis ones. Mn-toxicity affected C. grandis seedling growth, leaf CO2 assimilation, total soluble concentration, phosphorus (P) and magenisum (Mg) more than C. sinensis. Using cDNA-AFLP, we isolated 42 up-regulated and 80 down-regulated genes in Mn-toxicity C. grandis leaves. They were grouped into the following functional categories: biological regulation and signal transduction, carbohydrate and energy metabolism, nucleic acid metabolism, protein metabolism, lipid metabolism, cell wall metabolism, stress responses and cell transport. However, only 7 up-regulated and 8 down-regulated genes were identified in Mn-toxicity C. sinensis ones. The responses of C. grandis leaves to Mn-toxicity might include following several aspects: (1) accelerating leaf senescence; (2) activating the metabolic pathway related to ATPase synthesis and reducing power production; (3) decreasing cell transport; (4) inhibiting protein and nucleic acid metabolisms; (5) impairing the formation of cell wall; and (6) triggering multiple signal transduction pathways. We also identified many new Mn-toxicity-responsive genes involved in biological and signal transduction, carbohydrate and protein metabolisms, stress responses and cell transport. Conclusions: Our results demonstrated that C. sinensis was more tolerant to Mn-toxicity than C. grandis, and that Mn-toxicity affected gene expression far less in C. sinensis leaves. This might be associated with more Mn accumulation in roots and less Mn accumulation in leaves of Mn-toxicity C. sinensis seedlings than those of C. grandis seedlings. Our findings increase our understanding of the molecular mechanisms involved in the responses of plants to Mn-toxicity.

Description: Background: Two-partner secretion systems in Gram-negative bacteria consist of an outer membrane protein TpsB that mediates the secretion of a cognate TpsA protein into the extracellular milieu. TpsA proteins have diverse, often virulence-related functions, and some of them inhibit the growth of related bacteria. In Neisseria meningitidis, several functions have been attributed to the TpsA proteins. Downstream of the tpsB and tpsA genes, several shorter tpsA-related gene cassettes, called tpsC, are located interspersed with intervening open-reading frames (IORFs). It has been suggested that the tpsC cassettes may recombine with the tpsA gene as a mechanism of antigenic variation. Here, we investigated (i) whether TpsA of N. meningitidis also has growth-inhibitory properties, (ii) whether tpsC cassettes recombine with the tpsA gene, and (iii) what the consequences of such recombination events might be. Results: We demonstrate that meningococcal TpsA has growth-inhibitory properties and that the IORF located immediately downstream of tpsA confers immunity to the producing strain. Although bioinformatics analysis suggests that recombination between tpsC cassettes and tpsA occurs, detailed analysis of the tpsA gene in a large collection of disease isolates of three clonal complexes revealed that the frequency is very low and cannot be a mechanism of antigenic variation. However, recombination affected growth inhibition. In vitro experiments revealed that recombination can be mediated through acquirement of tpsC cassettes from the environment and it identified the regions involved in the recombination. Conclusions: Meningococcal TpsA has growth-inhibitory properties. Recombination between tpsA and tpsC cassettes occurs in vivo but is rare and has consequences for growth inhibition. A recombination model is proposed and we propose that the main goal of recombination is the collection of new IORFs for protection against a variety of TpsA proteins.

Description: Background: Sequencing of the genome of Propionibacterium acnes produced a catalogue of genes many of which enable this organism to colonise skin and survive exposure to the elements. Despite this platform, there was little understanding of the gene regulation that gives rise to an organism that has a major impact on human health and wellbeing and causes infections beyond the skin. To address this situation, we have undertaken a genome--wide study of gene regulation using a combination of improved differential and global RNA-sequencing and an analytical approach that takes into account the inherent noise within the data. Results: We have produced nucleotide-resolution transcriptome maps that identify and differentiate sites of transcription initiation from sites of stable RNA processing and mRNA cleavage. Moreover, analysis of these maps provides strong evidence for 'pervasive' transcription and shows that contrary to initial indications it is not biased towards the production of antisense RNAs. In addition, the maps reveal an extensive array of riboswitches, leaderless mRNAs and small non-protein-coding RNAs alongside vegetative promoters and post-transcriptional events, which includes unusual tRNA processing. The identification of such features will inform models of complex gene regulation, as illustrated here for ribonucleotide reductases and a potential quorum-sensing, two-component system. Conclusions: The approach described here, which is transferable to any bacterial species, has produced a step increase in whole-cell knowledge of gene regulation in P. acnes. Continued expansion of our maps to include transcription associated with different growth conditions and genetic backgrounds will provide a new platform from which to computationally model the gene expression that determines the physiology of P. acnes and its role in human disease.

Description: Background: Pancreatic cancer is a deadly disease with a five-year survival of less than 5%. A better understanding of the underlying biology may suggest novel therapeutic targets. Recent surveys of the pancreatic cancer genome have uncovered numerous new alterations; yet systematic functional characterization of candidate cancer genes has lagged behind. To address this challenge, here we have devised a highly-parallel RNA interference-based functional screen to evaluate many genomically-nominated candidate pancreatic cancer genes simultaneously. Results: For 185 candidate pancreatic cancer genes, selected from recurrently altered genomic loci, we performed a pooled shRNA library screen of cell growth/viability across 10 different cell lines. Knockdown-associated effects on cell growth were assessed by enrichment or depletion of shRNA hairpins, by hybridization to barcode microarrays. A novel analytical approach (COrrelated Phenotypes for On-Target Effects; COPOTE) was used to discern probable on-target knockdown, based on identifying different shRNAs targeting the same gene and displaying concordant phenotypes across cell lines. Knockdown data were integrated with genomic architecture and gene-expression profiles, and selected findings validated using individual shRNAs and/or independent siRNAs. The pooled shRNA library design delivered reproducible data. In all, COPOTE analysis identified 52 probable on-target gene-knockdowns. Knockdown of known oncogenes (KRAS, MYC, SMURF1 and CCNE1) and a tumor suppressor (CDKN2A) showed the expected contrasting effects on cell growth. In addition, the screen corroborated purported roles of PLEKHG2 and MED29 as 19q13 amplicon drivers. Most notably, the analysis also revealed novel possible oncogenic functions of nucleoporin NUP153 (ostensibly by modulating TGFbeta signaling) and Kruppel-like transcription factor KLF5 in pancreatic cancer. Conclusions: By integrating physical and functional genomic data, we were able to simultaneously evaluate many candidate pancreatic cancer genes. Our findings uncover new facets of pancreatic cancer biology, with possible therapeutic implications. More broadly, our study provides a general strategy for the efficient characterization of candidate genes emerging from cancer genome studies.

Description: Background: The specialisation of mammalian cells in time and space requires genes associated with specific pathways and functions to be co-ordinately expressed. Here we have combined a large number of publically available microarray datasets derived from human primary cells and analysed large correlation graphs of these data. Results: Using the network analysis tool BioLayout Express3D we identify robust co-associations of genes expressed in a wide variety of cell lineages. We discuss the biological significance of a number of these associations, in particular the coexpression of key transcription factors with the genes that they are likely to control. Conclusions: We consider the regulation of genes in human primary cells and specifically in the human mononuclear phagocyte system. Of particular note is the fact that these data do not support the identity of putative markers of antigen-presenting dendritic cells, nor classification of M1 and M2 activation states, a current subject of debate within immunological field. We have provided this data resource on the BioGPS web site (www.biogps.org) and on macrophages.com (www.macrophages.com).

Description: Background: Edwardsiella tarda is an enterobacterium which causes edwardsiellosis, a fatal disease of cultured fishes such as red sea bream, eel, and flounder. Preventing the occurrence of E. tarda infection has thus been an important issue in aquaculture. E. tarda has been isolated from other animals and from many environments; however, the relationship between the genotype and evolutionary process of this pathogen is not fully understood. To clarify this relationship, we sequenced and compared the genomes of pathogenic and non-pathogenic E. tarda strains isolated from fish, human, and eel pond using next-generation sequencing technology. Results: Eight strains of E. tarda were sequenced with high accuracy (〉99.9%) with coverages from 50- to 400-fold. The obtained reads were mapped to a public reference genome. By comparing single nucleotide and insertion/deletion polymorphisms, we found that an attenuated strain of E. tarda had a loss-of-function mutation in a gene related to the type III secretion system (T3SS), suggesting that this gene is involved in the virulence of E. tarda. A comprehensive gene comparison indicated that fish pathogenic strains possessed a type VI secretion system (T6SS) and pilus assembly genes in addition to the T3SS. Moreover, we found that an E. tarda strain isolated from red sea bream harbored two pathogenicity islands of T3SS and T6SS, which were absent in other strains. In particular, this T3SS was homologous to the locus of enterocyte effacement (LEE) in enteropathogenic and enterohemorrhagic Escherichia coli. Evolutionary analysis suggested that this locus, here named Et-LEE (E. tarda LEE), was introgressed into the E. tarda genome through horizontal transfer. Conclusions: We found significant differences in the presence/absence of virulence-related genes among E. tarda strains, reflecting their evolutionary relationship. In particular, a single genotype previously proposed for fish-pathogenic strains may be further divided into two subgroups. Furthermore, the current study demonstrated, for the first time, that a fish pathogenic bacterium carried a LEE-like pathogenicity island which was previously reported only in zoonotic pathogenic enterobacteria. These findings will contribute to the exploration of strain-specific drug targets against E. tarda in aquafarms, while also shedding light on the evolution of pathogenesis in enterobacteria.

Description: Background: Drought tolerance is an attribute maintained in plants by cross-talk between multiple and cascading metabolic pathways. Without a sequenced genome available for horse gram, it is difficult to comprehend such complex networks and intercalated genes associated with drought tolerance of horse gram (Macrotyloma uniflorum). Therefore, de novo transcriptome discovery and associated analyses was done for this highly drought tolerant yet under exploited legume to decipher its genetic makeup. Results: Eight samples comprising of shoot and root tissues of two horse gram genotypes (drought-sensitive; M-191 and drought-tolerant; M-249) were used for comparison under control and polyethylene glycol-induced drought stress conditions. Using Illumina sequencing technology, a total of 229,297,896 paired end read pairs were generated and utilized for de novo assembly of horse gram. Significant BLAST hits were obtained for 26,045 transcripts while, 3,558 transcripts had no hits but contained important conserved domains. A total of 21,887 unigenes were identified. SSRs containing sequences covered 16.25% of the transcriptome with predominant tri- and mono-nucleotides (43%). The total GC content of the transcriptome was found to be 43.44%. Under Gene Ontology response to stimulus, DNA binding and catalytic activity was highly expressed during drought stress conditions. Serine/threonine protein kinase was found to dominate in Enzyme Classification while pathways belonging to ribosome metabolism followed by plant pathogen interaction and plant hormone signal transduction were predominant in Kyoto Encyclopedia of Genes and Genomes analysis. Independent search on plant metabolic network pathways suggested valine degradation, gluconeogenesis and purine nucleotide degradation to be highly influenced under drought stress in horse gram. Transcription factors belonging to NAC, MYB-related, and WRKY families were found highly represented under drought stress. qRT-PCR validated the expression profile for 9 out of 10 genes analyzed in response to drought stress. Conclusions: De novo transcriptome discovery and analysis has generated enormous information over horse gram genomics. The genes and pathways identified suggest efficient regulation leading to active adaptation as a basal defense response against drought stress by horse gram. The knowledge generated can be further utilized for exploring other underexploited plants for stress responsive genes and improving plant tolerance.

Description: Background: In contrast to wild species, which have typically evolved phenotypes over long periods of natural selection, domesticates rapidly gained human-preferred agronomic traits in a relatively short-time frame via artificial selection. Under domesticated conditions, many traits can be observed that cannot only be due to environmental alteration. In the case of silkworms, aside from genetic divergence, whether epigenetic divergence played a role in domestication is an unanswered question. The silkworm is still an enigma in that it has two DNA methyltransferases (DNMT1 and DNMT2) but their functionality is unknown. Even in particular the functionality of the widely distributed DNMT1 remains unknown in insects in general. Results: By embryonic RNA interference, we reveal that knockdown of silkworm Dnmt1 caused decreased hatchability, providing the first direct experimental evidence of functional significance of insect Dnmt1. In the light of this fact and those that DNA methylation is correlated with gene expression in silkworms and some agronomic traits in domesticated organisms are not stable, we comprehensively compare silk gland methylomes of 3 domesticated (Bombyx mori) and 4 wild (Bombyx mandarina) silkworms to identify differentially methylated genes between the two. We observed 2-fold more differentiated methylated cytosinces (mCs) in domesticated silkworms as compared to their wild counterparts, suggesting a trend of increasing DNA methylation during domestication. Further study of more domesticated and wild silkworms narrowed down the domesticates' epimutations, and we were able to identify a number of differential genes. One such gene showing demethyaltion in domesticates correspondently displays lower gene expression, and more interestingly, has experienced selective sweep. A methylation-increased gene seems to result in higher expression in domesticates and the function of its Drosophila homolog was previously found to be essential for cell volume regulation, indicating a possible correlation with the enlargement of silk glands in domesticated silkworms. Conclusions: Our results imply epigenetic influences at work during domestication, which gives insight into long time historical controversies regarding acquired inheritance.

Description: Background: It was long assumed that proteins are at least 100 amino acids (AAs) long. Moreover, the detection of short translation products (e.g. coded from small Open Reading Frames, sORFs) is very difficult as the short length makes it hard to distinguish true coding ORFs from ORFs occurring by chance. Nevertheless, over the past few years many such non-canonical genes (with ORFs 〈 100 AAs) have been discovered in different organisms like Arabidopsis thaliana, Saccharomyces cerevisiae, and Drosophila melanogaster. Thanks to advances in sequencing, bioinformatics and computing power, it is now possible to scan the genome in unprecedented scrutiny, for example in a search of this type of small ORFs. Results: Using bioinformatics methods, we performed a systematic search for putatively functional sORFs in the Mus musculus genome. A genome-wide scan detected all sORFs which were subsequently analyzed for their coding potential, based on evolutionary conservation at the AA level, and ranked using a Support Vector Machine (SVM) learning model. The ranked sORFs are finally overlapped with ribosome profiling data, hinting to sORF translation. All candidates are visually inspected using an in-house developed genome browser. In this way dozens of highly conserved sORFs, targeted by ribosomes were identified in the mouse genome, putatively encoding micropeptides. Conclusion: Our combined genome-wide approach leads to the prediction of a comprehensive but manageable set of putatively coding sORFs, a very important first step towards the identification of a new class of bioactive peptides, called micropeptides.

Description: Background: The analysis of polyploid genomes is problematic because homeologous subgenome sequences are closely related. This relatedness makes it difficult to assign individual sequences to the specific subgenome from which they are derived, and hinders the development of polyploid whole genome assemblies. Results: We here present a next-generation sequencing (NGS)-based approach for assignment of subgenome-specific base-identity at sites containing homeolog-specific polymorphisms (HSPs): 'HSP base Assignment using NGS data through Diploid Similarity' (HANDS). We show that HANDS correctly predicts subgenome-specific base-identity at 〉90% of assayed HSPs in the hexaploid bread wheat (Triticum aestivum) transcriptome, thus providing a substantial increase in accuracy versus previous methods for homeolog-specific base assignment. Conclusion: We conclude that HANDS enables rapid and accurate genome-wide discovery of homeolog-specific base-identity, a capability having multiple applications in polyploid genomics.

Description: Background: The identification of the loci and specific alleles underlying variation in quantitative traits is an important goal for evolutionary biologists and breeders. Despite major advancements in genomics technology, moving from QTL to causal alleles remains a major challenge in genetics research. Near-isogenic lines are the ideal raw material for QTL validation, refinement of QTL location and, ultimately, gene discovery. Results: In this study, a population of 75 Arabidopsis thaliana near-isogenic lines was developed from an existing recombinant inbred line (RIL) population derived from a cross between physiologically divergent accessions Kas-1 and Tsu-1. First, a novel algorithm was developed to utilize genome-wide marker data in selecting RILs fully isogenic to Kas-1 for a single chromosome. Seven such RILs were used in 2 generations of crossing to Tsu-1 to create BC1 seed. BC1 plants were genotyped with SSR markers so that lines could be selected that carried Kas-1 introgressions, resulting in a population carrying chromosomal introgressions spanning the genome. BC1 lines were genotyped with 48 genome-wide SSRs to identify lines with a targeted Kas-1 introgression and the fewest genomic introgressions elsewhere. 75 such lines were selected and genotyped at an additional 41 SNP loci and another 930 tags using 2b-RAD genotyping by sequencing. The final population carried an average of 1.35 homozygous and 2.49 heterozygous introgressions per line with average introgression sizes of 5.32 and 5.16 Mb, respectively. In a simple case study, we demonstrate the advantage of maintaining heterozygotes in our library whereby fine-mapping efforts are conducted simply by self-pollination. Crossovers in the heterozygous interval during this single selfing generation break the introgression into smaller, homozygous fragments (sub-NILs). Additionally, we utilize a homozygous NIL for validation of a QTL underlying stomatal conductance, a low heritability trait. Conclusions: The present results introduce a new and valuable resource to the Brassicaceae research community that enables rapid fine-mapping of candidate loci in parallel with QTL validation. These attributes along with dense marker coverage and genome-wide chromosomal introgressions make this population an ideal starting point for discovery of genes underlying important complex traits of agricultural and ecological significance.

Description: Background: Advances in high-throughput sequencing technology have yielded a large number of publicly available vertebrate genomes, many of which are selected for inclusion in NCBI's RefSeq project and subsequently processed by NCBI's eukaryotic annotation pipeline. Genome annotation results are affected by differences in available support evidence and may be impacted by annotation pipeline software changes over time. The RefSeq project has not previously assessed annotation trends across organisms or over time. To address this deficiency, we have developed a comparative protocol which integrates analysis of annotated protein-coding regions across a data set of vertebrate orthologs in genomic sequence coordinates, protein sequences, and protein features. Results: We assessed an ortholog dataset that includes 34 annotated vertebrate RefSeq genomes including human. We confirm that RefSeq protein-coding gene annotations in mammals exhibit considerable similarity. Over 50% of the orthologous protein-coding genes in 20 organisms are supported at the level of splicing conservation with at least three selected reference genomes. Approximately 7,500 ortholog sets include at least half of the analyzed organisms, show highly similar sequence and conserved splicing, and may serve as a minimal set of mammalian "core proteins" for initial assessment of new mammalian genomes. Additionally, 80% of the proteins analyzed pass a suite of tests to detect proteins that lack splicing conservation and have unusual sequence or domain annotation. We use these tests to define an annotation quality metric that is based directly on the annotated proteins thus operates independently of other quality metrics such as availability of transcripts or assembly quality measures. Results are available on the RefSeq FTP site [1]. Conclusions: Our multi-factored analysis demonstrates a high level of consistency in RefSeq protein representation among vertebrates. We find that the majority of the RefSeq vertebrate proteins for which we have calculated orthology are good as measured by these metrics. The process flow described provides specific information on the scope and degree of conservation for the analyzed protein sequences and annotations and will be used to enrich the quality of RefSeq records by identifying targets for further improvement in the computational annotation pipeline, and by flagging specific genes for manual curation.

Description: Background: The tropical ascomycete Trichoderma reesei (Hypocrea jecorina) represents one of the most efficient plant cell wall degraders. Regulation of the enzymes required for this process is affected by nutritional signals as well as other environmental signals including light. Results: Our transcriptome analysis of strains lacking the photoreceptors BLR1 and BLR2 as well as ENV1 revealed a considerable increase in the number of genes showing significantly different transcript levels in light and darkness compared to wild-type. We show that members of all glycoside hydrolase families can be subject to light dependent regulation, hence confirming nutrient utilization including plant cell wall degradation as a major output pathway of light signalling. In contrast to N. crassa, photoreceptor mediated regulation of carbon metabolism in T. reesei occurs primarily by BLR1 and BLR2 via their positive effect on induction of env1 transcription, rather than by a presumed negative effect of ENV1 on the function of the BLR complex. Nevertheless, genes consistently regulated by photoreceptors in N. crassa and T. reesei are significantly enriched in carbon metabolic functions. Hence, different regulatory mechanisms are operative in these two fungi, while the light dependent regulation of plant cell wall degradation appears to be conserved.Analysis of growth on different carbon sources revealed that the oxidoreductive D-galactose and pentose catabolism is influenced by light and ENV1. Transcriptional regulation of the target enzymes in these pathways is enhanced by light and influenced by ENV1, BLR1 and/or BLR2. Additionally we detected an ENV1-regulated genomic cluster of 9 genes including the D-mannitol dehydrogenase gene lxr1, with two genes of this cluster showing consistent regulation in N. crassa. Conclusions: We show that one major output pathway of light signalling in Trichoderma reesei is regulation of glycoside hydrolase genes and the degradation of hemicellulose building blocks. Targets of ENV1 and BLR1/BLR2 are for the most part distinct and indicate individual functions for ENV1 and the BLR complex besides their postulated regulatory interrelationship.

Description: Background: Adult human fibroblasts grown in low oxygen and with FGF2 supplementation have the capacity to tip the healing outcome of skeletal muscle injury -- by favoring regeneration response in vivo over scar formation. Here, we compare the transcriptomes of control adult human dermal fibroblasts and induced regeneration-competent (iRC) fibroblasts to identify transcriptional changes that may be related to their regeneration competence. Results: We identified a unique gene-expression profile that characterizes FGF2-induced iRC fibroblast phenotype. Significantly differentially expressed genes due to FGF2 treatment were identified and analyzed to determine overrepresented Gene Ontology terms. Genes belonging to extracellular matrix components, adhesion molecules, matrix remodelling, cytoskeleton, and cytokines were determined to be affected by FGF2 treatment. Conclusions: Transcriptome analysis comparing control adult human fibroblasts with FGF2-treated fibroblasts identified functional groups of genes that reflect transcriptional changes potentially contributing to their regeneration competence. This comparative transcriptome analysis should contribute new insights into genes that characterize cells with greater regenerative potential.

Description: Background: Analyzing the RNA pool or transcription start sites requires effective means to convert RNA into cDNA libraries for digital expression counting. With current high-speed sequencers, it is necessary to flank the cDNAs with specific adapters. Adding template-switching oligonucleotides to reverse transcription reactions is the most commonly used approach when working with very small quantities of RNA even from single cells. Results: Here we compared the performance of DNA-RNA, DNA-LNA and DNA oligonucleotides in template-switching during nanoCAGE library preparation. Test libraries from rat muscle and HeLa cell RNA were prepared in technical triplicates and sequenced for comparison of the gene coverage and distribution of the reads within transcripts. The DNA-RNA oligonucleotide showed the highest specificity for capped 5[prime] ends of mRNA, whereas the DNA-LNA provided similar gene coverage with more reads falling within exons. Conclusions: While confirming the cap-specific preference of DNA-RNA oligonucleotides in template-switching reactions, our data indicate that DNA-LNA hybrid oligonucleotides could potentially find other applications in random RNA sequencing.

Description: Background: Agaricus bisporus is commercially grown on compost, in which the available carbon sources consist mainly of plant-derived polysaccharides that are built out of various different constituent monosaccharides. The major constituent monosaccharides of these polysaccharides are glucose, xylose, and arabinose, while smaller amounts of galactose, glucuronic acid, rhamnose and mannose are also present. Results: In this study, genes encoding putative enzymes from carbon metabolism were identified and their expression was studied in different growth stages of A. bisporus. We correlated the expression of genes encoding plant and fungal polysaccharide modifying enzymes identified in the A. bisporus genome to the soluble carbohydrates and the composition of mycelium grown compost, casing layer and fruiting bodies. Conclusions: The compost grown vegetative mycelium of A. bisporus consumes a wide variety of monosaccharides. However, in fruiting bodies only hexose catabolism occurs, and no accumulation of other sugars was observed. This suggests that only hexoses or their conversion products are transported from the vegetative mycelium to the fruiting body, while the other sugars likely provide energy for growth and maintenance of the vegetative mycelium. Clear correlations were found between expression of the genes and composition of carbohydrates. Genes encoding plant cell wall polysaccharide degrading enzymes were mainly expressed in compost-grown mycelium, and largely absent in fruiting bodies. In contrast, genes encoding fungal cell wall polysaccharide modifying enzymes were expressed in both fruiting bodies and vegetative mycelium, but different gene sets were expressed in these samples.

Description: Background: Chrysanthemum is one of the most important ornamental crops in the world and drought stress seriously limits its production and distribution. In order to generate a functional genomics resource and obtain a deeper understanding of the molecular mechanisms regarding chrysanthemum responses to dehydration stress, we performed large-scale transcriptome sequencing of chrysanthemum plants under dehydration stress using the Illumina sequencing technology. Results: Two cDNA libraries constructed from mRNAs of control and dehydration-treated seedlings were sequenced by Illumina technology. A total of more than 100 million reads were generated and de novo assembled into 98,180 unique transcripts which were further extensively annotated by comparing their sequencing to different protein databases. Biochemical pathways were predicted from these transcript sequences. Furthermore, we performed gene expression profiling analysis upon dehydration treatment in chrysanthemum and identified 8,558 dehydration-responsive unique transcripts, including 307 transcription factors and 229 protein kinases and many well-known stress responsive genes. Gene ontology (GO) term enrichment and biochemical pathway analyses showed that dehydration stress caused changes in hormone response, secondary and amino acid metabolism, and light and photoperiod response. These findings suggest that drought tolerance of chrysanthemum plants may be related to the regulation of hormone biosynthesis and signaling, reduction of oxidative damage, stabilization of cell proteins and structures, and maintenance of energy and carbon supply. Conclusions: Our transcriptome sequences can provide a valuable resource for chrysanthemum breeding and research and novel insights into chrysanthemum responses to dehydration stress and offer candidate genes or markers that can be used to guide future studies attempting to breed drought tolerant chrysanthemum cultivars.

Description: Background: The eastern mosquitofish (Gambusia holbrooki) has the potential to become a bioindicator organism of endocrine disrupting chemicals (EDCs) due to its androgen-driven secondary sexual characteristics. However, the lack of molecular information on G. holbrooki hinders its use as a bioindicator coupled with biomarker data. While traditional gene-by-gene approaches provide insight for biomarker development, a holistic analysis would provide more rapid and expansive determination of potential biomarkers. The objective of this study was to develop and utilize a mosquitofish microarray to determine potential biomarkers of subchronic androgen exposure. To achieve this objective, two specific aims were developed: 1) Sequence a G. holbrooki cDNA library, and 2) Use microarray analysis to determine genes that are differentially regulated by subchronic androgen exposure in hepatic tissues of 17beta-trenbolone (TB) exposed adult female G. holbrooki. Results: A normalized library of multiple organs of male and female G. holbrooki was prepared and sequenced by the Illumina GA IIx and Roche 454 XLR70. Over 30,000 genes with e-value 1.5 and 〈 -1.5). Fifteen gene ontology biological processes were enriched by TB exposure (Fisher's Exact Test, p 〈 0.05). The expression levels of 17beta-hydroxysteroid dehydrogenase 3 and zona pellucida glycoprotein 2 were validated by quantitative polymerase chain reaction (qPCR) (Student's t-test, p 〈 0.05). Conclusions: Coupling microarray data with phenotypic changes driven by androgen exposure in mosquitofish is key for developing this organism into a bioindicator for EDCs. Future studies using this array will enhance knowledge of the biology and toxicological response of this species. This work provides a foundation of molecular knowledge and tools that can be used to delve further into understanding the biology of G. holbrooki and how this organism can be used as a bioindicator organism for endocrine disrupting pollutants in the environment.

Description: Background: Small non-coding RNAs (ncRNAs) are important regulators of gene expression in eukaryotes. Previously, only microRNAs (miRNAs) and piRNAs have been identified in the silkworm, Bombyx mori. Furthermore, only ncRNAs (50-500nt) of intermediate size have been systematically identified in the silkworm. Results: Here, we performed a systematic identification and analysis of small RNAs (18-50nt) associated with the Bombyx mori argonaute2 (BmAgo2) protein. Using RIP-seq, we identified various types of small ncRNAs associated with BmAGO2. These ncRNAs showed a multimodal length distribution, with three peaks at ~20nt, ~27nt and ~33nt, which included tRNA-, transposable element (TE)-, rRNA-, snoRNA- and snRNA-derived small RNAs as well as miRNAs and piRNAs. The tRNA-derived fragments (tRFs) were found at an extremely high abundance and accounted for 69.90% of the BmAgo2-associated small RNAs. Northern blotting confirmed that many tRFs were expressed or up-regulated only in the BmNPV-infected cells, implying that the tRFs play a prominent role by binding to BmAgo2 during BmNPV infection. Additional evidence suggested that there are potential cleavage sites on the D, anti-codon and TpsiC loops of the tRNAs. TE-derived small RNAs and piRNAs also accounted for a significant proportion of the BmAgo2-associated small RNAs, suggesting that BmAgo2 could be involved in the maintenance of genome stability by suppressing the activities of transposons guided by these small RNAs. Finally, Northern blotting was also used to confirm the Bombyx 5.8 s rRNA-derived small RNAs, demonstrating that various novel small RNAs exist in the silkworm. Conclusions: Using an RIP-seq method in combination with Northern blotting, we identified various types of small RNAs associated with the BmAgo2 protein, including tRNA-, TE-, rRNA-, snoRNA- and snRNA-derived small RNAs as well as miRNAs and piRNAs. Our findings provide new clues for future functional studies of the role of small RNAs in insect development and evolution.

Description: Background: Release of domesticated strains of fish into nature may pose a threat to wild populations with respect to their evolved genetic structure and fitness. Understanding alterations that have occurred in both physiology and genetics as a consequence of domestication can assist in evaluating the risks posed by introgression of domesticated genomes into wild genetic backgrounds, however the molecular causes of these consequences are currently poorly defined. The present study has examined levels of mRNA in fast-growing pure domesticated (D), slow-growing age-matched pure wild (Wa), slow-growing size-matched pure wild (Ws), and first generation hybrid cross (W/D) rainbow trout (Oncorhynchus mykiss) to investigate the influence of genotype (domesticated vs. wild, and their interactions in hybrids) and developmental stage (age- or size-matched animals) on genetic responses (i.e. dominant vs. recessive) and specific physiological pathways. Results: Highly significant differences in mRNA levels were found between domesticated and wild-type rainbow trout genotypes (321 mRNAs), with many mRNAs in the wild-domesticated hybrid progeny showing intermediate levels. Differences were also found between age-matched and size-matched wild-type trout groups (64 mRNAs), with unique mRNA differences for each of the wild-type groups when compared to domesticated trout (Wa: 114 mRNAs, Ws: 88 mRNAs), illustrating an influence of fish developmental stage affecting findings when used as comparator groups to other genotypes. Analysis of differentially expressed mRNAs (found for both wild-type trout to domesticated comparisons) among the genotypes indicates that 34.8% are regulated consistent with an additive genetic model, whereas 39.1% and 26.1% show a recessive or dominant mode of regulation, respectively. These molecular data are largely consistent with phenotypic data (growth and behavioural assessments) assessed in domesticated and wild trout strains. Conclusions: The present molecular data are concordant with domestication having clearly altered rainbow trout genomes and consequent phenotype from that of native wild populations. Although mainly additive responses were noted in hybrid progeny, the prevalence of dominant and non-additive responses reveals that introgression of domesticated and wild genotypes alters the type of genetic control of mRNA levels from that of wild-type, which may lead to disruption of gene regulation systems important for developing phenotypes for optimal fitness in nature. A clear influence of both fish age and size (developmental stage) on mRNA levels was also noted in this study, which highlights the importance of examining multiple control samples to provide a comprehensive understanding of changes observed between strains possessing differences in growth rate.

Description: Background: Citrus variegated chlorosis (CVC), caused by Xylella fastidiosa, is one the most important citrus diseases, and affects all varieties of sweet orange (Citrus sinensis L. Osb). On the other hand, among the Citrus genus there are different sources of resistance against X. fastidiosa. For these species identifying these defense genes could be an important step towards obtaining sweet orange resistant varieties through breeding or genetic engineering. To assess these genes we made use of mandarin (C. reticulata Blanco) that is known to be resistant to CVC and shares agronomical characteristics with sweet orange. Thus, we investigated the gene expression in Ponkan mandarin at one day after infection with X. fastidiosa, using RNA-seq. A set of genes considered key elements in the resistance was used to confirm its regulation in mandarin compared with the susceptible sweet orange. Results: Gene expression analysis of mock inoculated and infected tissues of Ponkan mandarin identified 667 transcripts repressed and 724 significantly induced in the later. Among the induced transcripts, we identified genes encoding proteins similar to Pattern Recognition Receptors. Furthermore, many genes involved in secondary metabolism, biosynthesis and cell wall modification were upregulated as well as in synthesis of abscisic acid, jasmonic acid and auxin. Conclusions: This work demonstrated that the defense response to the perception of bacteria involves cell wall modification and activation of hormone pathways, which probably lead to the induction of other defense-related genes. We also hypothesized the induction of auxin-related genes indicates that resistant plants initially recognize X. fastidiosa as a necrotrophic pathogen.

Description: Background: Massively parallel sequencing offers the possibility of revolutionizing the study of viral populations by providing ultra deep sequencing (tens to hundreds of thousand fold coverage) of complete viral genomes. However, differentiation of true low frequency variants from sequencing errors remains challenging. Results: We developed a software package, V-Phaser 2, for inferring intrahost diversity within viral populations. This program adds three major new methodologies to the state of the art: a technique to efficiently utilize paired end read data for calling phased variants, a new strategy to represent and infer length polymorphisms, and an in line filter for erroneous calls arising from systematic sequencing artifacts. We have also heavily optimized memory and run time performance. This combination of algorithmic and technical advances allows V-Phaser 2 to fully utilize extremely deep paired end sequencing data (such as generated by Illumina sequencers) to accurately infer low frequency intrahost variants in viral populations in reasonable time on a standard desktop computer. V-Phaser 2 was validated and compared to both QuRe and the original V-Phaser on three datasets obtained from two viral populations: a mixture of eight known strains ofWest Nile Virus (WNV) sequenced on both 454 Titanium and Illumina MiSeq and a mixture of twenty-four known strains of WNV sequenced only on 454 Titanium. V-Phaser 2 outperformed the other two programs in both sensitivity and specificity while using more than five fold less time and memory. Conclusions: We developed V-Phaser 2, a publicly available software tool (V-Phaser 2 can be accessed via: http://www.broadinstitute.org/scientific-community/science/projects/viral-genomics/v-phaser-2 and is freely available for academic use) that enables the efficient analysis of ultra-deep sequencing data produced by common next generation sequencing platforms for viral populations.

Description: Background: Blood--brain barrier (BBB) disruption is an integral feature of numerous neurological disorders. However, there is a relative lack of knowledge regarding the underlying molecular mechanisms of immune-mediated BBB disruption. We have previously shown that CD8 T cells and perforin play critical roles in initiating altered permeability of the BBB in the peptide-induced fatal syndrome (PIFS) model developed by our laboratory. Additionally, despite having indistinguishable CD8 T cell responses, C57BL/6J (B6) mice are highly susceptible to PIFS, exhibiting functional motor deficits, increased astrocyte activation, and severe CNS vascular permeability, while 129S1/SvImJ (129S1) mice remain resistant. Therefore, to investigate the potential role of genetic factors, we performed a comprehensive genetic analysis of (B6 x 129S1) F2 progeny to define quantitative trait loci (QTL) linked to the phenotypic characteristics stated above that mediate CD8 T cell-initiated BBB disruption. Results: Using single nucleotide polymorphism (SNP) markers and a 95% confidence interval, we identified one QTL (PIFS1) on chromosome 12 linked to deficits in motor function (SNP markers rs6292954, rs13481303, rs3655057, and rs13481324, LOD score = 3.3). In addition we identified a second QTL (PIFS2) on chromosome 17 linked to changes in CNS vascular permeability (SNP markers rs6196216 and rs3672065, LOD score = 3.7). Conclusions: The QTL critical intervals discovered have allowed for compilation of a list of candidate genes implicated in regulating functional deficit and CNS vascular permeability. These genes encode for factors that may be potential targets for therapeutic approaches to treat disorders characterized by CD8 T cell-mediated BBB disruption.

Description: Background: Shigella flexneri is the major cause of shigellosis in the developing countries. The O-antigen component of the lipopolysaccharide is one of the key virulence determinants required for the pathogenesis of S. flexneri. The glucosyltransferase and/or acetyltransferase genes responsible for the modification of the O-antigen are encoded by temperate serotype converting bacteriophages present in the S. flexneri genome. Several serotype converting phages have previously been isolated and characterized, however, attempts to isolate a serotype converting phage which encodes the modification genes of serotypes 4a strain have not been successful. Results: In this study, a novel temperate serotype converting bacteriophage SfIV was isolated. Lysogenisation of phage SfIV converted serotype Y strain to serotype 4a. Electron microscopy indicated that SfIV belongs to Myoviridae family. The 39,758 bp genome of phage SfIV encompasses 54 open reading frames (orfs). Protein level comparison of SfIV with other serotype converting phages of S. flexneri revealed that SfIV is similar to phage SfII and SfV. The comparative analysis also revealed that SfIV phage contained five proteins which were not found in any other phages of S. flexneri. These proteins were: a tail fiber assembly protein, two hypothetical proteins with no clear function, and two other unknown proteins which were encoded by orfs present on a moron, that presumably got introduced in SfIV genome from another species via a transposon. These unique proteins of SfIV may play a role in the pathogenesis of the host. Conclusions: This study reports the isolation and complete genome sequence analysis of bacteriophage SfIV. The SfIV phage has a host range significantly different from the other phages of Shigella. Comparative genome analysis identified several proteins unique to SfIV, which may potentially be involved in the survival and pathogenesis of its host. These findings will further our understanding on the evolution of these phages, and will also facilitate studies on development of new phage vectors and therapeutic agents to control infections caused by S. flexneri.

Description: Background: With the price of next generation sequencing steadily decreasing, bacterial genome assembly is now accessible to a wide range of researchers. It is therefore necessary to understand the best methods for generating a genome assembly, specifically, which combination of sequencing and bioinformatics strategies result in the most accurate assemblies. Here, we sequence three E. coli strains on the Illumina MiSeq, Life Technologies Ion Torrent PGM, and Pacific Biosciences RS. We then perform genome assemblies on all three datasets alone or in combination to determine the best methods for the assembly of bacterial genomes. Results: Three E. coli strains -- BL21(DE3), Bal225, and DH5alpha -- were sequenced to a depth of 100x on the MiSeq and Ion Torrent machines and to at least 125x on the PacBio RS. Four assembly methods were examined and compared. The previously published BL21(DE3) genome [GenBank:AM946981.2], allowed us to evaluate the accuracy of each of the BL21(DE3) assemblies. BL21(DE3) PacBio-only assemblies resulted in a 90% reduction in contigs versus short read only assemblies, while N50 numbers increased by over 7-fold. Strikingly, the number of SNPs in PacBio-only assemblies were less than half that seen with short read assemblies (~20 SNPs vs. ~50 SNPs) and indels also saw dramatic reductions (~2 indel 〉5 bp in PacBio-only assemblies vs. ~12 for short-read only assemblies). Assemblies that used a mixture of PacBio and short read data generally fell in between these two extremes. Use of PacBio sequencing reads also allowed us to call covalent base modifications for the three strains. Each of the strains used here had a known covalent base modification genotype, which was confirmed by PacBio sequencing. Conclusion: Using data generated solely from the Pacific Biosciences RS, we were able to generate the most complete and accurate de novo assemblies of E. coli strains. We found that the addition of other sequencing technology data offered no improvements over use of PacBio data alone. In addition, the sequencing data from the PacBio RS allowed for sensitive and specific calling of covalent base modifications.

Description: Background: The advent of next-generation high-throughput technologies has revolutionized whole genome sequencing, yet some experiments require sequencing only of targeted regions of the genome from a very large number of samples. These regions can be amplified by PCR and sequenced by next-generation methods using a multidimensional pooling strategy. However, there is at present no available generalized tool for the computational analysis of target-enriched NGS data from multidimensional pools. Results: Here we present InsertionMapper, a pipeline tool for the identification of targeted sequences from multidimensional high throughput sequencing data. InsertionMapper consists of four independently working modules: Data Preprocessing, Database Modeling, Dimension Deconvolution and Element Mapping. We illustrate InsertionMapper with an example from our project 'New reverse genetics resources for maize', which aims to sequence-index a collection of 15,000 independent insertion sites of the transposon Ds in maize. Identified sequences are validated by PCR assays. This pipeline tool is applicable to similar scenarios requiring analysis of the tremendous output of short reads produced in NGS sequencing experiments of targeted genome sequences. Conclusions: InsertionMapper is proven efficacious for the identification of target-enriched sequences from multidimensional high throughput sequencing data. With adjustable parameters and experiment configurations, this tool can save great computational effort to biologists interested in identifying their sequences of interest within the huge output of modern DNA sequencers. InsertionMapper is freely accessible at https://sourceforge.net/p/insertionmapper and http://bo.csam.montclair.edu/du/insertionmapper.

Description: Background: Generation of large mate-pair libraries is necessary for de novo genome assembly but the procedure is complex and time-consuming. Furthermore, in some complex genomes, it is hard to increase the N50 length even with large mate-pair libraries, which leads to low transcript coverage. Thus, it is necessary to develop other simple scaffolding approaches, to at least solve the elongation of transcribed fragments. Results: We describe L_RNA_scaffolder, a novel genome scaffolding method that uses long transcriptome reads to order, orient and combine genomic fragments into larger sequences. To demonstrate the accuracy of the method, the zebrafish genome was scaffolded. With expanded human transcriptome data, the N50 of human genome was doubled and L_RNA_scaffolder out-performed most scaffolding results by existing scaffolders which employ mate-pair libraries. In these two examples, the transcript coverage was almost complete, especially for long transcripts. We applied L_RNA_scaffolder to the highly polymorphic pearl oyster draft genome and the gene model length significantly increased. Conclusions: The simplicity and high-throughput of RNA-seq data makes this approach suitable for genome scaffolding. L_RNA_scaffolder is available at http://www.fishbrowser.org/software/L_RNA_scaffolder.

Description: Background: The availability of mitochondrial genomes has allowed for the resolution of numerous questions regarding the evolutionary history of fungi and other eukaryotes. In the Gibberella fujikuroi species complex, the exact relationships among the so-called "African", "Asian" and "American" Clades remain largely unresolved, irrespective of the markers employed. In this study, we considered the feasibility of using mitochondrial genes to infer the phylogenetic relationships among Fusarium species in this complex. The mitochondrial genomes of representatives of the three Clades (Fusarium circinatum, F. verticillioides and F. fujikuroi) were characterized and we determined whether or not the mitochondrial genomes of these fungi have value in resolving the higher level evolutionary relationships in the complex. Results: Overall, the mitochondrial genomes of the three species displayed a high degree of synteny, with all the genes (protein coding genes, unique ORFs, ribosomal RNA and tRNA genes) in identical order and orientation, as well as introns that share similar positions within genes. The intergenic regions and introns generally contributed significantly to the size differences and diversity observed among these genomes. Phylogenetic analysis of the concatenated protein-coding dataset separated members of the Gibberella fujikuroi complex from other Fusarium species and suggested that F. fujikuroi ("Asian" Clade) is basal in the complex. However, individual mitochondrial gene trees were largely incongruent with one another and with the concatenated gene tree, because six distinct phylogenetic trees were recovered from the various single gene datasets. Conclusion: The mitochondrial genomes of Fusarium species in the Gibberella fujikuroi complex are remarkably similar to those of the previously characterized Fusarium species and Sordariomycetes. Despite apparently representing a single replicative unit, all of the genes encoded on the mitochondrial genomes of these fungi do not share the same evolutionary history. This incongruence could be due to biased selection on some genes or recombination among mitochondrial genomes. The results thus suggest that the use of individual mitochondrial genes for phylogenetic inference could mask the true relationships between species in this complex.

Description: Background: Buffalograss [Buchloe dactyloides (Nutt.) Engel. syn. Bouteloua dactyloides (Nutt.) Columbus] is a United States native turfgrass species that requires less irrigation, fungicides and pesticides compared to more commonly used turfgrass species. In areas where water is limited, interest in this grass species for lawns is increasing. While several buffalograss cultivars have been developed through buffalograss breeding, the timeframe for new cultivar development is long and is limited by a lack of useful genetic resources. Two high throughput next-generation sequencing techniques were used to increase the genomic resources available for buffalograss. Results: Total RNA was extracted and purified from leaf samples of two buffalograss cultivars. '378' and 'Prestige' cDNA libraries were subjected to high throughput sequencing on the Illumina GA and Roche 454 Titanium FLX sequencing platforms. The 454 platform (3 samples) produced 1,300,885 reads and the Illumina platform (12 samples) generated approximately 332 million reads. The multiple k-mer technique for de novo assembly using Velvet and Oases was applied. A total of 121,288 contigs were assembled that were similar to previously reported Ensembl commelinid sequences. Original Illumina reads were also mapped to the high quality assembly to estimate expression levels of buffalograss transcripts. There were a total of 325 differentially expressed genes between the two buffalograss cultivars. A glycosyl transferase, serine threonine kinase, and nb-arc domain containing transcripts were among those differentially expressed between the two cultivars. These genes have been previously implicated in defense response pathways and may in part explain some of the performance differences between 'Prestige' and '378'. Conclusions: To date, this is the first high throughput sequencing experiment conducted on buffalograss. In total, 121,288 high quality transcripts were assembled, significantly expanding the limited genetic resources available for buffalograss genetic studies. Additionally, 325 differentially expressed sequences were identified which may contribute to performance or morphological differences between 'Prestige' and '378' buffalograss cultivars.

Description: Background: Polyadenylation is a key regulatory step in eukaryotic gene expression and one of the major contributors of transcriptome diversity. Aberrant polyadenylation often associates with expression defects and leads to human diseases. Results: To better understand global polyadenylation regulation, we have developed a polyadenylation sequencing (PA-seq) approach. By profiling polyadenylation events in 13 human tissues, we found that alternative cleavage and polyadenylation (APA) is prevalent in both protein-coding and noncoding genes. In addition, APA usage, similar to gene expression profiling, exhibits tissue-specific signatures and is sufficient for determining tissue origin. A 3[prime] untranslated region shortening index (USI) was further developed for genes with tandem APA sites. Strikingly, the results showed that different tissues exhibit distinct patterns of shortening and/or lengthening of 3[prime] untranslated regions, suggesting the intimate involvement of APA in establishing tissue or cell identity. Conclusions: This study provides a comprehensive resource to uncover regulated polyadenylation events in human tissues and to characterize the underlying regulatory mechanism.

Description: Background: Phakopsora pachyrhizi, the causal agent responsible for soybean rust, is among the top hundred most virulent plant pathogens and can cause soybean yield losses of up to 80% when appropriate conditions are met. We used mRNA-Seq by Illumina to analyze pathogen transcript abundance at 15 seconds (s), 7 hours (h), 48 h, and 10 days (d) after inoculation (ai) of susceptible soybean leaves with P. pachyrhizi to gain new insights into transcript abundance in soybean and the pathogen at specific time-points during the infection including the uredinial stage. Results: Over three million five hundred thousand sequences were obtained for each time-point. Energy, nucleotide metabolism, and protein synthesis are major priorities for the fungus during infection and development as indicated by our transcript abundance studies. At all time-points, energy production is a necessity for P. pachyrhizi, as indicated by expression of many transcripts encoding enzymes involved in oxidative phosphorylation and carbohydrate metabolism (glycolysis, glyoxylate and dicarboxylate, pentose phosphate, pyruvate). However, at 15 sai, transcripts encoding enzymes involved in ATP production were highly abundant in order to provide enough energy for the spore to germinate, as observed by the expression of many transcripts encoding proteins involved in electron transport. At this early time-point, transcripts encoding proteins involved in RNA synthesis were also highly abundant, more so than transcripts encoding genes involved in DNA and protein synthesis. At 7 hai, shortly after germination during tube elongation and penetration, transcripts encoding enzymes involved in deoxyribonucleotide and DNA synthesis were highly abundant. At 48 hai, transcripts encoding enzymes involved in amino acid metabolism were highly abundant to provide for increased protein synthesis during haustoria maturation. During sporulation at 10 dai, the fungus still required carbohydrate metabolism, but there also was increased expression of transcripts encoding enzymes involved in fatty acid metabolism. Conclusion: This information provides insight into molecular events and their timing throughout the life cycle of the P. pachyrhizi, and it may be useful in the development of new methods of broadening resistance of soybean to soybean rust.

Description: Background: Microsatellite DNA is one of many powerful genetic markers used for the construction of genetic linkage maps and the study of population genetics. The biological databases in public domain hold vast numbers of microsatellite sequences for many organisms including fishes. The microsatellite data available in these data sources were extracted and managed into a database that facilitates sequences analysis and browsing relevant information. The system also helps to design primer sequences for flanking regions of repeat loci for PCR identification of polymorphism within populations.Description: FishMicrosat is a database of microsatellite sequences of fishes and shellfishes that includes important aquaculture species such as Lates calcarifer, Ctenopharyngodon idella, Hypophthalmichthys molitrix, Penaeus monodon, Labeo rohita, Oreochromis niloticus, Fenneropenaeus indicus and Macrobrachium rosenbergii. The database contains 4398 microsatellite sequences of 41 species belonging to 15 families from the Indian subcontinent. GenBank of NCBI was used as a prime data source for developing the database. The database presents information about simple and compound microsatellites, their clusters and locus orientation within sequences. The database has been integrated with different tools in a web interface such as primer designing, locus finding, mapping repeats, detecting similarities among sequences across species, and searching using motifs and keywords. In addition, the database has the ability to browse information on the top 10 families and the top 10 species, through record overview. Conclusions: FishMicrosat database is a useful resource for fish and shellfish microsatellite analyses and locus identification across species, which has important applications in population genetics, evolutionary studies and genetic relatedness among species. The database can be expanded further to include the microsatellite data of fishes and shellfishes from other regions and available information on genome sequencing project of species of aquaculture importance.

Description: Background: It has been demonstrated that the umbilical cord matrix, represented by the Wharton's Jelly (WJ), contains a great number of mesenchymal stem cells (MSCs), characterized by the expression of specific MSCs markers, shared by both human and animal models. The easy access to massive WJ amount makes it an attractive source of MSCs for cell-based therapies. However, as in other stem cell models, a deeper investigation of WJ-derived MSCs (WJ-MSCs) biological properties, probably modulated by their prolonged expansion and fast growth abilities, is required before their use in clinical settings. In this context, in order to analyze specific gene expression modifications occurring in WJ-MSCs, along with their culture prolongation, we investigated the transcriptomic profiles of WJ-MSCs after 4 and 12 passages of in vitro expansion by microarray analysis. Results: Hierarchical clustering analysis of the data set originated from a total of 6 experiments revealed that in vitro expansion of WJ-MSCs up to 12 passages promote selective over-expression of 157 genes and down-regulation of 440 genes compared to the 4th passage. IPA software analysis of the biological functions related to the identified sets of genes disclosed several transcripts related to inflammatory and cell stress response, cell proliferation and maturation, and apoptosis. Conclusions: Taken together, these modifications may lead to an impairment of both cell expansion ability and resistance to apoptosis, two hallmarks of aging cells. In conclusion, results provided by the present study suggest the need to develop novel culture protocols able to preserve stem cell plasticity.

Description: Background: Among the cold-adapted Antarctic notothenioid fishes, the high-latitude bald notothen Pagothenia borchgrevinki is particularly notable as the sole cryopelagic species, exploiting the coldest and iciest waters of the Southern Ocean. Because P. borchgrevinki is a frequent model for investigating notothenioid cold-adaptation and specialization, it is imperative that "omic" tools be developed for this species. In the absence of a sequenced genome, a well annotated reference transcriptome of the bald notothen will serve as a model of gene expression in the coldest and harshest of all polar marine environments, useful for future comparative studies of cold adaptation and thermal responses in polar teleosts and ectotherms. Results: We sequenced and annotated a reference transcriptome for P. borchgrevinki, with added attention to capturing the transcriptional responses to acute and chronic heat exposures. We sequenced by Roche 454 a normalized cDNA library constructed from pooled mRNA encompassing multiple tissues taken from environmental, warm acclimating, and acute heat stressed specimens. The resulting reads were assembled into 42,620 contigs, 17,951 of which could be annotated. We utilized this annotated portion of the reference transcriptome to map short Illumina reads sequenced from the gill and liver of environmental specimens, and also compared the gene expression profiles of these two tissue transcriptomes with those from the temperate model fish Danio rerio. From this, we identified a conserved group of 58 GO terms, in which terms related to transcription and its regulation, ubiquitin-protein ligase activity, protein ubiquitination, and protein binding among others are more prevalent in the bald notothen, suggesting the pertinent genes play essential roles in cold temperature functioning. Conclusion: We sequenced multiple tissue transcriptomes from native and heat-exposed experimental specimens of the high Antarctic, cryopelagic notothenioid P. borchgrevinki to construct a reference transcriptome. In a proof of concept, we utilized the annotated reference transcriptome to profile the gene expression patterns of gill and liver, and identified a suite of over and under-represented GO terms when compared to the tropical water zebrafish suggesting these functions may be important for surviving in freezing waters. The transcriptome resource from this study will aid future investigations of cold adaptation and thermal response of polar ectothermic species.

Description: Background: Mitochondrial genomic (mitogenomic) reorganizations are rarely found in closely-related animals, yet drastic reorganizations have been found in the Ranoides frogs. The phylogenetic relationships of the three major ranoid taxa (Natatanura, Microhylidae, and Afrobatrachia) have been problematic, and mitogenomic information for afrobatrachians has not been available. Several molecular models for mitochondrial (mt) gene rearrangements have been proposed, but observational evidence has been insufficient to evaluate them. Furthermore, evolutionary trends in rearranged mt genes have not been well understood. To gain molecular and phylogenetic insights into these issues, we analyzed the mt genomes of four afrobatrachian species (Breviceps adspersus, Hemisus marmoratus, Hyperolius marmoratus, and Trichobatrachus robustus) and performed molecular phylogenetic analyses. Furthermore we searched for two evolutionary patterns expected in the rearranged mt genes of ranoids. Results: Extensively reorganized mt genomes having many duplicated and rearranged genes were found in three of the four afrobatrachians analyzed. In fact, Breviceps has the largest known mt genome among vertebrates. Although the kinds of duplicated and rearranged genes differed among these species, a remarkable gene rearrangement pattern of non-tandemly copied genes situated within tandemly-copied regions was commonly found. Furthermore, the existence of concerted evolution was observed between non-neighboring copies of triplicated 12S and 16S ribosomal RNA regions. Conclusions: Phylogenetic analyses based on mitogenomic data support a close relationship between Afrobatrachia and Microhylidae, with their estimated divergence 100 million years ago consistent with present-day endemism of afrobatrachians on the African continent. The afrobatrachian mt data supported the first tandem and second non-tandem duplication model for mt gene rearrangements and the recombination-based model for concerted evolution of duplicated mt regions. We also showed that specific nucleotide substitution and compositional patterns expected in duplicated and rearranged mt genes did not occur, suggesting no disadvantage in employing these genes for phylogenetic inference.

Description: Background: Comparison and classification of metagenome samples is one of the major tasks in the study of microbial communities of natural environments or niches on human bodies. Bioinformatics methods play important roles on this task, including 16S rRNA gene analysis and some alignment-based or alignment-free methods on metagenomic data. Alignment-free methods have the advantage of not depending on known genome annotations and therefore have high potential in studying complicated microbiomes. However, the existing alignment-free methods are all based on unsupervised learning strategy (e.g., PCA or hierarchical clustering). These types of methods are powerful in revealing major similarities and grouping relations between microbiome samples, but cannot be applied for discriminating predefined classes of interest which might not be the dominating assortment in the data. Supervised classification is needed in the latter scenario, with the goal of classifying samples into predefined classes and finding the features that can discriminate the classes. The effectiveness of supervised classification with alignment-based features on metagenomic data have been shown in some recent studies. The application of alignment-free supervised classification methods on metagenome data has not been well explored yet. Results: We developed a method for this task using k-tuple frequencies as features counted directly from metagenome short reads and the R-SVM (Recursive SVM) for feature selection and classification. We tested our method on a simulation dataset, a real dataset composed of several known genomes, and a real metagenome NGS short reads dataset. Experiments on simulated data showed that the method can classify the classes almost perfectly and can recover major sequence signatures that distinguish the two classes. On the real human gut metagenome data, the method can discriminate samples of inflammatory bowel disease (IBD) patients from control samples with high accuracy, which cannot be separated when comparing the samples with unsupervised clustering approaches. Conclusions: The proposed alignment-free supervised classification method can perform well in discriminating of metagenomic samples of predefined classes and in selecting characteristic sequence features for the discrimination. This study shows as an example on the feasibility of using metagenome sequence features of microbiomes on human bodies to study specific human health conditions using supervised machine learning methods.

Description: Background: Lymphatic filariasis is a neglected tropical disease afflicting more than 120 million people, while another 1.3 billion people are at risk of infection. The nematode worm Brugia malayi is one of the causative agents of the disease and exists in a mutualistic symbiosis with Wolbachia bacteria. Since extensive lateral gene transfer occurs frequently between Wolbachia and its hosts, we sought to measure the extent of such LGT in B. malayi by whole genome sequencing of Wolbachia-depleted worms. Results: A considerable fraction (at least 115.4-kbp, or 10.6%) of the 1.08-Mbp Wolbachia wBm genome has been transferred to its nematode host and retains high levels of similarity, including 227 wBm genes and gene fragments. Complete open reading frames were transferred for 32 of these genes, meaning they have the potential to produce functional proteins. Moreover, four transfers have evidence of life stage-specific regulation of transcription at levels similar to other nematode transcripts, strengthening the possibility that they are functional. Conclusions: There is extensive and ongoing transfer of Wolbachia DNA to the worm genome and some transfers are transcribed in a stage-specific manner at biologically relevant levels.

Description: Background: Studies of protein association with DNA on a genome wide scale are possible through methods like ChIP-Chip or ChIP-Seq. Massive problems with false positive signals in our own experiments motivated us to revise the standard ChIP-Chip protocol. Analysis of chromosome wide binding of the alternative sigma factor sigma32 in Escherichia coli with this new protocol resulted in detection of only a subset of binding sites found in a previous study by Wade and colleagues. We suggested that the remainder of binding sites detected in the previous study are likely to be false positives. In a recent article the Wade group claimed that our conclusion is wrong and that the disputed sites are genuine sigma32 binding sites. They further claimed that the non-detection of these sites in our study was due to low data quality.Results/discussionWe respond to the criticism of Wade and colleagues and discuss some general questions of ChIP-based studies. We outline why the quality of our data is sufficient to derive meaningful results. Specific points are: (i) the modifications we introduced into the standard ChIP-Chip protocol do not necessarily result in a low dynamic range, (ii) correlation between ChIP-Chip replicates should not be calculated based on the whole data set as done in transcript analysis, (iii) control experiments are essential for identifying false positives. Suggestions are made how ChIP-based methods could be further optimized and which alternative approaches can be used to strengthen conclusions. Conclusion: We appreciate the ongoing discussion about the ChIP-Chip method and hope that it helps other scientist to analyze and interpret their results. The modifications we introduced into the ChIP-Chip protocol are a first step towards reducing false positive signals but there is certainly potential for further optimization. The discussion about the sigma32 binding sites in question highlights the need for alternative approaches and further investigation of appropriate methods for verification.

Description: Background: Verticillium wilt, caused by the fungal pathogen Verticillium dahliae, is the most severe disease in cotton (Gossypium spp.), causing great lint losses worldwide. Disease management could be achieved in the field if genetically improved, resistant plants were used. However, the interaction between V. dahliae and cotton is a complicated process, and its molecular mechanism remains obscure. To understand better the defense response to this pathogen as a means for obtaining more tolerant cultivars, we monitored the transcriptome profiles of roots from resistant plants of G. barbadense cv. Pima90-53 that were challenged with V. dahliae. Results: In all, 46,192 high-quality expressed sequence tags (ESTs) were generated from a full-length cDNA library of G. barbadense. They were clustered and assembled into 23126 unigenes that comprised 2661 contigs and 20465 singletons. Those unigenes were assigned Gene Ontology terms and mapped to 289 KEGG pathways. A total of 3027 unigenes were found to be homologous to known defense-related genes in other plants. They were assigned to the functional classification of plant--pathogen interactions, including disease defenses and signal transduction. The branch of "SA[rightwards arrow]NPR1[rightwards arrow]TGA[rightwards arrow]PR-1[rightwards arrow]Disease resistance" was first discovered in the interaction of cotton--V. dahliae, indicating that this wilt process includes both biotrophic and necrotrophic stages. In all, 4936 genes coding for putative transcription factors (TF) were identified in our library. The most abundant TF family was the NAC group (527), followed by G2-like (440), MYB (372), BHLH (331), bZIP (271) ERF, C3H, and WRKY. We also analyzed the expression of genes involved in pathogen-associated molecular pattern (PAMP) recognition, the activation of effector-triggered immunity, TFs, and hormone biosynthesis, as well as genes that are pathogenesis-related, or have roles in signaling/regulatory functions and cell wall modification. Their differential expression patterns were compared among mock-/inoculated- and resistant/susceptible cotton. Our results suggest that the cotton defense response has significant transcriptional complexity and that large accumulations of defense-related transcripts may contribute to V. dahliae resistance in cotton. Therefore, these data provide a resource for cotton improvement through molecular breeding approaches. Conclusions: This study generated a substantial amount of cotton transcript sequences that are related to defense responses against V. dahliae. These genomics resources and knowledge of important related genes contribute to our understanding of host--pathogen interactions and the defense mechanisms utilized by G. barbadense, a non-model plant system. These tools can be applied in establishing a modern breeding program that uses marker-assisted selections and oligonucleotide arrays to identify candidate genes that can be linked to valuable agronomic traits in cotton, including disease resistance.

Description: Background: One of the challenges in insect chemical ecology is to understand how insect pheromones are synthesised, detected and degraded. Genome wide survey by comparative sequencing and gene specific expression profiling provide rich resources for this challenge. A. ipsilon is a destructive pest of many crops and further characterization of the genes involved in pheromone biosynthesis and transport could offer potential targets for disruption of their chemical communication and for crop protection. Results: Here we report 454 next-generation sequencing of the A. ipsilon pheromone gland transcriptome, identification and expression profiling of genes putatively involved in pheromone production, transport and degradation. A total of 23473 unigenes were obtained from the transcriptome analysis, 86% of which were A. ipsilon specific. 42 transcripts encoded enzymes putatively involved in pheromone biosynthesis, of which 15 were specifically, or mainly, expressed in the pheromone glands at 5 to 120-fold higher levels than in the body. Two transcripts encoding for a fatty acid synthase and a desaturase were highly abundant in the transcriptome and expressed more than 40-fold higher in the glands than in the body. The transcripts encoding for 2 acetyl-CoA carboxylases, 1 fatty acid synthase, 2 desaturases, 3 acyl-CoA reductases, 2 alcohol oxidases, 2 aldehyde reductases and 3 acetyltransferases were expressed at a significantly higher level in the pheromone glands than in the body. 17 esterase transcripts were not gland-specific and 7 of these were expressed highly in the antennae. Seven transcripts encoding odorant binding proteins (OBPs) and 8 encoding chemosensory proteins (CSPs) were identified. Two CSP transcripts (AipsCSP2, AipsCSP8) were highly abundant in the pheromone gland transcriptome and this was confirmed by qRT-PCR. One OBP (AipsOBP6) were pheromone gland-enriched and three OBPs (AipsOBP1, AipsOBP2 and AipsOBP4) were antennal-enriched. Based on these studies we proposed possible A. ipsilon biosynthesis pathways for major and minor sex pheromone components. Conclusions: Our study identified genes potentially involved in sex pheromone biosynthesis and transport in A. ipsilon. The identified genes are likely to play essential roles in sex pheromone production, transport and degradation and could serve as targets to interfere with pheromone release. The identification of highly expressed CSPs and OBPs in the pheromone gland suggests that they may play a role in the binding, transport and release of sex pheromones during sex pheromone production in A. ipsilon and other Lepidoptera insects.

Description: Background: Though most of the transcripts are long non-coding RNAs (lncRNAs), little is known about their functions. lncRNAs usually function through interactions with proteins, which implies the importance of identifying the binding proteins of lncRNAs in understanding the molecular mechanisms underlying the functions of lncRNAs. Only a few approaches are available for predicting interactions between lncRNAs and proteins. In this study, we introduce a new method lncPro. Results: By encoding RNA and protein sequences into numeric vectors, we used matrix multiplication to score each RNA--protein pair. This score can be used to measure the interactions between an RNA--protein pair. This method effectively discriminates interacting and non-interacting RNA--protein pairs and predicts RNA--protein interactions within a given complex. Applying this method on all human proteins, we found that the long non-coding RNAs we collected tend to interact with nuclear proteins and RNA-binding proteins. Conclusions: Compared with the existing approaches, our method shortens the time for training matrix and obtains optimal results based on the model being used. The ability of predicting the associations between lncRNAs and proteins has also been enhanced. Our method provides an idea on how to integrate different information into the prediction process.

Description: Background: Brassica oleracea is a morphologically diverse species in the family Brassicaceae and contains a group of nutrition-rich vegetable crops, including common heading cabbage, cauliflower, broccoli, kohlrabi, kale, Brussels sprouts. This diversity along with its phylogenetic membership in a group of three diploid and three tetraploid species, and the recent availability of genome sequences within Brassica provide an unprecedented opportunity to study intra- and inter-species divergence and evolution in this species and its close relatives.DescriptionWe have developed a comprehensive database, Bolbase, which provides access to the B. oleracea genome data and comparative genomics information. The whole genome of B. oleracea is available, including nine fully assembled chromosomes and 1,848 scaffolds, with 45,758 predicted genes, 13,382 transposable elements, and 3,581 non-coding RNAs. Comparative genomics information is available, including syntenic regions among B. oleracea, Brassica rapa and Arabidopsis thaliana, synonymous (Ks) and non-synonymous (Ka) substitution rates between orthologous gene pairs, gene families or clusters, and differences in quantity, category, and distribution of transposable elements on chromosomes. Bolbase provides useful search and data mining tools, including a keyword search, a local BLAST server, and a customized GBrowse tool, which can be used to extract annotations of genome components, identify similar sequences and visualize syntenic regions among species. Users can download all genomic data and explore comparative genomics in a highly visual setting. Conclusions: Bolbase is the first resource platform for the B. oleracea genome and for genomic comparisons with its relatives, and thus it will help the research community to better study the function and evolution of Brassica genomes as well as enhance molecular breeding research. This database will be updated regularly with new features, improvements to genome annotation, and new genomic sequences as they become available. Bolbase is freely available at http://ocri-genomics.org/bolbase.

Description: Background: Previous whole-genome shotgun bisulfite sequencing experiments showed that DNA cytosine methylation in the honey bee (Apis mellifera) is almost exclusively at CG dinucleotides in exons. However, the most commonly used method, bisulfite sequencing, cannot distinguish 5-methylcytosine from 5-hydroxymethylcytosine, an oxidized form of 5-methylcytosine that is catalyzed by the TET family of dioxygenases. Furthermore, some analysis software programs under-represent non-CG DNA methylation and hydryoxymethylation for a variety of reasons. Therefore, we used an unbiased analysis of bisulfite sequencing data combined with molecular and bioinformatics approaches to distinguish 5-methylcytosine from 5-hydroxymethylcytosine. By doing this, we have performed the first whole genome analyses of DNA modifications at non-CG sites in honey bees and correlated the effects of these DNA modifications on gene expression and alternative mRNA splicing. Results: We confirmed, using unbiased analyses of whole-genome shotgun bisulfite sequencing (BS-seq) data, with both new data and published data, the previous finding that CG DNA methylation is enriched in exons in honey bees. However, we also found evidence that cytosine methylation and hydroxymethylation at non-CG sites is enriched in introns. Using antibodies against 5-hydroxmethylcytosine, we confirmed that DNA hydroxymethylation at non-CG sites is enriched in introns. Additionally, using a new technique, Pvu-seq (which employs the enzyme PvuRts1l to digest DNA at 5-hydroxymethylcytosine sites followed by next-generation DNA sequencing), we further confirmed that hydroxymethylation is enriched in introns at non-CG sites. Conclusions: Cytosine hydroxymethylation at non-CG sites might have more functional significance than previously appreciated, and in honey bees these modifications might be related to the regulation of alternative mRNA splicing by defining the locations of the introns.

Description: Background: The production of enzymes by an industrial strain requires a complex adaption of the bacterial metabolism to the conditions within the fermenter. Regulatory events within the process result in a dynamic change of the transcriptional activity of the genome. This complex network of genes is orchestrated by proteins as well as regulatory RNA elements. Here we present an RNA-Seq based study considering selected phases of an industry-oriented fermentation of Bacillus licheniformis. Results: A detailed analysis of 20 strand-specific RNA-Seq datasets revealed a multitude of transcriptionally active genomic regions. 3314 RNA features encoded by such active loci have been identified and sorted into ten functional classes. The identified sequences include the expected RNA features like housekeeping sRNAs, metabolic riboswitches and RNA switches well known from studies on Bacillus subtilis as well as a multitude of completely new candidates for regulatory RNAs. An unexpectedly high number of 855 RNA features are encoded antisense to annotated protein and RNA genes, in addition to 461 independently transcribed small RNAs. These antisense transcripts contain molecules with a remarkable size range variation from 38 to 6348 base pairs in length. The genome of the type strain B. licheniformis DSM13 was completely reannotated using data obtained from RNA-Seq analyses and from public databases. Conclusion: The hereby generated data-sets represent a solid amount of knowledge on the dynamic transcriptional activities during the investigated fermentation stages. The identified regulatory elements enable research on the understanding and the optimization of crucial metabolic activities during a productive fermentation of Bacillus licheniformis strains.

Description: Background: Crown gall (CG) (Agrobacterium tumefaciens) and the root lesion nematodes (RLNs) (Pratylenchus vulnus) are major challenges faced by the California walnut industry, reducing productivity and increasing the cost of establishing and maintaining orchards. Current nematode control strategies include nematicides, crop rotation, and tolerant cultivars, but these methods have limits. Developing genetic resistance through novel approaches like RNA interference (RNAi) can address these problems. RNAi-mediated silencing of CG disease in walnut (Juglans regia L.) has been achieved previously. We sought to place both CG and nematode resistance into a single walnut rootstock genotype using co-transformation to stack the resistance genes. A. tumefaciens, carrying self-complimentary iaaM and ipt transgenes, and Agrobacterium rhizogenes, carrying a self-complimentary Pv010 gene from P. vulnus, were used as co-transformation vectors. RolABC genes were introduced by the resident T-DNA in the A. rhizogenes Ri-plasmid used as a vector for plant transformation. Pv010 and Pv194 (transgenic control) genes were also transferred separately using A. tumefaciens. To test for resistance, transformed walnut roots were challenged with P. vulnus and microshoots were challenged with a virulent strain of A. tumefaciens. Results: Combining the two bacterial strains at a 1:1 rather than 1:3 ratio increased the co-transformation efficiency. Although complete immunity to nematode infection was not observed, transgenic lines yielded up to 79% fewer nematodes per root following in vitro co-culture than untransformed controls. Transgenic line 33-3-1 exhibited complete crown gall control and 32% fewer nematodes. The transgenic plants had thicker, longer roots than untransformed controls possibly due to insertion of rolABC genes. When the Pv010 gene was present in roots with or without rolABC genes there was partial or complete control of RLNs. Transformation using only one vector showed 100% control in some lines. Conclusions: CG and nematode resistance gene stacking controlled CG and RLNs simultaneously in walnuts. Silencing genes encoding iaaM, ipt, and Pv010 decrease CG formation and RLNs populations in walnut. Beneficial plant genotype and phenotype changes are caused by co-transformation using A. tumefaciens and A. rhizogenes strains. Viable resistance against root lesion nematodes in walnut plants may be accomplished in the future using this gene stacking technology.

Description: Background: Second generation sequencing has permitted detailed sequence characterisation at the whole genome level of a growing number of non-model organisms, but the data produced have short read-lengths and biased genome coverage leading to fragmented genome assemblies. The PacBio RS long-read sequencing platform offers the promise of increased read length and unbiased genome coverage and thus the potential to produce genome sequence data of a finished quality containing fewer gaps and longer contigs. However, these advantages come at a much greater cost per nucleotide and with a perceived increase in error-rate. In this investigation, we evaluated the performance of the PacBio RS sequencing platform through the sequencing and de novo assembly of the Potentilla micrantha chloroplast genome. Results: Following error-correction, a total of 28,638 PacBio RS reads were recovered with a mean read length of 1,902bp totalling 54,492,250 nucleotides and representing an average depth of coverage of 320x the chloroplast genome. The dataset covered the entire 154,959bp of the chloroplast genome in a single contig (100% coverage) compared to seven contigs (90.59% coverage) recovered from an Illumina data, and revealed no bias in coverage of GC rich regions. Post-assembly the data were largely concordant with the Illumina data generated and allowed 187 ambiguities in the Illumina data to be resolved. The additional read length also permitted small differences in the two inverted repeat regions to be assigned unambiguously. Conclusions: This is the first report to our knowledge of a chloroplast genome assembled de novo using PacBio sequence data. The PacBio RS data generated here were assembled into a single large contig spanning the P. micrantha chloroplast genome, with a higher degree of accuracy than an Illumina dataset generated at a much greater depth of coverage, due to longer read lengths and lower GC bias in the data. The results we present suggest PacBio data will be of immense utility for the development of genome sequence assemblies containing fewer unresolved gaps and ambiguities and a significantly smaller number of contigs than could be produced using short-read sequence data alone.

Description: Background: Elucidating gut microbiota among gallstone patients as well as the complex bacterial colonization of cholesterol gallstones may help in both the prediction and subsequent lowered risk of cholelithiasis. To this end, we studied the composition of bacterial communities of gut, bile, and gallstones from 29 gallstone patients as well as the gut of 38 normal individuals, examining and analyzing some 299, 217 bacterial 16S rRNA gene sequences from 120 samples. Results: First, as compared with normal individuals, in gallstone patients there were significant (P 〈 0.001) increases of gut bacterial phylum Proteobacteria and decreases of three gut bacterial genera, Faecalibacterium, Lachnospira, and Roseburia. Second, about 70% of gut bacterial operational taxonomic units (OTUs) from gallstone patients were detectable in the biliary tract and bacteria diversity of biliary tract was significantly (P 〈 0.001) higher than that of gut. Third, analysis of the biliary tract core microbiome (represented by 106 bacteria OTUs) among gallstone patients showed that 33.96% (36/106) of constituents can be matched to known bacterial species (15 of which have publicly available genomes). A genome-wide search of MDR, BSH, bG, and phL genes purpotedly associated with the formation of cholesterol gallstones showed that all 15 species with known genomes (e.g., Propionibacterium acnes, Bacteroides vulgates, and Pseudomonas putida) contained at least contained one of the four genes. This finding could potentially provide underlying information needed to explain the association between biliary tract microbiota and the formation of cholesterol gallstones. Conclusions: To the best of our knowledge, this is the first study to discover gut microbiota dysbiosis among gallstone patients, the presence of which may be a key contributor to the complex bacteria community assembly linked with the presence of cholesterol gallstones. Likewise, this study also provides the first large-scale glimpse of biliary tract microbiota potentially associated with cholesterol gallstones. Such a characterization of the biliary tract core microbiome has potentially important biological and medical implications regarding the role of bacteria in the formation cholesterol gallstones.

Description: Background: Multiple myeloma (MM) is a malignant proliferation of plasma B cells. Based on recurrent aneuploidy such as copy number alterations (CNAs), myeloma is divided into two subtypes with different CNA patterns and patient survival outcomes. How aneuploidy events arise, and whether they contribute to cancer cell evolution are actively studied. The large amount of transcriptomic changes resultant of CNAs (dosage effect) pose big challenges for identifying functional consequences of CNAs in myeloma in terms of specific driver genes and pathways. In this study, we hypothesize that gene-wise dosage effect varies as a result from complex regulatory networks that translate the impact of CNAs to gene expression, and studying this variation can provide insights into functional effects of CNAs. Results: We propose gene-wise dosage effect score and genome-wide karyotype plot as tools to measure and visualize concordant copy number and expression changes across cancer samples. We find that dosage effect in myeloma is widespread yet variable, and it is correlated with gene expression level and CNA frequencies in different chromosomes. Our analysis suggests that despite the enrichment of differentially expressed genes between hyperdiploid MM and non-hyperdiploid MM in the trisomy chromosomes, the chromosomal proportion of dosage sensitive genes is higher in the non-trisomy chromosomes. Dosage-sensitive genes are enriched by genes with protein translation and localization functions, and dosage resistant genes are enriched by apoptosis genes. These results point to future studies on differential dosage sensitivity and resistance of pro- and anti-proliferation pathways and their variation across patients as therapeutic targets and prognosis markers. Conclusions: Our findings support the hypothesis that recurrent CNAs in myeloma are selected by their functional consequences. The novel dosage effect score defined in this work will facilitate integration of copy number and expression data for identifying driver genes in cancer genomics studies. The accompanying R code is available at http://www.canevolve.org/dosageEffect/.

Description: Background: The epithelial-to-mesenchymal transition is an important mechanism in cancer metastasis. Although transcription factors including SNAIL, SLUG, and TWIST1 regulate the epithelial-to-mesenchymal transition, other unknown transcription factors could also be involved. Identification of the full complement of transcription factors is essential for a more complete understanding of gene regulation in this process. Chromatin immunoprecipitation-sequencing (ChIP-Seq) technologies have been used to detect genome-wide binding of transcription factors; here, we developed a systematic approach to integrate existing ChIP-Seq and transcriptome data. We scanned multiple transcription factors to investigate their functional impact on the epithelial-to-mesenchymal transition in the human A549 lung adenocarcinoma cell line. Results: Among the transcription factors tested, impact scores identified the forkhead box protein A1 (FOXA1) as the most significant transcription factor in the epithelial-to-mesenchymal transition. FOXA1 physically associates with the promoters of its predicted target genes. Several critical epithelial-to-mesenchymal transition effectors involved in cellular adhesion and cellular communication were identified in the regulatory network of FOXA1, including FOXA2, FGA, FGB, FGG, and FGL1. The implication of FOXA1 in the epithelial-to-mesenchymal transition via its regulatory network indicates that FOXA1 may play an important role in the initiation of lung cancer metastasis. Conclusions: We identified FOXA1 as a potentially important transcription factor and negative regulator in the initial stages of lung cancer metastasis. FOXA1 may modulate the epithelial-to-mesenchymal transition via its transcriptional regulatory network. Further, this study demonstrates how ChIP-Seq and expression data could be integrated to delineate the impact of transcription factors on a specific biological process.

Description: Background: Variation of gene expression can lead to phenotypic variation and have therefore been assumed to contribute the diversity of wine yeast (Saccharomyces cerevisiae) properties. However, the molecular bases of this variation of gene expression are unknown. We addressed these questions by carrying out an integrated genetical-genomic study in fermentation conditions. We report here quantitative trait loci (QTL) mapping based on expression profiling in a segregating population generated by a cross between a derivative of the popular wine strain EC1118 and the laboratory strain S288c. Results: Most of the fermentation traits studied appeared to be under multi-allelic control. We mapped five phenotypic QTLs and 1465 expression QTLs. Several expression QTLs overlapped in hotspots. Among the linkages unraveled here, several were associated with metabolic processes essential for wine fermentation such as glucose sensing or nitrogen and vitamin metabolism. Variations affecting the regulation of drug detoxification and export (TPO1, PDR12 or QDR2) were linked to variation in four genes encoding transcription factors (PDR8, WAR1, YRR1 and HAP1). We demonstrated that the allelic variation of WAR1 and TPO1 affected sorbic and octanoic acid resistance, respectively. Moreover, analysis of the transcription factors phylogeny suggests they evolved with a specific adaptation of the strains to wine fermentation conditions. Unexpectedly, we found that the variation of fermentation rates was associated with a partial disomy of chromosome 16. This disomy resulted from the well known 8--16 translocation. Conclusions: This large data set made it possible to decipher the effects of genetic variation on gene expression during fermentation and certain wine fermentation properties Our findings shed a new light on the adaptation mechanisms required by yeast to cope with the multiple stresses generated by wine fermentation. In this context, the detoxification and export systems appear to be of particular importance, probably due to nitrogen starvation. Furthermore, we show that the well characterized 8--16 translocation located in SSU1, which is associated with sulfite resistance, can lead to a partial chromosomic amplification in the progeny of strains that carry it, greatly improving fermentation kinetics. This amplification has been detected among other wine yeasts.

Description: Background: Amur ide (Leuciscus waleckii) is an economically and ecologically important cyprinid species in Northern Asia. The Dali Nor population living in the soda lake Dali Nor can adapt the extremely high alkalinity, providing us a valuable material to understand the adaptation mechanism against extreme environmental stress in teleost. Results: In this study, we generated high-throughput RNA-Seq data from three tissues gill, liver and kidney of L. waleckii living in the soda lake Dali Nor and the fresh water lake Ganggeng Nor, then performed parallel comparisons of three tissues. Our results showed that out of assembled 64,603 transcript contigs, 28,391 contigs had been assigned with a known function, corresponding to 20,371 unique protein accessions. We found 477, 2,761 and 3,376 differentially expressed genes (DEGs) in the gill, kidney, and liver, respectively, of Dali Nor population compared to Ganggeng Nor population with FDR = 2. Further analysis revealed that well-known functional categories of genes and signaling pathway, which are associated with stress response and extreme environment adaptation, have been significantly enriched, including the functional categories of "response to stimulus", "transferase activity", "transporter activity" and "oxidoreductase activity", and signaling pathways of "mTOR signaling", "EIF2 signaling", "superpathway of cholesterol biosynthesis". We also identified significantly DEGs encoding important modulators on stress adaptation and tolerance, including carbonic anhydrases, heat shock proteins, superoxide dismutase, glutathione S-transferases, aminopeptidase N, and aminotransferases. Conclusions: Overall, this study demonstrated that transcriptome changes in L. waleckii played a role in adaptation to complicated environmental stress in the highly alkalized Dali Nor lake. The results set a foundation for further analyses on alkaline-responsive candidate genes, which help us understand teleost adaptation under extreme environmental stress and ultimately benefit future breeding for alkaline-tolerant fish strains.

Description: Background: Differential processing of the amyloid precursor protein liberates either amyloid-Ss, a causative agent of Alzheimer's disease, or secreted amyloid precursor protein-alpha (sAPPalpha), which promotes neuroprotection, neurotrophism, neurogenesis and synaptic plasticity. The underlying molecular mechanisms recruited by sAPPalpha that underpin these considerable cellular effects are not well elucidated. As these effects are enduring, we hypothesised that regulation of gene expression may be of importance and examined temporally specific gene networks and pathways induced by sAPPalpha in rat hippocampal organotypic slice cultures. Slices were exposed to 1 nM sAPPalpha or phosphate buffered saline for 15 min, 2 h or 24 h and sAPPalpha-associated gene expression profiles were produced for each time-point using Affymetrix Rat Gene 1.0 ST arrays (moderated t-test using Limma: p 〈 0.05, and fold change +/- 1.15). Results: Treatment of organotypic hippocampal slice cultures with 1 nM sAPPalpha induced temporally distinct gene expression profiles, including mRNA and microRNA associated with Alzheimer's disease. Having demonstrated that treatment with human recombinant sAPPalpha was protective against N-methyl d-aspartate-induced toxicity, we next explored the sAPPalpha-induced gene expression profiles. Ingenuity Pathway Analysis predicted that short-term exposure to sAPPalpha elicited a multi-level transcriptional response, including upregulation of immediate early gene transcription factors (AP-1, Egr1), modulation of the chromatin environment, and apparent activation of the constitutive transcription factors CREB and NF-kappaB. Importantly, dynamic regulation of NF-kappaB appears to be integral to the transcriptional response across all time-points. In contrast, medium and long exposure to sAPPalpha resulted in an overall downregulation of gene expression. While these results suggest commonality between sAPPalpha and our previously reported analysis of plasticity-related gene expression, we found little crossover between these datasets. The gene networks formed following medium and long exposure to sAPPalpha were associated with inflammatory response, apoptosis, neurogenesis and cell survival; functions likely to be the basis of the neuroprotective effects of sAPPalpha. Conclusions: Our results demonstrate that sAPPalpha rapidly and persistently regulates gene expression in rat hippocampus. This regulation is multi-level, temporally specific and is likely to underpin the neuroprotective effects of sAPPalpha.

Description: Background: Comparative genomics is a powerful tool to transfer genomic information from model species to related non-model species. Channel catfish (Ictalurus punctatus) is the primary aquaculture species in the United States. Its existing genome resources such as genomic sequences generated from next generation sequencing, BAC end sequences (BES), physical maps, linkage maps, and integrated linkage and physical maps using BES-associated markers provide a platform for comparative genomic analysis between catfish and other model teleost fish species. This study aimed to gain understanding of genome organizations and similarities among catfish and several sequenced teleost genomes using linkage group 8 (LG8) as a pilot study. Results: With existing genome resources, 287 unique genes were identified in LG8. Comparative genome analysis indicated that most of these 287 genes on catfish LG8 are located on two homologous chromosomes of zebrafish, medaka, stickleback, and three chromosomes of green-spotted pufferfish. Large numbers of conserved syntenies were identified. Detailed analysis of the conserved syntenies in relation to chromosome level similarities revealed extensive inter-chromosomal and intra-chromosomal rearrangements during evolution. Of the 287 genes, 35 genes were found to be duplicated in the catfish genome, with the vast majority of the duplications being interchromosomal. Conclusions: Comparative genome analysis is a powerful tool even in the absence of a well-assembled whole genome sequence. In spite of sequence stacking due to low resolution of the linkage and physical maps, conserved syntenies can be identified although the exact gene order and orientation are unknown at present. Through chromosome-level comparative analysis, homologous chromosomes among teleosts can be identified. Syntenic analysis should facilitate annotation of the catfish genome, which in turn, should facilitate functional inference of genes based on their orthology.

Description: Background: Repetitive extragenic palindromic elements (REPs) constitute a group of bacterial genomic repeats known for their high abundance and several roles in host cells physiology. We analyzed the phylogenetic distribution of particular REP classes in genomic sequences of sixty-three bacterial strains belonging to the Pseudomonas fluorescens species complex and ten strains of Stenotrophomonas sp., in order to assess intraspecific REP diversity and to gain insight into long-term REP evolution. Results: Based on proximity to RAYT (REP-associated tyrosine transposase) genes, twenty-two and thirteen unique REP classes were determined in fluorescent pseudomonads and stenotrophomonads, respectively. In stenotrophomonads, REP elements were typically found in tens or a few hundred copies per genome. REPs of fluorescent pseudomonads were generally more numerous, occurring in hundreds or even over a thousand perfect copies of particular REP class per genome. REP sequences showed highly heterogeneous distribution. The abundances of REP classes roughly followed host strains phylogeny, differing markedly among individual clades. High abundances of particular REP classes appeared to depend on the presence of the cognate RAYT gene, and deviations from this state could be attributed to recent or ancient mutations of rayt-flanking REPs, or RAYT loss. RAYTs of both studied bacterial groups are monophyletic, and their cognate REPs show species-specific characteristics, suggesting shared evolutionary history of REPs, RAYTs and their hosts. Conclusions: The results of our large-scale analysis show that REP elements constitute intriguingly dynamic components of genomes of fluorescent pseudomonads and stenotrophomonads, and indicate that REP diversification and proliferation are ongoing processes. High numbers of REPs have apparently been retained during the entire evolutionary time since the establishment of these two bacterial lineages, probably because of their beneficial effect on host long-term fitness. REP elements in these bacteria represent a suitable platform to study the interplay between repeated elements, their mobilizers and host bacterial cells.

Description: Background: Obesity-associated organ-specific pathological states can be ensued from the dysregulation of the functions of the adipose tissues, liver and muscle. However, the influence of genetic differences underlying gross-compositional differences in these tissues is largely unknown. In the present study, the analytical method of ATR-FTIR spectroscopy has been combined with a genetic approach to identify genetic differences responsible for phenotypic alterations in adipose, liver and muscle tissues. Results: Mice from 29 BXD recombinant inbred mouse strains were put on high fat diet and gross-compositional changes in adipose, liver and muscle tissues were measured by ATR-FTIR spectroscopy. The analysis of genotype-phenotype correlations revealed significant quantitative trait loci (QTL) on chromosome 12 for the content of fat and collagen, collagen integrity, and the lipid to protein ratio in adipose tissue and on chromosome 17 for lipid to protein ratio in liver. Using gene expression and sequence information, we suggest Rsad2 (viperin) and Colec11 (collectin-11) on chromosome 12 as potential quantitative trait candidate genes. Rsad2 may act as a modulator of lipid droplet contents and lipid biosynthesis; Colec11 might play a role in apoptopic cell clearance and maintenance of adipose tissue. An increased level of Rsad2 transcripts in adipose tissue of DBA/2J compared to C57BL/6J mice suggests a cis-acting genetic variant leading to differential gene activation. Conclusion: The results demonstrate that the analytical method of ATR-FTIR spectroscopy effectively contributed to decompose the macromolecular composition of tissues that accumulate fat and to link this information with genetic determinants. The candidate genes in the QTL regions may contribute to obesity-related diseases in humans, in particular if the results can be verified in a bigger BXD cohort.

Description: Background: Identifying genes of adaptive significance in a changing environment is a major focus of ecological genomics. Such efforts were restricted, until recently, to researchers studying a small group of model organisms or closely related taxa. With the advent of next generation sequencing (NGS), genomes and transcriptomes of virtually any species are now available for studies of adaptive evolution. We experimentally manipulated temperature conditions for two groups of crimson spotted rainbowfish (Melanotaenia duboulayi) and measured differences in RNA transcription between them. This non-migratory species is found across a latitudinal thermal gradient in eastern Australia and is predicted to be negatively impacted by ongoing environmental and climatic change. Results: Using next generation RNA-seq technologies on an Illumina HiSeq2000 platform, we assembled a de novo transcriptome and tested for differential expression across the treatment groups. Quality of the assembly was high with a N50 length of 1856 bases. Of the 107,749 assembled contigs, we identified 4251 that were differentially expressed according to a consensus of four different mapping and significance testing approaches. Once duplicate isoforms were removed, we were able to annotate 614 up-regulated transfrags and 349 that showed reduced expression in the higher temperature group. Conclusions: Annotated blast matches reveal that differentially expressed genes correspond to critical metabolic pathways previously shown to be important for temperature tolerance in other fish species. Our results indicate that rainbowfish exhibit predictable plastic regulatory responses to temperature stress and the genes we identified provide excellent candidates for further investigations of population adaptation to increasing temperatures.

Description: Background: The Fra-1/AP-1 transcription factor regulates the expression of genes controlling various processes including migration, invasion, and survival as well as extracellular remodeling. We recently demonstrated that loss of Fra-1 leads to exacerbated bleomycin-induced pulmonary fibrosis, accompanied by enhanced expression of various inflammatory and fibrotic genes. To better understand the molecular mechanisms by which Fra-1 confers protection during bleomycin-induced lung injury, genome-wide mRNA expression profiling was performed. Results: We found that Fra-1 regulates gene expression programs that include: 1) several cytokines and chemokines involved in inflammation, 2) several genes involved in the extracellular remodeling and cell adhesion, and 3) several genes involved in programmed cell death. Conclusion: Loss of Fra-1 leads to the enhanced expression of genes regulating inflammation and immune responses and decreased the expression of genes involved in apoptosis, suggesting that this transcription factor distinctly modulates early pro-fibrotic cellular responses.

Description: Background: Viral upper respiratory tract infections are associated with increased colonization by Streptococcus pneumoniae but the mechanisms underlying this relationship are unclear. The objective of this study is to describe a comprehensive picture of the cellular interaction between the adhering bacteria and host cells in the presence or absence of a viral co-infection. Results: Gene expression profiles of Detroit-562 pharyngeal cells, which were either mock infected or infected with human respiratory syncytial virus (RSV) or human parainfluenza virus 3 (HPIV3), were analyzed using human microarrays. Transcription response of S. pneumoniae strain TIGR4 (serotype 4) in the presence of either mock- or viral-infected cells was analyzed by pneumococcal microarray. Significantly regulated genes were identified by both significance analysis of microarray (SAM) and a 〉= 2-fold change ratio cut-off. The adherence of S. pneumoniae to human pharyngeal cells was significantly augmented in the presence of RSV or HPIV3 infection. Global gene expression profiling of the host cells during infection with RSV or HPIV3 revealed increased transcription of carcinoembryonic antigen-related cell adhesion molecules (CEACAM1), CD47, fibronectin, interferon-stimulated genes and many other host cell adhesion molecules. Pneumococci increased transcription of several genes involved in adhesive functions (psaA, pilus islet), choline uptake and incorporation (lic operon), as well as transport and binding. Conclusions: We have identified a core transcriptome that represents the basic machinery required for adherence of pneumococci to D562 cells infected or not infected with a virus. These bacterial genes and cell adhesion molecules can potentially be used to control pneumococcal adherence occurring secondary to a viral infection.

Description: Background: The domestication and subsequent selection by humans to create breeds and biological types of cattle undoubtedly altered the patterning of variation within their genomes. Strong selection to fix advantageous large-effect mutations underlying domesticability, breed characteristics or productivity created selective sweeps in which variation was lost in the chromosomal region flanking the selected allele. Selective sweeps have now been identified in the genomes of many animal species including humans, dogs, horses, and chickens. Here, we attempt to identify and characterise regions of the bovine genome that have been subjected to selective sweeps. Results: Two datasets were used for the discovery and validation of selective sweeps via the fixation of alleles at a series of contiguous SNP loci. BovineSNP50 data were used to identify 28 putative sweep regions among 14 diverse cattle breeds. Affymetrix BOS 1 prescreening assay data for five breeds were used to identify 85 regions and validate 5 regions identified using the BovineSNP50 data. Many genes are located within these regions and the lack of sequence data for the analysed breeds precludes the nomination of selected genes or variants and limits the prediction of the selected phenotypes. However, phenotypes that we predict to have historically been under strong selection include horned-polled, coat colour, stature, ear morphology, and behaviour. Conclusions: The bias towards common SNPs in the design of the BovineSNP50 assay led to the identification of recent selective sweeps associated with breed formation and common to only a small number of breeds rather than ancient events associated with domestication which could potentially be common to all European taurines. The limited SNP density, or marker resolution, of the BovineSNP50 assay significantly impacted the rate of false discovery of selective sweeps, however, we found sweeps in common between breeds which were confirmed using an ultra-high-density assay scored in a small number of animals from a subset of the breeds. No sweep regions were shared between indicine and taurine breeds reflecting their divergent selection histories and the very different environmental habitats to which these sub-species have adapted.

Description: Background: Elucidating the biogeography of bacterial communities on the human body is critical for establishing healthy baselines from which to detect differences associated with disease; however, little is known about the baseline bacterial profiles from various human habitats of healthy Chinese undergraduates. Results: Using parallel barcoded 454 pyrosequencing targeting on the 16S rRNA gene V3 region, the bacterial diversity of the nasopharynx, saliva, dominant hands, and feces were investigated from 10 healthy Chinese junior boarding undergraduates at Zhejiang University. The participants were 21--24 years of age with a body mass index (BMI) 〈 24 kg/m2. A total of 156,717 high-quality pyrosequencing reads were obtained for evaluating bacterial diversity, which represented 29,887 unique phylotypes. The overall taxonomic distribution of the 16S rRNA gene-based amplicons demonstrated that these 4 habitats of the human body harbored distinct microbiota and could be divided into different clusters according to anatomic site, while the established patterns of bacterial diversity followed the human body habitat (feces, hands, saliva, and nasopharynx). Although significant inter-individual variation was observed, the healthy microbiota still shared a large number of phylotypes in each habitat, but not among the four habitats, indicating that a core microbiome existed in each healthy habitat. The vast majority of sequences from these different habitats were classified into different taxonmies that became the predominant bacteria of the healthy microbiota. Conclusions: We first established the framework of microbial communities from four healthy human habitats of the same participants with similar living environments for the Chinese undergraduates. Our data represent an important step for determining the diversity of Chinese healthy microbiota, and can be used for more large-scale studies that focus on the interactions between healthy and diseases states for young Chinese adults in the same age range.

Description: Background: White clover (Trifolium repens L.) is a temperate forage legume with an allotetraploid genome (2n=4x=32) estimated at 1093 Mb. Several linkage maps of various sizes, marker sources and completeness are available, however, no integrated map and marker set has explored consistency of linkage analysis among unrelated mapping populations. Such integrative analysis requires tools for homoeologue matching among populations. Development of these tools provides for a consistent framework map of the white clover genome, and facilitates in silico alignment with the model forage legume, Medicago truncatula. Results: This is the first report of integration of independent linkage maps in white clover, and adds to the literature on methyl filtered GeneThresher(R)-derived microsatellite (simple sequence repeat; SSR) markers for linkage mapping. Gene-targeted SSR markers were discovered in a GeneThresher(R) (TrGT) methyl-filtered database of 364,539 sequences, which yielded 15,647 SSR arrays. Primers were designed for 4,038 arrays and of these, 465 TrGT-SSR markers were used for parental consensus genetic linkage analysis in an F1 mapping population (MP2). This was merged with an EST-SSR consensus genetic map of an independent population (MP1), using markers to match homoeologues and develop a multi-population integrated map of the white clover genome. This integrated map (IM) includes 1109 loci based on 804 SSRs over 1274 cM, covering 97% of the genome at a moderate density of one locus per 1.2 cM. Eighteen candidate genes and one morphological marker were also placed on the IM. Despite being derived from disparate populations and marker sources, the component maps and the derived IM had consistent representations of the white clover genome for marker order and genetic length. In silico analysis at an E-value threshold of 1e-20 revealed substantial co-linearity with the Medicago truncatula genome, and indicates a translocation between T. repens groups 2 and 6 relative to M. truncatula. Conclusions: This integrated genetic linkage analysis provides a consistent and comprehensive linkage analysis of the white clover genome, with alignment to a model forage legume. Associated marker locus information, particularly the homoeologue-specific markers, offers a new resource for forage legume research to enable genetic analysis and improvement of this forage and grassland species.

Description: Background: There is a lack of understanding the evolutionary forces driving niche segregation of closely related organisms. In addition, pinpointing the genes driving ecological divergence is a key goal in molecular ecology. Here, larval transcriptome sequences obtained by next-generation-sequencing are used to address these issues in a morphologically cryptic sister species pair of non-biting midges (Chironomus riparius and C. piger). Results: More than eight thousand orthologous open reading frames were screened for interspecific divergence and intraspecific polymorphisms. Despite a small mean sequence divergence of 1.53% between the sister species, 25.1% of 18,115 observed amino acid substitutions were inferred by alpha statistics to be driven by positive selection. Applying McDonald-Kreitman tests to 715 alignments of gene orthologues identified eleven (1.5%) genes driven by positive selection. Conclusions: Three candidate genes were identified as potentially responsible for the observed niche segregation concerning nitrite concentration, habitat temperature and water conductivity. Additionally, signs of positive selection in the hydrogen sulfide detoxification pathway were detected, providing a new plausible hypothesis for the species' ecological differentiation. Finally, a divergently selected, nuclear encoded mitochondrial ribosomal protein may contribute to reproductive isolation due to cytonuclear coevolution.

Description: Background: Histone wrapping of DNA into nucleosomes almost certainly evolved in the Archaea, and predates Eukaryotes. In Eukaryotes, nucleosome positioning plays a central role in regulating gene expression and is directed by primary sequence motifs that together form a nucleosome positioning code. The experiments reported were undertaken to determine if archaeal histone assembly conforms to the nucleosome positioning code. Results: Eukaryotic nucleosome positioning is favored and directed by phased helical repeats of AA/TT/AT/TA and CC/GG/CG/GC dinucleotides, and disfavored by longer AT-rich oligonucleotides. Deep sequencing of genomic DNA protected from micrococcal nuclease digestion by assembly into archaeal nucleosomes has established that archaeal nucleosome assembly is also directed and positioned by these sequence motifs, both in vivo in Methanothermobacter thermautotrophicus and Thermococcus kodakarensis and in vitro in reaction mixtures containing only one purified archaeal histone and genomic DNA. Archaeal nucleosomes assembled at the same locations in vivo and in vitro, with much reduced assembly immediately upstream of open reading frames and throughout the ribosomal rDNA operons. Providing further support for a common positioning code, archaeal histones assembled into nucleosomes on eukaryotic DNA and eukaryotic histones into nucleosomes on archaeal DNA at the same locations. T. kodakarensis has two histones, designated HTkA and HTkB, and strains with either but not both histones deleted grow normally but do exhibit transcriptome differences. Comparisons of the archaeal nucleosome profiles in the intergenic regions immediately upstream of genes that exhibited increased or decreased transcription in the absence of HTkA or HTkB revealed substantial differences but no consistent pattern of changes that would correlate directly with archaeal nucleosome positioning inhibiting or stimulating transcription. Conclusions: The results obtained establish that an archaeal histone and a genome sequence together are sufficient to determine where archaeal nucleosomes preferentially assemble and where they avoid assembly. We confirm that the same nucleosome positioning code operates in Archaea as in Eukaryotes and presumably therefore evolved with the histone-fold mechanism of DNA binding and compaction early in the archaeal lineage, before the divergence of Eukaryotes.

Description: Background: Streptococcus equi ssp. zooepidemicus (S. zooepidemicus) is an important pathogen causing swine streptococcosis in China. Pathogenicity islands (PAIs) of S. zooepidemicus have been transferred among bacteria through horizontal gene transfer (HGT) and play important roles in the adaptation and increased virulence of S. zooepidemicus. The present study used comparative genomics to examine the different pathogenicities of S. zooepidemicus. Results: Genome of S. zooepidemicus ATCC35246 (Sz35246) comprises 2,167,264-bp of a single circular chromosome, with a GC content of 41.65%. Comparative genome analysis of Sz35246, S. zooepidemicus MGCS10565 (Sz10565), Streptococcus equi. ssp. equi. 4047 (Se4047) and S. zooepidemicus H70 (Sz70) identified 320 Sz35246-specific genes, clustered into three toxin-antitoxin (TA) systems PAIs and one restriction modification system (RM system) PAI. These four acquired PAIs encode proteins that may contribute to the overall pathogenic capacity and fitness of this bacterium to adapt to different hosts. Analysis of the in vivo and in vitro transcriptomes of this bacterium revealed differentially expressed PAI genes and non-PAI genes, suggesting that Sz35246 possess mechanisms for infecting animals and adapting to a wide range of host environments. Analysis of the genome identified potential Sz35246 virulence genes. Genes of the Fim III operon were presumed to be involved in breaking the host-restriction of Sz35246. Conclusion: Genome wide comparisons of Sz35246 with three other strains and transcriptome analysis revealed novel genes related to bacterial virulence and breaking the host-restriction. Four specific PAIs, which were judged to have been transferred into Sz35246 genome through HGT, were identified for the first time. Further analysis of the TA and RM systems in the PAIs will improve our understanding of the pathogenicity of this bacterium and could lead to the development of diagnostics and vaccines.

Description: Background: Streptococcus pneumoniae is a leading cause of childhood morbidity and mortality worldwide, despite the availability of effective pneumococcal vaccines. Understanding the molecular interactions between the bacterium and the host will contribute to the control and prevention of pneumococcal disease. Results: We used a combination of adherence assays, mutagenesis and functional genomics to identify novel factors involved in adherence. By contrasting these processes in two pneumococcal strains, TIGR4 and G54, we showed that adherence and invasion capacities vary markedly by strain. Electron microscopy showed more adherent bacteria in association with membranous pseudopodia in the TIGR4 strain. Operons for cell wall phosphorylcholine incorporation (lic), manganese transport (psa) and phosphate utilization (phn) were up-regulated in both strains on exposure to epithelial cells. Pneumolysin, pili, stress protection genes (adhC-czcD) and genes of the type II fatty acid synthesis pathway were highly expressed in the naturally more invasive strain, TIGR4. Deletion mutagenesis of five gene regions identified as regulated in this study revealed attenuation in adherence. Most strikingly, [increment]SP_1922 which was predicted to contain a B-cell epitope and revealed significant attenuation in adherence, appeared to be expressed as a part of an operon that includes the gene encoding the cytoplasmic pore-forming toxin and vaccine candidate, pneumolysin. Conclusion: This work identifies a list of novel potential pneumococcal adherence determinants.

Description: Background: Previous molecular genetic studies of physiology and pigmentation of sheep skin have focused primarily on a limited number of genes and proteins. To identify additional genes that may play important roles in coat color regulation, Illumina sequencing technology was used to catalog global gene expression profiles in skin of sheep with white versus black coat color. Results: There were 90,006 and 74,533 unigenes assembled from the reads obtained from white and black sheep skin, respectively. Genes encoding for the ribosomal proteins and keratin associated proteins were most highly expressed. A total of 2,235 known genes were differentially expressed in black versus white sheep skin, with 479 genes up-regulated and 1,756 genes down-regulated. A total of 845 novel genes were differentially expressed in black versus white sheep skin, consisting of 107 genes which were up-regulated (including 2 highly expressed genes exclusively expressed in black sheep skin) and 738 genes that were down-regulated. There was also a total of 49 known coat color genes expressed in sheep skin, from which 13 genes showed higher expression in black sheep skin. Many of these up-regulated genes, such as DCT, MATP, TYR and TYRP1, are members of the components of melanosomes and their precursor ontology category. Conclusion: The white and black sheep skin transcriptome profiles obtained provide a valuable resource for future research to understand the network of gene expression controlling skin physiology and melanogenesis in sheep.

Description: Background: Non-coding small RNAs, ranging from 20 to 30 nucleotides in length, mediate the regulation of gene expression and play important roles in many biological processes. One class of small RNAs, microRNAs (miRNAs), are highly conserved across taxa and mediate the regulation of the chromatin state and the post-transcriptional regulation of messenger RNA (mRNA). Another class of small RNAs is the Piwi-interacting RNAs, which play important roles in the silencing of transposons and other functional genes. Although the biological functions of the different small RNAs have been elucidated in several laboratory animals, little is known regarding naturally occurring variation in small RNA transcriptomes among closely related species. Results: We employed next-generation sequencing technology to compare the expression profiles of brain small RNAs between sympatric species of the Japanese threespine stickleback (Gasterosteus aculeatus). We identified several small RNAs that were differentially expressed between sympatric Pacific Ocean and Japan Sea sticklebacks. Potential targets of several small RNAs were identified as repetitive sequences. Female-biased miRNA expression from the old X chromosome was also observed, and it was attributed to the degeneration of the Y chromosome. Conclusions: Our results suggest that expression patterns of small RNA can differ between incipient species and may be a potential mechanism underlying differential mRNA expression and transposon activity.

Description: Background; In all Metazoa, transcription is inactive during the first mitotic cycles after fertilisation. In Drosophilamelanogaster, Zygotic Genome Activation (ZGA) occurs in two waves, starting respectively at mi-totic cycles 8 (approximately 60 genes) and 14 (over a thousand genes). The regulatory mechanismsunderlying these drastic transcriptional changes remain largely unknown.Results; We developed an original gene clustering method based on discretized transition profiles, and appliedit to datasets from three landmark early embryonic transcriptome studies. We identified 417 genes sig-nificantly up-regulated during ZGA. De novo motif discovery returned nine motifs over-represented intheir non-coding sequences (upstream, introns, UTR), three of which correspond to previously knowntranscription factors: Zelda, Tramtrack and Trithorax-like (Trl). The nine discovered motifs werecombined to scan ZGA-associated regions and predict about 1300 putative cis-regulatory modules.The fact that Trl is known to act as chromatin remodelling factor suggests that epigenetic regulationmight play an important role in zygotic genome activation. We thus systematically compared thelocations of predicted CRMs with ChIP-seq profiles for various transcription factors, 38 epigeneticmarks from ModENCODE, and DNAse1 accessibility profiles. This analysis highlighted a strong andspecific enrichment of predicted ZGA-associated CRMs for Zelda, CBP, Trl binding sites, as well asfor histone marks associated with active enhancers (H3K4me1) and for open chromatin regions.Conclusion;Based on the results of our computational analyses, we suggest a temporal model explaining the onsetof zygotic genome activation by the combined action of transcription factors and epigenetic signals.Although this study is mainly based on the analysis of publicly available transcriptome and ChiP-seqdatasets, the resulting model suggests novel mechanisms that underly the coordinated activation ofseveral hundreds genes at a precise time point during embryonic development.

Description: Background: Nicotianamine (NA), a ubiquitous molecule in plants, is an important metal ion chelator and the main precursor for phytosiderophores biosynthesis. Considerable progress has been achieved in cloning and characterizing the functions of nicotianamine synthase (NAS) in plants including barley, Arabidopsis and rice. Maize is not only an important cereal crop, but also a model plant for genetics and evolutionary study. The genome sequencing of maize was completed, and many gene families were identified. Although three NAS genes have been characterized in maize, there is still no systematic identification of maize NAS family by genomic mining. Results: In this study, nine NAS genes in maize were identified and their expression patterns in different organs including developing seeds were determined. According to the evolutionary relationship and tissue specific expression profiles of ZmNAS genes, they can be subgrouped into two classes. Moreover, the expression patterns of ZmNAS genes in response to fluctuating metal status were analysed. The class I ZmNAS genes were induced under Fe deficiency and were suppressed under Fe excessive conditions, while the expression pattern of class II genes were opposite to class I. The complementary expression patterns of class I and class II ZmNAS genes confirmed the classification of this family. Furthermore, the histochemical localization of ZmNAS1;1/1;2 and ZmNAS3 were determined using in situ hybridization. It was revealed that ZmNAS1;1/1;2, representing the class I genes, mainly expressed in cortex and stele of roots with sufficient Fe, and its expression can expanded in epidermis, as well as shoot apices under Fe deficient conditions. On the contrary, ZmNAS3, one of the class II genes, was accumulated in axillary meristems, leaf primordia and mesophyll cells. These results suggest that the two classes of ZmNAS genes may be regulated on transcriptional level when responds to various demands for iron uptake, translocation and homeostasis. Conclusion: These results provide significant insights into the molecular bases of ZmNAS in balancing iron uptake, translocation and homeostasis in response to fluctuating environmental Fe status.

Description: Background: Understanding the processes that drive the evolution of snake venom is a topic of great research interest in molecular and evolutionary toxinology. Recent studies suggest that ontogenetic changes in venom composition are genetically controlled rather than environmentally induced. However, the molecular mechanisms underlying these changes remain elusive. Here we have explored the basis and level of regulation of the ontogenetic shift in the venom composition of the Central American rattlesnake, Crotalus s. simus using a combined proteomics and transcriptomics approach. Results: Proteomic analysis showed that the ontogenetic shift in the venom composition of C. s. simus is essentially characterized by a gradual reduction in the expression of serine proteinases and PLA2 molecules, particularly crotoxin, a beta-neurotoxic heterodimeric PLA2, concominantly with an increment of PI and PIII metalloproteinases at age 9--18 months. Comparison of the transcriptional activity of the venom glands of neonate and adult C. s. simus specimens indicated that their transcriptomes exhibit indistinguisable toxin family profiles, suggesting that the elusive mechanism by which shared transcriptomes generate divergent venom phenotypes may operate post-transcriptionally. Specifically, miRNAs with frequency count of 1000 or greater exhibited an uneven distribution between the newborn and adult datasets. Of note, 590 copies of a miRNA targeting crotoxin B-subunit was exclusively found in the transcriptome of the adult snake, whereas 1185 copies of a miRNA complementary to a PIII-SVMP mRNA was uniquely present in the newborn dataset. These results support the view that age-dependent changes in the concentration of miRNA modulating the transition from a crotoxin-rich to a SVMP-rich venom from birth through adulhood can potentially explain what is observed in the proteomic analysis of the ontogenetic changes in the venom composition of C. s. simus. Conclusions: Existing snake venom toxins are the result of early recruitment events in the Toxicofera clade of reptiles by which ordinary genes were duplicated, and the new genes selectively expressed in the venom gland and amplified to multigene families with extensive neofunctionalization throughout the approximately 112--125 million years of ophidian evolution. Our findings support the view that understanding the phenotypic diversity of snake venoms requires a deep knowledge of the mechanisms regulating the transcriptional and translational activity of the venom gland. Our results suggest a functional role for miRNAs. The impact of specific miRNAs in the modulation of venom composition, and the integration of the mechanisms responsible for the generation of these miRNAs in the evolutionary landscape of the snake's venom gland, are further challenges for future research.

Description: Background: In temperate regions, the time lag between vegetative bud burst and bud set determines the duration of the growing season of trees (i.e. the duration of wood biomass production). Dormancy, the period during which the plant is not growing, allows trees to avoid cold injury resulting from exposure to low temperatures. An understanding of the molecular machinery controlling the shift between these two phenological states is of key importance in the context of climatic change. The objective of this study was to identify genes upregulated during endo- and ecodormancy, the two main stages of bud dormancy. Sessile oak is a widely distributed European white oak species. A forcing test on young trees was first carried out to identify the period most likely to correspond to these two stages. Total RNA was then extracted from apical buds displaying endo- and ecodormancy. This RNA was used for the generation of cDNA libraries, and in-depth transcriptome characterization was performed with 454 FLX pyrosequencing technology. Results: Pyrosequencing produced a total of 495,915 reads. The data were cleaned, duplicated reads removed, and sequences were mapped onto the oak UniGene data. Digital gene expression analysis was performed, with both R statistics and the R-Bioconductor packages (edgeR and DESeq), on 6,471 contigs with read numbers 〉= 5 within any contigs. The number of sequences displaying significant differences in expression level (read abundance) between endo- and ecodormancy conditions ranged from 75 to 161, depending on the algorithm used. 13 genes displaying significant differences between conditions were selected for further analysis, and 11 of these genes, including those for glutathione-S-transferase (GST) and dehydrin xero2 (XERO2) were validated by quantitative PCR. Conclusions: The identification and functional annotation of differentially expressed genes involved in the "response to abscisic acid", "response to cold stress" and "response to oxidative stress" categories constitutes a major step towards characterization of the molecular network underlying vegetative bud dormancy, an important life history trait of long-lived organisms.

Description: Background: The recA/RAD51 gene family encodes a diverse set of recombinase proteins that affect homologous recombination, DNA-repair, and genome stability. The recA gene family is expressed across all three domains of life - Eubacteria, Archaea, and Eukaryotes - and even in some viruses. To date, efforts to resolve the deep evolutionary origins of this ancient protein family have been hindered by the high sequence divergence between paralogous groups (i.e. ~30% average pairwise identity). Results: Through large taxon sampling and the use of a phylogenetic algorithm designed for inferring evolutionary events in highly divergent paralogs, we obtained a robust, parsimonious and more refined phylogenetic history of the recA/RAD51 superfamily. Conclusions: In summary, our model for the evolution of recA/RAD51 family provides a better understanding of the ancient origin of recA proteins and the multiple events that lead to the diversification of recA homologs in eukaryotes, including the discovery of additional RAD51 sub-families.

Description: Background: The Estrogen Receptor alpha (ERalpha) is the key transcriptional regulator in luminal breast cancer and is therefore the main target for adjuvant treatment of this subtype. Luminal gene signatures are dictated by the transcriptional capacities of ERalpha, which are a direct consequence of the receptors binding preference at specific sites on the chromatin. The identification of ERalpha binding signatures on a genome-wide level has greatly enhanced our understanding of Estrogen Receptor biology in cell lines and tumours, but the technique has its limitations with respect to its applicability in limited amounts of tumour tissue. Results: Here, we present a refinement of the ChIP-seq procedures to enable transcription factor mapping on limited amounts of tissue culture cells as well as from a limited amount of tumor tissue derived from core needle biopsies. Our approach uses a carrier that can be removed prior to DNA amplification and sequencing. Conclusion: We illustrate the applicability of this refined technology by mapping the ERalpha genome-wide chromatin binding landscape in core needle biopsy material from primary breast tumours. With this, our refined technology permits for a high-resolution transcription factor mapping even from clinical samples.

Description: Background: Mycobacterium bovis, the causative agent of bovine tuberculosis, is an intracellular pathogen that can persist inside host macrophages during infection via a diverse range of mechanisms that subvert the host immune response. In the current study, we have analysed and compared the transcriptomes of M. bovis-infected monocyte-derived macrophages (MDM) purified from six Holstein-Friesian females with the transcriptomes of non-infected control MDM from the same animals over a 24 h period using strand-specific RNA sequencing (RNA-seq). In addition, we compare gene expression profiles generated using RNA-seq with those previously generated by us using the high-density Affymetrix(R) GeneChip(R) Bovine Genome Array platform from the same MDM-extracted RNA. Results: A mean of 7.2 million reads from each MDM sample mapped uniquely and unambiguously to single Bos taurus reference genome locations. Analysis of these mapped reads showed 2,584 genes (1,392 upregulated; 1,192 downregulated) and 757 putative natural antisense transcripts (558 upregulated; 119 downregulated) that were differentially expressed based on sense and antisense strand data, respectively (adjusted P-value

Description: Background: Agarwood is an expensive resinous heartwood derived from Aquilaria plants that is widely used in traditional medicines, incense and perfume. Only wounded trees can produce agarwood, and the huge demand for the agarwood products has led all Aquilaria spp. being endangered and listed in the Appendix II of the CITES (http://www.cites.org). The major components of agarwood are sesquiterpenes and phenylethyl chromones. Owing to a lack of genomic information, the molecular basis of wound-induced sesquiterpenes biosynthesis and agarwood formation remains unknown. Results: To identify the primary genes that maybe related to agarwood formation, we sequenced 2 cDNA libraries generated from healthy and wounded A. sinensis (Lour.) Gilg. A total of 89,137 unigenes with an average length of 678.65 bp were obtained, and they were annotated in detail at bioinformatics levels. Of those associated with agarwood formation, 30 putatively encoded enzymes in the sesquiterpene biosynthesis pathway, and a handful of transcription factors and protein kinases were related to wound signal transduction. Three full-length cDNAs of sesquiterpene synthases (ASS1-3) were cloned and expressed in Escherichia coli, and enzyme assays revealed that they are active enzymes, with the major products being delta-guaiene. A methyl jasmonate (MJ) induction experiment revealed that the expression of ASS was significantly induced by MJ, and the production of sesquiterpenes was elevated accordingly. The expression of some transcription factors and protein kinases, especially MYB4, WRKY4, MPKK2 and MAPK2, was also induced by MJ and coordinated with ASS expression, suggesting they maybe positive regulators of ASS. Conclusions: This study provides extensive transcriptome information for Aquilaria spp. and valuable clues for elucidating the mechanism of wound-induced agarwood sesquiterpenes biosynthesis and their regulation.

Description: Background: In recent years, genome-wide association studies have successfully uncovered single-nucleotide polymorphisms (SNPs) associated with complex traits such as diseases and quantitative phenotypes. These variations account for a small proportion of heritability. With the development of high throughput techniques, abundant submicroscopic structural variations have been found in organisms, of which the main variations are copy number variations (CNVs). Therefore, CNVs are increasingly recognized as an important and abundant source of genetic variation and phenotypic diversity. Results: Analyses of CNVs in the genomes of three sheep breeds were performed using the Ovine SNP50 BeadChip array. A total of 238 CNV regions (CNVRs) were identified, including 219 losses, 13 gains, and six with both events (losses and gains), which cover 60.35 Mb of the sheep genomic sequence and correspond to 2.27% of the autosomal genome sequence. The length of the CNVRs on autosomes range from 13.66 kb to 1.30 Mb with a mean size of 253.57 kb, and 75 CNVRs events had a frequency 〉 3%. Among these CNVRs, 47 CNVRs identified by the PennCNV overlapped with the CNVpartition. Functional analysis indicated that most genes in the CNVRs were significantly enriched for involvement in the environmental response. Furthermore, 10 CNVRs were selected for validation and 6 CNVRs were further experimentally confirmed by qPCR. In addition, there were 57 CNVRs overlapped in our new dataset and other published ruminant CNV studies. Conclusions: In this study, we firstly constructed a sheep CNV map based on the Ovine SNP50 array. Our results demonstrated the differences of two detection tools and integration of multiple algorithms can enhance the detection of sheep genomic structure variations. Furthermore, our findings would be of help for understanding the sheep genome and provide preliminary foundation for carrying out the CNVs association studies with economically important phenotypes of sheep in the future.

Description: Background: Transcription factors (TF) play a crucial role in regulating gene expression and are fit to regulate diverse cellular processes by interacting with other proteins. A TF named calmodulin binding transcription activator (CAMTA) was identified in Arabidopsis thaliana (AtCAMTA1-6). To explore the role of CAMTA1 in drought response, the phenotypic differences and gene expression was studied between camta1 and Col-0 under drought condition. Results: In camta1, root development was abolished showing high-susceptibility to induced osmotic stress resulting in small wrinkled rosette leaves and stunted primary root. In camta1 under drought condition, we identified growth retardation, poor WUE, low photosystem II efficiency, decline in RWC and higher sensitivity to drought with reduced survivability. The microarray analysis of drought treated camta1 revealed that CAMTA1 regulates "drought recovery" as most indicative pathway along with other stress response, osmotic balance, apoptosis, DNA methylation and photosynthesis. Interestingly, majority of positively regulated genes were related to plasma membrane and chloroplast. Further, our analysis indicates that CAMTA1 regulates several stress responsive genes including RD26, ERD7, RAB18, LTPs, COR78, CBF1, HSPs etc. and promoter of these genes were enriched with CAMTA recognition cis-element. CAMTA1 probably regulate drought recovery by regulating expression of AP2-EREBP transcription factors and Abscisic acid response. Conclusion: CAMTA1 rapidly changes broad spectrum of responsive genes of membrane integrity and photosynthetic machinery by generating ABA response for challenging drought stress. Our results demonstrate the important role of CAMTA1 in regulating drought response in Arabidopsis, thus could be genetically engineered for improving drought tolerance in crop.

Description: Background: Due to the constantly growing number of sequenced microbial genomes, comparative genomics has been playing a major role in the investigation of regulatory interactions in bacteria. Regulon inference mostly remains a field of semi-manual examination since absence of a knowledgebase and informatics platform for automated and systematic investigation restricts opportunities for computational prediction. Additionally, confirming computationally inferred regulons by experimental data is critically important.Description: RegTransBase is an open-access platform with a user-friendly web interface that is publicly available at http://regtransbase.lbl.gov. It consists of two databases -- a manually collected hierarchical regulatory interactions database based on more than 7000 scientific papers which can serve as a knowledgebase for verification of predictions, and a large set of curated by experts transcription factor binding sites used in regulon inference by a variety of tools. RegTransBase captures the knowledge from published scientific literature using controlled vocabularies and contains various types of experimental data, such as: the activation or repression of transcription by an identified direct regulator; determination of the transcriptional regulatory function of a protein (or RNA) directly binding to DNA or RNA; mapping of binding sites for a regulatory protein; characterization of regulatory mutations. Analysis of the data collected from literature resulted in the creation of Putative Regulons from Experimental Data that are also available in RegTransBase. Conclusions: RegTransBase is a powerful user-friendly platform for the investigation of regulation in prokaryotes. It uses a collection of validated regulatory sequences that can be easily extracted and used to infer regulatory interactions by comparative genomics techniques thus assisting researchers in the interpretation of transcriptional regulation data.

Description: Background: Vascular system formation has been studied through molecular and genetic approaches in Arabidopsis, a herbaceous dicot that is used as a model system. Different vascular systems have developed in other plants such as crops and trees. Uncovering shared mechanisms underlying vascular development by transcriptome analysis of different vascular systems may help to transfer knowledge acquired from Arabidopsis to other economically important species. Results: Conserved vascular genes and biological processes fundamental to vascular development were explored across various plants. Through comparative transcriptome analysis, 226 genes from Arabidopsis, 217 genes from poplar and 281 genes from rice were identified as constituting 107 conserved vascular gene groups. These gene groups are expressed mainly in vascular tissues and form a complex coexpression network with multiple functional connections. To date, only half of the groups have been experimentally investigated. The conserved vascular gene groups were classified into 9 essential processes for vascular development. 18 groups (17%) lack of annotations were classified as having unknown functions. Conclusion: The study provides a map of fundamental biological processes conserved across different vascular systems. It identifies gaps in the experimental investigation of pathways active in vascular formation, which if explored, could lead to a more complete understanding of vascular development.

Description: Background: Although Candida albicans and Candida dubliniensis are most closely related, both species behave significantly different with respect to morphogenesis and virulence. In order to gain further insight into the divergent routes for morphogenetic adaptation in both species, we investigated qualitative along with quantitative differences in the transcriptomes of both organisms by cDNA deep sequencing. Results: Following genome-associated assembly of sequence reads we were able to generate experimentally verified databases containing 6016 and 5972 genes for C. albicans and C. dubliniensis, respectively. About 95% of the transcriptionally active regions (TARs) contain open reading frames while the remaining TARs most likely represent non-coding RNAs. Comparison of our annotations with publically available gene models for C. albicans and C. dubliniensis confirmed approximately 95% of already predicted genes, but also revealed so far unknown novel TARs in both species. Qualitative cross-species analysis of these databases revealed in addition to 5802 orthologs also 399 and 49 species-specific protein coding genes for C. albicans and C. dubliniensis, respectively. Furthermore, quantitative transcriptional profiling using RNA-Seq revealed significant differences in the expression of orthologs across both species. We defined a core subset of 84 hyphal-specific genes required for both species, as well as a set of 42 genes that seem to be specifically induced during hyphal morphogenesis in C. albicans. Conclusions: Species-specific adaptation in C. albicans and C. dubliniensis is governed by individual genetic repertoires but also by altered regulation of conserved orthologs on the transcriptional level.

Description: Background: Mammalian genes are regulated through the action of multiple regulatory elements, often distributed across large regions. The mechanisms that control the integration of these diverse inputs into specific gene expression patterns are still poorly understood. New approaches enabling the dissection of these mechanisms in vivo are needed. Results: Here, we describe TRACER (http://tracerdatabase.embl.de), a resource that centralizes information from a large on-going functional exploration of the mouse genome with different transposon-associated regulatory sensors. Hundreds of insertions have been mapped to specific genomic positions, and their corresponding regulatory potential has been documented by analysis of the expression of the reporter sensor gene in mouse embryos. The data can be easily accessed and provides information on the regulatory activities present in a large number of genomic regions, notably in gene-poor intervals that have been associated with human diseases. Conclusions: TRACER data enables comparisons with the expression pattern of neighbouring genes, activity of surrounding regulatory elements or with other genomic features, revealing the underlying regulatory architecture of these loci. TRACER mouse lines can also be requested for in vivo transposition and chromosomal engineering, to analyse further regions of interest.

Description: Background: Olive cDNA libraries to isolate candidate genes that can help enlightening the molecular mechanism of periodicity and / or fruit production were constructed and analyzed. For this purpose, cDNA libraries from the leaves of trees in "on year" and in "off year" in July (when fruits start to appear) and in November (harvest time) were constructed. Randomly selected 100 positive clones from each library were analyzed with respect to sequence and size. A fruit-flesh cDNA library was also constructed and characterized to confirm the reliability of each library's temporal and spatial properties. Results: Quantitative real-time RT-PCR (qRT-PCR) analyses of the cDNA libraries confirmed cDNA molecules that are associated with different developmental stages (e. g. "on year" leaves in July, "off year" leaves in July, leaves in November) and fruits. Hence, a number of candidate cDNAs associated with "on year" and "off year" were isolated. Comparison of the detected cDNAs to the current EST database of GenBank along with other non - redundant databases of NCBI revealed homologs of previously described genes along with several unknown cDNAs. Of around 500 screened cDNAs, 48 cDNA elements were obtained after eliminating ribosomal RNA sequences. These independent transcripts were analyzed using BLAST searches (cutoff E-value of 1.0E-5) against the KEGG and GenBank nucleotide databases and 37 putative transcripts corresponding to known gene functions were annotated with gene names and Gene Ontology (GO) terms. Transcripts in the biological process were found to be related with metabolic process (27%), cellular process (23%), response to stimulus (17%), localization process (8.5%), multicellular organismal process (6.25%), developmental process (6.25%) and reproduction (4.2 %). Conclusions: A putative P450 monooxigenase expressed fivefold more in the "on year" than that of "off year" leaves in July. Two putative dehydrins expressed significantly more in "on year" leaves than that of "off year" leaves in November. Homologs of UDP -- glucose epimerase, acyl - CoA binding protein, triose phosphate isomerase and a putative nuclear core anchor protein were significant in fruits only, while a homolog of an embryo binding protein / small GTPase regulator was detected in "on year" leaves only. One of the two unknown cDNAs was specific to leaves in July while the other was detected in all of the libraries except fruits. KEGG pathway analyses for the obtained sequences correlated with essential metabolisms such as galactose metabolism, amino sugar and nucleotide sugar metabolisms and photosynthesis. Detailed analysis of the results presents candidate cDNAs that can be used to dissect further the genetic basis of fruit production and / or alternate bearing which causes significant economical loss for olive growers.

Description: Background: The gene doublesex (dsx) is known as a key factor regulating genetic sex determination in many organisms. We previously identified two dsx genes (DapmaDsx1 and DapmaDsx2) from a freshwater branchiopod crustacean, Daphnia magna, which are expressed in males but not in females. D. magna produces males by parthenogenesis in response to environmental cues (environmental sex determination) and we showed that DapmaDsx1 expression during embryonic stages is responsible for the male trait development. The D. magna dsx genes are thought to have arisen by a cladoceran-specific duplication; therefore, to investigate evolutionary conservation of sex specific expression of dsx genes and to further assess their functions in the environmental sex determination, we searched for dsx homologs in four closely related cladoceran species. Results: We identified homologs of both dsx genes from, D. pulex, D. galeata, and Ceriodaphnia dubia, yet only a single dsx gene was found from Moina macrocopa. The deduced amino acid sequences of all 9 dsx homologs contained the DM and oligomerization domains, which are characteristic for all arthropod DSX family members. Molecular phylogenetic analysis suggested that the dsx gene duplication likely occurred prior to the divergence of these cladoceran species, because that of the giant tiger prawn Penaeus monodon is rooted ancestrally to both DSX1 and DSX2 of cladocerans. Therefore, this result also suggested that M. macrocopa lost dsx2 gene secondarily. Furthermore, all dsx genes identified in this study showed male-biased expression levels, yet only half of the putative 5' upstream regulatory elements are preserved in D. magna and D. pulex. Conclusions: The all dsx genes of five cladoceran species examined had similar amino acid structure containing highly conserved DM and oligomerization domains, and exhibited sexually dimorphic expression patterns, suggesting that these genes may have similar functions for environmental sex determination in cladocerans.

Description: Background: Nutritional symbioses play a central role in insects' adaptation to specialized diets and in their evolutionary success. The obligatory symbiosis between the pea aphid, Acyrthosiphon pisum, and the bacterium, Buchnera aphidicola, is no exception as it enables this important agricultural pest insect to develop on a diet exclusively based on plant phloem sap. The symbiotic bacteria provide the host with essential amino acids lacking in its diet but necessary for the rapid embryonic growth seen in the parthenogenetic viviparous reproduction of aphids. The aphid furnishes, in exchange, non-essential amino acids and other important metabolites. Understanding the regulations acting on this integrated metabolic system during the development of this insect is essential in elucidating aphid biology. Results: We used a microarray-based approach to analyse gene expression in the late embryonic and the early larval stages of the pea aphid, characterizing, for the first time, the transcriptional profiles in these developmental phases. Our analyses allowed us to identify key genes in the phenylalanine, tyrosine and dopamine pathways and we identified ACYPI004243, one of the four genes encoding for the aspartate transaminase (E.C. 2.6.1.1), as specifically regulated during development. Indeed, the tyrosine biosynthetic pathway is crucial for the symbiotic metabolism as it is shared between the two partners, all the precursors being produced by B. aphidicola. Our microarray data are supported by HPLC amino acid analyses demonstrating an accumulation of tyrosine at the same developmental stages, with an up-regulation of the tyrosine biosynthetic genes. Tyrosine is also essential for the synthesis of cuticular proteins and it is an important precursor for cuticle maturation: together with the up-regulation of tyrosine biosynthesis, we observed an up-regulation of cuticular genes expression. We were also able to identify some amino acid transporter genes which are essential for the switch over to the late embryonic stages in pea aphid development. Conclusions: Our data show that, in the development of A. pisum, a specific host gene set regulates the biosynthetic pathways of amino acids, demonstrating how the regulation of gene expression enables an insect to control the production of metabolites crucial for its own development and symbiotic metabolism.

Description: Background: Efficient screening of bacterial artificial chromosome (BAC) libraries with polymerase chain reaction (PCR)-based markers are feasible provided that a multidimensional pooling strategy is implemented. Single nucleotide polymorphisms (SNPs) can be screened in multiplexed format, therefore this marker type lends itself particularly well for medium- to high-throughput applications. Combining the power of multiplex-PCR assays with a multidimensional pooling system may prove to be especially challenging in a polyploid genome. In polyploid genomes two classes of SNPs need to be distinguished, polymorphisms between accessions (intragenomic SNPs) and those differentiating between homoeologous genomes (intergenomic SNPs). We have assessed whether the highly parallel Illumina GoldenGate(R) Genotyping Assay is suitable for the screening of a BAC library of the polyploid Brassica napus genome. Results: A multidimensional screening platform was developed for a Brassica napus BAC library which is composed of almost 83,000 clones. Intragenomic and intergenomic SNPs were included in Illumina's GoldenGate(R) Genotyping Assay and both SNP classes were used successfully for screening of the multidimensional BAC pools of the Brassica napus library. An optimized scoring method is proposed which is especially valuable for SNP calling of intergenomic SNPs. Validation of the genotyping results by independent methods revealed a success of approximately 80% for the multiplex PCR-based screening regardless of whether intra- or intergenomic SNPs were evaluated. Conclusions: Illumina's GoldenGate(R) Genotyping Assay can be efficiently used for screening of multidimensional Brassica napus BAC pools. SNP calling was specifically tailored for the evaluation of BAC pool screening data. The developed scoring method can be implemented independently of plant reference samples. It is demonstrated that intergenomic SNPs represent a powerful tool for BAC library screening of a polyploid genome.

Description: Background: The soilborne fungus, Verticillium dahliae, causes Verticillium wilt disease in plants. Verticillium wilt is difficult to control since V. dahliae is capable of persisting in the soil for 10 to 15 years as melanized microsclerotia, rendering crop rotation strategies for disease control ineffective. Microsclerotia of V. dahliae overwinter and germinate to produce infectious hyphae that give rise to primary infections. Consequently, microsclerotia formation, maintenance, and germination are critically important processes in the disease cycle of V. dahliae. Results: To shed additional light on the molecular processes that contribute to microsclerotia biogenesis and melanin synthesis in V. dahliae, three replicate RNA-seq libraries were prepared from 10 day-old microsclerotia (MS)-producing cultures of V. dahliae, strain VdLs.17 (average = 52.23 million reads), and those not producing microsclerotia (NoMS, average = 50.58 million reads). Analyses of these libraries for differential gene expression revealed over 200 differentially expressed genes, including up-regulation of melanogenesis-associated genes tetrahydroxynaphthalene reductase (344-fold increase) and scytalone dehydratase (231-fold increase), and additional genes located in a 48.8 kilobase melanin biosynthetic gene cluster of strain VdLs.17. Nearly 50% of the genes identified as differentially expressed in the MS library encode hypothetical proteins. Additional comparative analyses of gene expression in V. dahliae, under growth conditions that promote or preclude microsclerotial development, were conducted using a microarray approach with RNA derived from V. dahliae strain Dvd-T5, and from the amicrosclerotial vdh1 strain. Differential expression of selected genes observed by RNA-seq or microarray analysis was confirmed using RT-qPCR or Northern hybridizations. Conclusion: Collectively, the data acquired from these investigations provide additional insight into gene expression and molecular processes that occur during MS biogenesis and maturation in V. dahliae. The identified gene products could therefore potentially represent new targets for disease control through prevention of survival structure development.

Description: Background: The large-scale identification of physical protein-protein interactions (PPIs) is an important step toward understanding how biological networks evolve and generate emergent phenotypes. However, experimental identification of PPIs is a laborious and error-prone process, and current methods of PPI prediction tend to be highly conservative or require large amounts of functional data that may not be available for newly-sequenced organisms. Results: In this study we demonstrate a random-forest based technique, ENTS, for the computational prediction of protein-protein interactions based only on primary sequence data. Our approach is able to efficiently predict interactions on a whole-genome scale for any eukaryotic organism, using pairwise combinations of conserved domains and predicted subcellular localization of proteins as input features. We present the first predicted interactome for the forest tree Populus trichocarpa in addition to the predicted interactomes for Saccharomyces cerevisiae, Homo sapiens, Mus musculus, and Arabidopsis thaliana. Comparing our approach to other PPI predictors, we find that ENTS performs comparably to or better than a number of existing approaches, including several that utilize a variety of functional information for their predictions. We also find that the predicted interactions are biologically meaningful, as indicated by similarity in functional annotations and enrichment of co-expressed genes in public microarray datasets. Furthermore, we demonstrate some of the biological insights that can be gained from these predicted interaction networks. We show that the predicted interactions yield informative groupings of P. trichocarpa metabolic pathways, literature-supported associations among human disease states, and theory-supported insight into the evolutionary dynamics of duplicated genes in paleopolyploid plants. Conclusion: We conclude that the ENTS classifier will be a valuable tool for the de novo annotation of genome sequences, providing initial clues about regulatory and metabolic network topology, and revealing relationships that are not immediately obvious from traditional homology-based annotations.

Description: Background: Understanding how species-specific microRNAs (miRNAs) contribute to species-specific phenotypes is a central topic in biology. This study aimed to elucidate the role of ruminant-specific miRNAs in shaping mRNA expression divergence between ruminant and non-ruminant species. Results: We analyzed miRNA and mRNA transcriptomes generated by Illumina sequencing from whole blood samples of cattle and a closely related non-ruminant species, pig. We found evidence of expansion of cattle-specific miRNAs by analyzing miRNA conservation among 57 vertebrate species. The emergence of cattle-specific miRNAs was accompanied by accelerated sequence evolution at their target sites. Further, the target genes of cattle-specific miRNAs show markedly reduced expression compared to their pig and human orthologues. We found that target genes with conserved or non-conserved target sites of cattle-specific miRNAs exhibit reduced expression. One of the significantly enriched KEGG pathway terms for the target genes of the cattle-specific miRNAs is the insulin signalling pathway, raising the possibility that some of these miRNAs may modulate insulin resistance in ruminants. Conclusions: We provide evidence of rapid miRNA-mediated regulatory evolution in the ruminant lineage. Cattle-specific miRNAs play an important role in shaping gene expression divergence between ruminant and non-ruminant species, by influencing the expression of targets genes through both conserved and cattle-specific target sites.

Description: Background: Actinomycetes are a diverse group of medically, industrially and ecologically important bacteria, studied as much for the diseases they cause as for the cures they hold. The genomes of actinomycetes revealed that these bacteria have a large number of natural product gene clusters, although many of these are difficult to tie to products in the laboratory. Large scale comparisons of these clusters are difficult to perform due to the presence of highly similar repeated domains in the most common biosynthetic machinery: polyketide synthases (PKSs) and nonribosomal peptide synthetases (NRPSs). Results: We have used comparative genomics to provide an overview of the genomic features of a set of 102 closed genomes from this important group of bacteria with a focus on natural product biosynthetic genes. We have focused on well-represented genera and determine the occurrence of gene cluster families therein. Conservation of natural product gene clusters within Mycobacterium, Streptomyces and Frankia suggest crucial roles for natural products in the biology of each genus. The abundance of natural product classes is also found to vary greatly between genera, revealing underlying patterns that are not yet understood. Conclusions: A large-scale analysis of natural product gene clusters presents a useful foundation for hypothesis formulation that is currently underutilized in the field. Such studies will be increasingly necessary to study the diversity and ecology of natural products as the number of genome sequences available continues to grow.

Description: Background: The basidomycete Wallemia ichthyophaga from the phylogenetically distinct class Wallemiomycetes is the most halophilic fungus known to date. It requires at least 10% NaCl and thrives in saturated salt solution. To investigate the genomic basis of this exceptional phenotype, we obtained a de-novo genome sequence of the species type-strain and analysed its transcriptomic response to conditions close to the limits of its lower and upper salinity range. Results: The unusually compact genome is 9.6 Mb large and contains 1.67% repetitive sequences. Only 4884 predicted protein coding genes cover almost three quarters of the sequence. Of 639 differentially expressed genes, two thirds are more expressed at lower salinity. Phylogenomic analysis based on the largest dataset used to date (whole proteomes) positions Wallemiomycetes as a 250-million-year-old sister group of Agaricomycotina. Contrary to the closely related species Wallemia sebi, W. ichthyophaga appears to have lost the ability for sexual reproduction. Several protein families are significantly expanded or contracted in the genome. Among these, there are the P-type ATPase cation transporters, but not the sodium/ hydrogen exchanger family. Transcription of all but three cation transporters is not salt dependent. The analysis also reveals a significant enrichment in hydrophobins, which are cell-wall proteins with multiple cellular functions. Half of these are differentially expressed, and most contain an unusually large number of acidic amino acids. This discovery is of particular interest due to the numerous applications of hydrophobines from other fungi in industry, pharmaceutics and medicine. Conclusions: W. ichthyophaga is an extremophilic specialist that shows only low levels of adaptability and genetic recombination. This is reflected in the characteristics of its genome and its transcriptomic response to salt. No unusual traits were observed in common salt-tolerance mechanisms, such as transport of inorganic ions or synthesis of compatible solutes. Instead, various data indicate a role of the cell wall of W. ichthyophaga in its response to salt. Availability of the genomic sequence is expected to facilitate further research into this unique species, and shed more light on adaptations that allow it to thrive in conditions lethal to most other eukaryotes.

Description: Background: The ParS/ParR two component regulatory system plays critical roles for multidrug resistance in Pseudomonas aeruginosa. It was demonstrated that in the presence of antimicrobials, ParR enhances bacterial survival by distinct mechanisms including activation of the mexXY efflux genes, enhancement of lipopolysaccharide modification through the arn operon, and reduction of the expression of oprD porin. Results: In this study, we report on transcriptomic analyses of P. aeruginosa PAO1 wild type and parS and parR mutants growing in a defined minimal medium. Our transcriptomic analysis provides the first estimates of transcript abundance for the 5570 coding genes in P. aeruginosa PAO1. Comparative transcriptomics of P. aeruginosa PAO1 and par mutants identified a total of 464 genes regulated by ParS and ParR. Results also showed that mutations in the parS/parR system abolished expression of the mexEF-oprN operon by down-regulating the regulatory gene mexS. In addition to the known effects on drug resistance genes, transcript abundances of the quorum sensing genes (rhlIR and pqsABCDE-phnAB) were higher in both parS and parR mutants. In accordance with these results, a significant portion of the ParS/ParR regulated genes belonged to the MexEF-OprN and quorum sensing regulons. Deletion of the par genes also led to increased phenazine production and swarming motility, consistent with the up-regulation of the phenazine and rhamnolipid biosynthetic genes, respectively. Conclusion: Our results link the ParS/ParR two component signal transduction system to MexEF-OprN and quorum sensing systems in P. aeruginosa. These results expand our understanding of the roles of the ParS/ParR system in the regulation of gene expression in P. aeruginosa, especially in the absence of antimicrobials.

Description: Background: Candida glabrata follows C. albicans as the second or third most prevalent cause of candidemia worldwide. These two pathogenic yeasts are distantly related, C. glabrata being part of the Nakaseomyces, a group more closely related to Saccharomyces cerevisiae. Although C. glabrata was thought to be the only pathogenic Nakaseomyces, two new pathogens have recently been described within this group: C. nivariensis and C. bracarensis. To gain insight into the genomic changes underlying the emergence of virulence, we sequenced the genomes of these two, and three other non-pathogenic Nakaseomyces, and compared them to other sequenced yeasts. Results: Our results indicate that the two new pathogens are more closely related to the non-pathogenic N. delphensis than to C. glabrata. We uncover duplications and accelerated evolution that specifically affected genes in the lineage preceding the group containing N. delphensis and the three pathogens, which may provide clues to the higher propensity of this group to infect humans. Finally, the number of Epa-like adhesins is specifically enriched in the pathogens, particularly in C. glabrata. Conclusions: Remarkably, some features thought to be the result of adaptation of C. glabrata to a pathogenic lifestyle, are present throughout the Nakaseomyces, indicating these are rather ancient adaptations to other environments. Phylogeny suggests that human pathogenesis evolved several times, independently within the clade. The expansion of the EPA gene family in pathogens establishes an evolutionary link between adhesion and virulence phenotypes. Our analyses thus shed light onto the relationships between virulence and the recent genomic changes that occurred within the Nakaseomyces.Sequence Accession Numbers: Nakaseomyces delphensis: CAPT01000001 to CAPT01000179Candida bracarensis: CAPU01000001 to CAPU01000251Candida nivariensis: CAPV01000001 to CAPV01000123Candida castellii: CAPW01000001 to CAPW01000101Nakaseomyces bacillisporus: CAPX01000001 to CAPX01000186

Description: Background: MicroRNAs (miRNAs) and other types of small regulatory RNAs play critical roles in the regulation of gene expression at the post-transcriptional level in plants. Cotton is one of the most economically important crops, but little is known about the roles of miRNAs during cotton fiber elongation. Results: Here, we combined high-throughput sequencing with computational analysis to identify small RNAs (sRNAs) related to cotton fiber elongation in Gossypium hirsutum L. (G. hirsutum). The sequence analysis confirmed the expression of 79 known miRNA families in elongating fiber cells and identified 257 novel miRNAs, primarily derived from corresponding specific loci in the Gossypium raimondii Ulbr. (G. raimondii) genome. Furthermore, a comparison of the miRNAomes revealed that 46 miRNA families were differentially expressed throughout the elongation period. Importantly, the predicted and experimentally validated targets of eight miRNAs were associated with fiber elongation, with obvious functional relationships with calcium and auxin signal transduction, fatty acid metabolism, anthocyanin synthesis and the xylem tissue differentiation. Moreover, one tasiRNA was also identified, and its target, ARF4, was experimentally validated in vivo. Conclusion: This study not only facilitated the discovery of 257 novel low-abundance miRNAs in elongating cotton fiber cells but also revealed a potential regulatory network of nine sRNAs important for fiber elongation. The identification and characterization of miRNAs in elongating cotton fiber cells might promote the further study of fiber miRNA regulation mechanisms and provide insight into the importance of miRNAs in cotton.

Description: Background: Sunflower belongs to the largest plant family on earth, the genomically poorly explored Compositae. Downy mildew Plasmopara halstedii (Farlow) Berlese & de Toni is one of the major diseases of cultivated sunflower (Helianthus annuus L.). In the search for new sources of downy mildew resistance, the locus PlARG on linkage group 1 (LG1) originating from H. argophyllus is promising since it confers resistance against all known races of the pathogen. However, the mapping resolution in the PlARG region is hampered by significantly suppressed recombination and by limited availability of polymorphic markers. Here we examined a strategy developed for the enrichment of molecular markers linked to this specific genomic region. We combined bulked segregant analysis (BSA) with next-generation sequencing (NGS) and de novo assembly of the sunflower transcriptome for single nucleotide polymorphism (SNP) discovery in a sequence resource combining reads originating from two sunflower species, H. annuus and H. argophyllus. Results: A computational pipeline developed for SNP calling and pattern detection identified 219 candidate genes. For a proof of concept, 42 resistance gene-like sequences were subjected to experimental SNP validation. Using a high-resolution mapping population, 12 SNP markers were mapped to LG1. We successfully verified candidate sequences either co-segregating with or closely flanking PlARG. Conclusions: This study is the first successful example to improve bulked segregant analysis with de novo transcriptome assembly using next generation sequencing. The BSTA pipeline we developed provides a useful guide for similar studies in other non-model organisms. Our results demonstrate this method is an efficient way to enrich molecular markers and to identify candidate genes in a specific mapping interval.

Description: Background: Hydrogen peroxide (H2O2) is a reactive oxygen species (ROS), which is part of the oxidative burst encountered upon internalization of Salmonella enterica serovar Typhimurium (S. Typhimurium) by phagocytic cells. It has previously been established that, the ArcAB two-component system plays a critical role in ROS resistance, but the genes regulated by the system remained undetermined to date. We therefore investigated the ArcA regulon in aerobically growing S. Typhimurium before and after exposure to H2O2 by querying gene expression and other physiological changes in wild type and DeltaarcA strains. Results: In the DeltaarcA strain, expression of 292 genes showed direct or indirect regulation by ArcA in response to H2O2, of which 141were also regulated in aerobiosis, but in the opposite direction. Gene set enrichment analysis (GSEA) of the expression data from WT and DeltaarcA strains, revealed that, in response to H2O2 challenge in aerobically grown cells, ArcA down regulated multiple PEP-PTS and ABC transporters, while up regulating genes involved in glutathione and glycerolipid metabolism and nucleotide transport. Further biochemical analysis guided by GSEA results showed that deletion of arcA during aerobic growth lead to increased reactive oxygen species (ROS) production which was concomitant with an increased NADH/NAD+ ratio. In absence of ArcA under aerobic conditions, H2O2 exposure resulted in lower levels of glutathione reductase activity, leading to a decreased GSH (reduced glutathione)/ GSSG (oxidized glutathione) ratio. Conclusion: The ArcA regulon was defined in 2 conditions, aerobic growth and the combination of peroxide treatment and aerobic growth in S. Typhimurium. ArcA coordinates a response that involves multiple aspects of the carbon flux through central metabolism, which ultimately modulates the reducing potential of the cell.