Springer Science and Business Media LLC

Genes that belong to the same network are frequently co-expressed, but collectively, how the coordination of the whole transcriptome is perturbed during aging remains unclear. To explore this, we calculated the correlation of each gene in the transcriptome with every other, in the brain of young and older outbred deer mice (P. leucopus and P. maniculatus). In about 25 % of the genes, coordination was inversed during aging. Gene Ontology analysis in both species, for the genes that exhibited inverse transcriptomic coordination during aging pointed to alterations in the perception of smell, a known impairment occurring during aging. In P. leucopus, alterations in genes related to cholesterol metabolism were also identified. Among the genes that exhibited the most pronounced inversion in their coordination profiles during aging was THBS4, that encodes for thrombospondin-4, a protein that was recently identified as rejuvenation factor in mice. Relatively to its breadth, abolishment of coordination was more prominent in the long-living P. leucopus than in P. maniculatus but in the latter, the intensity of de-coordination was higher. There sults suggest that aging is associated with more stringent retention of expression profiles for some genes and more abrupt changes in others, while more subtle but widespread changes in gene expression appear protective. Our findings shed light in the mode of the transcriptional changes occurring in the brain during aging and suggest that strategies aiming to broader but more modest changes in gene expression may be preferrable to correct aging-associated deregulation in gene expression.

MicroRNAs play an important role in many fundamental biological and pathological processes. Defining the microRNAs profile underlying the processes by beneficial and detrimental lifestyles, including caloric restriction (CR), exercise and high-fat diet (HF), is necessary for understanding both normal physiology and the pathogenesis of metabolic disease. We used the microarray to detect microRNAs expression in livers from CR, EX and HF mice models. After predicted potential target genes of differentially expressed microRNAs with four algorithms, we applied GO and KEGG to analyze the function of predicted microRNA targets. We describe the overall microRNAs expression pattern, and identified 84 differentially expressed microRNAs changed by one or two or even all the three lifestyle modifications. The common and different enriched categories of gene function and main biochemical and signal transduction pathways were presented. We provided for the first time a comprehensive and thorough comparison of microRNAs expression profiles in liver among these lifestyle modifications. With this knowledge, our findings provide us with an overall vision of microRNAs in the molecular impact of lifestyle on health as well as useful clues for future and thorough research of the role of microRNAs.

Citrus blight is a citrus tree overall decline disease and causes serious losses in the citrus industry worldwide. Although it was described more than one hundred years ago, its causal agent remains unknown and its pathophysiology is not well determined, which hampers our understanding of the disease and design of suitable disease management. In this study, we sequenced and assembled the draft genome for Swingle citrumelo, one important citrus rootstock. The draft genome is approximately 280 Mb, which covers 74 % of the estimated Swingle citrumelo genome and the average coverage is around 15X. The draft genome of Swingle citrumelo enabled us to conduct transcriptome analysis of roots of blight and healthy Swingle citrumelo using RNA-seq. The RNA-seq was reliable as evidenced by the high consistence of RNA-seq analysis and quantitative reverse transcription PCR results (R2 = 0.966). Comparison of the gene expression profiles between blight and healthy root samples revealed the molecular mechanism underneath the characteristic blight phenotypes including decline, starch accumulation, and drought stress. The JA and ET biosynthesis and signaling pathways showed decreased transcript abundance, whereas SA-mediated defense-related genes showed increased transcript abundance in blight trees, suggesting unclassified biotrophic pathogen was involved in this disease. Overall, the Swingle citrumelo draft genome generated in this study will advance our understanding of plant biology and contribute to the citrus breeding. Transcriptome analysis of blight and healthy trees deepened our understanding of the pathophysiology of citrus blight.

Transcription promoters are fundamental genomic cis-elements controlling gene expression. They can be classified into two types by the degree of imprecision of their transcription start sites: peak promoters, which initiate transcription from a narrow genomic region; and broad promoters, which initiate transcription from a wide-ranging region. Eukaryotic transcription initiation is suggested to be associated with the genomic positions and modifications of nucleosomes. For instance, it has been recently shown that histone with H3K9 acetylation (H3K9ac) is more likely to be distributed around broad promoters rather than peak promoters; it can thus be inferred that there is an association between histone H3K9 and promoter architecture. Here, we performed a systematic analysis of transcription promoters and gene expression, as well as of epigenetic histone behaviors, including genomic position, stability within the chromatin, and several modifications. We found that, in humans, broad promoters, but not peak promoters, generally had significant associations with nucleosome positioning and modification. Specifically, around broad promoters histones were highly distributed and aligned in an orderly fashion. This feature was more evident with histones that were methylated or acetylated; moreover, the nucleosome positions around the broad promoters were more stable than those around the peak ones. More strikingly, the overall expression levels of genes associated with broad promoters (but not peak promoters) with modified histones were significantly higher than the levels of genes associated with broad promoters with unmodified histones. These results shed light on how epigenetic regulatory networks of histone modifications are associated with promoter architecture.

