Modeling naturalistic face processing in humans with deep convolutional neural networks Deep convolutional neural networks (DCNNs) trained for face identification can rival and even exceed human-level performance. The relationships between internal representations learned by DCNNs and those of the primate face processing system are not well understood, especially in naturalistic settings. We developed the largest naturalistic dynamic face stimulus set in human neuroimaging research (700+ naturalistic video clips of unfamiliar faces) to investigate this problem. DCNN representational geometries were weakly but significantly correlated with neural response geometries across the human face processing system. Intermediate layers better matched visual, face-selective cortices, and behavioral similarity judgments than the final fully-connected layers. Our results showed DCNNs captured only a small amount of the rich information in the neural representations during naturalistic face viewing. Future artificial neural networks trained with more ecological objective functions may help advance artificial intelligence toward the ultimate goal of mimicking human intelligence in naturalistic, real-world scenarios.
Neural basis for anxiety and anxiety-related physiological responses during a driving situation: An fMRI study While the exteroceptive and interoceptive prediction of a negative event increases a persons anxiety in daily life situations, the relationship between the brain mechanism of anxiety and the anxiety-related autonomic response has not been fully understood. In this fMRI study, we examined the neural basis of anxiety and anxiety-related autonomic responses in a daily driving situation. Participants viewed a driving video clip in the first-person perspective. During the video clip, participants were presented with a cue to indicate whether a subsequent crash could occur (attention condition) or not (safe condition). Enhanced activities in the anterior insula, bed nucleus of the stria terminalis, thalamus, and periaqueductal gray, and higher sympathetic nerve responses (pupil dilation and peripheral arterial stiffness) were triggered by the attention condition but not with the safe condition. Autonomic response-related functional connectivity was detected in the visual cortex, cerebellum, brainstem, and MCC/PCC with the right anterior insula and its adjacent regions as seed regions. Thus, the right anterior insula and adjacent regions, in collaboration with other regions play a role in eliciting anxiety based on the prediction of negative events, by mediating anxiety-related autonomic responses according to interoceptive information.
Towards a molecular mechanism underlying mitochondrial protein import through the TOM and TIM23 complexes Nearly all mitochondrial proteins need to be targeted for import from the cytosol. For the majority, the first port of call is the translocase of the outer membrane (TOM complex), followed by a procession of alternative molecular machines, conducting transport to their final destination. The pre-sequence translocase of the inner-membrane (TIM23-complex) imports proteins with cleavable pre-sequences, and comes in two flavours: the TIM23SORT complex mediates inner mitochondrial membrane (IMM) protein insertion and the TIM23MOTOR complex delivers proteins to the matrix. Progress in understanding these transport mechanisms has been hampered by the poor sensitivity and time-resolution of import assays. However, with the development of an assay based on split NanoLuc luciferase, we can now explore this process in greater detail. Here, we apply this new methodology to understand how {Delta}{psi} and ATP hydrolysis, the two main driving forces for import through the TIM23MOTOR complex, contribute to the import of pre-sequence-containing precursors (PCPs) with varying properties. Notably, we found that two major rate-limiting steps define PCP import time: passage of PCP across the OMM and initiation of IMM transport by the pre-sequence. The rates of these steps are influenced by PCP properties such as size and net charge, but not total amount of PCP imported - emphasising the importance of collecting rapid kinetic data, achieved here, to elucidate mechanistic detail. The apparent distinction between transport through the two membranes (passage through TOM is substantially complete before PCP-TIM engagement) is in contrast with the current view that import occurs through TOM and TIM in a single continuous step. Our results also indicate that PCPs spend very little time in the TIM23 channel - presumably rapid success or failure of import is critical for maintaining mitochondrial fitness. Genetic interactions drive heterogeneity in causal variant effect sizes for gene expression and complex traits Despite the growing number of genome-wide association studies (GWAS), it remains unclear to what extent gene-by-gene and gene-by-environment interactions influence complex traits in humans. The magnitude of genetic interactions in complex traits has been difficult to quantify because GWAS are generally underpowered to detect individual interactions of small effect. Here, we develop a method to test for genetic interactions that aggregates information across all trait-associated loci. Specifically, we test whether SNPs in regions of European ancestry shared between European American and admixed African American individuals have the same causal effect sizes. We hypothesize that in African Americans, the presence of genetic interactions will drive the causal effect sizes of SNPs in regions of European ancestry to be more similar to those of SNPs in regions of African ancestry. We apply our method to two traits: gene expression in 296 African Americans and 482 European Americans in the Multi-Ethnic Study of Atherosclerosis (MESA) and low-density lipoprotein cholesterol (LDL-C) in 74K African Americans and 296K European Americans in the Million Veteran Program (MVP). We find significant evidence for genetic interactions in our analysis of gene expression; for LDL-C, we observe a similar point estimate although this is not significant, likely due to lower statistical power. These results suggest that gene-by-gene or gene-by-environment interactions modify the effect sizes of causal variants in human complex traits. Chromosome counting in the mouse zygote using low-invasive super-resolution live-cell imaging In preimplantation embryos, an abnormal chromosome number causes developmental failure and a reduction in the pregnancy rate. Conventional chromosome testing methods requiring biopsy reduce the risk of associated genetic diseases; nevertheless, the reduction in cell number also reduces the pregnancy rate. Therefore, we attempted to count the chromosomes in mouse (Slc:ICR) embryos using super-resolution live-cell imaging as a new method of chromosome counting that does not reduce the cell number or viability. We counted the forty chromosomes at the first mitosis by injecting embryos with histone H2B-mCherry mRNA under conditions by which pups could be obtained; however, the results were often an underestimation of chromosome number and varied by embryo and time point. Therefore, we developed a method to count the chromosomes via CRISPR/dCas-mediated live-cell fluorescence in situ hybridization targeting the sequence of the centromere region, enabling us to count the chromosomes more accurately in mouse embryos. The methodology presented here may provide useful information for assisted reproductive technologies, such as those used in livestock animals/humans, as a technique for assessing the chromosomal integrity of embryos prior to transfer.

