Penalized Logistic Regression Analysis for Genetic Association Studies of Binary Phenotypes IntroductionIncreasingly, logistic regression methods for genetic association studies of binary phenotypes must be able to accommodate data sparsity, which arises from unbalanced case-control ratios and/or rare genetic variants. Sparseness leads to maximum likelihood estimators (MLEs) of log-OR parameters that are biased away from their null value of zero and tests with inflated type 1 errors. Different penalized-likelihood methods have been developed to mitigate sparse-data bias. We study penalized logistic regression using a class of log-F priors indexed by a shrinkage parameter m to shrink the biased MLE towards zero.

MethodsWe propose a two-step approach to the analysis of a genetic association study: first, a set of variants that show evidence of association with the trait is used to estimate m; and second, the estimated m is used for log-F -penalized logistic regression analyses of all variants using data augmentation with standard software. Our estimate of m is the maximizer of a marginal likelihood obtained by integrating the latent log-ORs out of the joint distribution of the parameters and observed data. We consider two approximate approaches to maximizing the marginal likelihood: (i) a Monte Carlo EM algorithm (MCEM) and (ii) a Laplace approximation (LA) to each integral, followed by derivative-free optimization of the approximation.

ResultsWe evaluate the statistical properties of our proposed two-step method and compared its performance to other shrinkage methods by a simulation study. Our simulation studies suggest that the proposed log-F -penalized approach has lower bias and mean squared error than other methods considered. We also illustrate the approach on data from a study of genetic associations with "super senior" cases and middle aged controls.

Discussion/ConclusionWe have proposed a method for single rare variant analysis with binary phenotypes by logistic regression penalized by log-F priors. Our method has the advantage of being easily extended to correct for confounding due to population structure and genetic relatedness through a data augmentation approach.
RapidoPGS: A rapid polygenic score calculator for summary GWAS data without a test dataset MotivationPolygenic scores (PGS) aim to genetically predict complex traits at an individual level. PGS are typically trained on genome-wide association summary statistics and require an independent test dataset to tune parameters. More recent methods allow parameters to be tuned on the training data, removing the need for independent test data, but approaches are computationally intensive. Based on fine-mapping principles, we present RapidoPGS, a flexible and fast method to compute PGS requiring summary-level GWAS datasets only, with little computational requirements and no test data required for parameter tuning.

ResultsWe show that RapidoPGS performs slightly less well than two out of three other widely-used PGS methods (LDpred2, PRScs, and SBayesR) for case-control datasets, with median r2 difference: -0.0092, -0.0042, and 0.0064, respectively, but up to 17,000-fold faster with reduced computational requirements. RapidoPGS is implemented in R and can work with user-supplied summary statistics or download them from the GWAS catalog.

