ConVarT: a search engine for matching human genetic variants with variants from non-human species The availability of genetic variants, together with phenotypic annotations from model organisms, facilitates comparing these variants with equivalent variants in humans. However, existing databases and search tools do not make it easy to scan for equivalent variants, namely "matching variants" (MatchVars) between humans and other organisms. Therefore, we developed an integrated search engine called ConVarT (http://www.convart.org/) for matching variants between humans, mice, and C. elegans. ConVarT incorporates annotations (including phenotypic and pathogenic) into variants, and these previously unexploited phenotypic MatchVars from mice and C. elegans can give clues about the functional consequence of human genetic variants. Our analysis shows that many phenotypic variants in different genes from mice and C. elegans, so far, have no counterparts in humans, and thus, can be useful resources when evaluating a relationship between a new human mutation and a disease.
Improved detection of evolutionary selection highlights potential bias from different sequencing strategies in complex genomic-regions Balancing selection occurs when multiple alleles are kept at elevated frequencies in equilibrium due to opposing evolutionary pressures. A new statistical method was developed to test for selection using efficient Bayesian techniques. Selection signals in three different data sets, generated with variable sequencing technologies, were compared: clinical trios, HLA NGS typed samples, and whole-genome long-read samples. Genome-wide, selection was observed across multiple gene families whose biological functions favor diversification, revealing established targets as well as 45 novel genes under selection. Using high-resolution HLA typing and long-read sequencing data, for the characterization of the MHC, revealed strong selection in expected peptide-binding domains as well as previously understudied intronic and intergenic regions of the MHC. Surprisingly, SIRPA, demonstrated dramatic selection signal, second only to the MHC in most settings. In conclusion, employing novel statistical approaches and improved sequencing technologies is critical to properly analyze complex genomic regions.
Activation of the type 3 secretion system of enteropathogenic E. coli leads to remodeling of its membrane composition and function The cell envelope of gram-negative bacteria is a complex structure, essential for bacterial survival and for resistance to many antibiotics. Channels that cross the bacterial envelope and the host cell membrane form secretion systems that are activated upon attachment to host, enabling bacteria to inject effector molecules into the host cell, required for bacterial-host interaction. The type III secretion system (T3SS) is critical for the virulence of several pathogenic bacteria, including enteropathogenic E. coli (EPEC). The EPEC T3SS activation is associated with repression of carbon storage regulator (CsrA), resulting in gene expression remodeling, which is known to affect EPEC central carbon metabolism and contributes to the adaptation to a cell-adherent lifestyle in a poorly understood manner. We reasoned that the changes in bacterial envelope upon attachment to host and the activation of a secretion system may involve a modification of the lipid composition of bacterial envelope. Accordingly, we performed a lipidomics analysis on mutant strains that simulate T3SS activation. We saw a shift in glycerophospholipid metabolism towards the formation of lysophospholipids, attributed to corresponding upregulation of the phospholipase pldA and the acyltransferase ygiH upon T3SS activation in EPEC. We also detected a shift from menaquinones and ubiquinones to undecaprenyl lipids, concomitant to abnormal synthesis of O-antigen. The remodeling of lipid metabolism is mediated by CsrA and associated with increased bacteria cell size and Zeta potential, and a corresponding alteration in EPEC permeability to vancomycin, increasing the sensitivity of T3SS-activated strains and of adherent wild type EPEC to the antibiotic.

ImportanceThe characterization of EPEC membrane lipid metabolism upon attachment to host is an important step towards a better understanding the shift of EPEC, a notable human pathogen, from a planktonic to adherent life style. It may also apply to other pathogenic bacteria that use this secretion system. We predict that upon attachment to host cells the lipid remodeling upon T3SS activation contributes to bacterial fitness and promotes host colonization, and show that it is associated with increased cell permeability and higher sensitivity to vancomycin.

To the best of our knowledge, this is the first demonstration of a bacterial lipid remodeling due to activation of a secretion system. Nanocrown electrodes for reliable and robust intracellular recording of cardiomyocytes and cardiotoxicity screening Drug-induced cardiotoxicity arises primarily when a compound alters the electrophysiological properties of cardiomyocytes. Features of intracellular action potentials (iAPs) are powerful biomarkers that predict proarrhythmic risks. However, the conventional patch clamp techniques for measuring iAPs are either laborious and low throughput or not suitable for measuring electrically connected cardiomyocytes. In the last decade, a number of vertical nanoelectrodes have been demonstrated to achieve parallel and minimally-invasive iAP recordings. Nanoelectrodes show great promise, but the large variability in success rate, signal strength, and the low throughput of device fabrication have hindered them from being broadly adopted for proarrhythmia drug assessment. In this work, we developed vertically-aligned and semi-hollow nanocrown electrodes that are mechanically robust and made through a scalable fabrication process. Nanocrown electrodes achieve >99% success rates in obtaining intracellular access through electroporation, allowing reliable and simultaneous iAP recordings from up to 57 human pluripotent stem-cell-derived cardiomyocytes (hPSC-CMs). The accuracy of nanocrown electrode recordings is validated by simultaneous patch clamp recording from the same cell. Nanocrown electrodes enable prolonged iAP recording for continual monitoring of the same cells upon the sequential addition of four to five incremental drug doses. In this way, the dose-response data is self-referencing, which avoids the cell-to-cell variations inherent to hPSC-CMs. We are hopeful that this technology development is a step towards establishing an iAP screening assay for preclinical evaluation of drug-induced arrhythmogenicity. Thermal preconditioning in a reef-building coral alleviates oxidative damage through a BI-1 mediated antioxidant response Global coral reef decline is largely driven by the breakdown of the coral-algal symbiosis during temperature stress. Corals can acclimatize to higher temperatures, but the cellular processes underlying this ability are poorly understood. We show that preconditioning-based improvements in thermal tolerance in Pocillopora acuta are accompanied by increases in host glutathione reductase (GR) activity and expression, which prevents DNA damage. A strong correlation between GR and BI-1 expression in heat-stressed preconditioned corals and the presence of an antioxidant response element (ARE) in the GR promoter suggest BI-1 could regulate GR expression through Nrf2/ARE pathway. To fortify this link, we developed an siRNA- mediated gene knockdown protocol and targeted the coral BI-1 gene. BI-1 knock-down specifically decreased GR expression and activity and increased oxidative DNA damage in heat- stressed preconditioned corals, showing that BI-1-mediated enhanced antioxidant response during acute heat stress is a key mechanism that prevents oxidative DNA damage after preconditioning.

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

ImportanceTransition between life cycle stages is a common feature among parasitic pathogens and its regulation must be optimized to balance persistence of infection with transmission. The early transcription factors (TFs) regulating commitment to differentiate are totally unknown in Giardia. In this work, we identified GLP4, a previously uncharacterized GARP-like TF, as an early-acting transcriptional repressor that inhibits G1+S cells from entering the encystation pathway. GLP4 is therefore a key regulator controlling the balance between proliferative growth and terminal differentiation into infective cysts. Hotspots of pest-induced US urban tree death, 2020-2050 O_LIUrban trees are important nature-based solutions for future wellbeing and livability but are at high risk of mortality from insect pests. In the United States (US), 82% of the population live in urban settings and this number is growing, making urban tree mortality a matter of concern for most of its population. Until now, the magnitudes and spatial distributions of risks were unknown.
C_LIO_LIHere, we combine new models of street tree populations in [~]30,000 US communities, species-specific spread predictions for 57 invasive insect species, and estimates of tree death due to insect exposure for 48 host tree genera.
C_LIO_LIWe estimate that 1.4 million street trees will be killed by invasive insects from 2020 through 2050, costing an annualized average of US$ 30M. However, these estimates hide substantial variation: 23% of urban centers will experience 95% of all insect-induced mortality. Further, 90% of all mortality will be due to emerald ash borer (Agrilus planipennis, EAB), which is expected to kill virtually all ash trees (Fraxinus spp.) in >6000 communities.
C_LIO_LIWe define an EAB high-impact zone spanning 902,500km2, largely within the southern and central US, within which we predict the death of 98.8% of all ash trees. "Mortality hotspot cities" include Milwaukee, WI; Chicago, IL; and New York, NY.
C_LIO_LIWe identify Asian wood borers of maple and oak trees as the highest risk future invaders, where a new establishment could cost US$ 4.9B over 30 years.
C_LIO_LIPolicy implications: To plan effective mitigation, managers must know which tree species in which communities will be at the greatest risk, as well as the highest-risk insects. We provide the first country-wide, spatial forecast of urban tree mortality due to invasive insect pests. This framework identifies dominant pest insects and spatial impact hotspots, which can provide the basis for spatial prioritization of spread control efforts such as quarantines and biological control release sites. Our results highlight the need for EAB early-detection efforts as far from current infestations as Seattle, WA. Further, these findings produce a list of biotic and spatiotemporal risk factors for future high-impact US urban forest insect pests.
C_LI A CTP-dependent gating mechanism enables ParB spreading on DNA Proper chromosome segregation is essential in all living organisms. The ParA-ParB-parS system is widely employed for chromosome segregation in bacteria. Previously, we showed that Caulobacter crescentus ParB requires cytidine triphosphate to escape the nucleation site parS and spread by sliding to the neighboring DNA 1. Here, we provide the structural basis for this transition from nucleation to spreading by solving co-crystal structures of a C-terminal domain truncated C. crescentus ParB with parS and with a CTP analog. Nucleating ParB is an open clamp, in which parS is captured at the DNA-binding domain (the DNA-gate). Upon binding CTP, the N-terminal domain (NTD) self-dimerizes to close the NTD-gate of the clamp. The DNA-gate also closes, thus driving parS into a compartment between the DNA-gate and the C-terminal domain. CTP hydrolysis and/or the release of hydrolytic products are likely associated with re-opening of the gates to release DNA and to recycle ParB. Overall, we suggest a CTP-operated gating mechanism that regulates ParB nucleation, spreading, and recycling. The compound of baicalein, wogonin, and oroxylin-A inhibits EMT in A549 cell line via PI3K/AKT-TWIST1-glycolysis pathway Non-small cell lung cancer (NSCLC) is a worldwide disease with high morbidity and mortality, which is most derived from its metastasis. Some researches show that the epithelial-mesenchymal transition (EMT) process promotes lung cancer cells migration and invasion, leading to NSCLC metastasis. Total Flavonoid Aglycones Extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds (baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)), the reconstructed TFAE (reTFAE) present inhibitory effect on the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC. 3CAC: improving the classification of phages and plasmids in metagenomic assemblies using assembly graphs MotivationBacteriophages and plasmids usually coexist with their host bacteria in microbial communities and play important roles in microbial evolution. Accurately identifying sequence contigs as phages, plasmids, and bacterial chromosomes in mixed metagenomic assemblies is critical for further unravelling their functions. Many classification tools have been developed for identifying either phages or plasmids in metagenomic assemblies. However, only two classifiers, PPR-Meta and viralVerify, were proposed to simultaneously identify phages and plasmids in mixed metagenomic assemblies. Due to the very high fraction of chromosome contigs in the assemblies, both tools achieve high precision in the classification of chromosomes but perform poorly in classifying phages and plasmids. Short contigs in these assemblies are often wrongly classified or classified as uncertain.

