GPI-anchored FGF directs cytoneme-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.
A cell cycle-dependent GARP-like transcriptional repressor regulates the initiation of differentiation in Giardia lamblia Transcriptional regulation of differentiation is critical for parasitic pathogens to adapt to environmental changes and regulate transmission. In response to encystation stimuli, Giardia lamblia cells shift from G1+S phase of the cell cycle to G2+M. By 2-4 hours, cyst wall proteins are upregulated, indicating that key regulatory steps occur within the first four hours of encystation. However, all characterized transcription factors (TFs) in Giardia have only been investigated at later time points of encystation. How TFs initiate encystation and link it to the cell cycle remains enigmatic. Here, we systematically screened six putative early upregulated TFs for nuclear localization, established their dynamic expression profiles, and determined their functional role in regulating encystation. We found a critical repressor, GLP4 that increases rapidly after 30 min of encystation stimuli and downregulates encystation specific markers, including Cyst Wall Proteins (CWPs) and enzymes in the cyst N-acetyl galactosamine (GalNAc) pathway. Depletion of GLP4 increases cyst production. Importantly, we observe that G2+M cells exhibit higher levels of CWP1 resulting from the activation of MYB2, a TF previously linked to encystation in Giardia. GLP4 upregulation occurs in G1+S cells, suggesting a role in repressing MYB2 and encystation specific genes in the G1+S phase of the cell cycle. Furthermore, we demonstrate that depletion of GLP4 upregulates MYB2 and promotes encystation while overexpression of GLP4 downregulates MYB2 and represses encystation. Together, these results suggest that Giardia employs a dose-dependent transcriptional response that involves the cell cycle regulated repressor GLP4 to orchestrate MYB2 and entry into the encystation pathway.