Availability and implementationOur method is available with a GPL license as an R package from GitHub.
Quantifying the Morphology and Mechanisms of Cancer Progression in 3D in-vitro environments: Integrating Experiments and Multiscale Models Mathematical models of cancer growth have become increasingly more accurate both in the space and time domains. However, the limited amount of data typically available has resulted in a larger number of qualitative rather than quantitative studies. In the present study, we provide an integrated experimental-computational framework for the quantification of the morphological characteristics and the mechanistic modelling of cancer progression in 3D environments. The proposed framework allows for the calibration of multiscale, spatiotemporal models of cancer growth using state-of-the-art 3D cell culture data, and their validation based on the resulting experimental morphological patterns. Its implementation enabled us to pursue two goals; first, the quantitative description of the morphology of cancer progression in 3D cultures, and second, the relation of tumour morphology with underlying biophysical mechanisms that govern cancer growth and migration. We applied this framework to the study of the spatiotemporal progression of Triple Negative Breast Cancer cells cultured in Matrigel scaffolds, and validated the hypothesis of chemotactic migration using a multiscale Keller-Segel model. The results revealed transient, non-random spatial distributions of cancer cells that consist of clustered, and dispersion patterns. The proposed model was able to describe the general characteristics of the experimental observations and suggests that cancer cells exhibited chemotactic migration and accumulation, as well as random motion during the examined time period of development. To our knowledge, this is the first time that a multiscale model is used to quantify the relationship between the spatial patterns and the underlying mechanisms of cancer growth in 3D environments. 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. Can single-component protein condensates form multiphase architectures? Phase-separated biomolecular condensates that contain multiple coexisting phases are widespread in vitro and in cells. Multiphase condensates emerge readily within multi-component mixtures of biomolecules (e.g. proteins and nucleic acids) when the different components present sufficient physicochemical diversity (e.g. in inter-molecular forces, structure, and chemical composition) to sustain separate coexisting phases. Because such diversity is highly coupled to the solution conditions (e.g. temperature, pH, salt, composition), it can manifest itself immediately from the nucleation and growth stages of condensate formation, develop spontaneously due to external stimuli, or progressively as the condensates age. Here, we investigate thermodynamic factors that can explain the intrinsic transformation of single-component condensates into multiphase architectures during the nonequilibrium process of ageing. We develop a multiscale model that integrates atomistic simulations of proteins, sequence-dependent coarse-grained simulations of condensates, and a minimal model of dynamically ageing condensates with non-conservative inter-molecular forces. Our nonequilibrium simulations of condensate ageing predict that single-component condensates that are initially homogeneous and liquid-like can transform into gel-core/liquid-shell or liquid-core/gel-shell multiphase condensates as they age, due to gradual and irreversible enhancement of inter-protein interactions. The type of multiphase architecture is determined by the ageing mechanism, the molecular organization of the gel and liquid phases, and the chemical make up of the protein. Notably, we predict that inter-protein disorder-to-order transitions within the prion-like domains of intracellular proteins could lead to the required non-conservative enhancement of inter-molecular interactions. Our study, therefore, predicts a potential mechanism