ResultsHere we present 3CAC, a new three-class classifier that improves the precision of phage and plasmid classification. 3CAC starts with an initial three-class classification generated by existing classifiers and improves the classification of short contigs and contigs with low confidence classification by using proximity in the assembly graph. Evaluation on simulated metagenomes and on real human gut microbiome samples showed that 3CAC outperformed PPR-Meta and viralVerify in both precision and recall, and increased F1-score by 10-60 percentage points.

AvailabilityThe 3CAC software is available on https://github.com/Shamir-Lab/3CAC.

Contactrshamir@tau.ac.il

Supplementary informationSupplementary data are available at Bioinformatics online. 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. Active enhancers strengthen insulation by RNA-mediated CTCF binding at TAD boundaries Vertebrate genomes are partitioned into Topologically Associating Domains (TADs), which are typically bound by head-to-head pairs of CTCF binding sites. Transcription at TAD boundaries correlates with better insulation, however, it is not known whether the boundary transcripts themselves contribute to boundary function. Here we characterize boundary-associated RNAs genome-wide, focusing on the disease-relevant INK4a/ARF TAD. Using CTCF site deletions and boundary-associated RNA knockdowns, we observe that boundary-associated RNAs facilitate recruitment and clustering of CTCF at TAD borders. The resulting CTCF enrichment enhances TAD insulation, enhancer:promoter interactions and TAD gene expression. Importantly, knockdown of boundary-associated RNAs results in loss of boundary insulation function. Using enhancer deletions and CRISPRi of promoters we show that active TAD enhancers but not promoters induce boundary-associated RNA transcription, thus defining a novel class of regulatory enhancer RNAs. Thermal modulation of Zebrafish exploratory statistics reveals constraints on individual behavioral variability BackgroundVariability is a hallmark of animal behavior. It contributes to survival by endowing individuals and populations with the capacity to adapt to ever-changing environmental conditions. Intra-individual variability is thought to reflect both endogenous and exogenous modulations of the neural dynamics of the central nervous system. However, how variability is internally regulated and modulated by external cues remains elusive. Here we address this question by analyzing the statistics of spontaneous exploration of freely swimming zebrafish larvae, and by probing how these locomotor patterns are impacted when changing the water temperatures within an ethologically relevant range.

