MetaLogo: a heterogeneity-aware sequence logo generator and aligner Sequence logos are used to visually display conservations and variations in short sequences. They can indicate the fixed patterns or conserved motifs in a batch of DNA or protein sequences. However, most of the popular sequence logo generators are based on the assumption that all the input sequences are from the same homologous group, which will lead to an overlook of the heterogeneity among the sequences during the sequence logo making process. Heterogeneous groups of sequences may represent clades of different evolutionary origins, or genes families with different functions. Therefore, it is essential to divide the sequences into different phylogenetic or functional groups to reveal their specific sequence motifs and conservation patterns. To solve these problems, we developed MetaLogo, which can automatically cluster the input sequences after multiple sequence alignment and phylogenetic tree construction, and then output sequence logos for multiple groups and aligned them in one figure. User-defined grouping is also supported by MetaLogo to allow users to investigate functional motifs in a more delicate and dynamic perspective. MetaLogo can highlight both the homologous and nonhomologous sites among sequences. MetaLogo can also be used to annotate the evolutionary positions and gene functions of unknown sequences, together with their local sequence characteristics. We provide users a public MetaLogo web server (http://metalogo.omicsnet.org), a standalone Python package (https://github.com/labomics/MetaLogo), and also a built-in web server available for local deployment. Using MetaLogo, users can draw informative, customized and publishable sequence logos without any programming experience to present and investigate new knowledge on specific sequence sets.
A novel, computationally tractable algorithm flags in big matrices every column associated in any way with others or a dependent variable, with much higher power when columns are linked like mutations in chromosomes When a data matrix DM has many independent variables IVs, it is not computationally tractable to assess the association of every distinct IV subset with the dependent variable DV of the DM, because the number of subsets explodes combinatorially as IVs increase. But model selection and correcting for multiple tests is complex even with few IVs.

DMs in genomics will soon summarize millions of mutation markers and genomes. Searching exhaustively in such DMs for markers that alone or synergistically with others are associated with a trait is therefore computationally tractable only for 1- and 2-marker effects. Also population geneticists study mainly 2-marker combinations.

I present a computationally tractable, fully parallelizable Participation in Association Score (PAS) that in a DM with markers detects one by one every column that is strongly associated in any way with others. PAS does not examine column subsets and its computational cost grows linearly with the number of columns, remaining reasonable even when DMs have millions of columns.

PAS exploits how associations of markers in the rows of a DM cause associations of matches in the rows pairwise comparisons. For every such comparison with a match at a tested column, PAS computes the matches at other columns by modifying the comparisons total matches (scored once per DM), yielding a distribution of conditional matches that reacts diagnostically to the associations of the tested column. Equally computationally tractable is dvPAS that flags DV-associated IVs by also probing the matches at the DV.

P values for the scores are readily obtained by permutation and accurately Sidak-corrected for multiple tests, bypassing model selection. The P values of a columns PASs for different orders of association are i.i.d. and readily turned into a single P value.

Simulations show that i) PAS and dvPAS generate uniform-(0,1)-distributed type I error in null DMs and ii) detect randomly encountered binary and trinary models of significant n-column association and n-IV association with a binary DV, respectively, with power in the order of magnitude of exhaustive evaluations and false positives that are uniform-(0,1)-distributed or straightforwardly tuned to be so. Power to detect 2-way associations that extend over 100+ columns is non-parametrically ultimate but that to detect pure n-column associations and pure n-IV DV associations sinks exponentially as n increases.

Important for geneticists, dvPAS power increases about twofold in trinary vs. binary DMs and by orders of magnitude with markers linked like mutations in chromosomes, specially in trinary DMs where furthermore dvPAS fine-maps with highest resolution.
CHRM1 is a Druggable Melanoma Target Whose Endogenous Activity is Determined by Inherited Genetic Variation in DOPA Production Melanoma risk is 30 times higher in people with lightly pigmented skin compared to those with darkly pigmented skin. Here we show that this difference results from more than melanin pigment and its ultraviolet radiation (UVR) shielding effect. Using primary human melanocytes representing the full human skin pigment continuum and several preclinical melanoma models, we show that cell-intrinsic differences between dark and light melanocytes regulate melanocyte proliferative capacity, overall cellular differentiation state, and susceptibility to malignant transformation, independently of melanin and UV exposure. We determined that these differences result from dihydroxyphenylalanine (DOPA), a melanin precursor synthesized at higher levels in melanocytes from dark skin. Although DOPA was not previously known to have specific signaling activity, we used both high throughput pharmacologic and genetic in vivo CRISPR screens to determine that DOPA limits melanocyte and melanoma cell proliferation by directly inhibiting the muscarinic acetylcholine receptor M1 (CHRM1), a G Protein-Coupled Receptor (GPCR) not previously known to bind DOPA, nor to affect melanoma pathobiology. Pharmacologic CHRM1 antagonism in melanoma leads to depletion of c-Myc and FOXM1, both of which are proliferation drivers associated with aggressive melanoma. In preclinical mouse melanoma models using both immune deficient and syngeneic immune competent mice, pharmacologic inhibition of CHRM1 or FOXM1 inhibited tumor growth. CHRM1 and FOXM1 may be new therapeutic targets for melanoma. Modular biosynthesis of plant hemicellulose and its impact on yeast cells BackgroundThe carbohydrate polymers that encapsulate plants cells have benefited humans for centuries and have valuable biotechnological uses. In the past five years, exciting possibilities have emerged in the engineering of polysaccharide-based biomaterials. Despite impressive advances on bacterial cellulose-based hydrogels, comparatively little is known about how plant hemicelluloses can be reconstituted and modulated in cells suitable for biotechnological purposes.

ResultsHere, we assembled cellulose synthase-like A (CSLA) enzymes using an optimized Pichia pastoris platform to produce tunable heteromannan (HM) polysaccharides in yeast. By swapping the domains of plant mannan and glucomannan synthases, we engineered chimeric CSLA proteins that made {beta}-1,4-linked mannan in quantities surpassing those of the native enzymes while minimizing the burden on yeast growth. Prolonged expression of a glucomannan synthase from Amorphophallus konjac was toxic to yeast cells: reducing biomass accumulation and ultimately leading to compromised cell viability. However, an engineered glucomannan synthase as well as CSLA pure mannan synthases and a CSLC glucan synthase did not inhibit growth. Interestingly, Pichia cell size could be increased or decreased depending on the composition of the CSLA protein sequence. HM yield and glucose incorporation could be further increased by co-expressing chimeric CSLA proteins with a MANNAN-SYNTHESIS-RELATED (MSR) co-factor from Arabidopsis thaliana.

ConclusionThe results provide novel routes for the engineering of polysaccharide-based biomaterials that are needed for a sustainable bioeconomy. The characterization of chimeric cellulose synthase-like enzymes in yeast offers an exciting avenue to produce plant polysaccharides in a tunable manner. Furthermore, cells modified with non-toxic plant polysaccharides such as {beta}-mannan offer a modular chassis to produce and encapsulate sensitive cargo such as therapeutic proteins. Heat-Dependent Hairpin Melting Drives Thermo-gated TRPV1 Opening The capsaicin receptor TRPV1 is a bio-thermometer in a primary afferent sensory neuron for noxious heat detection. However, the underlying molecular mechanism is unclear. Here, a hairpin thermodynamic model, together with graph theory, was to examine a role of temperature-dependent hairpin sizes in controlling the temperature threshold and sensitivity and gating states of TRPV1. The results showed that heat-dependent melting of the largest hairpin initiates a temperature threshold for channel opening while the smaller hairpins stabilize the open state and increase the heat efficacy or sensitivity. Therefore, this model may provide a broad structural basis for the thermo-gated TRP channels. SufB intein splicing in Mycobacterium tuberculosis is influenced by two remote conserved N-extein Histidines Inteins are auto-processing domains that implement a multi-step biochemical reaction termed protein splicing, marked by cleavage and formation of peptide bonds. They excise from a precursor protein, generating a functional protein via covalent bonding of flanking exteins. We report the kinetic study of splicing and cleavage reaction in [Fe-S] cluster assembly protein SufB from Mycobacterium tuberculosis. Although it follows a canonical intein splicing pathway, distinct features are added by extein residues present in the active site. Sequence analysis identified two conserved histidines in the N-extein region; His-5 and His-38. Kinetic analyses of His-5Ala and His-38Ala SufB mutants exhibited significant reductions in splicing and cleavage rates relative to the SufB wild-type precursor protein. Structural analysis and molecular dynamics simulations suggested that Mtu SufB displays a unique mechanism where two remote histidines work concurrently to facilitate N-terminal cleavage reaction. His-38 is stabilized by the solvent-exposed His-5, and can impact N-S acyl shift by direct interaction with the catalytic Cys1. Development of inteins as biotechnological tools or as pathogen specific novel antimicrobial targets requires a more complete understanding of such unexpected roles of conserved extein residues in protein splicing. Spatial and Texture Analysis of Root System Distribution with Earth Mover's Distance (STARSEED) PurposeRoot system architectures are complex, multidimensional, and challenging to characterize effectively for agronomic and ecological discovery.

MethodsWe propose a new method, Spatial and Texture Analysis of Root SystEm distribution with Earth movers Distance (STARSEED), for comparing root system distributions that incorporates spatial information through a novel application of the Earth Movers Distance (EMD).

ResultsWe illustrate that the approach captures the response of sesame root systems for different genotypes and soil moisture levels. STARSEED provides quantitative and visual insights into changes that occur in root architectures across experimental treatments.

ConclusionSTARSEED can be generalized to other plants and provides insight into root system architecture development and response to varying growth conditions not captured by existing root architecture metrics and models. The code and data for our experiments are publicly available: https://github.com/GatorSense/STARSEED. Predicting the zoonotic capacity of mammals to transmit SARS-CoV-2 Back and forth transmission of SARS-CoV-2 between humans and animals may lead to wild reservoirs of virus that can endanger efforts toward long-term control of COVID-19 in people, and protecting vulnerable animal populations that are particularly susceptible to lethal disease. Predicting high risk host species is key to targeting field surveillance and lab experiments that validate host zoonotic potential. A major bottleneck to predicting animal hosts is the small number of species with available molecular information about the structure of ACE2, a key cellular receptor required for viral cell entry. We overcome this bottleneck by combining species ecological and biological traits with 3D modeling of virus and host cell protein interactions using machine learning methods. This approach enables predictions about the zoonotic capacity of SARS-CoV-2 for over 5,000 mammals -- an order of magnitude more species than previously possible. The high accuracy predictions achieved by this approach are strongly corroborated by in vivo empirical studies. We identify numerous common mammal species whose predicted zoonotic capacity and close proximity to humans may further enhance the risk of spillover and spillback transmission of SARS-CoV-2. Our results reveal high priority areas of geographic overlap between global COVID-19 hotspots and potential new mammal hosts of SARS-CoV-2. With molecular sequence data available for only a small fraction of potential host species, predictive modeling integrating data across multiple biological scales offers a conceptual advance that may expand our predictive capacity for zoonotic viruses with similarly unknown and potentially broad host ranges. 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. Lmo7 recruits myosin II heavy chain to induce apical constriction in Xenopus ectoderm Apical constriction, or reduction of the apical domain, underlies many morphogenetic events during development, such as furrow or tube formation. Actomyosin complexes play an essential role in apical constriction, however the detailed analysis of molecular mechanisms is still pending. Here we show that Lim domain only protein 7 (Lmo7), a multidomain adaptor at apical junctions, promotes apical constriction in the Xenopus superficial ectoderm, whereas apical domain size increases in Lmo7-depleted cells. Lmo7 is localized to and promotes the formation of the circumferential actomyosin belt adjacent to apical junctions. We find that Lmo7-dependent apical constriction requires the RhoA-ROCK-Non-muscle myosin II (NMII) pathway. Strikingly, Lmo7 binds and recruits NMII heavy chains to apical junctions. Lmo7 overexpression altered the subcellular distribution of Wtip, a sensor of mechanical tension. Our findings suggest that Lmo7 serves as a scaffold organizing the actomyosin network to regulate contractility at apical junctions. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Development and Characterization of PAR-Trackers: New Tools for Detecting Poly(ADP-ribose) In Vitro and In Vivo AO_SCPLOWBSTRACTC_SCPLOWADP-ribosylation (ADPRylation) is a reversible post-translation modification resulting in the covalent attachment of ADP-ribose (ADPR) moieties on substrate proteins. Naturally-occurring protein motifs and domains, including WWEs, PBZs, and macrodomains, act as "readers" for protein-linked ADPR. Although recombinant, antibody-like ADPR detection reagents containing these readers have facilitated the detection of ADPR, they are limited in their ability to capture the dynamic nature of ADPRylation. Herein, we describe and characterize a set of poly(ADP-ribose) (PAR) Trackers (PAR-Ts) - optimized dimerization-dependent or split-protein reassembly PAR sensors in which a naturally occurring PAR binding domain, WWE, was fused to both halves of dimerization-dependent GFP (ddGFP) or split Nano Luciferase (NanoLuc), respectively. We demonstrate that these new tools allow the detection and quantification of PAR levels in extracts, living cells, and living tissues with greater sensitivity, as well as temporal and spatial precision. Importantly, these sensors detect changes in cellular ADPR levels in response to physiological cues (e.g., hormone-dependent induction of adipogenesis without DNA damage), as well as xenograft tumor tissues in living mice. Our results indicate that PAR Trackers have broad utility for detecting ADPR in many different experimental and biological systems. Swine H1N1 influenza virus variants with enhanced polymerase activity and HA stability promote airborne transmission in ferrets Understanding how animal influenza A viruses (IAVs) acquire airborne transmissibility in humans and ferrets is needed to prepare for and respond to pandemics. Here, we investigated in ferrets the replication and transmission of swine H1N1 isolates P4 and G15, whose majority population had decreased polymerase activity and poor HA stability, respectively. For both isolates, a minor variant was selected and transmitted in ferrets. Polymerase-enhancing variant PA-S321 airborne-transmitted and propagated in one ferret. HA-stabilizing variant HA1-S210 was selected in all G15-inoculated ferrets and was transmitted by contact and airborne routes. With an efficient polymerase and a stable HA, the purified minor variant G15-HA1-S210 had earlier and higher peak titers in inoculated ferrets and was recovered at a higher frequency after airborne transmission than P4 and G15. Overall, HA stabilization played a more prominent role than polymerase enhancement in the replication and transmission of these viruses in ferrets. The results suggest pandemic risk-assessment studies may benefit from deep sequencing to identify minor variants with human-adapted traits.

IMPORTANCEDiverse IAVs circulate in animals, yet few acquire the viral traits needed to start a human pandemic. A stabilized HA and mammalian-adapted polymerase have been shown to promote the adaptation of IAVs to humans and ferrets (the gold-standard model for IAV replication, pathogenicity, and transmissibility). Here, we used swine IAV isolates of the gamma lineage a model to investigate the importance of HA stability and polymerase activity in promoting replication and transmission in ferrets. These are emerging viruses that bind to both -2,6- and -2,3-linked receptors. Using isolates containing mixed populations, a stabilized HA was selected within days in inoculated ferrets. An enhanced polymerase was also selected and propagated after airborne transmission to a ferret. Thus, HA stabilization was a stricter requirement, yet both traits promoted transmissibility. Knowing the viral traits needed for pandemic potential, and the relative importance of each, will help identify emerging viruses of greatest concern. Polycystins are upregulated in fibrotic kidneys and promote renal fibrosis Autosomal polycystic kidney disease (ADPKD) is caused by mutations in PKD1or PKD2 gene, which encodes polycystin-1 or polycystin-2 respectively. Paradoxically, polycystins are up-regulated in human ADPKD kidneys and overexpression of polycystin-1 even induced renal cyst formation in mice, suggesting that up-regulation of polycystin is pathogenic. In this study, we determined the role of polycystin in renal fibrosis and explored its underlying mechanisms. We found that the expression of polycystin-1 or polycystin-2 was positively correlated with the expression of collagen-I in the kidney of three different fibrotic mouse models. Pharmaceutical inhibition or genetic knockout of polycystin-2 reduced the expression of pro-fibrotic markers in fibrotic kidneys of two animal models. Moreover, conditional knockout of Pkd1 gene attenuated renal fibrosis in two mouse models. Enhancer of Zeste homolog 2 (EZH2) is a methyltransferase which plays important role in renal fibrosis. We found that EZH2 was positively correlated with the expression of polycystins in fibrotic kidneys. Conditional knockout of the Ezh2 attenuated the anti-fibrotic responses induced by Pkd1 deletion in fibrotic kidneys.

In conclusion, polycystins are up-regulated in injured kidneys and polycystins promote renal fibrosis through up-regulation of EZH2 in fibrotic kidneys, suggesting that renal fibrosis in polycystic kidneys is not the direct effect of PKD gene mutations.

SignificanceMutations of polycystins cause autosomal polycystic kidney disease (ADPKD) which is associated with increased renal fibrosis in cystic kidneys. However, the role of polycystin in renal fibrosis remains unclear. Here, we showed that the expression of polycystins is positively correlated with the development of renal fibrosis in injured kidneys. Pharmaceutical inhibition or genetic deletion of polycystin-2 attenuated renal fibrosis in fibrotic kidneys. Moreover, genetic deletion of polycystin-1 inhibited renal fibrosis in two different animal models. Finally, we showed that ablation of the methyltransferase EZH2 mitigated the anti-fibrotic effect of polycystin-1 deletion in fibrotic kidneys. We conclude that polycystin promotes through EZH2 and renal fibrosis in cystic kidney is not the direct effect of PKD gene mutations. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. Mouse lemur, a new primate model to investigate cardiac pacemaker activity. BackgroundGray mouse lemurs (Microcebus murinus) are one of the smallest and most ancestral primates known. Their size falls in between that of mice and rats, and their genetic proximity to humans makes them an emerging model to study age-related neurodegeneration. Since Microcebus murinus replicates similar senescence processes of humans, it constitutes a useful model also for studying cardiovascular dysfunctions. However, its cardiac physiology is unknown. Thus, we investigated sinoatrial node (SAN) pacemaker activity of mouse lemurs, presenting the first characterization of heart automaticity in nonhuman primates.

Methods and ResultsWe recorded the heart automaticity in freely moving mouse lemurs and in their SAN tissues and pacemaker myocytes. Mouse lemurs have a heart rate (HR) in between those of mice and rats and a similar generation of electrical impulses in the SAN tissue. Conversely, their ventricular depolarization is more similar to large mammals. Mouse lemur SAN myocytes express the main pacemaker currents at densities similar to mice: the hyperpolarization-activated current (If) and the L-type (ICa,L) and T-type (ICa,T) calcium currents. Moreover, mouse lemur SAN myocytes display spontaneous ryanodine receptor (RyR)-dependent Ca2+ transients. Using muscle-derived stem cells (MDSCs) from the Microcebus murinus skeletal muscle, we also differentiated pacemaker-like (PML) cells showing spontaneous automaticity and expressing markers of native SAN myocytes (HCN4 and connexin-45).

ConclusionsIn conclusion, our characterization of heart automaticity in Microcebus murinus provides new opportunities for comparative cardio-physiology studies in primates and with humans, and for pharmacologic testing of molecules that could assist in modulating age-related dysfunctions of automaticity. An epigenomic shift in amygdala marks the transition to maternal behaviors in alloparenting virgin female mice In many species, adults care for young offspring that are not their own, a phenomenon called alloparenting. However, most nonparental adults must be sensitized by repeated or extended exposures to newborns before they will robustly display parental-like behaviors. To capture neurogenomic events underlying the transition to active parental caring behaviors, we analyzed brain gene expression and chromatin profiles of virgin female mice co-housed with mothers during pregnancy and after birth. After an initial display of antagonistic behaviors and a surge of defense-related gene expression, we observed a dramatic shift in the chromatin landscape specifically in amygdala of the pup-exposed virgin females, accompanied by a dampening of anxiety-related gene expression. This epigenetic shift coincided with hypothalamic expression of the oxytocin gene and the emergence of behaviors and gene expression patterns classically associated with maternal care. The results outline a neurogenomic program associated with dramatic behavioral changes and suggest molecular networks relevant to human postpartum mental health. EMBED: a low dimensional reconstruction of subject-specific gut microbiome dynamics based on ecological normal modes Dimensionality reduction can offer unique insights into high dimensional microbiome dynamics by leveraging collective abundance fluctuations of multiple bacteria driven by similar ecological perturbations. However, methods providing lower-dimensional representations of microbiome dynamics both at the community and individual taxa level are not currently available. To that end, we present EMBED: Essential MicroBiomE Dynamics, a probabilistic non-linear tensor factorization approach. Similar to normal mode analysis in structural biophysics, EMBED infers ecological normal modes (ECNs), which represent the unique orthogonal modes capturing the collective behavior of microbial communities. A very small number of ECNs can accurately approximate microbiome dynamics across multiple data sets. Inferred ECNs reflect specific ecological behaviors, providing natural templates along which the dynamics of individual bacteria may be partitioned. Moreover, the multi-subject treatment in EMBED systematically identifies subject-specific and universal abundance dynamics that are not detected by traditional approaches. Collectively, these results highlight the utility of EMBED as a versatile dimensionality reduction tool for studies of microbiome dynamics. A multiscale lipid and cellular atlas of the human kidney The human kidney is composed of over 26 cell types that actively coordinate with each other to form higher-order structures, such as the nephron. It is not yet understood how these structures vary throughout a single organ or amongst the same organs within the human population. We have developed an extensive lipid and cellular atlas of the human kidney consisting of over 3 million cells comprising 75,000 functional tissue units (i.e., glomeruli, proximal tubules, distal tubules, and collecting ducts) from 13 human subjects. This atlas was developed using spatially registered and integrated technologies consisting of imaging mass spectrometry, multiplexed immunofluorescence, stained microscopy, and autofluorescence microscopy, to comprehensively probe large (i.e., centimeter-sized) areas of tissue. The cellular organization and lipid profiles of glomeruli, proximal tubules, distal tubules, and collecting ducts were discovered through these multimodal imaging data as well as their intra- and inter-subject variance. Relating the lipid profiles obtained from imaging mass spectrometry to distinct cell types obtained from immunofluorescence allowed us to hypothesize the functional role of specific phospholipids that have not previously been described. These hypotheses include subject characteristics, such as BMI and sex. The integrated data from the aforementioned datasets provide a valuable reference for kidney researchers, are publicly available through the NIH Human Biomolecular Atlas Program (https://portal.hubmapconsortium.org/), and discussed below. Migration Restores Hybrid Incompatibility Driven By Mitochondrial-Nuclear Sexual Conflict In the mitochondrial genome, sexual asymmetry in transmission allows the accumulation of male-harming mutations since selection acts only on the effect of the mutation in females. Called the "Mothers Curse", this phenomenon induces a selective pressure for nuclear variants that compensate for this reduction in male fitness. Previous work has demonstrated the existence of these interactions and their potential to act as Dobzhansky-Muller incompatibilities, contributing to reproductive isolation between populations. However, it is not clear how readily they would give rise to and sustain hybrid incompatibilities. Here, we use computer simulations in SLiM 3 to investigate the consequences of sexually antagonistic mitochondrial-nuclear interactions in a subdivided population. We consider distinct migration schemes and vary the chromosomal location, and consequently the transmission pattern, of nuclear restorers. Disrupting these co-evolved interactions results in less-fit males, skewing the sex ratio toward females. Restoration of male fitness depends on both the chromosomal location of nuclear restorer loci and the migration scheme. Our results show that these interactions may act as Dobzhansky-Muller incompatibilities, but their strength is not enough to drive population isolation. Overall, this model shows the varied ways in which populations can respond to migrations disruption of co-evolved mitochondrial-nuclear interactions. Multi-layered regulation of neuroectoderm differentiation by retinoic acid in a primitive streak-like context The formation of the primitive streak (PS) and the subsequent induction of neuroectoderm are hallmarks of gastrulation. Combining an in vitro reconstitution of this process based on mouse embryonic stem cells (mESCs) with a collection of knockouts in reporter mESC lines, we reassessed the contribution of retinoic acid (RA) signaling at early stages of neural commitment and its cross-talk with TGF{beta} and Wnt signaling inhibition. Single-cell RNA sequencing analysis captured the temporal unfolding of cell type diversification from epiblast and primitive streak-like cells up to the emergence of anterior and posterior neural fates. In conditions thought to lack RA synthesis, we discovered a hitherto unidentified residual RA production via a sensitive RA reporter. Genetic perturbations proved that the RA-degrading enzyme Cyp26a1 safeguard the developmental capabilities of the PS-like cells, limiting neural differentiation in wild type and in Chrd-/-Nog-/- or Dkk1-/- cells. Finally, the knockout of the three RAR receptors highlighted their function as negative regulators of loci critical for neural induction. Overall, we identified two mechanisms whereby components of the RA pathway can control the formation of neural progenitors in our PS-like context: a RA-dependent neural induction gated by RARs, and a receptor-mediated repression in the absence of ligand. An ALE meta-analytic review of top-down and bottom-up processing of music in the brain A remarkable feature of the human brain is its ability to integrate information from the environment with internally generated content. The integration of top-down and bottom-up processes during complex multi-modal human activities, however, is yet to be fully understood. Music provides an excellent model for understanding this since music listening leads to the urge to move, and music making entails both playing and listening at the same time (i.e. audio-motor coupling). Here, we conducted activation likelihood estimation (ALE) meta-analyses of 130 neuroimaging studies of music perception, production and imagery, with 2660 foci, 139 experiments, and 2516 participants. We found that music perception and production rely on auditory cortices and sensorimotor cortices, while music imagery recruits distinct parietal regions. This indicates that the brain requires different structures to process similar information which is made available either by an interaction with the environment (i.e. bottom-up) or by internally generated content (i.e. top-down). Comparative immunogenicity of bacterially expressed soluble trimers and nanoparticle displayed influenza hemagglutinin stem immunogens Current influenza vaccines need to be updated annually due to mutations in the globular head of the viral surface protein, hemagglutinin (HA). To address this, vaccine candidates have been designed based on the relatively conserved HA stem domain and have shown protective efficacy in animal models. Oligomerization of the antigens either by fusion to oligomerization motifs or display on self-assembling nanoparticle scaffolds, can induce more potent immune responses compared to the corresponding monomeric antigen due to multivalent engagement of B-cells. Since nanoparticle display can increase manufacturing complexity, and often involves one or more mammalian cell expressed components, it is important to characterize and compare various display and oligomerization scaffolds. Using a structure guided approach, we successfully displayed multiple copies of a previously designed soluble, trimeric influenza stem domain immunogen, pH1HA10, on the ferritin like protein, MsDps2 (12 copies), Ferritin (24 copies) and Encapsulin (180 copies). All proteins were expressed in Escherichia coli. The nanoparticle fusion immunogens were found to be well folded and bound to the influenza stem directed broadly neutralizing antibodies with high affinity. An 8.5 [A] Cryo-EM map of Msdps2-pH1HA10 confirmed the successful design of the nanoparticle fusion immunogen. Mice immunization studies with the soluble trimeric stem and nanoparticle fusion constructs revealed that all of them were immunogenic, and protected mice against homologous (A/Belgium/145-MA/2009) and heterologous (A/Puerto Rico/8/1934) challenge with 10MLD50 mouse adapted virus. Although nanoparticle display conferred a small but statistically significant improvement in protection relative to the soluble trimer in a homologous challenge, heterologous protection was similar in both nanoparticle-stem immunized and trimeric stem immunized groups. Such rapidly producible, bacterially expressed antigens and nanoparticle scaffolds are useful modalities to tackle future influenza pandemics. Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin BackgroundEvolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation at many enhancers, distal regulatory elements that control spatio-temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific, code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

ResultsWe present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that is able to predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Then, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity and liver enhancers associated with liver regeneration have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

ConclusionThe framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements. Synthetic Cooling Agents in US-marketed E-cigarette Refill Liquids and Popular Disposable E-cigarettes: Chemical Analysis and Risk Assessment. BackgroundMenthol, through its cooling and pleasant sensory effects, facilitates smoking and tobacco product initiation, resulting in the high popularity of mint/menthol-flavored E-cigarettes. More recently, E-cigarette vendors started marketing synthetic cooling agents as additives that impart a cooling effect but lack a characteristic minty odor. Knowledge about content of synthetic coolants in US-marketed E-cigarette products and associated health risks is limited.

MethodsE-liquid vendor sites were searched with the terms "koolada", "kool/cool", "ice", or WS-3/WS-23, denoting individual cooling agents, and relevant refill E-liquids were purchased. "Ice" flavor varieties of Puffbar, the most popular disposable E-cigarette brand, were compared with non-"Ice" varieties. E-liquids were characterized, and synthetic coolants quantified using GC/MS. Margin of exposure (MOE), a risk assessment parameter, was calculated to assess the risk associated with synthetic coolant exposure from E-cigarette use.

ResultsWS-3 was detected in 24/25 refill E-liquids analyzed. All Puffbar flavor varieties contained either WS-23 (13/14) or WS-3 (5/14), in both "Ice"- and non-"Ice" flavors. Modeling consumption of WS-3 from vaped E-liquids, resulted in MOEs below the safe margin of 100 for most daily use scenarios. MOEs for WS-23 were <100 for 10/13 Puffbar flavors in all use scenarios. Puffbar power specifications are identical to Juul devices.

ConclusionsSynthetic cooling agents (WS-3/WS-23) were present in US-marketed E-cigarettes, at levels that may result in consumer exposures exceeding safety thresholds set by regulatory agencies. Synthetic coolants are not only found in mint-or menthol-flavored products, but also in fruit- and candy-flavored products, including popular disposable E-cigarette products such as Puffbar.

ImplicationsSynthetic cooling agents are widely used in "kool/cool"- and "ice"-flavored E-liquids and in E-liquids without these labels, both as a potential replacement for menthol or to add cooling  notes to non-menthol flavors. These agents may be used to bypass current and future regulatory limits on menthol content in tobacco products, and not just E-cigarettes. Since synthetic cooling agents are odorless, they may not fall under the category of "characterizing flavor", potentially circumventing regulatory measures based on this concept. Regulators need to consider the additional health risks associated with exposure to synthetic cooling agents. Skeletal muscle methylome and transcriptome integration reveals profound sex differences related to muscle function and substrate metabolism Nearly all human complex traits and diseases exhibit some degree of sex differences, with epigenetics being one of the main contributing factors. Various tissues display sex differences in DNA methylation, however this has not yet been explored in skeletal muscle, despite skeletal muscle being among the tissues with the most transcriptomic sex differences. For the first time, we investigated the effect of sex on autosomal DNA methylation in human skeletal muscle across three independent cohorts (Gene SMART, FUSION, and GSE38291) using a meta-analysis approach, totalling 369 human muscle samples (222 males, 147 females), and integrated this with known sex-biased transcriptomics. We found 10,240 differentially methylated regions (DMRs) at FDR < 0.005, 94% of which were hypomethylated in males, and gene set enrichment analysis revealed that differentially methylated genes were involved in muscle contraction and substrate metabolism. We then investigated biological factors underlying DNA methylation sex differences and found that circulating hormones were not associated with differential methylation at sex-biased DNA methylation loci, however these sex-specific loci were enriched for binding sites of hormone-related transcription factors (with top TFs including androgen (AR), estrogen (ESR1), and glucocorticoid (NR3C1) receptors). Fibre type proportions were associated with differential methylation across the genome, as well as across 16 % of sex-biased DNA methylation loci (FDR < 0.005). Integration of DNA methylomic results with transcriptomic data from the GTEx database and the FUSION cohort revealed 326 autosomal genes that display sex differences at both the epigenome and transcriptome levels. Importantly, transcriptional sex-biased genes were overrepresented among epigenetic sex-biased genes (p-value = 4.6e-13), suggesting differential DNA methylation and gene expression between male and female muscle are functionally linked. Finally, we validated expression of three genes with large effect sizes (FOXO3A, ALDH1A1, and GGT7) in the Gene SMART cohort with qPCR. GGT7, involved in antioxidant metabolism, displays male-biased expression as well as lower methylation in males across the three cohorts. In conclusion, we uncovered 8,420 genes that exhibit DNA methylation differences between males and females in human skeletal muscle that may modulate mechanisms controlling muscle metabolism and health.

SignificanceThe importance of uncovering biological sex differences and their translation to physiology has become increasingly evident. Using a large-scale meta-analysis of three cohorts, we perform the first comparison of genome-wide skeletal muscle DNA methylation between males and females, and identify thousands of genes that display sex-differential methylation. We then explore intrinsic biological factors that may be underlying the DNA methylation sex differences, such as fibre type proportions and sex hormones. Leveraging the GTEx database, we identify hundreds of genes with both sex-differential expression and DNA methylation in skeletal muscle. We further confirm the sex-biased genes with gene expression data from two cohorts included in the methylation meta-analysis. Our study integrates genomewide sex-biased DNA methylation and expression in skeletal muscle, shedding light on distinct sex differences in skeletal muscle.