Estimating the fundamental niche: accounting for the uneven availability of existing climates Studies that question important conceptual and methodological aspects of the field of ecological niche modeling (and species distribution modeling) have cast doubts on whether it is possible to estimate the fundamental niche of a species using presence-only data. The main limitation in niche estimation is that presence data come from the realized niche, which is only a subset of the fundamental niche. Most of the existing methods lack the ability to overcome this limitation and therefore fit niches that resemble the realized niche. To obtain a more accurate estimate of the fundamental niche, we propose using the geographic region that is accessible to a species (based on its dispersal ability) to determine a sampling distribution, in environmental space, from which we can quantify the likelihood of observing a particular environmental combination in a sample of presence points. We incorporate this sampling distribution into a multivariate normal model (i.e., a Mahalanobis distance model) by creating a weight function that modifies the probabilities of observing different environmental combinations in a sample of presences. This modification accounts for the uneven availability of environmental conditions. We show that the parameters of this weighted-normal model can be approximated with a maximum likelihood estimation approach; and then used to draw ellipsoids (confidence regions) that represent the fundamental niche of the species. We illustrate the application of our model with two worked examples. First, we use presence data for an invasive species and an accessible area that includes only its native range to evaluate whether the fitted model predicts confirmed establishments of the species outside its native range. Second, we use presence data for closely related species with known accessible areas to demonstrate how the different dispersal abilities of the species constrain a classic Mahalanobis distance model. Overall, we show that accounting for the distribution of environmental conditions that are accessible to a species indeed affects the estimation of the ellipsoids used to model its fundamental niche.
Hot and cold spots 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
Action at a distance: Defects in division plane positioning in the root meristematic zone affect cell organization in the differentiation zone Cell division plane orientation is critical for plant and animal development and growth. TANGLED1 (TAN1) and AUXIN-INDUCED-IN-ROOT-CULTURES9 (AIR9) are division-site localized microtubule-binding proteins required for division plane positioning. tan1 and air9 Arabidopsis thaliana single mutants have minor or no noticeable phenotypes but the tan1 air9 double mutant has synthetic phenotypes including stunted growth, misoriented divisions, and aberrant cell-file rotation in the root differentiation zone. These data suggest that TAN1 plays a role in nondividing cells. To determine whether TAN1 is required in elongating and differentiating cells in the tan1 air9 double mutant, we limited its expression to actively dividing cells using the G2/M-specific promoter of the syntaxin KNOLLE (pKN:TAN1-YFP). Unexpectedly, in addition to rescuing division plane defects, pKN:TAN1-YFP rescued root growth and the root differentiation zone cell file rotation defects in the tan1 air9 double mutant. This suggests that defects that occur in the meristematic zone later affect the organization of elongating and differentiating cells.

