A muscle-related contractile tissue specified by myocardin-related transcription factor activity in Porifera Muscle-based movement is a hallmark of animal biology, but the evolutionary origins of myocytes - the cells that comprise muscle tissues - are unknown. Sponges (Porifera) provide an opportunity to reconstruct the earliest periods of myocyte evolution. Although sponges are believed to lack muscle, they are capable of coordinated whole-body contractions that purge debris from internal water canals. This behavior has been observed for decades, but their contractile tissues remain uncharacterized; it is an open question whether they have affinity to muscle or non-muscle contractile tissues in other animals. Here, we characterize the endothelial-like lining of water canals (the endopinacoderm) as a primary contractile tissue in the sponge Ephydatia muelleri. We find tissue-wide organization of contractile actin-bundles that contain striated-muscle myosin II and transgelin, and that contractions are regulated by the release of internal Ca2+ stores upstream of the myosin-light-chain-kinase (MLCK) pathway. Further, we show that the endopinacoderm is developmentally specified by myocardin-related transcription factor (MRTF) as part of an environmentally-inducible transcriptional complex that otherwise is known to function in muscle development, plasticity, and regeneration in other animals. We conclude that both muscle tissues and the endopinacoderm evolved from a common, multifunctional contractile tissue in the animal stem-lineage. Furthermore, as an actin-regulated force-sensor, MRTF-activity offers a mechanism for how water canals dynamically remodel in response to flow and can re-form normally from stem-cells in the absence of the intrinsic positional cues characteristic of embryogenesis in other animals.
Positive Correlation between Bacterial GC Content and Growth Temperature BackgroundGC pairs are generally more stable than AT pairs; GC-rich genomes were proposed to be more adapted to high temperatures than AT-rich genomes. Previous studies consistently showed positive correlations between growth temperature and the GC contents of structural RNA genes. However, for the whole genome sequences and the silent sites of the codons in protein-coding genes, the relationship between GC content and growth temperature is in a long-lasting debate.

ResultsWith a dataset much larger than previous studies (681 bacteria and 155 archaea with completely assembled genomes), our phylogenetic comparative analyses showed positive correlations between optimal growth temperature (Topt) and GC content both in bacterial and archaeal structural RNA genes and in bacterial whole genome sequences, chromosomal sequences, plasmid sequences, core genes, and accessory genes. However, in the 155 archaea, we did not observe a significant positive correlation of Topt with whole-genome GC content (GCw) or GC content at four-fold degenerate sites. We randomly drew 155 samples from the 681 bacteria for 1000 rounds. In most cases (> 95%), the positive correlations between Topt and genomic GC contents became statistically nonsignificant (P > 0.05). This result suggested that the small sample sizes might account for the lack of positive correlations between growth temperature and genomic GC content in the 155 archaea and the bacterial samples of previous studies. Comparing the GC content among four categories (psychrophiles/psychrotrophiles, mesophiles, thermophiles, and hyperthermophiles) also revealed a positive correlation between GCw and growth temperature in bacteria. By including the GCw of incompletely assembled genomes, we expanded the sample size of archaea to 303. Positive correlations between GCw and Topt appear especially after excluding the halophilic archaea whose GC contents might be are strongly shaped by intense UV radiation.

ConclusionsThis study explains the previous contradictory observations and ends a long debate. Prokaryotes growing in high temperatures have higher GC contents. Thermal adaptation is one possible explanation for the positive association. Meanwhile, we propose that the elevated efficiency of DNA repair in response to heat mutagenesis might have the by-product of increasing GC content like that happens in intracellular symbionts and marine bacterioplankton.
A flux-based machine learning model to simulate the impact of pathogen metabolic heterogeneity on drug interactions AO_SCPLOWBSTRACTC_SCPLOWDrug combinations are a promising strategy to counter antibiotic resistance. However, current experimental and computational approaches do not account for the entire complexity involved in combination therapy design, such as the effect of pathogen metabolic heterogeneity, changes in the growth environment, drug treatment order and time interval. To address these limitations, we present a comprehensive approach that uses genome-scale metabolic modeling and machine learning to guide combination therapy design. Our mechanistic approach (a) accommodates diverse data types, (b) accounts for time- and order-specific interactions, and (c) accurately predicts drug interactions in various growth conditions and their robustness to pathogen metabolic heterogeneity. Our approach achieved high accuracy (AUC = 0.83 for synergy, AUC = 0.98 for antagonism) in predicting drug interactions for E. coli cultured in 57 metabolic conditions based on experimental validation. The entropy in bacterial metabolic response was predictive of combination therapy outcomes across time scales and growth conditions. Simulation of metabolic heterogeneity using population-FBA identified two sub-populations of E. coli cells defined by regulation of cis- to trans-isomerization of unsaturated fatty acids with distinct effects on drug interactions. Analysis of the vast landscape of condition-specific drug interactions revealed a set of 24 robustly synergistic drug combinations with potential for clinical use.

