The quantitative genetics of the prevalence of infectious diseases: hidden genetic variation due to Indirect Genetic Effects dominates heritable variation and response to selection Infectious diseases have profound effects on life, both in nature and agriculture. Despite the availability of well-established epidemiological theory, however, a quantitative genetic theory of the host population for the endemic prevalence of infectious diseases is almost entirely lacking. While several studies have demonstrated the relevance of the transmission dynamics of infectious diseases for heritable variation and response to selection of the host population, our current theoretical framework of quantitative genetics excludes these dynamics. As a consequence, we do not know which genetic effects of the host population determine the prevalence of an infection, and have no concepts of breeding value and heritable variation for endemic prevalence.

Here we integrate quantitative genetics and epidemiology, and propose a quantitative genetic theory for R0 and for the endemic prevalence of an infectious disease. We first identify the genetic factors that determine the prevalence of an infection, using an approach founded in epidemiological theory. Subsequently we investigate the population level consequences of individual genetic variation, both for R0 and for the endemic prevalence. Next, we present expressions for the breeding value and heritable variation, for both prevalence and individual binary disease status, and show that these parameters depend strongly on the level of the prevalence. Results show that heritable variation for endemic prevalence is substantially greater than currently believed, and increases when prevalence approaches zero, while heritability of individual disease status goes to zero. As a consequence, response of prevalence to selection accelerates considerably when prevalence goes down, in contrast to predictions from classical theory. Finally, we show that most of the heritable variation for the endemic prevalence of an infection is hidden due to indirect genetic effects, suggesting a key role for kin-group selection both in the evolutionary history of current populations and for genetic improvement strategies in animals and plants.
Birth Order Differences in Education Originate in Post-Natal Environments Siblings share many environments and much of their genetics. Yet, siblings turn out different. Intelligence and education are influenced by birth order, with earlier-born siblings outperforming later-borns. We investigate whether birth order differences in education are caused by biological differences present at birth, that is, genetic differences or in-utero differences. Using family data that spans two generations, combining registry, survey, and genotype information, this study is based on the Norwegian Mother, Father and Child Cohort Study (MoBa). We show that there are no genetic differences by birth order as captured by polygenic scores (PGSs) for educational attainment. Furthermore, we show that earlier-born have lower birth weight than later-born, indicating worse uterine environments. Educational outcomes are still higher for earlier-born children when we adjust for PGSs and in utero variables, indicating that birth order differences arise post-natally. Finally, we consider potential environmental influences, such as differences according to maternal age, parental educational attainment, and sibling genetic nurture. We show that birth order differences are not biological in origin, but pinning down their specific causes remains elusive.
Effect of a low concentration in chloride, sodium and potassium on oocyte maturation, oocyte hydration, ovulation and egg quality in rainbow trout. Water salinity is an important environmental factor known to have detrimental effects on salmonid reproduction, mostly when migrating female broodfish are held in sea water. In contrast, data obtained in freshwater are scarce and the impact of low water salinity, and more specifically of low water concentrations in sodium, chloride and potassium, during reproduction in freshwater is currently unknown. For this reason, and because ion and water fluxes are critical for the final steps of the female gamete formation, including oocyte hydration and ovulation, the aim of the present study was to investigate the impact of low salinity water on final oocyte maturation, ovulation and, ultimately, on egg quality, using rainbow trout as a physiological model and relevant aquaculture species.

Fish from the same commercial strain were raised either in a site characterized by low concentrations of Na+, K+, and Cl- ions in the water or in a closely located control site exhibiting higher concentration in these elements. Egg quality and duration of final oocyte maturation were investigated using innovative phenotyping tools such as automatic assessment of egg viability using the VisEgg system and non-invasive echograph-based monitoring of final oocyte maturation duration, respectively. Oocyte hydration during final oocyte maturation and after ovulation was also investigated. Finally, molecular phenotyping was performed using real-time PCR-based monitoring of several key players of final oocyte maturation and ovulation associated with ion and water transport, inflammation, proteolytic activity, and coagulation. Oocyte hydration and gene expression data were analyzed in the light of the duration of final oocyte maturation.

