Comparative pangenome analysis of Aspergillus flavus and Aspergillus oryzae reveals their phylogenetic, genomic, and metabolic homogeneity

Aspergillus flavus and Aspergillus oryzae are closely related fungal species with contrasting roles in food safety and fermentation. To comprehensively investigate their phylogenetic, genomic, and metabolic characteristics, we conducted an extensive comparative pangenome analysis using complete, dereplicated genome sets for both species. Phylogenetic analyses, employing both the entirety of the identified single-copy orthologous genes and six housekeeping genes commonly used for fungal classification, did not reveal clear differentiation between A. flavus and A. oryzae genomes. Upon analyzing the aflatoxin biosynthesis gene clusters within the genomes, we observed that non-aflatoxin-producing strains were dispersed throughout the phylogenetic tree, encompassing both A. flavus and A. oryzae strains. This suggests that aflatoxin production is not a distinguishing trait between the two species. Furthermore, A. oryzae and A. flavus strains displayed remarkably similar genomic attributes, including genome sizes, gene contents, and G + C contents, as well as metabolic features and pathways. The profiles of CAZyme genes and secondary metabolite biosynthesis gene clusters within the genomes of both species further highlight their similarity. Collectively, these findings challenge the conventional differentiation of A. flavus and A. oryzae as distinct species and highlight their phylogenetic, genomic, and metabolic homogeneity, potentially indicating that they may indeed belong to the same species.

BRD2-specific inhibitor, BBC0403, inhibits the progression of osteoarthritis pathogenesis in osteoarthritis-induced C57BL/6 male mice

BACKGROUND AND PURPOSE

The discovery of new bromo- and extra-terminal inhibitors presents new drugs to treat osteoarthritis (OA).

EXPERIMENTAL APPROACH

The new drug, BBC0403, was identified in the DNA-encoded library screening system by searching for compounds that target BRD (bromodomain-containing) proteins. The binding force with BRD proteins was evaluated using time-resolved fluorescence energy transfer (TR-FRET) and binding kinetics assays. Subsequently, in vitro and ex vivo analyses demonstrated the effects of the BRD2 inhibitor, BBC0403, on OA. For animal experiments, medial meniscus destabilization was performed to create a 12-week-old male C57BL/6 mouse model, and intra-articular (i.a.) injections were administered. Histological and immunohistochemical analyses were then performed. The underlying mechanism was confirmed by gene set enrichment analysis (GSEA) using RNA-seq.

KEY RESULTS

TR-FRET and binding kinetics assays revealed that BBC0403 exhibited higher binding specificity for BRD2 compared to BRD3 and BRD4. The anti-OA effects of BBC0403 were tested at concentrations of 5, 10 and 20 μM (no cell toxicity in the range tested). The expression of catabolic factors, prostaglandin E2 (PGE2) production and extracellular matrix (ECM) degradation was reduced. Additionally, the i.a. injection of BBC0403 prevented OA cartilage degradation in mice. Finally, BBC0403 was demonstrated to suppress NF-κB and MAPK signalling pathways.

CONCLUSION AND IMPLICATIONS

This study demonstrated that BBC0403 is a novel BRD2-specific inhibitor and a potential i.a.-injectable therapeutic agent to treat OA.

AMPK-HIF-1α signaling enhances glucose-derived de novo serine biosynthesis to promote glioblastoma growth

BACKGROUND

Cancer cells undergo cellular adaptation through metabolic reprogramming to sustain survival and rapid growth under various stress conditions. However, how brain tumors modulate their metabolic flexibility in the naturally serine/glycine (S/G)-deficient brain microenvironment remain unknown.

METHODS

We used a range of primary/stem-like and established glioblastoma (GBM) cell models in vitro and in vivo. To identify the regulatory mechanisms of S/G deprivation-induced metabolic flexibility, we employed high-throughput RNA-sequencing, transcriptomic analysis, metabolic flux analysis, metabolites analysis, chromatin immunoprecipitation (ChIP), luciferase reporter, nuclear fractionation, cycloheximide-chase, and glucose consumption. The clinical significances were analyzed in the genomic database (GSE4290) and in human GBM specimens.

RESULTS

The high-throughput RNA-sequencing and transcriptomic analysis demonstrate that the de novo serine synthesis pathway (SSP) and glycolysis are highly activated in GBM cells under S/G deprivation conditions. Mechanistically, S/G deprivation rapidly induces reactive oxygen species (ROS)-mediated AMP-activated protein kinase (AMPK) activation and AMPK-dependent hypoxia-inducible factor (HIF)-1α stabilization and transactivation. Activated HIF-1α in turn promotes the expression of SSP enzymes phosphoglycerate dehydrogenase (PHGDH), phosphoserine aminotransferase 1 (PSAT1), and phosphoserine phosphatase (PSPH). In addition, the HIF-1α-induced expression of glycolytic genes (GLUT1, GLUT3, HK2, and PFKFB2) promotes glucose uptake, glycolysis, and glycolytic flux to fuel SSP, leading to elevated de novo serine and glycine biosynthesis, NADPH/NADP+ ratio, and the proliferation and survival of GBM cells. Analyses of human GBM specimens reveal that the levels of overexpressed PHGDH, PSAT1, and PSPH are positively correlated with levels of AMPK T172 phosphorylation and HIF-1α expression and the poor prognosis of GBM patients.

