Crystal structures of alphavirus nonstructural protein 4 (nsP4) reveal an intrinsically dynamic RNA-dependent RNA polymerase fold. Alphaviruses such as Ross River virus (RRV), chikungunya virus, Sindbis virus (SINV), and Venezuelan equine encephalitis virus are mosquito-borne pathogens that can cause arthritis or encephalitis diseases. Nonstructural protein 4 (nsP4) of alphaviruses possesses RNA-dependent RNA polymerase (RdRp) activity essential for viral RNA replication. No 3D structure has been available for nsP4 of any alphaviruses despite its importance for understanding alphaviral RNA replication and for the design of antiviral drugs. Here, we report crystal structures of the RdRp domain of nsP4 from both RRV and SINV determined at resolutions of 2.6 and 1.9 [A]. The structure of the alphavirus RdRp domain appears most closely related to RdRps from pestiviruses, noroviruses, and picornaviruses. Hydrogendeuterium exchange mass spectrometry (HDX-MS) and nuclear magnetic resonance (NMR) methods, showed that in solution, nsP4 is highly dynamic with an intrinsically disordered N-terminal domain. Both full-length nsP4 and the RdRp domain were capable to catalyze RNA polymerization. Structure-guided mutagenesis using a trans-replicase system identified nsP4 regions critical for viral RNA replication.

