Clusters Of Orthologous Groups Categories Research Articles

Exploring the impact of paracetamol on the transcriptome of Streptococcus pneumoniae D39

IntroductionParacetamol (acetaminophen) is a medicine used for the treatment of fever, pain, and inflammation during pneumococcal infection. ObjectivesTo see how paracetamol affects the transcriptional profile of Streptococcus pneumonia across the genome. MethodsIn this study, microarray analysis was performed for transcriptional profiling. ResultsTranscriptome data showed differential expression of several genes in S. pneumoniae D39 wild-type incorporated with paracetamol in the growth medium. Furthermore, these genes were categorized using Clusters of Orthologous Groups (COG) functional categorization on the basis of the suspected functions of the respective proteins. The majority of differentially expressed genes are in COG categories E (Amino acid transport and metabolism) and I (Lipid transport and metabolism). Analysis of protein–protein interaction networks exhibited compactly connected networks between fatty acid transport/biosynthesis and antibiotic biosynthesis genes. Moreover, pathways enrichment analysis revealed that fatty acid metabolism and biosynthesis genes were significantly affected under the criteria we've established. ConclusionThese results suggest the fatty acid biosynthesis and metabolism genes to be potential target of paracetamol in S. pneumoniae D39 wild-type.

Intrauterine Hypoxia Changed the Colonization of the Gut Microbiota in Newborn Rats.

Background: Accumulating evidence suggests a connection between the gut microbiota and neonatal diseases. Hypoxia may play an important role in the intestinal lesions in neonates.Objective: This study aims to determine whether the gut microbiota differs between intrauterine hypoxic rats and healthy controls and to identify the factors that influence the changes in the gut microbiota.Methods: We constructed an intrauterine hypoxia model in rats and collected the intestinal contents of intrauterine hypoxic newborn rats and normal newborn rats within 4 h and on the seventh day after birth. They were divided them into the intrauterine hypoxia first-day group (INH1), intrauterine hypoxia seventh-day group (INH7), normal first-day group (NOR1), and normal seventh-day group (NOR7). The contents of the intestines were sequenced with 16S rRNA sequencing, the sequencing results were analyzed for biological information, and the differences in the diversity, richness, and individual taxa among the groups were analyzed.Results: The abundance of the gut microbiota of neonatal rats with intrauterine hypoxia was higher than that of the control group rats. Intrauterine hypoxia altered the structural composition of the gut microbiota in neonatal rats. The INH1 group showed increased species richness, phylogenetic diversity, and β-diversity, and altered relative abundance in several taxa compared to those in the control group. The differences in the microbiota among the four groups were significantly higher than those within the group, and the differences in the abundance and diversity of the INH7 and NOR7 groups decreased after 7 days of suckling. Functional analysis based on the Cluster of Orthologous Groups (COG) suggested that 23 functional COG categories. There was no significant difference in the functional categories between the hypoxia group and the normal group.Conclusion: Intrauterine hypoxia changed the initial colonization of the gut microbiota in neonatal rats. It could increase the species richness and β-diversity of the gut microbiota, and altered relative abundances of several taxa.

Liraglutide-induced structural modulation of the gut microbiota in patients with type 2 diabetes mellitus.

Accumulating evidence has suggested the importance of gut microbiota in the development of type 2 diabetes mellitus (T2DM). In the present study, 40 patients with T2DM were treated with liraglutide for 4 months. Feces samples and clinical characteristics were collected from these 40 T2DM patients before and after the liraglutide treatment. The diversity and composition of gut microbiota in the two groups were determined by sequencing the V4 region of bacterial 16S rRNA genes. Meanwhile, blood glucose, insulin, hemoglobin A1c (HbA1c), and lipid metabolism were also measured in the pre- and post-liraglutide-treatment groups. We find that Baseline HbA1c was associated with liraglutide treatment response (R2 = 0.527, β = − 0.726, p < 0.0001). After adjusted for baseline HbA1c, blood urea nitrogen was associated with liraglutide treatment response. Besides, our results showed reduced gut microbial alpha diversity, different community structure distribution and altered microbial interaction network in patients treated with liraglutide. The liner discriminant analysis (LDA) effect size (LEfSe) analysis showed that 21 species of bacteria were abundant in the pre-liraglutide-treatment group and 15 species were abundant in the post-liraglutide-treatment group. In addition, we also find that Megamonas were significantly correlated with older age, diabetes duration and diabetic retinopathy, Clostridum were significantly correlated with family history of diabetes and Oscillospira were significantly correlated with both diabetic retinopathy and diabetic peripheral neuropathy. Functional analysis based on Kyoto Encyclopedia of Genes and Genomes (KEGG) and cluster of orthologous groups (COG) annotations enriched three KEGG metabolic pathways and six functional COG categories in the post-liraglutide-treatment group. In conclusion, our research suggests that baseline HbA1c, blood urea nitrogen and gut microbiota are associated with the liraglutide treatment applied on patients with T2DM. These findings may contribute to the beneficial effects of liraglutide against diabetes.

