COWID: an efficient cloud-based genomics workflow for scalable identification of SARS-COV-2.

Implementing a specific cloud resource to analyze extensive genomic data on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) poses a challenge when resources are limited. To overcome this, we repurposed a cloud platform initially designed for use in research on cancer genomics (https://cgc.sbgenomics.com) to enable its use in research on SARS-CoV-2 to build Cloud Workflow for Viral and Variant Identification (COWID). COWID is a workflow based on the Common Workflow Language that realizes the full potential of sequencing technology for use in reliable SARS-CoV-2 identification and leverages cloud computing to achieve efficient parallelization. COWID outperformed other contemporary methods for identification by offering scalable identification and reliable variant findings with no false-positive results. COWID typically processed each sample of raw sequencing data within 5 min at a cost of only US$0.01. The COWID source code is publicly available (https://github.com/hendrick0403/COWID) and can be accessed on any computer with Internet access. COWID is designed to be user-friendly; it can be implemented without prior programming knowledge. Therefore, COWID is a time-efficient tool that can be used during a pandemic.

De novo assembly and annotation of the Patagonian toothfish (Dissostichus eleginoides) genome.

BACKGROUND: Patagonian toothfish (Dissostichus eleginoides) is an economically and ecologically important fish species in the family Nototheniidae. Juveniles occupy progressively deeper waters as they mature and grow, and adults have been caught as deep as 2500 m, living on or in just above the southern shelves and slopes around the sub-Antarctic islands of the Southern Ocean. As apex predators, they are a key part of the food web, feeding on a variety of prey, including krill, squid, and other fish. Despite its importance, genomic sequence data, which could be used for more accurate dating of the divergence between Patagonian and Antarctic toothfish, or establish whether it shares adaptations to temperature with fish living in more polar or equatorial climes, has so far been limited. RESULTS: A high-quality D. eleginoides genome was generated using a combination of Illumina, PacBio and Omni-C sequencing technologies. To aid the genome annotation, the transcriptome derived from a variety of toothfish tissues was also generated using both short and long read sequencing methods. The final genome assembly was 797.8 Mb with a N50 scaffold length of 3.5 Mb. Approximately 31.7% of the genome consisted of repetitive elements. A total of 35,543 putative protein-coding regions were identified, of which 50% have been functionally annotated. Transcriptomics analysis showed that approximately 64% of the predicted genes (22,617 genes) were found to be expressed in the tissues sampled. Comparative genomics analysis revealed that the anti-freeze glycoprotein (AFGP) locus of D. eleginoides does not contain any AFGP proteins compared to the same locus in the Antarctic toothfish (Dissostichus mawsoni). This is in agreement with previously published results looking at hybridization signals and confirms that Patagonian toothfish do not possess AFGP coding sequences in their genome. CONCLUSIONS: We have assembled and annotated the Patagonian toothfish genome, which will provide a valuable genetic resource for ecological and evolutionary studies on this and other closely related species.

NSPlex: an efficient method to analyze non-specific peaks amplified using commercial STR kits.

Commercial short tandem repeat (STR) kits exclusively contain human-specific primers; however, various non-human organisms with high homology to the STR kit's primer sequences can cause cross-reactivity. Owing to the proprietary nature of the primers in STR kits, the origins and sequences of most non-specific peaks (NSPs) remain unclear. Such NSPs can complicate data interpretation between the casework and reference samples; thus, we developed "NSPlex", an efficient method to discover the biological origins of NSPs. We used leftover STR kit amplicons after capillary electrophoresis and performed advanced bioinformatics analyses using next-generation sequencing followed by BLAST nucleotide searches. Using our method, we could successfully identify NSP generated from PCR amplicons of a sample mixture of human DNA and DNA extracted from matcha powder (finely ground powder of green tea leaves and previously known as a potential source of NSP). Our results showed our method is efficient for NSP analysis without the need for the primer information as in commercial STR kits.

The Effect of Enzymatic Treatment with Mutanase, Beta-Glucanase, and DNase on a Saliva-Derived Biofilm Model.

INTRODUCTION: The dental biofilm matrix is an important determinant of virulence for caries development and comprises a variety of extracellular polymeric substances that contribute to biofilm stability. Enzymes that break down matrix components may be a promising approach to caries control, and in light of the compositional complexity of the dental biofilm matrix, treatment with multiple enzymes may enhance the reduction of biofilm formation compared to single enzyme therapy. The present study investigated the effect of the three matrix-degrading enzymes mutanase, beta-glucanase, and DNase, applied separately or in combinations, on biofilm prevention and removal in a saliva-derived in vitro-grown model. METHODS: Biofilms were treated during growth to assess biofilm prevention or after 24 h of growth to assess biofilm removal by the enzymes. Biofilms were quantified by crystal violet staining and impedance-based real-time cell analysis, and the biofilm structure was visualized by confocal microscopy and staining of extracellular DNA (eDNA) and polysaccharides. RESULTS: The in vitro model was dominated by Streptococcus spp., as determined by 16S rRNA gene amplicon sequencing. All tested enzymes and combinations had a significant effect on biofilm prevention, with reductions of >90% for mutanase and all combinations including mutanase. Combined application of DNase and beta-glucanase resulted in an additive effect (81.0% ± 1.3% SD vs. 36.9% ± 21.9% SD and 48.2% ± 14.9% SD). For biofilm removal, significant reductions of up to 73.2% ± 5.5% SD were achieved for combinations including mutanase, whereas treatment with DNase had no effect. Glucans, but not eDNA decreased in abundance upon treatment with all three enzymes. CONCLUSION: Multi-enzyme treatment is a promising approach to dental biofilm control that needs to be validated in more diverse biofilms.

An updated LSU database and pipeline for environmental DNA identification of arbuscular mycorrhizal fungi.

Recent work established a backbone reference tree and phylogenetic placement pipeline for identification of arbuscular mycorrhizal fungal (AMF) large subunit (LSU) rDNA environmental sequences. Our previously published pipeline allowed any environmental sequence to be identified as putative AMF or within one of the major families. Despite this contribution, difficulties in implementation of the pipeline remain. Here, we present an updated database and pipeline with (1) an expanded backbone tree to include four newly described genera and (2) several changes to improve ease and consistency of implementation. In particular, packages required for the pipeline are now installed as a single folder (conda environment) and the pipeline has been tested across three university computing clusters. This updated backbone tree and pipeline will enable broadened adoption by the community, advancing our understanding of these ubiquitous and ecologically important fungi.

Identification of plasmids in avian-associated Escherichia coli using nanopore and illumina sequencing.

BACKGROUND: Avian pathogenic Escherichia coli (APEC) are the causative agents of colibacillosis in chickens, a disease which has significant economic impact on the poultry industry. Large plasmids detected in APEC are known to contribute to strain diversity for pathogenicity and antimicrobial resistance, but there could be other plasmids that are missed in standard analysis. In this study, we determined the impact of sequencing and assembly factors for the detection of plasmids in an E. coli whole genome sequencing project. RESULTS: Hybrid assembly (Illumina and Nanopore) combined with plasmid DNA extractions allowed for detection of the greatest number of plasmids in E. coli, as detected by MOB-suite software. In total, 79 plasmids were identified in 19 E. coli isolates. Hybrid assemblies were robust and consistent in quality regardless of sequencing kit used or if long reads were filtered or not. In contrast, long read only assemblies were more variable and influenced by sequencing and assembly parameters. Plasmid DNA extractions allowed for the detection of physically smaller plasmids, but when averaged over 19 isolates did not significantly change the overall number of plasmids detected. CONCLUSIONS: Hybrid assembly can be reliably used to detect plasmids in E. coli, especially if researchers are focused on large plasmids containing antimicrobial resistance genes and virulence factors. If the goal is comprehensive detection of all plasmids, particularly if smaller sized vectors are desired for biotechnology applications, the addition of plasmid DNA extractions to hybrid assemblies is prudent. Long read sequencing is sufficient to detect many plasmids in E. coli, however, it is more prone to errors when expanded to analyze a large number of isolates.

Transcriptome analysis of the endangered dung beetle Copris tripartitus (Coleoptera: Scarabaeidae) and characterization of genes associated to immunity, growth, and reproduction.

BACKGROUND: Dung beetles recycle organic matter through the decomposition of feces and support ecological balance. However, these insects are threatened by the indiscriminate use of agrochemicals and habitat destruction. Copris tripartitus Waterhouse (Coleoptera: Scarabaeidae), a dung beetle, is listed as a class-II Korean endangered species. Although the genetic diversity of C. tripartitus populations has been investigated through analysis of mitochondrial genes, genomic resources for this species remain limited. In this study, we analyzed the transcriptome of C. tripartitus to elucidate functions related to growth, immunity and reproduction for the purpose of informed conservation planning. RESULTS: The transcriptome of C. tripartitus was generated using next-generation Illumina sequencing and assembled de novo using a Trinity-based platform. In total, 98.59% of the raw sequence reads were processed as clean reads. These reads were assembled into 151,177 contigs, 101,352 transcripts, and 25,106 unigenes. A total of 23,450 unigenes (93.40%) were annotated to at least one database. The largest proportion of unigenes (92.76%) were annotated to the locally curated PANM-DB. A maximum of 5,512 unigenes had homologous sequences in Tribolium castaneum. Gene Ontology (GO) analysis revealed a maximum of 5,174 unigenes in the Molecular function category. Further, in Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, a total of 462 enzymes were associated with established biological pathways. Based on sequence homology to known proteins in PANM-DB, representative immunity, growth, and reproduction-related genes were screened. Potential immunity-related genes were categorized into pattern recognition receptors (PRRs), the Toll-like receptor signaling pathway, the MyD88- dependent pathway, endogenous ligands, immune effectors, antimicrobial peptides, apoptosis, and adaptation-related transcripts. Among PRRs, we conducted detailed in silico characterization of TLR-2, CTL, and PGRP_SC2-like. Repetitive elements such as long terminal repeats, short interspersed nuclear elements, long interspersed nuclear elements and DNA elements were enriched in the unigene sequences. A total of 1,493 SSRs were identified among all unigenes of C. tripartitus. CONCLUSIONS: This study provides a comprehensive resource for analysis of the genomic topography of the beetle C. tripartitus. The data presented here clarify the fitness phenotypes of this species in the wild and provide insight to support informed conservation planning.

Exploring microbial worlds: a review of whole genome sequencing and its application in characterizing the microbial communities.

The classical microbiology techniques have inherent limitations in unraveling the complexity of microbial communities, necessitating the pivotal role of sequencing in studying the diversity of microbial communities. Whole genome sequencing (WGS) enables researchers to uncover the metabolic capabilities of the microbial community, providing valuable insights into the microbiome. Herein, we present an overview of the rapid advancements achieved thus far in the use of WGS in microbiome research. There was an upsurge in publications, particularly in 2021 and 2022 with the United States, China, and India leading the metagenomics research landscape. The Illumina platform has emerged as the widely adopted sequencing technology, whereas a significant focus of metagenomics has been on understanding the relationship between the gut microbiome and human health where distinct bacterial species have been linked to various diseases. Additionally, studies have explored the impact of human activities on microbial communities, including the potential spread of pathogenic bacteria and antimicrobial resistance genes in different ecosystems. Furthermore, WGS is used in investigating the microbiome of various animal species and plant tissues such as the rhizosphere microbiome. Overall, this review reflects the importance of WGS in metagenomics studies and underscores its remarkable power in illuminating the variety and intricacy of the microbiome in different environments.

Root exudates and rhizosphere microbiomes jointly determine temporal shifts in plant-soil feedbacks.

Plants influence numerous soil biotic factors that can alter the performance of later growing plants-defined as plant-soil feedback (PSF). Here, we investigate whether PSF effects are linked with the temporal changes in root exudate diversity and the rhizosphere microbiome of two common grassland species (Holcus lanatus and Jacobaea vulgaris). Both plant species were grown separately establishing conspecific and heterospecific soils. In the feedback phase, we determined plant biomass, measured root exudate composition, and characterised rhizosphere microbial communities weekly (eight time points). Over time, we found a strong negative conspecific PSF on J. vulgaris in its early growth phase which changed into a neutral PSF, whereas H. lanatus exhibited a more persistent negative PSF. Root exudate diversity increased considerably over time for both plant species. Rhizosphere microbial communities were distinct in conspecific and heterospecific soils and showed strong temporal patterns. Bacterial communities converged over time. Using path models, PSF effects could be linked to the temporal dynamics of root exudate diversity, whereby shifts in rhizosphere microbial diversity contributed to temporal variation in PSF to a lesser extent. Our results highlight the importance of root exudates and rhizosphere microbial communities in driving temporal changes in the strength of PSF effects.

Assessment of the dibenzothiophene desulfurization potential of indigenously isolated bacterial consortium IQMJ-5: a different approach to safeguard the environment.

Biodesulfurization is emerging as a valuable technology for the desulfurization of dibenzothiophene (DBT) and its alkylated substitutes, which are otherwise regarded as refractory to other physical and chemical desulfurizing techniques. The inability of the currently identified pure cultures and artificial microbial consortia due to lower desulfurization rate and product inhibition issues has compelled the researcher to look for an alternative solution. Thus, in the present study, an indigenously isolated microbial consortium was employed to tackle the desulfurization issue. Herein, we isolated several kinds of DBT desulfurizing natural microbial consortia from hydrocarbon-contaminated soil samples by conventional enrichment technique. The most effective desulfurizing microbial consortium was sequenced through illumine sequencing technique. Finally, the effect of the products of the desulfurizing pathway (such as 2-hydroxybiphenyl (2-HBP) and sulfate (SO4-2) was evaluated on the growth and desulfurization capability of the isolated consortium. The outcomes of Gibb's assay analysis showed that six isolates followed the "4S" pathway and converted DBT to 2-HBP. Among the isolates, I5 showed maximum growth rate (1.078 g/L dry cell weight) and desulfurization activity (about 77% as indicated by HPLC analysis) and was considered for further in-depth experimentation. The analysis of 16S rRNA by high-throughput sequencing approach of the I5 isolate revealed five types of bacterial phyla including Proteobacteria, Bacteroidetes, Firmicutes, Patescibacteria, and Actinobacteria (in order of abundance). The isolate showed significant tolerance to the inhibitory effect of both 2-HBP and SO4-2 and maintained growth in the presence of even about 1.0 mM initial concentration of both products. This clearly suggests that the isolate can be an efficient candidate for future in-depth desulfurization studies of coal and other fossil fuels.

Genomic profiling and characteristics of a C1 degrading heterotrophic fresh-water bacterium Paracoccus sp. strain DMF.

Paracoccus species are metabolically versatile gram-negative, aerobic facultative methylotrophic bacteria showing enormous promise for environmental and bioremediation studies. Here we report, the complete genome analysis of Paracoccus sp. strain DMF (P. DMF) that was isolated from a domestic wastewater treatment plant in Kanpur, India (26.4287 °N, 80.3891 °E) based on its ability to degrade a recalcitrant organic solvent N, N-dimethylformamide (DMF). The results reveal a genome size of 4,202,269 base pairs (bp) with a G + C content of 67.9%. The assembled genome comprises 4141 coding sequences (CDS), 46 RNA sequences, and 2 CRISPRs. Interestingly, catabolic operons related to the conventional marine-based methylated amines (MAs) degradation pathway were functionally annotated within the genome of an obligated aerobic heterotroph that is P. DMF. The genomic data-based characterization presented here for the novel heterotroph P. DMF aims to improve the understanding of the phenotypic gene products, enzymes, and pathways involved with greater emphasis on facultative methylotrophic motility-based latent pathogenicity.

Unmapped short reads from whole-genome sequencing indicate potential infectious pathogens in german black Pied cattle.

When resequencing animal genomes, some short reads cannot be mapped to the reference genome and are usually discarded. In this study, unmapped reads from 302 German Black Pied cattle were analyzed to identify potential pathogenic DNA. These unmapped reads were assembled and blasted against NCBI's database to identify bacterial and viral sequences. The results provided evidence for the presence of pathogens. We found sequences of Bovine parvovirus 3 and Mycoplasma species. These findings emphasize the information content of unmapped reads for gaining insight into bacterial and viral infections, which is important for veterinarians and epidemiologists.

The Variations of Microbial Diversity and Community Structure Along Different Stream Orders in Wuyi Mountains.

The surface water is an important habitat for freshwater microorganisms, but there is a lack of understanding of the pattern of microbial diversity and structure in stream continuums of small subtropical forest watersheds. Therefore, this study aimed to understand the variations in microbial diversity and community structure along stream orders (1-5) in the small subtropical forest catchments of the Wuyi Mountains. Using GIS software, 20 streams were chosen and classified into 5 orders. Illumina sequencing was used to analyze the dynamics of microbial communities, along with stream orders and hydro-chemical properties of stream water were also determined. Our results indicated that the bacterial and fungal richness (ACE index) was higher in low-order (1 and 2 orders) streams than in high-order (3, 4, and 5 orders) streams, with the highest value in the order 2 streams (P < 0.05). The water temperature and dissolved oxygen were positively correlated with fungal richness (P < 0.05). The bacterial rare taxa had a significant correlation with the abundance taxa (P < 0.05). The relative abundances of Bacteroidetes, Actinobacteria, and Chytridiomycota microbial phyla were significantly different among different order streams (P < 0.05). Using the neutral community model, we found that the fungal community structure was significantly shaped by hydro-chemical properties, while the bacterial community structure was largely regulated by stochastic processes. Our findings suggest that variations in microbial community structure in subtropical headwaters are largely shaped by the water temperature and dissolved oxygen.

Extensive ITR expansion of the 2022 Mpox virus genome through gene duplication and gene loss.

Poxviruses are known to evolve slower than RNA viruses with only 1-2 mutations/genome/year. Rather than single mutations, rearrangements such as gene gain and loss, which have been discussed as a possible driver for host adaption, were described in poxviruses. In 2022 and 2023 the world is being challenged by the largest global outbreak so far of Mpox virus, and the virus seems to have established itself in the human community for an extended period of time. Here, we report five Mpox virus genomes from Germany with extensive gene duplication and loss, leading to the expansion of the ITR regions from 6400 to up to 24,600 bp. We describe duplications of up to 18,200 bp to the opposed genome end, and deletions at the site of insertion of up to 16,900 bp. Deletions and duplications of genes with functions of supposed immune modulation, virulence and host adaption as B19R, B21R, B22R and D10L are described. In summary, we highlight the need for monitoring rearrangements of the Mpox virus genome rather than for monitoring single mutations only.

Benchmarking of two bioinformatic workflows for the analysis of whole-genome sequenced Staphylococcus aureus collected from patients with suspected sepsis.

BACKGROUND: The rapidly growing area of sequencing technologies, and more specifically bacterial whole-genome sequencing, could offer applications in clinical microbiology, including species identification of bacteria, prediction of genetic antibiotic susceptibility and virulence genes simultaneously. To accomplish the aforementioned points, the commercial cloud-based platform, 1928 platform (1928 Diagnostics, Gothenburg, Sweden) was benchmarked against an in-house developed bioinformatic pipeline as well as to reference methods in the clinical laboratory. METHODS: Whole-genome sequencing data retrieved from 264 Staphylococcus aureus isolates using the Illumina HiSeq X next-generation sequencing technology was used. The S. aureus isolates were collected during a prospective observational study of community-onset severe sepsis and septic shock in adults at Skaraborg Hospital, in the western region of Sweden. The collected isolates were characterized according to accredited laboratory methods i.e., species identification by MALDI-TOF MS analysis and phenotypic antibiotic susceptibility testing (AST) by following the EUCAST guidelines. Concordance between laboratory methods and bioinformatic tools, as well as concordance between the bioinformatic tools was assessed by calculating the percent of agreement. RESULTS: There was an overall high agreement between predicted genotypic AST and phenotypic AST results, 98.0% (989/1006, 95% CI 97.3-99.0). Nevertheless, the 1928 platform delivered predicted genotypic AST results with lower very major error rates but somewhat higher major error rates compared to the in-house pipeline. There were differences in processing times i.e., minutes versus hours, where the 1928 platform delivered the results faster. Furthermore, the bioinformatic workflows showed overall 99.4% (1267/1275, 95% CI 98.7-99.7) agreement in genetic prediction of the virulence gene characteristics and overall 97.9% (231/236, 95% CI 95.0-99.2%) agreement in predicting the sequence types (ST) of the S. aureus isolates. CONCLUSIONS: Altogether, the benchmarking disclosed that both bioinformatic workflows are able to deliver results with high accuracy aiding diagnostics of severe infections caused by S. aureus. It also illustrates the need of international agreement on quality control and metrics to facilitate standardization of analytical approaches for whole-genome sequencing based predictions.

Domestication shapes the endophytic microbiome and metabolome of Salicornia europaea.

AIMS: We aim at understanding the effect of domestication on the endophytic microbiome and metabolome of Salicornia europaea and collecting evidence on the potential role of microbial populations and metabolites in the adaptation of plants to different ecological contexts (wild vs crops). METHODS AND RESULTS: Samples were collected from a natural salt marsh (wild) and an intensive crop field (crop). High-throughput sequencing of the 16S rRNA gene, gas chromatography-mass spectrometry (GC-MS) and ultra-performance liquid chromatography-mass spectrometry (UPLC-MS) were used to analyze the endophytic bacterial communities and the metabolite profiles of S. europaea roots, respectively. The elemental analysis of the plant shoots was performed by Inductively Coupled Plasma-Mass Spectroscopy (ICP-MS).Overall, significant differences were found between the microbiome of wild and cultivated plants. The later showed a higher relative abundance of the genera Erythrobacter, Rhodomicrobium, and Ilumatobacter than wild plants. The microbiome of wild plants was enriched in Marinobacter, Marixanthomonas, and Thalassospira. The metabolite profile of crop plants revealed higher amounts of saturated and non-saturated fatty acids and acylglycerols. In contrast, wild plants contained comparatively more carbohydrates and most macroelements (i.e. Na, K, Mg, and Ca). CONCLUSIONS: There is a strong correlation between plant metabolites and the endosphere microbiome of S. europaea. In wild populations, plants were enriched in carbohydrates and the associated bacterial community was enriched in genes related to primary metabolic pathways such as nitrogen metabolism and carbon fixation. The endosphere microbiome of crop plants was predicted to have higher gene counts related to pathogenesis. Crop plants also exhibited higher amounts of azelaic acid, an indicator of exposure to phytopathogens.

Relic DNA effects on the estimates of bacterial community composition and taxa dynamics in soil.

DNA-based analyses have become routine methods in soil microbial research, for their high throughput and resolution in characterizing microbial communities. Yet, concerns arise regarding the interference of relic DNA in estimates of viable bacterial community composition and individual taxa dynamics in soils that recovered from post-gamma irradiation. In this study, different soil samples with varying bacterial diversity but similar soil properties were randomly selected. We split each sample into two parts: one part was treated with propidium monoazide (PMA) before DNA extraction, PMA can bind to relic DNA and inhibit PCR amplification by chemical modification; DNA of the other part was extracted following the same process but without PMA pretreatment. Then, soil bacterial abundance was quantified by quantitative polymerase chain reaction, and bacterial community structure was examined by Illumina metabarcoding sequencing of 16S rRNA gene. The results showed that the higher bacterial richness and evenness were estimated when relic DNA was present. The variation trends of bacterial abundance, alpha diversity, and beta diversity remained the same, as reflected by the significant correlations between PMA-treated and -untreated samples (P < 0.05). Moreover, as the mean abundance increased, the reproducibility of detecting individual taxa dynamics between relic DNA present and absent treatments increased. These findings provide empirical evidence that a more even distribution of species abundance derived from relic DNA would result in the overestimation of richness in the total DNA pools and also have crucial implications for guiding proper application of high-throughput sequencing to estimate bacterial community diversity and taxonomic population dynamic. KEY POINTS: • Relic DNA effects on the bacterial community in sterilized soils were assessed. • More even species abundance distribution in relic DNA overestimates true richness. • The reproducibility of individual taxa dynamics increased with their abundance.

Analysis on the salt tolerance of Nitraria sibirica Pall. based on Pacbio full-length transcriptome sequencing.

KEY MESSAGE: Nitraria sibirica Pall. regulates its tolerance to salt stress mainly by adjusting ion balance, modifying cell wall structure, and activating signal transduction pathways. N. sibirica, as a typical halophyte, can not only effectively restore saline-alkali land, but also has high economic value. However, studies on its salt tolerance at combining molecular and physiological levels were limited. In this study, the salt tolerance of N. sibirica was analyzed based on Pacbio full-length transcriptome sequencing, and the salt tolerance in the physiological level was verified by key genes. The results showed that 89,017 full-length transcripts were obtained, of which 84,632 sequences were annotated. A total of 86,482 coding sequences (CDS) were predicted and 6561 differentially expressed genes (DEGs) were identified. DEGs were significantly enriched in "sodium ion homeostasis", "response to osmotic stress", "reactive oxygen species metabolic process", "defense response by cell wall thickening", "signal transduction", etc. The expression levels for most of these DEGs increased under salt stress. A total of 69 key genes were screened based on weighted gene co-expression network analysis (WGCNA), of which 33 were first reported on salt tolerance. Moreover, NsRabE1c gene with the highest expression level was selected to verify its salt tolerance. Over-expression of NsRabE1c gene enhanced the germination potential and root length of transgenic Arabidopsis thaliana plants without salt treatment as compared to those of Col-0 and AtRabE1c mutant. The expression levels of NsRabE1c decreased in the growth stagnation phase, while significantly increased in the growth recovery phase under salt stress. We predicted that NsRabE1c gene help N. sibirica resist salt stress through the regulation of plant growth. The results of this study deepen the understanding of salinity resistance in N. sibirica.

Changes in fungal community during different phases in conventional and bioreactor composting systems according to metatranscriptomics analysis.

We determined the changes that occurred in fungal community structures and their functions in conventional and bioreactor composting systems. The Illumina MiSeq platform was employed to sequence cDNA by reverse transcription to conduct metatranscriptomics analysis of RNA, and the FUNGuild tool was applied. The α-diversity of fungi in the bioreactor composter increased throughout composting, especially in the initial three phases, but decreased in the conventional composting system. The three dominant phyla in the bioreactor system were Ascomycota (30.27%-68.50%), Mortierellomycota (3.81%-39.51%), and Basidiomycota (9.17%-30.86%). Ascomycota (76.96%-97.18%) was the main phylum in the conventional composting system. Mortierella, Guehomyces, Plectosphaerella, Chaetomium, Millerozyma, and Coprinopsis were the main genera in the bioreactor composter. In the same phase, significant differences in the fungal functions were found between the two composting methods. Available phosphorus was the main factor that affected the community structures and functions of fungi in the bioreactor composter.

Exploring miRNA-mRNA regulatory modules responding to tannic acid stress in Micromelalopha troglodyta (Graeser) (Lepidoptera: Notodontidae) via small RNA sequencing.

MicroRNAs (miRNAs) are small noncoding RNAs (sRNAs) that regulate gene expression by inhibiting translation or degrading mRNA. Although the functions of miRNAs in many biological processes have been reported, there is currently no research on the possible roles of miRNAs in Micromelalopha troglodyta (Graeser) involved in the response of plant allelochemicals. In this article, six sRNA libraries (three treated with tanic acid and three control) from M. troglodyta were constructed using Illumina sequencing. From the results, 312 known and 43 novel miRNAs were differentially expressed. Notably, some of the most abundant miRNAs, such as miR-432, miR-541-3p, and miR-4448, involved in important physiological processes were also identified. To better understand the function of the targeted genes, we performed Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The results indicated that differentially expressed miRNA targets were involved in metabolism, development, hormone biosynthesis, and immunity. Finally, we visualized a miRNA-mRNA regulatory module that supports the role of miRNAs in host-allelochemical interactions. To our knowledge, this is the first report on miRNAs responding to tannic acid in M. troglodyta. This study provides indispensable information for understanding the potential roles of miRNAs in M. troglodyta and the applications of these miRNAs in M. troglodyta management.