SO_SCPLOWIGNIFICANCEC_SCPLOWWorldwide, 700,000 people die each year from drug-resistant infections. Drug combinations have great potential to reduce the spread of drug-resistant bacteria. However, their potency is impacted by both the pathogen growth environment and the heterogeneity in pathogen metabolism. The metabolic heterogeneity in a pathogen population allows them to survive antibiotic treatment. Here we present a flexible machine-learning framework that utilizes diverse data types to effectively search through the large design space of both sequential and simultaneous combination therapies across hundreds of simulated growth conditions and pathogen metabolic states. Our approach can serve as a useful guide for the selection of robustly synergistic drug combinations. Vapor exposure to Δ9-tetrahydrocannabinol (THC) slows locomotion of the Maine Lobster (Homarus americanus) RationaleDespite a long history of use in synaptic physiology, the lobster has been a neglected model for behavioral pharmacology. A restaurateur proposed that exposing lobster to cannabis smoke reduces anxiety and pain during the cooking process. It is unknown if lobster gill respiration in air would result in significant {Delta}9-tetrahydrocannabinol (THC) uptake and whether this would have any detectable behavioral effects.

ObjectiveThe primary goal was to determine tissue THC levels in the lobster after exposure to THC vapor. Secondary goals were to determine if THC vapor altered locomotor behavior or nociception.

MethodsTissue samples were collected (including muscle, brain and hemolymph) from Homarus americanus (N=3 per group) following 30 or 60 minutes of exposure to vapor generated by an e-cigarette device using THC (100 mg/mL in a propylene glycol vehicle). Separate experiments assessed locomotor behavior and hot water nociceptive responses following THC vapor exposure.

ResultsTHC vapor produced duration-related THC levels in all tissues examined. Locomotor activity was decreased (distance, speed, time-mobile) by 30 min inhalation of THC. Lobsters exhibit a temperature-dependent withdrawal response to immersion of tail, antennae or claws in warm water; this is novel evidence of thermal nociception for this species. THC exposure for 60 minutes had only marginal effect on nociception under the conditions assessed.