Universal and differential transcriptional regulatory pathways involved in the preparation of summer and winter diapauses in Pieris melete.

Much progress has been made in understanding the environmental and hormonal systems regulating winter diapause. However, transcriptional regulation of summer diapause is still largely unknown, making it difficult to understand an all-around regulation profile of seasonal adaptation. To bridge this gap, comparison RNA-seq to profile the transcriptome and to examine differential gene expression profiles between non-diapause, summer diapause, and winter diapause groups were performed. A total number of 113 million reads were generated and assembled into 79,117 unigenes, with 37,492 unigenes categorized into 58 functional gene ontology groups, 25 clusters of orthologous group categories, and 256 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. KEGG analysis mapped 2108 differentially expressed genes to 48 and 67 pathways for summer and winter diapauses, respectively. Enrichment statistics showed that 11 identical pathways similarly overlapped in the top 20 enriched functional groups both related to summer and winter diapauses. We also identified 35 key candidate genes for universal and differential functions related to summer and winter diapause preparation. Furthermore, we identified some genes involved in the signaling and metabolic pathways that may be the key drivers to integrate environmental signals into the summer and winter diapause preparation. The current study provided valuable insights into global molecular mechanisms underpinning diapause preparation.

Elevated CO2 weakens the shift in bacterial community structure in response to 8-year soybean straw return in the same experiment

A potentially important source for soil organic carbon (SOC) in the agricultural ecosystem is straw, straw return has been shown to affect soil bacterial communities. Facing global climate change, the response of bacterial communities to straw return will change at CO2 enrichment. In this study, we investigate the changes of bacterial communities in response to straw return (+straw) at elevated CO2 (eCO2, 700 ppm) in a long-term field experiment of northeast China. Soil samples were taken in the eighth year and analyzed by high throughput sequencing. Soil bacterial communities exhibited distinct clustering according to straw return and eCO2, while eCO2 shortened the distance of clustering between straw return and not. The relative abundances of 10 genera (Acidobacteria_norank, Candidatus_Solibacter, Gaiella, Nocardioides, Streptomyces, C0119_norank, Roseiflexus, Gemmatimonas, Mizugakiibacter and Rhodanobacter) were significantly affected by the interaction of straw × eCO2. In addition, straw return significantly decreased the relative abundances of Gaiellales_norank, Blastococcus, Psedarthrobacter, and Bacillus and increased that of Geminatimonadaceae_norank, Tepidisphaeraceae, Nitrosomonadaceae_norank, and SC-I-84_norank. These differential responses of genera abundances are illustrative of the susceptibility of bacterial communities and indicate their importance in evaluating the fate of exogenous C. The Clusters of Orthologous Groups (COG) analysis showed that straw return had a greater effect on the relative abundances of COG categories than eCO2. The present results point to the need to focus more strongly on the turnover and storage of straw-C during a chronic straw return in the future.

Integrated 16S rRNA Sequencing, Metagenomics, and Metabolomics to Characterize Gut Microbial Composition, Function, and Fecal Metabolic Phenotype in Non-obese Type 2 Diabetic Goto-Kakizaki Rats.