Summary StatementExpression of TAN1 in the root meristematic zone rescues cell file rotation defects in tan1 air9 mutants, suggesting defects that occur in mitosis may influence organization of nondividing cells. The adapted Activity-By-Contact model for enhancer-gene assignment and its application to single-cell data Identifying regulatory regions in the genome is of great interest for understanding the epigenomic landscape in cells. One fundamental challenge in this context is to find the target genes whose expression the regulatory regions affect. A recent successful method is the Activity-By-Contact (ABC) model [10] which scores enhancer-gene interactions based on enhancer activity and the contact frequency of an enhancer to its target gene. It requires two types of assays to measure enhancer activity, which limits the applicability. The ABC score describes regulatory interactions entirely from a genes perspective, and does not account for all the candidate target genes of an enhancer. Moreover, there is no implementation available which would allow for an integration with transcription factor (TF) binding information or an efficient analysis of single-cell data. We show that the ABC-model is comparably accurate with only one assay to measure enhancer activity. Further, we demonstrate that it can yield a higher accuracy by adapting the enhancer activity according to the number of contacts the enhancer has to its candidate target genes. We combined our adapted ABC-model with TF binding information and illustrate an analysis of a single-cell ATAC-seq data set of the human heart. We were able to characterise cell type-specific regulatory interactions as well as to prioritise candidate TFs that drive cell type-specific expression. All executed processing steps are incorporated into our new computational pipeline STARE. The software is available at https://github.com/schulzlab/STARE. 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. 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. Structure-Functional-Selectivity Relationship Studies on A-86929 Analogs and Small biaryl Fragments toward Discovery of biased D1 agonists Dopamine regulates normal functions such as movement, reinforcement learning, and cognition, and its dysfunction has been implicated in multiple psychiatric and neurological disorders. Dopamine acts through the D1- (D1R and D5R) and D2-class (D2R, D3R and D4R) of seven transmembrane receptors, and activates both G-protein- and {beta}-arrestin-dependent signaling pathways, to mediate its physiological effects. Current dopamine receptor-based therapies are used to ameliorate motor deficits in Parkinsons disease, or as antipsychotic medications for schizophrenia. These drugs show efficacy for ameliorating only some symptoms caused by dopamine dysfunction and are plagued by debilitating side-effects. Studies in primates and rodents have shown that shifting the balance of dopamine receptor signaling towards the arrestin pathway can be beneficial for inducing normal movement, while reducing motor side-effects such as dyskinesias, and can be efficacious at enhancing cognitive function compared to balanced agonists. Several structure-activity-relationship (SAR) studies have embarked on discovering {beta}-arrestin-biased dopamine agonists, focused on D2 partial agonists, non-catechol D1 agonists, and mixed D1/D2R dopamine receptor agonists. Here, we describe an SAR study to identify novel D1R {beta}-arrestin biased ligands using A-86929, a high-affinity D1R catechol agonist, as a core scaffold. Previously described and novel analogs of A-86929 were synthesized and screened in vitro for structure-functional-selectivity relationships (SFSR) studies to identify chemical motifs responsible for {beta}-arrestin biased activity at both D1 and D2Rs. Most of the A-86929 analogs screened were G protein biased but none of them were exclusively arrestin-biased. Additionally, various catechol aryl fragments were designed and synthesized. Other compounds surveyed included hydroxyl and chloro analogs of dopamine to test for hydrogen bonding and ionic interactions. Some of these small molecular probes displayed weak bias towards the {beta}-arrestin pathway. Continued in-depth SFSR studies informed by structure determination, molecular modeling, and mutagenesis studies will facilitate discovery of potent and efficacious arrestin-biased dopamine receptor ligands. "Adopt a Microorganism" methodology enhances the discourse richness of high school students in both hybrid learning and remote education The Internet has changed the way teachers and students access information and build knowledge. The recent COVID-19 pandemic has created challenges for both teachers and students and a demand for new methodologies of remote learning. In the life sciences, mixing online content with practical activities represents an even greater challenge. In microbiology, the implementation of an active teaching methodology, the #Adopt project, based on the social network Facebook(R), represents an excellent option for connecting remote education with classroom activities. In 2020, the version applied in high school, "Adopt a Microorganism", was adapted to meet the demands of emergency remote education owing to the suppression of face-to-face activities caused by the pandemic. In the present study, we assessed how the change in methodology impacted the discourse richness of students from high school integrated with technical education in the Business Administration program of the Federal Institute of Sao Paulo, Sorocaba Campus. Three questionnaires related to the groups of microorganisms (Archaea, Bacteria, Virus, Fungi, and Protozoan) were applied. The students responses in the 2019 and 2020 classes were compared concerning content richness and multiplicity of concepts through the application of the Shannon diversity index, an approach that is generally used to assess biodiversity in different environments. The observed results suggest that remote learning provided students with a conceptual basis and richness of content equivalent to that achieved by students subjected to the hybrid teaching model. In conclusion, this study suggests that the #Adopt project methodology increases students discourse richness in microbiology even without face-to-face traditional classes. Higher phage virulence accelerates the evolution of host resistance Parasites and pathogens vary strikingly in their virulence and the resulting selection they impose on their hosts. While the evolution of different virulence levels is well studied, the evolution of host resistance in response to different virulence levels is less understood and as of now mainly based on observations and theoretical predictions with few experimental tests. Increased virulence can increase selection for host resistance evolution if resistance costs are outweighed by the benefits of avoiding infection. To test this, we experimentally evolved the bacterium Vibrio alginolyticus in the presence of two variants of the filamentous phage, VALG{Phi}8, that differ in their virulence. The bacterial host exhibited two alternative defence strategies against future phage infection: (1) super infection exclusion (SIE) whereby phage-infected cells were immune to subsequent infection at a cost of reduced growth, and (2) surface receptor mutations (SRM) in genes encoding the MSHA type-IV pilus providing resistance to infection by preventing phage attachment. While SIE emerged rapidly against both phages, SRM evolved faster against the high virulence compared to the low virulence phage. Using a mathematical model of our system we show that increasing virulence strengthens selection for SRM due to the higher costs of infection suffered by SIE immune hosts. In both the experiments and the model, higher levels of SRM in the host population drove more rapid phage extinction. Thus, by accelerating the evolution of host resistance, more virulent phages caused shorter epidemics. CandyCollect: At-home saliva sampling for capture of respiratory pathogens Streptococcus pyogenes is a major human-specific bacterial pathogen and a common cause of a wide range of symptoms from mild infection such as pharyngitis (commonly called strep throat) to lifethreatening invasive infection. Traditional methods for diagnosis include collecting a sample using a pharyngeal swab, which can cause discomfort and even discourage adults and children from seeking proper testing and treatment in the clinic. Saliva samples are an alternative to pharyngeal swabs. To improve the testing experience for strep throat, we developed a novel lollipop-inspired sampling platform (called CandyCollect) to capture bacteria in saliva. The device can be used in clinics or in the home and shipped back to a lab for analysis, integrating with telemedicine. CandyCollect is designed to capture bacteria on an oxygen plasma treated polystyrene surface embedded with flavoring substances to enhance the experience for children and inform the required time to complete the sampling process. In addition, the open channel structure prevents the tongue from scraping and removing the captured bacteria. Flavoring substances did not affect bacterial capture and the device has a shelf life of at least 2 months (with experiments ongoing to extend the shelf life). We performed a usability study with 17 participants who provided feedback on the device design and the dissolving time of the candy. This technology and advanced processing techniques, including polymerase chain reaction (PCR), will enable user-friendly and effective diagnosis of streptococcal pharyngitis. Global membrane tension is independent of polyacrylamide substrate stiffness Most animal cells are surrounded by a cell membrane and an underlying actomyosin cortex. Both structures are linked, and they are under tension. In-plane membrane tension and cortical tension both influence many cellular processes, including cell migration, division, and endocytosis. However, while actomyosin tension is regulated by substrate stiffness, how membrane tension responds to mechanical substrate properties is currently poorly understood. Here, we probed the effective membrane tension of neurons and fibroblasts cultured on glass and polyacrylamide substrates of varying stiffness using optical tweezers. In contrast to actomyosin-based traction forces, both peak forces and steady state tether forces of cells cultured on hydrogels were independent of substrate stiffness and did not change after blocking myosin II activity using blebbistatin, indicating that tether and traction forces are not directly linked. Peak forces in fibroblasts on hydrogels were about twice as high as those in neurons, indicating stronger membrane-cortex adhesion in fibroblasts. Steady state tether forces were generally higher in cells cultured on hydrogels than on glass, which we explain by a mechanical model. Our results provide new insights into the complex regulation of effective membrane tension and pave the way for a deeper understanding of the biological processes it instructs. DNA Strand Breaks and Gaps Direct Retroviral Integration Integration into a host genome is essential for retrovirus infection and is catalyzed by a nucleoprotein complex (Intasome) containing the viral integrase (IN) and reverse transcribed (RT) copy DNA (cDNA). Previous studies demonstrated DNA site recognition limited intasome integration. Using single molecule Forster resonance energy transfer (smFRET), we show Prototype Foamy Virus (PFV) intasomes pause at DNA strand breaks and gaps. The break/gap discontinuities are similar to base excision repair (BER) lesion-processing intermediates, which affect retrovirus integration in vivo. Pausing targeted site-directed integration at the break/gap without inducing intasome conformational alterations. An 8-oxo-guanine lesion normally processes by BER and a G/T mismatch or a +T nucleotide insertion that induce flexibility or a bend in the DNA backbone did not promote intasome pausing or targeted integration. These results suggest that repair intermediates can modulate dynamic intasome-DNA interactions which target retroviral integration. Changes in Behavior and Salivary Serum Amyloid A in cows with Subclinical Mastitis The aim of this study was to identify detailed changes in behavior, and in salivary serum amyloid A (SAA), associated with subclinical mastitis. This included standard sickness behaviors (which form part of an adaptive response to conserve energy, minimize heat loss and facilitate recovery following infection and inflammation) and lesser-studied behaviors, that are biologically predicted to change at low-level infection, and therefore particularly relevant for detecting subclinical mastitis (here labelled  luxury behaviors). SAA is a physiological marker of systemic inflammation, with levels in milk and serum already known to increase during subclinical mastitis. We investigated whether the same was true for SAA in cow saliva. Data were collected for n = 34 commercial barn-housed dairy cows, tested in 17 matched-pairs. Each pair comprised a cow with subclinical mastitis (SCM) and a healthy control (CTRL), identified using somatic cell count (SCC) (SCM: SCC >200 x1000 cells/ml; CTRL: SCC <100 x1000 cells/ml). SCM cows were selected for study ad-hoc, at which point they were paired with a CTRL cow, based upon parity and calving date; consequently, the full data set was accrued over several months. Data were collected for each pair over 3 days: SCC (l4:00-l5:00h) Day 1; behavior (24h from 00:00h) Day 2; salivary serum amyloid-A (SAA) Day 3. We report, for the first time, that an increase in salivary SAA occurs during subclinical mastitis; SAA was higher in SCM cows and demonstrated a positive (weak) correlation with SCC. The behavioral comparisons revealed that SCM cows had reductions in activity (behavioral transitions and distance moved), social exploration, social reactivity (here: likelihood to move away/be displaced following receipt of agonism), performance of social grooming and head butts, and the receipt of agonistic non-contact challenges. In addition, SCM cows received more head swipes, and spent a greater proportion of time lying with their head on their flank than CTRL cows. SCM cows also displayed an altered feeding pattern; they spent a greater proportion of feeding time in direct contact with two conspecifics, and a lower proportion of feeding time at self-locking feed barriers, than CTRL cows. Behavioral measures were found to correlate, albeit loosely, with serum SAA in a direction consistent with predictions for sickness behavior. These included positive correlations with both lying duration and the receipt of all agonistic behavior, and negative correlations with feeding, drinking, the performance of all social and all agonistic behavior, and social reactivity. We conclude that changes in salivary SAA, social behavior, and activity offer potential in the detection of subclinical mastitis and recommend further investigation to substantiate and refine our findings. The isoforms of pyruvate kinase act as nutrient sensors for the β-cell KATP channel Insulin release from pancreatic {beta}-cells is triggered by the metabolism-dependent closure of ATP-sensitive K+ (KATP) channels. Here, we used isoform-specific deletion of pyruvate kinase (PK) to identify the unique roles of PKm1 and PKm2 in the {beta}-cell nutrient response. In the absence of PK activity, cytosolic ATP/ADP was insufficient to close KATP channels. Both PK isoforms were found on the plasma membrane, where they endow KATP channels with nutrient sensitivity and specificity. Although PKm1 accounts for >90% of cellular activity, PKm2 controls glucosedependent KATP closure. Amino acids, which generate mitochondrial phosphoenolpyruvate (PEP) without allosterically activating PKm2, depend on PKm1 to raise ATP/ADP and close KATP channels. Shifting the {beta}-cell from PKm2 to PKm1 correspondingly increased amino acid sensitivity and boosted secretory output. In addition to these membrane depolarizing "on-switches," {beta}-cell deletion of PCK2 exposed a mitochondrial PEP-dependent "off-switch" that facilitates Ca2+ extrusion. Restricting PCK2 prolonged Ca2+ influx and increased insulin secretion. In summary, it is not the KATP channel alone, but the PK-KATP microcompartment that confers nutrient sensitivity to the {beta}-cell plasma membrane. The differential response of the PK isoforms to the glycolytic and mitochondrial sources of PEP influences {beta}-cell nutrient preference, and controls the oscillatory cycle regulating insulin secretion. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. Natural variations in reproductive timing of brood size and X-chromosome nondisjunction in Caenorhabditis elegans Reproductive timing underlies self-sperm production and X-chromosome nondisjunction rate in Caenorhabditis elegans. These traits can be measured as brood size and male-production rate, and brood size decreases and male-production rate increases as worms age. These phenotypic changes occur simultaneously; however, whether these changes share common underlying genetic architectures still remains unclear. To enhance our understanding of reproductive timing in C. elegans, we measured and compared these two reproductive timing-related traits of 96 wild strains during early, late and total reproductive periods. We found that the two traits exhibited natural phenotypic variation with few outliers and a similar reproductive timing pattern as previous reports. The brood size and the male-production rate were not genetically correlated, implying that they might have different genetic architectures and that meiosis quality can be maintained despite increased progeny-production in the late reproductive period caused by more self-sperm. We also identified loci and candidate genetic variants significantly associated with male-production rate in the late and total reproductive periods. Our results provide an insight into life history traits in wild C. elegans strains. Rapid genotyping of tilapia lake virus (TiLV) using Nanopore sequencing Infectious diseases represent one of the major challenges to sustainable aquaculture production. Rapid, accurate diagnosis and genotyping of emerging pathogens during early-suspected disease cases is critical to facilitate timely response to deploy adequate control measures and prevent or reduce spread. Currently, most laboratories use PCR to amplify partial pathogen genomic regions, occasionally combined with sequencing of PCR amplicon(s) using conventional Sanger sequencing services for confirmatory diagnosis. The main limitation of this approach is the lengthy turnaround time. Here, we report an innovative approach using a previously developed specific PCR assay for pathogen diagnosis combined with a new Oxford Nanopore Technologies (ONT)-based amplicon sequencing method for pathogen genotyping. Using fish clinical samples, we applied this approach for the rapid confirmation of PCR amplicon sequences identity and genotyping of tilapia lake virus (TiLV), a disease-causing virus affecting tilapia aquaculture globally. The consensus sequences obtained after polishing exhibit strikingly high identity to references derived by Illumina and Sanger methods (99.83-100%). This study suggests that ONT-based amplicon sequencing is a promising platform to deploy in regional aquatic animal health diagnostic laboratories in low and medium income countries, for fast identification and genotyping of emerging infectious pathogens from field samples within a single day. Chemotherapy coupled to macrophage inhibition leads to enhanced T and B cell infiltration and durable triple negative breast cancer regression Immunosuppressive elements within the tumor microenvironment such as Tumor Associated Macrophages (TAMs) can present a barrier to successful anti-tumor responses by cytolytic T cells. We employed preclinical syngeneic p53 null mouse models of triple negative breast cancer (TNBC) to develop a treatment regimen that harnessed the immunostimulatory effects of low-dose cyclophosphamide coupled with the pharmacologic inhibition of TAMs using either a small molecule CSF1R inhibitor or an anti-CSF1R antibody. This therapeutic combination was used to successfully treat several highly aggressive TNBC murine mammary tumors and lung metastasis. Using this regimen and single cell RNA sequencing we characterized tumor infiltrating lymphocytes (TILs) including helper T cells and antigen-presenting B cells that were highly enriched in good responders to combination therapy. Using high dimensional imaging techniques, we identified the close spatial localization of B220+ CD86+ activated B cells and CD4+ T cells in tertiary lymphoid structures that were present up to 6 weeks post-treatment in one model that also exhibited long-term tumor regression post-treatment. We also characterized the transcriptional and metabolic heterogeneity of TAMS in these two closely related claudin-low/mesenchymal subtype tumor models with differential treatment responses. A murine TAM signature derived from the T12 model is highly expressed and conserved in human claudin-low breast cancers, and high expression of the T12 signature correlated with reduced overall survival. This T12 tumor TAM signature may help identify human claudin-low breast cancer patients that will benefit from the combination of cyclophosphamide and anti-CSF1R therapy. These studies illustrate the complexity of the tumor immune microenvironment and highlight different immune responses that result from rationale combinations of immunotherapy.

