Tunable Cellular Localization and Extensive Cytoskeleton-Interplay of Reflectins Reflectin proteins are nature copolymers consist of repeated canonical domains, locating in a biophotonic system called Bragg lamellae and manipulating the dynamic structural coloration of iridocytes. Their biological functions are intriguing, but the underlying mechanism is not fully understood. In this study, reflectin A1, A2, B1 and C were found to present distinguished cyto-/nucleoplasmic localization preferences. Comparable intracellular localization can be reproduced by truncated reflectin variants, suggesting a conceivable evolutionary order among reflectin proteins. Secondly, the size-dependent access of reflectin variants into nucleus demonstrate a potential model of how reflectins get into Bragg lamellae. Thirdly, RfA1 was found to extensively interact with cytoskeleton, including its binding to actin and microtubule minus-end-directed movement. This implies cytoskeleton system plays fundamental role during the organization and transportation of reflectin proteins. Findings presented here recommend reflectins as programmable biomaterials which can be used to decipher their evolution processes, to delineate their biological mechanism, and to achieve tunable intracellular targeting as editable tags.
Viral manipulation of mechanoresponsive signaling disassembles processing bodies Processing bodies (PBs) are ribonucleoprotein granules that suppress cytokine mRNA translation that are targeted for disassembly by many viruses. Kaposis sarcoma-associated herpesvirus is the etiological agent of the inflammatory endothelial cancer, Kaposis sarcoma, and a PB-regulating virus. The virus encodes Kaposin B (KapB), which induces actin stress fibres (SFs) and cell spindling as well as PB disassembly. We now show that KapB-mediated PB disassembly requires actin rearrangements, RhoA effectors and the mechanoresponsive transcription activator, YAP. Moreover, ectopic expression of active YAP or exposure of ECs to mechanical forces caused PB disassembly in the absence of KapB and mechanoresponsive PB disassembly also required YAP. Using the viral protein KapB, we identified a new consequence of the exposure of cells to mechanical forces that alter actin dynamics and activate YAP, namely the disassembly of PBs.