Type 2 diabetes mellitus (T2DM) is one of the most prevalent endocrine diseases in the world. Recent studies have shown that dysbiosis of the gut microbiota may be an important contributor to T2DM pathogenesis. However, the mechanisms underlying the roles of the gut microbiome and fecal metabolome in T2DM have not been characterized. Recently, the Goto-Kakizaki (GK) rat model of T2DM was developed to study the clinical symptoms and characteristics of human T2DM. To further characterize T2DM pathogenesis, we combined multi-omics techniques, including 16S rRNA gene sequencing, metagenomic sequencing, and metabolomics, to analyze gut microbial compositions and functions, and further characterize fecal metabolomic profiles in GK rats. Our results showed that gut microbial compositions were significantly altered in GK rats, as evidenced by reduced microbial diversity, altered microbial taxa distribution, and alterations in the interaction network of the gut microbiome. Functional analysis based on the cluster of orthologous groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations suggested that 5 functional COG categories belonged to the metabolism cluster and 33 KEGG pathways related to metabolic pathways were significantly enriched in GK rats. Metabolomics profiling identified 53 significantly differentially abundant metabolites in GK rats, including lipids and lipid-like molecules. These lipids were enriched in the glycerophospholipid metabolic pathway. Moreover, functional correlation analysis showed that some altered gut microbiota families, such as Verrucomicrobiaceae and Bacteroidaceae, significantly correlated with alterations in fecal metabolites. Collectively, the results suggested that an altered gut microbiota is associated with T2DM pathogenesis.

Molecular Signatures of Non-typeable Haemophilus influenzae Lung Adaptation in Pediatric Chronic Lung Disease.

Non-typeable Haemophilus influenzae (NTHi), an opportunistic pathogen of the upper airways of healthy children, can infect the lower airways, driving chronic lung disease. However, the molecular basis underpinning NTHi transition from a commensal to a pathogen is not clearly understood. Here, we performed comparative genomic and transcriptomic analyses of 12 paired, isogenic NTHi strains, isolated from the nasopharynx (NP) and bronchoalveolar lavage (BAL) of 11 children with chronic lung disease, to identify convergent molecular signatures associated with lung adaptation. Comparative genomic analyses of the 12 NP-BAL pairs demonstrated that five were genetically identical, with the remaining seven differing by only 1 to 3 mutations. Within-patient transcriptomic analyses identified between 2 and 58 differentially expressed genes in 8 of the 12 NP-BAL pairs, including pairs with no observable genomic changes. Whilst no convergence was observed at the gene level, functional enrichment analysis revealed significant under-representation of differentially expressed genes belonging to Coenzyme metabolism, Function unknown, Translation, ribosomal structure, and biogenesis Cluster of Orthologous Groups categories. In contrast, Carbohydrate transport and metabolism, Cell motility and secretion, Intracellular trafficking and secretion, and Energy production categories were over-represented. This observed trend amongst genetically unrelated NTHi strains provides evidence of convergent transcriptional adaptation of NTHi to pediatric airways that deserves further exploration. Understanding the pathoadaptative mechanisms that NTHi employs to infect and persist in the lower pediatric airways is essential for devising targeted diagnostics and treatments aimed at minimizing disease severity, and ultimately, preventing NTHi lung infections and subsequent chronic lung disease in children.

Analyses of Livestock-Associated Staphylococcus aureus Pan-Genomes Suggest Virulence Is Not Primary Interest in Evolution of Its Genome.

Staphylococcus aureus is not only part of normal flora but also an opportunistic pathogen relevant to microbial genomics, public health, and veterinary medicine. In addition to being a well-known human pathogen, S. aureus causes various infections in economically important livestock animals such as cows, sheep, goats, and pigs. There are very few studies that have examined the pan-genome of S. aureus or the host-specific strains' pan-genomes. We report on livestock-associated S. aureus' (LA-SA) pan-genome and suggest that virulence is not the primary interest in evolution of its genome. LA-SA' complete genomes were retrieved from the NCBI and pan-genome was constructed by high-speed Roary pipeline. The pan-genome size was 4637 clusters, whereas 42.46% of the pan-genome was associated with the core genome. We found 1268 genes were associated with the strain-unique genome, and the remaining 1432 cluster with the accessory genome. COG (clusters of orthologous group of proteins) analysis of the core genes revealed 34% of clusters related to metabolism responsible for amino acid and inorganic ion transport (COG categories E and P), followed by carbohydrate metabolism (category G). Virulent gene analysis revealed the core genes responsible for antiphagocytosis and iron uptake. The fluidity of pan-genome was calculated as 0.082 ± 0.025. Importantly, the positive selection analysis suggested a slower rate of evolution among the LA-SA genomes. We call for comparative microbial and pan-genome research between human and LA-SA that can help further understand the evolution of virulence and thus inform future microbial diagnostics and drug discovery.

The developmental transcriptome of Trichopria drosophilae (Hymenoptera: Diapriidae) and insights into cuticular protein genes