ImportanceTransition between life cycle stages is a common feature among parasitic pathogens and its regulation must be optimized to balance persistence of infection with transmission. The early transcription factors (TFs) regulating commitment to differentiate are totally unknown in Giardia. In this work, we identified GLP4, a previously uncharacterized GARP-like TF, as an early-acting transcriptional repressor that inhibits G1+S cells from entering the encystation pathway. GLP4 is therefore a key regulator controlling the balance between proliferative growth and terminal differentiation into infective cysts.
Tobramycin-induced secretion of P. aeruginosa 5' tRNA-fMet halves suppresses lung inflammation via AGO2 gene silencing Although inhaled tobramycin increases lung function in people with cystic fibrosis (pwCF), the density of P. aeruginosa in the lungs is only modestly reduced by tobramycin; hence, the mechanism whereby tobramycin improves lung function is unclear. Here, we demonstrate that tobramycin increases the abundance of two 5' tRNA-fMet halves in outer membrane vesicles (OMVs) secreted by P. aeruginosa and that the 5' tRNA-fMet halves reduce IL-8 secretion by CF bronchial epithelial cells (CF-HBECs). In mouse lung, the 5' tRNA-fMet halves attenuate KC secretion and neutrophil recruitment. We also report that the 5' tRNA-fMet halves suppress pro-inflammatory network gene expression by an Argonaut 2 (AGO2)-mediated gene silencing mechanism, thereby reducing IL-8 secretion in CF-HBECs. Moreover, tobramycin reduces the IL-8 concentration and neutrophil content in bronchoalveolar lavage fluid of pwCF. Thus, we conclude that tobramycin improves lung function in part by reducing chronic inflammation and neutrophil-mediated lung damage in pwCF. Structural Analysis of Spike Protein Mutations in the SARS-CoV-2 P.3 Variant A SARS-CoV-2 lineage designated as P.3 with multiple signature mutations in the Spike protein region was recently reported with cases from the Central Visayas Region of the Philippines. Whole genome sequencing revealed that the 33 samples under this lineage all contain the E484K, N501Y, and P681H Spike mutations previously found in variants of concern (VOC) such as the B.1.351, the P.1 and B.1.1.7 variants first reported in South Africa, Brazil, and the United Kingdom, respectively. The possible implications of the mutations found in the Spike protein of P.3 were analyzed for their potential effects on structure, stability, and molecular surface character. The analysis suggests that these mutations could significantly impact the possible interactions of the Spike protein with the ACE2 receptor and neutralizing antibodies, and warrants further clinical investigation. Some of the mutations affecting the N and C terminal domains may have effects on Spike monomer and trimer stability. This report provides insights on relevant targets for the design of future diagnostics, therapeutics and vaccines against the evolving SARS-CoV-2 variants in the Philippines. Creatine transporter deficient rat model show motor dysfunction, cerebellar alterations and muscle creatine deficiency without muscle atrophy Creatine (Cr) is a nitrogenous organic acid and plays roles as fast phosphate energy buffer to replenish ATP, osmolyte, antioxidant, neuromodulator, and as a compound with anabolic and ergogenic properties in muscle. Cr is taken from the diet or endogenously synthetized by the enzymes AGAT and GAMT, and specifically taken up by the transporter SLC6A8. Loss-of-function mutations in the genes encoding for the enzymes or the transporter cause Cerebral Creatine Deficiency Syndromes (CCDS). CCDS are characterized by brain Cr deficiency, intellectual disability with severe speech delay, behavioral troubles, epilepsy and motor dysfunction. Among CCDS, the X-linked Cr transporter deficiency (CTD) is the most prevalent with no efficient treatment so far. Different animal models of CTD show reduced brain Cr levels, cognitive deficiencies and together they cover other traits similar to those of patients. However, motor function was poorly explored in CTD models and some controversies in the phenotype exist in comparison with CTD patients. Our recently described Slc6a8Y389C knock-in (KI) rat model of CTD showed mild impaired motor function, morphological alterations in cerebellum, reduced muscular mass, Cr deficiency and increased guanidinoacetate content in muscle, although no consistent signs of muscle atrophy. Our results indicate that such motor dysfunction co-occurred with both nervous and muscle dysfunction, suggesting that muscle strength and performance as well as neuronal connectivity might be affected by this Cr deficiency in muscle and brain. Network Analysis Reveals Different Cellulose Degradation Strategies across Trichoderma harzianum Strains Associated with XYR1 and CRE1 Trichoderma harzianum, whose gene expression is tightly controlled by the transcription factors (TFs) XYR1 and CRE1, is a potential candidate for hydrolytic enzyme production. Here, we performed a network analysis of T. harzianum IOC-3844 and T. harzianum CBMAI-0179 to explore how the regulation of these TFs varies between these strains. In addition, we explored the evolutionary relationships of XYR1 and CRE1 protein sequences among Trichoderma spp. The results of the T. harzianum strains were compared with those of Trichoderma atroviride CBMAI-0020, a mycoparasitic species. Although transcripts encoding carbohydrate-active enzymes (CAZymes), TFs, transporters, and proteins with unknown functions were coexpressed with cre1 or xyr1, other proteins indirectly related to cellulose degradation were identified. The enriched GO terms describing the transcripts of these groups differed across all strains, and several metabolic pathways with high similarity between both regulators but strain-specific differences were identified. In addition, the CRE1 and XYR1 subnetworks presented different topology profiles in each strain, likely indicating differences in the influences of these regulators according to the fungi. The hubs of the cre1 and xyr1 groups included transcripts not yet characterized or described as being related to cellulose degradation. The first-neighbor analyses confirmed the results of the profile of the coexpressed transcripts in cre1 and xyr1. The analyses of the shortest paths revealed that CAZymes upregulated under cellulose degradation conditions are most closely related to both regulators, and new targets between such signaling pathways were discovered. Although the evaluated T. harzianum strains are phylogenetically close and their amino acid sequences related to XYR1 and CRE1 are very similar, the set of transcripts related to xyr1 and cre1 differed, suggesting that each T. harzianum strain used a specific regulation strategy for cellulose degradation. More interestingly, our findings may suggest that XYR1 and CRE1 indirectly regulate genes encoding proteins related to cellulose degradation in the evaluated T. harzianum strains. An improved understanding of the basic biology of fungi during the cellulose degradation process can contribute to the use of their enzymes in several biotechnological applications and pave the way for further studies on the differences across strains of the same species. 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. Anchorage by seed mucilage prevents seed dislodgement in high surface flow: a mechanistic investigation O_LIBackground and Aims: Seed mucilage is a common and highly diverse trait shared among thousands of angiosperm species. While long recognized that mucilage allows seeds to anchor to substrates (antitelechory), resisting abiotic and biotic dislodgement, we still lack a mechanistic understanding of this process.
C_LIO_LIMethods: We propose a mechanistic model of how mucilage affects substrate anchorage and fluid resistance, ultimately contributing to dislodgement resistance. To test this model, we subjected mucilaginous seeds of 52 species, varying in eight measured seed traits, to seven days of continuous water flow at a range of dislodgement potentials.
C_LIO_LIKey Results: Supporting our model, mucilage mass increased force necessary to dislodge both dry and wet seeds; our measurement of the dislodgement force of dry mucilage explained time to dislodgement well. The effect size was remarkably large; increasing the standardized mucilage mass by one standard deviation resulted in a 280-fold increase in the time to dislodgement. Fluid resistance was largely dependent on speed of water flow and the seeds modified drag coefficient, but not seed traits. Neither mucilage expansion speed nor mucilage decay rate explained dislodgement potential well.
C_LIO_LIConclusions: Our results suggest that the degree of anchorage to substrate, measured with a simple dislodgement force assay, is highly predictive of mucilaginous seed retention in highly erosive environments. In contrast, we found that other seed and mucilage traits are of lesser importance to anchorage.
C_LI 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. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. Dynamic gain decomposition reveals functional effects of dendrites, ion channels and input statistics in population coding The information processing capabilities of large neuronal circuits form the basis of higher-level cognition. It is unknown how the constituent cells properties interact with the external stimuli to shape the networks collective activity patterns and the associated computations. The dynamic gain function, a spectrally resolved approximation of a populations coding capability, is sensitive to the properties of stimuli and cells and could provide a mechanistic link from cell-level to network-level performance. However, the dynamic gains shape is co-determined in a complex manner by axonal and somatodendritic parameters and the populations operating regime, which so far precluded an understanding of any individual parameters impact. Here, we decomposed the dynamic gain function into three components corresponding to separate signal transformations. This allowed attribution of network-level encoding features to specific cell-level parameters. Applying the method to data from real neurons and biophysically plausible models, we found: 1. The encoding bandwidth of real neurons, approximately 400 Hz, is constrained by the voltage dependence of axonal currents during an early stage of action potential initiation. 2. The state-of-the-art models only achieve encoding bandwidths around 100 Hz and are mostly limited by subthreshold processes instead. 3. Large dendrites and low-threshold potassium currents modulate the bandwidth by shaping the subthreshold stimulus-to-voltage transformation. Dynamic gain decomposition can physiologically interpret changes in the dynamic gain curve, which is particularly desirable in studies of cellular effects in spectrinopathies and neurodegeneration. By pinpointing shortcomings of current models, it also guides inference of neuron models best suited for large-scale network simulations.

