Analyses of infant vocal communication in mouse models of autism spectrum disorder using a novel unsupervised machine learning approach Infant crying is an innate communicative behavior that is frequently impaired in certain neurodevelopmental disorders (NDDs). Since advanced paternal age is a reported risk factor for NDDs in offspring, we evaluated the impact of a fathers age on early vocal development in C57BL/6J mice. We recorded and applied a unique combination of computational analyses to ultrasonic vocalizations (USVs) emitted by mouse pups sired by young and aged fathers. Our data showed that advanced paternal age reduced the number and duration of USVs, and altered the syllable composition in pups. Moreover, pups born to young fathers showed convergent vocal characteristics with a rich repertoire during postnatal development, while those born to aged fathers exhibited more divergent vocal patterns with limited repertoire. Principal component analysis in conjunction with clustering analysis demonstrated that pups from aged fathers deviated from typical trajectories of vocal development, which were considered as atypical individuals. Thus, our study indicates that advanced paternal age has a significant effect on offsprings early vocal development. It is suggested that the trajectories of vocal development could be a useful marker of the NDD-like phenotype associated with the advanced paternal age. In addition, our comprehensive computational analysis described here is an effective approach to characterize the altered individual diversity relevant to neurodevelopmental disorders.

One Sentence SummaryAdvanced paternal age affects vocal development in early postnatal mice, with more pups showing atypical developmental trajectories.
In Vivo Multi-Day Calcium Imaging of CA1 Hippocampus in Freely Moving Rats Reveals a High Preponderance of Place Cells with Consistent Place Fields Calcium imaging using GCaMP indicators and miniature microscopes has been used to image cellular populations during long timescales and in different task phases, as well as to determine neuronal circuit topology and organization. Because the hippocampus (HPC) is essential for tasks of memory, spatial navigation, and learning, calcium imaging of large populations of HPC neurons can provide new insight on cell changes over time during these tasks. All reported HPC in vivo calcium imaging experiments have been done in mouse. However, rats have many behavioral and physiological experimental advantages over mice. In this paper, we present the first (to our knowledge) in vivo calcium imaging from CA1 hippocampus in freely moving male rats. Using the UCLA Miniscope, we demonstrate that, in rat, hundreds of cells can be visualized and held across weeks. We show that calcium events in these cells are highly correlated with periods of movement, with few calcium events occurring during periods without movement. We additionally show that an extremely large percent of cells recorded during a navigational task are place cells (77.3{+/-}5.0%, surpassing the percent seen during mouse calcium imaging), and that these cells enable accurate decoding of animal position and can be held over days with consistent place fields in a consistent spatial map. A detailed protocol is included, and implications of these advancements on in vivo imaging and place field literature are discussed.