SignificanceA treatment regimen that harnessed the immunostimulatory effects of cyclophosphamide coupled with the inhibition of CSF1R was used to successfully treat several highly aggressive claudin-low TNBC murine mammary tumors and lung metastasis. A dysmorphic mouse model reveals developmental interactions of chondrocranium and dermatocranium The cranial endo- and dermal skeletons, which comprise the vertebrate skull, evolved independently over 470 million years ago and form separately during embryogenesis. In mammals, much of the cartilaginous chondrocranium is transient, undergoing endochondral ossification or disappearing, so its role in skull morphogenesis is not well studied and it remains an enigmatic structure. We provide complete three-dimensional (3D) reconstructions of the laboratory mouse chondrocranium from embryonic day 13.5 through 17.5 using a novel methodology of uncertainty-guided segmentation of phosphotungstic enhanced 3D microcomputed tomography images with sparse annotation. We evaluate the embryonic mouse chondrocranium and dermatocranium in 3D and delineate the effects of a Fgfr2 variant on embryonic chondrocranial cartilages and on the association with forming dermal bones using the Fgfr2cC342Y/+ Crouzon syndrome mouse. We show that the dermatocranium develops outside of and in shapes that conform to the chondrocranium. Results reveal direct effects of the Fgfr2 variant on embryonic cartilage, on chondrocranium morphology, and on the association between chondrocranium and dermatocranium development. Histologically we observe a trend of relatively more chondrocytes, larger chondrocytes, and/or more matrix in the Fgfr2cC342Y/+ embryos at all timepoints before the chondrocranium begins to disintegrate at E16.5. The chondrocrania and forming dermatocrania of Fgfr2cC342Y/+ embryos are relatively large, but a contrasting trend begins at E16.5 and continues into early postnatal (P0 and P2) timepoints, with the skulls of older Fgfr2cC342Y/+ mice reduced in most dimensions compared to Fgfr2c+/+ littermates. Our findings have implications for the study and treatment of human craniofacial disease, for understanding the impact of chondrocranial morphology on skull growth, and potentially on the evolution of skull morphology. Functional buffering via cell-specific gene expression promotes tissue homeostasis and cancer robustness Functional buffering that ensures biological robustness is critical for maintaining tissue homeostasis, organismal survival, and evolution of novelty. However, the mechanism underlying functional buffering, particularly in multicellular organisms, remains largely elusive. Here, we developed an inference index (C-score) for Cell-specific Expression- BUffering (CEBU), whereby functional buffering is mediated via expression of buffering genes in specific cells and tissues in humans. By computing C-scores across 684 human cell lines using genome-wide CRISPR screens and transcriptomic RNA-seq, we report that C- score-identified putative buffering gene pairs are enriched for members of the same duplicated gene family, pathway, and protein complex. Furthermore, CEBU is especially prevalent in tissues of low regenerative capacity (e.g., bone and neuronal tissues) and is weakest in highly regenerative blood cells, linking functional buffering to tissue regeneration. Clinically, the buffering capacity enabled by CEBU can help predict patient survival for multiple cancers. Our results reveal CEBU as a buffering mechanism contributing to tissue homeostasis and cancer robustness in humans.

