Common variants in breast cancer risk loci predispose to distinct tumor subtypes BackgroundGenome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER), but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear.

MethodsAmong 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes.

ResultsEighty-five of 173 variants were associated with at least one tumor feature (false discovery rate <5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at P<0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions.

ConclusionThis report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.
A multi-dimensional integrative scoring framework for predicting functional variants in the human genome Attempts to identify and prioritize functional DNA elements in coding and noncoding regions, particularly through use of in silico functional annotation data,continue to increase in popularity. However, specific functional roles may vary widely from one variant to another, making it challenging to summarize different aspects of variant function. Here we propose Multi-dimensional Annotation Class Integrative Estimation (MACIE), an unsupervised multivariate mixed model framework capable of integrating annotations of diverse origin to assess multi-dimensional functional roles for both coding and noncoding variants. Unlike existing one-dimensional scoring methods, MACIE views variant functionality as a composite attribute encompassing multiple characteristics, and estimates the joint posterior functional probability vector of each genomic position, a quantity that offers richer and more interpretable information in the presence of multiple aspects of functionality. Applied to a variety of independent coding and non-coding datasets, MACIE demonstrates powerful and robust performance in discriminating between functional and non-functional variants. We also show an application of MACIE to fine-mapping using lipids GWAS summary statistics data from the European Network for Genetic and Genomic Epidemiology Consortium.
Structural effects driven by rare point mutations in amylin hormone, the type II diabetes-associated peptide BackgroundAmylin is a 37-amino-acid peptide hormone co-secreted with insulin, which participates in glucose homeostasis. This hormone is able to aggregate in a {beta}-sheet conformation and deposit in islet amyloids, a hallmark in type II diabetes. Since amylin is a gene-encoded hormone, this peptide has variants caused by point mutations that can impact its functions.

MethodsHere, we analyzed the structural effects caused by S20G and G33R point mutations which, according to the 1000 Genomes Project, have frequency in East Asian and European populations, respectively. The analyses were performed by means of aggrescan server, SNP functional effect predictors, and molecular dynamics.

ResultsWe found that both mutations have aggregation potential and cause changes in the monomeric forms when compared with wild-type amylin. Furthermore, comparative analyses with pramlintide, an amylin drug analogue, allowed us to infer that second -helix maintenance may be related to the aggregation potential.

ConclusionsThe S20G mutation has been described as pathologically related, which is in agreement with our findings. In addition, our data suggest that the G33R mutation might have a deleterious effect. The data presented here also provide new therapy opportunities, whether for creating more effective drugs for diabetes or implementing specific treatment for patients with these mutations.

General SignificanceOur data could help to better understand the impact of mutations on the wild-type amylin sequence, as a starting point for the evaluation and characterization of other variations. Moreover, these findings could improve the health of patients with type II diabetes. Feline tooth resorption findings from a questionnaire survey on 8115 Finnish cats and a case-control study for a more defined population drawn from the survey Tooth resorption (TR) is one of the most common dental diseases in cats. It is a painful disease characterised by progressive dental destruction, which eventually results in loss of teeth. The aetiology of the disease remains unclear, but associations with old age, breed, other oral and dental diseases, and certain environmental factors have been suspected. We collected health and environmental information of 8115 Finnish cats from 41 breeds through an online survey for cat owners to study TR in Finland with an observational study. The frequency of TR that the owner reported as veterinarian-diagnosed was 3.9% in the entire data (316/8115) and 15% in cats that the owner reported as having veterinarian-diagnosed oral or dental disease (316/2070). We determined the risk factors and breed variation of feline TR by logistic regression with the entire data and a case-control study in a more defined population drawn from the entire data. In the entire data, we defined those cats that the owners had reported as being diagnosed with TR by a veterinarian as TR-cats and compared them to all other cats. In the case-control study, we included only cats with owner-reported veterinarian-diagnosed oral and dental diseases. We defined cats with TR-diagnosis made under sedation as TR-cases. Other veterinarian-diagnosed TR-cats were excluded. Cats with other oral and dental diseases were used as controls. The model with the entire population included 308 TR-cats compared with 7508 cats, and the case-control study included 198 TR-cats and 1683 controls. Although the diagnoses were reportedly made by a veterinarian, the owner answered the survey and might have understood, remembered, or interpreted the diagnoses in a way that could affect the results. The results indicated an increasing risk of TR with age. Dental calculus, gingivitis, and periodontitis were potential risk factors for TR. These findings and the interaction between dental calculus and gingivitis suggest that inflammatory changes caused by dental calculus increase the risk of TR. Keeping food available constantly was a potential protective factor. We found that Cornish Rex, European, and Ragdoll breeds were at higher risk for TR. The observed differences between breeds highlight a genetic contribution. First detection of herpesvirus and prevalence of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. Dissociating the Neural Correlates of Planning and Executing Tasks with Nested Task Sets Task processing (e.g., the preparation and execution of responses) and task representation (e.g., the activation and maintenance of stimulus-response and context information) are two facets of cognitive control supported by lateral frontal cortex (LFC). However, the mechanistic overlap (or distinction) between these two facets is unknown. We explored this by combining a complex task mapping with a pre-cueing procedure. Participants made match/non-match judgments on pairs of stimuli during fMRI recording. Pre-cues on each trial gave variable amounts of information to the participant in anticipation of the stimulus. Our results demonstrated that regions throughout LFC were more active at the stimulus (when responses could be executed) than at the cue (when they could only be prepared), indicating that they supported execution of the task agnostic to the specific task representation. A subset of regions in left caudal LFC showed increased activity with more cue information at the cue and the reverse at the stimulus, suggesting their involvement in reducing uncertainty within the task representation. These results suggest that one component of task processing is preparing and executing the task according to the relevant representation, confined to left caudal LFC, while non-representational functions that occur primarily during execution are supported by different regions throughout the rest of LFC. We further conducted an exploratory investigation of connectivity between the two groups of regions in this study and their potential relationship to the fronto-parietal and cingulo-opercular networks. Regions with both patterns of activity appear to be part of the fronto-parietal network. Naturally enhanced neutralizing breadth to SARS-CoV-2 after one year Over one year after its inception, the coronavirus disease-2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) remains difficult to control despite the availability of several excellent vaccines. Progress in controlling the pandemic is slowed by the emergence of variants that appear to be more transmissible and more resistant to antibodies1,2. Here we report on a cohort of 63 COVID-19-convalescent individuals assessed at 1.3, 6.2 and 12 months after infection, 41% of whom also received mRNA vaccines3,4. In the absence of vaccination antibody reactivity to the receptor binding domain (RBD) of SARS-CoV-2, neutralizing activity and the number of RBD-specific memory B cells remain relatively stable from 6 to 12 months. Vaccination increases all components of the humoral response, and as expected, results in serum neutralizing activities against variants of concern that are comparable to or greater than neutralizing activity against the original Wuhan Hu-1 achieved by vaccination of naive individuals2,5-8. The mechanism underlying these broad-based responses involves ongoing antibody somatic mutation, memory B cell clonal turnover, and development of monoclonal antibodies that are exceptionally resistant to SARS-CoV-2 RBD mutations, including those found in variants of concern4,9. In addition, B cell clones expressing broad and potent antibodies are selectively retained in the repertoire over time and expand dramatically after vaccination. The data suggest that immunity in convalescent individuals will be very long lasting and that convalescent individuals who receive available mRNA vaccines will produce antibodies and memory B cells that should be protective against circulating SARS-CoV-2 variants. The N terminus of adhesion G protein-coupled receptor GPR126/ADGRG6 as allosteric force integrator The adhesion G protein-coupled receptor (aGPCR) GPR126/ADGRG6 plays an important role in several physiological functions, such as myelination or peripheral nerve repair. This renders the receptor an attractive pharmacological target. GPR126 is a mechano-sensor that translates binding of extracellular matrix (ECM) molecules to its N terminus into a metabotropic intracellular signal. To date, the structural requirements and the character of the forces needed for this ECM-mediated receptor activation are largely unknown.