Key PointsO_LICrystal structures of alphavirus nsP4 RNA polymerase domain from RRV and SINV.
C_LIO_LInsP4 protein is highly dynamic with an intrinsically disordered N-terminal domain.
C_LIO_LIOptimized RNA elongation activity assay to facilitate antiviral discovery.
C_LI
T606-phosphorylation deprives the function of Kaiso as a transcription factor and cancer repressor It is well known that Kaiso protein encoded by ZBTB33 gene is a transcription repressor and that Kaiso-P120ctn interaction increases the shift of Kaiso from the nucleus into the cytoplasm. However, the regulatory mechanisms of Kaiso compartmentalization are far from clear. Here, we reported that AKT1 could phosphorylate 606-threonine residue (T606) within the RSSTIP motif of Kaiso in the cytoplasm. The T606-phosphorylated Kaiso (pT606-Kaiso) could directly bind to 14-3-3 family proteins and the depletion of T606 phosphorylation by T606A mutation abolished most of the Kaiso-14-3-3 binding. In addition, the Kaiso-P120ctn interaction was essential for the pT606-Kaiso accumulation in the cytoplasm. Notably, enforced 14-3-3{sigma} (SFN) overexpression could increase the pT606-Kaiso accumulation in the cytoplasm and de-repress the transcription of Kaiso target gene CDH1. Decreased amounts of both pT606-Kaiso and CDH1 proteins were frequently observed in human gastric cancer tissues relative to paired normal controls. The mRNA levels of 14-3-3{sigma} and Kaiso target gene CDH1 were positively and significantly correlated with each other in bioinformatics analyses using publicly available RNA-seq datasets for human normal tissues (n=11688, r=0.60, p<0.001) in the GTEx project and for cancer cell lines (n=1156, r=0.41, p<0.001) in the CCLE project. Furthermore, Kaiso T606A mutant (unable to be phosphorylated) significantly increased the migration and invasion of cancer cells in vitro as well as boosted the growth of these cells in vivo. In conclusion, Kaiso could be phosphorylated by AKT1 at the T606 and the pT606-Kaiso accumulates in the cytoplasm through binding to 14-3-3/P120ctn that de-represses the expression of Kaiso target gene CDH1 in normal tissues. Decreased Kaiso phosphorylation may contribute to the development of gastrointestinal cancer. The status of Kaiso phosphorylation is a determinant factor for the role of Kaiso in the development of cancer.
Convergent evolution of multiple mutations improves the viral fitness of SARS-CoV-2 variants by balancing positive and negative selection Multiple mutations have been seen to undergo convergent evolution in SARS-CoV-2 variants of concern. One such evolution occurs in Beta, Gamma, and Omicron variants at three amino acid positions K417, E484, and N501 in the receptor binding domain of the spike protein. We examined the physical mechanisms underlying the convergent evolution of three mutations K417T/E484K/N501Y by delineating the individual and collective effects of mutations on binding to angiotensin converting enzyme 2 receptor, immune escape from neutralizing antibodies, protein stability and expression. Our results show that each mutation serves a distinct function that improves virus fitness supporting its positive selection, even though individual mutations have deleterious effects that make them prone to negative selection. Compared to the wild-type, K417T escapes Class 1 antibodies, has increased stability and expression; however, it has decreased receptor binding. E484K escapes Class 2 antibodies; however, it has decreased receptor binding, stability and expression. N501Y increases receptor binding; however, has decreased stability and expression. When these mutations come together, the deleterious effects are mitigated due to the presence of compensatory effects. Triple mutant K417T/E484K/N501Y has increased receptor binding, escapes both Class 1 and Class 2 antibodies, and has similar stability and expression as that of the wild-type. These results show the implications of presence of multiple mutations on virus evolution that enhance viral fitness on different fronts by balancing both positive and negative selection and improves the chances of selection of mutations together. Algal photosystem I dimer and high resolution model of PSI:plastocyanin complex Photosystem I (PSI) enables photo-electron transfer and regulates photosynthesis in the bioenergetic membranes of cyanobacteria and chloroplasts. Being a multi-subunit complex, its macromolecular organization affects the dynamics of photosynthetic membranes. Here, we reveal a chloroplast PSI from the green alga Chlamydomonas reinhardtii that is organized as a homodimer, comprising 40 protein subunits with 118 transmembrane helices that provide scaffold for 568 pigments. Our cryo-EM structure identifies that the absence of PsaH and Lhca2 gives rise to a head-to-head relative orientation of the PSI-LHCI monomers in a way that is essentially different from the oligomer formation in cyanobacteria. The light-harvesting protein Lhca9 is the key element for mediating this dimerization. The interface between the monomers is lacking PsaH, and thus partially overlaps with the surface area that would bind one of the LHCII complexes in state transitions. We also define the most accurate available PSI-LHCI model at 2.3 [A] resolution, including a flexibly bound electron donor plastocyanin, and assign correct identities and orientations of all the pigments, as well as 621 water molecules that affect energy transfer pathways. A map of binary SARS-CoV-2 protein interactions implicates host immune regulation and ubiquitination Key steps in viral propagation, immune suppression, and pathology are mediated by direct, binary, physical interactions between viral and host proteins. To understand the biology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, we generated an unbiased systematic map of binary interactions between viral and host proteins, complementing previous co-complex association maps by conveying more direct mechanistic understanding and potentially enabling targeted disruption of direct interactions. To this end, we deployed two parallel strategies, identifying 205 virus-host and 27 intraviral binary interactions amongst 171 host and 19 viral proteins, and confirming high quality of these interactions via a calibrated orthogonal assay. Host proteins interacting with SARS-CoV-2 proteins are enriched in various cellular processes, including immune signaling and inflammation, protein ubiquitination, and membrane trafficking. Specific subnetworks provide new hypotheses related to viral modulation of host protein homeostasis and T-cell regulation. The binary virus-host protein interactions we identified can now be prioritized as targets for therapeutic intervention. More generally, we provide a resource of systematic maps describing which SARS-CoV-2 and human proteins interact directly. An individualized causal framework for learning intercellular communicationnetworks that define microenvironments of individual tumors Cells within a tumor microenvironment (TME) dynamically communicate and influence each others cellular states through an intercellular communication network (ICN). In cancers, intercellular communications underlie immune evasion mechanisms of individual tumors. We developed an individualized causal analysis framework for discovering tumor-specific ICNs. Using head and neck squamous cell carcinoma (HNSCC) tumors as a testbed, we first mined single-cell RNA-sequencing data to discover gene expression modules (GEMs) that reflect the states of transcriptomic processes within tumor and stromal single cells. By deconvoluting bulk transcriptomes of HNSCC tumors profiled by The Cancer Genome Atlas (TCGA), we estimated the activation states of these transcriptomic processes in individual tumors. Finally, we applied individualized causal network learning to discover an ICN within each tumor. Our results show that cellular states of cells in TMEs are coordinated through ICNs that enable multi-way communications among epithelial, fibroblast, endothelial, and immune cells. Further analyses of individual ICNs revealed structural patterns that were shared across subsets of tumors, leading to the discovery of 4 different subtypes of networks that underlie disparate TMEs of HNSCC. Patients with distinct TMEs exhibited significantly different clinical outcomes. Our results show that the capability of estimating individual ICNs reveals heterogeneity of ICNs and sheds light on the importance of intercellular communication in impacting disease development and progression. Frequency modulated timer regulates mammalian hibernation Mammalian hibernators decrease basal metabolism and body temperature (Tb) to minimize energy expenditure in harsh seasons. During hibernation, Tb drops to low temperature (<10 {degrees}C) and remains constant for days, known as deep torpor in small mammalian hibernators. Spontaneous interbout arousals interrupt torpor bouts, when Tb recovers to euthermic state ~37 {degrees}C. Torpor-interbout arousal event repeats during several months of hibernation. However, little is known about mechanisms governing Tb fluctuation across torpor-interbout arousal cycles during hibernation. Recent improvement in data-logging techniques enables us to monitor Tb for more than hundred days with high precision, opening up new avenues for quantitative analysis to address the principle governing Tb fluctuation. Here, we analyzed Tb fluctuation across torpor-interbout arousal cycle of Syrian hamster, which can hibernate in response to chronic cold and short photoperiod under a laboratory condition, using generalized harmonic analysis and discovered a model with frequency modulation quantitatively reproducing Tb fluctuation. This analysis also identified an unexpectedly longer period of 120-430 days as the period that modulates another period of several days, generating Tb fluctuation for Syrian hamster. We propose that concerted action of two endogenous periods governs torpor-interbout arousal cycles during hibernation. 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. Greater spear nosed bats in Panama do not use social proximity to improve foraging efficiency Group foraging animals may change their spatial relationships in response to the resource landscape. The strength of social network ties across foraging groups could then reflect the energy balance of individual group members through their effect on foraging efficiency. Greater spear-nosed bats (Phyllostomus hastatus) live in social groups where unrelated females develop group-specific social calls, defend each others offspring, and are hypothesized to coordinate social foraging to feed on balsa flowers. To test if social groups forage together, and if this increases foraging efficiency, we tracked 18 individuals from three social groups of P. hastatus on Isla Colon, Panama using high frequency GPS for 1 - 3 days each. We found that individuals commuted alone but then joined group members in small foraging patches with high densities of flowering balsas, 25 km from their roost on the mainland. Proximity-based social network measures did not predict foraging efficiency or energy expenditure, but individuals with higher social associations spent more time outside of the cave. Individuals did not coordinate movement, but close proximity among individuals increased resting bout duration and synchrony, particularly when in the same patch. These results show that stronger social proximity associations were more closely tied to resting behavior, and indicate that factors other than increased feeding efficiency may structure network relationships of group members while foraging.