ResultsWe show that, for this simple animal model, five short-term kinematic parameters - interbout interval, turn amplitude, travelled distance, turn probability and orientational flipping rate - together control the long-term exploratory dynamics. We establish that the bath temperature consistently impacts the means of these parameters, but leave their pairwise covariance unchanged. These results indicate that the temperature merely controls the sampling statistics within a well-defined kinematic space delineated by this robust statistical structure. At a given temperature, individual animals explore the behavioral space over a timescale of tens of minutes, suggestive of a slow internal state modulation that could be externally biased through the bath temperature. By combining these various observations into a minimal stochastic model of navigation, we show that this thermal modulation of locomotor kinematics results in a thermophobic behavior, complementing direct gradient-sensing mechanisms.

ConclusionsThis study establishes the existence of a well-defined locomotor space accessible to zebrafish larvae during spontaneous exploration, and quantifies self-generated modulation of locomotor patterns. Intra-individual variability reflects a slow diffusive-like probing of this space by the animal. The bath temperature in turn restricts the sampling statistics to sub-regions, endowing the animal with basic thermophobicity. This study suggests that in Zebrafish, as well as in other ectothermic animals, ambient temperature could be used to efficiently manipulate internal states in a simple and ethological way. Mechanistic basis for understanding the dual activities of the bifunctional Azotobacter vinelandii mannuronan C-5 epimerase and alginate lyase AlgE7 The functional properties of alginates are dictated by the monomer composition and molecular weight distribution. Mannuronan C-5 epimerases determine the former by epimerizing {beta}-D-mannuronic acid residues (M) into -L-guluronic acid residues (G). The molecular weight is affected by alginate lyases, which cleave alginate chains through {beta}-elimination. The reaction mechanisms for the epimerization and cleavage are similar and some enzymes can perform both. These dualistic enzymes share high sequence identity with mannuronan C-5 epimerases without lyase activity, and the mechanism behind their activity as well as the amino acids responsible for it are still unknown. In this study, we investigate mechanistic determinants of the bifunctional epimerase and lyase activity of AlgE7 from Azotobacter vinelandii. Based on sequence analyses, a range of AlgE7 variants were constructed and subjected to activity assays and product characterization by NMR. Our results show that the lyase activity of AlgE7 is regulated by the type of ion present: Calcium promotes it, whereas NaCl reduces it. By using defined poly-M and poly-MG substrates, the preferred cleavage sites of AlgE7 were found to be M{downarrow}XM and G{downarrow}XM, where X can be either M or G. By studying AlgE7 mutants, R148 was identified as an important residue for the lyase activity, and the point mutant R148G resulted in an enzyme with only epimerase activity. Based on the results obtained in the present study we suggest a unified catalytic reaction mechanism for both epimerase and lyase activities where H154 functions as the catalytic base and Y149 as the catalytic acid.

