BackgroundThe filamentous fungus Penicillium oxalicum is a potential alternative to Trichoderma reesei for industrial production of a complete cellulolytic enzyme system for a bio-refinery. Comparative omics approaches can support rational genetic engineering and/or breeding of filamentous fungi with improved cellulase production capacity. In this study, comparative genomic, transcriptomic and secretomic profiling of P. oxalicum HP7-1 and its cellulase and xylanase hyper-producing mutant EU2106 were employed to screen for novel regulators of cellulase and xylanase gene expression.ResultsThe 30.62 Mb P. oxalicum HP7-1 genome was sequenced, and 9834 protein-coding genes were annotated. Re-sequencing of the mutant EU2106 genome identified 274 single nucleotide variations and 12 insertion/deletions. Comparative genomic, transcriptomic and secretomic profiling of HP7-1 and EU2106 revealed four candidate regulators of cellulase and xylanase gene expression. Deletion of these candidate genes and measurement of the enzymatic activity of the resultant mutants confirmed the identity of three regulatory genes. POX02484 and POX08522, encoding a putative Zn(II)2Cys6 DNA-binding domain and forkhead protein, respectively, were found to be novel, while PoxClrB is an ortholog of ClrB, a key transcriptional regulator of cellulolytic enzyme gene expression in filamentous fungi. ΔPOX02484 and ΔPOX08522 mutants exhibited significantly reduced β-glucosidase activity, increased carboxymethylcellulose cellulase and xylanase activities, and altered transcription level of cellulase and xylanase genes compared with the parent strain ΔPoxKu70, with Avicel as the sole carbon source.ConclusionsTwo novel genes, POX02484 and POX08522, were found and characterized to regulate the expression of cellulase and xylanase genes in P. oxalicum. These findings are important for engineering filamentous fungi to improve cellulase and xylanase production.Electronic supplementary materialThe online version of this article (doi:10.1186/s13068-016-0616-9) contains supplementary material, which is available to authorized users.

The large yellow croaker Larimichthys crocea (L. crocea) is one of the most economically important marine fish in China and East Asian countries. It also exhibits peculiar behavioral and physiological characteristics, especially sensitive to various environmental stresses, such as hypoxia and air exposure. These traits may render L. crocea a good model for investigating the response mechanisms to environmental stress. To understand the molecular and genetic mechanisms underlying the adaptation and response of L. crocea to environmental stress, we sequenced and assembled the genome of L. crocea using a bacterial artificial chromosome and whole-genome shotgun hierarchical strategy. The final genome assembly was 679 Mb, with a contig N50 of 63.11 kb and a scaffold N50 of 1.03 Mb, containing 25,401 protein-coding genes. Gene families underlying adaptive behaviours, such as vision-related crystallins, olfactory receptors, and auditory sense-related genes, were significantly expanded in the genome of L. crocea relative to those of other vertebrates. Transcriptome analyses of the hypoxia-exposed L. crocea brain revealed new aspects of neuro-endocrine-immune/metabolism regulatory networks that may help the fish to avoid cerebral inflammatory injury and maintain energy balance under hypoxia. Proteomics data demonstrate that skin mucus of the air-exposed L. crocea had a complex composition, with an unexpectedly high number of proteins (3,209), suggesting its multiple protective mechanisms involved in antioxidant functions, oxygen transport, immune defence, and osmotic and ionic regulation. Our results reveal the molecular and genetic basis of fish adaptation and response to hypoxia and air exposure. The data generated by this study will provide valuable resources for the genetic improvement of stress resistance and yield potential in L. crocea.

Objective: Leaflet thickening, fibrosis, and hardening are early pathological features of calcific aortic valve disease (CAVD). An inadequate understanding of the resident aortic valve cells involved in the pathological process may compromise the development of therapeutic strategies. We aim to construct a pattern of the human aortic valve cell atlas in healthy and CAVD clinical specimens, providing insight into the cellular origins of CAVD and the complex cytopathological differentiation process. Approach and Results: We used unbiased single-cell RNA sequencing for the high-throughput evaluation of cell heterogeneity in 34 632 cells isolated from 6 different human aortic valve leaflets. Cellular experiments, in situ localization, and bulk sequencing were performed to verify the differences between normal, healthy valves and those with CAVD. By comparing healthy and CAVD specimens, we identified 14 cell subtypes, including 3 heterogeneous subpopulations of resident valve interstitial cells, 3 types of immune-derived cells, 2 types of valve endothelial cells, and 6 novel valve-derived stromal cells found particularly in CAVD leaflets. Combining additional verification experiments with single-cell transcriptome profiling provided evidence of endothelial to mesenchymal transition involved in lesion thickening of the aortic valve leaflet. Conclusions: Our findings deconstructed the aortic valve cell atlas and suggested novel functional interactions among resident cell subpopulations. Our findings may provide insight into future targeted therapies to prevent CAVD.