Here we show that low water salinity (i.e., low water concentration in chloride, sodium and potassium) negatively influences final oocyte maturation, ovulation and, ultimately, egg quality. Low water salinity triggered delayed ovulation and lower oocyte viability. When investigating the impact of low water salinity on final oocyte maturation duration, individuals presenting the most severe phenotypes exhibited impaired oocyte hydration and abnormally reduced gene expression levels of several key players of the ovulatory process. While the under expression of water (i.e., aquaporins) and ion (i.e., solute carriers) transporters is consistent with impaired oocyte hydration, our observations also indicate that the overall ovulatory gene expression program is disrupted. Our results raise the question of the mechanisms underlying the negative influence of low salinity water on the dynamics of the preovulatory phase, on the control of the oocyte homeostasis, including hydration, and on the overall success of the maturation-ovulation process.

HighlightsO_LILow water salinity impairs final oocyte maturation and egg quality in rainbow trout
C_LIO_LILow water salinity induces delayed ovulation and impaired oocyte hydration
C_LIO_LILow water salinity induces a dysregulation of several key ovulatory genes
C_LIO_LIMonitoring of final oocyte maturation can be performed using ultrasound staging
C_LI TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Protein cost minimization promotes the emergence of coenzyme redundancy Coenzymes distribute a variety of chemical moieties throughout cellular metabolism, participating in group (e.g., phosphate, acyl) and electron transfer. For a variety of reactions requiring acceptors or donors of specific resources, there often exist degenerate sets of molecules (e.g., NAD(H) and NADP(H)) that carry out similar functions. Although the physiological roles of various coenzyme systems are well established, it is unclear what selective pressures may have driven the emergence of coenzyme redundancy. Here we use genome-wide metabolic modeling approaches to decompose the selective pressures driving enzymatic specificity for either NAD(H) or NADP(H) in the metabolic network of Escherichia coli. We found that few enzymes are thermodynamically constrained to using a single coenzyme, and in principle, a metabolic network relying on only NAD(H) is feasible. However, structural and sequence analyses revealed widespread conservation of residues that retain selectivity for either NAD(H) or NADP(H), suggesting that additional forces may shape specificity. Using a model accounting for the cost of oxidoreductase enzyme expression, we found that coenzyme redundancy universally reduces the minimal amount of protein required to catalyze coenzyme-coupled reactions, inducing individual reactions to strongly prefer one coenzyme over another when reactions are near thermodynamic equilibrium. We propose that protein minimization generically promotes coenzyme redundancy, and that coenzymes typically thought to exist in a single pool (e.g., CoA) may exist in more than one form (e.g., dephospho-CoA).

Significance statementMetabolism relies on a small class of molecules (coenzymes) that serve as universal donors and acceptors of key chemical groups and electrons. Although metabolic networks crucially depend on structurally redundant coenzymes (e.g., NAD(H) and NADP(H)) associated with different enzymes, the criteria that led to the emergence of this redundancy remain poorly understood. Our combination of modeling, and structural and sequence analysis indicates that coenzyme redundancy is not essential for metabolism, but rather an evolved strategy promoting efficient usage of enzymes when biochemical reactions are near equilibrium. Our work suggests that early metabolism may have operated with fewer coenzymes, and that adaptation for metabolic efficiency may have driven the rise of coenzyme diversity in living systems. Improving GWAS discovery and genomic prediction accuracy in Biobank data Genetically informed and deep-phenotyped biobanks are an important research resource. The cost of phenotyping far outstrips that of genotyping, and therefore it is imperative that the most powerful, versatile and efficient analysis approaches are used. Here, we apply our recently developed Bayesian grouped mixture of regressions model (GMRM) in the UK and Estonian Biobanks and obtain the highest genomic prediction accuracy reported to date across 21 heritable traits. On average, GMRM accuracies were 15% (SE 7%) greater than prediction models run in the LDAK software with SNP annotation marker groups, 18% (SE 3%) greater than a baseline BayesR model without SNP markers grouped into MAF-LD-annotation categories, and 106% (SE 9%) greater than polygenic risk scores calculated from mixed-linear model association (MLMA) estimates. For height, the prediction accuracy R2 was 47% in a UK Biobank hold-out sample, which was 76% of the estimated [Formula]. We then extend our GMRM prediction model to provide MLMA SNP marker estimates for GWAS discovery, which increased the independent loci detected to 7,910 in unrelated UK Biobank individuals, as compared to 5,521 from BoltLMM and 5,727 from Regenie, a 43% and 38% increase respectively. The average{chi} 2 value of the leading markers was 34% (SE 5.11) higher for GMRM as compared to Regenie, and increased by 17% for every 1% increase in prediction accuracy gained over a baseline BayesR model across the traits. Thus, we show that modelling genetic associations accounting for MAF and LD differences among SNP markers, and incorporating prior knowledge of genomic function, is important for both genomic prediction and for discovery in large-scale individual-level biobank-scale studies. PRDM9 losses in vertebrates are coupled to those of paralogs ZCWPW1 and ZCWPW2 In most mammals and likely throughout vertebrates, the gene PRDM9 specifies the locations of meiotic double strand breaks; in mice and humans at least, it also aids in their repair. For both roles, many of the molecular partners remain unknown. Here, we take a phylogenetic approach to identify genes that may be interacting with PRDM9, by leveraging the fact that PRDM9 arose before the origin of vertebrates, but was lost many times, either partially or entirely--and with it, its role in recombination. As a first step, we characterize PRDM9 domain composition across 446 vertebrate species, inferring at least thirteen independent losses. We then use the interdigitation of PRDM9 orthologs across vertebrates to test whether it co-evolved with any of 241 candidate genes co-expressed with PRDM9 in mice or associated with recombination phenotypes in mammals. Accounting for the phylogenetic relationship among species, we find two genes whose presence and absence is unexpectedly coincident with that of PRDM9: ZCWPW1, which was recently shown to facilitate double strand break repair, and its paralog ZCWPW2, as well as more tentative evidence for TEX15 and FBXO47. ZCWPW2 is expected to be recruited to sites of PRDM9 binding; its tight coevolution with PRDM9 across vertebrates suggests that it is a key interactor, with a role either in recruiting the recombination machinery or in double strand break repair.