The pupal endoparasitoid wasp Trichopria drosophilae (Hymenoptera: Diapriidae) plays an important role in biological control of many frugivorous fruit fly species including Drosophila suzukii, a well-known invasive pest. Here, we report the transcriptomes of T. drosophilae among different developmental stages. A total of 601,148,438 high-quality reads were obtained and de novo assembled into 187,704 unigenes with an average length of 1096 bp. Among them, 21,735 unigenes were annotated into 52 Gene Ontology terms and 36,898 were assigned to 25 Cluster of Orthologous Groups categories, whereas 30,585 unigenes were mapped to 270 Kyoto Encyclopedia of Genes and Genomes different pathways. Numbers of differentially expressed genes were found through comparisons between different developmental stages. We further identified 137 cuticular protein genes (CPs) from T. drosophilae transcriptome, including 59 from CPR family, 2 from Tweedle family, 1 from CPF family, 46 from CPAP family, and 29 from other CP families. We analyzed expression patterns of the CPs at different developmental stages of T. drosophilae, and found some stage-specific CPs. Quantitative real-time PCR results confirmed RNA sequencing findings based on the relative expression levels of eight randomly selected CPs. This study provides a valuable transcriptomic resource for a comprehensive understanding of the development and physiology of T. drosophilae, and will help to improve their parasitism efficiency for biological control purposes.

English

Taro (Colocasia esculenta. (L.) Schott) is a genus of perennial plants that is widely distributed in the tropics or subtropics of Asia, Africa and America, which is the fourteenth most consumed vegetable of the world. However, molecular genetic research of Colocasia has been hindered by the insufficient genomic and transcriptome information. Here, the transcriptome of taro variety ‘Jingjiang Xiangsha’ from Jiangsu, China, was sequenced using the Illumina HiSeqTM 2000 platform in 2015. A total of 58,263,364 reads were generated, and assembly resolved into 65,878 unigenes with a N50 length of 1,357 bp. A total of 40,375 unigene sequences were successfully annotated based on searches against six public databases. Among the annotated unigenes, 14,753 were identified by gene Ontologyterms, 16,643 were classified to Clusters of Orthologous Groups categories, and 25,401 were mapped to 127 pathways in the Kyoto Encyclopedia of genes and genomes database. Also, 11,363 potential microsatellite loci were identified in 5,671 unigenes, and 150 primer pairs were randomly selected and amplified in 18 accessions of C. esculenta. A total of 100 primer pairs showed polymorphisms in repeat length. The number of alleles per locus ranged from 2 to 8. Across the 100 microsatellite loci, the polymorphism information content values ranged from 0.042 to 0.778. The transcriptomic data and microsatellite markers will play important roles in future functional gene analyses and genetic map construction of taro. Key words: Colocasia, genetic diversity, microsatellite markers, transcriptome.

Transcriptome analysis and development of simple sequence repeat (SSR) markers in Zingiber striolatum Diels.

Illumina-based paired-end sequencing technology was used for the high-throughput transcriptome sequencing of combined Zingiber striolatum Diels tissues (i.e., root, stem, leaf, flower, and fruit tissues). More than 130 million sequencing reads were generated, and a de novo assembly yielded 287,959 contigs and 112,107 unigenes with an average length of 1029 and 28,891bp, respectively. Similarity searches with known sequences led to the identification of 51,804 (46.21%) genes. Of the annotated unigenes, 6867 and 51,987 were assigned to Gene Ontology and Clusters of Orthologous Groups categories, respectively. Additionally, 8384 simple sequence repeats (SSRs) were identified as potential molecular markers in the unigenes. Thirty pairs of polymerase chain reaction primers were designed and used to validate the unigenes and assess the associated genomic polymorphism. The PCR amplification products for 25 primer pairs were ofthe expected size. These primers may represent usable molecular markers. The thousands of SSR markers identified in the present study may be useful for analyses of genetic diversity, genetic linkage mapping, and the identification and improvement of varieties during the breeding of Z. striolatum Diels. The unigene sequences and SSR markers described herein may serve as valuable resources for future investigations of Z. striolatum Diels.

Complete genome sequence of Citrobacter werkmanii strain BF-6 isolated from industrial putrefaction