ImportancePost-harvest valorisation and upgrading of algal constituents is a promising strategy in the development of a sustainable bioeconomy based on algal biomass. In this respect, alginate epimerases and lyases are valuable enzymes for tailoring of the functional properties of alginate, a polysaccharide extracted from brown seaweed with numerous applications in food, medicine and material industries. By providing a better understanding of the reaction mechanism and of how the two enzyme reactions can be altered by changes in reaction conditions, this study opens for further applications of bacterial epimerases and lyases in enzymatic tailoring of alginate polymers. Robust retrieval of data stored in DNA by de Bruijn graph-based de novo strand assembly DNA data storage is a rapidly developing technology with great potential due to its high density, long-term durability, and low maintenance cost. The major technical challenges include various errors, such as the strand breaks, rearrangements, and indels that frequently arise during DNA synthesis, amplification, sequencing, and preservation. Through a de novo assembly strategy, we developed an algorithm based on the de Bruijn graph and greedy path search (DBGPS) to address these issues. DBGPS shows distinct advantages in handling DNA breaks, rearrangements, and indels. The robustness of DBGPS is demonstrated by accelerated aging, multiple independent data retrievals, deep error-prone PCR, and large data scale simulations. Remarkably, 6.8 MB of data can be retrieved accurately from a seriously corrupted sample that has been treated at 70 {degrees}C for 70 days. With DBGPS, we were able to achieve a logical density of 1.30 bits/cycle and a physical density of 295 PB/g.