In this study we provide this information by combining classic second messenger detection with single cell atomic force microscopy. We establish a monoclonal antibody targeting the N terminus to stimulate GPR126 and compare it to the activation through its known ECM ligands collagen IV and laminin 211. As each ligand uses a distinct mode of action, the N terminus can be viewed as an allosteric module that can fine-tune receptor activation in a context-specific manner. A counter-intuitive method for slowing down an analyte in a nanopore and its application to polymer sequencing and analysis An unresolved problem in nanopore sensing is the high translocation speed ([~]10-100 monomers/s) of an analyte (nucleotide, DNA, amino acid (AA), peptide) through the pore. Here a method based on reversing the pore voltage and changing the solution pH is described. A simplified Fokker-Planck model shows mean translocation times of 1-10 ms in a nanopore of length 10 nm. Simulations show that a positive-negative voltage profile can trap an analyte for [~]1 ms. This method can be used for free nucleotides, single AAs, oligonucleotides, and oligopeptides. Its applicability to existing nanopore sensing and sequencing techniques and implementation issues are discussed. 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. A protease-activatable luminescent biosensor andreporter cell line for authentic SARS-CoV-2 infection Efforts to define serological correlates of protection against COVID-19 have been hampered by the lack of a simple, scalable, standardised assay for SARS-CoV-2 infection and antibody neutralisation. Plaque assays remain the gold standard, but are impractical for high-throughput screening. In this study, we show that expression of viral proteases may be used to quantitate infected cells. Our assays exploit the cleavage of specific oligopeptide linkers, leading to the activation of cell-based optical biosensors. First, we characterise these biosensors using recombinant SARS-CoV-2 proteases. Next, we confirm their ability to detect viral protease expression during replication of authentic virus. Finally, we generate reporter cells stably expressing an optimised luciferase-based biosensor, enabling viral infection to be measured within 24 h in a 96- or 384-well plate format, including variants of concern. We have therefore developed a luminescent SARS-CoV-2 reporter cell line, and demonstrated its utility for the relative quantitation of infectious virus and titration of neutralising antibodies.