Summary blurbWe unveil a genome-wide functional buffering mechanism, termed Cell-specific Expression Buffering (CEBU), whereby gene expression contributes to functional buffering in specific cell types and tissues. We link CEBU to genetic interactions, tissue homeostasis and cancer robustness. The Telomere Length Landscape of Prostate Cancer Replicative immortality is a hallmark of cancer, and can be achieved through telomere lengthening and maintenance. We report telomere lengths (TLs) of 392 localized prostate cancer tumours and characterize their relationship to genomic, transcriptomic and proteomic features. Shorter tumour TLs were associated with elevated genomic instability, including single-nucleotide variants, indels and structural variants. Genes involved in cell proliferation and signaling were correlated with tumour TL at all levels of the central dogma. TL was also associated with multiple clinical features of a tumour. Longer TLs in non-tumour samples were associated with a lower rate of biochemical relapse after definitive local therapy. Our analysis integrates multi-omics data to illuminate the relationship of specific genomic alterations in a tumour and TL in prostate cancer. Although the role of telomere length in cancer has been well studied, its association to genomic features is less well known. We describe the multi-level integration of telomere length, genomics, transcriptomics and proteomics in localized prostate cancer. Patient Summary We examined the association between telomere length and multiple omics-level data in prostate cancer. We observed that traditional telomere mutations are rare in prostate cancer and that telomere length is associated with multiple measure of genomic instability. The female pheromone (Z)-4-undecenal mediates flight attraction and courtship in Drosophila melanogaster Specific mate communication and recognition underlies reproduction and hence speciation. Mate communication evolves during adaptation to ecological niches and makes use of social signals and habitat cues.