One-Sentence SummaryA de Bruijn graph-based de novo assembly algorithm for DNA data storage enables fast and robust data readouts even with DNA samples that have been severely corrupted. Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. Small extracellular vesicles from young mice prevent frailty, improve healthspan and decrease epigenetic age in old mice. Aging is associated with an increased risk of frailty, disability, comorbidities, institutionalization, falls, fractures, hospitalization, and mortality. Searching for strategies to delay the degenerative changes associated with aging and frailty is interesting. We treated old animals intravenously with small extracellular vesicles (sEVs) derived from adipose mesenchymal stem cells (ADSCs) of young animals, and we found an improvement of several functional parameters usually altered with aging, such as motor coordination, grip strength, fatigue resistance, fur regeneration, and renal function. Frailty index analysis showed that 40% of old control mice were frail, whereas none of the old ADSCs-sEVs treated mice were. Molecular and structural benefits in muscle and kidney accompanied this functional improvement. ADSCs-sEVs induced pro-regenerative effects and a decrease in oxidative stress, inflammation, and senescence markers. Moreover, predicted epigenetic age was lower in tissues of old mice treated with ADSCs-sEVs and their metabolome changed to a youth-like pattern. Finally, we gained some insight into the miRNAs contained in sEVs that might be, at least in part, responsible for the effects observed. We propose that young sEVs treatment can be beneficial against frailty and therefore can promote healthy aging. Differential impact of mutations in cell cycle checkpoint kinases on breast cancer presentation and prognosis. Cell cycle checkpoint kinases ATM/Chk2 and ATR/Chk1 are important tumor suppressors that are activated by DNA damage and other stressors to arrest replication. Germline variants in ATM and CHEK2 associate with increased incidence of estrogen receptor positive (ER+)/HER2- breast cancer and some evidence suggests that somatic inactivation of ATM/Chk2 induces therapy resistance in ER+/HER2- breast cancer patients. However, the relative contributions of germline and somatic inactivation of each of the four cell cycle checkpoint kinases to ER+/HER2- and ER-/PR-/HER2- (or triple negative breast cancer; TNBC) breast cancer incidence and outcomes has not yet been systematically assessed. Here, we comprehensively analyze associations between mutations in each cell cycle checkpoint kinase gene with age at diagnosis, tumor characteristics, and patient survival outcomes using six breast cancer datasets. Overall, we find that of all the cell cycle checkpoint kinases, CHEK2 (p=0.001) is the only gene that is preferentially mutated in ER+/HER2- breast cancer samples relative to TNBCs. CHEK2 mutations also associate with younger age at diagnosis of metastatic ER+/HER2- breast cancer (median age 41.5 years compared to 50 years for CHEK2 wildtype patients, p=0.03) that is highly PR positive (10% increase in PR positivity in metastatic vs primary CHEK2 mutant tumors, as opposed to a 25% decrease in PR positivity in CHEK2 wildtype tumors, p=0.03). These associations appear to be driven by germline variants, rather than somatic mutations. Whether germline or somatic, deleterious (truncating, frameshift or splice site) mutations in CHEK2 associate with rapid disease progression in response to standard first line endocrine therapy when compared to tumors with wildtype CHEK2 in two independent datasets (HR=6.15, p=0.01). Conversely, neither ATM nor ATR mutations are enriched in ER+/HER2- breast cancers relative to TNBC. However, co-incident mutation of ATR and TP53 is 2-fold enriched in this subset (p=0.002), and further enriched in metastatic (3.3%) vs primary (1%) ER+/HER2- breast cancers (p=0.002). This co-incident mutation of ATR and TP53 associates with decreased PR positivity in primary ER+/HER2- breast cancer diagnoses (p<0.0001) and significantly worse metastasis-free survival than mutations in either ATR or TP53 alone in two independent datasets (HR=2.01, p=0.007). This association is not specific to ER+/HER2- breast cancer but is observable in ER- breast cancer patients as well. These results provide the first systematic analysis of the contribution of germline and somatic mutations in cell cycle checkpoint kinases to tumor characteristics affecting patient prognosis and treatment response. They demonstrate a unique role for CHK2 in ER signaling and response to endocrine therapy, and a more universal role for ATR in combination with TP53 in metastatic progression across breast cancer subtypes. These findings provide rationale for streamlined testing of checkpoint kinase mutations as prognostic markers of metastatic potential and/or predictive markers of responsiveness to endocrine treatment in breast cancer patients. Conserved organ-specific microbial assemblages in different populations of a terrestrial crab Brachyuran crabs originated in the oceans and evolved specific morphological and physiological adaptations to live in freshwater, intertidal and even terrestrial habitats but the role of a selection mechanism involving symbiotic microorganisms long these colonization processes are not known. In this work we investigated the associated microbiota of three populations of a terrestrial brachyuran crab, Chiromantes haematocheir, to find evidence of a conserved crab-specific microbiome unrelated to the population of origin and dissimilar from environmental microbial assemblages. Bacterial 16S rRNA gene and fungal ITS sequences were obtained from selected crab organs and environmental matrices to profile microbial communities. In spite of the presence of truly marine larval stages and no gregarious behaviour, favouring microbiota exchanges, we found common, organ-specific microbiota, associated to the gut and the gills of the crabs (with more than 15% of the genera detected specifically enriched only in one organ). Our results suggest an early establishment of a new common, stable microbiota in the transition from water to land. Wrinkle force microscopy: a new machine learning based approach to predict cell mechanics from images Combining experiments with artificial intelligence algorithms, we propose a new machine learning based approach to extract the cellular force distributions from the microscope images. The full process can be divided into three steps. First, we culture the cells on a special substrate allowing to measure both the cellular traction force on the substrate and the corresponding substrate wrinkles simultaneously. The cellular forces are obtained using the traction force microscopy (TFM), at the same time that cell-generated contractile forces wrinkle their underlying substrate. Second, the wrinkle positions are extracted from the microscope images. Third, we train the machine learning system with GAN (generative adversarial network) by using sets of corresponding two images, the traction field and the input images (raw microscope images or extracted wrinkle images), as the training data. The network understands the way to convert the input images of the substrate wrinkles to the traction distribution from the training. After sufficient training, the network is utilized to predict the cellular forces just from the input images. Our system provides a powerful tool to evaluate the cellular forces efficiently because the forces can be predicted just by observing the cells under the microscope, which is a way simpler method compared to the TFM experiment. Additionally, the machine learning based approach presented here has the profound potential for being applied to diverse cellular assays for studying mechanobiology of cells.