CONCLUSION

Our findings reveal that metabolic stress-enhanced glucose-derived de novo serine biosynthesis is a critical metabolic feature of GBM cells, and highlight the potential to target SSP for treating human GBM.

Novel molecule BBC0901 inhibits BRD4 and acts as a catabolic regulator in the pathogenesis of osteoarthritis

Osteoarthritis (OA) is induced by matrix degradation and inflammation mediated by bromo-domain-containing protein 4 (BRD4)-dependent catabolic factors. BRD4 acts as both a transcriptional regulator and an epigenetic reader. BBC0901 was identified as an inhibitor of BRD4 using a DNA-encoded library screening system. We aimed to demonstrate the effects of BBC0901 on OA pathogenesis by in vitro, ex vivo, and in vivo analyses. BBC0901 inhibited the expression of catabolic factors that degrade cartilage without significantly affecting the viability of mouse articular chondrocytes. Additionally, ex vivo experiments under conditions mimicking OA showed that BBC0901 suppressed extracellular matrix degradation. RNA sequencing analysis of gene expression patterns showed that BBC0901 inhibited the expression of catabolic factors, such as matrix metalloproteinases (MMPs) and cyclooxygenase (COX)2, along with reactive oxygen species (ROS) production. Furthermore, intra-articular (IA) injection of BBC0901 into the knee joint blocked osteoarthritic cartilage destruction by inhibition of MMP3, MMP13, COX2, interleukin (IL)6, and ROS production, thereby obstructing the nuclear factor kappa-light-chain-enhancer of activated B cell and mitogen activated protein kinase signaling. In conclusion, BBC0901-mediated BRD4 inhibition prevented OA development by attenuating catabolic signaling and hence, can be considered a promising IA therapeutic for OA.

Blockade of Activin Receptor IIB Protects Arthritis Pathogenesis by Non-Amplification of Activin A-ACVR2B-NOX4 Axis Pathway

Although activin receptor IIB (ACVR2B) is emerging as a novel pathogenic receptor, its ligand and assembled components (or assembly) are totally unknown in the context of osteoarthritis (OA) pathogenesis. The present results suggest that upregulation of ACVR2B and its assembly could affect osteoarthritic cartilage destruction. It is shown that the ACVR2B ligand, activin A, regulates catabolic factor expression through ACVR2B in OA development. Activin A Tg mice (Col2a1-Inhba) exhibit enhanced cartilage destruction, whereas heterozygous activin A KO mice (Inhba^+/- ) show protection from cartilage destruction. In silico analysis suggests that the Activin A-ACVR2B axis is involved in Nox4-dependent ROS production. Activin A Tg:Nox4 KO (Col2a1-Inhba:Nox4^-/- ) mice show inhibition of experimental OA pathogenesis. NOX4 directly binds to the C-terminal binding site on ACVR2B-ACVR1B and amplifies the pathogenic signal for cartilage destruction through SMAD2/3 signaling. Together, the findings reveal that the ACVR2B assembly, which comprises Activin A, ACVR2B, ACVR1B, Nox4, and AP-1-induced HIF-2α, accelerates OA development. Furthermore, it is shown that shRNA-mediated ACVR2B knockdown or trapping ligands of ACVR2B abrogate OA development by competitively disrupting the ACVR2B-Activin A interaction. These results suggest that the ACVR2B assembly is required to amplify osteoarthritic cartilage destruction and could be a potential therapeutic target in efforts to treat OA.

A practical comparison of the next-generation sequencing platform and assemblers using yeast genome

Assembling fragmented whole-genomic information from the sequencing data is an inevitable process for further genome-wide research. However, it is intricate to select the appropriate assembly pipeline for unknown species because of the species-specific genomic properties. Therefore, our study focused on relatively more static proclivities of sequencing platforms and assembly algorithms than the fickle genome sequences. A total of 212 draft and polished de novo assemblies were constructed under the different sequencing platforms and assembly algorithms with the repetitive yeast genome. Our comprehensive data indicated that sequencing reads from Oxford Nanopore with R7.3 flow cells generated more continuous assemblies than those derived from the PacBio Sequel, although the homopolymer-based assembly errors and chimeric contigs exist. In addition, the comparison between two second-generation sequencing platforms showed that Illumina NovaSeq 6000 provides more accurate and continuous assembly in the second-generation-sequencing-first pipeline, but MGI DNBSEQ-T7 provides a cheap and accurate read in the polishing process. Furthermore, our insight into the relationship among the computational time, read length, and coverage depth provided clues to the optimal pipelines of yeast assembly.