BackgroundIn our previous study, Citrobacter werkmanii BF-6 was isolated from an industrial spoilage sample and demonstrated an excellent ability to form biofilms, which could be affected by various environmental factors. However, the genome sequence of this organism has not been reported so far.ResultsWe report the complete genome sequence of C. werkmanii BF-6 together with the description of the genome features and its annotation. The size of the complete chromosome is 4,929,789 bp with an average coverage of 137×. The chromosome exhibits an average G + C content of 52.0%, and encodes 4570 protein coding genes, 84 tRNA genes, 25 rRNA operons, 3 microsatellite sequences and 34 minisatellite sequences. A previously unknown circular plasmid designated as pCW001 was also found with a length of 212,549 bp and a G + C content of 48.2%. 73.5%, 75.6% and 92.6% of the protein coding genes could be assigned to GO Ontology, KEGG Pathway, and COG (Clusters of Orthologous Groups) categories respectively. C. werkmanii BF-6 and C. werkmanii NRBC 105721 exhibited the closest evolutionary relationships based on 16S ribosomal RNA and core-pan genome assay. Furthermore, C. werkmanii BF-6 exhibits typical bacterial biofilm formation and development. In the RT-PCR experiments, we found that a great number of biofilm related genes, such as bsmA, bssR, bssS, hmsP, tabA, csgA, csgB, csgC, csgD, csgE, and csgG, were involved in C. werkmanii BF-6 biofilm formation.ConclusionsThis is the first complete genome of C. werkmanii. Our work highlights the potential genetic mechanisms involved in biofilm formation and paves a way for further application of C. werkmanii in biofilms research.

Global effect of CsrA on gene expression in enterohemorrhagic Escherichia coli O157:H7

The post-transcriptional regulator CsrA regulates multiple unrelated processes such as central carbon metabolism, motility, biofilm formation and bacterial virulence in different bacteria. However, regulation by CsrA in enterohemorrhagic Escherichia coli (EHEC) O157:H7 is still largely unknown. In this study, we performed a detailed analysis of gene expression differences between the EHEC O157:H7 wild-type strain and a corresponding csrA::kan mutant using RNA-seq technology. Genes whose expression was affected by CsrA were identified and grouped into different clusters of orthologous group categories. Genes located in the locus of enterocyte effacement (LEE) pathogenicity island were significantly upregulated, whereas expression of flagella-related genes was significantly reduced in the csrA::kan mutant. Subsequent bacterial adherence and motility assays showed that inactivation of CsrA in EHEC O157:H7 resulted in a significant increase in bacterial adherence to host epithelial cells, with a concomitant loss of swimming motility on semi-solid agar plates. Furthermore, we also found that CsrA regulates genes not previously identified in other bacterial species, including genes encoding cytochrome oxidases and those required for nitrogen metabolism. Our results provide essential insight into the regulatory function of CsrA.

Potential factors involved in virulence of Cronobacter sakazakii isolates by comparative transcriptome analysis

Cronobacter species are important foodborne pathogens causing severe infections in neonates through consumption of contaminated powdered infant formula. However, the virulence-associated factors in Cronobacter are largely unknown. In this study, the transcriptome analysis between highly virulent Cronobacter sakazakii G362 and attenuated L3101 strains was used to reveal the potential factors involved in virulence. The total transcripts were grouped into 20 clusters of orthologous group categories and summarized in 3 gene ontology categories (biological process, cellular component, and molecular function). In addition, the differentially expressed genes (DEG) between these isolates were analyzed using Volcano plots and gene ontology enrichment. The predominant DEG were flagella-associated genes such as flhD, motA, flgM, flgB, and fliC. Furthermore, the expression abundance of outer membrane protein or lipoprotein genes (ompW, slyB, blc, tolC, and lolA), potential virulence-related factors (hlyIII and hha), and regulation factors (sdiA, cheY, Bss, fliZ) was also significantly different between G362 and L3101. Interestingly, 3 hypothetical protein genes (ESA_01022, ESA_01609, and ESA_00609) were found to be expressed only in G362. Our findings provide valuable transcriptomic information about potential virulence factor genes, which will be needed in future molecular biology studies designed to understand the pathogenic mechanism of Cronobacter.

De novo transcriptome analysis of Fraxinus velutina using Illumina platform and development of EST-SSR markers