Bread wheat (Triticum aestivum) has a large, complex and hexaploid genome consisting of A, B and D homoeologous chromosome sets. Therefore each wheat gene potentially exists as a trio of A, B and D homoeoloci, each of which may contribute differentially to wheat phenotypes. We describe a novel approach combining wheat cytogenetic resources (chromosome substitution ‘nullisomic-tetrasomic’ lines) with next generation deep sequencing of gene transcripts (RNA-Seq), to directly and accurately identify homoeologue-specific single nucleotide variants and quantify the relative contribution of individual homoeoloci to gene expression. We discover, based on a sample comprising ~5-10% of the total wheat gene content, that at least 45% of wheat genes are expressed from all three distinct homoeoloci. Most of these genes show strikingly biased expression patterns in which expression is dominated by a single homoeolocus. The remaining ~55% of wheat genes are expressed from either one or two homoeoloci only, through a combination of extensive transcriptional silencing and homoeolocus loss. We conclude that wheat is tending towards functional diploidy, through a variety of mechanisms causing single homoeoloci to become the predominant source of gene transcripts. This discovery has profound consequences for wheat breeding and our understanding of wheat evolution.

Cork oak (Quercus suber) is one of the rare trees with the ability to produce cork, a material widely used to make wine bottle stoppers, flooring and insulation materials, among many other uses. The molecular mechanisms of cork formation are still poorly understood, in great part due to the difficulty in studying a species with a long life-cycle and for which there is scarce molecular/genomic information. Cork oak forests are of great ecological importance and represent a major economic and social resource in Southern Europe and Northern Africa. However, global warming is threatening the cork oak forests by imposing thermal, hydric and many types of novel biotic stresses. Despite the economic and social value of the Q. suber species, few genomic resources have been developed, useful for biotechnological applications and improved forest management. We generated in excess of 7 million sequence reads, by pyrosequencing 21 normalized cDNA libraries derived from multiple Q. suber tissues and organs, developmental stages and physiological conditions. We deployed a stringent sequence processing and assembly pipeline that resulted in the identification of ~159,000 unigenes. These were annotated according to their similarity to known plant genes, to known Interpro domains, GO classes and E.C. numbers. The phylogenetic extent of this ESTs set was investigated, and we found that cork oak revealed a significant new gene space that is not covered by other model species or EST sequencing projects. The raw data, as well as the full annotated assembly, are now available to the community in a dedicated web portal at http://www.corkoakdb.org . This genomic resource represents the first trancriptome study in a cork producing species. It can be explored to develop new tools and approaches to understand stress responses and developmental processes in forest trees, as well as the molecular cascades underlying cork differentiation and disease response.

Heat stress (HS) is a major stress event in the life of an animal, with detrimental upshots in production and health. Long-non-coding RNAs (lncRNAs) play an important role in many biological processes by transcriptional regulation. However, no research has been reported on the characterization and functionality of lncRNAs in heat-stressed rats. We studied expression levels of lncRNAs in rats during HS, using strand-specific RNA sequencing. Six rats, three in each of Control (22 ± 1 °C) and H120 (42 °C for 120 min) experimental groups, were used to screen for lncRNAs in their liver and adrenal glands. Totally, 4498 and 7627 putative lncRNAs were identified in liver and adrenal glands of the Control and H120 groups, respectively. The majority of lncRNAs were relatively shorter and contained fewer exons than protein-coding transcripts. In total, 482 (174 up-regulated and 308 down-regulated) and 271 (126 up-regulated and 145 down-regulated) differentially-expressed lncRNAs (DElncRNAs, P < 0.05) were identified in the liver and adrenal glands of the Control and H120 groups, respectively. Furthermore, 1274, 121, and 73 target differentially-expressed genes (DEGs) in the liver were predicted to interact with DElncRNAs based on trans−/cis- and sequence similarity regulatory modes. Functional annotation analyses indicated that these DEGs were mostly significantly enriched in insulin signalling, myeloid leukaemia, and glucagon signalling pathways. Similarly, 437, 73 and 41 target DEGs in the adrenal glands were mostly significantly enriched in the cell cycle (trans-prediction) and lysosome pathways (cis-prediction). The DElncRNAs interacting with DEGs that encode heat shock proteins (HSPs) may play an important role in HS response, which include Hsf4, Dnaja1, Dnajb4, Hsph1 and Hspb1 in the liver, and Dnajb13 and Hspb8 in the adrenal glands. The strand-specific RNA sequencing findings were also further verified through RT-qPCR. This study is the first to provide a detailed characterization and functional analysis of expression levels of lncRNAs in liver and adrenal glands of heat-stressed rats, which provides basis for further studies on the biological functions of lncRNAs under heat stress in rats and other mammalian species.