Galectin-4 increases the ability of M2 macrophages to enhance antiviral CD4+ T-cell responses

Galectin-4 (Gal-4) is a β-galactoside-binding protein belonging to the galectin family. Although Gal-4 is known to be involved in several physiologic processes of the gastrointestinal tract, its immunomodulatory roles remain unclear. In this study, we investigated whether Gal-4 influences the function of M1 and M2 macrophages. Gal-4 treatment drove more robust changes in the gene expression of M2 macrophages compared to M1 macrophages. Antiviral immune response-related genes were significantly upregulated in Gal-4-treated M2 macrophages. Gal-4 significantly enhanced the immunostimulatory activity of M2 macrophages upon Toll-like receptor 7 stimulation or infection with lymphocytic choriomeningitis virus (LCMV). Moreover, the antibody production against LCMV infection and the antiviral CD4+ T-cell responses, but not the antiviral CD8+ T-cell responses, were greatly increased by Gal-4-treated M2 macrophages in vivo. The present results indicate that Gal-4 enhances the ability of M2 macrophages to promote antiviral CD4+ T-cell responses. Thus, Gal-4 could be used to boost antiviral immune responses.

Genomic and functional features of yeast species in Korean traditional fermented alcoholic beverage and soybean products

In this review, we describe the genomic and physiological features of the yeast species predominantly isolated from Nuruk, a starter for traditional Korean rice wines, and Jang, a traditional Korean fermented soy product. Nuruk and Jang have several prevalent yeast species, including Saccharomycopsis fibuligera, Hyphopichia burtonii, and Debaryomyces hansenii complex, which belong to the CUG clade showing high osmotic tolerance. Comparative genomics revealed that the interspecies hybridization within yeast species for generating heterozygous diploid genomes occurs frequently as an evolutional strategy in the fermentation environment of Nuruk and Jang. Through gene inventory analysis based on the high-quality reference genome of S. fibuligera, new genes involved in cellulose degradation and volatile aroma biosynthesis and applicable to the production of novel valuable enzymes and chemicals can be discovered. The integrated genomic and transcriptomic analysis of Hyphopichia yeasts, which exhibit strong halotolerance, provides insights into the novel mechanisms of salt and osmo-stress tolerance for survival in fermentation environments with a low-water activity and high-concentration salts. In addition, Jang yeast isolates, such as D. hansenii, show probiotic potential for the industrial application of yeast species beyond fermentation starters to diverse human health sectors.

Genomic features, aroma profiles, and probiotic potential of the Debaryomyces hansenii species complex strains isolated from Korean soybean fermented food

Fermented soybean products are gaining attention in the food industry owing to their nutritive value and health benefits. In this study, we performed genomic analysis and physiological characterization of two Debaryomyces spp. yeast isolates obtained from a Korean traditional fermented soy sauce "ganjang". Both Debaryomyces hansenii ganjang isolates KD2 and C11 showed halotolerance to concentrations of up to 15% NaCl and improved growth in the presence of salt. Ploidy and whole-genome sequencing analyses indicated that the KD2 genome is haploid, whereas the C11 genome is heterozygous diploid with two distinctive subgenomes. Interestingly, phylogenetic analysis using intron sequences indicated that the C11 strain was generated via hybridization between D. hansenii and D. tyrocola ancestor strains. The D. hansenii KD2 and D. hansenii-hybrid C11 produced various volatile flavor compounds associated with butter, caramel, cheese, and fruits, and showed high bioconversion activity from ferulic acid to 4-vinylguaiacol, a characteristic flavor compound of soybean products. Both KD2 and C11 exhibited viability in the presence of bile salts and at low pH and showed immunomodulatory activity to induce high levels of the anti-inflammatory cytokine IL-10. The safety of the yeast isolates was confirmed by analyzing virulence and acute oral toxicity. Together, the D. hansenii ganjang isolates possess physiological properties beneficial for improving the flavor and nutritional value of fermented products.

Development of DNA aptamers specific for small therapeutic peptides using a modified SELEX method

