Multi-omic meta-analysis of transcriptomes and the bibliome uncovers novel hypoxia-inducible genes Hypoxia is a condition in which cells, tissues, or organisms are deprived of sufficient oxygen supply. Aerobic organisms have a hypoxic response system, represented by hypoxia-inducible factor 1- (HIF1A), to adapt to this condition. Due to publication bias, there has been little focus on genes other than well-known signature hypoxia-inducible genes. Therefore, in this study, we performed a meta-analysis to identify novel hypoxia-inducible genes. We searched publicly available transcriptome databases to obtain hypoxia-related experimental data, retrieved the metadata, and manually curated it. We selected the genes that are differentially expressed by hypoxic stimulation, and evaluated their relevance in hypoxia by performing enrichment analyses. Next, we performed a bibliometric analysis using gene2pubmed data to examine genes that have not been well studied in relation to hypoxia. Gene2pubmed data that provides information about the relationship between gene and publication. We calculated and evaluated the number of reports and similarity coefficients of each gene to HIF1A, which is a representative gene in hypoxia studies. In this data-driven study, we report that several genes that were not known to be associated with hypoxia, including the G protein-coupled receptor 146 gene, are upregulated by hypoxic stimulation.
MEDICC2: whole-genome doubling aware copy-number phylogenies for cancer evolution Chromosomal instability (CIN) and somatic copy-number alterations (SCNA) play a key role in the evolutionary process that shapes cancer genomes. SC-NAs comprise many classes of clinically relevant events, such as localised amplifications, gains, losses, loss-of-heterozygosity (LOH) events, and recently discovered parallel evolutionary events revealed by multi-sample phasing. These events frequently appear jointly with whole genome doubling (WGD), a transformative event in tumour evolution involving tetraploidization of genomes preceded or followed by individual chromosomal copy-number changes and associated with an overall increase in structural CIN.

While SCNAs have been leveraged for phylogeny reconstruction in the past, existing methods do not take WGD events into account and cannot model parallel evolution. They frequently make use of the infinite sites assumption, do not model horizontal dependencies between adjacent genomic loci and can not infer ancestral genomes. Here we present MEDICC2, a new phylogeny inference algorithm for allele-specific SCNA data that addresses these shortcomings. MEDICC2 dispenses with the infinite sites assumption, models parallel evolution and accurately identifies clonal and subclonal WGD events. It times SCNAs relative to each other, quantifies SCNA burden in single-sample studies and infers phylogenetic trees and ancestral genomes in multi-sample or single-cell sequencing scenarios with thousands of cells.