To facilitate functional genomic analysis and molecular breeding of velvet ash (Fraxinus velutina Torr), the de novo sequencing was carried out by Illumina sequencing technology. The cDNA samples were prepared from eleven different tissues of velvet ash and sequenced by using the Illumina genome analyzer. Subsequently, de novo assemebly, functional annotation databases, and the screening of expressed sequence tag-simple sequence repeats (EST-SSRs) were performed by comparing with corresponding databases using BLASTx and software tools. We obtained 51 698 unigenes with an average length of 661 bp and an N50 length of 980 bp. Among all these unigenes, 41 267 (79.8 %) were annotated in the NCBI non-redundant protein database and 25 236 (48.8 %) were annotated in the Swiss-Prot database. A total of 31 546 (61.0 %) and 13 281 (25.7 %) unigenes were successfully categorized to 59 and 25 functional groups, respectively, by gene ontology categories and clusters of orthologous group categories. A total of 22 323 (43.2 %) unigenes were assigned to 128 pathways using the Kyoto encyclopedia of genes and genomes pathway database. Additionally, 3 249 EST-SSRs markers were detected in 51 698 unigenes from velvet ash. Based on 3 249 EST-SSRs markers, 1 800 primer pairs were successfully designed using Primer 3. In the 50 randomly selected primers, 48 successfully amplified fragments, and 42 showed polymorphisms. We completed a successful application of the Illumina platform to de novo transcriptome assembly of velvet ash, which has the potential to be used for discovering novel genes and further characterization of gene expression profiles.

Complete genome sequence and comparative genome analysis of the Paenibacillus mucilaginosus K02

AimPaenibacillus mucilaginosus (P. mucilaginosus) K02 is implicated in mineral weathering. However, relevant molecular mechanisms remain obscure. The study aims to uncover the bacterium's physiological processes using genomic approaches. Methods and resultsGenomic DNA from P. mucilaginosus K02 was sequenced using high-throughput Solexa sequencing technology and then conducted for Clusters of Orthologous Group (COG) annotation. Thereafter, genome sequences of K02 were compared with two strains, 3016 and KNP414. Mummer was applied for collinearity analysis of three P. mucilaginosus genomes. BLAST was used to identify pan and core genes in these strains. Finally, a phylogenetic tree was constructed using the maximum likelihood method by TreeBeST. Complete genome sequence of P. mucilaginosus K02 indicated the strain comprises one circular chromosome with 8,819,200 bases containing 58.3% GC content and 84.75% coding regions. A total of 7299 predicted ORFs were identified in the genome, among them, several genes were related to carbonic anhydrase (CA), and exopolysaccharide biosynthesis and secretion. Moreover, proteins of the predicted genes were annotated in COG categories such as “Carbohydrate transport and metabolism” and “Inorganic ion transport and metabolism.” In comparison with KNP414 and 3016, K02 exhibited chromosomal recombination or transposition. A total of 6662 core genes were identified among three P. mucilaginosus strains. The phylogenomic study indicated that P. mucilaginosus K02 was clustered with P. mucilaginosus strains 3016 and KNP414. ConclusionsIn P. mucilaginosus K02, genes related to CA and exopolysaccharide biosynthesis and secretion, and that involved in metabolism-related processes might play significant roles in mineral weathering.

De novo assembly and characterization of foot transcriptome and microsatellite marker development for Paphia textile

Paphia textile is an important, aquaculture bivalve clam species distributed mainly in China, Philippines, and Malaysia. Recent studies of P. textile have focused mainly on artificial breeding and nutrition analysis, and the transcriptome and genome of P. textile have rarely been reported. In this work, the transcriptome of P. textile foot tissue was sequenced on an Illumina HiSeq™ 2000 platform. A total of 20,219,795 reads were generated, resulting in 4.08Gb of raw data. The raw reads were cleaned and assembled into 54,852 unigenes with an N50 length of 829bp. Of these unigenes, 38.92% were successfully annotated based on their matches to sequences in seven public databases. Among the annotated unigenes, 14,571 were assigned Gene Ontology terms, 5448 were classified to Clusters of Orthologous Groups categories, and 6738 were mapped to 228 pathways in the Kyoto Encyclopedia of Genes and Genomes database. For functional marker development, 5605 candidate simple sequence repeats were identified in the transcriptome and 80 primer pairs were selected randomly and amplified in a wild population of P. textile. A total of 36 loci that exhibited obvious repeat length polymorphisms were detected. The transcriptomic data and microsatellite markers will provide valuable resources for future functional gene analyses, genetic map construction, and quantitative trait loci mapping in P. textile.