Significance statementIn vivo calcium imaging in freely moving animals allows the visualization of cellular activity across days. In this paper, we present the first in vivo Ca2+ recording from CA1 hippocampus in freely moving rats. We demonstrate that hundreds of cells can be visualized and held across weeks, and that calcium activity corresponds to periods of movement. We show that a high percentage (77.3{+/-}5.0%) of imaged cells are place cells, and that these place cells enable accurate decoding and can be held stably over days with little change in field location. Because the hippocampus is essential for many tasks involving memory, navigation, and learning, imaging of large populations of HPC neurons can shed new insight on cellular activity changes and organization.
Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

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

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

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. Somatic uniparental disomy mitigates the most damaging EFL1 allele combination in Shwachman-Diamond syndrome Shwachman-Diamond syndrome (SDS; OMIM: #260400) is caused by variants in SBDS (Shwachman-Bodian-Diamond syndrome gene), which encodes a protein that plays an important role in ribosome assembly. Recent reports suggest that recessive variants in EFL1 are also responsible for SDS. However, the precise genetic mechanism that leads to EFL1-induced SDS remains incompletely understood. Here we present three unrelated Korean SDS patients that carry biallelic pathogenic variants in EFL1 with biased allele frequencies, resulting from a bone marrow-specific somatic uniparental disomy (UPD) in chromosome 15. The recombination events generated cells that were homozygous for the relatively milder variant, allowing for the evasion of catastrophic physiological consequences. Still, the milder EFL1 variant was solely able to impair 80S ribosome assembly and induce SDS features in cell line, zebrafish, and mouse models. The loss of EFL1 resulted in a pronounced inhibition of terminal oligo-pyrimidine element-containing ribosomal protein transcript 80S assembly. Therefore, we propose a more accurate pathogenesis mechanism of EFL1 dysfunction that eventually leads to aberrant translational control and ribosomopathy. Age-associated B cells are long-lasting effectors that restrain reactivation of latent γHV68 Age-associated B cells (ABCs; CD19+CD11c+T-bet+) are a unique population that are increased in an array of viral infections, though their role during latent infection is largely unexplored. Here, we use murine gammaherpesvirus 68 ({gamma}HV68) to demonstrate that ABCs remain elevated long-term during latent infection and express IFN{gamma} and TNF. Using a strain of {gamma}HV68 that is cleared following acute infection, we show that ABCs persist in the absence of latent virus, though their expression of IFN{gamma} and TNF is decreased. With a fluorescent virus we demonstrate that ABCs are infected with {gamma}HV68 at similar rates to other previously activated B cells. We find that mice without ABCs display defects in anti-viral IgG2a/c antibodies and are more susceptible to {gamma}HV68 reactivation when challenged with heterologous infection. Together, these results indicate that ABCs are a persistent effector subset during latent viral infection that restrains {gamma}HV68 reactivation. Embryonic exposure to cannabidiol disrupts active neural circuits, an effect increased by Δ-9-tetrahydrocannabinol and involving CB1R and CB2R in zebrafish In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. A NANOG-pERK reciprocal regulatory circuit mediates Nanog autoregulation and ERK signaling dynamics. The self-renewal and differentiation potential of Embryonic stem cells (ESCs) is maintained by the regulated expression of core pluripotency factors. The expression level of core pluripotency factor Nanog is tightly regulated by a negative feedback autorepression loop. However, it remains unclear how the ESCs perceive the NANOG levels and execute autorepression. Here, we show that a dose-dependent induction of Fgfbp1 and Fgfr2 by NANOG activates an autocrine mediated ERK signaling in high-Nanog cells to trigger autorepression. pERK recruits NONO to Nanog locus to repress transcription by preventing POL2 loading. The Nanog autorepression process establishes a self-perpetuating NANOG-pERK reciprocal regulatory circuit. We further demonstrate that the reciprocal regulatory circuit induces the pERK heterogeneity and ERK signaling dynamics in pluripotent stem cells. 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. "Neighbourhood watch" model: embryonic epiblast cells assess positional information in relation to their neighbours In many developing and regenerating systems, tissue pattern is established through gradients of informative morphogens, but we know little about how cells interpret these. Using experimental manipulation of early chick embryos including misexpression of an inducer (VG1 or ACTIVIN) and an inhibitor (BMP4), we test two alternative models for their ability to explain how the site of primitive streak formation is positioned relative to the rest of the embryo. In one model, cells read morphogen concentrations cell-autonomously. In the other, cells sense changes in morphogen status relative to their neighbourhood. We find that only the latter model can account for the experimental results, including some counter-intuitive predictions. This mechanism (which we name "neighbourhood watch" model) illuminates the classic "French Flag Problem" and how positional information is interpreted by a sheet of cells in a large developing system.

Summary statementIn a large developing system, the chick embryo before gastrulation, cells interpret gradients of positional signals relative to their neighbours to position the primitive streak, establishing bilateral symmetry. First Report of Lactoria cornuta from Karnataka, West coast of India, Role of tropical cyclones in Range Extension India is known for its rich biodiversity and is fortunate to have several endemic species from different classes of vertebrates. India is home to 7.5% of the global fish diversity, with 91 endemic species of ray finned fishes, the actinopterygians. Some fish species have never been reported until recently, one such example is the long-horned Cowfish (Ostracidae), best known for inhabiting only certain areas of the Indo-Pacific oceans. However, this fish has been reported from Bay of Bengal, Costs of Goa and Lakshadweep. This is the first time this fish has been found on the Coasts of Karnataka, west coast of India. Given its morphology, the migration seems highly improbable since it can only endure lethargic swimming. In this article, we have explored the role of other external forces that could have contributed to its range extension in the east and west coast of India. In this paper we are exploring the role of Super Cyclones as a vital force in determining the expansion of its range. The fact of collection longhorn cowfish for the first time from the coast of Karwar coincides with the passage of Cyclone Ockhi on the west coast of India this suggests the possible role of super cyclone in its range extension. A year later, another specimen was recovered from a location very close to the first, which indicates the successful establishment of Lacturia cornuta in its new environment. Related morphometric and meristic studies of our specimens are consistent with previous studies reported in the Bay of Bengal. Modulation of P2X4 pore closure by magnesium,potassium, and ATP The P2X4 receptor plays a prominent role in cellular responses to extracellular ATP. Through classical all-atom molecular dynamics (MD) simulations totaling 24 {micro}s we have investigated how metal-complexed ATP stabilizes the channels open state and prevents its closing. We have identified two metal-binding sites, magnesium (Mg2+) and potassium (K+), one at the intersection of the three subunits in the ectodomain (MBS1) and the second one near the ATP binding site (MBS2), similar to those characterized in Gulf coast P2X. Our data indicate that when Mg2+ and K+ ions are complexed with ATP, the channel is locked into an open state. Interestingly, irrespective of the number of bound ATP molecules, Mg2+ ions bound to the MBS2 resisted collapsing of the open state protein to a closed state by stabilizing the ATP-protein interactions. However, when Mg2+ in the MBS2 was replaced with K+ ions, as might be expected when in equilibrium with an extracellular solution, the interactions between the subunits were weakened and we found evidence of pore collapse. This collapse was apparent when fewer than two ATP were bound to MBS2 in the presence of K+. Therefore, the different capacities of common cations to stabilize the channel may underlie a mechanism governing P2X4 channel gating in physiological systems. This study provides structural insights into the differential modulation of ATP activation of P2X4 by Mg2+ and K+. Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

Here, we designed, prepared, and characterized two combinatorial protein libraries consisting of randomized proteins, each 105 residues in length. The first library constructed proteins from the entire canonical alphabet of 20 amino acids. The second library used a subset of only 10 residues (A,S,D,G,L,I,P,T,E,V) that represent a consensus view of plausibly available amino acids through prebiotic chemistry. Our study shows that compact conformations resistant to proteolysis are (i) abundant (up to 40%) in random sequence space, (ii) independent of general Hsp70 chaperone system activity, and (iii) not granted solely by "late" and complex amino acid additions. The Hsp70 chaperone system effectively increases solubility and refoldability of the canonical alphabet but has only a minor impact on the "early" library. The early alphabet proteins are inherently more soluble and refoldable, possibly assisted by the cell-like environment in which these assays were performed.

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. Transfer learning framework for cell segmentation with incorporation of geometric features With recent advances in multiplexed imaging and spatial transcriptomic and proteomic technologies, cell segmentation is becoming a crucial step in biomedical image analysis. In recent years, Fully Convolutional Networks (FCN) have achieved great success in nuclei segmentation in in vitro imaging. Nevertheless, it remains challenging to perform similar tasks on in situ tissue images with more cluttered cells of diverse shapes. To address this issue, we propose a novel transfer learning, cell segmentation framework incorporating shape-aware features in a deep learning model, with multi-level watershed and morphological post-processing steps. Our results show that incorporation of geometric features improves generalizability to segmenting cells in in situ tissue images, using solely in vitro images as training data. Targeted biallelic integration of an inducible Caspase 9 suicide gene for safer cellular therapies prevents development of drug-resistant escapees in human iPSCs Drug-inducible suicide systems may help to minimize risks of cellular therapies due to the tumor forming potential of human induced pluripotent stem cells (hiPSCs). Recent research challenged the usefulness of such systems since rare drug-resistant subclones were observed that showed elimination or silencing of the transgene.

We have introduced a drug-inducible Caspase9 suicide system (iCASP9) into the AAVS1 safe harbor locus of hiPSCs. In these cells, apoptosis could be efficiently induced in vitro. In mice, drug treatment generally led to rapid elimination of teratomas, but individual animals subsequently formed tumor tissue from monoallelic iCASP9 hiPSCs. Very rare drug-resistant subclones of monoallelic iCASP9 hiPSCs appeared in vitro with frequencies of ~ 3x10-8. Transgene elimination, presumably via Loss of Heterozygosity (LoH), or methylation of the CAG promoter but not methylation of the O_SCPLOWPPPC_SCPLOWO_SCPCAP1C_SCPCAPO_SCPLOWRC_SCPLOWO_SCPCAP12C_SCPCAPO_SCPLOWCC_SCPLOW locus were identified as underlying mechanisms. In contrast, we never observed any escapees from biallelic iCASP9 hiPSCs, even after treatment of up to 0.8 billion hiPSCs.

In conclusion, biallelic integration of an iCASP9 system in the AAVS1 locus may substantially contribute to the safety level of iPSC-based therapies, which should be calculated by relating clonal escapee frequencies to the cell number in tumors of a size that is readily detectable during routine screening procedures. On the incongruence of genotype-phenotype and fitness landscapes The mapping from genotype to phenotype to fitness typically involves multiple nonlinearities that can transform the effects of mutations. For example, mutations may contribute additively to a phenotype, but their effects on fitness may combine non-additively because selection favors a low or intermediate value of that phenotype. This can cause incongruence between the topographical properties of a fitness landscape and its underlying genotype-phenotype landscape. Yet, genotype-phenotype landscapes are often used as a proxy for fitness landscapes to study the dynamics and predictability of evolution. Here, we use theoretical models and empirical data on transcription factor-DNA interactions to systematically study the incongruence of genotype-phenotype and fitness landscapes when selection favors a low or intermediate phenotypic value. Using the theoretical models, we prove a number of fundamental results. For example, selection for low or intermediate phenotypic values does not change simple sign epistasis into reciprocal sign epistasis, yet it changes reciprocal sign epistasis into simple sign epistasis and no sign epistasis with equal probability. More broadly, we show that such selection tends to create fitness landscapes that are more rugged than the underlying genotype-phenotype landscape, but this increased ruggedness typically does not frustrate adaptive evolution because the local adaptive peaks in the fitness landscape tend to be nearly as tall as the global peak. Many of these results carry forward to the empirical genotype-phenotype landscapes, which may help to explain why low- and intermediate-affinity transcription factor-DNA interactions are so prevalent in eukaryotic gene regulation.

AUTHOR SUMMARYHow do mutations change phenotypic traits and organismal fitness? This question is often addressed in the context of a classic metaphor of evolutionary theory -- the fitness landscape. A fitness landscape is akin to a physical landscape, in which genotypes define spatial coordinates, and fitness defines the elevation of each coordinate. Evolution then acts like a hill-climbing process, in which populations ascend fitness peaks as a consequence of mutation and selection. It is becoming increasingly common to construct such landscapes using experimental data from high-throughput sequencing technologies and phenotypic assays, in systems such as macromolecules and gene regulatory circuits. Although these landscapes are typically defined by molecular phenotypes, and are therefore more appropriately referred to as genotype-phenotype landscapes, they are often used to study evolutionary dynamics. This requires the assumption that the molecular phenotype is a reasonable proxy for fitness, which need not be the case. For example, selection may favor a low or intermediate phenotypic value, causing incongruence between a fitness landscape and its underlying genotype-phenotype landscape. Here, we study such incongruence using a diversity of theoretical models and experimental data from gene regulatory systems. We regularly find incongruence, in that fitness landscapes tend to comprise more peaks than their underlying genotype-phenotype landscapes. However, using evolutionary simulations, we show that this increased ruggedness need not impede adaptation. Material-mediated histogenesis using mechano-chemically microstructured cell niches Stem-cell derived tissue models are commonly cultured under globally-delivered stimuli that trigger histogenesis via self-organizing activity. However, the culture of such tissue models is prone to stochastic behavior, limiting the reproducibility of cellular composition and resulting in non-physiological architectures. To overcome these shortcomings, we developed a method for printing cell niche microenvironments with microstructured cues that mediate local histogenic processes, including mechanosensing and differentiation of selected cell types. Microstructured cues include independently tunable mechano-chemical properties, with conjugated peptides, proteins, and morphogens across a range of Youngs moduli. By rationally designing niches, we mediate the structure of tissues derived from stem-cell-progenitor sources, including a bone-fat assembly from stromal mesenchyme, and embryonic tissues derived from hiPSC. We show that microstructured cues can recapitulate mechano-chemical signals resembling early embryonic histogenesis. This outcome includes a role for niche mechanics in human embryonic organization, where soft niche mechanics bias markers of mesendodermal differentiation and epithelial-to-mesenchymal-transition (EMT), as well as a demonstration of a material-mediated morphogen signaling centers able to induce foci of mesenchymal and EMT differentiation. Thus, microstructured materials can mediate local histogenic processes to enhance the structure and composition of tissue models. Flexible reprogramming of Pristionchus pacificus motivation for attacking Caenorhabditis elegans in predator-prey competition Animals with diverse diets must adapt their food priorities to a wide variety of environmental conditions. This diet optimization problem is especially complex for predators that compete with prey for food. Although predator-prey competition is widespread and ecologically critical, it remains difficult to disentangle predatory and competitive motivations for attacking competing prey. Here, we dissect the foraging decisions of the omnivorous nematode Pristionchus pacificus to reveal that its seemingly failed predatory attempts against Caenorhabditis elegans are actually motivated acts of efficacious territorial aggression. While P. pacificus easily kills and eats larval C. elegans with a single bite, adult C. elegans typically survives and escapes from bites. However, nonfatal biting can provide competitive benefits by reducing access of adult C. elegans and its progeny to bacterial food that P. pacificus also eats. We show that P. pacificus considers costs and benefits of both predatory and territorial outcomes to decide which food goal, prey or bacteria, should guide its motivation for biting. These predatory and territorial motivations impose different sets of rules for adjusting willingness to bite in response to changes in bacterial abundance. In addition to biting, predatory and territorial motivations also influence which search tactic P. pacificus uses to increase encounters with C. elegans. When treated with an octopamine receptor antagonist, P. pacificus switches from territorial to predatory motivation for both biting and search. Overall, we demonstrate that P. pacificus assesses alternate outcomes of attacking C. elegans and flexibly reprograms its foraging strategy to prioritize either prey or bacterial food. Limited selection towards lighter seeds in Danish grasslands over an eight-year period Time-series analysis of cover data from 236 Danish grassland sites demonstrated community selection on seed mass. Across four grassland habitat types and during an eight-year period there was significant community selection towards lighter seeds. In the analysis, the continuous seed size variable was used to group plant species into functional types. This method was chosen in order to account for the sampling process of the abundance data, and is in contrast to most other analyses of trait selection, where the community weighted mean of the traits is used as the dependent variable. Generally, little is known on the underlying causes of possible plant trait community selection in semi-natural grasslands that are not subjected to successional processes. Membrane remodeling due to mixture of multiple curvature proteins We present an extension of the Monte Carlo based mesoscopic membrane model, where the membrane is represented as a dynamically triangulated surface and the proteins are modeled as anisotropic inclusions formulated as in-plane nematic field variables adhering to the deformable elastic sheet. In the extended model, we have augmented the Hamiltonian to study membrane deformation due to a mixture of multiple types of curvature generating proteins. This feature opens the door for understanding how multiple kinds of curvature-generating proteins may be working in a coordinated manner to induce desired membrane morphologies. For example, among other things, we study membrane deformations and tubulation due to a mixture of positive and negative curvature proteins as mimics of various proteins from BAR domain family working together for curvature formation and stabilization. We also study the effect of membrane anisotropy, which manifests as membrane localization and differential binding affinity of a given curvature protein, leading to insights into the tightly regulated cargo sorting and transport processes. Our simulation results show different morphology of deformed vesicles that depend on the curvatures and densities of the participating proteins as well as on the protein-protein and membrane-proteins interactions.

SignificanceMembrane remodeling requires highly orchestrated interactions between several types of lipids and curvature proteins. Experimentally probing the membrane deformation intermediates is non-trivial. For example, it is not known exactly how two or more different kinds of BAR domain proteins work in concert to induce and stabilize certain membrane curvatures. In this work, we use mesoscopic continuum modeling and explore the factors that induce and stabilize a range of membrane deformation due to a mixture of positive and negative curvature proteins. Marburg and Ebola virus infections elicit a complex, muted inflammatory state in bats The Marburg and Ebola filoviruses cause a severe, often fatal, disease in humans and nonhuman primates but have only subclinical effects in bats, including Egyptian rousettes, which are a natural reservoir of Marburg virus. A fundamental question is why these viruses are highly pathogenic in humans but fail to cause disease in bats. To understand how bats resist the disease caused by filoviruses, we infected one cohort of Egyptian rousette bats with Marburg virus and another cohort with Ebola virus and harvested multiple tissues for mRNA expression analysis. While virus transcripts were found primarily in the liver, Principal component analysis (PCA) revealed coordinated changes across multiple tissues. Gene signatures in kidney and liver pointed at induction of vasodilation, reduction of coagulation and changes in the regulation of iron metabolism. Signatures of immune response detected in spleen and liver indicated a robust anti-inflammatory state signified by macrophages in the M2 state and an active T cell response. Many of the responsive genes were found to be evolutionarily divergent, providing a framework for understanding the differences in outcomes of filovirus infections between bats and humans. In this study, we outline multiple interconnected pathways that respond to infection by MARV and EBOV, providing insights into the complexity of the mechanisms that enable bats to resist the disease caused by filoviral infections. The results have the potential to aid in the development of new strategies to effectively mitigate and treat the disease caused by these viruses in humans. 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. Improving GWAS discovery and genomic prediction accuracy in Biobank data Genetically informed and deep-phenotyped biobanks are an important research resource. The cost of phenotyping far outstrips that of genotyping, and therefore it is imperative that the most powerful, versatile and efficient analysis approaches are used. Here, we apply our recently developed Bayesian grouped mixture of regressions model (GMRM) in the UK and Estonian Biobanks and obtain the highest genomic prediction accuracy reported to date across 21 heritable traits. On average, GMRM accuracies were 15% (SE 7%) greater than prediction models run in the LDAK software with SNP annotation marker groups, 18% (SE 3%) greater than a baseline BayesR model without SNP markers grouped into MAF-LD-annotation categories, and 106% (SE 9%) greater than polygenic risk scores calculated from mixed-linear model association (MLMA) estimates. For height, the prediction accuracy R2 was 47% in a UK Biobank hold-out sample, which was 76% of the estimated [Formula]. We then extend our GMRM prediction model to provide MLMA SNP marker estimates for GWAS discovery, which increased the independent loci detected to 7,910 in unrelated UK Biobank individuals, as compared to 5,521 from BoltLMM and 5,727 from Regenie, a 43% and 38% increase respectively. The average{chi} 2 value of the leading markers was 34% (SE 5.11) higher for GMRM as compared to Regenie, and increased by 17% for every 1% increase in prediction accuracy gained over a baseline BayesR model across the traits. Thus, we show that modelling genetic associations accounting for MAF and LD differences among SNP markers, and incorporating prior knowledge of genomic function, is important for both genomic prediction and for discovery in large-scale individual-level biobank-scale studies. Native proline-rich motifs exploit sequence context to target actin-remodeling Ena/VASP proteins The human proteome is replete with short linear motifs (SLiMs) of 4-6 residues that are critical for protein-protein interactions, yet the importance of the sequence surrounding such motifs is underexplored. We devised a proteomic screen to systematically examine the influence of SLiM sequence context on protein-protein interactions. Focusing on the EVH1 domain of ENAH, an actin regulator that is upregulated in invasive cancers, we screened 36-residue proteome-derived peptides for binding. We discovered a pocket on the ENAH EVH1 domain that diverged from its orthologs to recognize extended SLiMs, and we found that proteins with two EVH1-binding SLiMs can wrap around a single domain. We also found that the ciliary protein PCARE uses an extended 23-residue region to obtain higher affinity than any known ENAH EVH1-binding motif. Our screen provides a way to uncover the effects of broader proteomic context on motif-mediated interactions, revealing diverse mechanisms of contextual control over EVH1 interactions and establishing that SLiMs cant be fully understood outside of their native context. Genetic variation associated with plastic and homeostatic growth responses to drought in Arabidopsis O_LINatural genetic variation influences plant responses to environmental stressors. However, the extent to which such variation underlies plastic versus homeostatic response phenotypes has been little studied.
C_LIO_LIWe quantified the extent of drought-induced changes in leaf area in a set of Iberian Arabidopsis accessions and then performed association studies correlating variation in this phenotype with genomic and transcriptomic variation.
C_LIO_LIDrought-induced plastic reductions in relative leaf area typified accessions originating from productive environments. In contrast, homeostasis in relative leaf area typified accessions originating from unproductive environments. Genome-Wide Association Studies (GWAS), Transcriptome Wide Association Studies (TWAS), and expression GWAS (eGWAS) highlighted the importance of auxin-related processes and, in particular, the potential role of the SMALL AUXIN UP RNA 26 (SAUR26) gene in conferring leaf area plasticity. Homeostatic responses in relative leaf area were associated with a diverse gene set and positively associated with a higher intrinsic water use efficiency (WUEi) as confirmed in a TWAS metanalysis of previously published {delta}13C measurements.
C_LIO_LIWe have identified not only candidate "plasticity genes" but also candidate "homeostasis genes" controlling leaf area. Our results demonstrate the value of a combined GWAS, TWAS, and eGWAS approach to identify mechanisms underlying phenotypic responses to stress.
C_LI

One-sentence summaryInformation on phenotype, genotype, and transcript abundance is integrated to identify both plasticity and homeostasis genes and processes associated with local adaptation to drought stress in Arabidopsis accessions of the Iberian Peninsula. Human papillomavirus integration transforms chromatin to drive oncogenesis Human papillomavirus (HPV) drives almost all cervical cancers and up to [~]70% of head and neck cancers. Frequent integration into the host genome occurs only for tumourigenic strains of HPV. We hypothesized that changes in the epigenome and transcriptome contribute to the tumourigenicity of HPV. We found that viral integration events often occurred along with changes in chromatin state and expression of genes near the integration site. We investigated whether introduction of new transcription factor binding sites due to HPV integration could invoke these changes. Some regions within the HPV genome, particularly the position of a conserved CTCF binding site, showed enriched chromatin accessibility signal. ChIP-seq revealed that the conserved CTCF binding site within the HPV genome bound CTCF in 4 HPV+ cancer cell lines. Significant changes in CTCF binding pattern and increases in chromatin accessibility occurred exclusively within 100 kbp of HPV integration sites. The chromatin changes co-occurred with out-sized changes in transcription and alternative splicing of local genes. We analyzed the essentiality of genes upregulated around HPV integration sites of The Cancer Genome Atlas (TCGA) HPV+ tumours. HPV integration upregulated genes which had significantly higher essentiality scores compared to randomly selected upregulated genes from the same tumours. Our results suggest that introduction of a new CTCF binding site due to HPV integration reorganizes chromatin and upregulates genes essential for tumour viability in some HPV+ tumours. These findings emphasize a newly recognized role of HPV integration in oncogenesis. Identifying complex motifs in massive omics data with a variable-convolutional layer in deep neural network Motif identification is among the most common and essential computational tasks for bioinformatics and genomics. Here we proposed a novel convolutional layer for deep neural network, named Variable Convolutional (vConv) layer, for effective motif identification in high-throughput omics data by learning kernel length from data adaptively. Empirical evaluations on DNA-protein binding and DNase footprinting cases well demonstrated that vConv-based networks have superior performance to their convolutional counterparts regardless of model complexity. Meanwhile, vConv could be readily integrated into multi-layer neural networks as an "in-place replacement" of canonical convolutional layer. All source codes are freely available on GitHub for academic usage.