Impact StatementLow-toxic super-resolution observation enables chromosome counting in preimplantation embryos without cell collection. Transformations, Lineage Comparisons, and Analysis of Down to Up Protomer States of Variants of the SARS-CoV-2 Prefusion Spike Protein Including the UK Variant B.1.1.7 Monitoring and strategic response to variants in SARS-CoV-2 represents a considerable challenge in the current pandemic, as well as potentially future viral outbreaks of similar magnitude. In particular mutations and deletions involving the virions prefusion Spike protein have significant potential impact on vaccines and therapeutics that utilize this key structural viral protein in their mitigation strategies. In this study, we have demonstrated how dominant energetic landscape mappings ("glue points") coupled with sequence alignment information can potentially identify or flag key residue mutations and deletions associated with variants. Surprisingly, we also found excellent homology of stabilizing residue glue points across the lineage of {beta} coronavirus Spike proteins, and we have termed this as "sequence homologous glue points". In general, these flagged residue mutations and/or deletions are then computationally studied in detail using all-atom biocomputational molecular dynamics over approximately one microsecond in order to ascertain structural and energetic changes in the Spike protein associated variants. Specifically, we examined both a theoretically-based triple mutant and the so-called UK or B.1.1.7 variant. For the theoretical triple mutant, we demonstrated through Alanine mutations, which help "unglue" key residue-residue interactions, that these three key stabilizing residues could cause the transition of Down to Up protomer states, where the Up protomer state allows binding of the prefusion Spike protein to hACE2 host cell receptors, whereas the Down state is believed inaccessible. Thus, we are able to demonstrate the importance of glue point residue identification in the overall stability of the prefusion Spike protein. For the B.1.1.7 variant, we demonstrated the critical importance of D614G and N5017 on the structure and binding, respectively, of the Spike protein. Notably, we had previously identified D614 as a key glue point in the inter-protomer stabilization of the Spike protein prior to the emergence of its mutation. The mutant D614G is a structure breaking Glycine mutation demonstrating a relatively more distal Down state RBD and a more stable conformation in general. In addition, we demonstrate that the mutation N501Y may significantly increase the Spike protein binding to hACE2 cell receptors through its interaction with Y41 of hACE2 forming a potentially strong hydrophobic residue binding pair. We note that these two key mutations, D614G and N501Y, are also found in the so-called South African (SA; B.1.351) variant of SARS-CoV-2. Future studies along these lines are, therefore, aimed at mapping glue points to residue mutations and deletions of associated prefusion Spike protein variants in order to help identify and analyze possible "variants of interest" and optimize efforts aimed at the mitigation of this current and future virions. Citation needed? Wikipedia and the COVID-19 pandemic BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. 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. Force profile of the two-handed hardstyle kettlebell swing performed by an RKC-certified Instructor BackgroundThe effects of hardstyle kettlebell training are increasingly cited in strength and conditioning research, yet reference data from a proficient swing is scarce. The aim of this exploratory study was to investigate the force profile of a two-handed hardstyle swing performed by a Russian Kettlebell Challenge (RKC) instructor.

MethodsThe subject is a 44-year-old male, body mass 75.6 kg, height 173.5 cm, with six years of regular hardstyle training experience. Two-handed hardstyle swings were performed with a series of incremental mass kettlebells (8-68 kg). Ground reaction force (GRFs) was obtained from a floor-mounted force platform. Force-time curves (FTCs), peak force, forward force, rate of force development (RFD) and swing cadence were investigated.

ResultsData revealed the FTC of a proficient swing is highly consistent and dominated by a single force peak (mean SD = 47 N), with a profile that remained largely unchanged to 24 kg. Pearson correlation analysis revealed a very strong positive correlation in peak force with kettlebell mass (r = 0.95), which increased disproportionately from the lightest to heaviest kettlebells; net peak force increased from 8.36 {+/-} 0.75 N.kg-1 (0.85 x BW) to 12.82 {+/-} 0.39 N.kg-1 (1.3x BW). There was a strong negative correlation between RFD and kettlebell mass (r = 0.82) that decreased from 39.2 N.s-1.kg-1 to 21.5 N.s-1.kg-1. There was a very strong positive correlation in forward ground reaction force with kettlebell mass (r = 0.99), expressed as a ratio of vertical ground reaction, that increased from 0.092 (9.2%) to 0.205 (20.5%). Swing cadence exceeded 40 swings per minute (SPM) with all kettlebells.