Author SummaryOur understanding of meiotic recombination in mammals has seen great progress over the past 15 years, spurred in part by the convergence of lines of evidence from molecular biology, statistical genetics and evolutionary biology. We now know that in most mammals and likely in many vertebrates, the gene PRDM9 specifies the location of meiotic double strand breaks and that in mice and humans at least, it also aids in their repair. For both roles, however, many of the molecular partners remain unknown. To search for these, we take a phylogenetic approach, leveraging the fact that the complete PRDM9 has been lost at least thirteen times in vertebrates and thus that its presence is interdigitated across species. By this approach, we identify two genes whose presence or absence across vertebrates is coupled to the presence or absence of PRDM9, ZCWPW1 and ZCWPW2, as well as two genes for which the evidence is weaker, TEX15 and FBXO47. ZCWPW1 was recently shown to be recruited to sites of PRDM9 binding and to aid in the repair of double strand breaks. ZCWPW2 is likely recruited to sites of PRDM9 binding as well; its tight coevolution with PRDM9 across vertebrates suggests that it plays an important role either in double strand break formation, potentially as the missing link that recruits the recombination machinery to sites of PRDM9 binding, or in double strand break repair. Cellular-level phenotyping of tumor-immune microenvironment (TiME) in patients in vivo reveals distinct inflammation and endothelial anergy signatures Immunotherapies have shown unprecedented clinical benefits in several malignancies1-3. However, clinical responses remain variable and unpredictable, indicating the need to develop predictive platforms that can improve patient stratification4. Phenotyping of tumors into hot, altered, or cold5 based on assessment of only T-cell infiltration in static tumor biopsies provides suboptimal prediction of immunotherapy response6,7. In vivo dynamic mechanisms within the tumor microenvironment such as tumor angiogenesis and leukocyte trafficking5,8,9 also play a central role in modulating anti-tumor immunity and therefore immunotherapy response. Here, we report novel tumor immune microenvironment (TiME) phenotyping in vivo in patients with non-invasive spatially-resolved cellular-level imaging based on endogenous contrast. Investigating skin cancers as a model, with reflectance confocal microscopy (RCM) imaging10, we determined four major phenotypes with variable prevalence of vasculature (Vasc) and inflammation (Inf) features: VaschiInfhi, VaschiInflo, VascloInfhi and Vascmed/hiInflo. The VaschiInfhi phenotype correlates with high immune activation, exhaustion, and vascular signatures while VaschiInflo with endothelial anergy and immune exclusion. Automated quantification of TiME features demonstrates moderate-high accuracy and correlation with corresponding gene expression. Prospectively analyzed response to topical immunotherapy show highest response in VascloInfhi, and reveals the added value of vascular features in predicting treatment response. Our novel in vivo cellular-level imaging and phenotyping approach can potentially advance our fundamental understanding of TiME, develop robust predictors for immunotherapy outcomes and identify novel targetable pathways in future. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. A network approach to genetic circuit designs As genetic circuits become more sophisticated, the size and complexity of data about their designs increases. This data captured goes beyond monolithic genetic sequences and towards circuit modularity and functional details, which are beneficial for analyzing circuit performance and establishing design automation techniques. However, the accessibility, visualisation and usability of design data (and metadata) have received relatively little attention to date. Here, we present a method to turn circuit designs into networks and showcase its potential to enhance the utility of design data. Since networks are dynamic structures, initial graphs can be interactively shaped into sub-networks of relevant information based on requirements such as abstraction, hierarchy and protein interactions. Additionally, several visual changes can be applied, such as colouring or clustering nodes based on types (e.g., genes or promoters), resulting in easier comprehension from a user perspective. This approach allows circuit designs to be coupled to other networks, such as metabolic pathways or implementation protocols captured in graph-like formats. Therefore, we advocate using networks to structure, access and improve synthetic biology information. Widespread inhibitory projections from the cerebellar interposed nucleus The cerebellum consists of parallel parasagittal modules that contribute to diverse behaviors, spanning motor to cognitive. Recent work employing cell-type specific tracing has identified circumscribed output channels of the cerebellar nuclei that could confer tight functional specificity. These studies have largely focused on excitatory projections of the cerebellar nuclei, however, leaving open the question of whether inhibitory neurons also constitute multiple output modules. We mapped output and input patterns to intersectionally restricted cell types of the interposed and adjacent interstitial nuclei. In contrast to the widespread assumption of primarily excitatory outputs and restricted inferior olive-targeting inhibitory output, we found that inhibitory neurons from this region ramified widely within the brainstem, targeting both motor- and sensory-related nuclei, distinct from excitatory output targets. Despite differences in output targeting, monosynaptic rabies tracing revealed largely shared afferents to both cell classes. We discuss the potential novel functional roles for inhibitory outputs in the context of cerebellar theory. Microvesicles Transfer Mitochondria and Increase Mitochondrial Function in Brain Endothelial Cells We have demonstrated, for the first time that microvesicles, a sub-type of extracellular vesicles (EVs) derived from hCMEC/D3: a human brain endothelial cell (BEC) line transfer polarized mitochondria to recipient BECs in culture and to neurons in mice acute brain cortical and hippocampal slices. This mitochondrial transfer increased ATP levels by 100 to 200-fold (relative to untreated cells) in the recipient BECs exposed to oxygen-glucose deprivation, an in vitro model of cerebral ischemia. We have also demonstrated that transfer of microvesicles, the larger EV fraction, but not exosomes resulted in increased mitochondrial function in hypoxic endothelial cultures. Gene ontology and pathway enrichment analysis of EVs revealed a very high association to glycolysis-related processes. In comparison to heterotypic macrophage- derived EVs, BEC-derived EVs demonstrated a greater selectivity to transfer mitochondria and increase endothelial cell survival under ischemic conditions.