ImpactThe dynamic gain function quantifies how neurons can engage in collective, network-level activity, shaping brain rhythms and corresponding information encoding capabilities. Its shape results from a complex interaction between properties of different molecules (ion channels) and cell compartments (morphology, resistance), and is so far only understood for the simplest neuron models. Here we provide an interpretable analysis, decomposing dynamic gain based on the stimulus transformation steps in individual neurons. It is applicable to models of any complexity and real neurons, and allows attribution of dynamic gain changes to specific sub- and suprathreshold processes. Using this decomposition method, we reveal the relevance of subthreshold potassium channels for ultrafast information encoding and expose the shortcomings of even the state-of-the-art neuron models. Hematopoietic differentiation is characterized by a transient peak of entropy at a single cell level. Hematopoietic differentiation has been metaphorically represented as linear trajectories with discrete steps from hematopoietic stem cells to mature cells. While the transcriptional state of cells at the beginning or at the end of these trajectories are well described from bulk analysis, what happens in the intermediate states has remained elusive until the use of single cell approaches. Applying Shannon entropy to measure cell-to-cell variability among cells at the same stage of differentiation, we observed a transient peak of gene expression variability in all the hematopoietic differentiation pathways. Strikingly, genes with the highest entropy variation in a given differentiation pathway matched genes known as pathway-specific, whereas genes with the highest expression variation were common to all pathways. Finally, we showed that the level of cell-to-cell variation is increased in the most immature compartment of hematopoiesis in myelodysplastic syndromes. These data suggest that differentiation could be better conceptualized as a dynamical stochastic process with a transient stage of cellular indetermination. Effects of phenotypic variation on consumer coexistence and prey community structure A popular idea in ecology is that trait variation among individuals from the same species may promote the coexistence of competing species. However, theoretical and empirical tests of this idea have yielded inconsistent findings. We manipulated intraspecific trait diversity in a ciliate competing with a nematode for bacterial prey in experimental microcosms. We found that intraspecific trait variation inverted the original competitive hierarchy to favor the consumer with variable traits, ultimately resulting in competitive exclusion. This competitive outcome was driven by foraging traits (size, speed, and directionality) that increased the ciliates fitness ratio and niche overlap with the nematode. The interplay between consumer trait variation and competition resulted in non-additive cascading effects - mediated through prey defense traits - on prey community assembly. Our results suggest that predicting consumer competitive population dynamics and the assembly of prey communities will require understanding the complexities of trait variation within consumer species. Depletion of Sun1/2 induces heterochromatin accrual in mesenchymal stem cells during adipogenesis Critical to the mechano-regulation of mesenchymal stem cells (MSC), Linker of the Nucleoskeleton and Cytoskeleton (LINC) complex transduces cytoskeletal forces to the nuclei. The LINC complex contains outer nuclear membrane Nesprin proteins that associate with the cytoskeleton and their inner nuclear membrane couplers, Sun proteins. In addition to coupling Nesprin-associated cytoskeletal elements to inner nuclear membrane, Sun proteins also function in regulating nuclear mechanics and chromatin tethering to inner nuclear membrane. This suggests that release of LINC-mediated cytoskeletal connections from cell nuclei may have different effects on chromatin organization and MSC differentiation than those due to ablation of intra-nuclear Sun proteins. To test this hypothesis we compared Sun1/2 depletion with expression of the dominant-negative KASH domain (dnKASH) that inhibits Nesprin-Sun association. In cells cultured under adipogenic conditions, disconnecting LINC from the cytoskeleton through dnKASH expression, increased adipogenic gene expression and fat droplet formation; heterochromatin H3K27 and H3K9 methylation was unaltered. In contrast, Sun1/2 depletion inhibited adipogenic gene expression and fat droplet formation; as well the anti-adipogenic effect of Sun1/2 depletion was accompanied by increased H3K9me3, which was enriched on Adipoq, silencing this fat locus. We conclude that releasing the nucleus from cytoskeletal constraints via dnKASH accelerates adipogenesis while depletion of Sun1/2 increases heterochromatin accrual on adipogenic genes in a fashion independent of LINC complex function. Therefore, while these two approaches both disable LINC functions, their divergent effects on the epigenetic landscape indicate they cannot be used interchangeably to study mechanical regulation of cell differentiation. Interpreting ciliopathy-associated missense variants of uncertain significance (VUS) in an animal model PurposeBetter methods are required to interpret the pathogenicity of disease-associated variants of uncertain significance (VUS), which cannot be actioned clinically. In this study, we explore the use of a tractable animal model (Caenorhabditis elegans) for in vivo interpretation of missense VUS alleles of TMEM67, a cilia gene associated with ciliopathies.