Lay SummarySocial foraging can increase foraging efficiency on hard-to-find food sources. Miniaturized high-frequency GPS allowed us to test this on groups of greater spear-nosed bats. We estimated food intake rates, energy expenditure, and the strength of bat social relationships. In contrast to our expectations, we did not find a relationship between bat social relationships and energy expenditure or foraging efficiency. However, bats took longer resting breaks when group members were close. Why these bats commute 25 km individually across the ocean to then forage near each other thus requires further investigation. Diversity of SARS-CoV-2 genome among various strains identified in Lucknow, Uttar Pradesh Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has emerged as a significant challenge worldwide. Rapid genome sequencing of SARS-CoV-2 is going on across the globe to detect mutations and genomic modifications in SARS-CoV-2. In this study, we have sequenced twenty-three SARS-CoV-2 positive samples collected during the first pandemic from the state of Uttar Pradesh, India. We observed a range of already reported mutations (2-22), including; D614G, L452R, Q613H, Q677H, T1027I in the S gene; S194L in the N gene; Q57H, L106F, T175I in the ORF3. Few unreported mutations such as P309S in the ORF1ab gene; T379I in the N gene; and L52F, V77I in the ORF3a gene were also detected. Phylogenetic genome analysis showed similarity with other SARS-CoV-2 viruses reported from Uttar Pradesh. The observed mutations may be associated with SARS-CoV-2 virus pathogenicity or disease severity. Ultrahigh-Throughput Ambient MS: Direct Analysis at 22 Samples per Second by Infrared Matrix-Assisted Desorption Electrospray Ionization Mass Spectrometry Infrared Matrix-Assisted Laser Desorption Electrospray Ionization (IR-MALDESI) mass spectrometry is an ambient-direct sampling method being developed for high-throughput, label-free, biochemical screening of large-scale compound libraries. Here, we report the development of an ultrahigh-throughput continuous motion IR-MALDESI sampling approach capable of acquiring data at rates up to 22.7 samples per second in a 384-well microtiter plate. At top speed, less than 1% analyte carryover is observed from well-to-well and signal intensity relative standard deviations (RSD) of 11.5% and 20.9% for 3 M 1-hydroxymidazolam and 12 M dextrorphan, respectively, are achieved. The ability to perform parallel kinetics studies on 384 samples with ~30s time resolution using an isocitrate dehydrogenase 1 (IDH1) enzyme assay is shown. Finally, we demonstrate the repeatability and throughput of our approach by measuring 115,200 samples from 300 microtiter plate reads consecutively over 5.54 hours with RSDs under 8.14% for each freshly introduced plate. Taken together, these results demonstrate the use of IR-MALDESI at sample acquisition rates that surpass other currently reported direct sampling mass spectrometry approaches used for high throughput compound screening.