HighlightsO_LIMicrovesicles transfer mitochondria to endothelial cells and brain slice neurons
C_LIO_LIMitochondrial transfer increased ATP in ischemic brain endothelial cells (BECs)
C_LIO_LITransfer of microvesicles increased mitochondrial function in hypoxic BECs
C_LIO_LITransfer of exosomes did not affect mitochondrial function in hypoxic BECs
C_LIO_LIHomotypic BEC-derived EVs result in greater ATP levels in the recipient BECs
C_LI Groping in the fog: soaring migrants exhibit wider scatter in flight directions and respond differently to wind under low visibility conditions Atmospheric conditions are known to affect flight propensity, behaviour during flight, and migration route in birds. Yet, the effects of fog have only rarely been studied although they could disrupt orientation and hamper flight. Fog could limit the visibility of migrating birds such that they might not be able to detect landmarks that guide them during their journey. Soaring migrants modulate their flight speed and direction in relation to the wind vector to optimize the cost of transport. Consequently, landmark-based orientation, as well as adjustments of flight speed and direction in relation to wind conditions, could be jeopardized when flying in fog. Using a radar system operated in a migration bottleneck (Strait of Messina, Italy), we studied the behaviour of soaring birds under variable wind and fog conditions over two consecutive springs (2016 and 2017), discovering that migrating birds exhibited a wider scatter of flight directions and responded differently to wind under fog conditions. Birds flying through fog deviated more from the mean migration direction and increased their speed with increasing crosswinds. In addition, airspeed and groundspeed increased in the direction of the crosswind, causing the individuals to drift laterally. Our findings represent the first quantitative empirical evidence of flight behaviour changes when birds migrate through fog and explain why low visibility conditions could risk their migration journey. 'rbioacc': an R-package to analyse toxicokinetic data O_LI rbioacc is an  R-package dedicated to the analysis of experimental data collected from bioaccumulation tests during which organisms are exposed to a chemical (exposure phase) and then put into a clean media (depuration phase). Internal concentrations are regularly measured over time all along the experiment.
C_LIO_LI rbioacc provides ready-to-use functions to visualize and fully analyze such data. Under a Bayesian framework, this package fits a generic one-compartment toxicokinetic (TK) model automatically built from the data. It provides TK parameter estimates (appropriate uptake and elimination rates) and bioaccumulation metrics (e.g., BCF, BSAF, BMF). All parameter estimates, bioaccumulation metrics as well as predictions of internal concentrations into organisms are delivered with their uncertainty.
C_LIO_LIBioaccumulation metrics are primarily provided in support of environmental risk assessment, in full compliance with regulatory requirements required to approve marketing applications of chemical substances.
C_LIO_LIThis paper gives brief worked examples of the use of  rbioacc from data collected through standard bioaccumulation tests, and publicly available within the scientific literature. These examples constitute step-by-step user-guides to analyze any new data set, uploaded in the right format.
C_LI A General Biphasic-Bodyweight Model for Scaling Basal Metabolic Rate, Glomerular Filtration Rate and Drug Clearance from Birth to Adulthood Understanding the maturation process of human physiology and metabolism has broad medical and pharmaceutical implications. Age and bodyweight are frequently considered as separate variables in modeling the dynamical changes of human organ functions and of drug clearance from birth to adulthood. The objective of this study is to propose a unified, continuous and bodyweight-only equation to quantify the changes of human basal metabolic rate (BMR), glomerular filtration rate (GFR) and drug clearance (CL) from infancy to adulthood. The BMR datasets were retrieved from a comprehensive historical database of male and female subjects (0.02 to 64 years). The CL datasets for 17 drugs and the GFR dataset were generated from age-incorporated maturation-and-growth models with reported parameter values. The model used in the simulation is independent of the proposed model. A statistical approach was used to simulate the model generated CL and GFR data for a hypothetical population with 26 age groups (ranging from 0 to 20 years). Besides, individual CL data for one drug, and sparse PBPK-modeled CL values for two drugs were also included for analysis. A 4-parameter, mixed-allometry equation with two power-law functions of bodyweight was proposed and evaluated as a general model using nonlinear regression and dimensionless analysis. All datasets universally reveal biphasic curves with two distinct linear segments on log-log plots. Compared with simple allometry, the biphasic model fits satisfactorily to all datasets (based on Akaikes Information Criterion and residual plots). The biphasic equation consists of two reciprocal allometric terms that asymptotically determine the overall curvature. The fitting results show a superlinear scaling phase (slope >1; ca. 1.5 - 3.5) below the characteristic bodyweight at the phase transition; and above which, a sublinear scaling phase (slope <1; ca. 0.5 - 0.7) is evident. The phase-transition bodyweight is ranging from 5 to 20 kg (corresponding to 0.5 - 9 years) and the mean value is around 10 kg ([~]2 years) for all data sets. The dimensionless analysis generalizes, and offers quantitative realization of, the maturation and growth process. In conclusion, the proposed mixed-allometry equation is a generic model that quantitatively describes the phase transition occurring in the human maturation process of BMR, GFR and drug CL. Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. Glypiated FGF directs contact-mediated bidirectional signaling to self-regulate tissue-specific dispersion How signaling proteins generate a multitude of information to organize tissue patterns is critical to understanding morphogenesis. In Drosophila, FGF produced in wing-disc cells regulates the development of the disc-associated air-sac-primordium (ASP). Here, we show that FGF is Glycosylphosphatidylinositol-anchored to the producing cell surface and that this modification both inhibits free FGF secretion and activates target-specific bidirectional FGF-FGFR signaling through cytonemes. FGF-source and recipient ASP cells extend cytonemes that present FGF and FGFR on their surfaces and reciprocally recognize each other over distance by contacting through CAM-like FGF-FGFR binding. Contact-mediated FGF-FGFR binding induces bidirectional signaling, which, in turn, promotes ASP and source cells to polarize cytonemes toward each other and reinforce signaling contacts. Subsequent un-anchoring of FGFR-bound-FGF from the source cell membrane dissociates cytoneme contacts and delivers FGF target-specifically to ASP cytonemes for paracrine functions. Thus, GPI-anchored FGF organizes both source and recipient cells and self-regulates its cytoneme-mediated tissue-specific dispersion and signaling. Standardized genome-wide function prediction enables comparative functional genomics: a new application area for Gene Ontologies in plants BackgroundGenome-wide gene function annotations are useful for hypothesis generation and for prioritizing candidate genes potentially responsible for phenotypes of interest. We functionally annotated the genes of 18 crop plant genomes across 14 species using the GOMAP pipeline.

