SSAGES : Software Suite for Advanced General Ensemble Simulations

Molecular simulation has emerged as an essential tool for modern-day research, but obtaining proper results and making reliable conclusions from simulations requires adequate sampling of the system under consideration. To this end, a variety of methods exist in the literature that can enhance sampling considerably, and increasingly sophisticated, effective algorithms continue to be developed at a rapid pace. Implementation of these techniques, however, can be challenging for experts and non-experts alike. There is a clear need for software that provides rapid, reliable, and easy access to a wide range of advanced sampling methods and that facilitates implementation of new techniques as they emerge. Here we present SSAGES, a publicly available Software Suite for Advanced General Ensemble Simulations designed to interface with multiple widely used molecular dynamics simulations packages. SSAGES allows facile application of a variety of enhanced sampling techniques—including adaptive biasing force, string methods, and forward flux sampling—that extract meaningful free energy and transition path data from all-atom and coarse-grained simulations. A noteworthy feature of SSAGES is a user-friendly framework that facilitates further development and implementation of new methods and collective variables. In this work, the use of SSAGES is illustrated in the context of simple representative applications involving distinct methods and different collective variables that are available in the current release of the suite. The code may be found at: https://github.com/MICCoM/SSAGES-public

Ensembl Genomes: an integrative resource for genome-scale data from non-vertebrate species

Ensembl Genomes (http://www.ensemblgenomes.org) is an integrative resource for genome-scale data from non-vertebrate species. The project exploits and extends technology (for genome annotation, analysis and dissemination) developed in the context of the (vertebrate-focused) Ensembl project and provides a complementary set of resources for non-vertebrate species through a consistent set of programmatic and interactive interfaces. These provide access to data including reference sequence, gene models, transcriptional data, polymorphisms and comparative analysis. Since its launch in 2009, Ensembl Genomes has undergone rapid expansion, with the goal of providing coverage of all major experimental organisms, and additionally including taxonomic reference points to provide the evolutionary context in which genes can be understood. Against the backdrop of a continuing increase in genome sequencing activities in all parts of the tree of life, we seek to work, wherever possible, with the communities actively generating and using data, and are participants in a growing range of collaborations involved in the annotation and analysis of genomes

Silencing of Vlaro2 for chorismate synthase revealed that the phytopathogen Verticillium longisporum induces the cross-pathway control in the xylem

The first leaky auxotrophic mutant for aromatic amino acids of the near-diploid fungal plant pathogen Verticillium longisporum (VL) has been generated. VL enters its host Brassica napus through the roots and colonizes the xylem vessels. The xylem contains little nutrients including low concentrations of amino acids. We isolated the gene Vlaro2 encoding chorismate synthase by complementation of the corresponding yeast mutant strain. Chorismate synthase produces the first branch point intermediate of aromatic amino acid biosynthesis. A novel RNA-mediated gene silencing method reduced gene expression of both isogenes by 80% and resulted in a bradytrophic mutant, which is a leaky auxotroph due to impaired expression of chorismate synthase. In contrast to the wild type, silencing resulted in increased expression of the cross-pathway regulatory gene VlcpcA (similar to cpcA/GCN4) during saprotrophic life. The mutant fungus is still able to infect the host plant B. napus and the model Arabidopsis thaliana with reduced efficiency. VlcpcA expression is increased in planta in the mutant and the wild-type fungus. We assume that xylem colonization requires induction of the cross-pathway control, presumably because the fungus has to overcome imbalanced amino acid supply in the xylem

Evolutionary genomics of a cold-adapted diatom: Fragilariopsis cylindrus

The Southern Ocean houses a diverse and productive community of organisms1, 2. Unicellular eukaryotic diatoms are the main primary producers in this environment, where photosynthesis is limited by low concentrations of dissolved iron and large seasonal fluctuations in light, temperature and the extent of sea ice3, 4, 5, 6, 7. How diatoms have adapted to this extreme environment is largely unknown. Here we present insights into the genome evolution of a cold-adapted diatom from the Southern Ocean, Fragilariopsis cylindrus8, 9, based on a comparison with temperate diatoms. We find that approximately 24.7 per cent of the diploid F. cylindrus genome consists of genetic loci with alleles that are highly divergent (15.1 megabases of the total genome size of 61.1 megabases). These divergent alleles were differentially expressed across environmental conditions, including darkness, low iron, freezing, elevated temperature and increased CO2. Alleles with the largest ratio of non-synonymous to synonymous nucleotide substitutions also show the most pronounced condition-dependent expression, suggesting a correlation between diversifying selection and allelic differentiation. Divergent alleles may be involved in adaptation to environmental fluctuations in the Southern Ocean

Genome Characterization of the Oleaginous Fungus Mortierella alpina

Mortierella alpina is an oleaginous fungus which can produce lipids accounting for up to 50% of its dry weight in the form of triacylglycerols. It is used commercially for the production of arachidonic acid. Using a combination of high throughput sequencing and lipid profiling, we have assembled the M. alpina genome, mapped its lipogenesis pathway and determined its major lipid species. The 38.38 Mb M. alpina genome shows a high degree of gene duplications. Approximately 50% of its 12,796 gene models, and 60% of genes in the predicted lipogenesis pathway, belong to multigene families. Notably, M. alpina has 18 lipase genes, of which 11 contain the class 2 lipase domain and may share a similar function. M. alpina's fatty acid synthase is a single polypeptide containing all of the catalytic domains required for fatty acid synthesis from acetyl-CoA and malonyl-CoA, whereas in many fungi this enzyme is comprised of two polypeptides. Major lipids were profiled to confirm the products predicted in the lipogenesis pathway. M. alpina produces a complex mixture of glycerolipids, glycerophospholipids and sphingolipids. In contrast, only two major sterol lipids, desmosterol and 24(28)-methylene-cholesterol, were detected. Phylogenetic analysis based on genes involved in lipid metabolism suggests that oleaginous fungi may have acquired their lipogenic capacity during evolution after the divergence of Ascomycota, Basidiomycota, Chytridiomycota and Mucoromycota. Our study provides the first draft genome and comprehensive lipid profile for M. alpina, and lays the foundation for possible genetic engineering of M. alpina to produce higher levels and diverse contents of dietary lipids

Correction: Quevillon, M.J., et al. Charge Transport and Phase Behavior of Imidazolium-Based Ionic Liquid Crystals from Fully Atomistic Simulations. Materials 2018, 11, 64

In our publication [...

The Backfolded Odijk Regime for Wormlike Chains Confined in Rectangular Nanochannels

We confirm Odijk’s scaling laws for (i) the average chain extension; (ii) the variance about the average extension; and (iii) the confinement free energy of a wormlike chain confined in a rectangular nanochannel smaller than its chain persistence length through pruned-enriched Rosenbluth method (PERM) simulations of asymptotically long, discrete wormlike chains. In the course of this analysis, we also computed the global persistence length of ideal wormlike chains for the modestly rectangular channels that are used in many experimental systems. The results are relevant to genomic mapping systems that confine DNA in channel sizes around 50 nm, since fabrication constraints generally lead to rectangular cross-sections