MethodsCRISPR/Cas9 gene editing was used to generate homozygous knock-in C. elegans worm strains carrying TMEM67 patient variants engineered into the orthologous gene (mks-3). Quantitative phenotypic assays of sensory cilia structure and function measured if the variants affect mks-3 gene function. Results from worms were validated by a genetic complementation assay in a human TMEM67 knock-out hTERT-RPE1 cell line that tests a TMEM67 signaling function.

ResultsAssays in C. elegans accurately distinguished between known benign (Asp359Glu, Thr360Ala) and known pathogenic (Glu361Ter, Gln376Pro) variants. Analysis of eight missense VUS generated evidence that three are benign (Cys173Arg, Thr176Ile, Gly979Arg) and five are pathogenic (Cys170Tyr, His782Arg, Gly786Glu, His790Arg, Ser961Tyr).

ConclusionEfficient genome editing and quantitative functional assays in C. elegans make it a tractable in vivo animal model that allows rapid, cost-effective interpretation of ciliopathy-associated missense VUS alleles. Darwin's small and medium ground finches might have taste preferences, but not for human foods Urbanization is rapidly changing ecological niches. On the inhabited Galapagos Islands, Darwins finches consume human-introduced foods preferentially; however, it remains unclear why. Here we presented pastry with flavour profiles typical of human foods (oily, salty, sweet) to small ground finches (Geospiza fuliginosa) and medium ground finches (Geospiza fortis) to test if latent taste preferences might drive selection of human foods. If human-food flavours were consumed more than a neutral or bitter control only at sites with human foods, then we predicted tastes were acquired after urbanization; however, if no site-differences were found then this would indicate latent taste preferences. Contrary to both predictions, we found little evidence that human-food flavours were preferred compared to control flavours at any site. Instead, finches showed a weak aversion to oily foods, but only at remote (no human foods present) sites. This was further supported by behavioural responses, with beak wiping occurring more often at remote sites after finches tasted flavours associated with human foods. Our results suggest, therefore, that while Darwins finches may have acquired a tolerance to the flavours of human food, latent taste preferences are unlikely to have played a major role in their dietary response to increased urbanization. Behavioral and postural analyses establish sleep-like states for mosquitoes that can impact host landing and blood feeding Sleep is an evolutionarily conserved process that has been described in different animal systems. For insects, sleep characterization has been primarily achieved using behavioral and electrophysiological correlates in a few systems. Sleep in mosquitoes, which are important vectors of disease-causing pathogens, has not been directly examined. This is surprising as circadian rhythms, which have been well studied in mosquitoes, influence sleep in other systems. In this study, we characterized sleep in mosquitoes using body posture analysis and behavioral correlates, and quantified the effect of sleep deprivation on sleep rebound, host landing and blood-feeding propensity. Body and appendage position metrics revealed a clear distinction between the posture of mosquitoes in their putative sleep and awake states for multiple species, which correlate with a reduction in responsiveness to host cues. Sleep assessment informed by these posture analyses indicated significantly more sleep during periods of low activity. Nighttime and daytime sleep deprivation resulting from the delivery of vibration stimuli induced sleep rebound in the subsequent phase in day and night active mosquitoes, respectively. Lastly, sleep deprivation suppressed host landing in both laboratory and field settings and also impaired blood feeding of a human host when mosquitoes would normally be active. These results suggest that quantifiable sleep states occur in mosquitoes, and highlight the potential epidemiological importance of mosquito sleep. Mayaro Virus Infection Elicits an Innate Immune Response in Anopheles stephensi Mayaro virus (MAYV) is an arboviral pathogen in the genus Alphavirus that is circulating in South America with potential to spread to naive regions. MAYV is also one of the few viruses with the ability to be transmitted by mosquitoes in the genus Anopheles, as well as the typical arboviral transmitting mosquitoes in the genus Aedes. Few studies have investigated the infection response of Anopheles mosquitoes. In this study we detail the transcriptomic and small RNA responses of An. stephensi to infection with MAYV via infectious bloodmeal at 2, 7, and 14 days post infection (dpi). 487 unique transcripts were significantly regulated, 78 putative novel miRNAs were identified, and an siRNA response is observed targeting the MAYV genome. Gene ontology analysis of transcripts regulated at each timepoint suggested activation of the Toll pathway at 7 dpi and repression of pathways related to autophagy and apoptosis at 14 dpi. These findings provide a basic understanding of the infection response of An. stephensi to MAYV and help to identify host factors which might be useful to target to inhibit viral replication in Anopheles mosquitoes.