We demonstrate MEDICC2s ability on simulated data, real-world data of 2,778 single sample tumours from the Pan-cancer analysis of whole genomes (PCAWG), 10 bulk multi-region prostate cancer patients and two recent single-cell datasets of triple-negative breast cancer comprising several thousands of single cells.
A General Biphasic-Bodyweight Model for Scaling Basal Metabolic Rate, Glomerular Filtration Rate and Drug Clearance from Birth to Adulthood Understanding the maturation process of human physiology and metabolism has broad medical and pharmaceutical implications. Age and bodyweight are frequently considered as separate variables in modeling the dynamical changes of human organ functions and of drug clearance from birth to adulthood. The objective of this study is to propose a unified, continuous and bodyweight-only equation to quantify the changes of human basal metabolic rate (BMR), glomerular filtration rate (GFR) and drug clearance (CL) from infancy to adulthood. The BMR datasets were retrieved from a comprehensive historical database of male and female subjects (0.02 to 64 years). The CL datasets for 17 drugs and the GFR dataset were generated from age-incorporated maturation-and-growth models with reported parameter values. The model used in the simulation is independent of the proposed model. A statistical approach was used to simulate the model generated CL and GFR data for a hypothetical population with 26 age groups (ranging from 0 to 20 years). Besides, individual CL data for one drug, and sparse PBPK-modeled CL values for two drugs were also included for analysis. A 4-parameter, mixed-allometry equation with two power-law functions of bodyweight was proposed and evaluated as a general model using nonlinear regression and dimensionless analysis. All datasets universally reveal biphasic curves with two distinct linear segments on log-log plots. Compared with simple allometry, the biphasic model fits satisfactorily to all datasets (based on Akaikes Information Criterion and residual plots). The biphasic equation consists of two reciprocal allometric terms that asymptotically determine the overall curvature. The fitting results show a superlinear scaling phase (slope >1; ca. 1.5 - 3.5) below the characteristic bodyweight at the phase transition; and above which, a sublinear scaling phase (slope <1; ca. 0.5 - 0.7) is evident. The phase-transition bodyweight is ranging from 5 to 20 kg (corresponding to 0.5 - 9 years) and the mean value is around 10 kg ([~]2 years) for all data sets. The dimensionless analysis generalizes, and offers quantitative realization of, the maturation and growth process. In conclusion, the proposed mixed-allometry equation is a generic model that quantitatively describes the phase transition occurring in the human maturation process of BMR, GFR and drug CL. Role of posterodorsal medial amygdala urocortin-3 in pubertal timing in female mice Post-traumatic stress disorder impedes pubertal development and disrupts pulsatile LH secretion in humans and rodents. The posterodorsal sub-nucleus of the medial amygdala (MePD) is an upstream modulator of the hypothalamic gonadotropin-releasing hormone (GnRH) pulse generator, pubertal timing, as well as emotional processing and anxiety. Psychosocial stress exposure alters neuronal activity within the MePD increasing the expression of Urocortin3 (Ucn3) and its receptor corticotropin-releasing factor type-2 receptor (CRFR2) while enhancing the inhibitory output from the MePD to key hypothalamic reproductive centres. We test the hypothesis that psychosocial stress, processed by the MePD, is relayed to the hypothalamic GnRH pulse generator to delay puberty in female mice. We exposed C57Bl6/J female mice to the predator odor, 2,4,5-Trimethylthiazole (TMT), during pubertal transition and examined the effect on pubertal timing, pre-pubertal LH pulses and anxiety-like behaviour. Subsequently, we virally infected Ucn3-cre-tdTomato female mice with stimulatory DREADDs targeting MePD Ucn3 neurons and determined the effect on pubertal timing and pre-pubertal LH pulse frequency. Exposure to TMT during pubertal development delayed puberty, suppressed pre-pubertal LH pulsatility and enhanced anxiety-like behaviour, while activation of MePD Ucn3 neurons reduced LH pulse frequency and delayed puberty. Early psychosocial stress exposure decreases GnRH pulse generator frequency delaying puberty while inducing anxiety-behaviour in female mice, an effect potentially involving Ucn3 neurons in the MePD. 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. Learning and motor inhibitory control in crows and domestic chickens Cognitive abilities allow animals to navigate through complex, fluctuating environments. In the present study, we tested the performance of a captive group of eight crows, Corvus corone, and ten domestic chickens, Gallus gallus domesticus, in the cylinder task, as a test of motor inhibitory control, and reversal learning as a measure of learning ability and behavioural flexibility. Four crows and nine chickens completed the cylinder task, eight crows and six chickens completed the reversal learning experiment. Crows performed better in the cylinder task compared to chickens. In the reversal learning experiment, species did not significantly differ in the number of trials until learning criterion was reached. Performance in the reversal learning experiment did not correlate with performance in the cylinder task in chickens. Our results suggest crows to possess better motor inhibitory control compared to chickens. In contrast, learning performance in a reversal learning task did not differ between the species, indicating similar levels of behavioural flexibility. Interestingly, we describe notable individual differences in performance. We stress the importance not only to compare cognitive performance between species but also between individuals of the same species when investigating the evolution of cognitive skills. Synergy of chemotherapy and macrophage depletion leads to B cell-mediated T cell memory activation 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. Evaluation of US state pollinator plans using 3 evidence-based policymaking frameworks After the US federal government created a national pollinator protection plan in 2015, many states followed with their own. Since their goal is to promote pollinating insect conservation, we wanted to know whether the state plans are using best practices for evidence-based science policy. In early 2019 we found and downloaded every existing, publicly available US state pollinator protection plan. We then used content analysis to assess the goals, scope, and implementation of state-level pollinator protection plans across the US. This analysis was conducted using three distinct frameworks for evidence-based policymaking: US Department of Interior Adaptive Resources Management (ARM), US Environmental Protection Agency management pollinator protection plan (MP3) guidance, and Pew Trusts Pew-MacAthur Results First Project elements of evidence-based state policymaking (PEW) framework. Then we scored them using the framework criteria, to assess whether the plans were using known best practices for evidence based policymaking. Of the 31 states with a state pollinator plan, Connecticut was the state with the lowest total score across the three evaluation frameworks. The state with the highest overall scores, across the three frameworks, was Missouri. Most states did not score highly on the majority of the frameworks. Overall, many state plans were lacking policy elements that address monitoring, evaluation, and adjustment. These missing elements impact states ability to achieve their conservation goals. Our results indicate that states can improve their pollinator conservation policies to better match evidence-based science policy guidance, regardless of which framework is used. Micro topographical instruction of bacterial attachment, biofilm formation and in vivo host response Bio-instructive materials that prevent bacterial biofilm formation and drive an appropriate host immune response have the potential to significantly reduce the burden of medical device-associated infections. As bacterial surface attachment is sensitive to surface topography, we experimentally surveyed 2,176 combinatorially generated shapes using an unbiased high-throughput micro topographical polystyrene polymer chip screen. This identified topographies that reduced colonization in vitro by up to 15-fold compared with a flat surface for both motile and non-motile bacterial pathogens. Equivalent reductions were achieved on polyurethane, a polymer commonly used in medical devices. Using machine learning methods, a set of design rules based on generalisable topo-descriptors was established for predicting bacteria-resistant micro topographies. In a murine foreign body infection model, anti-attachment topographies were shown to be refractory to P. aeruginosa colonization and serendipitously, the fibrotic response to the implant was reduced, offering control of foreign body response using simple topographical patterning of non-eluting implants. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. Structure-function analysis of the nsp14 N7-guanine methyltransferase reveals an essential role in Betacoronavirus replication As coronaviruses (CoVs) replicate in the host cell cytoplasm, they rely on their own capping machinery to ensure the efficient translation of their mRNAs, protect them from degradation by cellular 5 exoribonucleases, and escape innate immune sensing. The CoV nonstructural protein 14 (nsp14) is a bi-functional replicase subunit harboring an N-terminal 3'-to-5' exoribonuclease (ExoN) domain and a C-terminal (N7-guanine)-methyltransferase (N7-MTase) domain that is presumably involved in viral mRNA capping. Here, we aimed to integrate structural, biochemical, and virological data to assess the importance of conserved N7-MTase residues for nsp14s enzymatic activities and virus viability. We revisited the crystal structure of severe acute respiratory syndrome (SARS)-CoV nsp14 to perform an in silico comparative analysis between betacoronaviruses. We identified several residues likely involved in the formation of the N7-MTase catalytic pocket, which presents a fold distinct from the Rossmann fold observed in most known MTases. Next, for SARS-CoV and Middle East respiratory syndrome-CoV, site-directed mutagenesis of selected residues was used to assess their importance for in vitro enzymatic activity. Most of the engineered mutations abolished N7-MTase activity, while not affecting nsp14-ExoN activity. Upon reverse engineering of these mutations into different betacoronavirus genomes, we identified two substitutions (R310A and F426A in SARS-CoV nsp14) abrogating virus viability and one mutation (H424A) yielding a crippled phenotype across all viruses tested. Our results identify the N7-MTase as a critical enzyme for betacoronavirus replication and define key residues of its catalytic pocket that can be targeted to design inhibitors with a potential pan-coronaviral activity spectrum.