ConclusionsVapor exposure of lobsters, using an e-cigarette based model, produces dose-dependent THC levels in all tissues and reduces locomotor activity. Hot water nociception was temperature dependent, but only minimal anti-nociceptive effect of THC exposure was confirmed. TRF2 rescues telomere attrition and prolongs cell survival in Duchenne muscular dystrophy cardiomyocytes derived from human iPSCs Duchenne muscular dystrophy (DMD) is a severe muscle wasting disease caused by the lack of dystrophin. Heart failure, driven by cardiomyocyte death, fibrosis, and the development of dilated cardiomyopathy, is the leading cause of death in DMD patients. Current treatments decrease the mechanical load on the heart, but do not address the root cause of dilated cardiomyopathy: cardiomyocyte death. Previously, we showed that telomere shortening is a hallmark of cardiomyocytes in DMD patient cardiac tissues and in a mouse model of DMD we generated that manifests dilated cardiomyopathy. We also found that this telomere shortening is recapitulated in DMD cardiomyocytes differentiated from patient-derived induced pluripotent stem cells (iPSCs). Here we further characterize that DMD cardiomyocytes exhibit reduced cell size, nuclear size, and sarcomere density compared to healthy isogenic controls. The telomere-binding protein, TRF2, is a core component of the shelterin complex, which protects chromosome ends and is expressed at lower levels in DMD cardiomyocytes compared to controls. We investigated whether preservation of telomere length could be achieved with upregulation of TRF2. We found that TRF2 upregulation in DMD cardiomyocytes increased telomere lengths, cell size, nuclear size, sarcomere density, and cell survival. These data suggest TRF2 gene therapy has the potential to delay the onset of dilated cardiomyopathy. Plant host shift: The viability and persistence of a polyphagous, non-native insect Assessing the effects of a plant-host shift is important for monitoring insect populations over long time periods and for interventions in a conservation or pest management framework. In a heterogeneous environment, individuals may disperse between sources and sinks in order to persist. Here we propose a single-species two-patch model that aims to capture the generational movement of an insect that exhibits density-dependent dispersal, to see how shifting between hosts could alter its viability and asymptotic dynamics. We then analyse the stability and persistence properties of the model and further validate it using parameter estimates derived from laboratory experiments. In order to evaluate the potential of this model, we applied it to Drosophila suzukii (Diptera: Drosophilidae), which has become a harmful pest in several countries around the world. Although many studies have investigated the preference and attractiveness of potential hosts on this invasive drosophilid, no studies thus far have investigated whether a shift of fruit host could affect such a species ecological viability or spatiotemporal persistence. The model results show that a shift in host choice can significantly affect the growth potential and fecundity of a species such as D. suzukii, which ultimately could aid such invasive populations in their ability to persist within a changing environment. Naive human B cells engage the receptor binding domain of SARS-CoV-2, variants of concern, and related sarbecoviruses Exposure to a pathogen elicits an adaptive immune response aimed to control and eradicate. Interrogating the abundance and specificity of the naive B cell repertoire contributes to understanding how to potentially elicit protective responses. Here, we isolated naive B cells from 8 seronegative human donors targeting the SARS-CoV-2 receptor-binding domain (RBD). Single B cell analysis showed diverse gene usage with no restricted complementarity determining region lengths. We show that recombinant antibodies engage SARS-CoV-2 RBD, circulating variants, and pre-emergent coronaviruses. Representative antibodies signal in a B cell activation assay and can be affinity matured through directed evolution. Structural analysis of a naive antibody in complex with spike shows a conserved mode of recognition shared with infection-induced antibodies. Lastly, both naive and affinity-matured antibodies can neutralize SARS-CoV-2. Understanding the naive repertoire may inform potential responses recognizing variants or emerging coronaviruses enabling the development of pan-coronavirus vaccines aimed at engaging germline responses.

One Sentence SummaryIsolation of antibody germline precursors targeting the receptor binding domain of coronaviruses. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. IRGM/Irgm1 Aggravates Progression of Atherosclerosis by Inducing Macrophage Apoptosis through the MAPK Signaling Pathway AimsAtherosclerosis underlies most cardiovascular diseases, among which acute coronary syndrome (ACS) caused by plaque rupture (PR) often leads to death. Immune-related GTPases (IRGM/Irgm1) have been extensively studied in inflammatory diseases, but their role in atherosclerosis is unclear. Determining how IRGM/Irgm1 promotes atherosclerotic plaque vulnerability will provide information for new biomarkers and/or therapeutic targets.

Methods and resultsWe identified ruptured and unruptured plaques by optical coherence tomography, and found that serum IRGM was highly expressed in patients with ST-segment elevation myocardial infarction in PR. We used ApoE-/-Irgm1+/+, ApoE-/-Irgm1+/- mice and chimeric mice to establish a model of advanced atherosclerosis. The results of pathological experiments showed that Irgm1 caused plaque necrosis. The ratio of neutral lipids and cholesterol crystals increases, while the content of collagen fibers decreases, aggravating the destabilization of atherosclerotic plaques. In vitro, we used multiple approaches to confirm that Irgm1 promotes macrophage apoptosis by promoting the production of reactive oxygen species and activating the MAPK signaling pathway.

ConclusionsIRGM may be a potential risk factor for PR. Mechanistic studies have shown that IRGM/Irgm1 contributes to the formation and rupture of fragile plaques. This is partly mediated by the induction of macrophage apoptosis via the MAPK signaling pathway. IRGM may offer new strategies for early treatment of ACS.