ConclusionOur findings challenge some of the popular beliefs of the hardstyle kettlebell swing. Consistent with hardstyle practice, and previous kinematic analysis of expert and novice, force-time curves show a characteristic single large force peak, differentiating passive from active shoulder flexion. Ground reaction force did not increase proportionate to kettlebell mass, with a magnitude of forward force smaller than described in practice. These results could be useful for coaches and trainers wanting to improve athletic performance, and healthcare providers using the kettlebell swing for therapeutic purposes. Findings from this study were used to inform the BELL Trial, a pragmatic controlled trial of kettlebell training with older adults. www.anzctr.org.au ACTRN12619001177145. Non-canonical Wnt/Ror2 signaling regulates cell-matrix crosstalk to prompt directional tumor cell invasion and dissemination in breast cancer Bidirectional cell-extracellular matrix (ECM) interactions represent fundamental exchanges during tumor progression. We demonstrate the noncanonical Wnt receptor, Ror2, regulates tumor cell-driven matrix remodeling in models of breast cancer. Wnt/Ror2 loss-of-function triggers tumor cell invasion, accompanied by changes in actin cytoskeleton, adhesion, and collagen crosslinking gene expression programs. E-cadherin levels decline upon Ror2 depletion, and spatially, we pinpoint the upregulation and redistribution of 5 and {beta}3 integrins together with the production of fibronectin in areas of invasion. Wnt/{beta}-catenin-dependent and Wnt/Ror2 alternative Wnt signaling appear to regulate distinct functions for tumor cells regarding their ability to modify cell-ECM exchanges during invasion. Furthermore, blocking either integrin or focal adhesion kinase (FAK), a downstream mediator of integrin-mediated signal transduction, abrogates the enhanced migration observed upon Ror2 loss. These results reveal a critical function for the alternative Wnt receptor, Ror2, as a determinant of reciprocal communication between tumor cells and ECM during cancer invasion and metastasis. Risky choice: probability weighting explains Independence Axiom violations in monkeys Expected Utility Theory (EUT) provides axioms for maximizing utility in risky choice. The Independence Axiom (IA) is its most demanding axiom: preferences between two options should not change when altering both options equally by mixing them with a common gamble. We tested common consequence (CC) and common ratio (CR) violations of the IA over several months in thousands of stochastic choices using a large variety of binary option sets. Three monkeys showed consistently few outright Preference Reversals (8%) but substantial graded Preference Changes (46%) between the initial preferred gamble and the corresponding altered gamble. Linear Discriminant Analysis (LDA) indicated that gamble probabilities predicted most Preference Changes in CC (72%) and CR (88%) tests. The Akaike Information Criterion indicated that probability weighting within Cumulative Prospect Theory (CPT) explained choices better than models using Expected Value (EV) or EUT. Fitting by utility and probability weighting functions of CPT resulted in nonlinear and non-parallel indifference curves (IC) in the Marschak-Machina triangle and suggested IA non-compliance of models using EV or EUT. Indeed, CPT models predicted Preference Changes better than EV and EUT models. Indifference points in out-of-sample tests were closer to CPT-estimated ICs than EV and EUT ICs. Finally, while the few outright Preference Reversals may reflect the long experience of our monkeys, their more graded Preference Changes corresponded to those reported for humans. In benefitting from the wide testing possibilities in monkeys, our stringent axiomatic tests contribute critical information about risky decision-making and serve as basis for investigating neuronal decision mechanisms. Extreme Heterogeneity in Genomic Differentiation between Phenotypically Divergent Songbirds: A Test of Mitonuclear Co-introgression Comparisons of genomic variation among closely related species often show more differentiation in mitochondrial DNA (mtDNA) and sex chromosomes than in autosomes, a pattern expected due to the differing effective population sizes and evolutionary dynamics of these genomic components. Yet, introgression can cause species pairs to deviate dramatically from general differentiation trends. The yellowhammer (Emberiza citrinella) and the pine bunting (E. leucocephalos) are hybridizing avian sister species that differ greatly in appearance and moderately in nuclear DNA, but that show no mtDNA differentiation. This mitonuclear discordance is best explained by adaptive mtDNA introgression--a process that can select for co-introgression at nuclear genes with mitochondrial functions (mitonuclear genes). To better understand the extent of mitonuclear discordance and characterize nuclear differentiation patterns in this system, we investigated genome-wide differentiation between allopatric yellowhammers and pine buntings and compared it to what was seen previously in mtDNA. We found significant nuclear differentiation that was highly heterogeneous across the genome, with a particularly wide differentiation peak on the sex chromosome Z. We further tested for preferential introgression of mitonuclear genes and found statistical support for this process in yellowhammers. A role for mitonuclear coevolution in this system is supported by a stronger signal of co-introgression in genes coding for subunits of the mitoribosome and electron transport chain complexes. Altogether, our study emphasizes the extreme variation seen in differentiation across genomic components and study systems as well as highlights the ramifications of this variation in species evolution. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. Symbiodiniaceae conduct under natural bleaching stress during advanced gametogenesis stages of a mesophotic coral The mesophotic coral Alveopora allingi from the northern Gulf of Eilat/Aqaba, Red Sea, is affected by year-round partial coral-bleaching events. During these events, the migration of Symbiodiniaceae takes place from the coral-host mesoglea to the developed oocytes in bleached parts of colonies of A. allingi but not in the non-bleached parts. Additionally, these oocytes are abnormal, missing part of the structural material of the peripheral areas and are also significantly larger in the bleached areas of the colonies. Hence, we suggest a parasitic behavior of the symbionts or a commensalism relationship which enhance symbionts needs during bleaching periods and may boost the gametogenesis development in these corals. We propose that evolutionarily, this behavior may greatly contribute to the symbiont community survival throughout the bleaching period, and it can also be beneficial for the hosts persistence and adaptation to bleaching through the acquisition of a specific symbiont community following the bleaching event. Synthetic periphyton as a model system to understand species dynamics in complex microbial freshwater communities Phototrophic biofilms, also known as periphyton, are microbial freshwater communities that drive crucial ecological processes in streams and lakes. Gaining a deep mechanistic understanding of the biological processes occurring in natural periphyton remains challenging due to the high complexity and variability of such communities. To address this challenge, we rationally developed a workflow to construct a synthetic community by co-culturing 26 phototrophic species (i.e., diatoms, green algae and cyanobacteria) that were inoculated in a successional sequence to create a periphytic biofilm on glass slides. We show that this community is diverse, stable and highly reproducible in terms of microbial composition, function and 3D spatial structure of the biofilm. We also demonstrate the ability to monitor microbial dynamics at the single species level during periphyton development and how their abundances are impacted by stressors such as increased temperature and a herbicide, singly and in combination. Overall, such a synthetic periphyton, grown under controlled conditions, can be used as a model system for theory testing through targeted manipulation. Ustilago maydis telomere protein Pot1 harbors an extra N-terminal OB fold and regulates homology-directed DNA repair factors in a dichotomous and context-dependent manner The telomere G-strand binding protein Pot1 plays multifaceted roles in telomere maintenance and protection. We examined the structure and activities of Pot1 in Ustilago maydis, a fungal model that recapitulates key features of mammalian telomere regulation. Compared to the well-characterized primate and fission yeast Pot1 orthologs, UmPot1 harbors an extra N-terminal OB-fold domain (OB-N), which was recently shown to be present in most metazoans. UmPot1 binds directly to Rad51 and regulates the latters strand exchange activity. Deleting the OB-N domain, which is implicated in Rad51-binding, caused telomere shortening, suggesting that Pot1-Rad51 interaction facilitates telomere replication. Depleting Pot1 through transcriptional repression triggered growth arrest as well as rampant recombination, leading to multiple telomere aberrations. In addition, telomere repeat RNAs transcribed from both the G- and C-strand were dramatically up-regulated, and this was accompanied by elevated levels of telomere RNA-DNA hybrids. Telomere abnormalities of pot1-deficient cells were suppressed, and cell viability was restored by the deletion of genes encoding Rad51 or Brh2 (the BRCA2 ortholog), indicating that homology-directed repair (HDR) proteins are key mediators of telomere aberrations and cellular toxicity. Together, these observations underscore the complex physical and functional interactions between Pot1 and DNA repair factors, leading to context-dependent and dichotomous effects of HDR proteins on telomere maintenance and protection. A parametric logistic equation with light flux and medium concentration for cultivation planning of microalgae Microalgae are considered to be promising producers of bioactive chemicals, feeds, and fuels from carbon dioxide by photosynthesis. Thus, the prediction of microalgal growth profiles is important for the planning of cost-effective and sustainable cultivation-harvest cycles. This paper proposes a mathematical model capable of predicting the effect of light flux into culture and medium concentration on the growth profiles of microalgae by incorporating these growth-limiting factors into the logistic equation. The specific form of the equation is derived based on the experimentally measured growth profiles of Monoraphidium sp., a microalgal strain isolated by the authors, under 16 conditions consisting of the combinations of incident light fluxes into culture and initial medium concentrations. Using a cross-validation method, it is shown that the proposed model has the ability to predict necessary incident light flux into culture and initial medium concentration for harvesting target biomass at target time. Finally, the model-guided cultivation planning is performed and is evaluated by comparing the result with the experimental data. Chromosomal rearrangements with stable repertoires of genes and transposable elements in an invasive forest-pathogenic fungus AO_SCPLOWBSTRACTC_SCPLOWChromosomal rearrangements have been largely described among eukaryotes, and may have important consequences on evolution of species. High genome plasticity has been often reported in Fungi, which may explain their apparent ability to quickly adapt to new environments. Cryphonectria parasitica, causing the Chestnut blight disease, is an invasive fungal pathogen species associated with several recent host shifts during its successive introductions from Asia to North America and Europe. Previous cytological karyotyping and genomic studies suggested several chromosomal rearrangements which remains to be described in detail for this species. A serious limitation for valid genome comparisons is the access to robust genome assemblies that usually contain genomic regions of low complexity. We present a new de novo whole-genome assembly obtained from a new method of DNA extraction and long-reads sequencing Nanopore technology obtained from a Japanese isolate sampled in the native area of the species. The comparison with a recently published reference genome revealed stable gene and transposable elements (TEs) repertoires. We also showed that the C. parasitica genome is lowly compartmentalized, with a poor association between TEs and genes, such as those potentially involved in host interactions (i.e., genes coding for small secreted proteins or for secondary metabolites). This genome comparison, however, detected several large chromosomal rearrangements that may have important consequences in gene regulations and sexual mating in this invasive species. This study opens the way for more comparisons of high-quality assembled genomes, and questions the role of structural variations in the invasive success of this fungal pathogen species. Biomechanical property of limbal niche maintains stemness through YAP Stem cells (SCs) decision to self-renew or differentiate largely depends on the external control of their niche. However, the complex mechanisms that underlie this crosstalk are poorly understood. To address this question, we focused on the corneal epithelial SC model in which the SC niche, known as the limbus, is spatially segregated from the differentiation compartment. We report that the unique biomechanical property of the limbus supports the nuclear localization and function of Yes-associated protein (YAP), a putative mediator of the mechanotransduction pathway. Perturbation of tissue stiffness or YAP activity affects SC function as well as tissue integrity under homeostasis and significantly inhibited the regeneration of the SC population following SC depletion. In vitro experiments revealed that substrates with the rigidity of the corneal differentiation compartment inhibit YAP localization and induce differentiation, a mechanism that is mediated by the TGFbeta-SMAD2/3 pathway. Taken together, these results indicate that SC sense biomechanical niche signals and that manipulation of mechano-sensory machinery or its downstream biochemical output may bear fruits in SC expansion for regenerative therapy. Identifying and testing marker-trait associations for growth and phenology in three pine species: implications for genomic prediction In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 SNPs from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments. Material-mediated histogenesis using mechano-chemically microstructured cell niches Stem-cell derived tissue models are commonly cultured under globally-delivered stimuli that trigger histogenesis via self-organizing activity. However, the culture of such tissue models is prone to stochastic behavior, limiting the reproducibility of cellular composition and resulting in non-physiological architectures. To overcome these shortcomings, we developed a method for printing cell niche microenvironments with microstructured cues that mediate local histogenic processes, including mechanosensing and differentiation of selected cell types. Microstructured cues include independently tunable mechano-chemical properties, with conjugated peptides, proteins, and morphogens across a range of Youngs moduli. By rationally designing niches, we mediate the structure of tissues derived from stem-cell-progenitor sources, including a bone-fat assembly from stromal mesenchyme, and embryonic tissues derived from hiPSC. We show that microstructured cues can recapitulate mechano-chemical signals resembling early embryonic histogenesis. This outcome includes a role for niche mechanics in human embryonic organization, where soft niche mechanics bias markers of mesendodermal differentiation and epithelial-to-mesenchymal-transition (EMT), as well as a demonstration of a material-mediated morphogen signaling centers able to induce foci of mesenchymal and EMT differentiation. Thus, microstructured materials can mediate local histogenic processes to enhance the structure and composition of tissue models. Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: the Multi-Targeting Drug DREAM Challenge A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ( polypharmacology). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Author SummaryMany modern drugs are developed with the goal of modulating a single cellular pathway or target. However, many drugs are, in fact,  dirty; they target multiple cellular pathways or targets. This phenomenon is known as multi-targeting or polypharmacology. While some strive to develop  cleaner therapeutics that eliminate secondary targets, recent work has shown that multi-targeting therapeutics have key advantages for a variety of diseases. However, while multi-targeting drugs that affect a precisely-defined profile of targets may be more effective, it is difficult to computationally predict which molecules have desirable target profiles. Here, we report the results of a competitive crowdsourcing project (the Multi-Targeting Drug DREAM Challenge), where we challenged participants to predict chemicals that have desired target profiles for cancer and neurodegenerative disease. Evaluation of critical data processing steps for reliable prediction of gene co-expression from large collections of RNA-seq data MotivationGene co-expression analysis is an attractive tool for leveraging enormous amounts of public RNA-seq datasets for the prediction of gene functions and regulatory mechanisms. However, the optimal data processing steps for the accurate prediction of gene co-expression from such large datasets remain unclear. Especially the importance of batch effect correction is understudied.