Aptamers are short single-stranded DNA or RNA oligonucleotides capable of binding with high affinity and specificity to target molecules. Because of their durability and ease of synthesis, aptamers are used in a wide range of biomedical fields, including the diagnosis of diseases and targeted delivery of therapeutic agents. The aptamers were selected using a process called systematic evolution of ligands by exponential enrichment (SELEX), which has been improved for various research purposes since its development in 1990. In this protocol, we describe a modified SELEX method that rapidly produces high aptamer screening yields using two types of magnetic beads. Using this method, we isolated an aptamer that specifically binds to an antimicrobial peptide. We suggest that by conjugating a small therapeutic-specific aptamer to a gold nanoparticle-based delivery system, which enhances the stability and intracellular delivery of peptides, aptamers selected by our method can be used for the development of therapeutic agents utilizing small therapeutic peptides.

A New Branchipolynoe (Aphroditiformia: Polynoidae) Scale Worm from the Onnuri Deep-sea Hydrothermal Vent Field, Northern Central Indian Ridge

Deep-sea hydrothermal vents are dynamic environments with exotic fauna, including bathymodiolin mussels and scale worm annelids that are often in close association. In this study, we found a new species of Branchipolynoe (Aphroditiformia: Polynoidae) living in the recently discovered mussel Gigantidas vrijenhoeki in deep-sea hydrothermal vents and methane seeps at 2,014-2,023 m depth. Based on the morphology and full mitochondrial genome sequences of specimens of Branchipolynoe from the Onnuri vent field (OVF) on the northern Central Indian Ridge, we describe them as a new species: Branchipolynoe onnuriensis sp. nov. This species resembles B. longqiensis and B. tjiasmantoi, but can be distinguished from these species by the shape of the notopodial acicular lobe and the tips of the subacicular neurochaetae. This identity is well-supported by genetic distance and phylogenetic analyses based on the mitochondrial c oxidase subunit I (COI) gene, with the new species being closest to the Western Pacific species B. tjiasmantoi. Phylogenetic analyses support close relationships between the Indian Ocean and Western Pacific hydrothermal polychaetes. Our data provide a foundation for exploring the evolutionary relationship between scale worms and bathymodiolin mussels.

Construction of a tunable promoter library to optimize gene expression in Methylomonas sp. DH-1, a methanotroph, and its application to cadaverine production

BACKGROUND

As methane is 84 times more potent than carbon dioxide in exacerbating the greenhouse effect, there is an increasing interest in the utilization of methanotrophic bacteria that can convert harmful methane into various value-added compounds. A recently isolated methanotroph, Methylomonas sp. DH-1, is a promising biofactory platform because of its relatively fast growth. However, the lack of genetic engineering tools hampers its wide use in the bioindustry.

RESULTS

Through three different approaches, we constructed a tunable promoter library comprising 33 promoters that can be used for the metabolic engineering of Methylomonas sp. DH-1. The library had an expression level of 0.24-410% when compared with the strength of the lac promoter. For practical application of the promoter library, we fine-tuned the expressions of cadA and cadB genes, required for cadaverine synthesis and export, respectively. The strain with P_rpmB-cadA and P_DnaA-cadB produced the highest cadaverine titre (18.12 ± 1.06 mg/L) in Methylomonas sp. DH-1, which was up to 2.8-fold higher than that obtained from a non-optimized strain. In addition, cell growth and lysine (a precursor of cadaverine) production assays suggested that gene expression optimization through transcription tuning can afford a balance between the growth and precursor supply.

CONCLUSIONS

The tunable promoter library provides standard and tunable components for gene expression, thereby facilitating the use of methanotrophs, specifically Methylomonas sp. DH-1, as a sustainable cell factory.

Evidence for reciprocal evolution of the global repressor Mlc and its cognate phosphotransferase system sugar transporter

Because the bacterial phosphoenolpyruvate:carbohydrate phosphotransferase system (PTS) is involved in the regulation of various physiological processes in addition to carbohydrate transport, its expression is precisely regulated in response to the availability of PTS sugars. The PTS consists of enzyme I and histidine phosphocarrier protein, and several sugar-specific enzymes II. In Escherichia coli, genes for enzymes II specific for glucose and related sugars are co-regulated by the global repressor Mlc, and glucose induction of the Mlc regulon genes is achieved by its interaction with glucose-specific enzyme II (EII^Glc ). In this study, we revealed that, in Vibrio species, which are phylogenetically older than Enterobacteriaceae, the membrane sequestration of Mlc and thereby the induction of its regulon genes is mediated by N-acetylglucosamine (NAG)-specific EII. While Vibrio Mlc interacts only with the EIIB domain of EII^Nag , E. coli Mlc interacts with the EIIB domain of both EII^Glc and EII^Nag . The present data suggest that EII^Nag may be the primordial regulator of Mlc, and EII^Glc has evolved to interact with Mlc since an EIIA domain was fused to EII^Nag in Enterobacteriaceae. Our findings provide insight into the coevolutionary dynamics between a transcription factor and its cognate regulator according to long-term resource availability in the bacterial natural habitat.

Synthetically engineered microbial scavengers for enhanced bioremediation