Author summaryTechniques for measuring infection with SARS-CoV-2 in the laboratory are laborious and time-consuming, and different laboratories use different approaches. There is therefore no generally agreed way to quantitate neutralising antibodies against SARS-CoV-2, which block infection with the virus and protect people from COVID-19. In this study, we describe a new way to measure SARS-CoV-2 infection, which is much simpler and faster than existing methods. It relies on the production of a specific protease enzyme by the virus, which is able to cleave and activate an engineered protein biosensor in infected cells. This biosensor emits light in the presence of viral infection, and the amount of light released is used as a readout for the amount of infectious SARS-CoV-2 present. The signal is very sensitive, so the number of infected cells required is very small, and the method can be scaled-up to test many samples at once. In particular, we demonstrate how it can be used to detect different variants of SARS-CoV-2, and quantitate neutralising antibodies against these viruses. Uncovering Differences in Hydration Free Energies and Structures for Model Compound Mimics of Charged Sidechains of Amino Acids Free energies of hydration are of fundamental interest for modeling and understanding conformational and phase equilibria of macromolecular solutes in aqueous phases. Of particular relevance to systems such as intrinsically disordered proteins are the free energies of hydration and hydration structures of model compounds that mimic charged sidechains of Arg, Lys, Asp, and Glu. Here, we deploy a Thermodynamic Cycle based Proton Dissociation (TCPD) approach in conjunction with data from direct measurements to obtain estimates for the free energies of hydration for model compounds that mimic the sidechains of Arg+, Lys+, Asp-, and Glu-. Irrespective of the choice made for the hydration free energy of the proton, the TCPD approach reveals clear trends regarding the free energies of hydration for Arg+, Lys+, Asp-, and Glu-. These trends include asymmetries between the hydration free energies of acidic (Asp- and Glu-) and basic (Arg+ and Lys+) residues. Further, the TCPD analysis, which relies on a combination of experimental data, shows that the free energy of hydration of Arg+ is less favorable than that of Lys+. We sought a physical explanation for the TCPD derived trends free energies of hydration. To this end, we performed temperature dependent calculations of free energies of hydration and analyzed hydration structures from simulations that use the polarizable AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications) forcefield and water model. At 298 K, the AMOEBA model generates estimates of free energies of hydration that are consistent with TCPD values with a free energy of hydration for the proton of {approx} -259 kcal / mol. Analysis of temperature dependent simulations leads to a structural explanation for the observed differences in free energies of hydration of ionizable residues and reveals that the heat capacity of hydration is positive for Arg+ and Lys+ and negative for Asp- and Glu-. Dynamic Phosphorylation of RGS Provides Spatial Regulation of MAP Kinase And Promotes Completion of Cytokinesis during the Yeast Pheromone Response Yeast use a G-protein coupled receptor (GPCR) signaling pathway to detect mating pheromone, arrest in G1, and direct polarized growth towards the potential mating partner. The primary negative regulator of this pathway is the regulator of G-protein signaling (RGS), Sst2, which induces G GTPase activity and subsequent inactivation of all downstream signaling including a MAPK cascade. The MAPK Fus3 phosphorylates the RGS in response to pheromone, but the role of this modification is unknown. We set out to examine the role of RGS phosphorylation during the pheromone response. We found that RGS phosphorylation peaks early in the pheromone response and diminishes RGS localization to the polarity site and focuses MAPK complexes there. At later time points, RGS is predominantly unphosphorylated, which promotes RGS localization to the polar cap and broadens the distribution of MAPK complexes relative to the Cdc42 polarity machinery. Surprisingly, we found that phosphorylation of the RGS is required for the completion of cytokinesis prior to pheromone induced growth. The completion of cytokinesis in the presence of pheromone is promoted by the formin Bnr1 and the kelch-repeat protein, Kel1, both proteins previously found to interact with the RGS. Methylome dynamic upon proteasome inhibition by the Pseudomonas syringae 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. Inhibition of BCMA signaling alters protein translation and Multiple Myeloma growth Recent T Cell therapies have been effective in the treatment of hematological cancers. However, immunotoxicity and treatment relapse pose significant clinical challenges. Here, we revealed distinctive requirement for neutralizing TNF receptor ligands APRIL and BAFF in MM and DLBCL. Furthermore, we investigated the use of BCMA decoy receptor (sBCMA-Fc) as a therapeutic inhibitor of ARPIL and BAFF. While wild-type sBCMA-Fc successfully blocked APRIL signaling with picomolar binding affinity, inhibiting tumor growth in MM models, it lacked efficacy in inhibiting DLBCL progression due to its weak binding for BAFF. To expand the therapeutic utility of sBCMA-Fc, using a protein engineering approach, we generated an affinity enhanced mutant sBCMA-Fc fusion molecule (sBCMA-Fc V3) with 4-folds and 500-folds enhancement in binding to APRIL and BAFF respectively. The sBCMA-Fc V3 clone significantly enhanced antitumor activity against both MM and DLBCL. Importantly, sBCMA-Fc V3 was proven to be a viable clinical candidate by showing adequate toxicity profile and on-target mechanism of action in nonhuman primates.