Significance StatementBiomolecular condensates are highly diverse systems spanning not only homogeneous liquid droplets, but also gels, glasses, and even multiphase architectures that contain various coexisting liquid-like and/or gel-like inner phases. Multiphase architectures form when the different biomolecular components in a multi-component condensate establish sufficiently imbalanced inter-molecular forces to sustain different coexisting phases. While such a requirement seems, at first glance, impossible to fulfil for a condensate formed exclusively of chemically-identical proteins (i.e., single-component), our simulations predict conditions under which this may be possible. During condensate ageing, a sufficiently large imbalance in inter-molecular interactions can emerge intrinsically from the accumulation of protein structural transitions--driving even single-component condensates into nonequilibrium liquid-core/gel-shell or gel-core/liquid-shell multiphase architectures. An exon-intron split framework to prioritize miRNA-driven post-transcriptional regulatory signals and its application to study energy homeostasis in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Social behavioral profiling by unsupervised deep learning reveals a stimulative effect of dopamine D3 agonists on zebrafish sociality It has been a major challenge to systematically evaluate and compare how pharmacological perturbations influence social behavioral outcomes. Although some pharmacological agents are known to alter social behavior, precise description and quantification of such effects have proven difficult. The complexity of brain functions regulating sociality makes it challenging to predict drug effects on social behavior without testing in live animals, and most existing behavioral assays are low-throughput and provide only unidimensional readouts of social function. To achieve richer characterization of drug effects on sociality, we developed a scalable social behavioral assay for zebrafish named ZeChat based on unsupervised deep learning. High-dimensional and dynamic social behavioral phenotypes are automatically classified using this method. By screening a neuroactive compound library, we found that different classes of chemicals evoke distinct patterns of social behavioral fingerprints. By examining these patterns, we discovered that dopamine D3 agonists possess a social stimulative effect on zebrafish. The D3 agonists pramipexole, piribedil, and 7-hydroxy-DPAT-HBr rescued social deficits in a valproic acid-induced zebrafish autism model. The ZeChat platform provides a promising approach for dissecting the pharmacology of social behavior and discovering novel social-modulatory compounds. Quantifying the Morphology and Mechanisms of Cancer Progression in 3D in-vitro environments: Integrating Experiments and Multiscale Models Mathematical models of cancer growth have become increasingly more accurate both in the space and time domains. However, the limited amount of data typically available has resulted in a larger number of qualitative rather than quantitative studies. In the present study, we provide an integrated experimental-computational framework for the quantification of the morphological characteristics and the mechanistic modelling of cancer progression in 3D environments. The proposed framework allows for the calibration of multiscale, spatiotemporal models of cancer growth using state-of-the-art 3D cell culture data, and their validation based on the resulting experimental morphological patterns. Its implementation enabled us to pursue two goals; first, the quantitative description of the morphology of cancer progression in 3D cultures, and second, the relation of tumour morphology with underlying biophysical mechanisms that govern cancer growth and migration. We applied this framework to the study of the spatiotemporal progression of Triple Negative Breast Cancer cells cultured in Matrigel scaffolds, and validated the hypothesis of chemotactic migration using a multiscale Keller-Segel model. The results revealed transient, non-random spatial distributions of cancer cells that consist of clustered, and dispersion patterns. The proposed model was able to describe the general characteristics of the experimental observations and suggests that cancer cells exhibited chemotactic migration and accumulation, as well as random motion during the examined time period of development. To our knowledge, this is the first time that a multiscale model is used to quantify the relationship between the spatial patterns and the underlying mechanisms of cancer growth in 3D environments. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. Neurological mechanism of sensory deficits after developmental exposure to non-dioxin-like polychlorinated biphenyls (PCBs) in zebrafish larvae The most abundant polychlorinated biphenyl (PCB) congeners found in the environment and in humans are neurotoxic. This is of particular concern for early life stages because the exposure of the more vulnerable developing nervous system to chemicals can result in neurobehavioral disorders. To uncover currently unknown links between PCB target mechanisms and neurobehavioral deficits, we investigated the effects of the non-dioxin-like (NDL) congener PCB153 on neuronal morphology and synaptic transmission linked to the proper execution of a sensorimotor response using zebrafish as a vertebrate model. Zebrafish that were exposed during development to concentrations similar to those found in human cord blood and PCB contaminated sites showed a delay in startle response. A similar delay was observed for other NDL congeners but not for the potent dioxin-like congener PCB126. Morphological and biochemical data demonstrate that while exposure to PCB153 induced swelling of afferent sensory neurons, the disruption of dopaminergic and GABAergic signaling is the dominant mechanism associated with motor movement. The effects on important and broadly conserved signaling mechanisms in vertebrates suggest that NDL PCBs may contribute to neurodevelopmental abnormalities in humans and, with the startle response being critical for the survival of fish, to evolutionary adaptation in wildlife. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. Scaling effects of temperature on parasitism from individuals to host-parasite systems Parasitism is expected to change in a warmer future, but whether warming leads to substantial increases in parasitism remains unclear. Understanding how warming effects on parasitism in individual hosts (e.g., parasite load) translate to effects on population-level parasitism (e.g., prevalence, R0) remains a major knowledge gap. We analyzed the temperature dependence of parasitism at both host and population levels in thirteen empirical vector-borne host-parasite systems and found a strong, significant positive correlation between the thermal optima of individual- and population-level parasitism. We also found a significant, positive correlation in eleven environmentally-transmitted parasite systems, though several of these systems exhibited thermal optima >5{o}C apart between individual and population levels. Similarly, parasitism thermal optima were close to host performance thermal optima in vector-borne systems but not in environmentally-transmitted systems. We then adapted and simulated simple models for both transmission modes and found a similar pattern to the empirical systems: thermal optima in vector-borne systems were more strongly correlated across scales compared to environmentally-transmitted systems. Generally, our results suggest that information on the temperature-dependence, and specifically the thermal optimum, at either the individual- or population-level should provide a useful--though not quantitatively exact--baseline for predicting temperature dependence at the other level, especially in vector-borne parasite systems. Environmentally-transmitted parasitism may operate by a different set of rules, in which temperature-dependence is decoupled in some systems, requiring the need for trait-based studies of temperature-dependence at individual and population levels. Principles of mRNA targeting via the Arabidopsism6A-binding protein ECT2 Specific recognition of N6-methyladenosine (m6A) in mRNA by RNA-binding proteins containing a YT521-B homology (YTH) domain is important in eukaryotic gene regulation. The Arabidopsis YTH-domain protein ECT2 is thought to bind to mRNA at URU(m6A)Y sites, yet RR(m6A)CH is the canonical m6A consensus site in all eukaryotes and ECT2 functions require m6A binding activity. Here, we apply iCLIP (individual-nucleotide resolution cross-linking and immunoprecipitation) and HyperTRIBE (targets of RNA-binding proteins identified by editing) to define high-quality target sets of ECT2, and analyze the patterns of enriched sequence motifs around ECT2 crosslink sites. Our analyses show that ECT2 does in fact bind to RR(m6A)CH. Pyrimidine-rich motifs are enriched around, but not at m6A-sites, reflecting a preference for N6-adenosine methylation of RRACH/GGAU islands in pyrimidine-rich regions. Such motifs, particularly oligo-U and UNUNU upstream of m6A sites, are also implicated in ECT2 binding via its intrinsically disordered region (IDR). Finally, URUAY-type motifs are enriched at ECT2 crosslink sites, but their distinct properties suggest function as sites of competition between binding of ECT2 and as yet unidentified RNA-binding proteins. Our study provides coherence between genetic and molecular studies of m6A-YTH function in plants, and reveals new insight into the mode of RNA recognition by YTH-domain-containing proteins. Genome editing in the unicellular holozoan Capsaspora owczarzaki suggests a premetazoan role for the Hippo pathway in multicellular morphogenesis Animal development is mediated by a surprisingly small set of canonical signaling pathways such as Wnt, Hedgehog, TGF-beta, Notch and Hippo pathways. Although once thought to be present only in animals, recent genome sequencing has revealed components of these pathways in the closest unicellular relatives of animals. These findings raise questions about the ancestral functions of these developmental pathways and their potential role in the emergence of animal multicellularity. Here, we provide the first functional characterization of any of these developmental pathways in unicellular organisms by developing techniques for genetic manipulation in Capsaspora owczarzaki, a close unicellular relative of animals that displays aggregative multicellularity. We then use these tools to characterize the Capsaspora ortholog of the Hippo signaling nuclear effector YAP/TAZ/Yorkie (coYki), a key regulator of tissue size in animals. In contrast to what might be expected based on studies in animals, we show that coYki is dispensable for cell proliferation but regulates cytoskeletal dynamics and the three-dimensional (3D) shape of multicellular structures. We further demonstrate that the cytoskeletal abnormalities of individual coYki mutant cells underlie the abnormal 3D shape of coYki mutant aggregates. Taken together, these findings implicate an ancestral role for the Hippo pathway in cytoskeletal dynamics and multicellular morphogenesis predating the origin of animal multicellularity, which was co-opted during evolution to regulate cell proliferation. Expression of the cyanobacterial F0F1 ATP synthase inhibitor Atp{Theta} depends on small basic DNA-binding proteins and differential mRNA stability F0F1 ATP synthases produce ATP, the universal biological energy source. ATP synthase complexes on cyanobacterial thylakoid membranes use proton gradients generated either by photosynthesis or respiration. Atp{Theta} is an ATP synthase regulator in cyanobacteria which is encoded by the gene atpT. Atp{Theta} inhibits the hydrolysis of ATP (reverse reaction) that otherwise would occur under unfavorable conditions. In the cyanobacterium Synechocystis sp. PCC 6803, Atp{Theta} is expressed maximum in darkness but at very low levels under optimum phototrophic growth conditions or in the presence of glucose. DNA coimmunoprecipitation experiments followed by mass spectrometry identified the binding of the two transcriptional regulators cyAbrB1 and cyAbrB2 to the promoter and the histone-like protein HU to the 5UTR of atpT. Analyses of nucleotide substitutions in the promoter and GFP reporter assays identified a functionally relevant sequence motif resembling the HLR1 element bound by the RpaB transcription factor. Electrophoretic mobility shift assays confirmed interaction of cyAbrB1, cyAbrB2 and RpaB with the promoter DNA. However, overall the effect of transcriptional regulation was comparatively low. In contrast, atpT transcript stabilities differed dramatically, half-lives were 1.6 min in the light, 33 min in the dark and substantial changes were observed if glucose or DCMU were added. These findings show that basic transcriptional control of atpT involves nucleoid-associated DNA-binding proteins, positive regulation through RpaB, while the major effect on the condition-dependent regulation of atpT expression is mediated by controlling mRNA stability, which is related to the cellular redox and energy status.