AUTHOR SUMMARYMayaro virus (MAYV) is a mosquito-borne Alphavirus responsible for outbreaks in South America and the Caribbean. In this study we infected Anopheles stephensi with MAYV and sequenced mRNA and small RNA to understand how MAYV infection impacts gene transcription and the expression of small RNAs in the mosquito vector. Genes involved with innate immunity and signaling pathways related to cell death are regulated in response to MAYV infection of An. stephensi, we also discover novel miRNAs and describe the expression patterns of miRNAs, siRNAs, and piRNAs following bloodmeal ingestion. These results suggest that MAYV does induce a molecular response to infection in its mosquito vector species. Local adaptation to novel habitats and climatic changes in a wild tomato species via selective sweeps. O_LIPositive Darwinian selection is the driving force underpinning local adaptation, and leaves footprints of selective sweeps at the underlying major genes. Quantifying the timing of selection and revealing the genetic bases of local adaptation in plants species occurring in steep and varying environmental gradients is crucial to predict a species ability to respond to climate changes and explore new niches.
C_LIO_LIWe use whole genome sequence data from six populations across three different habitats of the wild tomato species Solanum chilense, to 1) infer the past demographic history, and 2) search for genes under strong positive selection. We then correlate current and past climatic projections with the demographic history, allele frequencies in space and time, the age of selection events, and changes in habitat suitability.
C_LIO_LIWe find evidence for several selective sweeps targeting regulatory networks involved in root hair development in low altitude, and response to photoperiod and vernalization in high altitude populations. These sweeps occur in a concerted fashion in a given regulatory gene network at particular periods of substantial climatic change.
C_LIO_LIWe decipher the genetic bases and the timing of local adaptation during plant colonization of semi-arid habitats using a unique combination of genome scans for selection and modelling of past climatic data.
C_LI Critical re-evaluation of the slope factor of the operational model of agonism Background and PurposeAlthough being a relative term, agonist efficacy is a cornerstone in the proper assessment of agonist selectivity and signalling bias. The operational model of agonism (OMA) has become successful in the determination of agonist efficacies and ranking them. In 1985, Black et al. introduced the slope factor to the OMA to make it more flexible and allow for fitting steep as well as flat concentration-response curves.