Microbial bioremediation has gained attention as a cheap, efficient, and sustainable technology to manage the increasing environmental pollution. Since microorganisms in nature are not evolved to degrade pollutants, there is an increasing demand for developing safer and more efficient pollutant-scavengers for enhanced bioremediation. In this review, we introduce the strategies and technologies developed in the field of synthetic biology and their applications to the construction of microbial scavengers with improved efficiency of biodegradation while minimizing the impact of genetically engineered microbial scavengers on ecosystems. In addition, we discuss recent achievements in the biodegradation of fastidious pollutants, greenhouse gases, and microplastics using engineered microbial scavengers. Using synthetic microbial scavengers and multidisciplinary technologies, toxic pollutants could be more easily eliminated, and the environment could be more efficiently recovered.

Estrogen-related receptor-gamma influences Helicobacter pylori infection by regulating TFF1 in gastric cancer

Helicobacter pylori infection is a crucial factor in the development of gastric cancer (GC). Molecular therapeutic targets and mechanisms contributing to H. pylori infection-associated GC induction are poorly understood and this study aimed to fill that research gap. We found that the nuclear receptor estrogen-related receptor gamma (ESRRG) is a candidate factor influencing H. pylori infection-driven GC. ESRRG suppressed H. pylori infection and cell growth induced by H. pylori infection in GC cells and organoid models In addition, H. pylori infection downregulates ESRRG expression. Gene expression profiling revealed that trefoil factor 1 (TFF1), a well-known tumor suppressor in GC, is a downstream target of ESRRG. Mechanistically, ESRRG directly binds to the TFF1 promoter and induces TFF1 gene expression. Furthermore, TFF1 activation by ESRRG was inhibited by nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB)/p65, which is induced by inflammation, such as by H. pylori infection. Our current study provides new molecular insights into how ESRRG regulates H. pylori infection, contributing to GC development. We suggest that modulation of ESRRG-suppressing H. pylori infection could be a therapeutic target for the treatment of GC patients.

Nonlabens ponticola sp. nov., isolated from seawater and reclassification of Nonlabens sediminis as a later heterotypic synonym of Nonlabens tegetincola

A Gram-stain-negative, orange-pigmented and strictly aerobic bacterium, designated strain MJ115^T, was isolated from seawater in Pohang, South Korea. Cells were non-motile rods and showed positive reactions for catalase and oxidase tests. Growth of strain MJ115^T was observed at 4-35 °C (optimum, 30 °C), pH 6.0-7.0 (optimum, pH 6.5) and in the presence of 0-8.0 % (w/v) NaCl (optimum, 2.0%). Strain MJ115^T contained iso-C_15 : 0, anteiso-C_15 : 0, anteiso-C_17 : 1  ω9c, C_17 : 0 2-OH, iso-C_16 : 0 3-OH, iso-C_17 : 0 3-OH and summed feature 3 (C_16 : 1  ω7c and/or C_16 : 1  ω6c) as major cellular fatty acids and menaquinone-6 as the major respiratory quinone. Phosphatidylethanolamine, two unidentified aminolipids and four unidentified lipids were detected as major polar lipids. The G+C content of the genomic DNA was 40.7 mol%. Phylogenetic analysis based on 16S rRNA gene sequences showed that strain MJ115^T formed a phyletic lineage with Nonlabens marinus S1-08^T, Nonlabens agnitus JC2678^T and Nonlabens antarcticus AKS 622^T within the genus Nonlabens. Strain MJ115^T was also most closely related to N. marinus S1-08^T, N. agnitus JC2678^T and N. antarcticus AKS 622^T with 96.5, 96.4 and 96.0 % 16S rRNA sequence similarities, respectively. Here it is proposed that strain MJ115^T represents a new species of the genus Nonlabens, for which the name Nonlabens ponticola sp. nov. is proposed. The type strain is MJ115^T (=KCTC 72237^T=NBRC 113963^T). In addition, the comparison of the whole genome sequences and phenotypic features suggested that Nonlabens tegetincola and Nonlabens sediminis belong to the same species. Therefore, it is proposed that N. sediminis is reclassified as a later heterotypic synonym of N. tegetincola.

TRIM24-RIP3 axis perturbation accelerates osteoarthritis pathogenesis

OBJECTIVES

Recently, necroptosis has attracted increasing attention in arthritis research; however, it remains unclear whether its regulation is involved in osteoarthritis (OA) pathogenesis. Since receptor-interacting protein kinase-3 (RIP3) plays a pivotal role in necroptosis and its dysregulation is involved in various pathological processes, we investigated the role of the RIP3 axis in OA pathogenesis.

METHODS