Significance StatementCell-generated forces are indispensable determinants of fundamental cell functions such as motility and cell division. As such, quantifying how the forces change upon perturbations to the cells such as gene mutations and drug administration is of profound importance. Here we present a novel machine learning based system that allows for efficient estimations of the forces that are determined only by "observing" microscope images. Given that the cellular traction forces are regulated downstream of diverse signaling pathways, our system - that helps significantly improve the throughput of the measurements - presents a new, high throughput platform for real time analysis of the effects of a massive number of genetic and molecular perturbations on the forces and resulting cell mechanics. Sleep spindle frequency: overnight dynamics, afternoon nap effects, and possible circadian modulation Homeostatic and circadian processes play a pivotal role in determining sleep structure, timing and quality. In sharp contrast with the wide accessibility of the EEG index of sleep homeostasis, an electrophysological measure of the circadian modulation of sleep is still non-available. Evidence suggests that sleep spindle frequencies decelerate during biological night. In order to test the feasibility of measuring this marker in common polysomnographic protocols, the Budapest-Munich database of sleep records (N = 251 healthy subjects, 122 females, age range: 4-69 years), as well as an afternoon nap sleep record database (N = 112 healthy subjects, 30 females, age range: 18-30 years) were analysed by the Individual Adjustment Method of sleep spindle analysis. Slow and fast sleep spindle frequencies were characterized by U-shaped overnight dynamics, with highest values in the first and the fourth-to-fifth sleep cycle and the lowest values in the middle of the sleeping period (cycles 2-3). Age-related attenuation of sleep spindle deceleration was evident. Estimated phases of the nadirs in sleep spindle frequencies were advanced in children as compared to other age groups. Additionally, nap sleep spindles were faster than night sleep spindles (0.57 and 0.39 Hz difference for slow and fast types, respectively). The fine frequency resolution analysis of sleep spindles is a feasible method of measuring the assumed circadian modulation of sleep. Moreover, age-related attenuation of circadian sleep modulation might be measurable by assessing the overnight dynamics in sleep spindle frequency. Phase of the minimal sleep spindle frequency is a putative biomarker of chronotype. 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. Secretory cells are the primary source of pIgR in small airways BackgroundLoss of secretory immunoglobulin A (SIgA) is common in COPD small airways and likely contributes to disease progression. We hypothesized loss of SIgA results from reduced expression of pIgR, a chaperone protein needed for SIgA transcytosis, in the COPD small airway epithelium.