Our study provides new insights in Drosophila melanogaster premating olfactory communication, showing that female pheromone (Z)-4-undecenal (Z4-11Al) and male pheromone cVA interact with food odour in a sex-specific manner. Furthermore, Z4-11Al, which mediates upwind flight attraction in both sexes, also elicits courtship in experienced males.

Twin variants of the olfactory receptor Or69a are co-expressed in the same olfactory sensory neurons, and feed into the same glomerulus in the antennal lobe. Z4-11Al is perceived via Or69aB, while the food odorant (R)-linalool is a main ligand for the other variant, Or69aA. That Z4-11Al mediates courtship in experienced males, not (R)- linalool, is probably due to courtship learning. Behavioural discrimination is reflected by calcium imaging of the antennal lobe, showing distinct glomerular activation patterns by these two compounds.

Male sex pheromone cVA is known to affect male and female courtship at close range, but does not elicit upwind flight attraction as a single compound, in to contrast to Z4-11Al. A blend of cVA and the food odour vinegar attracted females, while a blend of female pheromone Z4-11Al and vinegar attracted males instead.

Sex-specific upwind flight attraction to blends of food volatiles and male and female pheromone, respectively, adds a new element to Drosophila olfactory premating communication and is an unambiguous paradigm for identifying the behaviourally active components, towards a more complete concept of food-pheromone odour objects.