For Table of Contents Only

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=125 SRC="FIGDIR/small/465730v3_ufig1.gif" ALT="Figure 1">
View larger version (49K):
org.highwire.dtl.DTLVardef@5f916corg.highwire.dtl.DTLVardef@4473b4org.highwire.dtl.DTLVardef@82e8forg.highwire.dtl.DTLVardef@128b57b_HPS_FORMAT_FIGEXP  M_FIG C_FIG Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. Bacterial and fungal communities experience rapid secondary succession during the first year following a wildfire in a California chaparral. The rise in wildfire frequency and severity across the globe has increased interest in secondary succession. However, despite the role of soil microbial communities in the regeneration of post-fire vegetation and biogeochemical cycling, the lack of measurements immediately post-fire and at high temporal resolution has resulted in a limited understanding of microbial secondary succession. To fill this knowledge gap, we sampled soils at 17, 25, 34, 67, 95, 131, 187, 286, and 376 days after a Southern California wildfire in fire-adapted chaparral shrublands. We assessed bacterial and fungal biomass with qPCR of 16S and 18S and richness and composition with Illumina MiSeq sequencing of 16S and ITS2 amplicons. Fire severely reduced bacterial biomass by 47% and richness by 46%, but the impacts were stronger for fungi, with biomass decreasing by 86% and richness by 68%. While bacterial biomass and richness oscillated in response to precipitation, the declines in biomass and richness for both groups persisted for the entire year. In contrast to the unburned plots, burned fungal and bacterial communities experienced rapid succession, with 5-6 compositional turnover periods. Analogous to plants, "fire-loving" pyrophilous microbes increased post-fire and drove successional dynamics. Interestingly, many of the same pyrophilous taxa found in forests worldwide also dominated post-fire chaparral and traded off markedly in abundance over time. For fungi, an initial peak in the Basidiomycete yeast Geminibasidium traded off with the filamentous Ascomycetes Pyronema, Aspergillus, and Penicillium. For bacteria, the Proteobacteria genus Massilia dominated burned plots throughout the year, but the following most dominant bacteria, the Firmicute Bacillus and Proteobacteria Noviherbaspirillum, traded off in abundance over time. Our unprecedented high-resolution temporal sampling allowed us to capture post-fire microbial secondary successional dynamics and show that putative tradeoffs in thermotolerance, colonization, and competition among dominant pyrophilous microbes control microbial succession with possible implications for ecosystem function. Single-domain antibodies efficiently neutralize SARS-CoV-2 variants of concern, including Omicron variant Virus-neutralizing antibodies are one of the few treatment options for COVID-19. The evolution of SARS-CoV-2 virus has led to the emergence of virus variants with reduced sensitivity to some antibody-based therapies. The development of potent antibodies with a broad spectrum of neutralizing activity is urgently needed. Here we isolated a panel of single-domain antibodies that specifically bind to the receptor-binding domain of SARS-CoV-2 S glycoprotein. Three of the selected antibodies exhibiting most robust neutralization potency were used to generate dimeric molecules. We observed that these modifications resulted in up to a 200-fold increase in neutralizing activity. The most potent heterodimeric molecule efficiently neutralized each of SARS-CoV-2 variant of concern, including Alpha, Beta, Gamma, Delta and Omicron variants. This heterodimeric molecule could be a promising drug candidate for a treatment for COVID-19 caused by virus variants of concern. Impedance scaling for gold and platinum microelectrodes ObjectiveElectrical measurement of the activity of individual neurons is a primary goal for many invasive neural electrodes. Making these "single unit" measurements requires that we fabricate electrodes small enough so that only a few neurons contribute to the signal, but not so small that the impedance of the electrode creates overwhelming noise or signal attenuation. Thus, neuroelectrode design often must strike a balance between electrode size and electrode impedance, where the impedance is often assumed to scale linearly with electrode area.