ResultsBy comparison to existing GO annotation datasets, GOMAP-generated datasets cover more genes, contain more GO terms, and produce datasets similar in quality (based on precision and recall metrics using existing gold standards as the basis for comparison). From there, we sought to determine whether the datasets across multiple species could be used together to carry out comparative functional genomics analyses in plants. To test the idea and as a proof of concept, we created dendrograms of functional relatedness based on terms assigned for all 18 genomes. These dendrograms were compared to well-established species-level evolutionary phylogenies to determine whether trees derived were in agreement with known evolutionary relationships, which they largely are. Where discrepancies were observed, we determined branch support based on jack-knifing then removed individual annotation sets by genome to identify the annotation sets causing unexpected relationships.

ConclusionsGOMAP-derived functional annotations used together across multiple species generally retain sufficient biological signal to recover known phylogenetic relationships based on genome-wide functional similarities, indicating that comparative functional genomics across species based on GO data hold promise for generating novel hypotheses about comparative gene function and traits. Dilute phase oligomerization can oppose phase separation and modulate material properties of a ribonucleoprotein condensate Ribonucleoprotein bodies are exemplars of membraneless biomolecular condensates that can form via spontaneous or driven phase transitions. The fungal protein Whi3 forms ribonucleoprotein condensates with different RNA molecules, and these condensates are implicated in key processes such as cell-cycle control and generating cell polarity. Whi3 has a modular architecture that includes a Q-rich intrinsically disordered region (IDR) and a tandem RNA recognition module. Here, we demonstrate that a 21-residue stretch within the Q-rich IDR has a weak intrinsic preference for forming alpha-helical conformations. Through mutagenesis, we find that increased alpha helicity enhances oligomerization in the dilute phase. One consequence of enhanced oligomerization is a dilution of Whi3 in the dense phase. The opposite behavior is observed when helicity within the 21-residue stretch of the Q-rich region is abrogated. Thus, the formation of dilute phase oligomers, driven by a specific sequence motif and potential synergies with the rest of the IDR, opposes incorporation of the Whi3 protein into the dense phase, thereby altering the dense phase stoichiometry of protein to RNA. Our findings, which stand in contrast to other systems where oligomerization has been shown to enhance the drive for phase separation, point to a novel mechanism that might be operative for influencing compositions of condensates. Our work also points to routes for designing synthetic ribonucleoprotein condensates whereby modulation of protein oligomerization via homotypic interactions can impact dense phase concentrations, stoichiometries, and material properties.

