Meta-review of protein network regulating obesity between validated obesity candidate genes in the white adipose tissue of high-fat diet-induced obese C57BL/6J mice.

Worldwide obesity and related comorbidities are increasing, but identifying new therapeutic targets remains a challenge. A plethora of microarray studies in diet-induced obesity models has provided large datasets of obesity associated genes. In this review, we describe an approach to examine the underlying molecular network regulating obesity, and we discuss interactions between obesity candidate genes. We conducted network analysis on functional protein-protein interactions associated with 25 obesity candidate genes identified in a literature-driven approach based on published microarray studies of diet-induced obesity. The obesity candidate genes were closely associated with lipid metabolism and inflammation. Peroxisome proliferator activated receptor gamma (Pparg) appeared to be a core obesity gene, and obesity candidate genes were highly interconnected, suggesting a coordinately regulated molecular network in adipose tissue. In conclusion, the current network analysis approach may help elucidate the underlying molecular network regulating obesity and identify anti-obesity targets for therapeutic intervention.

Biological features of core networks that result from a high-fat diet in hepatic and pulmonary tissues in mammary tumour-bearing, obesity-resistant mice.

We previously demonstrated that the chronic consumption of a high-fat diet (HFD) promotes lung and liver metastases of 4T1 mammary carcinoma cells in obesity-resistant BALB/c mice. To examine early transcriptional responses to tumour progression in the liver and lungs of HFD-fed mice, 4-week-old female BALB/c mice were divided into four groups: sham-injected, control diet (CD)-fed; sham-injected, HFD-fed (SH); 4T1 cell-injected, CD-fed (TC); 4T1 cell-injected, HFD-fed (TH). Following 16 weeks of either a CD or HFD, 4T1 cells were injected into the mammary fat pads of mice in the TC and TH groups and all mice were continuously fed identical diets. At 14 d post-injection, RNA was isolated from hepatic and pulmonary tissues for microarray analysis of mRNA expression. Functional annotation and core network analyses were conducted for the TH/SH Unique gene set. Inflammation in hepatic tissues and cell mitosis in pulmonary tissues were the most significant biological functions in the TH/SH Unique gene set. The biological core networks of the hepatic TH/SH Unique gene set were characterised as those genes involved in the activation of acute inflammatory responses (Orm1, Lbp, Hp and Cfb), disordered lipid metabolism and deregulated cell cycle progression. Networks of the pulmonary Unique gene set displayed the deregulation of cell cycle progression (Cdc20, Cdk1 and Bub1b). These HFD-influenced alterations may have led to favourable conditions for the formation of both pro-inflammatory and pro-mitotic microenvironments in the target organs that promote immune cell infiltration and differentiation, as well as the infiltration and proliferation of metastatic tumour cells.

Small RNA and transcriptome deep sequencing proffers insight into floral gene regulation in Rosa cultivars.

BACKGROUND: Roses (Rosa sp.), which belong to the family Rosaceae, are the most economically important ornamental plants--making up 30% of the floriculture market. However, given high demand for roses, rose breeding programs are limited in molecular resources which can greatly enhance and speed breeding efforts. A better understanding of important genes that contribute to important floral development and desired phenotypes will lead to improved rose cultivars. For this study, we analyzed rose miRNAs and the rose flower transcriptome in order to generate a database to expound upon current knowledge regarding regulation of important floral characteristics. A rose genetic database will enable comprehensive analysis of gene expression and regulation via miRNA among different Rosa cultivars. RESULTS: We produced more than 0.5 million reads from expressed sequences, totalling more than 110 million bp. From these, we generated 35,657, 31,434, 34,725, and 39,722 flower unigenes from Rosa hybrid: 'Vital', 'Maroussia', and 'Sympathy' and Rosa rugosa Thunb., respectively. The unigenes were assigned functional annotations, domains, metabolic pathways, Gene Ontology (GO) terms, Plant Ontology (PO) terms, and MIPS Functional Catalogue (FunCat) terms. Rose flower transcripts were compared with genes from whole genome sequences of Rosaceae members (apple, strawberry, and peach) and grape. We also produced approximately 40 million small RNA reads from flower tissue for Rosa, representing 267 unique miRNA tags. Among identified miRNAs, 25 of them were novel and 242 of them were conserved miRNAs. Statistical analyses of miRNA profiles revealed both shared and species-specific miRNAs, which presumably effect flower development and phenotypes. CONCLUSIONS: In this study, we constructed a Rose miRNA and transcriptome database, and we analyzed the miRNAs and transcriptome generated from the flower tissues of four Rosa cultivars. The database provides a comprehensive genetic resource which can be used to better understand rose flower development and to identify candidate genes for important phenotypes.