Experimental OA was induced in wild-type or Rip3 knockout mice by surgery to destabilise the medial meniscus (DMM) or the intra-articular injection of adenovirus carrying a target gene (Ad-Rip3 and Ad-Trim24 shRNA). RIP3 expression was examined in OA cartilage from human patients; Trim24, a negative regulator of RIP3, was identified by microarray and in silico analysis. Connectivity map (CMap) and in silico binding approaches were used to identify RIP3 inhibitors and to examine their direct regulation of RIP3 activation in OA pathogenesis.

RESULTS

RIP3 expression was markedly higher in damaged cartilage from patients with OA than in undamaged cartilage. In the mouse model, adenoviral RIP3 overexpression accelerated cartilage disruption, whereas Rip3 depletion reduced DMM-induced OA pathogenesis. Additionally, TRIM24 knockdown upregulated RIP3 expression; its downregulation promoted OA pathogenesis in knee joint tissues. The CMap approach and in silico binding assay identified AZ-628 as a potent RIP3 inhibitor and demonstrated that it abolished RIP3-mediated OA pathogenesis by inhibiting RIP3 kinase activity.

CONCLUSIONS

TRIM24-RIP3 axis perturbation promotes OA chronicity by activating RIP3 kinase, suggesting that the therapeutic manipulation of this pathway could provide new avenues for treating OA.

Comparative Evaluation of Genome Assemblers from Long-Read Sequencing for Plants and Crops

The availability of recent state-of-the-art long-read sequencing technologies has significantly increased the ease and speed of producing high-quality plant genome assemblies. A wide variety of genome-related software tools are now available and they are typically benchmarked using microbial or model eukaryotic genomes such as Arabidopsis and rice. However, many plant species have much larger and more complex genomes than these, and the choice of tools, parameters, and/or strategies that can be used is not always obvious. Thus, we have compared the metrics of assemblies generated by various pipelines to discuss how assembly quality can be affected by two different assembly strategies. First, we focused on optimizing read preprocessing and assembler variables using eight different de novo assemblers on five different Pacific Biosciences long-read datasets of diploid and tetraploid species. Then, we examined a single scaffolding tool (quickmerge) that has been employed for the postprocessing step. We then merged the outputs from multiple assemblies to produce a higher quality consensus assembly. Then, we benchmarked the assemblies for completeness and accuracy (assembly metrics and BUSCO), computer memory, and CPU times. Two lightweight assemblers, Miniasm/Minimap/Racon and WTDBG, were deemed good for novice users because they involved smaller required learning curves and light computational resources. However, two heavyweight tools, CANU and Flye, should be the first choice when the goal is to achieve accurate and complete assemblies. Our results will provide valuable guidance in future plant genome projects and beyond.

3'-Sialyllactose prebiotics prevents skin inflammation via regulatory T cell differentiation in atopic dermatitis mouse models

3'-Sialyllactose (3'-SL), a natural prebiotic, maintains immune homeostasis and exerts anti-inflammatory and anti-arthritic effects. Although regulatory T cells (Tregs) prevent excessive inflammation and maintain immune tolerance, the effect of 3'-SL on Treg regulation is unclear. This study aimed to investigate the effect of 3'-SL on Treg responses in atopic dermatitis (AD) pathogenesis. Oral administration of 3'-SL reduced AD-like symptoms such as ear, epidermal, and dermal thickness in repeated topical application of house dust mites (HDM) and 2,4-dinitrochlorobenzene (DNCB). 3'-SL inhibited IgE, IL-1β, IL-6, and TNF-α secretion and markedly downregulated AD-related cytokines including IL-4, IL-5, IL-6, IL-13, IL-17, IFN-γ, TNF-α, and Tslp through regulation of NF-κB in ear tissue. Additionally, in vitro assessment of Treg differentiation revealed that 3'-SL directly induced TGF-β-mediated Treg differentiation. Furthermore, 3'-SL administration also ameliorated sensitization and elicitation of AD pathogenesis by suppressing mast cell infiltration and production of IgE and pro-inflammatory cytokines in mouse serum by mediating the Treg response. Furthermore, Bifidobacterium population was also increased by 3'-SL administration as prebiotics. Our data collectively show that 3'-SL has therapeutic effects against AD progression by inducing Treg differentiation, downregulating AD-related cytokines, and increasing the Bifidobacterium population.

Upregulation of cellulase activity and mRNA levels by bacterial challenge in the earthworm Eisenia andrei, supporting the involvement of cellulases in innate immunity