Significance StatementThe ongoing SARS-CoV-2 pandemic emphasizes the urgent need to develop efficient broad-spectrum anti-CoV drugs. The structure-function characterization of conserved CoV replicative enzymes is key to identifying the most suitable drug targets. Using a multidisciplinary comparative approach and different betacoronaviruses, we characterized the key conserved residues of the nsp14 (N7-guanine)-methyltransferase, a poorly defined subunit of the CoV mRNA-synthesizing machinery. Our study highlights the unique structural features of this enzyme and establishes its essential role in betacoronavirus replication, while identifying two residues that are critical for the replication of the four betacoronaviruses tested, including SARS-CoV-2. Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector. Artificial seawater based long-term culture of colonial ascidians Tunicates are highly diverse marine invertebrate filter-feeders that are vertebrates closest relatives. These organisms, despite a drastically different body plan during their adulthood, have a tissue complexity related to that of vertebrates. Ascidians, which compose most of the Tunicata, are benthic sessile hermaphrodites that reproduce sexually through a motile tadpole larval stage. Over half of the known ascidians species are able to reproduce asexually by budding, typically leading to the formation of colonies where animals, called zooids, are interconnected through an external vascular system. In addition, colonial ascidians are established models for important biological processes including allorecognition, immunobiology, aging, angiogenesis and whole-body regeneration. However, the current paucity in breeding infrastructures limits the study of these animals to coastal regions.