SummaryThe female-produced, species-specific volatile pheromone of D. melanogaster attracts both sexes from a distance, alone and in concert with food odorants, and elicits courtship in experienced males. Is Vaccination a Viable Method to Control Johne's Disease Caused by Mycobacterium avium subsp. paratuberculosis? Data from 12 Million Ovine Vaccinations and 7.6 Million Carcass Examinations in New South Wales, Australia from 1999-2009 BackgroundMycobacterium avium subsp. paratuberculosis (MAP) causes Johnes disease (or paratuberculosis), a chronic wasting disease of ruminants and other animals resulting from granulomatous enteritis. There are increasing concerns that MAP is zoonotic. The prevalence of Johnes disease is increasing worldwide. In an attempt to control an epidemic of ovine Johnes disease (OJD) in New South Wales (NSW), a government/industry sponsored voluntary vaccination/on-farm management program commenced in 2000. We report herein an observational study of changes in disease prevalence as vaccination progressed, based on abattoir surveillance data for OJD from 1999 to 2009. We also discuss the epidemiological, policy, regulatory, research, economic and sociological elements that contributed to the development of a mature control program, whose aim was to halt the epidemic spread of OJD in a naive sheep population.

MethodsNSW was divided into areas of "High" (HPA), "Medium" (MPA) and "Low" (LPA) OJD prevalence. A killed whole cell vaccine (Gudair(R)) was administered to sheep from 2000 to 2009. Trained examiners evaluated the viscera of adult sheep carcasses at slaughter for gross evidence of OJD. MAP infection was confirmed by histopathology.