Approach and main resultsHere we study how impedance scales with neural electrode area and find that the 1 kHz impedance of Pt electrodes (but not Au electrodes) transitions from scaling with area (r-2) to scaling with perimeter (r-1) when the electrode radius falls below 10 microns. This effect can be explained by the transition from planar to spherical diffusion behavior previously reported for electrochemical microelectrodes. Significance: These results provide important intuition for designing small, single unit recording electrodes. Specifically, for materials where the impedance is dominated by a pseudo-capacitance that is associated with a diffusion limited process, the total impedance will scale with perimeter rather than area when the electrode size becomes comparable with the diffusion layer thickness. For Pt electrodes this transition occurs around 10 um radius electrodes. At even lower frequencies (1 Hz) impedance approaches a constant. This transition to r-1 scaling implies that electrodes with a pseudo-capacitance can be made smaller than one might expect before thermal noise or voltage division limits the ability to acquire high-quality single-unit recordings. The Resting-State Causal Human Connectome is Characterized by Hub Connectivity of Executive and Attentional Networks We demonstrate a data-driven approach for calculating a "causal connectome" of directed connectivity from resting-state fMRI data using a greedy adjacency search and pairwise non-Gaussian edge orientations. We used this approach to construct n=442 causal connectomes. These connectomes were very sparse in comparison to typical Pearson correlation-based graphs (roughly 2.25% edge density) yet were fully connected in nearly all cases. Prominent highly connected hubs of the causal connectome were situated in attentional (dorsal attention) and executive (frontoparietal and cingulo-opercular) networks. These hub networks had distinctly different connectivity profiles: attentional networks shared incoming connections with sensory regions and outgoing connections with higher cognitive networks, while executive networks primarily connected to other higher cognitive networks and had a high degree of bidirected connectivity. Virtual lesion analyses accentuated these findings, demonstrating that attentional and executive hub networks are points of critical vulnerability in the human causal connectome. These data highlight the central role of attention and executive control networks in the human cortical connectome and set the stage for future applications of data-driven causal connectivity analysis in psychiatry. Re-convolving the compositional landscape of primary and recurrent glioblastoma reveals prognostic andtargetable tissue states Glioblastoma is an aggressive diffusely infiltrating neoplasm that spreads beyond surgical resection margins, where it intermingles with non-neoplastic brain cells. This complex microenvironment harbouring infiltrating glioma and non-neoplastic brain cells is the origin of tumor recurrence. Thus, understanding the cellular and molecular features of the glioma microenvironment is therapeutically and prognostically important. Here, we used single-nucleus RNA sequencing (snRNAseq) of primary and recurrent glioma to define three compositional  tissue-states rooted in cohabitation of cell-types and transcriptional states. These comprise a state featuring A) abundant normal brain microenvironment cells, B) reactive/inflamed brain microenvironment, and C) cellular/proliferative tumor. All these states exhibited variable degrees of infiltration by glioma cells. Spatial transcriptomics confirmed that the cell-types and transcriptomics states which compositionally cohabitate tissue states also colocalize in space. Tissue states are clinically significant because they correlate with radiographic, histopathologic, and prognostic features. We found that enrichment of tissue state B signature correlated with shorter survival. Importantly, we found that our compositionally defined tissue states are enriched in distinct metabolic pathways. One such pathway is fatty acid biosynthesis, which was enriched in tissue state B - a state enriched in recurrent glioblastoma and composed of astrocyte-like/mesenchymal glioma cells, reactive astrocytes resembling those seen in neurodegeneration, and monocyte-like myeloid cells. We showed that treating acute slices of GBM with a fatty acid synthesis inhibitor is sufficient to deplete the transcriptional signature of tissue state B. Our findings define a novel compositional approach to the glioma-infiltrated tissue which allows us to discover prognostic and targetable features, paving the way to new mechanistic and therapeutic discoveries. Quantifying the potential for red blood cell β-adrenergic sodium-proton exchangers to protect oxygen transport in hypoxic and hypercapnic white seabass White seabass (Atractoscion nobilis) are increasingly experiencing periods of low oxygen (O2; hypoxia) and high carbon dioxide (CO2, hypercapnia) due to climate change and eutrophication of the coastal waters of California. Hemoglobin (Hb) is the principal O2 carrier in the blood and in many teleost fishes Hb-O2 binding is compromised at low pH; however, the red blood cells (RBC) of some species regulate intracellular pH with adrenergically-stimulated sodium-proton-exchangers ({beta}-NHE). We hypothesized that RBC {beta}-NHEs in white seabass are an important mechanism that can protect the blood O2-carrying capacity during hypoxia and hypercapnia. We determined the O2-binding characteristics of white seabass blood, the response of RBCs to adrenergic stimulation, and quantified the protective effect of {beta}-NHE activity on Hb-O2 saturation. White seabass had typical teleost Hb characteristics, with a moderate O2 affinity (PO2 at half-saturation; P50 2.9 kPa) that was highly pH-sensitive (Bohr coefficient -0.92; Root effect 52%). The presence of RBC {beta}-NHEs was confirmed by functional, molecular and bioinformatic data and super-resolution imaging revealed, for the first time, the subcellular location of {beta}-NHE protein in vesicle-like structures and on the RBC membrane, and its translocation after adrenergic stimulation. The activation of RBC {beta}-NHEs increased Hb-O2 saturation by [~]8% in normoxia at 1 kPa PCO2, and by up to 20% in hypoxia. Our results confirm that RBC {beta}-NHE activity in white seabass has great potential to protect arterial O2 transport in environmentally relevant conditions of hypoxia and hypercapnia, but also reveal a potential vulnerability of fish to combinations of these stressors. 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. Inhibitory proteins block substrate access by occupying the active site cleft of Bacillus subtilis intramembrane metalloprotease SpoIVFB Intramembrane proteases function in numerous signaling pathways that impact health, but knowledge about regulation of intramembrane proteolysis is limited. We examined inhibition of intramembrane metalloprotease SpoIVFB by proteins BofA and SpoIVFA. We found that BofA residues in and near a predicted transmembrane segment are required for SpoIVFB inhibition, and cross-linking experiments indicated that this transmembrane segment occupies the SpoIVFB active site cleft. SpoIVFA is also required for SpoIVFB inhibition. The inhibitory proteins block access of the substrate N-terminal Proregion to the membrane-embedded SpoIVFB active site, based on additional cross-linking experiments; however, the inhibitory proteins did not prevent interaction between the substrate C-terminal region and the SpoIVFB soluble domain. A structural model was built of SpoIVFB in complex with BofA and parts of SpoIVFA and substrate, using partial homology and constraints from cross-linking and co-evolutionary analyses. The model predicts that conserved BofA residues interact to stabilize a transmembrane segment and a membrane-embedded C-terminal region. SpoIVFA is predicted to bridge the BofA C-terminal region and SpoIVFB, forming a membrane-embedded inhibition complex. Our results reveal a novel mechanism of intramembrane protease inhibition with clear implications for relief from inhibition in vivo and design of inhibitors as potential therapeutics. Tunable Cellular Localization and Extensive Cytoskeleton-Interplay of Reflectins Reflectin proteins are nature copolymers consist of repeated canonical domains, locating in a biophotonic system called Bragg lamellae and manipulating the dynamic structural coloration of iridocytes. Their biological functions are intriguing, but the underlying mechanism is not fully understood. In this study, reflectin A1, A2, B1 and C were found to present distinguished cyto-/nucleoplasmic localization preferences. Comparable intracellular localization can be reproduced by truncated reflectin variants, suggesting a conceivable evolutionary order among reflectin proteins. Secondly, the size-dependent access of reflectin variants into nucleus demonstrate a potential model of how reflectins get into Bragg lamellae. Thirdly, RfA1 was found to extensively interact with cytoskeleton, including its binding to actin and microtubule minus-end-directed movement. This implies cytoskeleton system plays fundamental role during the organization and transportation of reflectin proteins. Findings presented here recommend reflectins as programmable biomaterials which can be used to decipher their evolution processes, to delineate their biological mechanism, and to achieve tunable intracellular targeting as editable tags. Parallel evolution of urban-rural clines in melanism in a widespread mammal Urbanization is the dominant trend of global land use change. The replicated nature of environmental change associated with urbanization should drive parallel evolution, yet insight into the repeatability of evolutionary processes in urban areas has been limited by a lack of multi-city studies. Here we leverage community science data on coat color in >60,000 eastern gray squirrels (Sciurus carolinensis) across 43 North American cities to test for parallel clines in melanism, a genetically based trait associated with thermoregulation and crypsis. We show the prevalence of melanism was positively associated with urbanization as measured by impervious cover. Urban-rural clines in melanism were strongest in the largest cities with extensive forest cover and weakest or absent in cities with warmer winter temperatures, where thermal selection likely limits the prevalence of melanism. Our results suggest that novel traits can evolve in a highly repeatable manner among urban areas, modified by factors intrinsic to individual cities, including their size, land cover, and climate. Long-term conservation effects of protected areas in stochastic population dynamics Terrestrial and marine protected areas are essential tools in mitigating anthropogenic impacts and promoting population persistence and resource sustainability. Adequately implemented protected areas (PAs) promote long-term conservation effects. Stochasticity causes fluctuations in the conservation effects of PAs, and so it is important to investigate the variabilities of these conservation effects to inform their long-term conservation effects. To investigate long-term conservation effects, I develop and analyze new models of stochastic processes that encompass the fluctuations generated by demographic or environmental stochasticity in PAs management. The stochastic model is built upon individual processes. In the model, density-independent mortality, migration between PAs and non-PAs, and site preferences characterize the features of the PA. The effect of PAs size is also examined. The long-term conservation effects are quantified using the coefficient of variation (CV) of population size in PAs, where a lower CV indicates higher robustness in stochastic variations. Typically, the results from this study demonstrate that sufficiently reduced density-independent mortality in PAs and high site preference and immigration rate of PA are likely to decrease the CV. However, different types of stochasticity induce rather different consequences: under demographic stochasticity, the CV is always reduced when PAs increase the population size therein, but an increased population size by PAs does not always decrease the CV under environmental stochasticity. The deterministic dynamics of the model are investigated, facilitating effective management decisions. 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. Risk modelling of single cell populations revealed the heterogeneity of immune infiltration in hepatocellular carcinoma Hepatocellular carcinoma (HCC) is one of the most common primary hepatic malignancies. E2F transcription factors play an important role in the tumorigenesis and progression of HCC, mainly through the RB/E2F pathway. Prognostic models for HCC based on gene signatures have been developed rapidly in recent years; however, their discriminating ability at the single cell level remains elusive, which could reflect the underlying mechanisms driving the sample bifurcation. In this study, we constructed and validated a predictive model based on E2F expression, successfully stratifying HCC patients into two groups with different survival risks. Then we used a single-cell dataset to test the discriminating ability of the predictive model on infiltrating T cells, demonstrating remarkable cellular heterogeneity as well as altered cell fates. We identified distinct cell subpopulations with diverse molecular characteristics. We also found that the distribution of cell subpopulations varied considerably across onset stages among patients, providing fundamental basis for patient-oriented precision evaluation. Moreover, single-sample gene set enrichment analysis revealed that subsets of CD8+ T cells with significantly different cell adhesion level could be associated with different patterns of tumor cell dissemination. Therefore, our findings linked the conventional prognostic gene signature to the immune microenvironment and cellular heterogeneity at the single-cell level, thereby providing deeper insights into the understanding of HCC tumorigenesis. Sex differences and risk factors for bleeding in Alagille syndrome Spontaneous lethal bleeds are major cause of death in the pediatric liver disease Alagille syndrome (ALGS), yet risk factors and screening methods have not been established. We performed a systematic review and identified significantly more female than male patients with idiopathic intracranial hemorrhage (10:1). We investigated bleeding and vasculature in patients and a mouse model for ALGS (Jag1Ndr/Ndr mice) and asked whether phenotypes identified in mice could be detected in patients non-invasively. Jag1Ndr/Ndr mice bled spontaneously, exhibiting a thin skull and vascular defects including artery-vein crossings, tortuous vessels, capillary breakdown and CADASIL-like sparse vascular smooth muscle cell coverage which was aggravated by hypertension. Retinographs from patients confirmed tortuous blood vessels and artery-vein crossings in ALGS. In conclusion, Jag1Ndr/Ndr mice could be used to develop interventions for vascular defects in ALGS, and retinography could provide a non-invasive method for vascular analysis in these pediatric patients.