De novo assembly of the Japanese lawngrass (Zoysia japonica Steud.) root transcriptome and identification of candidate unigenes related to early responses under salt stress.

Japanese lawngrass (Zoysia japonica Steud.) is an important warm-season turfgrass that is able to survive in a range of soils, from infertile sands to clays, and to grow well under saline conditions. However, little is known about the molecular mechanisms involved in its resistance to salt stress. Here, we used high-throughput RNA sequencing (RNA-seq) to investigate the changes in gene expression of Zoysia grass at high NaCl concentrations. We first constructed two sequencing libraries, including control and NaCl-treated samples, and sequenced them using the Illumina HiSeq™ 2000 platform. Approximately 157.20 million paired-end reads with a total length of 68.68 Mb were obtained. Subsequently, 32,849 unigenes with an N50 length of 1781 bp were assembled using Trinity. Furthermore, three public databases, the Kyoto Encyclopedia of Genes and Genomes (KEGG), Swiss-prot, and Clusters of Orthologous Groups (COGs), were used for gene function analysis and enrichment. The annotated genes included 57 Gene Ontology (GO) terms, 120 KEGG pathways, and 24 COGs. Compared with the control, 1455 genes were significantly different (false discovery rate ≤0.01, |log2Ratio |≥1) in the NaCl-treated samples. These genes were enriched in 10 KEGG pathways and 73 GO terms, and subjected to 25 COG categories. Using high-throughput next-generation sequencing, we built a database as a global transcript resource for Z. japonica Steud. roots. The results of this study will advance our understanding of the early salt response in Japanese lawngrass roots.

Draft genome sequence of Methylibium sp. strain T29, a novel fuel oxygenate-degrading bacterial isolate from Hungary.

Methylibium sp. strain T29 was isolated from a gasoline-contaminated aquifer and proved to have excellent capabilities in degrading some common fuel oxygenates like methyl tert-butyl ether, tert-amyl methyl ether and tert-butyl alcohol along with other organic compounds. Here, we report the draft genome sequence of M. sp. strain T29 together with the description of the genome properties and its annotation. The draft genome consists of 608 contigs with a total size of 4,449,424 bp and an average coverage of 150×. The genome exhibits an average G + C content of 68.7 %, and contains 4754 protein coding and 52 RNA genes, including 48 tRNA genes. 71 % of the protein coding genes could be assigned to COG (Clusters of Orthologous Groups) categories. A formerly unknown circular plasmid designated as pT29A was isolated and sequenced separately and found to be 86,856 bp long.

A global transcriptional analysis of Megalobrama amblycephala revealing the molecular determinants of diet-induced hepatic steatosis

Blunt snout bream (Megalobrama amblycephala), a prevalent species in China's intensive polyculture systems, is highly susceptible to hepatic steatosis, resulting in considerable losses to the fish farming industry. Due to a lack of genomic resources, the molecular mechanisms of lipid metabolism in M. amblycephala are poorly understood. Here, a hepatic cDNA library was generated from equal amounts of mRNAs isolated from M. amblycephala fed normal-fat and high-fat diets. Sequencing of this library using the Illumina/Solexa platform produced approximately 51.87 million clean reads, which were assembled into 48,439 unigenes with an average length of 596bp and an N50 value of 800bp. These unigenes were searched against the nucleotide (NT), non-redundant (NR), Swiss-Prot, Cluster of Orthologous Groups (COG), and Kyoto Encyclopedia of Genes and Genome (KEGG) databases using the BLASTn or BLASTx algorithms (E-value≤10−5). A total of 8602 unigenes and 22,155 unigenes were functionally classified into 25 COG categories and 259 KEGG pathways, respectively. Furthermore, 22,072 unigenes were grouped into 62 sub-categories belonging to three main Gene Ontology (GO) terms. Using a digital gene expression analysis and the M. amblycephala transcriptome as a reference, 477 genes (134 up-regulated and 343 down-regulated) were identified as differentially expressed in fish fed a high-fat diet versus a normal-fat diet. KEGG and GO functional enrichment analyses of the differentially expressed unigenes were performed and 12 candidate genes related to lipid metabolism were identified. This study provides a global survey of hepatic transcriptome profiles and identifies candidate genes that may be related to lipid metabolism in M. amblycephala. These findings will facilitate further investigations of the mechanisms underlying hepatic steatosis in M. amblycephala.