To promote a wider scientific spreading and popularity of colonial ascidians, we have developed a flexible recirculating husbandry setup for their long-term in-lab culture. Our system is inspired both by the flow-through aquariums used by coastal ascidian labs, as well as by the recirculating in-lab systems used for zebrafish research. Our hybrid system thus combines colony breeding, water filtering and food culturing in a semi-automated system where specimens develop on hanging microscopy glass slides. Temperature, light/dark cycles, flow speed and feeding rates can be controlled independently in four different breeding environments to provide room for species-specific optimization as well as for running experiments. This setup is complemented with a quarantine for the acclimatization of wild isolates.

Herein we present our success in breeding Botrylloides diegensis, a species of colonial ascidians, for more than 3 years in recirculating artificial seawater over 600 km away from their natural habitat. We show that colonies adapt well to in-lab culturing provided that a suitable marine microbiome is present, and that a specific strain can be isolated, propagated and efficiently used for research over prolonged periods of time. The flexible and modular structure of our system can be scaled and adapted to the needs of specific species, such as Botryllus schlosseri, as well as of particular laboratory spaces. Overall, we show that Botrylloides diegensis can be proficiently bred in-land and suggest that our results can be extended to other species of colonial ascidians to promote research on these fascinating animals.

HighlightsO_LIFirst artificial seawater based recirculating aquaculture for colonial ascidians
C_LIO_LIOver 3 years of continuous breeding
C_LIO_LISemi-automated setup with minimized maintenance
C_LIO_LIGood biomass production for strain propagation
C_LIO_LI4 different culture conditions for optimized breeding for species of interest
C_LI Fecal immunoglobulin A (IgA) and its subclasses in systemic lupus erythematosus patients are nuclear antigen reactive and this feature correlates with gut permeability marker levels Systemic lupus erythematosus (SLE) is characterized by the production of anti-nuclear autoantibodies. Here, for the first time, we show that the abundances of gut permeability marker Zonulin and IgA1- and IgA2-subclasses are significantly higher in the fecal samples of SLE patients compared to HCs. Importantly, IgA-total, IgA1- and IgA2-subclasses from SLE patients showed higher nAg reactivity titers. Notably, we found that not only the nuclear antigen (nAg) reactive fecal IgA1:IgA2 ratio is higher in SLE patients, but also the abundance and nAg reactivity of fecal IgA and subclasses, IgA1 particularly, correlate with the fecal levels of Zonulin, which is produced primarily in the small intestine. These observations suggest that higher amounts of nAg-reactive IgA are produced, particularly, in the proximal gut as indicated by higher IgA1:IgA2 abundance and nAg reactivity ratios and correlation of these features with the Zonulin levels. Nuclear Progestin Receptor Mediated Linkage of Blood Coagulation and Ovulation Ovulation is a dramatic remodeling process that includes rupture of blood capillaries and clotting, but coagulation is not thought to directly regulate this process. Herein, we report remarkable increases of coagulation factors V (f5, [~]3145-fold) and tissue factor (f3a, [~]120-fold) in zebrafish ovarian follicle cells during ovulation. This increase was mediated through the nuclear progestin receptor (Pgr), which is essential for ovulation in zebrafish, and was totally abolished in ovarian follicular cells from pgr-/- mutants. In addition, promoter activities of f5 and f3a were significantly enhanced by progestin (DHP) via Pgr. Similar regulation of human F5 promoter activity was induced via human PGRB, suggesting a conserved mechanism. Site-directed mutagenesis of the zebrafish f5 promoter further demonstrated a direct regulation of coagulation factors via progestin response elements. Moreover, a stark increase of erythrocytes occurred in capillaries meshed in wildtype preovulatory follicles but was absent in pgr-/- mutants. Interestingly, anticoagulants significantly inhibited ovulation both in vitro and in vivo, respectively. Furthermore, reduced fecundity was observed in f5+/- female zebrafish. Taken together, our study provides plausible evidence for steroid regulation of coagulation factors, and a new hypothesis for blood clotting triggered ovulation in vertebrates. High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. Structural basis for the inhibition of the Bacillus subtilis c-di-AMP cyclase CdaA by the phosphoglucomutase GlmM Cyclic-di-adenosine monophosphate (c-di-AMP) is an important nucleotide signalling molecule, which plays a key role in osmotic regulation in bacteria. Cellular c-di-AMP levels are tightly regulated, as both high and low levels have a negative impact on bacterial growth. Here, we investigated how the activity of the main Bacillus subtilis c-di-AMP cyclase CdaA is regulated by the phosphoglucomutase GlmM. c-di-AMP is produced from two molecules of ATP by proteins containing a deadenylate cyclase (DAC) domain. CdaA is a membrane-linked cyclase with an N-terminal transmembrane domain followed by the cytoplasmic DAC domain. Here we show, using the soluble catalytic B. subtilis CdaACD domain and purified full-length GlmM or the GlmMF369 variant lacking the C-terminal flexible domain 4, that the cyclase and phosphoglucomutase form a stable complex in vitro and that GlmM is a potent cyclase inhibitor. We determined the crystal structure of the individual B. subtilis CdaACD and GlmM proteins, both of which form dimers in the structures, and of the CdaACD:GlmMF369 complex. In the complex structure, a CdaACD dimer is bound to a GlmMF369 dimer in such a manner that GlmM blocks the oligomerization of CdaACD and formation of active head-to-head cyclase oligomers, thus providing molecular details on how GlmM acts as cyclase inhibitor. The function of a key amino acid residue in CdaACD in complex formation was confirmed by mutagenesis analysis. As the amino acids at the CdaACD:GlmM interphase are conserved, we propose that the observed inhibition mechanism of CdaA by GlmM is conserved among Firmicutes. High Spatial Resolution MALDI Imaging Mass Spectrometry ofFresh-Frozen Bone Bone and bone marrow are vital to mammalian structure, movement, and immunity. These tissues are also commonly subjected to molecular alterations giving rise to debilitating diseases like rheumatoid arthritis, osteoporosis, osteomyelitis, and cancer. Technologies such as matrix-assisted laser desorption/ionization (MALDI) imaging mass spectrometry (IMS) enable the discovery of spatially resolved chemical information in biological tissue samples to help elucidate the complex molecular processes underlying pathology. Traditionally, preparation of native osseous tissue for MALDI IMS has been difficult due to the mineralized composition and heterogenous morphology of the tissue, and compensation for these challenges with decalcification and fixation protocols can remove or delocalize molecular species. Here, sample preparation methods were advanced to enable multimodal MALDI IMS of undecalcified, fresh-frozen murine femurs allowing the distribution of endogenous lipids to be linked to specific tissue structures and cell types. Adhesive-bound bone sections were mounted onto conductive glass slides with a microscopy-compatible glue and freeze-dried to minimize artificial bone marrow damage. Subliming matrix does not induce further bone marrow cracking, and recrystallizing the deposited matrix improves lipid signal. High spatial resolution (10 m) MALDI IMS was employed to characterize lipid distributions in fresh-frozen bone, and use of complementary microscopy modalities aided tissue and cell assignments. For example, various phosphatidylcholines localize to bone marrow, adipose tissue, marrow adipose tissue, and muscle. Further, sphingomyelin(42:1) was abundant in megakaryocytes, whereas sphingomyelin(42:2) was diminished in this cell type. These data reflect the vast molecular and cellular heterogeneity indicative of the bone marrow and the soft tissue surrounding the femur. Multimodal MALDI IMS has the potential to advance bone-related biomedical research by offering deep molecular coverage with spatial relevance in a preserved native bone microenvironment. Define protein variant functions with high-complexity mutagenesis libraries and enhanced mutation detection software ASMv1.0 Open reading frame (ORF) variant libraries have advanced our ability to query the functions of a large number of variants of a protein simultaneously in a single experiment. Variant libraries targeting full-length ORFs typically consists of all possible single-amino-acid substitutions and a stop codon at each amino-acid position. Because a variant differs from the template ORF by a single codon variation, variant quantification presents the most profound challenge to this technology. Efforts such as dividing a library into sub-libraries for direct sequencing, or tag-directed subassembly are practical only for small ORFs. Our approach, however, features an enhanced variant-calling tool that processes massive shotgun sequencing data, allowing one to screen a single variant library for ORFs sized up to [~]3600 bases. The variant-calling tool detects variants reliably, and also presents variant annotations in datafiles enabling analyses that have reshaped our strategies governing library design, screen deconvolution, sequencing and sequencing data analysis. Tracheal motile cilia in mice require CAMSAP3 for formation of central microtubule pair and coordinated beating Motile cilia of multiciliated epithelial cells undergo synchronized beating to produce fluid flow along the luminal surface of various organs. Each motile cilium consists of an axoneme and a basal body, which are linked by a  transition zone. The axoneme exhibits a characteristic 9+2 microtubule arrangement important for ciliary motion, but how this microtubule system is generated is not yet fully understood. Here we show that CAMSAP3, a protein that can stabilize the minus end of a microtubule, concentrates at multiple sites of the cilium-basal body complex, including the upper region of the transition zone or the axonemal basal plate where the central pair of microtubules (CP) initiates. CAMSAP3 dysfunction resulted in loss of the CP and partial distortion of the basal plate, as well as the failure of multicilia to undergo synchronized beating. These findings suggest that CAMSAP3 plays pivotal roles in the formation or stabilization of the CP by localizing at the basal region of the axoneme, and thereby supports the coordinated motion of multicilia in airway epithelial cells. Persistent cell migration emerges from a coupling between protrusion dynamics and polarized trafficking Migrating cells present a variety of paths, from random to highly directional ones. While random movement can be explained by basal intrinsic activity, persistent movement requires stable polarization. Here, we quantitatively address emergence of persistent migration in RPE1 cells over long timescales. By live-cell imaging and dynamic micropatterning, we demonstrate that the Nucleus-Golgi axis aligns with direction of migration leading to efficient cell movement. We show that polarized trafficking is directed towards protrusions with a 20 min delay, and that migration becomes random after disrupting internal cell organization. Eventually, we prove that localized optogenetic Cdc42 activation orients the Nucleus-Golgi axis. Our work suggests that polarized trafficking stabilizes the protrusive activity of the cell, while protrusive activity orients this polarity axis, leading to persistent cell migration. Using a minimal physical model, we show that this feedback is sufficient to recapitulate the quantitative properties of cell migration in the timescale of hours. Modeling microRNA-driven post-transcriptional regulation by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. An Evolutionary Process Without Variation and Selection Natural selection successfully explains how organisms accumulate adaptive change despite that traits acquired over a lifetime are eliminated at the end of each generation. However, in some domains that exhibit cumulative, adaptive change--e.g., cultural evolution, and earliest life--acquired traits are retained; these domains do not face the problem that Darwins theory was designed to solve. Lack of transmission of acquired traits occurs when germ cells are protected from environmental change, due to a self-assembly code used in two distinct ways: (i) actively interpreted during development to generate a soma, and (ii) passively copied without interpretation during reproduction to generate germ cells. Early life and cultural evolution appear not to involve a self-assembly code used in these two ways. We suggest that cumulative, adaptive change in these domains is due to a lower-fidelity evolutionary process, and model it using Reflexively Autocatalytic and Foodset-generated networks. We refer to this more primitive evolutionary process as Self-Other Reorganisation (SOR) because it involves internal self-organising and self-maintaining processes within entities, as well as interaction between entities. SOR encompasses learning but in general operates across groups. We discuss the relationship between SOR and Lamarckism, and illustrate a special case of SOR without variation. Teosinte mexicana introgression modulates phosphatidylcholine levels and maize flowering time Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp.parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability, and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern US, Canada and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time. Actomyosin activity-dependent apical targeting of Rab11 vesicles reinforces apical constriction During tissue morphogenesis, cell shape changes resulting from cell-generated forces often require active regulation of intracellular trafficking. How mechanical stimuli influence intracellular trafficking and how such regulation impacts tissue mechanics are not fully understood. In this study, we identify an actomyosin dependent mechanism involving Rab11- mediated trafficking in regulating apical constriction in the Drosophila embryo. During Drosophila mesoderm invagination, apical actin and Myosin II (actomyosin) contractility induces apical accumulation of Rab11-marked vesicle-like structures ("Rab11 vesicles") by promoting a directional bias in dynein mediated vesicle transport. At the apical domain, Rab11 vesicles are enriched near the adherens junctions (AJs). The apical accumulation of Rab11 vesicles is essential to prevent fragmented apical AJs, breaks in the supracellular actomyosin network and a reduction in the apical constriction rate. This Rab11 function is separate from its role in promoting apical Myosin II accumulation. These findings suggest a feedback mechanism between actomyosin activity and Rab11-mediated intracellular trafficking that regulates the force generation machinery during tissue folding. A Bayesian model selection approach to mediation analysis Mediation analysis is a powerful tool for discovery of causal relationships. We describe a Bayesian model selection approach to mediation analysis that is implemented in our bmediatR software. Using simulations, we show that bmediatR performs as well or better than established methods including the Sobel test, while allowing greater flexibility in both model specification and in the types of inference that are possible. We applied bmediatR to genetic data from mice and human cell lines to demonstrate its ability to derive biologically meaningful findings. The Bayesian model selection framework is extensible to support a wide variety of mediation models.