SignificanceA large sub-class of biomolecular condensates are linked to RNA regulation and are known as ribonucleoprotein (RNP) bodies. While extensive work has identified driving forces for biomolecular condensate formation, relatively little is known about forces that oppose assembly. Here, using a fungal RNP protein, Whi3, we show that a portion of its intrinsically disordered, glutamine-rich region modulates phase separation by forming transient alpha helical structures that promote the assembly of dilute phase oligomers. These oligomers detour Whi3 proteins from condensates, thereby impacting the driving forces for phase separation, the protein-to-RNA ratio in condensates, and the material properties of condensates. Our findings show how nanoscale conformational and oligomerization equilibria can influence mesoscale phase equilibria. Fission yeast polycystin Pkd2p promotes the cell expansion and antagonizes the Hippo pathway SIN Polycystins are conserved mechanosensitive channels whose mutations lead to the common human renal disorder ADPKD. Previously we discovered that the plasma membrane-localized fission yeast homologue Pkd2p is an essential protein required for cytokinesis, but the mechanism remains unclear. Here, we isolated a novel temperature-sensitive mutant pkd2-B42. Among its strong growth defects, the most unique was that many mutant cells often lost significant portion of their volume in just 5 minutes followed by a gradual recovery, a process that we termed Deflation. Unlike cell lysis, deflation did not result in the plasma membrane rupture and it occurred independently from the cell cycle progression. The tip extension of pkd2-B42 cells was 80% slower than the wild type and their turgor pressure was 50% lower. Both pkd2-B42 and the other mutant pkd2-81KD partially rescued the mutants of the yeast Hippo signaling pathway Septation Initiation Network, by preventing cell lysis, enhancing septum formation, and doubling the number of Sid2/Mob1 molecules at the spindle pole bodies. We conclude that Pkd2p promotes cell size expansion during interphase by regulating turgor pressure and antagonizes SIN during cytokinesis.