To investigate whether earthworm cellulases contribute to the innate immune system, the responsiveness of cellulase activity and mRNA expression to bacterial challenge was examined by zymography and RNA sequencing. A zymographic analysis revealed that the activity levels of earthworm cellulases were upregulated in response to either a bacterial (Bacillus subtilis or Escherichia coli) or LPS challenge. After the challenge, significant increases in cellulase 1 and cellulase 2 activity levels were observed within 8-16 and 16-24 h, respectively. In the coelomic fluid, both activities were significantly upregulated at 8 h post-injection with B. subtilis. Based on RNA sequencing, cellulase-related mRNAs encoding beta-1,4-endoglucanases were upregulated by 3-fold within 6 h after B. subtilis injection. Our results clearly demonstrated that earthworm cellulases are upregulated by bacterial challenge at the mRNA and protein levels. These results support the view that earthworm cellulases act as inducible humoral effectors of innate immunity against bacterial infection.

Cirsium japonicum var. maackii and apigenin block Hif-2α-induced osteoarthritic cartilage destruction

Although Hif-2α is a master regulator of catabolic factor expression in osteoarthritis development, Hif-2α inhibitors remain undeveloped. The aim of this study was to determine whether Cirsium japonicum var. maackii (CJM) extract and one of its constituents, apigenin, could attenuate the Hif-2α-induced cartilage destruction implicated in osteoarthritis progression. In vitro and in vivo studies demonstrated that CJM reduced the IL-1β-, IL-6, IL-17- and TNF-α-induced up-regulation of MMP3, MMP13, ADAMTS4, ADAMTS5 and COX-2 and blocked osteoarthritis development in a destabilization of the medial meniscus mouse model. Activation of Hif-2α, which directly up-regulates MMP3, MMP13, ADAMTS4, IL-6 and COX-2 expression, is inhibited by CJM extract. Although cirsimarin, cirsimaritin and apigenin are components of CJM and can reduce inflammation, only apigenin effectively reduced Hif-2α expression and inhibited Hif-2α-induced MMP3, MMP13, ADAMTS4, IL-6 and COX-2 expression in articular chondrocytes. IL-1β induction of JNK phosphorylation and IκB degradation, representing a critical pathway for Hif-2α expression, was completely blocked by apigenin in a concentration-dependent manner. Collectively, these effects indicate that CJM and one of its most potent constituents, apigenin, can lead to the development of therapeutic agents for blocking osteoarthritis development as novel Hif-2α inhibitors.

Accelerated pseudogenization of trace amine-associated receptor genes in primates

Trace amines (TAs) in the mammalian brain have been investigated for four decades. Trace amine‐associated receptors (TAARs) were discovered during the search for receptors activated by TAs. TAARs are considered a second class of vertebrate olfactory receptors and successfully proliferated in conjunction with adaptation to living on the ground to detect carnivore odors. Thus, therian mammals have a high number of TAAR genes due to rapid species‐specific gene duplications. In primate lineages, however, their genomes have significantly smaller numbers of TAAR genes than do other mammals. To elucidate the evolutionary force driving these patterns, exhaustive data mining of TAAR genes was performed for 13 primate genomes (covering all four infraorders) and two nonprimate euarchontan genomes. This study identified a large number of pseudogenes in many of these primate genomes and thus investigated the pseudogenization event process for the TAAR repertoires. The degeneration of TAARs is likely associated with arboreal inhabitants reducing their exposure to carnivores, and this was accelerated by the change in the nose shape of haplorhines after their divergence from strepsirrhines. Arboreal life may have decreased the reliance on the chemosensing of predators, suggestive of leading to the depauperation of TAAR subfamilies. The evolutionary deterioration of TAARs in primates has been reestablished in recently derived primates due to high selection pressure and probably functional diversity.

Patterns of Gene Expression in Western Corn Rootworm (Diabrotica virgifera virgifera) Neonates, Challenged with Cry34Ab1, Cry35Ab1 and Cry34/35Ab1, Based on Next-Generation Sequencing

With Next Generation Sequencing technologies, high-throughput RNA sequencing (RNAseq) was conducted to examine gene expression in neonates of Diabrotica virgifera virgifera (LeConte) (Western Corn Rootworm, WCR) challenged with individual proteins of the binary Bacillus thuringiensis insecticidal proteins, Cry34Ab1 and Cry35Ab1, and the combination of Cry34/Cry35Ab1, which together are active against rootworm larvae. Integrated results of three different statistical comparisons identified 114 and 1300 differentially expressed transcripts (DETs) in the Cry34Ab1 and Cry34/35Ab1 treatment, respectively, as compared to the control. No DETs were identified in the Cry35Ab1 treatment. Putative Bt binding receptors previously identified in other insect species were not identified in DETs in this study. The majority of DETs (75% with Cry34Ab1 and 68.3% with Cry34/35Ab1 treatments) had no significant hits in the NCBI nr database. In addition, 92 DETs were shared between Cry34Ab1 and Cry34/35Ab1 treatments. Further analysis revealed that the most abundant DETs in both Cry34Ab1 and Cry34/35Ab1 treatments were associated with binding and catalytic activity. Results from this study confirmed the nature of these binary toxins against WCR larvae and provide a fundamental profile of expression pattern of genes in response to challenge of the Cry34/35Ab1 toxin, which may provide insight into potential resistance mechanisms.