The 2010 annual conference of the Asia Pacific Bioinformatics Network (APBioNet), Asia’s oldest bioinformatics organisation formed in 1998, was organized as the 9th International Conference on Bioinformatics (InCoB), Sept. 26-28, 2010 in Tokyo, Japan. Initially, APBioNet created InCoB as forum to foster bioinformatics in the Asia Pacific region. Given the growing importance of interdisciplinary research, InCoB2010 included topics targeting scientists in the fields of genomic medicine, immunology and chemoinformatics, supporting translational research. Peer-reviewed manuscripts that were accepted for publication in this supplement, represent key areas of research interests that have emerged in our region. We also highlight some of the current challenges bioinformatics is facing in the Asia Pacific region and conclude our report with the announcement of APBioNet’s 100 BioDatabases (BioDB100) initiative. BioDB100 will comply with the database criteria set out earlier in our proposal for Minimum Information about a Bioinformatics and Investigation (MIABi), setting the standards for biocuration and bioinformatics research, on which we will report at the next InCoB, Nov. 27 – Dec. 2, 2011 at Kuala Lumpur, Malaysia.

Infants born prematurely, particularly extremely low birth weight infants (ELBW) have altered gut microbial communities. Factors such as maternal health, gut immaturity, delivery mode, and antibiotic treatments are associated with microbiota disturbances, and are linked to an increased risk of certain diseases such as necrotising enterocolitis. Therefore, there is a requirement to optimally characterise microbial profiles in this at-risk cohort, via standardisation of methods, particularly for studying the influence of microbiota therapies (e.g. probiotic supplementation) on community profiles and health outcomes. Profiling of faecal samples using the 16S rRNA gene is a cost-efficient method for large-scale clinical studies to gain insights into the gut microbiota and additionally allows characterisation of cohorts were sample quantities are compromised (e.g. ELBW infants). However, DNA extraction method, and the 16S rRNA region targeted can significantly change bacterial community profiles obtained, and so confound comparisons between studies. Thus, we sought to optimise a 16S rRNA profiling protocol to allow standardisation for studying ELBW infant faecal samples, with or without probiotic supplementation. Using ELBW faecal samples, we compared three different DNA extraction methods, and subsequently PCR amplified and sequenced three hypervariable regions of the 16S rRNA gene (V1 + V2 + V3), (V4 + V5) and (V6 + V7 + V8), and compared two bioinformatics approaches to analyse results (OTU and paired end). Paired shotgun metagenomics was used as a ‘gold-standard’. Results indicated a longer bead-beating step was required for optimal bacterial DNA extraction and that sequencing regions (V1 + V2 + V3) and (V6 + V7 + V8) provided the most representative taxonomic profiles, which was confirmed via shotgun analysis. Samples sequenced using the (V4 + V5) region were found to be underrepresented in specific taxa including Bifidobacterium, and had altered diversity profiles. Both bioinformatics 16S rRNA pipelines used in this study (OTU and paired end) presented similar taxonomic profiles at genus level. We determined that DNA extraction from ELBW faecal samples, particularly those infants receiving probiotic supplementation, should include a prolonged beat-beating step. Furthermore, use of the 16S rRNA (V1 + V2 + V3) and (V6 + V7 + V8) regions provides reliable representation of ELBW microbiota profiles, while inclusion of the (V4 + V5) region may not be appropriate for studies where Bifidobacterium constitutes a resident microbiota member.

Duplications of large genomic segments provide genetic diversity in genome evolution. Despite their importance, how these duplications are generated remains uncertain, particularly for distant duplicated genomic segments. Here we provide evidence of the participation of circular DNA intermediates in the single generation of some large human segmental duplications. A specific reversion of sequence order from A-B/C-D to B-A/D-C between duplicated segments and the presence of only microhomologies and short indels at the evolutionary breakpoints suggest a circularization of the donor ancestral locus and an accidental replicative interaction with the acceptor locus. This novel mechanism of random genomic mutation could explain several distant genomic duplications including some of the ones that took place during recent human evolution.