Principal findingsFrom 2000-2009, 12 million vaccine doses were administered in NSW (91%; 10.9 million in the HPA). Many of the vaccinated flocks were suffering > 5% annual mortality in adult sheep, with some individual flocks with 10-15% losses attributable to OJD. A total of 7.6 million carcasses were examined (38%; 2.9 million from the HPA). Overall, 16% of slaughter consignments (sheep consigned to the abattoir from a single vendor) were positive for OJD, of which 94% were from the HPA. In the HPA, the percentage of animals with lesions attributable to OJD at slaughter fell progressively from 2.4% (10,406/432,860) at commencement of vaccination in 2000 to 0.8% (1,573/189,564) by 2009. Herd immunity from vaccination in the HPA was estimated at 70% by 2009. This coincided with a progressive decrease in reports of clinical disease and mortalities in vaccinated flocks.

SignificanceWe show a decrease in the prevalence of lesions attributable to OJD in NSW concomitant with initiation of voluntary vaccination, on-farm management plans, abattoir monitoring and feedback of animal prevalence data to sheep producers. We conclude that a target of [&le;] 1% regional prevalence of OJD affected sheep at slaughter is achievable using these interventions. ZIKV disrupts placental ultrastructure and drug transporter expression in mice Congenital Zika virus (ZIKV) infection can induce fetal brain abnormalities. Here, we investigated whether maternal ZIKV infection affects placental physiology and metabolic transport potential and impacts the fetal outcome, regardless of viral presence in the fetus at term. Low (103 PFU-ZIKVPE243; low ZIKV) and high (5x107 PFU-ZIKVPE243; high ZIKV) virus titers were injected into immunocompetent (ICompetent C57BL/6) and immunocompromised (ICompromised A129) mice at gestational day (GD) 12.5 for tissue collection at GD18.5 (term). High ZIKV elicited fetal death rates of 66% and 100%, whereas low ZIKV induced fetal death rates of 0% and 60% in C57BL/6 and A129 dams, respectively. All surviving fetuses exhibited intrauterine growth restriction (IUGR) and decreased placental efficiency. High-ZIKV infection in C57BL/6 and A129 mice resulted in virus detection in maternal spleens and placenta, but only A129 fetuses presented virus RNA in the brain. Nevertheless, pregnancies in both strains produced fetuses with decreased head sizes (p<0.05). Low-ZIKV-A129 dams had higher IL-6 and CXCL1 levels (p<0.05), and their placentas showed increased CCL-2 and CXCL-1 contents (p<0.05). In contrast, low-ZIKV-C57BL/6 dams had an elevated CCL2 serum level and increased type I and II IFN expression in the placenta. Notably, less abundant microvilli and mitochondrial degeneration were evidenced in the placental labyrinth zone (Lz) of ICompromised and high-ZIKV-ICompetent mice but not in low-ZIKV-C57BL/6 mice. In addition, decreased placental expression of the drug transporters P-glycoprotein (P-gp) and breast cancer resistance protein (Bcrp) and the lipid transporter Abca1 was detected in all ZIKV-infected groups, but Bcrp and Abca1 were only reduced in ICompromised and high-ZIKV ICompetent mice. Our data indicate that gestational ZIKV infection triggers specific proinflammatory responses and affects placental turnover and transporter expression in a manner dependent on virus concentration and maternal immune status. Placental damage may impair proper fetal-maternal exchange function and fetal growth/survival, likely contributing to congenital Zika syndrome. 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. Morphogenesis and cell ordering in confined bacterial biofilms Biofilms are aggregates of bacterial cells surrounded by an extracellular matrix. Much progress has been made in studying biofilm growth on solid substrates; however, little is known about the biophysical mechanisms underlying biofilm development in three-dimensional confined environments, in which the biofilm-dwelling cells must push against and even damage the surrounding environment to proliferate. Here, combining single-cell imaging, mutagenesis, and rheological measurement, we reveal the key morphogenesis steps of Vibrio cholerae biofilms embedded in hydrogels as they grow by four orders of magnitude from their initial size. We show that the morphodynamics and cell ordering in embedded biofilms are fundamentally different from those of biofilms on flat surfaces. Treating embedded biofilms as inclusions growing in an elastic medium, we quantitatively show that the stiffness contrast between the biofilm and its environment determines biofilm morphology and internal architecture, selecting between spherical biofilms with no cell ordering and oblate ellipsoidal biofilms with high cell ordering. When embedded in stiff gels, cells self-organize into a bipolar structure that resembles the molecular ordering in nematic liquid crystal droplets. In vitro biomechanical analysis shows that cell ordering arises from stress transmission across the biofilm-environment interface, mediated by specific matrix components. Our imaging technique and theoretical approach are generalizable to other biofilm-forming species, and potentially to biofilms embedded in mucus or host tissues as during infection. Our results open an avenue to understand how confined cell communities grow by means of a compromise between their inherent developmental program and the mechanical constraints imposed by the environment.

Significance StatementBiofilms are microbial cities in which bacterial cells reside in a polymeric matrix. They are commonly found inside soft confining environments such as food matrices and host tissues, against which bacteria must push to proliferate. Here, by combining single-cell live imaging and mechanical characterization, we show that the confining environment determines the dynamics of biofilm shape and internal structure. This self-organized evolution of biofilm architecture is caused by force transmission between the environment and the biofilm, mediated by the extracellular matrix secreted by the cells. Our findings lead to new ways to understand how bacterial communities develop under mechanical constraints, and potentially to new strategies for preventing and controlling biofilm growth in three-dimensional environments.