ResultsWe processed RNA-seq data of 68 human and 76 mouse cell types and tissues using 50 different workflows into 7,200 genome-wide gene co-expression networks. We then conducted a systematic analysis of the factors that result in high-quality co-expression predictions, focusing on normalization, batch effect correction, and measure of correlation. We confirmed the key importance of high sample counts for high-quality predictions. However, choosing a suitable normalization approach and applying batch effect correction can further improve the quality of co-expression estimates, equivalent to a >80% and >40% increase in samples. In larger datasets, batch effect removal was equivalent to a more than doubling of the sample size. Finally, Pearson correlation appears more suitable than Spearman correlation, except for smaller datasets.

ConclusionA key point for accurate prediction of gene co-expression is the collection of many samples. However, paying attention to data normalization, batch effects, and the measure of correlation can significantly improve the quality of co-expression estimates. Creation of a single cell RNA seq meta-atlas to define human liver immune homeostasis The liver is unique in both its ability to maintain immune homeostasis and in its potential for immune tolerance following solid organ transplantation. Single-cell RNA sequencing (scRNA seq) is a powerful approach to generate highly dimensional transcriptome data to understand cellular phenotypes. However, when scRNA data is produced by different groups, with different data models, different standards, and samples processed in different ways, it can be challenging to draw meaningful conclusions from the aggregated data. The goal of this study was to establish a method to combine  human liver scRNA seq datasets by 1) characterizing the heterogeneity between studies and 2) using the meta-atlas to define the dominant phenotypes across immune cell subpopulations in healthy human liver. Publicly available scRNA seq data generated from liver samples obtained from a combined total of 17 patients and [~]32,000 cells were analyzed. Liver-specific immune cells (CD45+) were extracted from each dataset, and immune cell subpopulations (myeloid cells, NK and T cells, plasma cells, and B cells) were examined using dimensionality reduction (UMAP), differential gene expression, and ingenuity pathway analysis. All datasets co-clustered, but cell proportions differed between studies. Gene expression correlation demonstrated similarity across all studies, and canonical pathways that differed between datasets were related to cell stress and oxidative phosphorylation rather than immune-related function. Next, a meta-atlas was generated via data integration and compared against PBMC data to define gene signatures for each hepatic immune subpopulation. This analysis defined key features of hepatic immune homeostasis, with decreased expression across immunologic pathways and enhancement of pathways involved with cell death. This method for meta-analysis of scRNA seq data provides a novel approach to broadly define the features of human liver immune homeostasis. Specific pathways and cellular phenotypes described in this human liver immune meta-atlas provide a critical reference point for further study of immune mediated disease processes within the liver. 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. Deletion of Rap-Phr systems in Bacillus subtilis influences in vitro biofilm formation and plant root colonization Natural isolates of the soil-dwelling bacterium Bacillus subtilis form robust biofilms under laboratory conditions and colonize plant roots. B. subtilis biofilm gene expression displays phenotypic heterogeneity that is influenced by a family of Rap-Phr regulatory systems. Most Rap-Phr systems in B. subtilis have been studied independently, in different genetic backgrounds and under distinct conditions, hampering true comparison of the Rap-Phr systems impact on bacterial differentiation. Here, we investigated each of the 12 Rap-Phr systems of B. subtilis NCIB 3610 for their effect on biofilm formation. By studying single {Delta}rap-phr mutants, we show that despite redundancy between the cell-cell communication systems, deletion of each of the 12 Rap-Phr systems influences matrix gene expression, which could possibly enable fine-tuning of the timing and level of matrix production in response to specific conditions. Furthermore, some of the {Delta}rap-phr mutants demonstrated altered biofilm formation in vitro and colonization of Arabidopsis thaliana roots, but not necessarily similarly in both processes, indicating that the pathways regulating matrix gene expression and other factors important for biofilm formation may be differently regulated under these distinct conditions.

Significance StatementNatural isolates of Bacillus subtilis form robust biofilms in vitro and on plant roots. The formation of these heterogeneous populations is regulated by diverse Rap-Phr systems. However, most Rap-Phr systems in B. subtilis have been studied independently and in different genetic backgrounds. Here, we report that all 12 Rap-Phr systems affect matrix gene expression, while some of them affect development of in vitro biofilms and plant root colonization. Our study highlights the importance of the Rap-Phr systems in environmental adaptation of B. subtilis, specifically during biofilm formation in the rhizosphere.