IMPORTANCEF0F1 ATP synthases are protein complexes that produce ATP, the universal biological energy source in all kinds of organisms. Under unfavorable conditions, ATP synthases can operate in a futile reverse reaction, pumping protons while ATP is used up. Cyanobacteria perform plant-like photosynthesis but they cannot use the same mechanism as plants that inhibit chloroplast ATP synthases entirely during the night because respiratory and photosynthetic complexes are both located in the same membrane system. Atp{Theta} is a small peptide inhibitor of the reverse ATPase function in cyanobacteria encoded by the gene atpT. The production of Atp{Theta} is highly regulated to ensure that it is only synthetized when it is needed. In the here presented work we found that three transcription factors contribute to the regulation of atpT expression. However, we identified the control of mRNA stability as the major regulatory process governing atpT expression. Thus, it is the interplay between transcriptional and posttranscriptional regulation that position the Atp{Theta}-based inhibitory mechanism within the context of the cellular redox and energy balance. Programmed knockout mutation of liver fluke granulin, Ov-grn-1, impedes malignant transformation during chronic opisthorchiasis Infection with the food-borne liver fluke Opisthorchis viverrini is the principal risk factor for cholangiocarcinoma in the Mekong Basin countries of Thailand, Lao PDR, Vietnam, Myanmar and Cambodia. Using a novel model of CCA, involving infection with gene-edited liver flukes in the hamster during concurrent exposure to a nitrosamine, we explored the role of the fluke granulin-like growth factor Ov-GRN-1 in malignancy. We produced programmed gene knockout flukes ({Delta}Ov-grn-1) by delivery of a CRISPR/Cas9/gRNA system by electroporation. Genome sequencing confirmed Cas9-catalyzed mutations in the targeted genes, which was accompanied by rapid depletion of transcripts and the cognate proteins. Whereas Ov-grn-1 gene-edited parasites colonized the biliary tract and developed into adult flukes, less hepatobiliary tract disease manifested during chronic infection with{Delta} Ov-grn-1 worms in comparison to hamsters infected with control parasites. Specifically, immunohistochemical analysis of thin sections of livers revealed markedly less periductal fibrosis surrounding the flukes and less liver fibrosis globally during infection with{Delta} Ov-grn-1 genotype worms, minimal biliary epithelial cell proliferation, and markedly fewer mutations of TP53 in biliary epithelial cells. Moreover, fewer hamsters developed high-grade cholangiocarcinoma when infected with the{Delta} Ov-grn-1 flukes compared to controls. The clinically-relevant, pathophysiological phenotype of the hepatobiliary tract confirmed a role for this secreted growth factor in malignancy and morbidity during opisthorchiasis. The coiled coil of FtsLB is a structurally detuned element critical for modulating activation in bacterial cell division The FtsLB complex is a key regulator of bacterial cell division, existing in either an off or on state which supports the activation of septal peptidoglycan synthesis. In Escherichia coli, residues known to be critical for this activation are located in a region near the C-terminal end of the periplasmic coiled-coil domain of FtsLB, raising questions about the precise role of this conserved domain in the activation mechanism. Here, we investigate an unusual cluster of polar amino acids found within the core of the FtsLB coiled coil. We hypothesized that these amino acids likely reduce the structural stability of the domain and thus may be important for governing conformational changes. We found that mutating these positions to hydrophobic residues increased the thermal stability of FtsLB but caused cell division defects, suggesting that the coiled-coil domain is a "detuned" structural element. In addition, we identified suppressor mutations within the polar cluster, indicating that the precise identity of the polar amino acids is important for fine-tuning the structural balance between the off and on states. We propose a revised structural model of the tetrameric FtsLB (named the "Y-model") in which the periplasmic domain splits into a pair of coiled-coil branches. In this configuration, the hydrophilic terminal moieties of the polar amino acids remain more favorably exposed to water than in the original four-helix bundle model ("I-model"). We propose that a shift in this architecture, dependent on its marginal stability, is involved in activating the FtsLB complex and triggering septal cell wall reconstruction. Cytokinin oxidase/dehydrogenase family genes exhibit functional divergence and overlap in rice growth and development, especially in control of tillering Cytokinins play key roles in plant growth and development; hence, cytokinin biosynthesis and degradation have been extensively studied. Cytokinin oxidase/dehydrogenases (CKXs) are a group of enzymes that regulate oxidative cleavage to maintain cytokinin homeostasis. In rice, 11 OsCKX genes have been identified to date; however, most of their functions remain unknown. Here, we comprehensively analyzed the expression patterns and functions of OsCKX genes. Using CRISPR/Cas9 technology, we constructed mutants of all OsCKX genes to determine the functions of OsCKXs in rice development. The results revealed that the single osckx and higher-order osckx4 osckx9 mutant lines showed functional overlap and subfunctionalization. Notably, the osckx1 osckx2 and osckx4 osckx9 double mutants displayed contrasting phenotypic changes in tiller number and panicle size compared to the wild type. Moreover, we identified several genes with significantly altered expression in osckx4 and osckx9 single and double mutant plants. Many differentially expressed genes were found to be associated with auxin and cytokinin pathways. Additionally, the cytokinins in osckx4 osckx9 mutants were increased compared to the wild type. Overall, our findings provide new insights into the functions of OsCKX genes in rice growth and may be used as a foundation for future studies aimed at improving rice yield. Key words: Cytokinin, expression pattern, OsCKX, panicle, phenotype, rice, tiller