Translation viewOur findings indicate that IRGM/Irgm1 contributes to formation and rupture of vulnerable plaques. It suggests that IRGM may provide a new target for the early treatment of ACS. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. Analysis of Drosophila cardiac hypertrophy by micro-computerized tomography for genetic dissection of heart growth mechanisms Heart failure is often preceded by pathological cardiac hypertrophy, a thickening of the heart musculature driven by complex gene regulatory and signaling processes. The Drosophila heart has great potential as a genetic model for deciphering the underlying mechanisms of cardiac hypertrophy. However, current methods for evaluating hypertrophy of the Drosophila heart are laborious and difficult to carry out reproducibly. Here we demonstrate that micro-computerized tomography (microCT) is an accessible, highly reproducible method for non-destructive, quantitative analysis of Drosophila heart morphology and size. To validate our microCT approach for analyzing Drosophila cardiac hypertrophy, we show that expression of constitutively active Ras (Ras85DV12), previously shown to cause hypertrophy of the fly heart, results in significant thickening of both adult and larval heart walls when measured from microCT images. We then show using microCT analysis that genetic upregulation of store-operated Ca2+ entry (SOCE) driven by expression of constitutively active Stim (StimCA) or Orai (OraiCA) proteins also results in significant hypertrophy of the Drosophila heart, through a process that specifically depends on Orai Ca2+ influx channels. Intravital imaging of heart contractility revealed significantly reduced end diastolic dimensions in StimCA and OraiCA expressing hearts, consistent with the hypertrophic phenotype. These results demonstrate that increased SOCE activity is an important driver of hypertrophic cardiomyocyte growth, and demonstrate how microCT analysis combined with tractable genetic tools in Drosophila can be used to delineate molecular signaling processes that underlie cardiac hypertrophy and heart failure.

NEW AND NOTEWORTHYGenetic analysis of cardiac hypertrophy in Drosophila holds immense potential for the discovery of new therapeutic targets to prevent and treat heart failure. However, this potential has been hindered by a lack of rapid and effective methods for analysis of heart size in flies. Here we demonstrate that analysis of the Drosophila heart with micro-computerized tomography yields accurate and highly reproducible heart size measurements that can be used to efficiently analyze heart growth and cardiac hypertrophy in Drosophila. Analysis of alternative mRNA splicing in vemurafenib-resistant melanoma cells Alternative mRNA splicing is common in cancers. In BRAF V600E mutated malignant melanoma a frequent mechanism of acquired resistance to BRAF inhibitors involves alternative splicing (AS) of BRAF. The resulting shortened BRAF protein constitutively dimerizes and conveys drug resistance. Here, we have analysed AS in SKMEL-239 melanoma cells and a BRAF inhibitor (vemurafenib) resistant derivative that expresses an AS, shortened BRAF V600E transcript. Transcriptome analysis showed differential expression of spliceosome components between the two cell lines. As there is no consensus approach to analysing AS events, we used and compared four common AS softwares based on different principles, DEXSeq, rMATS, ASpli, and LeafCutter. Two of them correctly identified the BRAF V600E AS in the vemurafenib resistant cells. Only 12 AS events were identified by all four softwares. Testing the AS predictions experimentally showed that these overlapping predictions are highly accurate. Interestingly, they identified AS caused alterations in the expression of melanin synthesis and cell migration genes in the vemurafenib resistant cells. This analysis shows that combining different AS analysis approaches produce reliable results and meaningful, biologically testable hypotheses.

Graphic Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=70 SRC="FIGDIR/small/484656v2_ufig1.gif" ALT="Figure 1">
View larger version (24K):
org.highwire.dtl.DTLVardef@9f4920org.highwire.dtl.DTLVardef@179e9bcorg.highwire.dtl.DTLVardef@e05632org.highwire.dtl.DTLVardef@3b1c41_HPS_FORMAT_FIGEXP  M_FIG C_FIG Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

