This study demonstrated the strengths of invivo molecular staining coupled with automated imaging analysis in Daphnia magna. A multiwell plate protocol was developed to assess mitochondrial membrane potential using the JC-1 dye. The suitability of five common anesthetics was initially tested, and 5% ethanol performed best in terms of anesthetic effects and healthy recovery. The staining conditions were optimized to 30min staining with 2 μM JC-1 for best J-aggregate formation. The protocol was validated with the model compound carbonyl cyanide 3-chlorophenylhydrazone (CCCP) and used to measure the effect of four environmental contaminants, 2,4-dinitrophenol, triclosan, n-(1,3-dimethylbutyl)-N′-phenyl-p-phenylenediamine (6PPD), and ibuprofen, on mitochondrial health. Test organisms were imaged using anautomated confocal microscope, and fluorescence intensities were automatically quantified. The effect concentrations for CCCP were lower by a factor of 30 compared with the traditional OECD 202 acute toxicity test. Mitochondrial effects were also detected at lower concentrations for all tested environmental contaminants compared to the OCED 202 test. For 2,4-dinitrophenol, mitochondria effects were detectable after 2h exposure to environmentally relevant concentrations and predicted organism death was observed after 24h. The high sensitivity and time efficiency of this novel automated imaging method make it a valuable tool for advancing ecotoxicological testing.