Experimental ApproachFunctional analysis of OMA demonstrates that the slope factor implemented by Black et al. affects relationships among parameters of the OMA, leading to wrong estimates of operational efficacy and affinity of an agonist. In contrast, functional analysis of OMA with Hill coefficient as a slope factor has not this flaw.

Key ResultsFitting of the OMA with Black et al. slope factor to concentration-response curves of several models (homotropic allosteric modulation, substrate inhibition and non-competitive auto-inhibition) gives wrong estimates of operational efficacy and affinity. In contrast, fitting of the OMA modified by the Hill coefficient to concentration-response curves of these several models gives correct estimates of operational efficacy and affinity.

Conclusion and ImplicationsThe OMA modified by the Hill coefficient can be used for correct ranking of agonism and subsequent analysis, like quantification of signalling bias. Force-dependent remodeling of cytoplasmic ZO-1 condensates contributes to cell-cell adhesion Although the physiological importance of biomolecular condensates is widely recognized, how it is controlled in time and space during development is largely unknown. Here we show that a tight junction protein ZO-1 forms cytoplasmic condensates in the trophectoderm (TE) of the mouse embryo before E4.0. These disappear via dissolution, and ZO-1 accumulates at the cell junction as the blastocyst cavity grows and internal pressure on TE cells increases. In contrast, this dissolution was less evident in TE cells attached to the inner cell mass, as they receive weaker tensile forces. Furthermore, analyses using MDCK cells have demonstrated that the ZO-1 condensates are generated and maintained by liquid-liquid phase separation. Our study also highlights that the dynamics of these condensates depends on the physical environment via an interaction between ZO-1 and F-actin. We propose that the force-dependent regulation of ZO-1 condensation contributes to establishing robust cell-cell adhesion during early development. 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. 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. Early functional mismatch between breast cancer cells and their tumour microenvironment suppresses long term growth Cancer cells thrive embedded in a fine-tuned cellular and extracellular environment or tumour microenvironment (TME). There is a general understanding of a co-evolution between cancer cells and their surrounding TME, pointing at a functional connection between cancer cells characteristics and the perturbations induced in their surrounding tissue. However, whether this functional connection needs to be set from the start or if aggressive cancer cells can always be dominating their microenvironment has never been formally proven with a dedicated experimental setting where malignant cells can be challenged to grow in a different TME from the one they would naturally create. Here we generated an experimental setting where we transiently perturb the secretory profile of aggressive breast cancer cells without affecting their intrinsic growth ability. This led to the initial establishment of an atypical TME. Interestingly, even if initially tumours are formed, this atypical TME evolves to impair long term in vivo cancer growth. Using a combination of in vivo transcriptomics, protein arrays and in vitro co-cultures, we found that the atypical TME culminates in the infiltration of macrophages with STAT1high activity. These macrophages show strong anti-tumoural functions which reduce long-term tumour growth, despite lacking canonical M1 markers. Importantly, gene signatures of the mesenchymal compartment of the TME, as well as the anti-tumoural macrophages show striking prognostic power-correlating with less aggressive human breast cancers. Development of a Novel PROTAC using the Nucleic Acid Aptamer as a Targeting Ligand for Tumor Selective Degradation of Nucleolin PROteolysis TArgeting Chimeras (PROTACs) induce targeted protein degradation by hijacking the intracellular ubiquitin proteasome system, thus emerging as a new strategy for drug development. However, most PROTACs generated lack cell-type selectivity and have poorly soluble in water. To address this drawback, we developed a novel PROTAC ZL216 using aptamer AS1411 as a targeting ligand of nucleolin to conjugate with a small molecule ligand of E3 ligase VHL, which shows high aqueous solubility and serum stability. Based on the differential expression of nucleolin on the cell surface, ZL216 could bind to and internalize into breast cancer cells, but not normal breast cells. Furthermore, we revealed that ZL216 promoted the formation of a nucleolin-ZL216-VHL ternary complex in breast cancer cells and potently induced nucleolin degradation in vitro and in vivo. As a result, ZL216 inhibited the proliferation and migration of breast cancer cells. These studies demonstrate that in addition to peptides and small molecule compounds, nuclei acid aptamers can also be used to generate PROTACs, which broadens the toolbox constructing PROTACs and provides a promising strategy for development of tumor-selective PROTACs. 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.