Here, we designed, prepared, and characterized two combinatorial protein libraries consisting of randomized proteins, each 105 residues in length. The first library constructed proteins from the entire canonical alphabet of 20 amino acids. The second library used a subset of only 10 residues (A,S,D,G,L,I,P,T,E,V) that represent a consensus view of plausibly available amino acids through prebiotic chemistry. Our study shows that compact conformations resistant to proteolysis are (i) abundant (up to 40%) in random sequence space, (ii) independent of general Hsp70 chaperone system activity, and (iii) not granted solely by "late" and complex amino acid additions. The Hsp70 chaperone system effectively increases solubility and refoldability of the canonical alphabet but has only a minor impact on the "early" library. The early alphabet proteins are inherently more soluble and refoldable, possibly assisted by the cell-like environment in which these assays were performed.

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. Evaluation of methods for estimating coalescence times using ancestral recombination graphs The ancestral recombination graph (ARG) is a structure that describes the joint genealogies of sampled DNA sequences along the genome. Recent computational methods have made impressive progress towards scalably estimating whole-genome genealogies. In addition to inferring the ARG, some of these methods can also provide ARGs sampled from a defined posterior distribution. Obtaining good samples of ARGs is crucial for quantifying statistical uncertainty and for estimating population genetic parameters such as effective population size, mutation rate, and allele age. Here, we use standard neutral coalescent simulations to benchmark the estimates of pairwise coalescence times from three popular ARG inference programs: ARGweaver, Relate, and tsinfer+tsdate. We compare 1) the true coalescence times to the inferred times at each locus; 2) the distribution of coalescence times across all loci to the expected exponential distribution; 3) whether the sampled coalescence times have the properties expected of a valid posterior distribution. We find that inferred coalescence times at each locus are most accurate in ARGweaver, and often more accurate in Relate than in tsinfer+tsdate. However, all three methods tend to overestimate small coalescence times and underestimate large ones. Lastly, the posterior distribution of ARGweaver is closer to the expected posterior distribution than Relates, but this higher accuracy comes at a substantial trade-off in scalability. The best choice of method will depend on the number and length of input sequences and on the goal of downstream analyses, and we provide guidelines for the best practices. Differed cell division angle, position of cell proliferative area, and localization of ANGUSTIFOLIA3 in lateral organs Leaf meristem is a cell proliferative zone present in the lateral organ primordia. In this study, we investigated how the proliferative zone affects the final morphology of the lateral organs. We examined how cell proliferative zones differ in the primordia of planar floral organs and polar auxin transport inhibitor (PATI)-treated leaves from normal foliage leaf primordia of Arabidopsis thaliana with a focus on the spatial accumulation pattern of ANGUSTIFOLIA3 (AN3), a key element for leaf meristem positioning. We found that organ shape changes by PATI treatment were correlated to the angle of the cell division plane relative to the leaf primordia axis in the leaf meristem (cell division angle), but not with leaf-meristem positioning, size of the leaf meristem, or the localization pattern of AN3 protein. In contrast, different shapes between sepals and petals compared with foliage leaves were associated with both altered meristem position associated with altered AN3 expression patterns and different distributions of cell division angles. These results suggest that lateral organ shapes are regulated via two aspects: position of meristem and cell division angles