HighlightThe osckx4 osckx9 double mutant had a significantly greater number of tillers, whereas the osckx1 osckx2 double mutant showed the opposite phenotypic change, compared to the wild type SERCA function is impaired in skeletal and cardiac muscles from young DBA/2J mdx mice The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice. Species-specific developmental timing dictates expansion of the avian wing skeletal pattern A fundamental question in biology is how species-specific rates of embryonic development are controlled in differently sized species. To address this problem, we compared wing development in the quail and the larger chick. We reveal that pattern formation (humerus to digits) is faster in the quail wing as determined by the earlier activation of 5 Hox genes, termination of developmental organisers (Shh and Fgf8) and the laying down of the skeleton (Sox9). Using interspecies tissue grafts, we show that developmental timing can be reset during a critical window in which retinoic acid signalling is active. Accordingly, retinoic acid can switch developmental timing between the quail and chick, and also between the chick and the larger turkey. We show that the incremental growth rate is equivalent between all three species, suggesting that differences in the expansion of the skeletal pattern are governed primarily by the tempo of development. Our findings have implications for how development is timed both within and between different species. Host phenology can drive the evolution of intermediate virulence strategies in some obligate-killer parasites AO_SCPLOWBSTRACTC_SCPLOWThe traditional mechanistic trade-offs resulting in a negative correlation between transmission and virulence are the foundation of nearly all current theory on the evolution of parasite virulence. Several ecological factors have been shown to modulate the optimal virulence strategies predicted from mechanistic trade-off models, but these ecological factors have not yet been shown to be sufficient to explain the intermediate virulence strategies observed in any natural system. The timing of seasonal activity, or phenology, is a common factor that influences the types and impact of many ecological interactions but is difficult to incorporate into virulence evolution studies. We develop a mathematical model of a disease system with seasonal host activity to study the evolutionary consequences of host phenology on the virulence of obligate-killer parasite. Results from this model demonstrated that seasonal host activity is sufficient to drive the evolution of intermediate parasite virulence in some types of natural disease systems, even when a traditional mechanistic trade-off between transmission and virulence is not assumed in the modeling framework. The optimal virulence strategy in these systems can be determined by both the duration of the host activity period as well as the variation in the host emergence timing. Parasites with low virulence strategies are favored in environments with long host activity periods and in environments in which hosts emerge synchronously. The results demonstrate that host phenology can be sufficient to select for intermediate optimal virulence strategies, providing an alternative mechanism to account for virulence evolution in some natural systems. Macroscopic control of synchronous electrical signaling with chemically-excited gene expression Excitable cells can convert electrical signals into chemical outputs to facilitate the active transport of information across larger distances. This electrical-to-chemical conversion requires a tightly regulated expression of ion channels. Alterations of ion channel expression provide landmarks of numerous pathological diseases, such as cardiac arrhythmia, epilepsy, or cancer. Although the activity of ion channels can be locally regulated by external light or chemical stimulus, it remains challenging to coordinate the expression of ion channels on extended spatial-temporal scales in a non-invasive manner. Here, we have engineered yeast S. cerevisiae to read and convert local chemical concentrations into a dynamic electrical field distributed across cell populations. The core mechanism encodes a chemically-excitable dual-feedback gene circuit that precisely tunes the expression domain of potassium channels, globally coordinating cyclic firing of the plasma membrane potential (PMP). We demonstrate that this mechanism leverages an engineered constitutively open bacterial potassium channel KcsA to directly couple chemical stimuli with ion flux through gene expression and it can interface with the host ion channels through the pulsatile production of toxins. Our study provides a robust synthetic transcriptional toolbox underlying the conversion of local chemical environments into spatiotemporally organized electrical impulses for various cellular engineering, synthetic biology, and potential therapeutic applications. Integrated quality control of allele-specific copy numbers, mutations and tumour purity from cancer whole genome sequencing assays Cancer genomes contain thousands of somatic point mutations, chromosome copy alterations and more complex structural variants, which contribute to tumour growth and therapy response. Whole genome sequencing is a well established approach for somatic variant identification, but its broad application comes with complications, particularly in how proposed calls are quality assessed. To address this issue, we present CNAqc, a quantitative framework to quality control somatic mutations and allele-specific copy numbers, both in clonal and subclonal settings while accounting for variations in tumour purity, as commonly seen in bulk sampling. We test the model via extensive simulations, validate it using low-pass single-cell data, and apply it to 2778 single-sample PCAWG whole-genomes, 10 in-house multi-region whole-genomes and 48 TCGA whole-exomes. CNAqc is compatible with common bioinformatic pipelines and designed to support automated parameterization processes that are crucial in the era of large-scale whole genome sequencing. Analysis of Drosophila cardiac hypertrophy by micro-computerized tomography for genetic dissection of heart growth mechanisms Heart failure is often preceded by pathological cardiac hypertrophy, a thickening of the heart musculature driven by complex gene regulatory and signaling processes. The Drosophila heart has great potential as a genetic model for deciphering the underlying mechanisms of cardiac hypertrophy. However, current methods for evaluating hypertrophy of the Drosophila heart are laborious and difficult to carry out reproducibly. Here we demonstrate that micro-computerized tomography (microCT) is an accessible, highly reproducible method for non-destructive, quantitative analysis of Drosophila heart morphology and size. To validate our microCT approach for analyzing Drosophila cardiac hypertrophy, we show that expression of constitutively active Ras (Ras85DV12), previously shown to cause hypertrophy of the fly heart, results in significant thickening of both adult and larval heart walls when measured from microCT images. We then show using microCT analysis that genetic upregulation of store-operated Ca2+ entry (SOCE) driven by expression of constitutively active Stim (StimCA) or Orai (OraiCA) proteins also results in significant hypertrophy of the Drosophila heart, through a process that specifically depends on Orai Ca2+ influx channels. Intravital imaging of heart contractility revealed significantly reduced end diastolic dimensions in StimCA and OraiCA expressing hearts, consistent with the hypertrophic phenotype. These results demonstrate that increased SOCE activity is an important driver of hypertrophic cardiomyocyte growth, and demonstrate how microCT analysis combined with tractable genetic tools in Drosophila can be used to delineate molecular signaling processes that underlie cardiac hypertrophy and heart failure.