Complete genome sequence of the metabolically versatile plant growth-promoting endophyte Variovorax paradoxus S110.

Variovorax paradoxus is a microorganism of special interest due to its diverse metabolic capabilities, including the biodegradation of both biogenic compounds and anthropogenic contaminants. V. paradoxus also engages in mutually beneficial interactions with both bacteria and plants. The complete genome sequence of V. paradoxus S110 is composed of 6,754,997 bp with 6,279 predicted protein-coding sequences within two circular chromosomes. Genomic analysis has revealed multiple metabolic features for autotrophic and heterotrophic lifestyles. These metabolic diversities enable independent survival, as well as a symbiotic lifestyle. Consequently, S110 appears to have evolved into a superbly adaptable microorganism that is able to survive in ever-changing environmental conditions. Based on our findings, we suggest V. paradoxus S110 as a potential candidate for agrobiotechnological applications, such as biofertilizer and biopesticide. Because it has many associations with other biota, it is also suited to serve as an additional model system for studies of microbe-plant and microbe-microbe interactions.

Genome-wide identification of the subcellular localization of the Escherichia coli B proteome using experimental and computational methods.

Escherichia coli K-12 and B strains have most widely been employed for scientific studies as well as industrial applications. Recently, the complete genome sequences of two representative descendants of E. coli B strains, REL606 and BL21(DE3), have been determined. Here, we report the subproteome reference maps of E. coli B REL606 by analyzing cytoplasmic, periplasmic, inner and outer membrane, and extracellular proteomes based on the genome information using experimental and computational approaches. Among the total of 3487 spots, 651 proteins including 410 non-redundant proteins were identified and characterized by 2-DE and LC-MS/MS; they include 440 cytoplasmic, 45 periplasmic, 50 inner membrane, 61 outer membrane, and 55 extracellular proteins. In addition, subcellular localizations of all 4205 ORFs of E. coli B were predicted by combined computational prediction methods. The subcellular localizations of 1812 (43.09%) proteins of currently unknown function were newly assigned. The results of computational prediction were also compared with the experimental results, showing that overall precision and recall were 92.16 and 92.16%, respectively. This work represents the most comprehensive analyses of the subproteomes of E. coli B, and will be useful as a reference for proteome profiling studies under various conditions. The complete proteome data are available online (http://ecolib.kaist.ac.kr).

Comparative analysis of pepper and tomato reveals euchromatin expansion of pepper genome caused by differential accumulation of Ty3/Gypsy-like elements.

BACKGROUND: Among the Solanaceae plants, the pepper genome is three times larger than that of tomato. Although the gene repertoire and gene order of both species are well conserved, the cause of the genome-size difference is not known. To determine the causes for the expansion of pepper euchromatic regions, we compared the pepper genome to that of tomato. RESULTS: For sequence-level analysis, we generated 35.6 Mb of pepper genomic sequences from euchromatin enriched 1,245 pepper BAC clones. The comparative analysis of orthologous gene-rich regions between both species revealed insertion of transposons exclusively in the pepper sequences, maintaining the gene order and content. The most common type of the transposon found was the LTR retrotransposon. Phylogenetic comparison of the LTR retrotransposons revealed that two groups of Ty3/Gypsy-like elements (Tat and Athila) were overly accumulated in the pepper genome. The FISH analysis of the pepper Tat elements showed a random distribution in heterochromatic and euchromatic regions, whereas the tomato Tat elements showed heterochromatin-preferential accumulation. CONCLUSIONS: Compared to tomato pepper euchromatin doubled its size by differential accumulation of a specific group of Ty3/Gypsy-like elements. Our results could provide an insight on the mechanism of genome evolution in the Solanaceae family.

In silico identification and characterization of microRNAs and their putative target genes in Solanaceae plants.

MicroRNAs (miRNAs) are a class of small, single-stranded, noncoding RNAs ranging from 19 to 25 nucleotides. The miRNA control various cellular functions by negatively regulating gene expression at the post-transcriptional level. The miRNA regulation over their target genes has a central role in regulating plant growth and development; however, only a few reports have been published on the function of miRNAs in the family Solanaceae. We identified Solanaceae miRNAs and their target genes by analyzing expressed sequence tag (EST) data from five different Solanaceae species. A comprehensive bioinformatic analysis of EST data of Solanaceae species revealed the presence of at least 11 miRNAs and 54 target genes in pepper (Capsicum annuum L.), 22 miRNAs and 221 target genes in potato (Solanum tuberosum L.), 12 miRNAs and 417 target genes in tomato (Solanum lycopersicum L.), 46 miRNAs and 60 target genes in tobacco (Nicotiana tabacum L.), and 7 miRNAs and 28 target genes in Nicotiana benthamiana. The identified Solanaceae miRNAs and their target genes were deposited in the SolmiRNA database, which is freely available for academic research only at http://genepool.kribb.re.kr/SolmiRNA. Our data indicate that the Solanaceae family has both conserved and specific miRNAs and that their target genes may play important roles in growth and development of Solanaceae plants.

Complete sequencing and comparative analyses of the pepper (Capsicum annuum L.) plastome revealed high frequency of tandem repeats and large insertion/deletions on pepper plastome.

Plants in the family Solanaceae are used as model systems in comparative and evolutionary genomics. The complete chloroplast genomes of seven solanaceous species have been sequenced, including tobacco, potato and tomato, but not peppers. We analyzed the complete chloroplast genome sequence of the hot pepper, Capsicum annuum. The pepper chloroplast genome was 156,781 bp in length, including a pair of inverted repeats (IR) of 25,783 bp. The content and the order of 133 genes in the pepper chloroplast genome were identical to those of other solanaceous plastomes. To characterize pepper plastome sequence, we performed comparative analysis using complete plastome sequences of pepper and seven solanaceous plastomes. Frequency and contents of large indels and tandem repeat sequences and distribution pattern of genome-wide sequence variations were investigated. In addition, a phylogenetic analysis using concatenated alignments of coding sequences was performed to determine evolutionary position of pepper in Solanaceae. Our results revealed two distinct features of pepper plastome compared to other solanaceous plastomes. Firstly, large indels, including insertions on accD and rpl20 gene sequences, were predominantly detected in the pepper plastome compared to other solanaceous plastomes. Secondly, tandem repeat sequences were particularly frequent in the pepper plastome. Taken together, our study represents unique features of evolution of pepper plastome among solanaceous plastomes.

Identification of functional BrFAD2-1 gene encoding microsomal delta-12 fatty acid desaturase from Brassica rapa and development of Brassica napus containing high oleic acid contents.

Microsomal delta-12 fatty acid desaturase (FAD2) functions in the first committed step of the biosynthesis of polyunsaturated fatty acids via the desaturation of oleic acid to linoleic acid. In this study, two FAD2 genes were identified through genome-wide analysis of Brassica rapa. One BrFAD2-1 gene harbors functional sequence information, but another BrFAD2-2 gene has mutations that generated a premature stop codon, rendering it nonfunctional. From a database of 120,000 B. rapa expressed sequence tags, we determined that all sequences coding for FAD2 corresponded to the BrFAD2-1 gene. The BrFAD2-1 protein was shown to share high sequence homology (71-99%) with FAD2 proteins from other plant species. An intron in the 5'-untranslated region and three histidine boxes in the protein, which are characteristic of plant FAD2 genes, have been well-conserved. BrFAD2-1 transcripts were detected in various organs of B. rapa. When a pBrFAD2-1:mRFP construct was introduced into tobacco epidermal cells, the fluorescent signal was noted in the endoplasmic reticulum. Ectopic expression of BrFAD2-1:mRFP complemented the Arabidopsis fad2-2 mutant. Finally, transgenic Korean rapeseed Tammi containing high oleic acid contents (78 mol%) was developed via the expression of the BrFAD2-1 gene in an antisense orientation. The data demonstrate that B. rapa harbors only one functional FAD2 that can be utilized for the development of the high-oleic acid Korean rapeseed cultivar Tammi, which might be useful for both human consumption and industrial applications.

Genome sequence of the polymyxin-producing plant-probiotic rhizobacterium Paenibacillus polymyxa E681.

Paenibacillus polymyxa E681, a spore-forming, low-G+C, Gram-positive bacterium isolated from the rhizosphere of winter barley grown in South Korea, has great potential for agricultural applications due to its ability to promote plant growth and suppress plant diseases. Here we present the complete genome sequence of P. polymyxa E681. Its 5.4-Mb genome encodes functions specialized to the plant-associated lifestyle and characteristics that are beneficial to plants, such as the production of a plant growth hormone, antibiotics, and hydrolytic enzymes.

Genome sequence of the probiotic bacterium Bifidobacterium animalis subsp. lactis AD011.

Bifidobacterium animalis subsp. lactis is a probiotic bacterium that naturally inhabits the guts of most mammals, including humans. Here we report the complete genome sequence of B. animalis subsp. lactis AD011 that was isolated from an infant fecal sample. Biological functions encoded in a single circular chromosome of 1,933,695 bp, smallest among the completely sequenced bifidobacterial genomes, are suggestive of their probiotic functions, such as utilization of bifidogenic factors and a variety of glycosidic enzymes and biosynthesis of polysaccharides.

Evolution of ribosomal DNA-derived satellite repeat in tomato genome.

BACKGROUND: Tandemly repeated DNA, also called as satellite DNA, is a common feature of eukaryotic genomes. Satellite repeats can expand and contract dramatically, which may cause genome size variation among genetically-related species. However, the origin and expansion mechanism are not clear yet and needed to be elucidated. RESULTS: FISH analysis revealed that the satellite repeat showing homology with intergenic spacer (IGS) of rDNA present in the tomato genome. By comparing the sequences representing distinct stages in the divergence of rDNA repeat with those of canonical rDNA arrays, the molecular mechanism of the evolution of satellite repeat is described. Comprehensive sequence analysis and phylogenetic analysis demonstrated that a long terminal repeat retrotransposon was interrupted into each copy of the 18S rDNA and polymerized by recombination rather than transposition via an RNA intermediate. The repeat was expanded through doubling the number of IGS into the 25S rRNA gene, and also greatly increasing the copy number of type I subrepeat in the IGS of 25-18S rDNA by segmental duplication. Homogenization to a single type of subrepeat in the satellite repeat was achieved as the result of amplifying copy number of the type I subrepeat but eliminating neighboring sequences including the type II subrepeat and rRNA coding sequence from the array. FISH analysis revealed that the satellite repeats are commonly present in closely-related Solanum species, but vary in their distribution and abundance among species. CONCLUSION: These results represent that the dynamic satellite repeats were originated from intergenic spacer of rDNA unit in the tomato genome. This result could serve as an example towards understanding the initiation and the expansion of the satellite repeats in complex eukaryotic genome.

Towards pathogenomics: a web-based resource for pathogenicity islands.

Pathogenicity islands (PAIs) are genetic elements whose products are essential to the process of disease development. They have been horizontally (laterally) transferred from other microbes and are important in evolution of pathogenesis. In this study, a comprehensive database and search engines specialized for PAIs were established. The pathogenicity island database (PAIDB) is a comprehensive relational database of all the reported PAIs and potential PAI regions which were predicted by a method that combines feature-based analysis and similarity-based analysis. Also, using the PAI Finder search application, a multi-sequence query can be analyzed onsite for the presence of potential PAIs. As of April 2006, PAIDB contains 112 types of PAIs and 889 GenBank accessions containing either partial or all PAI loci previously reported in the literature, which are present in 497 strains of pathogenic bacteria. The database also offers 310 candidate PAIs predicted from 118 sequenced prokaryotic genomes. With the increasing number of prokaryotic genomes without functional inference and sequenced genetic regions of suspected involvement in diseases, this web-based, user-friendly resource has the potential to be of significant use in pathogenomics. PAIDB is freely accessible at http://www.gem.re.kr/paidb.

Complete genome sequence of Leuconostoc citreum KM20.

Leuconostoc citreum is one of the most prevalent lactic acid bacteria during the manufacturing process of kimchi, the best-known Korean traditional dish. We have determined the complete genome sequence of L. citreum KM20. It consists of a 1.80-Mb chromosome and four circular plasmids and reveals genes likely involved in kimchi fermentation and its probiotic effects.

ASPMF: a new approach for identifying alternative splicing isoforms using peptide mass fingerprinting.

Alternative splicing is generally accepted as a mechanism that explains the discrepancy between the number of genes and proteins. We used peptide mass fingerprinting with a theoretical database and scoring method to discover and identify alternative splicing isoforms. Our theoretical database was built using published alternative splicing databases such as ECgene, H-DBAS, and TISA. According to our theoretical database of 190,529 isoforms, 37% of human genes have multiple isoforms. The isoforms produced from a gene partially share common peptide fragments because they have common exons, making it difficult to distinguish isoforms. Therefore, we developed a new method that effectively distinguishes a true isoform among multiple isoforms in a gene. In order to evaluate our algorithm, we made test sets for 4226 protein isoforms extracted from our theoretical database randomly. Consequently, 94% of true isoforms were identified by our scoring algorithm.

Pepper EST database: comprehensive in silico tool for analyzing the chili pepper (Capsicum annuum) transcriptome.

BACKGROUND: There is no dedicated database available for Expressed Sequence Tags (EST) of the chili pepper (Capsicum annuum), although the interest in a chili pepper EST database is increasing internationally due to the nutritional, economic, and pharmaceutical value of the plant. Recent advances in high-throughput sequencing of the ESTs of chili pepper cv. Bukang have produced hundreds of thousands of complementary DNA (cDNA) sequences. Therefore, a chili pepper EST database was designed and constructed to enable comprehensive analysis of chili pepper gene expression in response to biotic and abiotic stresses. RESULTS: We built the Pepper EST database to mine the complexity of chili pepper ESTs. The database was built on 122,582 sequenced ESTs and 116,412 refined ESTs from 21 pepper EST libraries. The ESTs were clustered and assembled into virtual consensus cDNAs and the cDNAs were assigned to metabolic pathway, Gene Ontology (GO), and MIPS Functional Catalogue (FunCat). The Pepper EST database is designed to provide a workbench for (i) identifying unigenes in pepper plants, (ii) analyzing expression patterns in different developmental tissues and under conditions of stress, and (iii) comparing the ESTs with those of other members of the Solanaceae family. The Pepper EST database is freely available at http://genepool.kribb.re.kr/pepper/. CONCLUSION: The Pepper EST database is expected to provide a high-quality resource, which will contribute to gaining a systemic understanding of plant diseases and facilitate genetics-based population studies. The database is also expected to contribute to analysis of gene synteny as part of the chili pepper sequencing project by mapping ESTs to the genome.

Isolation and characterization of multiple abundant lipid transfer protein isoforms in developing sesame (Sesamum indicum L.) seeds.

Sesame (Sesamum indicum) is an important oilseed crop; approximately 50% of the seed dry weight is storage oil. In a previous report, developing sesame seed expressed sequence tags (ESTs) revealed that ESTs encoding lipid transfer protein (LTPs) were one of the most abundant groups of sesame ESTs. LTP functions in the transfer of wax or cutin monomers and in the defense response against pathogen attack. To study the biological role of the abundant LTP isoforms in developing seeds, 122 ESTs out of 3328 sesame ESTs were analyzed against Arabidopsis and rice proteome databases. LTP fraction, which was partially purified from developing sesame seeds, actively transferred fluorescent phospholipids and bound to fatty acids. Full-length cDNAs of five out of 21 LTP isoforms were isolated and named SiLTP1-SiLTP5. The predicted amino acid sequences of the five SiLTPs harbor typical characteristics of LTPs, including conserved arrangement of cysteine residues. Northern blot analysis revealed that the five SiLTP isoforms were most abundantly expressed in developing seeds, but were also detected in flower tissues. Also, SiLTP3 and SiLTP4 transcripts were expressed in leaves and seed-pot walls, respectively. In addition, SiLTP2 and SiLTP4 transcripts were significantly induced in 6-day-old sesame seedlings by application of NaCl, mannitol, and abscisic acid (ABA). Transient expression of green fluorescent protein (GFP)-fusion constructs in Arabidopsis protoplasts revealed that SiLTP1 and SiLTP2 were secreted by different pathways. Taken together, the abundant LTPs in developing sesame seeds are involved in lipid transfer into the extracellular matrix. Possible biological roles of SiLTPs related to organ-specific expression and abiotic stresses are discussed.

A combined approach for the classification of G protein-coupled receptors and its application to detect GPCR splice variants.

G protein-coupled receptors (GPCRs) constitute the largest family of cell surface receptors and play a central role in cellular signaling pathways. The importance of GPCRs has led to their becoming the targets of more than 50% of prescription drugs. However, drug compounds that do not differentiate between receptor subtypes can have considerable side effects and efficacy problems. An accurate classification of GPCRs can solve the side effect problems and raise the efficacy of drugs. Here, we introduce an approach that combines a fingerprint method, statistical profiles and physicochemical properties of transmembrane (TM) domains for a highly accurate classification of the receptors. The approach allows both the recognition and classification for GPCRs at the subfamily and subtype level, and allows the identification of splice variants. We found that the approach demonstrates an overall accuracy of 97.88% for subfamily classification, and 94.57% for subtype classification.

The genome sequence of the capnophilic rumen bacterium Mannheimia succiniciproducens.

The rumen represents the first section of a ruminant animal's stomach, where feed is collected and mixed with microorganisms for initial digestion. The major gas produced in the rumen is CO(2) (65.5 mol%), yet the metabolic characteristics of capnophilic (CO(2)-loving) microorganisms are not well understood. Here we report the 2,314,078 base pair genome sequence of Mannheimia succiniciproducens MBEL55E, a recently isolated capnophilic Gram-negative bacterium from bovine rumen, and analyze its genome contents and metabolic characteristics. The metabolism of M. succiniciproducens was found to be well adapted to the oxygen-free rumen by using fumarate as a major electron acceptor. Genome-scale metabolic flux analysis indicated that CO(2) is important for the carboxylation of phosphoenolpyruvate to oxaloacetate, which is converted to succinic acid by the reductive tricarboxylic acid cycle and menaquinone systems. This characteristic metabolism allows highly efficient production of succinic acid, an important four-carbon industrial chemical.

Genomic blueprint of Hahella chejuensis, a marine microbe producing an algicidal agent.

Harmful algal blooms, caused by rapid growth and accumulation of certain microalgae in the ocean, pose considerable impacts on marine environments, aquatic industries and even public health. Here, we present the 7.2-megabase genome of the marine bacterium Hahella chejuensis including genes responsible for the biosynthesis of a pigment which has the lytic activity against a red-tide dinoflagellate. H.chejuensis is the first sequenced species in the Oceanospiralles clade, and sequence analysis revealed its distant relationship to the Pseudomonas group. The genome was well equipped with genes for basic metabolic capabilities and contained a large number of genes involved in regulation or transport as well as with characteristics as a marine heterotroph. Sequence analysis also revealed a multitude of genes of functional equivalence or of possible foreign origin. Functions encoded in the genomic islands include biosynthesis of exopolysacchrides, toxins, polyketides or non-ribosomal peptides, iron utilization, motility, type III protein secretion and pigmentation. Molecular structure of the algicidal pigment, which was determined through LC-ESI-MS/MS and NMR analyses, indicated that it is prodigiosin. In conclusion, our work provides new insights into mitigating algal blooms in addition to genetic make-up, physiology, biotic interactions and biological roles in the community of a marine bacterium.