Summary statementDifferent lateral organs with different morphology possess different properties of meristems; cell division angles, position of cell proliferative area and AN3 localization patterns. Frequency modulated timer regulates mammalian hibernation Mammalian hibernators decrease basal metabolism and body temperature (Tb) to minimize energy expenditure in harsh seasons. During hibernation, Tb drops to low temperature (<10 {degrees}C) and remains constant for days, known as deep torpor in small mammalian hibernators. Spontaneous interbout arousals interrupt torpor bouts, when Tb recovers to euthermic state ~37 {degrees}C. Torpor-interbout arousal event repeats during several months of hibernation. However, little is known about mechanisms governing Tb fluctuation across torpor-interbout arousal cycles during hibernation. Recent improvement in data-logging techniques enables us to monitor Tb for more than hundred days with high precision, opening up new avenues for quantitative analysis to address the principle governing Tb fluctuation. Here, we analyzed Tb fluctuation across torpor-interbout arousal cycle of Syrian hamster, which can hibernate in response to chronic cold and short photoperiod under a laboratory condition, using generalized harmonic analysis and discovered a model with frequency modulation quantitatively reproducing Tb fluctuation. This analysis also identified an unexpectedly longer period of 120-430 days as the period that modulates another period of several days, generating Tb fluctuation for Syrian hamster. We propose that concerted action of two endogenous periods governs torpor-interbout arousal cycles during hibernation. Notch signaling via Hey1 and Id2b regulates Muller glia's regenerative response to retinal injury Zebrafish Muller glia (MG) respond to retinal injury by suppressing Notch signaling and producing progenitors for retinal repair. A certain threshold of injury-derived signal must be exceeded in order to engage MG in a regenerative response (MGs injury-response threshold). Pan-retinal Notch inhibition expands the zone of injury-responsive MG at the site of focal injury, suggesting that Notch signaling regulates MGs injury-response threshold. We found that Notch signaling enhanced chromatin accessibility and gene expression at a subset of regeneration-associated genes in the uninjured retina. Two Notch effector genes, hey1 and id2b, were identified that reflect bifurcation of the Notch signaling pathway, and differentially regulate MGs injury-response threshold and proliferation of MG-derived progenitors. Furthermore, Notch signaling component gene repression in the injured retina suggests a role for Dll4, Dlb, and Notch3 in regulating Notch signaling in MG and epistasis experiments confirm that the Dll4/Dlb-Notch3-Hey1/Id2b signaling pathway regulates MGs injury-response threshold and proliferation. A high throughput lipidomics method using scheduled multiple reaction monitoring Lipid compositions of cells, tissues and bio-fluids are complex, with varying concentrations and structural diversity, which makes their identification challenging. Newer methods for comprehensive analysis of lipids are thus necessary. Herein, we propose a targeted-mass spectrometry based method for large-scale lipidomics using a combination of variable retention time window and relative dwell time weightage. Using this, we detected more than 1000 lipid species, including structural isomers. The limit of detection varied from femtomolar to nanomolar range and the coefficient of variance <30% for 849 lipid species. We used this method to identify lipids altered due to Vitamin B12 deficiency and found that the levels of lipids with {omega}-3 fatty acid chains decreased while those with {omega}-6 increased. This method enables identification of by far the largest number of lipid species with structural isomers in a single experiment and would significantly advance our understanding of the role of lipids in biological processes. Post-translational flavinylation is associated with diverse extracytosolic redox functionalities throughout bacterial life Disparate redox activities that take place beyond the bounds of the prokaryotic cell cytosol must connect to membrane or cytosolic electron pools. Proteins post-translationally flavinylated by the enzyme ApbE mediate electron transfer in several characterized extracytosolic redox systems but the breadth of functions of this modification remains unknown. Here we present a comprehensive bioinformatic analysis of 31,910 prokaryotic genomes that provides evidence of extracytosolic ApbEs within ~50% of bacteria and the involvement of flavinylation in numerous uncharacterized biochemical processes. By mining flavinylation-associated gene clusters, we identify five protein classes responsible for transmembrane electron transfer and two domains of unknown function (DUF2271 and DUF3570) that are flavinylated by ApbE. We observe flavinylation/iron transporter gene colocalization patterns that implicate functions in iron reduction and assimilation. We find associations with characterized and uncharacterized respiratory oxidoreductases that highlight roles of flavinylation in respiratory electron transport chains. Finally, we identify interspecies gene cluster variability consistent with flavinylation/cytochrome functional redundancies and discover a class of "multi-flavinylated proteins'' that may resemble multiheme cytochromes in facilitating longer distance electron transfer. These findings provide key mechanistic insight into an important facet of bacterial physiology and establish flavinylation as a functionally diverse mediator of extracytosolic electron transfer.

Impact StatementBioinformatic and biochemical studies find that covalently bound flavins are common and participate in wide-ranging extracytosolic redox activities throughout bacterial life. A pilot study on the home range and movement patterns of the Andean Fox Lycalopex culpaeus (Molina 1782) in Cotopaxi National Park, Ecuador This study reports movement patterns and home range estimates of an Andean fox (Lycalopex culpaeus) in Cotopaxi National Park in Ecuador, representing the first GPS-tagging of the species. The GPS functioned well during the 197-day tracking period. Home range sizes ranged between 4.9 - 8.1 km2, depending on the estimation method. Movement speeds averaged 0.17 km/hr at day vs. 0.87 km/hr. at night, and distance traveled averaged 0.23 km at day vs. 0.89 km at night. These preliminary results highlight the importance of collecting unbiased, high-quality data which enables an enhanced understanding on mammal behavior and human/animal interaction. The integrity and assay performance of tissue mitochondrial DNA is considerably affected by choice of isolation method The integrity of mitochondrial DNA (mtDNA) isolated from solid tissues is critical for analyses such as long-range PCR, but is typically assessed under conditions that fail to provide information on the individual mtDNA strands. Using denaturing gel electrophoresis, we show that commonly-used isolation procedures generate mtDNA containing several single-strand breaks per strand. Through systematic comparison of DNA isolation methods, we identify a procedure yielding the highest integrity of mtDNA that we demonstrate displays improved performance in downstream assays. Our results highlight the importance of isolation method choice, and serve as a resource to researchers requiring high-quality mtDNA from solid tissues. Identifying complex motifs in massive omics data with a variable-convolutional layer in deep neural network Motif identification is among the most common and essential computational tasks for bioinformatics and genomics. Here we proposed a novel convolutional layer for deep neural network, named Variable Convolutional (vConv) layer, for effective motif identification in high-throughput omics data by learning kernel length from data adaptively. Empirical evaluations on DNA-protein binding and DNase footprinting cases well demonstrated that vConv-based networks have superior performance to their convolutional counterparts regardless of model complexity. Meanwhile, vConv could be readily integrated into multi-layer neural networks as an "in-place replacement" of canonical convolutional layer. All source codes are freely available on GitHub for academic usage. A novel method for producing functionalized vesicles that efficiently deliver oligonucleotides in vitro in cancer cells and in vivo in mice Nano-based delivery systems have enhanced our ability to administer and target drugs and macromolecules to their targets. Oligonucleotide drugs have great therapeutic potential but often have off-target effects and stability issues. Therefore, they are often encapsulated in vesicles with targeting ligands such as antibodies (Ab) to deliver their cargo. Herein, we describe a novel, scalable and straightforward approach to producing functionalized vesicles called the "Functionalized Lipid Insertion Method," which differs from the older "Detergent-Dialysis Method." The latter method required excess detergent and extensive dialysis over many hours to produce the functionalized vesicles. With our method, only the functionalized lipid is detergent-solubilized during the engineering of the vesicle. The approach reduces the dialysis time, keeps the vesicles intact while orienting the targeting moieties of the functionalized lipid toward the outside of the vesicle. Pilot in vitro and in vivo experiments was performed to show the feasibility of our method. Dynamic light scattering (DLS) experiments suggested that the original vesicular structure was relatively unperturbed, and the functionalized lipid was inserted externally. Our approach efficiently delivered oligonucleotides and affected the function of liver cancer HepG2 cells. Furthermore, functionalized vesicles achieved targeted delivery of oligonucleotides in mice without inducing a significant innate immune response. The industrial and therapeutic significance and implications of functionalized vesicles produced by our method are also discussed. Additional experiments and analyses are recommended to bring out the full potential of this molecular delivery technology. Platform for Remote, Multifocal Brain Interventions Transcranial focused ultrasound brings precision healthcare to the human brain, delivering focal interventions to deep brain sites without the need for surgeries or systemic drug administrations. The approaches have included neuromodulation and localized delivery of drugs through an intact or temporarily opened blood brain barrier. Clinical translation of these approaches requires a system that can deliver ultrasound into specific brain targets on command. Here, we demonstrate a platform with these features, and validate its function in nonhuman primates. Specifically, we used the system to flexibly modulate two deep brain nuclei--the left and right lateral geniculate nucleus-during visual choice behavior. We found that ultrasound induced a behavioral bias whose polarity depended on the stimulation waveform. Moreover, repeated applications of ultrasound over 39 consecutive behavioral sessions produced no signs of harm to the stimulated circuits. The system can be be applied to diagnose and possibly treat disorders of brain function with a precision and flexibility previously unattainable. ConVarT: a search engine for matching human genetic variants with variants from non-human species The availability of genetic variants, together with phenotypic annotations from model organisms, facilitates comparing these variants with equivalent variants in humans. However, existing databases and search tools do not make it easy to scan for equivalent variants, namely "matching variants" (MatchVars) between humans and other organisms. Therefore, we developed an integrated search engine called ConVarT (http://www.convart.org/) for matching variants between humans, mice, and C. elegans. ConVarT incorporates annotations (including phenotypic and pathogenic) into variants, and these previously unexploited phenotypic MatchVars from mice and C. elegans can give clues about the functional consequence of human genetic variants. Our analysis shows that many phenotypic variants in different genes from mice and C. elegans, so far, have no counterparts in humans, and thus, can be useful resources when evaluating a relationship between a new human mutation and a disease. Reliable  in silico  ranking of engineered therapeutic TCR binding affinities with MMPBSA and MMGBSA. Accurate and efficient in silico ranking of protein-protein binding affinities is useful for protein design with applications in biological therapeutics. One popular approach to rank binding affinities is to apply the molecular mechanics Poisson Boltzmann/generalized Born surface area (MMPB/GBSA) method to molecular dynamics trajectories. Here, we identify protocols that enable the reliable evaluation of T-cell receptor (TCR) variants binding to their target, peptide-human leukocyte antigens (pHLAs). We suggest different protocols for variant sets with few ([&le;]4) or many mutations, with entropy corrections important for the latter. We demonstrate how potential outliers could be identified in advance and that just 5-10 replicas of short (4 ns) MD simulations may be sufficient for reproducible and accurate ranking of TCR variants. The protocols developed here can be applied towards in silico screening during the optimization of therapeutic TCRs, potentially reducing both the cost and time taken for biologic development.