Summary statementMutations of polycystins lead to human genetic disorder ADPKD. We discovered that the fission yeast homologue Pkd2p promotes the cell expansion during interphase growth and antagonizes the Hippo pathway SIN during cytokinesis. 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 combined human in silico and CRISPR/Cas9-mediated in vivo zebrafish based approach for supporting gene target validation in early drug discovery The clinical heterogeneity of heart failure has challenged our understanding of the underlying genetic mechanisms of this disease. In this respect, large-scale patient DNA sequencing studies have become an invaluable strategy for identifying potential genetic contributing factors. The complex aetiology of heart failure, however, also means that in vivo models are vital to understand the links between genetic perturbations and functional impacts. Traditional approaches (e.g. genetically-modified mice) are optimal for assessing small numbers of proposed target genes, but less practical when multiple targets are identified. The zebrafish, in contrast, offers great potential for higher throughput in vivo gene functional assessment to aid target prioritisation and support definitive studies undertaken in mice. Here we used whole-exome sequencing and bioinformatics on human patient data to identify 3 genes (API5, HSPB7, and LMO2) suggestively associated with heart failure that were also predicted to play a broader role in disease aetiology. The role of these genes in cardiovascular system development and function was then further investigated using in vivo CRISPR/Cas9-mediated gene mutation analysis in zebrafish. We observed multiple impacts in F0 knockout zebrafish embryos (crispants) following effective somatic mutation, including reductions in ventricle size, pericardial oedema, and chamber malformation. In the case of lmo2, there was also a significant impact on cardiovascular function as well as an expected reduction in erythropoiesis. The data generated from both the human in silico and zebrafish in vivo assessments undertaken supports roles for API5, HSPB7, and LMO2 in human cardiovascular disease and identifies them as potential drug targets for further investigation. The data presented also supports the use of human in silico genetic variant analysis, in combination with zebrafish crispant phenotyping, as a powerful approach for assessing gene function as part of an integrated multi-level drug target validation strategy. Antibody ligation of carcinoembryonic antigen-related cell adhesion molecule 1 (CEACAM1), CEACAM3, and CEACAM6, differentially enhance the cytokine release of human neutrophils in responses to Candida albicans Invasive candidiasis, often caused by Candida albicans, is an important healthcare-associated fungal infection that results in a high mortality rate of up to 40%. Neutrophils are the first line of defense during Candida infections. They can initiate various killing mechanisms and release cytokines to attract further immune cells to the site of infection. These responses are tightly controlled, since they can also lead to severe tissue/organ damage. We hypothesized that the regulation of C. albicans-specific neutrophil functions by the immunoregulatory C. albicans receptors CEACAM1, CEACAM3, and CEACAM6 are involved in the immune pathology of candidemia. Here, we analyzed the effects of the specific antibodies B3-17, 308/3-3, and 1H7-4B, respectively, targeting the three CEACAM receptors on C. albicans-induced neutrophil responses. We show that ligation of CEACAM6 significantly enhanced the response to C. albicans, as evidenced by the increased CXCL8/IL-8 secretion. By assessing the transcriptional responses, we found that CEACAM6 ligation and to some extent CEACAM1 ligation, but not CEACAM3 ligation, resulted in altered gene regulation of the C. albicans-stimulated neutrophils. Genes that were differentially regulated by the different CEACAM-targeting antibodies were analyzed for affected cellular processes and signaling pathways using various bioinformatics methods, including integrated network analyses and dynamic simulations of signaling cascades. Predicted changes in cellular pathways and cellular functions included CEACAM-specific alterations in apoptosis and cytokine secretion. In particular, we verified predicted changes in IL-1{beta}/IL-6 expression in response to the antibody ligation of all three targeted CEACAMs and apoptosis induction by anti-CEACAM6 antibody treatment in presence of C. albicans stimulation. Specifically, CEACAM6 ligation by 1H7-4B enhanced neutrophil apoptosis and increased immediate and long-term cytokine release in responses to C. albicans. CEACAM1 ligation by B3-17 mainly enhanced the immediate secretion of pro-inflammatory cytokines, and CEACAM3 ligation by 308/3-3 increased the long-term release of pro-inflammatory cytokines. Taken together, we demonstrated for the first time that CEACAM receptors have an important and differential impact on the regulation of C. albicans-induced immune functions in human neutrophils. Getting more out of FLAG-Tag co-immunoprecipitation mass spectrometry experiments using FAIMS. Co-immunoprecipitation of proteins coupled to mass spectrometry is critical for the understanding of protein interaction networks. In instances where a suitable antibody is not available, it is common to graft synthetic tags onto target protein sequences and allowing the use of commercially available antibodies for affinity purification. A common approach is through FLAG-Tag co-immunoprecipitation. To allow the selective elution of protein complexes, competitive displacement using a large molar excess of the tag peptides is often carried out. Yet, this creates downstream challenges for the mass spectrometry analysis due to the presence of large quantities of these peptides. Here, we demonstrate that Field Asymmetric Ion Mobility Spectrometry (FAIMS), a gas phase ion separation device prior to mass spectrometry analysis can be applied to FLAG-Tag co-immunoprecipitation experiment to increase the depth of protein coverage. By excluding these abundant tag peptides, we were able to observe deeper coverage of interacting proteins and as a result, deeper biological insights, without the need for additional sample handling or altering sample preparation protocols. Water levels primarily drive variation in photosynthesis and nutrient use of scrub Red Mangroves in the southeastern Florida Everglades We investigated how mangrove-island micro-elevation (i.e., habitat: center vs. edge) affects tree physiology in a scrub mangrove forest of the southeastern Everglades. We measured leaf gas exchange rates of scrub Rhizophora mangle L. trees monthly during 2019, hypothesizing that CO2 assimilation (Anet) and stomatal conductance (gsw) would decline with increasing water levels and salinity, expecting more-considerable differences at mangrove-island edges than centers, where physiological stress is greatest. Water levels varied between 0 and 60 cm from the soil surface, rising during the wet season (May-October) relative to the dry season (November-April). Porewater salinity ranged from 15 to 30 ppt, being higher at mangrove-island edges than centers. Anet maximized at 15.1 {micro}mol m-2 s-1, and gsw was typically <0.2 mol m-2 s-1, both of which were greater in the dry than the wet season and greater at island centers than edges, with seasonal variability being roughly equal to variation between habitats. After accounting for season and habitat, water level positively affected Anet in both seasons but did not affect gsw. Our findings suggest that inundation stress (i.e., water level) is the primary driver of variation in leaf gas exchange rates of scrub mangroves in the Florida Everglades, while also constraining Anet more than gsw. The interaction between inundation stress due to permanent flooding and habitat varies with season as physiological stress is alleviated at higher-elevation mangrove-island center habitats during the dry season. Freshwater inflows during the wet season, increase water levels and inundation stress at higher-elevation mangrove-island centers, but also potentially alleviate salt and sulfide stress in soils. Thus, habitat heterogeneity leads to differences in nutrient and water acquisition and use between trees growing in island centers versus edges, creating distinct physiological controls on photosynthesis, which likely affect carbon flux dynamics of scrub mangroves in the Everglades. Post-translational flavinylation is associated with diverse extracytosolic redox functionalities throughout bacterial life Disparate redox activities that take place beyond the bounds of the prokaryotic cell cytosol must connect to membrane or cytosolic electron pools. Proteins post-translationally flavinylated by the enzyme ApbE mediate electron transfer in several characterized extracytosolic redox systems but the breadth of functions of this modification remains unknown. Here we present a comprehensive bioinformatic analysis of 31,910 prokaryotic genomes that provides evidence of extracytosolic ApbEs within ~50% of bacteria and the involvement of flavinylation in numerous uncharacterized biochemical processes. By mining flavinylation-associated gene clusters, we identify five protein classes responsible for transmembrane electron transfer and two domains of unknown function (DUF2271 and DUF3570) that are flavinylated by ApbE. We observe flavinylation/iron transporter gene colocalization patterns that implicate functions in iron reduction and assimilation. We find associations with characterized and uncharacterized respiratory oxidoreductases that highlight roles of flavinylation in respiratory electron transport chains. Finally, we identify interspecies gene cluster variability consistent with flavinylation/cytochrome functional redundancies and discover a class of "multi-flavinylated proteins'' that may resemble multiheme cytochromes in facilitating longer distance electron transfer. These findings provide key mechanistic insight into an important facet of bacterial physiology and establish flavinylation as a functionally diverse mediator of extracytosolic electron transfer.

Impact StatementBioinformatic and biochemical studies find that covalently bound flavins are common and participate in wide-ranging extracytosolic redox activities throughout bacterial life.