NEW AND NOTEWORTHYGenetic analysis of cardiac hypertrophy in Drosophila holds immense potential for the discovery of new therapeutic targets to prevent and treat heart failure. However, this potential has been hindered by a lack of rapid and effective methods for analysis of heart size in flies. Here we demonstrate that analysis of the Drosophila heart with micro-computerized tomography yields accurate and highly reproducible heart size measurements that can be used to efficiently analyze heart growth and cardiac hypertrophy in Drosophila. CXCR6 by increasing retention of memory CD8 T cells in the ovarian tumor microenvironment promotes immunosurveillance and control of ovarian cancer PurposeResident memory CD8 T cells owing to their ability to reside and persist in peripheral tissues, impart adaptive sentinel activity and amplify local immune response, have beneficial implications for tumor surveillance and control. The current study aims to clarify the less known chemotactic mechanisms that govern the localization, retention, and residency of memory CD8 T cells in the ovarian tumor microenvironment.

Experimental DesignRNA/FACS based profiling of chemokine receptor expression in CD8+ resident memory T cells in human ovarian cancer and analyze their association with survival. Analyze chemokine receptor role in anti-tumor response and control by resident memory T cells using prophylactic mice models of ovarian cancer, treated with adoptive transfer of OT1 T cells and vaccination with maraba virus-OVA to target Ovalbumin expressing tumor.