This study was conducted to (1) characterize coagulation cascade and complement system in systemic lupus erythematosus (SLE); (2) evaluate the associations between coagulation cascade, complement system, inflammatory response and SLE disease severity; (3) test the diagnostic value of a combination of D-dimer and C4 for lupus activity. Transcriptomics, proteomics and metabolomics were performed in 24 SLE patients and 24 healthy controls. The levels of ten coagulations, seven complements and three cytokines were measured in 112 SLE patients. Clinical data were collected from 2025 SLE patients. The analysis of multi-omics data revealed the common links for the components of coagulation cascade and complement system. The results of ELISA showed coagulation cascade and complement system had an interaction effect on SLE disease severity, this effect was pronounced among patients with excess inflammation. The analysis of clinical data revealed a combination of D-dimer and C4 provided good diagnostic performance for lupus activity. This study suggested that coagulation cascade and complement system become ‘partners in crime’, contributing to SLE disease severity and identified the diagnostic value of D-dimer combined with C4for lupus activity.

The diamondback moth, Plutella xylostella (L.), is the major cosmopolitan pest of brassica and other cruciferous crops. Its larval midgut is a dynamic tissue that interfaces with a wide variety of toxicological and physiological processes. The draft sequence of the P. xylostella genome was recently released, but its annotation remains challenging because of the low sequence coverage of this branch of life and the poor description of exon/intron splicing rules for these insects. Peptide sequencing by computational assignment of tandem mass spectra to genome sequence information provides an experimental independent approach for confirming or refuting protein predictions, a concept that has been termed proteogenomics. In this study, we carried out an in-depth proteogenomic analysis to complement genome annotation of P. xylostella larval midgut based on shotgun HPLC-ESI-MS/MS data by means of a multialgorithm pipeline. A total of 876,341 tandem mass spectra were searched against the predicted P. xylostella protein sequences and a whole-genome six-frame translation database. Based on a data set comprising 2694 novel genome search specific peptides, we discovered 439 novel protein-coding genes and corrected 128 existing gene models. To get the most accurate data to seed further insect genome annotation, more than half of the novel protein-coding genes, i.e. 235 over 439, were further validated after RT-PCR amplification and sequencing of the corresponding transcripts. Furthermore, we validated 53 novel alternative splicings. Finally, a total of 6764 proteins were identified, resulting in one of the most comprehensive proteogenomic study of a nonmodel animal. As the first tissue-specific proteogenomics analysis of P. xylostella, this study provides the fundamental basis for high-throughput proteomics and functional genomics approaches aimed at deciphering the molecular mechanisms of resistance and controlling this pest. Thousands of eukaryotic genomes have been sequenced with the rapid development of massive parallel sequencing technologies. The annotation of their genomic DNA sequences is a prerequisite to delineate conserved elements and understand the specific encoded functions. Today, conventional approaches used to identify protein-coding genes are highly dependent on predictions based on computational algorithms and homology searches against known proteins. Nonmodel organisms, such as most arthropods, are by nature distantly related to well-studied organisms (1). As a result, their genomes encode a large majority of proteins with little similarity to known proteins, challenging the annotation process. Much recent focus on computational gene finding is using transcript evidence to complement this predictionbased approach. Indeed, the use of information from deep sequencing of mRNA-derived cDNA libraries (RNA-seq) or expressed sequence tag (EST) 1 libraries can dramatically improve the genome annotation confidence (2-4). However, this analysis remains at the transcript level and cannot make a crystal-clear distinct...

Abiotic stress adversely inhibits the growth and development of plants, by changing the expression of multiple genes. Circular RNAs (circRNAs), as a class of non-coding RNAs, function in transcriptional and posttranscriptional regulation. Yet, the involvement of circRNAs in abiotic stress response is rarely reported. In this study, the participation and function of circRNAs in low-temperature (LT)-induced stress response were investigated in tomato leaves. We generated genome-wide profiles of circRNAs and mRNAs in tomato leaves grown at 25 • C room temperature (RT) and 12 • C LT. Our results show that 1,830 circRNAs were identified in tomato leaves in both RT and LT treatments, among which 1,759 were differentially induced by the LT treatment. We find that the identified circRNAs are mainly located at exons of genes, but less distributed at introns of genes or intergenic regions. Our results suggest that there are 383 differentially expressed circRNAs predicted to function as putative sponges of 266 miRNAs to target 4,476 mRNAs in total. Moreover, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis assays indicate that multiple pathways were enriched in both differentially expressed genes induced by LT and parental genes of differentially expressed circRNAs induced by LT, revealing the key functions of circRNAs and the corresponding targeted genes in response to LT stress. Our results suggest that circRNAs may be involved in regulating metabolism (i.e., carbohydrate, amino acid, lipid, and energy), signal transduction, and environmental adaptation-related pathways and that these circRNAs were predicted to regulate the expression of transcription factors, genes in signal transduction pathways, and genes related to the Ca 2+ channel through targeting the corresponding proteins, such as WRKY, NAC, cytochrome P450, and calmodulin binding protein. Taken together, our study uncovers that multiple circRNAs are isolated and differently regulated in response to LT stress and provides the resource and potential networks of circRNA-miRNA-mRNA under LT stress for further investigations in tomato leaves.

1The large yellow croaker Larimichthys crocea (L. crocea) is one of the most economically