MethodspIgR-expressing cells were defined and quantified at single-cell resolution in human airways using RNA in-situ hybridization, immunostaining, and single-cell RNA sequencing. Complementary studies in mice utilized immunostaining, primary murine tracheal epithelial cell (MTEC) culture, and transgenic mice with secretory or ciliated cell-specific knockout of pIgR. SIgA degradation by human neutrophil elastase or secreted bacterial proteases from non-typeable Haemophilus influenzae (NTHi) was evaluated in vitro.

ResultsWe found that secretory cells are the predominant cell type responsible for pIgR expression in human and murine airways. Loss of SIgA in small airways was not associated with a reduction in secretory cells but rather a reduction in pIgR protein expression despite intact PIGR mRNA expression. Neutrophil elastase and NTHi-secreted proteases are both capable of degrading SIgA in vitro and may also contribute to a deficient SIgA immunobarrier in COPD.

InterpretationLoss of the SIgA immunobarrier in small airways of patients severe COPD is complex and likely results from both pIgR-dependent defects in IgA transcytosis and SIgA degradation.

Key MessagesWhat is the key question? Localized SIgA deficiency in small airways is an established driver of COPD pathogenesis, but the mechanism of loss remains unclear. We hypothesized loss of SIgA is due to reduced numbers of pIgR-expressing cells in SIgA-deficient small airways.