Phylogenomic analysis of Copepoda (Arthropoda, Crustacea) reveals unexpected similarities with earlier proposed morphological phylogenies

BACKGROUND

Copepods play a critical role in marine ecosystems but have been poorly investigated in phylogenetic studies. Morphological evidence supports the monophyly of copepods, whereas interordinal relationships continue to be debated. In particular, the phylogenetic position of the order Harpacticoida is still ambiguous and inconsistent among studies. Until now, a small number of molecular studies have been done using only a limited number or even partial genes and thus there is so far no consensus at the order-level.

RESULTS

This study attempted to resolve phylogenetic relationships among and within four major copepod orders including Harpacticoida and the phylogenetic position of Copepoda among five other crustacean groups (Anostraca, Cladocera, Sessilia, Amphipoda, and Decapoda) using 24 nuclear protein-coding genes. Phylogenomics has confirmed the monophyly of Copepoda and Podoplea. However, this study reveals surprising differences with the majority of the copepod phylogenies and unexpected similarities with postembryonic characters and earlier proposed morphological phylogenies; More precisely, Cyclopoida is more closely related to Siphonostomatoida than to Harpacticoida which is likely the most basally-branching group of Podoplea. Divergence time estimation suggests that the origin of Harpacticoida can be traced back to the Devonian, corresponding well with recently discovered fossil evidence. Copepoda has a close affinity to the clade of Malacostraca and Thecostraca but not to Branchiopoda. This result supports the hypothesis of the newly proposed clades, Communostraca, Multicrustacea, and Allotriocarida but further challenges the validity of Hexanauplia and Vericrustacea.

CONCLUSIONS

The first phylogenomic study of Copepoda provides new insights into taxonomic relationships and represents a valuable resource that improves our understanding of copepod evolution and their wide range of ecological adaptations.

Transcriptional analysis of susceptible and resistant European corn borer strains and their response to Cry1F protoxin

Background

Despite a number of recent reports of insect resistance to transgenic crops expressing insecticidal toxins from Bacillus thuringiensis (Bt), little is known about the mechanism of resistance to these toxins. The purpose of this study is to identify genes associated with the mechanism of Cry1F toxin resistance in European corn borer (Ostrinia nubilalis Hübner). For this, we compared the global transcriptomic response of laboratory selected resistant and susceptible O. nubilalis strain to Cry1F toxin. We further identified constitutive transcriptional differences between the two strains.

Results

An O. nubilalis midgut transcriptome of 36,125 transcripts was assembled de novo from 106 million Illumina HiSeq and Roche 454 reads and used as a reference for estimation of differential gene expression analysis. Evaluation of gene expression profiles of midgut tissues from the Cry1F susceptible and resistant strains after toxin exposure identified a suite of genes that responded to the toxin in the susceptible strain (n = 1,654), but almost 20-fold fewer in the resistant strain (n = 84). A total of 5,455 midgut transcripts showed significant constitutive expression differences between Cry1F susceptible and resistant strains. Transcripts coding for previously identified Cry toxin receptors, cadherin and alkaline phosphatase and proteases were also differentially expressed in the midgut of the susceptible and resistant strains.

Conclusions

Our current study provides a valuable resource for further molecular characterization of Bt resistance and insect response to Cry1F toxin in O. nubilalis and other pest species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1751-6) contains supplementary material, which is available to authorized users.

Genetic identification and molecular phylogeny of Pseudodiaptomus species (Calanoida, Pseudodiaptomidae) in Korean waters

Five species of the genus Pseudodiaptomus, P. inopinus, P. poplesia, P. marinus, P. nihonkaiensis, and P. sp. occur in Korea nearshore waters. Although the four species except for P. sp. have been classified into Lobus and Ramosus groups, two for each group, based on morphological characters, this classification had yet to be confirmed by molecular characters. Here, we determined molecular characters and phylogenetic relationships of the five species in order to evaluate the morphology-based groupings and the species identifications. For this, a 625-bp DNA region of the mitochondrial gene cytochrome oxidase subunit I (mtCOI) was sequenced and compared among the species. Intraspecific variation of the sequences is less than 0.6%, while interspecific variation ranges from 17.6-26.7%, indicating every species, including P. sp., is a genetically distinct, valid species. Phylogenetic trees of the mtCOI DNA reveal that the Lobus-group species including P. inopinus and P. poplesia form a well-supported clade and that P. sp. belongs to this group. On the other hand, the Ramosus group clade consisting of P. marinus and P. nihonkaiensis is not well supported by bootstrap analyses, suggesting that further evaluation of the validity of this group assignment is needed.