SUMMARYThis study demonstrates the dichotomous function of APRIL and BAFF in MM and DLBCL, that can be safely targeted by an engineered fusion protein designed to trap APRIL and BAFF with ultra-high binding affinity. Learning and motor inhibitory control in crows and domestic chickens Cognitive abilities allow animals to navigate through complex, fluctuating environments. In the present study, we tested the performance of a captive group of eight crows, Corvus corone, and ten domestic chickens, Gallus gallus domesticus, in the cylinder task, as a test of motor inhibitory control, and reversal learning as a measure of learning ability and behavioural flexibility. Four crows and nine chickens completed the cylinder task, eight crows and six chickens completed the reversal learning experiment. Crows performed better in the cylinder task compared to chickens. In the reversal learning experiment, species did not significantly differ in the number of trials until learning criterion was reached. Performance in the reversal learning experiment did not correlate with performance in the cylinder task in chickens. Our results suggest crows to possess better motor inhibitory control compared to chickens. In contrast, learning performance in a reversal learning task did not differ between the species, indicating similar levels of behavioural flexibility. Interestingly, we describe notable individual differences in performance. We stress the importance not only to compare cognitive performance between species but also between individuals of the same species when investigating the evolution of cognitive skills. Structural and Functional Analyses Explain Pea KAI2 Receptor Diversity and Reveal Stereoselective Catalysis During Signal Perception KAI2 are plant /{beta} hydrolase receptors, which perceive smoke-derived butenolide signals (karrikins) and putative endogenous, yet unidentified phytohormones (KAI2-ligands, KLs). The number of functional KAI2 receptors varies among plant species. It has been suggested that KAI2 gene duplication and sub-functionalization plays an adaptative role for diverse environments or ligand diversification by altering the receptor responsiveness to specific KLs. Legumes represent one of the largest families of flowering plants and contain many essential agronomic crops. Prior to legume diversification, KAI2 underwent duplication, resulting in KAI2A and KAI2B. Integrating plant genetics, ligand perception and enzymatic assays, and protein crystallography, we demonstrate that Pisum sativum KAI2A and KAI2B act as receptors and enzymes with divergent ligand stereoselectivity. KAI2B has a stronger affinity than KAI2A towards the KAI2-ligand (-)-GR24 and remarkably hydrolyses a broader range of substrates including the strigolactone-like isomer (+)-GR24. We determine the crystal structures of PsKAI2B in apo and butenolide-bound states. The biochemical and structural analyses as well as recorded mass spectra of KAI2s reveal a transient intermediate on the catalytic serine and a stable adduct on the catalytic histidine, further illuminating the role of KAI2 not only as receptors but also as bona fide enzymes. Our work uncovers the stereoselectivity of ligand perception and catalysis by evolutionarily diverged KAI2 receptors in KAR/KL signaling pathways and proposes adaptive sensitivity to KAR/KL and strigolactone phytohormones by KAI2B. A Wnt/β-catenin signaling driven adipocyte population required for beiging and glucose homeostasis Wnt/{beta}-catenin signaling has been well established as a potent inhibitor of adipogenesis. Here, we identified a population of adipocytes that exhibit persistent activity of Wnt/{beta}-catenin signaling, as revealed by the Tcf/Lef-Gfp reporter allele, in embryonic and adult mouse fat depots, named as Wnt+ adipocytes. We showed that the {beta}-catenin-mediated signaling activation in these cells is Wnt ligand- and receptor-independent but relies on AKT/mTOR pathway and is essential for cell survival. Such adipocytes are distinct from classical ones in transcriptomic and genomic signatures and can be induced from various sources of mesenchymal stromal cells including human cells. Genetic lineage-tracing and targeted cell ablation studies revealed that these adipocytes convert into beige adipocytes directly and are also required for beige fat recruitment under thermal challenge, demonstrating both cell autonomous and non-cell autonomous roles in adaptive thermogenesis. Furthermore, mice bearing targeted ablation of these adipocytes exhibited glucose intolerance, while mice receiving exogenously supplied such cells manifested enhanced glucose utilization. Our studies uncover a unique adipocyte population in regulating beiging in adipose tissues and systemic glucose homeostasis. A parametric logistic equation with light flux and medium concentration for cultivation planning of microalgae Microalgae are considered to be promising producers of bioactive chemicals, feeds, and fuels from carbon dioxide by photosynthesis. Thus, the prediction of microalgal growth profiles is important for the planning of cost-effective and sustainable cultivation-harvest cycles. This paper proposes a mathematical model capable of predicting the effect of light flux into culture and medium concentration on the growth profiles of microalgae by incorporating these growth-limiting factors into the logistic equation. The specific form of the equation is derived based on the experimentally measured growth profiles of Monoraphidium sp., a microalgal strain isolated by the authors, under 16 conditions consisting of the combinations of incident light fluxes into culture and initial medium concentrations. Using a cross-validation method, it is shown that the proposed model has the ability to predict necessary incident light flux into culture and initial medium concentration for harvesting target biomass at target time. Finally, the model-guided cultivation planning is performed and is evaluated by comparing the result with the experimental data. Integrated annotation and analysis of genomic features reveal new types of functional elements and large-scale epigenetic phenomena in the developing zebrafish Zebrafish, a popular model for embryonic development and for modelling human diseases, has so far lacked a systematic functional annotation programme akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created the first central repository to store and process zebrafish developmental functional genomic data. Our Data Coordination Center (https://danio-code.zfin.org) combines a total of 1,802 sets of unpublished and reanalysed published genomics data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements in development, including novel classes with distinct features dependent on their activity in time and space. We delineated the distinction between regulatory elements active during zygotic genome activation and those active during organogenesis, identifying new aspects of how they relate to each other. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predict functional relationships between them beyond sequence similarity, extending the utility of zebrafish developmental genomics to mammals. Bulk-assembly of monodisperse coacervates and giant unilamellar vesicles with programmable hierarchical complexity Compartmentalization is an essential step for the emergence of protocellular life on the early earth and the endeavor to bottom-up construct synthetic cells in the lab. Among manifold natural strategies for compartmentalization, assembly of lipid vesicles or coacervate droplets with resemblance to cell envelope or cytoplasm have emerged as two dominant paradigms. However, the spontaneous assembly of lipid vesicles or coacervates generally results in polydisperse microcompartments and lacks the ability to efficiently integrate functional building blocks towards defined complexity, which give rise to the difficulty to achieve the order and complexity towards the scope of that of life via spontaneous assembly. Herein, we show that the interplay of coacervates and colloidal particles provides a bulk-assembly approach to form monodisperse coacervates and giant unilamellar vesicles (GUVs) with programmable hierarchical complexity. The hierarchical complexity of the synthetic cell microcompartments can be delicately engineered from the control of the spatial or temporal organization of components in entity unit to manifold forms of interconnected synthetic cell consortia, achieving a compartment hierarchy analogous to the organelle-cell-tissue structure. The increase in complexity gives rise to emergent properties, for instance, collective morphology evolution of coacervate assemblages, gated permeability of GUVs without sophisticated protein machinery, and remarkable structural and functional stability of GUVs at extreme conditions. This work paves an unprecedented step to form monodisperse, hierarchical coacervates and GUVs via bulk-assembly, which provides a prebiotically plausible approach for the proper integration of inanimate matter towards the order and complexity of life and monodisperse, hierarchical microdroplets for follow-up scientific and technological applications. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. An exon-intron split framework to prioritize miRNA-driven post-transcriptional regulatory signals and its application to study energy homeostasis in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. The emergence of highly fit SARS-CoV-2 variants accelerated by epistasis and caused by recombination The SARS-CoV-2 pandemic recently entered an alarming new phase with the emergence of the variants of concern (VOC) and understanding their biology is paramount to predicting future ones. Current efforts mainly focus on mutations in the spike glycoprotein (S), but changes in other regions of the viral proteome are likely key. We analyzed more than 900,000 SARS-CoV-2 genomes with a computational systems biology approach including a haplotype network and protein structural analyses to reveal lineage-defining mutations and their critical functional attributes. Our results indicate that increased transmission is promoted by epistasis, i.e., combinations of mutations in S and other viral proteins. Mutations in the non-S proteins involve immune-antagonism and replication performance, suggesting convergent evolution. Furthermore, adaptive mutations appear in geographically disparate locations, suggesting that either independent, repeat mutation events or recombination among different strains are generating VOC. We demonstrate that recombination is a stronger hypothesis, and may be accelerating the emergence of VOC by bringing together cooperative mutations. This emphasizes the importance of a global response to stop the COVID-19 pandemic. Discovery of highly potent pancoronavirus fusion inhibitors that also effectively inhibit COVID-19 variants from the UK (Alpha), South Africa (Beta), and India (Delta) We report here the discovery of several highly potent small molecules that showed low nM potency against SARS-CoV (IC50: as low as 13 nM), SARS-CoV-2 (IC50: as low as 23 nM), and MERS-CoV (IC50: as low as 76 nM) in pseudovirus based assays with excellent selectivity indices (SI: as high as > 5000) demonstrating their pancoronavirus inhibition. Some compounds also show 100% inhibition of CPE (IC100) at 1.25 {micro}M against an authentic SARS-CoV-2 (US_WA-1/2020). Furthermore, the most active inhibitors also potently inhibited variants of concerns (VOCs), such as the UK (B.1.1.7), South Africa (B.1.351), and Delta variant (B.1.617.2), originated in India. We confirmed that one of the potent inhibitors binds to the prefusion spike protein trimer of SARS-CoV-2 by SPR. Besides, we showed that they inhibit virus-mediated cell-cell fusion. The ADME data of one of the most active inhibitors, NBCoV1, show drug-like properties. In vivo PK of NBCoV1 in rats demonstrated excellent half-life (t1/2) of 11.3 h, mean resident time (MRT) of 14.2 h, and oral bioavailability. We expect the lead inhibitors to pave the way for further development to preclinical and clinical candidates. Scaling effects of temperature on parasitism from individuals to host-parasite systems Parasitism is expected to change in a warmer future, but whether warming leads to substantial increases in parasitism remains unclear. Understanding how warming effects on parasitism in individual hosts (e.g., parasite load) translate to effects on population-level parasitism (e.g., prevalence, R0) remains a major knowledge gap. We analyzed the temperature dependence of parasitism at both host and population levels in thirteen empirical vector-borne host-parasite systems and found a strong, significant positive correlation between the thermal optima of individual- and population-level parasitism. We also found a significant, positive correlation in eleven environmentally-transmitted parasite systems, though several of these systems exhibited thermal optima >5{o}C apart between individual and population levels. Similarly, parasitism thermal optima were close to host performance thermal optima in vector-borne systems but not in environmentally-transmitted systems. We then adapted and simulated simple models for both transmission modes and found a similar pattern to the empirical systems: thermal optima in vector-borne systems were more strongly correlated across scales compared to environmentally-transmitted systems. Generally, our results suggest that information on the temperature-dependence, and specifically the thermal optimum, at either the individual- or population-level should provide a useful--though not quantitatively exact--baseline for predicting temperature dependence at the other level, especially in vector-borne parasite systems. Environmentally-transmitted parasitism may operate by a different set of rules, in which temperature-dependence is decoupled in some systems, requiring the need for trait-based studies of temperature-dependence at individual and population levels.