ResultsChemokine receptor profiling of CD8+CD103+ resident memory TILs in ovarian cancer patients revealed high expression of CXCR6. Analysis of the TCGA ovarian cancer database revealed CXCR6 to be associated with CD103 and increased patient survival. Functional studies in mouse models of ovarian cancer revealed that CXCR6 is a marker of resident, but not circulatory tumor-specific memory CD8 T cells. Knockout of CXCR6 in tumor-specific CD8 T cells showed reduced retention in tumor tissues leading to diminished resident memory responses and poor control of ovarian cancer

ConclusionsCXCR6 by promoting increased retention in tumor tissues serves a critical role in resident memory T cell-mediated immunosurveillance and control of ovarian cancer. Future studies warrant exploiting CXCR6 to promote resident memory response in cancers. A somatosensory computation that unifies limbs and tools It is often claimed that tools are embodied by the user, but whether the brain actually repurposes its body-based computations to perform similar tasks with tools is not known. A fundamental computation for localizing touch on the body is trilateration. Here, the location of touch on a limb is computed by integrating estimates of the distance between sensory input and its boundaries (e.g., elbow and wrist of the forearm). As evidence of this computational mechanism, tactile localization on a limb is most precise near its boundaries and lowest in the middle. We show that the brain repurposes trilateration to localize touch on a tool. In a large sample of participants, we found that localizing touch on a tool produced the signature of trilateration, with highest precision close to the base and tip of the tool. A computational model of trilateration provided a good fit to the observed localization behavior. To further demonstrate the computational plausibility of repurposing trilateration, we implemented it in a three-layer neural network that was based on principles of probabilistic population coding. This network determined hit location in tool-centered coordinates by using a tools unique pattern of vibrations when contacting an object. Simulations demonstrated the expected signature of trilateration, in line with the behavioral patterns. Our results have important implications for how trilateration may be implemented by somatosensory neural populations. We conclude that trilateration is a fundamental spatial computation that unifies limbs and tools.