Given the need for modern researchers to produce open, reproducible scientific output, the lack of standards and best practices for sharing data and workflows used to produce and analyze molecular dynamics (MD) simulations has become an important issue in the field. There are now multiple well-established packages to perform molecular dynamics simulations, often highly tuned for exploiting specific classes of hardware, each with strong communities surrounding them, but with very limited interoperability/transferability options. Thus, the choice of the software package often dictates the workflow for both simulation production and analysis. The level of detail in documenting the workflows and analysis code varies greatly in published work, hindering reproducibility of the reported results and the ability for other researchers to build on these studies. An increasing number of researchers are motivated to make their data available, but many challenges remain in order to effectively share and reuse simulation data. To discuss these and other issues related to best practices in the field in general, we organized a workshop in November 2018 (https://bioexcel.eu/events/workshop-on-sharing-data-from-molecular-simulations/). Here, we present a brief overview of this workshop and topics discussed. We hope this effort will spark further conversation in the MD community to pave the way toward more open, interoperable, and reproducible outputs coming from research studies using MD simulations.

While important advances have been made in high-resolution mass spectrometry (HRMS) and its applications in non-targeted analysis (NTA), the number of identified compounds in biological and environmental samples often does not exceed 5% of the detected chemical features. Our aim was to develop a computational pipeline that leverages data from HRMS but also incorporates physicochemical properties (equilibrium partition ratios between organic solvents and water; K_{solvent–water}) and can propose molecular structures for detected chemical features. As these physicochemical properties are often sufficiently different across isomers, when put together, they can form a unique profile for each isomer, which we describe as the “physicochemical fingerprint”. In our study, we used a comprehensive database of compounds that have been previously reported in human blood and collected their K_{solvent–water} values for 129 partitioning systems. We used RDKit to calculate the number of RDKit fragments and the number of RDKit bits per molecule. We then developed and trained an artificial neural network, which used as an input the physicochemical fingerprint of a chemical feature and predicted the number and types of RDKit fragments and RDKit bits present in that structure. These were then used to search the database and propose chemical structures. The average success rate of predicting the right chemical structure ranged from 60 to 86% for the training set and from 48 to 81% for the testing set. These observations suggest that physicochemical fingerprints can assist in the identification of compounds with NTA and substantially improve the number of identified compounds.

Technological advances in transmission electron microscopes and detectors have turned cryogenic electron microscopy (cryo-EM) into an essential tool for structural biology. A commonly used cryo-EM data analysis method, single particle analysis, averages hundreds of thousands of low-dose images of individual macromolecular complexes to determine a density map of the complex. The presence of symmetry in the complex is beneficial since each projection image can be assigned to multiple views of the complex. However, data processing that applies symmetry can average out asymmetric features and consequently data analysis methods are required to resolve asymmetric structural features. Scipion is a cryo-EM image processing framework that integrates functions from different image processing packages as plugins. To extend its functionality for handling symmetry mismatches, we present here a Scipion plugin termed LocalRec implementing the localized reconstruction method. When tested on an adenovirus data set, the plugin enables resolving the symmetry-mismatched trimeric fibre bound to the five-fold vertices of the capsid. Furthermore, it improves the structure determination of the icosahedral capsid by dealing with the defocus gradient across the particle. LocalRec is expected to be widely applicable in a range of cryo-EM investigations of flexible and symmetry mismatched complexes.

T cells expressing the cutaneous lymphocyte antigen (CLA) mediate pathogenic inflammation in atopic dermatitis (AD). The molecular alterations contributing to their dysregulation remain unclear. With the aim to elucidate putative altered pathways in AD we profiled DNA methylation levels and miRNA expression in sorted T cell populations -(CD4(+), -CD4(+)CD45RA(+) naive, -CD4(+)CLA(+), and -CD8(+)) from adult AD patients and healthy controls (HC). Skin homing -CD4(+)CLA(+) T cells from AD patients showed significant differences in DNA methylation in 40 genes compared to HC (p < 0.05). Reduced DNA methylation levels in the upstream region of the interleukin-13 gene (IL13) in -CD4(+)CLA(+) T cells from AD patients correlated with increased IL13 mRNA expression in these cells. Sixteen miRNAs showed differential expression in -CD4(+)CLA(+) T cells from AD patients targeting genes in 202 biological processes (p < 0.05). An integrated network analysis of miRNAs and CpG sites identified two communities of strongly interconnected regulatory elements with strong antagonistic behaviours that recapitulated the differences between AD patients and HC. Functional analysis of the genes linked to these communities revealed their association with key cytokine signaling pathways, MAP kinase signaling and protein ubiquitination. Our findings support that epigenetic mechanisms play a role in the pathogenesis of AD by affecting inflammatory signaling molecules in skin homing -CD4(+)CLA(+) T cells and uncover putative molecules participating in AD pathways.

Helminth infections and allergic diseases are associated with IgE hyperresponsiveness but the genetics of this phenotype remain to be defined. Susceptibility to Ascaris lumbricoides infection and antibody levels to this helminth are associated with polymorphisms in locus 13q33-34. We aimed to explore this and other genomic regions to identify genetic variants associated with the IgE responsiveness in humans. Forty-eight subjects from Cartagena, Colombia, with extreme values of specific IgE to Ascaris and ABA-1, a resistance marker of this nematode, were selected for targeted resequencing. Burden analyses were done comparing extreme groups for IgE values. One-hundred one SNPs were genotyped in 1258 individuals of two well-characterized populations from Colombia and Sweden. Two low-frequency coding variants in the gene encoding the Acidic Mammalian Chitinase (CHIA rs79500525, rs139812869, tagged by rs10494133) were found enriched in high IgE responders to ABA-1 and confirmed by genetic association analyses. The SNP rs4950928 in the Chitinase 3 Like 1 gene (CHI3L1) was associated with high IgE to ABA-1 in Colombians and with high IgE to Bet v 1 in the Swedish population. CHIA rs10494133 and ABDH13 rs3783118 were associated with IgE responses to Ascaris. SNPs in the Tumor Necrosis Factor Superfamily Member 13b gene (TNFSF13B) encoding the cytokine B cell activating Factor were associated with high levels of total IgE in both populations. This is the first report on the association between low-frequency and common variants in the chitinases- related genes CHIA and CHI3L1 with the intensity of specific IgE to ABA-1 in a population naturally exposed to Ascaris and with Bet v 1 in a Swedish population. Our results add new information about the genetic influences of human IgE responsiveness; since the genes encode for enzymes involved in the immune response to parasitic infections, they could be helpful for understanding helminth immunity and allergic responses. We also confirmed that TNFSF13B has an important and conserved role in the regulation of total IgE levels, which supports potential evolutionary links between helminth immunity and allergic response.

DNA methylation changes may predispose becoming IgE-sensitized to allergens. We analyzed whether DNA methylation in peripheral blood mononuclear cells (PBMC) is associated with IgE sensitization at 5 years of age (5Y). DNA methylation was measured in 288 PBMC samples from 74 mother/child pairs from the birth cohort ALADDIN (Assessment of Lifestyle and Allergic Disease During INfancy) using the HumanMethylation450BeadChip (Illumina). PBMCs were obtained from the mothers during pregnancy and from their children in cord blood, at 2 years and 5Y. DNA methylation levels at each time point were compared between children with and without IgE sensitization to allergens at 5Y. For replication, CpG sites associated with IgE sensitization in ALADDIN were evaluated in whole blood DNA of 256 children, 4 years old, from the BAMSE (Swedish abbreviation for Children, Allergy, Milieu, Stockholm, Epidemiology) cohort. We found 34 differentially methylated regions (DMRs) associated with IgE sensitization to airborne allergens and 38 DMRs associated with sensitization to food allergens in children at 5Y (Sidak p <= 0.05). Genes associated with airborne sensitization were enriched in the pathway of endocytosis, while genes associated with food sensitization were enriched in focal adhesion, the bacterial invasion of epithelial cells, and leukocyte migration. Furthermore, 25 DMRs in maternal PBMCs were associated with IgE sensitization to airborne allergens in their children at 5Y, which were functionally annotated to the mTOR (mammalian Target of Rapamycin) signaling pathway. This study supports that DNA methylation is associated with IgE sensitization early in life and revealed new candidate genes for atopy. Moreover, our study provides evidence that maternal DNA methylation levels are associated with IgE sensitization in the child supporting early in utero effects on atopy predisposition.

Blue mold, caused by the necrotrophic fungal pathogen Penicillium expansum, causes serious postharvest losses in apple, and threatens human health through production of the potent mycotoxin patulin. Recent studies indicate a quantitative control of resistance against this disease in apple cultivars. A whole genome apple microarray covering 60k transcripts was used to identify gene(s) that appear to be differentially regulated between resistant and susceptible cultivars in P. expansum-infected fruits. A number of potential candidates was encountered among defense- and oxidative stress-related genes, cell wall modification and lignification genes, and genes related to localization and transport. Induction of one cell wall-related gene and three genes involved in the 'down-stream' flavonoid biosynthesis pathway, demonstrates the fundamental role of the cell wall as an important barrier, and suggests that fruit flavonoids are involved in the resistance to blue mold. Moreover, exogenous application of the plant hormone methyl jasmonate (MeJA) reduced the symptoms resulting from inoculating apples with P. expansum. This is the first report linking MeJA and activation of cell wall and flavonoid pathway genes to resistance against blue mold in a study comparing different cultivars of domesticated apple. Our results provide an initial categorization of genes that are potentially involved in the resistance mechanism, and should be useful for developing tools for gene marker-assisted breeding of apple cultivars with an improved resistance to blue mold.

Tannery wastewater containing chromium (Cr) is one of the most serious problems in leather industry. In order to develop an effective and eco-friendly treatment technology, a combined chemical-biological treatment system was performed for Cr remediation and recycling. The aim of the present study is to design a laboratory scale system using chemical precipitation of Cr(III) combined with biological removal of Cr(VI) from tannery wastewater, and to investigate the possibility of recycling the recovered Cr(III) in the tanning industry. Chemical precipitation of Cr(III) was carried out using lime and cement dust. The actinomycete strain Kitasatosporia sp. was used in microcosm studies for Cr(VI) bioremoval. Moreover, parameters such as type of porous medium, inoculum size, flow rate and culture conditions were investigated. The precipitated Cr(III) that was recovered from the chemical precipitation stage was recycled in the leather tanning industry. Our findings indicate that the maximum Cr(III) precipitation (98%) was achieved using 2 g/100 mL of lime and 2 h of settling rate. On the other hand, microcosm columns using sand that was inoculated with induced culture (OD600 = 2.43) and flow rate (2 mL/min) gave the maximum recovery (99%) of Cr(VI). The experimental Cr(III) was successfully recycled in the tanning process and the experimental leathers showed comparable properties as same as the leathers tanned with commercial Cr(III). Thus, we concluded that using combined chemical-biological treatment system for Cr remediation from tanning wastewater together with recycling process for the recovered Cr(III) is a promising strategy for economic and environmental friendly tanning industry.

The Lateglacial/early Holocene sediments from the ancient lake at Hasseldala Port, southern Sweden provide an important archive for the environmental and climatic shifts at the end of the last ice age and the transition into the present Interglacial. The existing multi-proxy data set highlights the complex interplay of physical and ecological changes in response to climatic shifts and lake status changes. Yet, it remains unclear how microorganisms, such as Archaea, which do not leave microscopic features in the sedimentary record, were affected by these climatic shifts. Here we present the metagenomic data set of Hasseldala Port with a special focus on the abundance and biodiversity of Archaea. This allows reconstructing for the first time the temporal succession of major Archaea groups between 13.9 and 10.8 ka BP by using ancient environmental DNA metagenomics and fossil archaeal cell membrane lipids. We then evaluate to which extent these findings reflect physical changes of the lake system, due to changes in lake-water summer temperature and seasonal lake-ice cover. We show that variations in archaeal composition and diversity were related to a variety of factors (e.g., changes in lake water temperature, duration of lake ice cover, rapid sediment infilling), which influenced bottom water conditions and the sediment-water interface. Methanogenic Archaea dominated during the Allerod and Younger Dryas pollen zones, when the ancient lake was likely stratified and anoxic for large parts of the year. The increase in archaeal diversity at the Younger Dryas/Holocene transition is explained by sediment infilling and formation of a mire/peatbog.

We use functionalized iron oxide magnetic multi-core particles of 100 nm in size (hydrodynamic particle diameter) and AC susceptometry (ACS) methods to measure the binding reactions between the magnetic nanoparticles (MNPs) and bio-analyte products produced from DNA segments using the rolling circle amplification (RCA) method. We use sensitive induction detection techniques in order to measure the ACS response. The DNA is amplified via RCA to generate RCA coils with a specific size that is dependent on the amplification time. After about 75 min of amplification we obtain an average RCA coil diameter of about 1 mu m. We determine a theoretical limit of detection (LOD) in the range of 11 attomole (corresponding to an analyte concentration of 55 fM for a sample volume of 200 mu L) from the ACS dynamic response after the MNPs have bound to the RCA coils and the measured ACS readout noise. We also discuss further possible improvements of the LOD.

The human mitochondria possess a dedicated set of ribosomes (mitoribosomes) that translate 13 essential protein components of the oxidative phosphorylation complexes encoded by the mitochondria! genome. Since all proteins synthesized by human mitoribosomes are integral membrane proteins, human mitoribosomes are tethered to the mitochondrial inner membrane during translation. Compared to the cytosolic ribosome the mitoribosome has a sedimentation coefficient of 55S, half the rRNA content, no 5S rRNA and 36 additional proteins. Therefore, a higher protein-to-RNA ratio and an atypical structure make the human mitoribosome substantially distinct from its cytosolic counterpart. Despite the central importance of the mitoribosome to life, no protocols were available to purify the intact complex from human cell lines. Traditionally, mitoribosomes were isolated from mitochondria-rich animal tissues that required kilograms of starting material. We reasoned that mitochondria in dividing HEK293-derived human cells grown in nutrient-rich expression medium would have an active mitochondrial translation, and, therefore, could be a suitable source of material for the structural and biochemical studies of the mitoribosome. To investigate its structure, we developed a protocol for large-scale purification of intact mitoribosomes from HEK cells. Herein, we introduce nitrogen cavitation method as a faster, less labor-intensive and more efficient alternative to traditional mechanical shear-based methods for cell lysis. This resulted in preparations of the mitoribosome that allowed for its structural determination to high resolution, revealing the composition of the intact human mitoribosome and its assembly intermediates. Here, we follow up on this work and present an optimized and more cost-effective method requiring only similar to 10(10) cultured HEK cells. The method can be employed to purify human mitoribosomal translating complexes, mutants, quality control assemblies and mitoribosomal subunits intermediates. The purification can be linearly scaled up tenfold if needed, and also applied to other types of cells.

Translation of mitochondrial messenger RNA (mt-mRNA) is performed by distinct mitoribosomes comprising at least 36 mitochondria-specific proteins. How these mitoribosomal proteins assist in the binding of mt-mRNA and to what extent they are involved in the translocation of transfer RNA (mt-tRNA) is unclear. To visualize the process of translation in human mitochondria, we report similar to 3.0 angstrom resolution structure of the human mitoribosome, including the L7/L12 stalk, and eight structures of its functional complexes with mt-mRNA, mt-tRNAs, recycling factor and additional trans factors. The study reveals a transacting protein module LRPPRC-SLIRP that delivers mt-mRNA to the mitoribosomal small subunit through a dedicated platform formed by the mitochondria-specific protein mS39. Mitoribosomal proteins of the large subunit mL40, mL48, and mL64 coordinate translocation of mt-tRNA. The comparison between those structures shows dynamic interactions between the mitoribosome and its ligands, suggesting a sequential mechanism of conformational changes.

The Lon protease is widely conserved in both prokaryotic and eukaryotic organisms and fulfills important regulatory functions. Nevertheless, the number of identified Lon substrates is limited in most organisms, and the precise role of Lon in regulating these proteins is poorly understood. Here, we describe the α-proteobacterial general stress response sigma factor σ^T as a novel Lon substrate in Caulobacter crescentus. Based on previously published quantitative proteomics data, we find σ^T to be a promising putative Lon substrate and confirm a direct role of Lon in degrading σ^T. We show that Lon contributes to the downregulation of σ^T abundance under optimal conditions and during recovery from sucrose-induced osmotic stress. Furthermore, the presence of the Lon activity regulator LarA enhances Lon-mediated degradation of σ^T in vitro and reduces σ^T levels in vivo indicating a role of LarA in modulating Lon-mediated degradation of σ^T. Together, our results highlight the importance of Lon during the recovery phase following stress exposure by adjusting the concentrations of critical regulators of stress responses.

Introduction

Acute respiratory distress syndrome (ARDS) is a common complication of influenza virus (IV) infection. During ARDS, alveolar protein concentrations often reach 40-90% of plasma levels, causing severe impairment of gas exchange and promoting deleterious alveolar remodeling. Protein clearance from the alveolar space is at least in part facilitated by the multi-ligand receptor megalin through clathrin-mediated endocytosis.

Methods

To investigate whether IV infection impairs alveolar protein clearance, we examined albumin uptake and megalin expression in MLE-12 cells and alveolar epithelial cells (AEC) from murine precision-cut lung slices (PCLS) and in vivo, under IV infection conditions by flow cytometry and western blot. Transcriptional levels from AEC and broncho-alveolar lavage (BAL) cells were analyzed in an in-vivo mouse model by RNAseq.

Results

IV significantly downregulated albumin uptake, independently of activation of the TGF- β1/GSK3β axis that has been previously implicated in the regulation of megalin function. Decreased plasma membrane abundance, total protein levels, and mRNA expression of megalin were associated with this phenotype. In IV-infected mice, we identified a significant upregulation of matrix metalloproteinase (MMP)-14 in BAL fluid cells. Furthermore, the inhibition of this protease partially recovered total megalin levels and albumin uptake.

Discussion

Our results suggest that the previously described MMP-driven shedding mechanisms are potentially involved in downregulation of megalin cell surface abundance and clearance of excess alveolar protein. As lower alveolar edema protein concentrations are associated with better outcomes in respiratory failure, our findings highlight the therapeutic potential of a timely MMP inhibition in the treatment of IV-induced ARDS.

Cyanobacteria are widely distributed primary producers with significant implications for the global biogeochemical cycles of carbon and nitrogen. Diazotrophic cyanobacteria of subsection V (Order Stigonematales) are particularly ubiquitous in photoautotrophic microbial mats of hot springs. The Stigonematal cyanobacterium strain CHPI isolated from the Porcelana hot spring (Chile) was one of the major contributors of the new nitrogen through nitrogen fixation. Further morphological and genetic characterization verified that the strain CHP1 belongs to Stigonematales, and it formed a separate Glade together with other thermophiles of the genera Fischerella and Mastigocladus. Strain CHP1 fixed maximum N-2 in the light, independent of the temperature range. At 50 degrees C niJH gene transcripts showed high expression during the light period, whereas the nifH gene expression at 45 degrees C was arrhythmic. The strain displayed a high affinity for nitrate and a low tolerance for high ammonium concentrations, whereas the narB and glnA genes showed higher expression in light and at the beginning of the dark phase. It is proposed that Mastigocladus sp. strain CHPI would represent a good model for the study of subsection V thermophilic cyanobacteria, and for understanding the adaptations of these photoautotrophic organisms inhabiting microbial mats in hot springs globally.

RT-LAMP detection of SARS-CoV-2 has been shown to be a valuable approach to scale up COVID-19 diagnostics and thus contribute to limiting the spread of the disease. Here we present the optimization of highly cost-effective in-house produced enzymes, and we benchmark their performance against commercial alternatives. We explore the compatibility between multiple DNA polymerases with high strand-displacement activity and thermostable reverse transcriptases required for RT-LAMP. We optimize reaction conditions and demonstrate their applicability using both synthetic RNA and clinical patient samples. Finally, we validate the optimized RT-LAMP assay for the detection of SARS-CoV-2 in unextracted heat-inactivated nasopharyngeal samples from 184 patients. We anticipate that optimized and affordable reagents for RT-LAMP will facilitate the expansion of SARS-CoV-2 testing globally, especially in sites and settings where the need for large scale testing cannot be met by commercial alternatives.

Background: Sampling genomes with Fosmid vectors and sequencing of pooled Fosmid libraries on the Illumina platform for massive parallel sequencing is a novel and promising approach to optimizing the trade-off between sequencing costs and assembly quality. Results: In order to sequence the genome of Norway spruce, which is of great size and complexity, we developed and applied a new technology based on the massive production, sequencing, and assembly of Fosmid pools (FP). The spruce chromosomes were sampled with similar to 40,000 bp Fosmid inserts to obtain around two-fold genome coverage, in parallel with traditional whole genome shotgun sequencing (WGS) of haploid and diploid genomes. Compared to the WGS results, the contiguity and quality of the FP assemblies were high, and they allowed us to fill WGS gaps resulting from repeats, low coverage, and allelic differences. The FP contig sets were further merged with WGS data using a novel software package GAM-NGS. Conclusions: By exploiting FP technology, the first published assembly of a conifer genome was sequenced entirely with massively parallel sequencing. Here we provide a comprehensive report on the different features of the approach and the optimization of the process. We have made public the input data (FASTQ format) for the set of pools used in this study: ftp://congenie.org/congenie/Nystedt_2013/Assembly/ProcessedData/FosmidPools/.(alternatively accessible via http://congenie.org/downloads).The software used for running the assembly process is available at http://research.scilifelab.se/andrej_alexeyenko/downloads/fpools/.

FunCoup (http://FunCoup.sbc.su.se) is a database that maintains and visualizes global gene/protein networks of functional coupling that have been constructed by Bayesian integration of diverse high-throughput data. FunCoup achieves high coverage by orthology-based integration of data sources from different model organisms and from different platforms. We here present release 2.0 in which the data sources have been updated and the methodology has been refined. It contains a new data type Genetic Interaction, and three new species: chicken, dog and zebra fish. As FunCoup extensively transfers functional coupling information between species, the new input datasets have considerably improved both coverage and quality of the networks. The number of high-confidence network links has increased dramatically. For instance, the human network has more than eight times as many links above confidence 0.5 as the previous release. FunCoup provides facilities for analysing the conservation of subnetworks in multiple species. We here explain how to do comparative interactomics on the FunCoup website.

Background: MCMC-based methods are important for Bayesian inference of phylogeny and related parameters. Although being computationally expensive, MCMC yields estimates of posterior distributions that are useful for estimating parameter values and are easy to use in subsequent analysis. There are, however, sometimes practical difficulties with MCMC, relating to convergence assessment and determining burn-in, especially in large-scale analyses. Currently, multiple software are required to perform, e.g., convergence, mixing and interactive exploration of both continuous and tree parameters.

Results: We have written a software called VMCMC to simplify post-processing of MCMC traces with, for example, automatic burn-in estimation. VMCMC can also be used both as a GUI-based application, supporting interactive exploration, and as a command-line tool suitable for automated pipelines.

Conclusions: VMCMC is a free software available under the New BSD License. Executable jar files, tutorial manual and source code can be downloaded from https://bitbucket. org/rhali/visualmcmc/.

Background: Homology inference is pivotal to evolutionary biology and is primarily based on significant sequence similarity, which, in general, is a good indicator of homology. Algorithms have also been designed to utilize conservation in gene order as an indication of homologous regions. We have developed GenFamClust, a method based on quantification of both gene order conservation and sequence similarity. Results: In this study, we validate GenFamClust by comparing it to well known homology inference algorithms on a synthetic dataset. We applied several popular clustering algorithms on homologs inferred by GenFamClust and other algorithms on a metazoan dataset and studied the outcomes. Accuracy, similarity, dependence, and other characteristics were investigated for gene families yielded by the clustering algorithms. GenFamClust was also applied to genes from a set of complete fungal genomes and gene families were inferred using clustering. The resulting gene families were compared with a manually curated gold standard of pillars from the Yeast Gene Order Browser. We found that the gene-order component of GenFamClust is simple, yet biologically realistic, and captures local synteny information for homologs. Conclusions: The study shows that GenFamClust is a more accurate, informed, and comprehensive pipeline to infer homologs and gene families than other commonly used homology and gene-family inference methods.

Background

Clustering sequences into families has long been an important step in characterization of genes and proteins. There are many algorithms developed for this purpose, most of which are based on either direct similarity between gene pairs or some sort of network structure, where weights on edges of constructed graphs are based on similarity. However, conserved synteny is an important signal that can help distinguish homology and it has not been utilized to its fullest potential.

Results

Here, we present GenFamClust, a pipeline that combines the network properties of sequence similarity and synteny to assess homology relationship and merge known homologs into groups of gene families. GenFamClust identifies homologs in a more informed and accurate manner as compared to similarity based approaches. We tested our method against the Neighborhood Correlation method on two diverse datasets consisting of fully sequenced genomes of eukaryotes and synthetic data.

Conclusions

The results obtained from both datasets confirm that synteny helps determine homology and GenFamClust improves on Neighborhood Correlation method. The accuracy as well as the definition of synteny scores is the most valuable contribution of GenFamClust.

Metagenomic and metatranscriptomic data were generated from size-fractionated samples from 11 sites within the Baltic Sea and adjacent marine waters of Kattegat and freshwater Lake Tornetrask in order to investigate the diversity, distribution, and transcriptional activity of virioplankton. Such a transect, spanning a salinity gradient from freshwater to the open sea, facilitated a broad genome-enabled investigation of natural as well as impacted aspects of Baltic Sea viral communities. Taxonomic signatures representative of phages within the widely distributed order Caudovirales were identified with enrichments in lesser-known families such as Podoviridae and Siphoviridae. The distribution of phage reported to infect diverse and ubiquitous heterotrophic bacteria (SAR11 clades) and cyanobacteria (Synechococcus sp.) displayed population-level shifts in diversity. Samples from higher-salinity conditions (>14 practical salinity units [PSU]) had increased abundances of viruses for picoeukaryotes, i.e., Ostreococcus. These data, combined with host diversity estimates, suggest viral modulation of diversity on the whole-community scale, as well as in specific prokaryotic and eukaryotic lineages. RNA libraries revealed single-stranded DNA (ssDNA) and RNA viral populations throughout the Baltic Sea, with ssDNA phage highly represented in Lake Tornetrask. Further, our data suggest relatively high transcriptional activity of fish viruses within diverse families known to have broad host ranges, such as Nodoviridae (RNA), Iridoviridae (DNA), and predicted zoonotic viruses that can cause ecological and economic damage as well as impact human health. IMPORTANCE Inferred virus-host relationships, community structures of ubiquitous ecologically relevant groups, and identification of transcriptionally active populations have been achieved with our Baltic Sea study. Further, these data, highlighting the transcriptional activity of viruses, represent one of the more powerful uses of omics concerning ecosystem health. The use of omics-related data to assess ecosystem health holds great promise for rapid and relatively inexpensive determination of perturbations and risk, explicitly with regard to viral assemblages, as no single marker gene is suitable for widespread taxonomic coverage.

The advent of machine learning-based structure prediction algorithms such as AlphaFold2 (AF2) and RoseTTa Fold have moved the generation of accurate structural models for the entire cellular protein machinery into the reach of the scientific community. However, structure predictions of protein complexes are based on user-provided input and may require experimental validation. Mass spectrometry (MS) is a versatile, time-effective tool that provides information on post-translational modifications, ligand interactions, conformational changes, and higher-order oligomerization. Using three protein systems, we show that native MS experiments can uncover structural features of ligand interactions, homology models, and point mutations that are undetectable by AF2 alone. We conclude that machine learning can be complemented with MS to yield more accurate structural models on a small and large scale.

Signal peptides (SPs) are short amino acid sequences in the amino terminus of many newly synthesized proteins that target proteins into, or across, membranes. Bioinformatic tools can predict SPs from amino acid sequences, but most cannot distinguish between various types of signal peptides. We present a deep neural network-based approach that improves SP prediction across all domains of life and distinguishes between three types of prokaryotic SPs.

To improve the mechanistic screening of reproductive toxicants in  chemical-risk assessment and drug development, we have developed a three-dimensional (3D) heterogenous testicular co-culture model from neonatal mice. Di-n-butyl phthalate (DBP), an environmental contaminant that can affect reproductive health negatively, was used as a model compound to illustrate the utility of the in vitro model. The cells were treated with DBP (1 nM to 100 µM) for 7 days. Automated high-content imaging confirmed the presence of cell-specific markers of Leydig cells (CYP11A1 +), Sertoli cells (SOX9 +), and germ cells (DAZL +). Steroidogenic activity of Leydig cells was demonstrated by analyzing testosterone levels in the culture medium. DBP induced a concentration-dependent reduction in testosterone levels and decreased the number of Leydig cells compared to vehicle control. The levels of steroidogenic regulator StAR and the steroidogenic enzyme CYP11A1 were decreased already at the lowest DBP concentration (1 nM), demonstrating upstream effects in the testosterone biosynthesis pathway. Furthermore, exposure to 10 nM DBP decreased the levels of the germ cell-specific RNA binding protein DAZL, central for the spermatogenesis. The 3D model also captured the development of the Sertoli cell junction proteins, N-cadherin and Zonula occludens protein 1 (ZO-1), critical for the blood–testis barrier. However, DBP exposure did not significantly alter the cadherin and ZO-1 levels. Altogether, this 3D in vitro system models testicular cellular signaling and function, making it a powerful tool for mechanistic screening of developmental testicular toxicity. This can open a new avenue for high throughput screening of chemically-induced reproductive toxicity during sensitive developmental phases.

Di-n-butyl phthalate (DBP) is a ubiquitous environmental contaminant linked with various adverse health effects, including immune system dysfunction. Gut microbial dysbiosis can contribute to a wide range of pathogenesis, particularly immune disease. Here, we investigated the impact of DBP on the gut microbiome and examined correlations with immune system changes after five weeks oral exposure (10 or 100 mg/kg/day) in adult male mice. The fecal microbiome composition was characterized using 16S rRNA sequencing. DBP-treated mice displayed a significantly distinct microbial community composition, indicated by Bray-Curtis distance. Numerous amplicon sequence variants (ASVs) at the genus level were altered. Compared to the vehicle control group, the 10 mg/kg/day DBP group had 63 more abundant and 65 less abundant ASVs, while 60 ASVs were increased and 76 ASVs were decreased in the 100 mg/kg/day DBP group. Both DBP treatment groups showed higher abundances of ASVs assigned to Desulfovibrio (Proteobacteria phylum) and Enterorhabdus genera, while ASVs belonging to Parabacteroides, Lachnospiraceae UCG-006 and Lachnoclostridium were less common compared to the control group. Interestingly, an ASV belonging to Rumniniclostridium 6, which was less abundant in DBP-treated mice, demonstrated a negative correlation with the increased number of non-classical monocytes observed in the blood of DBP-treated animals. In addition, an ASV from Lachnospiraceae UCG-001, which was more abundant in the DBP-treated animals, showed a positive correlation with the non-classical monocyte increase. This study shows that DBP exposure greatly modifies the gut bacterial microbiome and indicates a potential contribution of microbial dysbiosis to DBP-induced immune system impairment, illustrating the importance of investigating how interactions between exposome components can affect health.

Background: Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

Results: Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

Conclusions: Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

The Baltic Sea is one of the world's largest brackish water bodies and is characterised by pronounced physicochemical gradients where microbes are the main biogeochemical catalysts. Meta-omic methods provide rich information on the composition of, and activities within, microbial ecosystems, but are computationally heavy to perform. We here present the Baltic Sea Reference Metagenome (BARM), complete with annotated genes to facilitate further studies with much less computational effort. The assembly is constructed using 2.6 billion metagenomic reads from 81 water samples, spanning both spatial and temporal dimensions, and contains 6.8 million genes that have been annotated for function and taxonomy. The assembly is useful as a reference, facilitating taxonomic and functional annotation of additional samples by simply mapping their reads against the assembly. This capability is demonstrated by the successful mapping and annotation of 24 external samples. In addition, we present a public web interface, BalticMicrobeDB, for interactive exploratory analysis of the dataset. [GRAPHICS] .

Achieving high accuracy in orthology inference is essential for many comparative, evolutionary and functional genomic analyses, yet the true evolutionary history of genes is generally unknown and orthologs are used for very different applications across phyla, requiring different precision-recall trade-offs. As a result, it is difficult to assess the performance of orthology inference methods. Here, we present a community effort to establish standards and an automated web-based service to facilitate orthology benchmarking. Using this service, we characterize 15 well-established inference methods and resources on a battery of 20 different benchmarks. Standardized benchmarking provides a way for users to identify the most effective methods for the problem at hand, sets a minimum requirement for new tools and resources, and guides the development of more accurate orthology inference methods.

The identification of orthologs-genes in different species which descended from the same gene in their last common ancestor-is a prerequisite for many analyses in comparative genomics and molecular evolution. Numerous algorithms and resources have been conceived to address this problem, but benchmarking and interpreting them is fraught with difficulties (need to compare them on a common input dataset, absence of ground truth, computational cost of calling orthologs). To address this, the Quest for Orthologs consortium maintains a reference set of proteomes and provides a web server for continuous orthology benchmarking (http://orthology.benchmarkservice.org). Furthermore, consensus ortholog calls derived from public benchmark submissions are provided on the Alliance of Genome Resources website, the joint portal of NIH-funded model organism databases.

Genome-scale metabolic models (GSMMs) can provide novel insights into metabolic reprogramming during disease progression and therapeutic interventions. We developed a context-specific system-level GSMM of people living with HIV (PLWH) using global RNA sequencing data from PBMCs with suppressive viremia either by natural (elite controllers, PLWHEC) or drug-induced (PLWHART) control. This GSMM was compared with HIV-negative controls (HC) to provide a comprehensive systems-level metabo-transcriptomic characterization. Transcriptomic analysis identified up-regulation of oxidative phosphorylation as a characteristic of PLWHART, differentiating them from PLWHEC with dysregulated complexes I, III, and IV. The flux balance analysis identified altered flux in several intermediates of glycolysis including pyruvate, a-ketoglutarate, and glutamate, among others, in PLWHART. The in vitro pharmacological inhibition of OXPHOS complexes in a latent lymphocytic cell model (J-Lat 10.6) suggested a role for complex IV in latency reversal and immunosenescence. Furthermore, inhibition of complexes I/III/IV induced apoptosis, collectively indicating their contribution to reservoir dynamics.

The clinical outcome and disease severity in coronavirus disease 2019 (COVID-19) are heterogeneous, and the progression or fatality of the disease cannot be explained by a single factor like age or comorbidities. In this study, we used system-wide network-based system biology analysis using whole blood RNA sequencing, immunophenotyping by flow cytometry, plasma metabolomics, and single-cell-type metabolo-mics of monocytes to identify the potential determinants of COVID-19 severity at personalized and group levels. Digital cell quantification and immunophenotyping of the mononuclear phagocytes indicated a sub-stantial role in coordinating the immune cells that mediate COVID-19 severity. Stratum-specific and person-alized genome-scale metabolic modeling indicated monocarboxylate transporter family genes (e.g., SLC16A6), nucleoside transporter genes (e.g., SLC29A1), and metabolites such as a-ketoglutarate, succi-nate, malate, and butyrate could play a crucial role in COVID-19 severity. Metabolic perturbations targeting the central metabolic pathway (TCA cycle) can be an alternate treatment strategy in severe COVID-19.

The current human reference sequence (GRCh38) is a foundation for large-scale sequencing projects. However, recent studies have suggested that GRCh38 may be incomplete and give a suboptimal representation of specific population groups. Here, we performed a de novo assembly of two Swedish genomes that revealed over 10 Mb of sequences absent from the human GRCh38 reference in each individual. Around 6 Mb of these novel sequences (NS) are shared with a Chinese personal genome. The NS are highly repetitive, have an elevated GC-content, and are primarily located in centromeric or telomeric regions. Up to 1 Mb of NS can be assigned to chromosome Y, and large segments are also missing from GRCh38 at chromosomes 14, 17, and 21. Inclusion of NS into the GRCh38 reference radically improves the alignment and variant calling from short-read whole-genome sequencing data at several genomic loci. A re-analysis of a Swedish population-scale sequencing project yields > 75,000 putative novel single nucleotide variants (SNVs) and removes > 10,000 false positive SNV calls per individual, some of which are located in protein coding regions. Our results highlight that the GRCh38 reference is not yet complete and demonstrate that personal genome assemblies from local populations can improve the analysis of short-read whole-genome sequencing data.

Here we describe the SweGen data set, a comprehensive map of genetic variation in the Swedish population. These data represent a basic resource for clinical genetics laboratories as well as for sequencing-based association studies by providing information on genetic variant frequencies in a cohort that is well matched to national patient cohorts. To select samples for this study, we first examined the genetic structure of the Swedish population using high-density SNP-array data from a nation-wide cohort of over 10 000 Swedish-born individuals included in the Swedish Twin Registry. A total of 1000 individuals, reflecting a cross-section of the population and capturing the main genetic structure, were selected for whole-genome sequencing. Analysis pipelines were developed for automated alignment, variant calling and quality control of the sequencing data. This resulted in a genome-wide collection of aggregated variant frequencies in the Swedish population that we have made available to the scientific community through the website https://swefreq.nbis.se. A total of 29.2 million single-nucleotide variants and 3.8 million indels were detected in the 1000 samples, with 9.9 million of these variants not present in current databases. Each sample contributed with an average of 7199 individual-specific variants. In addition, an average of 8645 larger structural variants (SVs) were detected per individual, and we demonstrate that the population frequencies of these SVs can be used for efficient filtering analyses. Finally, our results show that the genetic diversity within Sweden is substantial compared with the diversity among continental European populations, underscoring the relevance of establishing a local reference data set.

Ionic liquids (ILs) are a class of salts existing in the liquid state below 100 degrees C, possessing low volatility, high thermal stability as well as many highly attractive solvent and electrochemical capabilities, etc., making them highly tunable for a great variety of applications, such as lubricants, electrolytes, and soft functional materials. In many applications, ILs are first either physi- or chemisorbed on a solid surface to successively create more functional materials. The functions of ILs at solid surfaces can differ considerably from those of bulk ILs, mainly due to distinct interfacial layers with tunable structures resulting in new ionic liquid interface layer properties and enhanced performance. Due to an almost infinite number of possible combinations among the cations and anions to form ILs, the diversity of various solid surfaces, as well as different external conditions and stimuli, a detailed molecular-level understanding of their structure-property relationship is of utmost significance for a judicious design of IL-solid interfaces with appropriate properties for task-specific applications. Many experimental techniques, such as atomic force microscopy, surface force apparatus, and so on, have been used for studying the ion structuring of the IL interface layer. Molecular Dynamics simulations have been widely used to investigate the microscopic behavior of the IL interface layer. To interpret and clarify the IL structure and dynamics as well as to predict their properties, it is always beneficial to combine both experiments and simulations as close as possible. In another theoretical model development to bridge the structure and properties of the IL interface layer with performance, thermodynamic prediction & property modeling has been demonstrated as an effective tool to add the properties and function of the studied nanomaterials. Herein, we present recent findings from applying the multiscale triangle experiment-simulation-thermodynamic modeling in the studies of ion structuring of ILs in the vicinity of solid surfaces, as well as how it qualitatively and quantitatively correlates to the overall ILs properties, performance, and function. We introduce the most common techniques behind experiment-simulation-thermodynamic modeling and how they are applied for studying the IL interface layer structuring, and we highlight the possibilities of the IL interface layer structuring in applications such as lubrication and energy storage. Integrative experiment-simulation-thermodynamic modeling is highly demanded for qualitatively and quantitatively correlating the ionic liquids interface layer structuring to the overall properties, performance, and function.

Gut microbes play important roles for their hosts. Previous studies suggest that host-microbial systems can form long-term associations over evolutionary time and the dynamic changes of the intestinal system may represent major driving forces and contribute to insect dietary diversification and speciation. Our study system includes a set of six closely related leaf beetle species (Galerucella spp.) and our study aims to separate the roles of host phylogeny and ecology in determining the gut microbial community and to identify eventual relationship between host insects and gut bacteria. We collected adult beetles from their respective host plants and quantified their microbial community using 16S rRNA sequencing. The results showed that the gut bacteria community composition was structured by host beetle phylogeny, where more or less host-specific gut bacteria interact with the different Galerucella species. For example, the endosymbiotic bacteria Wolbachia was found almost exclusively in G. nymphaea and G. sagittariae. Diversity indicators also suggested that α- and β-diversities of gut bacteria communities varied among host beetle species. Overall, our results suggest a phylogenetically controlled co-occurrence pattern between the six closely related Galerucella beetles and their gut bacteria, indicating the potential of co-evolutionary processes occurring between hosts and their gut bacterial communities.

The malarial parasite Plasmodium exports its own proteins to the cell surfaces of red blood cells (RBCs) during infection. Examples of exported proteins include members of the repetitive interspersed family (RIFIN) and subtelomeric variable open reading frame (STEVOR) family of proteins from Plasmodium falciparum. The presence of these parasite-derived proteins on surfaces of infected RBCs triggers the adhesion of infected cells to uninfected cells (rosetting) and to the vascular endothelium potentially obstructing blood flow. While there is a fair amount of information on the localization of these proteins on the cell surfaces of RBCs, less is known about how they can be exported to the membrane and the topologies they can adopt during the process. The first step of export is plausibly the cotranslational insertion of proteins into the endoplasmic reticulum (ER) of the parasite, and here, we investigate the insertion of three RIFIN and two STEVOR proteins into the ER membrane. We employ a well-established experimental system that uses N-linked glycosylation of sites within the protein as a measure to assess the extent of membrane insertion and the topology it assumes when inserted into the ER membrane. Our results indicate that for all the proteins tested, transmembranes (TMs) 1 and 3 integrate into the membrane, so that the protein assumes an overall topology of Ncyt-Ccyt. We also show that the segment predicted to be TM2 for each of the proteins likely does not reside in the membrane, but is translocated to the lumen.

The discriminatory power of the noncoding control region (CR) of domestic dog mitochondrial DNA alone is relatively low. The extent to which the discriminatory power could be increased by analyzing additional highly variable coding regions of the mitochondrial genome (mtGenome) was therefore investigated. Genetic variability across the mtGenome was evaluated by phylogenetic analysis, and the three most variable similar to 1kb coding regions identified. We then sampled 100 Swedish dogs to represent breeds in accordance with their frequency in the Swedish population. A previously published dataset of 59 dog mtGenomes collected in the United States was also analyzed. Inclusion of the three coding regions increased the exclusion capacity considerably for the Swedish sample, from 0.920 for the CR alone to 0.964 for all four regions. The number of mtDNA types among all 159 dogs increased from 41 to 72, the four most frequent CR haplotypes being resolved into 22 different haplotypes.

The detection of recurrent somatic chromosomal rearrangements is standard of care for most leukemia types. Even though karyotype analysis-a low-resolution genome-wide chromosome analysis-is still the gold standard, it often needs to be complemented with other methods to increase resolution. To evaluate the feasibility and applicability of mate pair whole genome sequencing (MP-WGS) to detect structural chromosomal rearrangements in the diagnostic setting, we sequenced ten bone marrow samples from leukemia patients with recurrent rearrangements. Samples were selected based on cytogenetic and FISH results at leukemia diagnosis to include common rearrangements of prognostic relevance. Using MP-WGS and in-house bioinformatic analysis all sought rearrangements were successfully detected. In addition, unexpected complexity or additional, previously undetected rearrangements was unraveled in three samples. Finally, the MP-WGS analysis pinpointed the location of chromosome junctions at high resolution and we were able to identify the exact exons involved in the resulting fusion genes in all samples and the specific junction at the nucleotide level in half of the samples. The results show that our approach combines the screening character from karyotype analysis with the specificity and resolution of cytogenetic and molecular methods. As a result of the straightforward analysis and high-resolution detection of clinically relevant rearrangements, we conclude that MP-WGS is a feasible method for routine leukemia diagnostics of structural chromosomal rearrangements.

Since the original publication of sRNAtoolbox in 2015, small RNA research experienced notable advances in different directions. New protocols for small RNA sequencing have become available to address important issues such as adapter ligation bias, PCR amplification artefacts or to include internal controls such as spike-in sequences. New microRNA reference databases were developed with different foci, either prioritizing accuracy (low number of false positives) or completeness (low number of false negatives). Additionally, other small RNA molecules as well asmicroRNA sequence and length variants (isomiRs) have continued to gain importance. Finally, the number of microRNA sequencing studies deposited in GEO nearly triplicated from 2014 (280) to 2018 (764). These developments imply that fast and easy-to-use tools for expression profiling and subsequent downstream analysis of miRNAseq data are essential to many researchers. Key features in this sRNAtoolbox release include addition of all major RNA library preparation protocols to sRNAbench and improvements in sRNAde, a tool that summarizes several aspects of small RNA sequencing studies including the detection of consensus differential expression. A special emphasis was put on the user-friendliness of the tools, for instance sRNAbench now supports parallel launching of several jobs to improve reproducibility and user time efficiency.