What is the bottom line? pIgR is primarily expressed by secretory cells in human and murine airways. Although numbers of secretory cells are similar between SIgA-deficient and SIgA-replete airways in COPD, there is reduced expression of pIgR protein, but not mRNA, in SIgA-deficient airways. Additionally, host and bacterial proteases degrade SIgA in vitro, suggesting loss of SIgA may relate to both impaired transcytosis and increased degradation.

Why read on? This study highlights the complexity of SIgA immunobarrier maintenance and suggests that strategies aimed at restoring the SIgA immunobarrier will need to account for both impaired transcytosis and degradation by host and/or bacterial proteases. Feline tooth resorption findings from a questionnaire survey on 8115 Finnish cats and a case-control study for a more defined population drawn from the survey Tooth resorption (TR) is one of the most common dental diseases in cats. It is a painful disease characterised by progressive dental destruction, which eventually results in loss of teeth. The aetiology of the disease remains unclear, but associations with old age, breed, other oral and dental diseases, and certain environmental factors have been suspected. We collected health and environmental information of 8115 Finnish cats from 41 breeds through an online survey for cat owners to study TR in Finland with an observational study. The frequency of TR that the owner reported as veterinarian-diagnosed was 3.9% in the entire data (316/8115) and 15% in cats that the owner reported as having veterinarian-diagnosed oral or dental disease (316/2070). We determined the risk factors and breed variation of feline TR by logistic regression with the entire data and a case-control study in a more defined population drawn from the entire data. In the entire data, we defined those cats that the owners had reported as being diagnosed with TR by a veterinarian as TR-cats and compared them to all other cats. In the case-control study, we included only cats with owner-reported veterinarian-diagnosed oral and dental diseases. We defined cats with TR-diagnosis made under sedation as TR-cases. Other veterinarian-diagnosed TR-cats were excluded. Cats with other oral and dental diseases were used as controls. The model with the entire population included 308 TR-cats compared with 7508 cats, and the case-control study included 198 TR-cats and 1683 controls. Although the diagnoses were reportedly made by a veterinarian, the owner answered the survey and might have understood, remembered, or interpreted the diagnoses in a way that could affect the results. The results indicated an increasing risk of TR with age. Dental calculus, gingivitis, and periodontitis were potential risk factors for TR. These findings and the interaction between dental calculus and gingivitis suggest that inflammatory changes caused by dental calculus increase the risk of TR. Keeping food available constantly was a potential protective factor. We found that Cornish Rex, European, and Ragdoll breeds were at higher risk for TR. The observed differences between breeds highlight a genetic contribution. The RAG1 Ubiquitin Ligase Domain Enhances Assembly and Selection of T Cell Receptor Genes to Restrain the Autoimmune Hazard of Generating T Cell Receptor Diversity RAG1/RAG2 (RAG) endonuclease-mediated assembly of diverse lymphocyte antigen receptor genes by V(D)J recombination is critical for the development and immune function of T and B cells. However, this process creates highly self-reactive cells that must be properly selected to suppress autoimmunity. The RAG1 protein contains a ubiquitin ligase domain that stabilizes RAG1 and stimulates RAG endonuclease activity in vitro. We report that mice with a mutation that inactivates the RAG1 ubiquitin ligase in vitro exhibit modestly reduced thymic cellularity, decreased assembly and altered repertoires of T cell receptor (TCR) {beta} and  genes in thymocytes, and impaired thymocyte developmental transitions that require the assembly of TCR{beta} or  genes and signaling by their proteins. These RAG1 mutant mice also exhibit less efficient positive selection and superantigen-mediated negative selection of conventional {beta} T cells, 2) impaired differentiation of iNKT lineage {beta} T cells, and 3) CD4+ {beta} T cells with elevated autoimmune potential. Our findings demonstrate that the RAG1 ubiquitin ligase domain functions in vivo to stimulate the assembly and selection of TCR{beta} and TCR genes, thereby establishing replete diversity of {beta} TCRs and {beta} T cell lineages while restraining the inherent autoimmune hazard of generating diverse antigen specificities. Identifying and testing marker-trait associations for growth and phenology in three pine species: implications for genomic prediction In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 SNPs from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments.