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.
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.
Cryptic pathogen-sugar interactions revealed by universal saturation transfer analysis Many host pathogen interactions such as human viruses (including non-SARS-coronaviruses) rely on attachment to host cell-surface glycans. There are conflicting reports about whether the Spike protein of SARS-CoV-2 binds to sialic acid commonly found on host cell-surface N-linked glycans. In the absence of a biochemical assay, the ability to analyze the binding of glycans to heavily- modified proteins and resolve this issue is limited. Classical Saturation Transfer Difference (STD) NMR can be confounded by overlapping sugar resonances that compound with known experimental constraints. Here we present  universal saturation transfer analysis (uSTA), an NMR method that builds on existing approaches to provide a general and automated workflow for studying protein-ligand interactions. uSTA reveals that B-origin-lineage-SARS-CoV-2 spike trimer binds sialoside sugars in an  end on manner and modelling guided by uSTA localises binding to the spike N-terminal domain (NTD). The sialylated-polylactosamine motif is found on tetraantennary human N-linked-glycoproteins in deeper lung and may have played a role in zoonosis. Provocatively, sialic acid binding is abolished by mutations in some subsequent SARS- CoV-2 variants-of-concern. A very high resolution cryo-EM structure confirms the NTD location and  end on mode; it rationalises the effect of NTD mutations and the structure-activity relationship of sialic acid analogues. uSTA is demonstrated to be a robust, rapid and quantitative tool for analysis of binding, even in the most demanding systems.

Extended AbstractThe surface proteins found on both pathogens and host cells mediate entry (and exit) and influence disease progression and transmission. Both types can bear host-generated post- translational modifications such as glycosylation that are essential for function but can confound biophysical methods used for dissecting key interactions. Several human viruses (including non- SARS-coronaviruses) attach to host cell-surface N-linked glycans that include forms of sialic acid (sialosides). There remains, however, conflicting evidence as to if or how SARS-associated coronaviruses might use such a mechanism. Here, we demonstrate quantitative extension of  saturation transfer protein NMR methods to a complete mathematical model of the magnetization transfer caused by interactions between protein and ligand. The method couples objective resonance-identification via a deconvolution algorithm with Bloch-McConnell analysis to enable a structural, kinetic and thermodynamic analysis of ligand binding beyond previously-perceived limits of exchange rates, concentration or system. Using an automated and openly available workflow this  universal saturation transfer analysis (uSTA) can be readily-applied in a range of even heavily-modified systems in a general manner to now obtain quantitative binding interaction parameters (KD, kEx). uSTA proved critical in mapping direct interactions between natural sialoside sugar ligands and relevant virus-surface attachment glycoproteins - SARS-CoV-2-spike and influenza-H1N1-haemagglutinin variants - by quantitating ligand signal in spectral regions otherwise occluded by resonances from mobile protein glycans (that also include sialosides). In B- origin-lineage-SARS-CoV-2 spike trimer  end on-binding to sialoside sugars was revealed contrasting with  extended surface-binding for heparin sugar ligands; uSTA-derived constraints used in structural modelling suggested sialoside-glycan binding sites in a beta-sheet-rich region of spike N-terminal domain (NTD). Consistent with this, uSTA-glycan binding was minimally- perturbed by antibodies that neutralize the ACE2-binding domain (RBD) but strongly disrupted in spike from the B1.1.7/alpha and B1.351/beta variants-of-concern, which possess hotspot mutations in the NTD. Sialoside binding in B-origin-lineage-NTD was unequivocally pinpointed by cryo-EM to a site that is created from residues that are notably deleted in variants (e.g. H69,V70,Y145 in alpha). An analysis of beneficial genetic variances in cohorts of patients from early 2020 suggests a model in which this site in the NTD of B-origin-lineage-SARS-CoV-2 (but not in alpha/beta-variants) may have exploited a specific sialylated-polylactosamine motif found on tetraantennary human N-linked-glycoproteins in deeper lung. Together these confirm a novel binding mode mediated by the unusual NTD of SARS-CoV-2 and suggest how it may drive virulence and/or zoonosis via modulation of glycan attachment. Since cell-surface glycans are widely relevant to biology and pathology, uSTA can now provide ready, quantitative, widespread analysis of complex, host-derived and post-translationally modified proteins with putative ligands relevant to disease even in previously confounding complex systems. Enhancing CDK4/6 inhibitor therapy for medulloblastoma using nanoparticle delivery and scRNA-seq-guided combination with sapanisertib CDK4/6 inhibitors hold promise for brain tumor treatment, but efficacy has been limited by recurrence in both preclinical models and clinical trials. To address recurrence, we tested a nanoparticle formulation of the CDK4/6 inhibitor palbociclib (POx-palbo) in mice genetically-engineered to develop SHH-driven medulloblastoma. We then analyzed medulloblastomas in mice receiving palbociclib treatment, and compared the efficacy of combining palbociclib with specific inhibitors suggested by our analysis. POx-Palbo showed reduced toxicity compared to conventional palbociclib, was tolerable in parenteral administration, improved CNS pharmacokinetics, and extended survival of mice with medulloblastoma. Recurrence, however, remained problematic as fractions of tumor cells proliferated during therapy. ScRNA-seq identified a gene expression pattern unique to proliferating medulloblastoma cells in POx-Palbo-treated mice, marked by up-regulation of the glutamate transporter Slc1a2 and down-regulation of diverse ribosomal genes. Reduced mTORC1 signaling, suggested by ribosomal suppression in POx-Palbo-treated tumors was confirmed by decreased 4EBP1 phosphorylation (p4EBP1). Further reducing mTORC1 activity by combining POx-Palbo with the mTORC1 inhibitor sapanisertib produced mutually enhancing effects, with increased suppression of both pRB and p4EBP1, and prolonged mouse survival compared to either agent alone. In contrast, targeting cell cycle progression by combining POx-Palbo with the SHH-pathway inhibitor vismodegib, or with the replication-targeting agents gemcitabine or etoposide, failed to enhance efficacy. Our data show the potential of nanoparticle formulation and scRNA-seq analysis of resistance to improve brain tumor treatment, and identify POx-palbo plus sapanisertib as effective combinatorial therapy for SHH medulloblastoma. This combination may be appropriate for testing in patients with recurrence, who need new options. Cell cycle-independent integration of stress signals promotes Non-G1/G0 quiescence entry Cellular quiescence is a non-proliferative state required for cell survival under stress and during development. In most quiescent cells, proliferation is stopped in a reversible state of low Cdk1 kinase activity; in many organisms, however, quiescent states with high Cdk1 activity can also be established through still uncharacterized stress or developmental mechanisms. Here, we used a microfluidics approach coupled to phenotypic classification by machine learning to identify stress pathways associated with starvation-triggered high-Cdk1 quiescent states in Saccharomyces cerevisiae. We found that low- and high-Cdk1 quiescent states shared a core of stress-associated processes, such as autophagy, protein aggregation, and mitochondrial upregulation, but differed in the nuclear accumulation of the stress transcription factors Xbp1, Gln3, and Sfp1. The decision between low- or high-Cdk1 quiescence was controlled by cell cycle-independent accumulation of Xbp1, which acted as a time-delayed integrator of the duration of stress stimuli. Our results show how cell cycle-independent stress-activated factors promote cellular quiescence outside of G1/G0. Citation needed?Wikipedia bibliometrics during the first waveof the COVID-19 pandemic. BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. Thermal modulation of Zebrafish exploratory statistics reveals constraints on individual behavioral variability BackgroundVariability is a hallmark of animal behavior. It contributes to survival by endowing individuals and populations with the capacity to adapt to ever-changing environmental conditions. Intra-individual variability is thought to reflect both endogenous and exogenous modulations of the neural dynamics of the central nervous system. However, how variability is internally regulated and modulated by external cues remains elusive. Here we address this question by analyzing the statistics of spontaneous exploration of freely swimming zebrafish larvae, and by probing how these locomotor patterns are impacted when changing the water temperatures within an ethologically relevant range.

ResultsWe show that, for this simple animal model, five short-term kinematic parameters - interbout interval, turn amplitude, travelled distance, turn probability and orientational flipping rate - together control the long-term exploratory dynamics. We establish that the bath temperature consistently impacts the means of these parameters, but leave their pairwise covariance unchanged. These results indicate that the temperature merely controls the sampling statistics within a well-defined kinematic space delineated by this robust statistical structure. At a given temperature, individual animals explore the behavioral space over a timescale of tens of minutes, suggestive of a slow internal state modulation that could be externally biased through the bath temperature. By combining these various observations into a minimal stochastic model of navigation, we show that this thermal modulation of locomotor kinematics results in a thermophobic behavior, complementing direct gradient-sensing mechanisms.

ConclusionsThis study establishes the existence of a well-defined locomotor space accessible to zebrafish larvae during spontaneous exploration, and quantifies self-generated modulation of locomotor patterns. Intra-individual variability reflects a slow diffusive-like probing of this space by the animal. The bath temperature in turn restricts the sampling statistics to sub-regions, endowing the animal with basic thermophobicity. This study suggests that in Zebrafish, as well as in other ectothermic animals, ambient temperature could be used to efficiently manipulate internal states in a simple and ethological way. NPAS4 in the medial prefrontal cortex mediates chronic social defeat stress- induced anhedonia and dendritic spine loss Chronic stress can produce reward system deficits (i.e. anhedonia) and other common symptoms associated with depressive disorders, as well as neural circuit hypofunction in the medial prefrontal cortex (mPFC). However, the molecular mechanisms by which chronic stress promotes depressive-like behavior and hypofrontality remain unclear. We show here that the neuronal activity-regulated transcription factor, NPAS4, in the mPFC is regulated by chronic social defeat stress (CSDS), and its required in this brain region for CSDS-induced changes in sucrose preference and natural reward motivation. Interestingly, NPAS4 is not required for CSDS-induced social avoidance or anxiety-like behavior. We also find that mPFC NPAS4 is required for CSDS-induced reduction of pyramidal cell dendritic spine density, revealing a relationship between mPFC dendritic spine changes and anhedonia-like behavior, but not social avoidance behavior. Finally, transcriptomic analysis from the mPFC revealed that NPAS4 influences expression of numerous genes linked to glutamatergic synapses and ribosomal function, as well as many dysregulated genes observed in common neuropsychiatric disorders, including depression. Together our findings reveal an essential role for the activity-regulated transcription factor, NPAS4, in chronic stress-induced mPFC hypofunction and anhedonia. 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. Sperm Histone H3 Lysine 4 tri-methylation serves as a metabolic sensor of paternal obesity and is associated with the inheritance of metabolic dysfunction ObjectiveParental environmental exposures can strongly influence descendant risks for adult disease. How paternal obesity changes the sperm chromatin leading to the acquisition of metabolic disease in offspring remains controversial and ill-defined. The objective of this study was to assess: (1) whether obesity induced by a high-fat diet alters sperm histone methylation; (2) whether paternal obesity can induce metabolic disturbances across generations; (3) whether there could be cumulative damage to the sperm epigenome leading to enhanced metabolic dysfunction in descendants; and (4) whether obesity-sensitive regions associate with embryonic epigenetic and transcriptomic profiles. Using a genetic mouse model of epigenetic inheritance, we investigated the role of histone H3 lysine 4 methylation (H3K4me3) in the paternal transmission of metabolic dysfunction. This transgenic mouse overexpresses the histone demethylase enzyme KDM1A in the developing germline and has an altered sperm epigenome at the level of histone H3K4 methylation. We hypothesized that challenging transgenic sires with a high-fat diet would further erode the sperm epigenome and lead to enhanced metabolic disturbances in the next generations.

MethodsTo assess whether paternal obesity can have inter- or transgenerational impacts, and if so, to identify potential mechanisms of this non-genetic inheritance, we used wildtype C57BL/6NCrl and transgenic males with a pre-existing altered sperm epigenome. To induce obesity, sires were fed either a control or high-fat diet (10% or 60% kcal fat, respectively) for 10-12 weeks, then bred to wildtype C57BL/6NCrl female fed a regular diet. F1 and F2 descendants were characterized for metabolic phenotypes by examining the effects of paternal obesity by sex, on body weight, fat mass distribution, the liver transcriptome, intraperitoneal glucose and insulin tolerance tests. To determine whether obesity altered the F0 sperm chromatin, native chromatin immunoprecipitation-sequencing targeting H3K4me3 was performed. To gain insight into mechanisms of paternal transmission, we compared our sperm H3K4me3 profiles with embryonic and placental chromatin states, histone modification and gene expression profiles.

ResultsObesity-induced alterations in H3K4me3 occurred at genes implicated in metabolic, inflammatory, and developmental processes. These processes were associated with offspring metabolic dysfunction and corresponded to genes enriched for H3K4me3 in embryos, and overlapped embryonic and placenta gene expression profiles. Transgenerational susceptibility to metabolic disease was only observed when obese F0 had a pre-existing modified sperm epigenome. This coincided with increased H3K4me3 alterations in sperm and more severe phenotypes affecting their offspring.

ConclusionsOur data suggest sperm H3K4me3 might serve as a metabolic sensor that connects paternal diet with offspring phenotypes via the placenta. This non-DNA based knowledge of inheritance has the potential to improve our understanding of how environment shapes heritability and may lead to novel routes for the prevention of disease. This study highlights the need to further study the connection between the sperm epigenome, placental development and childrens health. High temporal resolution systems profiling reveals distinct patterns of interferon response after Covid-19 mRNA vaccination and SARS-CoV2 infection Knowledge of the mechanisms underpinning the development of protective immunity conferred by mRNA vaccines is fragmentary. Here we investigated responses to COVID-19 mRNA vaccination via ultra-low-volume sampling and high-temporal-resolution transcriptome profiling (23 subjects across 22 timepoints, and with 117 COVID-19 patients used as comparators). There were marked differences in the timing and amplitude of the responses to the priming and booster doses. Notably, we identified two distinct interferon signatures. The first signature (A28/S1) was robustly induced both post-prime and post-boost and in both cases correlated with the subsequent development of antibody responses. In contrast, the second interferon signature (A28/S2) was robustly induced only post-boost, where it coincided with a transient inflammation peak. In COVID19 patients, a distinct phenotype dominated by A28/S2 was associated with longer duration of intensive care. In summary, high-temporal-resolution transcriptomic permitted the identification of post- vaccination phenotypes that are determinants of the course of COVID-19 disease. Effects of substrate on shrimp growth, water quality and bacterial community in the biofloc system nursing Penaeus vannamei postlarvae This study aimed to investigate the effects of substrate on water quality, shrimp growth and bacterial community in the biofloc system with a salinity of 5{per thousand}. Two treatments, biofloc system with (sB) or without (nB) addition of elastic solid packing filler (nylon) as substrate, were set up. Penaeus vannamei postlarvae (PL, [~] stage 15) were stocked at a density of 4000 PL m-3 in each treatment with triplicates for a 28-days culture experiment, taking glucose as carbon source (C:N 15:1). Results showed that the survival rate (96.3{+/-}3.6%), FCR (0.76{+/-}0.06) and productivity (1.54{+/-}0.12 kg m-3) in sB treatment were significantly better than those in nB treatment (81.0{+/-}7.1%, 0.98{+/-}0.08 and 1.14{+/-}0.09 kg m-3, P<0.05). All water parameters were in the recommended ranges. Substrate showed significant effect on TAN, TSS, turbidity, biofloc volume, pH and carbonate alkalinity (P < 0.05). Actinobacteria (4.0-22.7%), Bacteroidetes (10.4-33.5%), Firmicutes (0.2-11.2%), Planctomycetes (4.0-14.9%) and Proteobacteria (29.4-59.0%) were the most dominant phyla for both treatments. However, the bacterial community in sB treatment showed to be significantly different from that in nB treatment (Jaccard distance 0.94{+/-}0.01, P=0.001). Substrate showed significant effects on Shannon, Heip, Pielou and Simpson index, as well as relative abundances of Actinobacteria, Bacteroidetes and Planctomycetes (P < 0.05). The results suggested that addition of substrate affected the shrimp growth, water quality and bacterial community in the biofloc system nursing P. vannamei PL with a 5{per thousand} salinity. METHYLOME DYNAMICS UPON PROTEASOME INHIBITION BY THE PSEUDOMONAS VIRULENCE FACTOR SYRINGOLIN A DNA methylation is an important epigenetic mark required for proper gene expression and silencing of transposable elements. DNA methylation patterns can be modified by environmental factors such as pathogen infection, where modification of DNA methylation can be associated with plant resistance. To counter the plant defense pathways, pathogens produce effectors molecule, several of which act as proteasome inhibitors. Here we investigated the effect of proteasome inhibition by the bacterial virulence factor Syringolin A on genome-wide DNA methylation. We show that Syringolin A treatment results in an increase of DNA methylation at centromeric and pericentromeric regions of Arabidopsis chromosomes. We identify several CHH DMRs that are enriched in the proximity of transcriptional start sites. Syringolin A treatment does not result in significant changes in small RNA composition. However, significant changes in genome transcriptional activity can be observed, including a strong upregulation of resistance genes that are located on chromosomal arms. We hypothesize that DNA methylation changes could be linked to the upregulation of some atypical members of the de novo DNA methylation pathway: AGO3, AGO9 and DRM1. Our data suggests that modification of genome-wide DNA methylation resulting from an inhibition of the proteasome by bacterial effectors could be part of a epi-genomic arms race against pathogens. Hierarchical Bayesian Integrated Model for Estimating Migratory Bird Harvest in Canada The Canadian Wildlife Service (CWS) requires reliable estimates of the harvest of migratory game birds, including waterfowl, to effectively manage populations of these hunted species. The National Harvest Survey is an annual survey of hunters who purchase Canadas mandatory migratory game bird hunting permit, integrating information from a survey of hunting activity with information from a separate survey of species composition in the harvest. We use these survey data to estimate the number of birds harvested for each species, as well as hunting activity metrics such as the number of active hunters and days spent hunting. The analytical methods used to generate these estimates have not changed since the survey was first designed in the early 1970s. Here we describe a new hierarchical Bayesian integrated model, which replaces the series of ratio estimators that comprised the old model. We are now using this new model to generate estimates for migratory bird harvests as of the 2019-2020 hunting season, and to generate updated estimates for all earlier years. The hierarchical Bayesian model uses over-dispersed Poisson distributions to model mean hunter activity and harvest (zero inflated Poisson and zero truncated Poisson, respectively). It also includes multinomial distributions to model some key components including, variation in total harvest across periods of the hunting season, the species composition of the harvest within each of those periods, and the age and sex composition in the harvests of a given species. We estimated the parameters of the Poisson and the multinomial distributions for each year as random effects using first-difference time-series. This time-series component allows the model to share information across years and reduces the sensitivity of the estimates to annual sampling noise. The new model estimates are generally very similar to those from the old model, particularly for the species that occur most commonly in the harvest, and so the results do not suggest any major changes to harvest management decisions and regulations. However, estimates for all species from the new model are more precise and less susceptible to annual sampling error, particularly for species that occur less commonly in the harvest (e.g., sea ducks and other species of conservation concern). This new model, with its hierarchical Bayesian framework, will also facilitate future improvements and elaborations, allowing the incorporation of prior information from the rich literature and knowledge in game bird management and biology. Clonal evolution in serially passaged Cryptococcus neoformans x deneoformans hybrids reveals a heterogenous landscape of genomic change Cryptococcus neoformans x deneoformans hybrids (also known as serotype AD hybrids) are basidiomycete yeasts that are common in a clinical setting. Like many hybrids, the AD hybrids are largely locked at the F1 stage and are mostly unable to undergo normal meiotic reproduction. However, these F1 hybrids, which display a high ([~]10%) sequence divergence are known to genetically diversify through mitotic recombination and aneuploidy, and this diversification may be adaptive. In this study, we evolved a single AD hybrid genotype in six diverse environments by serial passaging and then used genome resequencing of evolved clones to determine evolutionary mechanisms of adaptation. The evolved clones generally increased fitness after passaging, accompanied by an average of 3.3 point mutations, 2.9 loss of heterozygosity (LOH) events, and 0.7 trisomic chromosomes per clone. LOH occurred through nondisjunction of chromosomes, crossing over consistent with break-induced replication, and gene conversion, in that order of prevalence. The breakpoints of these recombination events were significantly associated with regions of the genome with lower sequence divergence between the parents and clustered in subtelomeric regions, notably in regions that had undergone introgression between the two parental species. Parallel evolution was observed, particularly through repeated homozygosity via nondisjunction, yet there was little evidence of environment-specific parallel change for either LOH, aneuploidy, or mutations. These data show that AD hybrids have both a remarkable genomic plasticity and yet are challenged in the ability to recombine through sequence divergence and chromosomal rearrangements, a scenario likely limiting the precision of adaptive evolution to novel environments. Etmopterus lantern sharks use coelenterazine as the substrate for their luciferin-luciferase bioluminescence system The lantern shark genus Etmopterus is a group of deep-sea bioluminescent fishes. They emit blue light mainly from the ventral body surface, and the primary biological function is considered to be for camouflage by counterillumination. In this study, we detected both coelenterazine and coelenterazine-specific luciferase activity in the ventral photophore tissues. The results suggested that bioluminescence in lantern sharks is produced using coelenterazine as the substrate for the luciferin-luciferase reaction. CyanoHABs: unavoidable evolutionary ecological consequence for low nutrient-requiring cyanobacteria in eutrophic water The mechanism of cyanobacterial harmful algal blooms (CyanoHABs) is complicated and confusing. One major reason is they are studied primarily from an ecological perspective and on bloom-forming species only. This narrow angle loses a broader evolutionary and ecological context in which CyanoHABs occur and fails to provide information on relevant components to achieve a wholistic understanding.

To derive a comprehensive mechanism of CyanoHABs, we examine CyanoHABs through the overlooked evolutionary and ecological lenses: evolutionary radiation, ecological comparison with co-living algae, and recently identified genomic functional repertoire between blooming and non-blooming species. We found four key factors: (1) elaborate diverse functional repertoire and low nutrient requirement in cyanobacteria molded by early adaptive evolution, (2) cyanobacteria having lower nutrient requirements than co-living algae, (3) there is no directed evolution in biological functions toward water eutrophication in cyanobacteria, (4) the CyanoHAB-associated functional repertoire are more abundant and complete in blooming than non-blooming species. These factors lead us to postulate a preliminary mechanism of CyanoHABs as a synergistic quad: superior functional repertoire, established with long adaptive radiation under nutrient-deficient conditions and not evolved toward eutrophic conditions, enables cyanobacteria to efficiently utilize elevated nutrients under current eutrophic levels for excess growth and CyanoHABs thereof, due to their lower nutrient requirements than co-living algae. This working comprehensive mechanism without doubt needs further empirical testing, which can be undertaken with more comparative studies of multiple species using integrated systems biology approaches. 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. scJoint: transfer learning for data integration of atlas-scale single-cell RNA-seq and ATAC-seq Single-cell multi-omics data continues to grow at an unprecedented pace, and effectively integrating different modalities holds the promise for better characterization of cell identities. Although a number of methods have demonstrated promising results in integrating multiple modalities from the same tissue, the complexity and scale of data compositions typically present in cell atlases still pose a significant challenge for existing methods. Here we present scJoint, a transfer learning method to integrate atlas-scale, heterogeneous collections of scRNA-seq and scATAC-seq data. scJoint leverages information from annotated scRNA-seq data in a semi-supervised framework and uses a neural network to simultaneously train labeled and unlabeled data, enabling label transfer and joint visualization in an integrative framework. Using multiple atlas data and a biologically varying multi-modal data, we demonstrate scJoint is computationally efficient and consistently achieves significantly higher cell type label accuracy than existing methods while providing meaningful joint visualizations. This suggests scJoint is effective in overcoming the heterogeneity in different modalities towards a more comprehensive understanding of cellular phenotypes. Glucose dynamics during ozone exposure measured using radiotelemetry: Stress drivers BackgroundStress-related neurobehavioral and metabolic disorders are associated with altered circulating adrenal-derived hormones and hyperglycemia. Temporal assessment of glucose and these hormones is critical for insights on an individuals health.

ObjectivesHere we used implantable-telemetry in rats to assess real-time changes in circulating glucose during and after exposure to the air pollutant ozone, and link responses to circulating neuroendocrine stress and metabolic hormones. We also compared rodent glucose and corticosterone (cortisol in humans) responses to humans exposed to ozone.

MethodsFirst, using a cross-over design, we monitored glucose levels during single or repeated ozone exposures (0.0, 0.2, 0.4 and 0.8-ppm) and non-exposure periods in male Wistar-Kyoto-rats implanted with glucose-telemeters. A second cohort of un-implanted rats was exposed to ozone (0.0, 0.4 or 0.8-ppm) for 30-min, 1-hour, 2-hour, or 4-hour with hormones measured immediately after exposure. Then we assessed glucose metabolism in sham and adrenalectomized rats with or without pharmacological interventions of adrenergic and glucocorticoid receptors. Finally, we assessed glucose and cortisol in serum samples from a clinical study involving exposure of human volunteers to air or 0.3 ppm ozone.

ResultsOzone (0.8-ppm) caused hyperglycemia and hypothermia beginning 90-min into exposure, with reversal of effects 4-6 hours post-exposure. Glucose monitoring during four daily 4-hour ozone exposures revealed duration of hyperglycemia, adaptation, and diurnal variations. Ozone-induced hyperglycemia was preceded by increased adrenocorticotropic hormone, corticosterone, and epinephrine, but depletion of thyroid-stimulating, prolactin, and luteinizing hormones. Hyperglycemia was inhibited in rats that were adrenalectomized and/or treated with glucocorticoid inhibitor. The depletion of cortisol was dampened in humans exposed to ozone during intermittent exercise.

DiscussionWe demonstrate for the first time the temporality of neuroendocrine-stress-mediated biological sequalae responsible for ozone-induced metabolic dysfunction as exposure occurs. Real-time glucose monitoring with stress hormones assessment may be useful in identifying interactions among pollutants and stress-related illnesses. Litter commensal bacteria can limit the horizontal gene transfer of antimicrobial resistance to Salmonella in chickens Fostering a balanced gut microbiome through the administration of beneficial microbes that can competitively exclude pathogens has gained a lot of attention and use in human and animal medicine. However, little is known about how microbes affect the horizontal gene transfer of antimicrobial resistance (AMR). To shed more light on this question, we challenged neonatal broiler chicks raised on reused broiler chicken litter - a complex environment made up of decomposing pine shavings, feces, uric acid, feathers, and feed, with Salmonella Heidelberg (S. Heidelberg), a model pathogen. Neonatal chicks challenged with S. Heidelberg and raised on reused litter were more resistant to S. Heidelberg cecal colonization than chicks grown on fresh litter. Furthermore, chicks grown on reused litter were at a lower risk of colonization with S. Heidelberg strains that encoded AMR on IncI1 plasmids. We used 16S rRNA gene sequencing and shotgun metagenomics to show that the major difference between chicks grown on fresh litter and reused litter was the microbiome harbored in the litter and ceca. The microbiome of reused litter samples was more uniform and enriched in functional pathways related to the biosynthesis of organic and antimicrobial molecules than fresh litter samples. We found that E. coli was the main reservoir of plasmids encoding AMR and that the IncI1 plasmid was maintained at a significantly lower copy per cell in reused litter compared to fresh litter. These findings support the notion that commensal bacteria play an integral role in the horizontal transfer of plasmids encoding AMR to pathogens like Salmonella.

Importance/SignificanceAntimicrobial resistance spread is a worldwide health challenge, stemming in large part, from the ability of microorganisms to share their genetic material through horizontal gene transfer. To address this issue, many countries and international organization have adopted a One health approach to curtail the proliferation of antimicrobial resistant bacteria. This includes the removal and reduction of antibiotics used in food animal production and the development of alternatives to antibiotics. However, there is still a significant knowledge gap in our understanding of how resistance spreads in the absence of antibiotic selection and the role commensal bacteria play in reducing antibiotic resistance transfer. In this study, we show that commensal bacteria play a key role in reducing the horizontal gene transfer of antibiotic resistance to Salmonella and provide the identity of the bacterial species that potentially perform this function in broiler chickens and also postulate the mechanism involved. Glycolysis/gluconeogenesis specialization in microbes is driven by biochemical constraints of flux sensing Central carbon metabolism is highly conserved across microbial species, but can catalyze very different pathways depending on the organism and their ecological niche. Here, we study the dynamic re-organization of central metabolism after switches between the two major opposing pathway configurations of central carbon metabolism, glycolysis and gluconeogenesis in Escherichia coli, Pseudomonas aeruginosa and Pseudomonas putida. We combined growth dynamics and dynamic changes of intracellular metabolite levels with a coarse-grained model that integrates fluxes, regulation, protein synthesis and growth and uncovered fundamental limitations of the regulatory network: after nutrient shifts, metabolite concentrations collapse to their equilibrium, rendering the cell unable to sense which direction the flux is supposed to flow through the metabolic network. The cell can partially alleviate this by picking a preferred direction of regulation at the expense of increasing lag times in the opposite direction. Moreover, decreasing both lag times simultaneously comes at the cost of reduced growth rate or higher futile cycling between metabolic enzymes. These three trade-offs can explain why microorganisms specialize for either glycolytic or gluconeogenic substrates and can help elucidate the complex growth patterns exhibited by different microbial species.

Graphical synopsis O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY

Standfirst textMicrobes face a series of fundamental trade-offs that limit their ability to optimize simultaneously for both glycolytic and gluconeogenic growth.

Bullet pointsO_LILag times between glycolysis and gluconeogenesis show asymmetry in many microbes: A long lag in one direction, but a short lag in the other.
C_LIO_LILong lag times are caused by an inability to sense fluxes after nutrient shifts.
C_LIO_LIWith existing regulation, lag time asymmetry can only be overcome by reducing either growth rate or increasing futile cycling in metabolism.
C_LI Kidney failure alters parathyroid Pin1 phosphorylation and parathyroid hormone mRNA binding proteins leading to secondary hyperparathyroidism Parathyroid hormone (PTH) regulates calcium metabolism and bone strength. Chronic kidney disease leads to secondary hyperparathyroidism (SHP) which increases morbidity and mortality. In experimental SHP, the increased PTH gene expression is due to enhanced PTH mRNA stability, mediated by changes in its interaction with stabilizing AUF1 and destabilizing KSRP. Pin1 isomerizes target proteins and leads to KSRP dephosphorylation. In SHP, Pin1 isomerase activity is decreased and phosphorylated KSRP fails to bind PTH mRNA, resulting in high PTH mRNA stability and levels. However, the up- and down-stream mechanisms by which kidney failure stimulates the parathyroid remain elusive. We now reveal a pathway where kidney failure induces parathyroid Pin1 phosphorylation, linking Pin1, KSRP and PTH mRNA stability as key players for the high PTH expression in SHP. We identified by mass-spectrometry, changes in rat parathyroid proteome and phosphoproteome profiles induced by impaired renal function, including KSRP phosphorylation at Pin1 target sites. Furthermore, both acute and chronic kidney failure led to parathyroid-specific Pin1 Ser16 and Ser71 phosphorylation, which disrupts Pin1 activity. Accordingly, pharmacologic Pin1 inhibition, that mimics the decreased Pin1 activity in SHP, increased PTH expression ex-vivo in parathyroid organ cultures and in transfected cells, through the PTH mRNA protein-interacting element and KSRP phosphorylation at potential Pin1-binding motifs. Therefore, kidney failure leads to loss of parathyroid Pin1 activity by inducing Pin1 phosphorylation. This predisposes parathyroids to increase PTH production through impaired PTH mRNA decay that is dependent on KSRP phosphorylation. Pin1 and KSRP phosphorylation and the Pin1-KSRP-PTH mRNA axis thus drive SHP. A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. Placental uptake and metabolism of 25(OH)Vitamin D determines its activity within the fetoplacental unit Pregnancy 25-hydroxyvitamin D (25(OH)D) concentrations are associated with maternal and fetal health outcomes. Using physiological human placental perfusion and villous explants, we investigate the role of the placenta in regulating the relationships between maternal 25(OH)D and fetal physiology. We demonstrate active placental uptake of 25(OH)D3 by endocytosis, placental metabolism of 25(OH)D3 into 24,25-dihydroxyvitamin D3 and active 1,25-dihydroxyvitamin D [1,25(OH)2D3], with subsequent release of these metabolites into both the maternal and fetal circulations. Active placental transport of 25(OH)D3 and synthesis of 1,25(OH)2D3 demonstrate that fetal supply is dependent on placental function rather than simply the availability of maternal 25(OH)D3. We demonstrate that 25(OH)D3 exposure induces rapid effects on the placental transcriptome and proteome. These map to multiple pathways central to placental function and thereby fetal development, independent of vitamin D transfer. Our data suggest that the underlying epigenetic landscape helps dictate the transcriptional response to vitamin D treatment. This is the first quantitative study demonstrating vitamin D transfer and metabolism by the human placenta, with widespread effects on the placenta itself. These data demonstrate a complex interplay between vitamin D and the placenta and will inform future interventions using vitamin D to support fetal development and maternal adaptations to pregnancy. New Segmentation and Feature Extraction Algorithm for Classification of White Blood Cells in Peripheral Smear Images This article addresses a new method for the classification of white blood cells (WBCs) using image processing techniques and machine learning methods. The proposed method consists of three steps: detecting the nucleus and cytoplasm, extracting features, and classification. At first, a new algorithm is designed to segment the nucleus. For the cytoplasm to be detected, only a part of it located inside the convex hull of the nucleus is involved in the process. This attitude helps us overcome the difficulties of segmenting the cytoplasm. In the second phase, three shapes and four novel color features are devised and extracted. Finally, by using an SVM model, the WBCs are classified. The segmentation algorithm can detect the nucleus with a dice similarity coefficient of 0.9675. The proposed method can categorize WBCs in Raabin-WBC, LISC, and BCCD datasets with accuracies of 94.65 %, 92.21 %, and 94.20 %, respectively. It is worth mentioning that the hyperparameters of the classifier are fixed only with the Raabin-WBC dataset, and these parameters are not readjusted for LISC and BCCD datasets. The obtained results demonstrate that the proposed method is robust, fast, and accurate.