ImportanceFor the first time, we demonstrate that processing bodies (PBs), cytoplasmic sites of RNA decay, are regulated by mechanical signaling events that alter actin dynamics and that this requires the mechanoresponsive transcription factor, YAP. Using the overexpression of a viral protein called KapB, known previously to mediate PB disassembly, we show that actin stress fibers (SFs) and the mechanoresponsive transcription factor, YAP, are required for PB loss. We also show that other established mechanical signals (shear stress or stiff extracellular matrix) that lead to the formation of SFs and activate YAP also cause PB disassembly. This is important because it means that KapB activates, from the inside out, a pathway that links cell shape to post-transcriptional gene regulation via cytoplasmic PBs.
Fatal COVID-19 and non-COVID-19 Acute Respiratory Distress Syndrome is Associated with Incomplete Alveolar Type 1 Epithelial Cell Differentiation from the Transitional State Without Fibrosis ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. Lipid membranes modulate the activity of RNA through sequence-dependent interactions RNA is a ubiquitous biomolecule that can serve as both catalyst and information carrier. Understanding how RNA bioactivity is controlled is crucial for elucidating its physiological roles and potential applications in synthetic biology. Here we show that lipid membranes can act as RNA organization platforms, introducing a novel mechanism for ribo-regulation. The activity of R3C ribozyme can be modified by the presence of lipid membranes, with direct RNA-lipid interactions dependent on RNA nucleotide content, base pairing and length. In particular, the presence of guanine in short RNAs is crucial for RNA-lipid interactions, and G-quadruplex formation further promotes lipid binding. Lastly, by artificially modifying the R3C substrate sequence to enhance membrane binding we generated a lipid-sensitive ribozyme reaction with riboswitch-like behavior. These findings introduce RNA-lipid interactions as a tool for developing synthetic riboswitches and novel RNA-based lipid biosensors, and bear significant implications for RNA World scenarios for the origin of life. Single-cell Transcriptomic Profiling Unveils Cardiac Cell-type Specific Response to Maternal Hyperglycemia Underlying the Risk of Congenital Heart Defects Congenital heart disease (CHD) is the most prevalent structural malformations of the heart affecting [~]1% of live births. To date, both damaging genetic variations and adverse environmental exposure such as maternal diabetes have been found to cause CHD. Clinical studies show [~]fivefold higher risk of CHD in the offspring of mothers with pregestational diabetes. Maternal pregestational diabetes affects the gene regulatory networks key to proper cardiac development in the fetus. However, the cell-type specificity of these gene regulatory responses to maternal diabetes and their association with the observed cardiac defects in the fetuses remains unknown. To uncover the transcriptional responses to maternal diabetes in the early embryonic heart, we used an established murine model of pregestational diabetes. In this model, we have previously demonstrated an increased incidence of CHD. Here, we show maternal hyperglycemia (matHG) elicits diverse cellular responses during heart development by single-cell RNA-sequencing in embryonic hearts exposed to control and matHG environment. Through differential gene-expression and pseudotime trajectory analyses of this data, we identified changes in lineage specifying transcription factors, predominantly affecting Isl1+ second heart field progenitors and Tnnt2+ cardiomyocytes with matHG. Using in vivo cell-lineage tracing studies, we confirmed that matHG exposure leads to impaired second heart field-derived cardiomyocyte differentiation. Finally, this work identifies matHG-mediated transcriptional determinants in cardiac cell lineages elevate CHD risk and show perturbations in Isl1-dependent gene-regulatory network (Isl1-GRN) affect cardiomyocyte differentiation. Functional analysis of this GRN in cardiac progenitor cells will provide further mechanistic insights into matHG-induced severity of CHD associated with diabetic pregnancies. EEG-based neurofeedback with network components extraction: a data-driven approach by multilayer ICA extension and simultaneous EEG-fMRI measurements Several studies have reported advanced treatments for depressive symptoms, such as real-time neurofeedback (NF) with functional MRI (fMRI) and/or electroencephalogram (EEG). NF focusing on a regularization of brain activity associated with the amygdala or functional connectivity (FC) between the executive control network (ECN) and default mode network (DMN) has been applied to reduce depressive symptoms. However, it is practically difficult to install the fMRI-NF system and to consistently provide treatment, because of high cost. Additionally, no practical signal processing techniques have been developed extracting FC-related features from EEG signals, particularly when no physical forward models are available. In this regard, stacked pooling and linear components estimation (SPLICE), recently proposed as a multilayer extension of independent component analysis (ICA) and related independent subspace analysis (ISA), can be a promising alternative. The resting-state EEG network features can be correlated with fMRI network activity corresponding to the DMN or ECN. This may enable the modulation of the target FC-related features in EEG-based NF.

In this study, we developed a real-time EEG NF system for improving depressive symptoms by using the SPLICE. Utilizing information from the fMRI biomarkers, we evaluated our paradigm for effectiveness with regard to upregulation of the dorsolateral prefrontal cortex /middle frontal gyrus or downregulation of the precuneus/posterior cingulate cortex. We conducted an NF experiment in participants with subclinical depression; the participants were divided into the NF group (n=8) and the sham group (n=9). We found a significant reduction and a large effect size in the rumination response scale (RRS) score (reflection) in the NF group, compared to the sham group.

However, we did not find a significant relationship between the training score and difference in symptoms. This suggests that increased controllability of the EEG signals did not directly reduce the RRS reflection score. This could be due to various reasons such as improper feature extraction, individual differences, and the targeted brain regions. In this paper, we also discuss the possible ways to modify our NF protocol including the design of the experiment, sample size, and online processing. We then discuss way to improve the NF training, based on our results. Foraging with friends, moving with anyone: effects of the social environment on movement-integrated habitat selection O_LIMovement links the distribution of habitats with the social environment of animals using those habitats; yet integrating movement, habitat selection, and socioecology remains an opportunity for further study.
C_LIO_LIHere, our objective was to disentangle the roles of habitat selection and social association as drivers of collective movement in a gregarious ungulate. To accomplish this objective, we (1) assessed whether socially familiar individuals form discrete social communities and whether social communities have high spatial, but not necessarily temporal, overlap; and (2) we modelled the relationship between collective movement and selection of foraging habitats using socially informed integrated step selection analysis.
C_LIO_LIWe used social network analysis to assign individuals to social communities and determine short and long-term social preference among individuals. Using integrated step selection functions (iSSF), we then modelled the effect of social processes, i.e., nearest neighbour distance and social preference, and movement behaviour on patterns of habitat selection.
C_LIO_LIBased on assignment of individuals to social communities and home range overlap analyses, individuals assorted into discrete social communities, and these communities had high spatial overlap. By unifying social network analysis with iSSF, we identified movement-dependent social association, where individuals foraged with more familiar individuals, but moved collectively with any between foraging patches.
C_LIO_LIOur study demonstrates that social behaviour and space use are inter-related based on spatial overlap of social communities and movement-dependent habitat selection. Movement, habitat selection, and social behaviour are linked in theory. Here, we put these concepts into practice to demonstrate that movement is the glue connecting individual habitat selection to the social environment.
C_LI A family of unusual immunoglobulin superfamily genes in an invertebrate histocompatibility complex Most colonial marine invertebrates are capable of allorecognition, the ability to distinguish between themselves and conspecifics. One long-standing question is whether invertebrate allorecognition genes are homologous to vertebrate histocompatibility genes. In the cnidarian Hydractinia symbiolongicarpus, allorecognition is controlled by at least two genes, Allorecognition 1 (Alr1) and Allorecognition 2 (Alr2), which encode highly polymorphic cell surface proteins that serve as markers of self. Here, we show that Alr1 and Alr2 are part of a family of 41 Alr genes, all of which reside a single genomic interval called the Allorecognition Complex (ARC). Using sensitive homology searches and highly accurate structural predictions, we demonstrate that the Alr proteins are members of the immunoglobulin superfamily (IgSF) with V-set and I-set Ig domains unlike any previously identified in animals. Specifically, their primary amino acid sequences lack many of the motifs considered diagnostic for V-set and I-set domains, yet they adopt secondary and tertiary structures nearly identical to canonical Ig domains. Thus, the V-set domain, which played a central role in the evolution of vertebrate adaptive immunity, was present in the last common ancestor of cnidarians and bilaterians. Unexpectedly, several Alr proteins also have immunoreceptor tyrosine-based activation motifs (ITAMs) and immunoreceptor tyrosine-based inhibitory motifs (ITIMs) in their cytoplasmic tails, suggesting they could participate in pathways homologous to those that regulate immunity in humans and flies. This work expands our definition of the IgSF with the addition of a family of unusual members, several of which play a role in invertebrate histocompatibility.

Significance StatementThe immunoglobulin superfamily (IgSF) is one of the largest and most functionally versatile domain families in animal genomes. Although their amino acid sequences can vary considerably, IgSF domains have been traditionally defined by conserved residues at several key positions in their fold. Here, we sequenced an invertebrate histocompatibility complex and discovered a family of IgSF genes with amino acid sequences that lack most of these residues yet are predicted to adopt folds virtually identical to canonical V-set and I-set IgSF domains. This work broadens the definition of the IgSF and shows that the V-set domain was present earlier in animal evolution than previously appreciated. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. Antibody Resistance of SARS-CoV-2 Variants B.1.351 and B.1.1.7 The COVID-19 pandemic has ravaged the globe, and its causative agent, SARS-CoV-2, continues to rage. Prospects of ending this pandemic rest on the development of effective interventions. Single and combination monoclonal antibody (mAb) therapeutics have received emergency use authorization1-3, with more in the pipeline4-7. Furthermore, multiple vaccine constructs have shown promise8, including two with ~95% protective efficacy against COVID-199,10. However, these interventions were directed toward the initial SARS-CoV-2 that emerged in 2019. The recent emergence of new SARS-CoV-2 variants B.1.1.7 in the UK11 and B.1.351 in South Africa12 is of concern because of their purported ease of transmission and extensive mutations in the spike protein. We now report that B.1.1.7 is refractory to neutralization by most mAbs to the N-terminal domain (NTD) of spike and relatively resistant to a few mAbs to the receptor-binding domain (RBD). It is not more resistant to convalescent plasma or vaccinee sera. Findings on B.1.351 are more worrisome in that this variant is not only refractory to neutralization by most NTD mAbs but also by multiple individual mAbs to the receptor-binding motif on RBD, largely due to an E484K mutation. Moreover, B.1.351 is markedly more resistant to neutralization by convalescent plasma (9.4 fold) and vaccinee sera (10.3-12.4 fold). B.1.351 and emergent variants13,14 with similar spike mutations present new challenges for mAb therapy and threaten the protective efficacy of current vaccines. Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data. Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. 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. 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. Production and characterization of virus-free, CRISPR-CAR T cells capable of inducing solid tumor regression BackgroundChimeric antigen receptor (CAR) T cells traditionally harbor viral vector-based sequences that encode the CAR transgene in the genome. These T cell products have yet to show consistent anti-tumor activity in patients with solid tumors. Further, viral vector manufacturing is resource intensive, suffers from batch-to-batch variability, and includes several animal components, adding regulatory and supply chain pressures.

MethodsAnti-GD2 CAR T cells were generated using CRISPR/Cas9 within nine days using recombinant Cas9 protein and nucleic acids, without any viral vectors or animal components. The CAR was specifically targeted to the T Cell Receptor Alpha Constant gene (TRAC). T cell products were characterized at the level of the genome, transcriptome, proteome, and secretome using CHANGE-seq, scRNA-seq, spectral cytometry, and ELISA assays. Functionality was evaluated in vivo in an NSG xenograft neuroblastoma model.

ResultsIn comparison to traditional retroviral CAR T cells, virus-free CRISPR CAR (VFC-CAR) T cells exhibit TRAC-targeted genomic integration of the CAR transgene, elevation of transcriptional and protein characteristics associated with a memory phenotype, and low tonic signaling prior to infusion arising in part from the the knockout of the TCR. Upon exposure to the GD2 target antigen, anti-GD2 VFC-CAR T cells exhibited specific cytotoxicity against GD2+ cells in vitro and induced solid tumor regression in vivo, with robust homing, persistence, and low exhaustion against a human neuroblastoma xenograft model.

ConclusionsThis proof-of-principle study leveraging virus-free genome editing technology could enable flexible manufacturing of clinically relevant, high-quality CAR T cells to treat cancers, including solid tumors. 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 genomic and morphometric analysis of alpine bumblebees: ongoing reductions in tongue length but no clear genetic component Over the last six decades, populations of the bumblebees Bombus sylvicola and Bombus balteatus in Colorado have experienced decreases in tongue length, a trait important for plant-pollinator mutualisms. It has been hypothesized that this observation reflects selection resulting from shifts in floral composition under climate change. Here we used morphometrics and population genomics to determine whether morphological change is ongoing, investigate the genetic basis of morphological variation, and analyze population structure in these populations. We generated a genome assembly of B. balteatus. We then analyzed whole-genome sequencing data and morphometric measurements of 580 samples of both species from seven high-altitude localities. Out of 281 samples originally identified as B. sylvicola, 67 formed a separate genetic cluster comprising a newly-discovered cryptic species ("incognitus"). However, an absence of genetic structure within species suggests that gene flow is common between mountains. We found a significant decrease in tongue length between bees collected between 2012-2014 and in 2017, indicating that morphological shifts are ongoing. We did not discover any genetic associations with tongue length, but a SNP related to production of a proteolytic digestive enzyme was implicated in body size variation. We identified evidence of covariance between kinship and both tongue length and body size, which is suggestive of a genetic component of these traits, although it is possible that shared environmental effects between colonies are responsible. Our results provide evidence for ongoing modification of a morphological trait important for pollination and indicate that this trait likely has a complex genetic and environmental basis. 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. Interpretable and Generalizable Attention-Based Model for Predicting Drug-Target Interaction Using 3D Structure of Protein Binding Sites: SARS-CoV-2 Case Study and in-Lab Validation In this study, we introduce and implement an interpretable graph-based deep learning prediction model, which utilizes protein binding sites along with self-attention to learn which protein binding sites interact with a given ligand. Our proposed model enables interpretability by identifying the protein binding sites that contribute the most towards the Drug-Target Interaction. Results on three benchmark datasets show improved performance compared to previous graph-based models. More significantly, unlike previous studies our model performance remains close to the optimal performance when tested with new proteins (ie., high generalizablity). Through multidisciplinary collaboration, we further experimentally evaluate the practical potential of our proposed approach. To achieve this, we first computationally predict binding interaction of some candidate compounds with a target protein, then experimentally validate the binding interactions for these pairs in the laboratory. The high agreement between the computationally-predicted and experimentally-observed (measured) DTIs illustrates the potential of our method as an effective pre-screening tool in drug re-purposing applications. Thermodynamic constraints on the assembly and diversity ofmicrobial ecosystems are different near to and far fromequilibrium Non-equilibrium thermodynamics has long been an area of substantial interest to ecologists because most fundamental biological processes, such as protein synthesis and respiration, are inherently energy-consuming. However, most of this interest has focused on developing coarse ecosystem-level maximisation principles, providing little insight into underlying mechanisms that lead to such emergent constraints. Microbial communities are a natural system to decipher this mechanistic basis because their interactions in the form of substrate consumption, metabolite production, and cross-feeding can be described explicitly in thermodynamic terms. Previous work has considered how thermodynamic constraints impact competition between pairs of species, but restrained from analysing how this manifests in complex dynamical systems. To address this gap, we develop a thermodynamic microbial community model with fully reversible reaction kinetics, which allows direct consideration of free-energy dissipation. This also allows species to interact via products rather than just substrates, increasing the dynamical complexity, and allowing a more nuanced classification of interaction types to emerge. Using this model, we find that community diversity increases with substrate lability, because greater free-energy availability allows for faster generation of niches. Thus, more niches are generated in the time frame of community establishment, leading to higher final species diversity. We also find that allowing species to make use of near-to-equilibrium reactions increases diversity in a low free-energy regime. In such a regime, two new thermodynamic interaction types that we identify here reach comparable strengths to the conventional (competition and facilitation) types, emphasising the key role that thermodynamics plays in community dynamics. Our results suggest that accounting for realistic thermodynamic constraints is vital for understanding the dynamics of real-world microbial communities.

Author summaryThere is a growing interest in microbial communities due to their important role in biogeochemical cycling as well as plant and animal health. Although our understanding of thermodynamic constraints on individual cells is rapidly improving, the impact of these constraints on complex microbial communities remains largely unexplored theoretically and empirically. Here, we develop a new microbial community model which allows thermodynamic efficiency and entropy production to be calculated directly. We find that availability of substrates with greater free-energy allows for a faster rate of niche generation, leading to higher final species diversity. We also show that when the free-energy availability is low, species with reactions close to thermodynamic equilibrium are favoured, leading to more diverse and efficient communities. In addition to the conventional interaction types (competition and facilitation), our model reveals the existence of two novel interaction types mediated by products rather than substrates. Though the conventional interactions are generally the strongest, the novel interaction types are significant when free-energy availability is low. Our results suggest that non-equilibrium thermodynamics need to be considered when studying microbial community dynamics. Embryo development in Carica papaya Linn. Papaya (Carica papaya Linn.) is a tropical plant whose draft genome has been sequenced. Papaya produces large fruits rich in vitamins A and C and is an important cash crop in developing countries. Nonetheless, little is known about how the female gametophyte develops, how it is fertilized and how it develops into a mature seed containing an embryo and an endosperm. The Papaya female gametophyte displays a Polygonum-type architecture consisting of two synergid cells, an egg cell, a central cell, and three antipodal cells. Reports are available of the presumed existence of varieties in which cross fertilization is bypassed and autonomous development of embryos occurs (e.g., apomixis). In this study, we analyzed the development of female gametophytes in a commercial Hawaiian parental line and in the presumed apomictic Costa Rican line L1. Samples were collected before and after anthesis to compare the overall structure, size and transcriptional patterns of several genes that may be involved in egg and endosperm cell fate and proliferation. These genes were the putative papaya homologs of ARGONAUTE9 (AGO9), MEDEA (MEA), RETINOBLASTOMA RELATED-1 (RBR1), and SLOW WALKER-1 (SWA1). Our results suggest that its feasibly to identify structural features of Polygonum-type development, and that in bagged female flowers of line L1 we might have observed autonomous development of embryo-like structures. Possible downregulation of papaya homologs for AGO9, MEA, RBR1 and SWA1 was observed in embryo sacs from line L1 before and after anthesis, which may suggest a link between suspected apomixis and transcriptional downregulation of genes for RNA-directed DNA methylation, histone remodelers, and rRNA processing. Most notably, the large size of the papaya embryo sac suggests that it could be a cytological alternative to Arabidopsis thaliana for study. Significant variation in embryo sac size was observed between the varieties under study, suggesting wide differences in the genetic regulation of anatomical features. Orientation dependent proton transverse relaxation in the human brain white matter: The magic angle effect on a cylindrical helix PurposeTo overcome some limitations of prior orientation-dependent proton transverse relaxation formalisms in white matter (WM) with a novel framework based on the generalized magic angle effect function.

MethodsA cylindrical helix model was developed embracing anisotropic rotational and translational diffusion of restricted molecules in human brain WM, with the former characterized by an axially symmetric system. Transverse relaxation rates R2 and [Formula] were divided into isotropic [Formula] and anisotropic parts, [Formula], with  denoting an open angle and{varepsilon} 0 an orientation ({Phi}) offset from DTI-derived primary diffusivity direction. The proposed framework (Fit A) was compared with prior models without{varepsilon} 0 on previously published water and methylene proton transverse relaxation rates from developing, healthy, and pathological WM at 3T. Goodness of fit was represented by root-mean-square error (RMSE). F-test and linear correlation were used with statistical significance set to P [&le;] 0.05.

ResultsFit A significantly (P<0.01) outperformed prior models as demonstrated by reduced RMSEs, e.g., 0.349 vs. 0.724 in myelin water. Fitted{varepsilon} 0 was in good agreement with calculated{varepsilon} 0 from directional diffusivities. Compared with those from healthy adult, the fitted [Formula], and  from neonates were substantially reduced but{varepsilon} 0 increased, consistent with incomplete myelination. Significant positive [Formula] and negative ( and [Formula]) correlations were found with aging (demyelination) in elderly.

ConclusionThe developed framework can better characterize orientation dependences from a wide range of proton transverse relaxation measurements in human brain WM, shedding new light on myelin microstructural alterations at the molecular level. C. elegans PEZO-1 is a Mechanosensitive Ion Channel Involved in Food Sensation PIEZO channels are force sensors essential for physiological processes, including baroreception and proprioception. The Caenorhabditis elegans genome encodes an ortholog gene of the Piezo family, pezo-1, which is expressed in several tissues, including the pharynx. This myogenic pump is an essential component of the C. elegans alimentary canal, whose contraction and relaxation are modulated by mechanical stimulation elicited by food content. Whether pezo-1 encodes a mechanosensitive ion channel and contributes to pharyngeal function remain unknown. Here, we leverage genome editing, genetics, microfluidics, and electropharyngeogram recording to establish that pezo-1 is expressed in the pharynx, including in a proprioceptive-like neuron, and regulates pharyngeal function. Knockout (KO) and gain-of-function (GOF) mutants reveal that pezo-1 is involved in fine-tuning pharyngeal pumping frequency, sensing osmolarity, and food mechanical properties. Using pressure-clamp experiments in primary C. elegans embryo cultures, we determine that pezo-1 KO cells do not display mechanosensitive currents, whereas cells expressing wild-type or GOF PEZO-1 exhibit mechanosensitivity. Moreover, infecting the Spodoptera frugiperda cell line with a baculovirus containing the G-isoform of pezo-1 (among the longest isoforms) demonstrates that pezo-1 encodes a mechanosensitive channel. Our findings reveal that pezo-1 is a mechanosensitive ion channel that regulates food sensation in worms. Contribution Evenness: A functional redundancy metric sensitive to trait stability in microbialcommunities The concept of functional redundancy has received considerable attention in both the macroecology and microbial ecology literature. As a result, multiple metrics of functional redundancy have been proposed. These vary in how they weight trait levels, species abundance, functional richness, and species richness. Here we present a new functional redundancy metric tailored for community-aggregated traits, which are traits that are quantified at the community level and can be quantitatively partitioned among species. We call this metric Contribution Evenness (CE) because it measures how evenly species contribute to a community-aggregated trait. As CE is an evenness measurement, it ranges from 0 and 1, where 0 corresponds to a single species contributing to a community-aggregated trait and 1 corresponds to all species contributing equally. Using in silico simulations of species extinctions, we demonstrate that CE reflects the stability of an ecosystem function to species extinction, a hypothesized ecological consequence of functional redundancy. As a positive control and to illustrate how CE can be used with sequence data, we analyzed the functional redundancy of eight nitrogen-transforming pathways using 2,631 metagenome-assembled genomes from 47 marine sites. CE for marine nitrogen cycle marker genes was consistent with our qualitative understanding of which nitrogen pathways are most functionally redundant in the ocean. We found that, on average, the NH4+ assimilation pathway was the most functionally redundant (0.44 {+/-} 0.08) while dissimilatory nitrate reduction was the least redundant (0.005 {+/-} 0.005). As demonstrated here, CE provides a promising framework for measuring trait stability in microbiomes. Phenotypic heterogeneity driven by plasticity of the intermediate EMT state governs disease progression and metastasis in breast cancer The Epithelial-to-Mesenchymal Transition (EMT) is a developmental cellular program frequently coopted by cancer cells1 and is a key contributor to both heterogeneity in solid tumors2-4 and later stage chemo-resistance and metastasis5,6. Rather than being a switch from an epithelial to a mesenchymal state, increasing evidence points to the existence of intermediate EMT states, wherein cells co-express both epithelial and mesenchymal traits7-13. Multiple stable intermediate EMT states possessing unique characteristics exist across the EMT spectrum7,8,14,15, contributing to the complex heterogeneity of tumors to promote metastasis16. While much work has been carried out identifying and characterizing EMT-inducing transcription factors17-20, the transcriptional and epigenetic networks responsible for the stability and maintenance of the midpoints along the EMT spectrum are poorly defined. In addition, there are currently no approaches to identifying and quantifying intermediate EMT subpopulations within patient tumors to evaluate their prognostic significance. Using clonally isolated derivatives of the SUM149PT breast cancer cell line, we systematically interrogate how each EMT state independently contributes to heterogeneity and influences metastatic progression, uncovering the role of RUNX2 in stabilizing certain intermediate states. Using SUM149PT-derived tumors as a training set, we develop an entropy-based model to quantify phenotypic heterogeneity and EMT status. Remarkably, the majority of cell states captured in the SUM149PT model are represented in a selection of patient tumors, laying the foundation for quantification of epithelial-mesenchymal heterogeneity and understanding the role of the intermediate EMT state in tumor progression. Genome-wide association study suggests an independent genetic basis of zinc and cadmium concentrations in fresh sweet corn kernels Despite being one of the most consumed vegetables in the United States, the elemental profile of sweet corn (Zea mays L.) is limited in its dietary contributions. To address this through genetic improvement, a genome-wide association study was conducted for the concentrations of 15 elements in fresh kernels of a sweet corn association panel. In concordance with mapping results from mature maize kernels, we detected a probable pleiotropic association of zinc and iron concentrations with nicotianamine synthase5 (nas5), which purportedly encodes an enzyme involved in synthesis of the metal chelator nicotianamine. Additionally, a pervasive association signal was identified for cadmium concentration within a recombination suppressed region on chromosome 2. The likely causal gene underlying this signal was heavy metal ATPase3 (hma3), whose counterpart in rice, OsHMA3, mediates vacuolar sequestration of cadmium and zinc in roots, whereby regulating zinc homeostasis and cadmium accumulation in grains. In our association panel, hma3 associated with cadmium but not zinc accumulation in fresh kernels. This finding implies that selection for low cadmium will not affect zinc levels in fresh kernels. Although less resolved association signals were detected for boron, nickel, and calcium, all 15 elements were shown to have moderate predictive abilities via whole-genome prediction. Collectively, these results help enhance our genomics-assisted breeding efforts centered on improving the elemental profile of fresh sweet corn kernels.