Gene duplications and phylogenomic conflict underlie major pulses of phenotypic evolution in gymnosperms Inferring the intrinsic and extrinsic drivers of species diversification and phenotypic disparity across the Tree of Life is a major challenge in evolutionary biology. In green plants, polyploidy (or whole-genome duplication, WGD) is known to play a major role in microevolution and speciation1, but the extent to which WGD has shaped macroevolutionary patterns of diversification and phenotypic innovation across plant phylogeny remains an open question. Here we examine the relationship of various facets of genomic evolution--including gene and genome duplication, genome size, and chromosome number--with macroevolutionary patterns of phenotypic innovation, species diversification, and climatic occupancy in gymnosperms. We show that genomic changes, such as WGD and genome-size shifts, underlie the origins of most major extant gymnosperm clades, and notably our results support an ancestral WGD in the gymnosperm lineage. Spikes of gene duplication typically coincide with major spikes of phenotypic innovation, while increased rates of phenotypic evolution are typically found at nodes with high gene-tree conflict, representing historic population-level dynamics during speciation. Most shifts in gymnosperm diversification since the rise of angiosperms are decoupled from putative WGDs and instead are associated with increased rates of climatic occupancy evolution, particularly in cooler and/or more arid climatic conditions, suggesting that ecological opportunity, especially in the later Cenozoic, and environmental heterogeneity have driven a resurgence of gymnosperm diversification. Our study provides critical insight on the processes underlying diversification and phenotypic evolution in gymnosperms, with important broader implications for the major drivers of both micro- and macroevolution in plants.
Factors behind intraspecific competition during post-glacial recolonization: a Norway spruce case study O_LITo a large extent, the current distribution and population structure of many species were shaped by the cycles of isolation in glacial refugia and subsequent population expansions. Isolation in, and post-glacial expansion through, heterogeneous environments led to either neutral or adaptive divergence. Norway spruce is no exception, and its current distribution is the consequence of a complex interplay between evolutionary forces.
C_LIO_LIWe investigated population differentiation and adaptations of Norway spruce for juvenile growth, diameter of the trunk at breast height and its components: wood density and tracheid traits. Data from 4,461 phenotyped and genotyped Norway spruce from 396 half-sib families in two progeny tests were used to test for divergent selection in the framework of QST vs FST, and to carry out multiple Genome-Wide Association Studies.
C_LIO_LIWe show that (i) the macroscopic resultant trait (trunk diameter), unlike its microscopic components (tracheid dimensions) and juvenile growth, was under divergent selection, and (ii) adaptation predated the Last Glaciation Maximum.
C_LIO_LIAltogether, the current variation in these phenotypic traits in Norway spruce is better explained by local adaptation to ancestral environments than to current ones, where populations were partly pre-adapted mainly through growth-related traits.
C_LI
Cell Cycle Delay, Pro-metaphase Arrest and C-metaphase Inducing Effects of Petroleum Ether Fraction of Leaf Aqueous Extract of Clerodendrum viscosum Vent. Clerodendrum viscosum is a traditionally used medicinal plant. The present study aimed to analyze cell cycle delay, pro-metaphase arrest, and c-metaphase inducing effects of the petroleum ether fraction (AQPEF) of leaf aqueous extract of C. Viscosum (LAECV). The LAECV was fractionated with petroleum ether and its metaphase arrest, cell cycle delay, and c-metaphase inducing activities were tested on A. cepa root tip cells. The AQPEF induced cell cycle delay, and colchicine like metaphase, c-metaphase, in A. cepa root tip cells. Thus, the present study explores AQPEF as an active fraction of LAECV having metaphase arresting, cell cycle delay, and c-metaphase inducing potentials. bMSI-CAST: a systematic method for next generation sequencing-based microsatellite instability detection in plasma cell-free DNA Microsatellite instability (MSI) is a well-established prognostic and predictive biomarker in certain types of cancers. MSI detection using tumour tissue is often limited by the availability of specimens. Next generation sequencing (NGS)-based MSI detection in plasma cell-free DNA (cfDNA) is challenged by a much lower signal-to-noise ratio. We developed a highly accurate cfDNA MSI detection method called bMSI-CAST (blood MSI Caller Adjusted with Sequence duplicaTes), with improvement on three features including a set of locus selection principles ensuring loci with high robustness and compatibility across sequencing platforms, an MSI-specific duplicate removal strategy, and a calling algorithm that dynamically matches baselines with a broad range of duplication levels. Analytical validation via MSI-high (MSI-H) cell gDNA showed an LOD of 0.15%. Furthermore, in an analysis of 95 evaluable cfDNA samples from patients with gastrointestinal cancers, bMSI-CAST exhibited a positive predictive agreement (PPA) of 92.9% (39/42) and negative predictive agreement (NPA) of 100% (53/53) with tissue MSI-PCR. In conclusion, bMSI-CAST provides novel and advanced solutions to key aspects fundamental to cfDNA MSI calling but not sufficiently addressed by existing methods, and it is a validated method ready to be applied to aid clinical decisions for cancer patients. Partial exogastrulation due to apical-basal polarity of F-actin distribution disruption in sea urchin embryo by omeprazole Gastrulation is a universal process in the morphogenesis of many animal embryos. Although morphological and molecular events in gastrulation have been well studied, the mechanical driving forces and underlying regulatory mechanisms are not fully understood. Here, we investigated the gastrulation of embryos of a sea urchin, Hemicentrotus pulcherrimus, which involves the invagination of a single-layered vegetal plate into the blastocoel. We observed that omeprazole, a proton pump inhibitor capable of perturbing the left-right asymmetry of sea urchin embryo, induced "partial exogastrulation" where the secondary invagination proceeds outward. During early gastrulation, intracellular apical-basal polarity of F-actin distribution in vegetal half were higher than those in animal half, while omeprazole treatment disturbed the apical-basal polarity of F-actin distribution in vegetal half. Furthermore, gastrulation stopped and even partial exogastrulation did not occur when F-actin polymerization or degradation in whole embryo was partially inhibited via RhoA or YAP1 knockout. A mathematical model of the early gastrulation reproduced the shapes of both normal and exogastrulating embryos using cell-dependent cytoskeletal features based on F-actin. Additionally, such cell position-dependent intracellular F-actin distributions might be regulated by intracellular pH distributions. Therefore, apical-basal polarity of F-actin distribution disrupted by omeprazole may induce the partial exogastrulation via anomalous secondary invagination. Natural tolerance to transposition is associated with double-strand break repair and germ-cell differentiation Transposable elements (TE) are mobile genetic parasites whose unregulated activity in the germline causes DNA damage and sterility. While the regulation of TE mobilization by hosts is studied extensively, little is known about mechanisms that could allow germline cells to persist in the face of genotoxic stress imposed by active transposition. Such tolerance mechanisms are predicted to be beneficial when new TEs invade and host repression has not yet evolved. Here we use hybrid dysgenesis--a sterility syndrome of Drosophila caused by transposition of invading DNA transposons--to uncover genetic variants that confer tolerance to transposition. Using a panel of highly recombinant inbred lines of Drosophila melanogaster, we identified two linked quantitative trait loci (QTL), that determine tolerance in young and old females, respectively. Through transcriptomic and phenotypic comparisons, we provide evidence that young tolerant females exhibit enhanced repair of double-stranded breaks, explaining their ability to withstand high germline transposition rates. We furthermore identify the germline differentiation factor brat as an independent tolerance factor, whose activity may promote germline maintenance in aging dysgenic females. Together, our work reveals the diversity of potential tolerance mechanisms across development, as well as tolerant variants that may be beneficial in the context of P-element transposition. 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. A map of binary SARS-CoV-2 protein interactions implicates host immune regulation and ubiquitination Key steps in viral propagation, immune suppression, and pathology are mediated by direct, binary, physical interactions between viral and host proteins. To understand the biology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, we generated an unbiased systematic map of binary interactions between viral and host proteins, complementing previous co-complex association maps by conveying more direct mechanistic understanding and potentially enabling targeted disruption of direct interactions. To this end, we deployed two parallel strategies, identifying 205 virus-host and 27 intraviral binary interactions amongst 171 host and 19 viral proteins, and confirming high quality of these interactions via a calibrated orthogonal assay. Host proteins interacting with SARS-CoV-2 proteins are enriched in various cellular processes, including immune signaling and inflammation, protein ubiquitination, and membrane trafficking. Specific subnetworks provide new hypotheses related to viral modulation of host protein homeostasis and T-cell regulation. The binary virus-host protein interactions we identified can now be prioritized as targets for therapeutic intervention. More generally, we provide a resource of systematic maps describing which SARS-CoV-2 and human proteins interact directly. First evidence of in vitro cytotoxic effects of marine microlitter on Merluccius merluccius and Mullus barbatus, two Mediterranean commercial fish species Marine litter is composed mainly of plastics and is recognized as a serious threats to marine ecosystems. Ecotoxicological approaches have started elucidating the potential severity of microplastics (MPs) in controlled laboratory studies with pristine materials but no information exist on marine environmental microlitter as a whole. Here, we characterized the litter in the coastal Northern Tyrrhenian sea and in the stomach of two fish species of socio-economic importance, and exposed primary cell cultures of mucosal and lymphoid organs to marine microlitter for evaluating possible cytotoxic effects. An average of 0.30 {+/-} 0.02 microlitter items m-3 was found in water samples. FT-IR analysis revealed that plastic particles, namely HDPE, polyamide and polypropylene were present in 100% and 83.3% of Merluccius merluccius and Mullus barbatus analyzed, which overall ingested 14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual, respectively. Moreover, microlitter was confirmed as a vector of microorganisms. Lastly, the apical end-point of viability was found to be significantly reduced in splenic cells exposed in vitro to two microlitter conditions. Considering the role of the spleen in the mounting of adaptive immune responses, our results warrant more in-depth investigations for clarifying the actual susceptibility of these two species to anthropogenic microlitter.

HighlightsO_LI0.30 {+/-} 0.02 microlitter items m-3 were found at the surface of coastal Northern Tyrrhenian sea
C_LIO_LI14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual were retrieved from the stomach of hakes and mullets
C_LIO_LIThe ingested microlitter contained plastic items
C_LIO_LIMicrolitter was validated as a carrier of bacteria, fungi and flagellates
C_LIO_LISplenic cells exposed to two microlitter conditions for 72 hours suffered cytotoxicity
C_LI Micro topographical instruction of bacterial attachment, biofilm formation and in vivo host response Bio-instructive materials that prevent bacterial biofilm formation and drive an appropriate host immune response have the potential to significantly reduce the burden of medical device-associated infections. As bacterial surface attachment is sensitive to surface topography, we experimentally surveyed 2,176 combinatorially generated shapes using an unbiased high-throughput micro topographical polystyrene polymer chip screen. This identified topographies that reduced colonization in vitro by up to 15-fold compared with a flat surface for both motile and non-motile bacterial pathogens. Equivalent reductions were achieved on polyurethane, a polymer commonly used in medical devices. Using machine learning methods, a set of design rules based on generalisable topo-descriptors was established for predicting bacteria-resistant micro topographies. In a murine foreign body infection model, anti-attachment topographies were shown to be refractory to P. aeruginosa colonization and serendipitously, the fibrotic response to the implant was reduced, offering control of foreign body response using simple topographical patterning of non-eluting implants. Rattus Economicus: A neural model for rational regulation of water-seeking effort in rodents In the laboratory, animals motivation to work tends to be positively correlated with reward magnitude. But in nature, rewards earned by work are essential to survival (e.g., working to find water), and the payoff of that work can vary on long timescales (e.g., seasonally). Under these constraints, the strategy of working less when rewards are small could be fatal. We found that instead, rats in a closed economy did more work for water rewards when the rewards were stably smaller, a phenomenon also observed in human labor supply curves. Like human consumers, rats also showed elasticity of demand, consuming far more water per day when its price in effort was lower. The neural mechanisms underlying such "rational" market behaviors remain largely unexplored. We propose a dynamic utility maximization model that can account for the dependence of rat labor supply (trials/day) on the wage rate (ml/trial), and also predict the temporal dynamics of when rats work. Based on data from mice, we hypothesize that SFOGLUT neurons in lamina terminalis continuously compute the instantaneous marginal utility of voluntary work for water reward, and causally determine the amount and timing of work. NGF-p75 signaling coordinates skeletal cell migration during bone repair Bone regeneration following injury is initiated by inflammatory signals and occurs in association with infiltration by sensory nerve fibers. Together, these events are believed to coordinate angiogenesis and tissue reprogramming, but the mechanism of coupling immune signals to re-innervation and osteogenesis is unknown. Here, we found that NGF is expressed following cranial bone injury and signals via p75 in resident mesenchymal osteogenic precursors to impact their migration into the damaged tissue. Mice lacking Ngf in myeloid cells demonstrated reduced migration of osteogenic precursors to the injury site with consequently delayed bone healing. These features were phenocopied by mice lacking p75 in Pdgfra+ osteoblast precursors. Single-cell transcriptomics identified mesenchymal subpopulations with potential roles in cell migration and immune response, altered in the context of p75 deletion. Together, these results identify the role of p75 signaling pathway in coordinating skeletal cell migration during early bone repair. 'rbioacc': an R-package to analyse toxicokinetic data O_LI rbioacc is an  R-package dedicated to the analysis of experimental data collected from bioaccumulation tests during which organisms are exposed to a chemical (exposure phase) and then put into a clean media (depuration phase). Internal concentrations are regularly measured over time all along the experiment.
C_LIO_LI rbioacc provides ready-to-use functions to visualize and fully analyze such data. Under a Bayesian framework, this package fits a generic one-compartment toxicokinetic (TK) model automatically built from the data. It provides TK parameter estimates (appropriate uptake and elimination rates) and bioaccumulation metrics (e.g., BCF, BSAF, BMF). All parameter estimates, bioaccumulation metrics as well as predictions of internal concentrations into organisms are delivered with their uncertainty.
C_LIO_LIBioaccumulation metrics are primarily provided in support of environmental risk assessment, in full compliance with regulatory requirements required to approve marketing applications of chemical substances.
C_LIO_LIThis paper gives brief worked examples of the use of  rbioacc from data collected through standard bioaccumulation tests, and publicly available within the scientific literature. These examples constitute step-by-step user-guides to analyze any new data set, uploaded in the right format.
C_LI Multivariate genome-wide association study on tissue-sensitive diffusion metrics identifies key molecular pathways for axonal growth, synaptogenesis, and astrocyte-mediated neuroinflammation The molecular determinants of tissue composition of the human brain remain largely unknown. Recent genome-wide association studies (GWAS) on this topic have had limited success due to methodological constraints. Here, we apply advanced whole-brain analyses on multi-shell diffusion imaging data and multivariate GWAS to two large scale imaging genetic datasets (UK Biobank and the Adolescent Brain Cognitive Development study) to identify and validate genetic association signals. We discovered 503 unique genetic loci that have impact on multiple regions of human brain. Among them, more than 79% were validated in either of two large-scale independent imaging datasets. Key molecular pathways involved in axonal growth, astrocyte-mediated neuroinflammation, and synaptogenesis during development were found to significantly impact the measured variations in tissue-specific imaging features. Our results shed new light on the biological determinants of brain tissue composition and their potential overlap with the genetic basis of neuropsychiatric disorders. SUMM4 biochemically complex couples insulator function and DNA replication timingcontrol The mechanisms that establish DNA replication timing programs in eukaryotes remain incompletely understood. Drosophila SNF2-related factor SUUR imparts under-replication (UR) of late-replicating intercalary heterochromatin (IH) in polytene chromosomes. We developed a proteomics technique termed MERCI to isolate a native complex SUMM4 comprising SUUR and chromatin boundary protein Mod(Mdg4)-67.2. Mod(Mdg4) stimulates the ATPase activity of SUUR and is required for its normal spatiotemporal distribution in vivo. Both SUMM4 subunits mediate the activities of gypsy insulator disrupting enhancer-promoter interactions and establishing chromatin barriers. Furthermore, SuUR or mod(mdg4) mutations reverse UR of IH. Our findings uncover a critical role for architectural proteins in attenuating replication fork progression and suggest an alternative mechanism for DNA replication timing that does not depend on an asynchronous firing of replication origins.

One-Sentence SummaryA stable protein complex comprising an insulator factor and a SNF2-like ATPase imparts late replication of heterochromatin. Magnetic wire active microrheology of human respiratory mucus Mucus is a viscoelastic gel secreted by the pulmonary epithelium in the tracheobronchial region of the lungs. The coordinated beating of cilia in contact with the gel layer moves mucus upwards towards pharynx, removing inhaled pathogens and particles from the airways. The efficacy of this clearance mechanism depends primarily on the rheological properties of mucus. Here we use a magnetic wire based microrheology technique to study the viscoelastic properties of human mucus collected from human bronchus tubes. The response of wires between 5 and 80 {micro}m in length to a magnetic rotating field is monitored by optical time-lapse microscopy and analyzed using constitutive equation models of rheology, including Maxwell and Kelvin-Voigt. The static shear viscosity and elastic modulus can be inferred from low frequency (10-3 - 10 rad s-1) measurements, leading to the evaluation of the mucin network relaxation time. This relaxation time is found to be widely distributed, from one to several hundred seconds. Mucus is identified as a viscoelastic liquid with an elastic modulus of 2.5 {+/-} 0.5 Pa and a static viscosity of 100 {+/-} 40 Pa s. Our work shows that beyond the established spatial variations in rheological properties due to microcavities, mucus exhibits secondary inhomogeneities associated with the relaxation time of the mucin network that may be important for its flow properties. Overexpression-based detection of translatable circular RNAs is vulnerable to coexistent linear RNA byproducts In RNA field, the demarcation between coding and non-coding has been negotiated by the recent discovery of occasionally translated circular RNAs (circRNAs). Although absent of 5 cap structure, circRNAs can be translated cap-independently. Complementary intron-mediated overexpression is one of the most utilized methodologies for circRNA research but not without bearing echoing skepticism for its poorly defined mechanism and latent coexistent side products. In this study, leveraging such circRNA overexpression system, we have interrogated the protein-coding potential of 30 human circRNAs containing infinite open reading frames in HEK293T cells. Surprisingly, pervasive translation signals are detected by immunoblotting. However, intensive mutagenesis reveals that numerous translation signals are generated independently of circRNA synthesis. We have developed a dual tag strategy to isolate translation noise and directly demonstrate that the fallacious translation signals originate from cryptically spliced linear transcripts. The concomitant linear RNA byproducts, presumably concatemers, can be translated to allow pseudo rolling circle translation signals, and can involve backsplicing junction (BSJ) to disqualify the BSJ-based evidence for circRNA translation. We also find non-AUG start codons may engage in the translation initiation of circRNAs. Taken together, our systematic evaluation sheds light on heterogeneous translational outputs from circRNA overexpression vector and comes with a caveat that ectopic overexpression technique necessitates extremely rigorous control setup in circRNA translation and functional investigation. Cell-type specific circadian bioluminescence rhythms in Dbp reporter mice Circadian rhythms are endogenously generated physiological and molecular rhythms with a cycle length of about 24 h. Bioluminescent reporters have been exceptionally useful for studying circadian rhythms in numerous species. Here, we report development of a reporter mouse generated by modification of a widely expressed and highly rhythmic gene encoding D-site albumin promoter binding protein (Dbp). In this line of mice, firefly luciferase is expressed from the Dbp locus in a Cre-recombinase-dependent manner, allowing assessment of bioluminescence rhythms in specific cellular populations. A mouse line in which luciferase expression was Cre-independent was also generated. The Dbp reporter alleles do not alter Dbp gene expression rhythms in liver or circadian locomotor activity rhythms. In vivo and ex vivo studies show the utility of the reporter alleles for monitoring rhythmicity. Our studies reveal cell-type specific characteristics of rhythms among neuronal populations within the suprachiasmatic nuclei ex vivo. In vivo studies show Dbp-driven bioluminescence rhythms in the liver of Albumin-Cre;DbpKI/+ "liver reporter" mice. After a shift of the lighting schedule, locomotor activity achieved the proper phase relationship with the new lighting cycle more rapidly than hepatic bioluminescence did. As previously shown, restricting food access to the daytime altered the phase of hepatic rhythmicity. Our model allowed assessment of the rate of recovery from misalignment once animals were provided with food ad libitum. These studies confirm the previously demonstrated circadian misalignment following environmental perturbations and reveal the utility of this model for minimally invasive, longitudinal monitoring of rhythmicity from specific mouse tissues. A network approach to genetic circuit designs As genetic circuits become more sophisticated, the size and complexity of data about their designs increases. This data captured goes beyond monolithic genetic sequences and towards circuit modularity and functional details, which are beneficial for analyzing circuit performance and establishing design automation techniques. However, the accessibility, visualisation and usability of design data (and metadata) have received relatively little attention to date. Here, we present a method to turn circuit designs into networks and showcase its potential to enhance the utility of design data. Since networks are dynamic structures, initial graphs can be interactively shaped into sub-networks of relevant information based on requirements such as abstraction, hierarchy and protein interactions. Additionally, several visual changes can be applied, such as colouring or clustering nodes based on types (e.g., genes or promoters), resulting in easier comprehension from a user perspective. This approach allows circuit designs to be coupled to other networks, such as metabolic pathways or implementation protocols captured in graph-like formats. Therefore, we advocate using networks to structure, access and improve synthetic biology information. Distinct SARS-CoV-2 antibody reactivity patterns elicited by natural infection and mRNA Vaccination We analyzed data from two ongoing COVID-19 longitudinal serological surveys in Orange County, CA., between April 2020 and March 2021. A total of 8,476 finger stick blood specimens were collected before and after an aggressive mRNA vaccination campaign. IgG levels were determined using a multiplex antigen microarray containing 10 SARS-CoV-2 antigens, 4 SARS, 3 MERS, 12 Common CoV, and 8 Influenza antigens. Twenty-six percent of 3,347 specimens from unvaccinated Orange County residents in December 2020 were SARS-CoV-2 seropositive. The Ab response was predominantly against nucleocapsid (NP), full length spike and the spike S2 domain. Anti-receptor binding domain (RBD) reactivity was low and there was no cross-reactivity against SARS S1 or SARS RBD. An aggressive mRNA vaccination campaign at the UCI Medical Center started on December 16, 2020 and 6,724 healthcare workers were vaccinated within 3 weeks. Seroprevalence increased from 13% in December to 79% in January, 93% in February and 99% in March. mRNA vaccination induced much higher Ab levels especially against the RBD domain and significant cross-reactivity against SARS RBD and S1 was also observed. Nucleocapsid protein Abs can be used to distinguish individuals in a population of vaccinees to classify those who have been previously infected and those who have not, because nucleocapsid is not in the vaccine. Previously infected individuals developed higher Ab titers to the vaccine than those who have not been previously exposed. These results indicate that mRNA vaccination rapidly induces a much stronger and broader Ab response than SARS-CoV-2 infection. Spontaneous differentiation leads to emergence of hybrid cell states relate to poor prognosis in oral cancer PurposeCellular dynamics between phenotypically heterogeneous subpopulations of cancer cells within individual tumor is shown to be responsible for drug tolerance and overall poor prognosis; however, evidences were largely missing in oral cancer. Therefore, this study was undertaken to describe the dynamic phenotypic states among oral cancer cells, its influence on transcriptomic heterogeneity as well as its clinical significance.

Experimental DesignWe multiplexed phenotypic markers of putative oral-stem-like cancer cells (SLCCs) and characterized diversity among CD44-positive oral cancer cell subpopulations with respect to distinct expression of CD24 and aldehyde dehydrogenase (ALDH)-activity in multiple cell lines. Population trajectories were characterized by Markov model and cell states were defined based on the population specific RNA sequencing (RNAseq). ssGSEA based gene expression signatures were explore for prognostic significance.

ResultsOral cancer cells followed two distinct patterns of spontaneous repopulation dynamics with stochastic inter-conversions on  ALDH-axis, however a strict non-interconvertible transition on  CD24-axis. Interestingly, plastic  ALDH-axis was harnessed to enrich ALDHHigh subpopulations in response to Cisplatin treatment, to adapt a drug tolerant state. Phenotype-specific RNAseq results suggested the organization of subpopulations into hierarchical structure with possible maintenance of intermediate states of stemness within the differentiating oral cancer cells. Further, survival analysis with each subpopulation-specific gene signature strongly suggested that the cell-state dynamics may act as possible mechanism to drive ITH, resulting in poor prognosis in patient.

ConclusionsOur results emphasized the prognostic power of the population dynamics in oral cancer. Importantly, we have described the phenotypic-composition of heterogeneous subpopulations critical for global tumor behaviour in oral cancer; which is a prerequisite knowledge important for precision treatment, however largely lacking for most solid tumors.

Graphical AbstractWe have characterized diversity among CD44-positive oral cancer cells lines with respect to distinct expression of CD24 and ALDH-activity. Subpopulations showed stochastic inter-conversions on ALDH-axis but a strict non-interconvertible transition of CD24Low to CD24High phenotype, even in response to chemotherapy-induced stress. RNAseq study suggested the organization of subpopulations into hierarchical structure with possible maintenance of intermediate alternate states of stemness within the differentiating oral cancer cells. The described population dynamics demonstrtaed influence tumor behaviour possibly by increasing intratumoral heterogeneity in aggressive oral tumors.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=112 SRC="FIGDIR/small/457509v3_ufig1.gif" ALT="Figure 1">
View larger version (20K):
org.highwire.dtl.DTLVardef@1dc6a7org.highwire.dtl.DTLVardef@dd0948org.highwire.dtl.DTLVardef@18c9603org.highwire.dtl.DTLVardef@cca8fb_HPS_FORMAT_FIGEXP  M_FIG C_FIG Translational relevanceIntratumoral heterogeneity (ITH) has been the clinically important factor, impacting aggressive cancer behaviour, drug tolerance and overall poor prognosis. Recent high-throughput studies have provided better cellular and molecular resolution of ITH; however, the prerequisite knowledge which defines the composition of subpopulations critical for global tumor behaviour is majorly lacking for most of the solid tumors. By combining phenotypic markers, we have defined four subpopulations of oral cancer cells. These subpopulations showed stochastic inter-conversions as well as a strict non-interconvertible transition among them to acheive heterogeneity. Importantly, transcriptional states of each subpopulations indicated a clinically relevant signatures for patient prognosis. Also, we observed interconversions of these subpopulations in response to Cisplatin to accumulate drug-tolerant cell state, as rapid and reversible strategy to respond to chemotherapy induced stress. Thus, the characteristics of described phenotypic subgroups may be translated to the clinic for estimating the extent of intratumoral heterogeneity in oral cancer patients. 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. DOMINANT AWN INHIBITOR encodes the ALOG protein originating from gene duplication and inhibits awn elongation by suppressing cell proliferation and elongation in sorghum The awn, a needle-like structure extending from the tip of the lemma in grass species, plays a role in environmental adaptation and fitness. In some crops, awns appear to have been eliminated during domestication. Although numerous genes involved in awn development have been identified, several dominant genes that eliminate awns are also known to exist. For example, in sorghum (Sorghum bicolor), the dominant awn inhibiting gene has been known since 1921; however, its molecular features remain uncharacterized. In this study, we conducted quantitative trait locus analysis and a genome-wide association study of awn-related traits in sorghum and identified DOMINANT AWN INHIBITOR (DAI), which encodes the ALOG family protein on chromosome 3. DAI appeared to be present in most awnless sorghum cultivars, likely because of its effectiveness. Detailed analysis of the ALOG protein family in cereals revealed that DAI originated from duplication of its twin paralog (DAIori) on chromosome 10. Observations of immature awns in near-isogenic lines revealed that DAI inhibits awn elongation by suppressing both cell proliferation and elongation. We also found that only DAI gained a novel function to inhibit awn elongation through an awn-specific expression pattern distinct from that of DAIori. Interestingly, heterologous expression of DAI with its own promoter in rice inhibited awn elongation in the awned cultivar Kasalath. We found that DAI originated from gene duplication, providing an interesting example of gain-of-function that occurs only in sorghum but shares its functionality with rice and sorghum. Mapping the microbial diversity and natural resistome of North Antarctica soils The rising of multiresistant bacterial pathogens is currently one of the most critical threats to global health, demanding a better understanding of the origin and spread of antibiotic resistance. In this regard, the resistome hosted by the microbiota from natural and remote environments remains poorly explored. Moreover, little is known about the availability of antimicrobial resistance genes (ARGs) from these environments to be disseminated through horizontal transfer, potentially mediating the rise of novel resistance factors among clinically relevant pathogens. In this context, the North Antarctica soils are attractive ecosystems to study due to the presence of a microbiota naturally adapted to thrive in harsh conditions, including potential factors to resist natural toxic substances. In this work, we evaluated the antibiotic resistance of bacteria isolated from soils collected in humanized and non-intervened areas of North Antarctica. We identified resistance to a wide array of antibiotics, with isolates harboring up to 10 simultaneous resistances, mainly native Pseudomonas. Genomic analysis revealed the presence of a wide array of genes encoding efflux pumps but the lack of genes explaining some of the resistance phenotypes, suggesting novel uncharacterized mechanisms. Also, using 16S rRNA amplicon and shotgun metagenome sequencing, we explored the microbial diversity in the sampled soils and evaluated the presence of ARGs and their host microbiota. High microbial diversity was found in all the sites, with Proteobacteria, Bacteroidota, Acidobacteriota, and Verrucomicrobiota being the most abundant Phyla, while Candidatus Udaeobacter, RB41, Polaromonas, and Ferruginibacter the most abundant genera. We identified hundreds of genes potentially conferring resistance to more than 15 drug classes, both by short reads analyses and ARG detection among assembled contigs and MAGs obtained combining short and long-read sequence data. Polaromonas, Pseudomonas, Streptomyces, Variovorax, Bhurkolderia, and Gemmatimonas were the main host taxa of the identified ARGs. Part of these ARGs was found inside predicted plasmids, including a putative OXA-like beta-lactamase from Polaromonas harboring the key conserved residues of this kind of enzyme and a conserved predicted protein structure. All this evidence indicates that microbial communities from North Antarctica soil have a highly diverse natural resistome, part of it located inside mobile genetic elements, which would act as a source of novel ARGs. Fatal neuroinvasion and SARS-CoV-2 tropism in K18-hACE2 mice is partially independent on hACE2 expression Animal models recapitulating distinctive features of severe COVID-19 are critical to enhance our understanding of SARS-CoV-2 pathogenesis. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. The precise mechanisms of lethality in this mouse model remain unclear. Here, we evaluated the spatiotemporal dynamics of SARS-CoV-2 infection for up to 14 days post-infection. Despite infection and moderate pneumonia, rapid clinical decline or death of mice was invariably associated with viral neuroinvasion and direct neuronal injury (including brain and spinal neurons). Neuroinvasion was observed as early as 4 dpi, with virus initially restricted to the olfactory bulb supporting axonal transport via the olfactory neuroepithelium as the earliest portal of entry. No evidence of viremia was detected suggesting neuroinvasion occurs independently of entry across the blood brain barrier. SARS-CoV-2 tropism was not restricted to ACE2-expressing cells (e.g., AT1 pneumocytes), and some ACE2-positive lineages were not associated with the presence of viral antigen (e.g., bronchiolar epithelium and brain capillaries). Detectable ACE2 expression was not observed in neurons, supporting overexpression of ACE2 in the nasal passages and neuroepithelium as more likely determinants of neuroinvasion in the K18-hACE2 model. Although our work incites caution in the utility of the K18-hACE2 model to study global aspects of SARS-CoV-2 pathogenesis, it underscores this model as a unique platform for exploring the mechanisms of SARS-CoV-2 neuropathogenesis that may have clinical relevance acknowledging the growing body of evidence that suggests COVID-19 may result in long-standing neurologic consequences.

IMPORTANCECOVID-19 is predominantly a respiratory disease caused by SARS-CoV-2 that has infected more than 191 million people with over 4 million fatalities (2021-07-20). The development of animal models recapitulating distinctive features of severe COVID-19 is critical to enhancing our understanding of SARS-CoV-2 pathogenesis and in the evaluation of vaccine and therapeutic efficacy. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. Here, we show lethality of this model is invariably associated with viral neuroinvasion linked with viral replication and assembly. Importantly, pneumonia albeit invariably present was generally moderate with the absence of culturable infectious virus at peak neuroinvasion. The dynamics of viral neuroinvasion and pneumonia were only partially dependent on hACE2. Overall, this study provides an in-depth sequential characterization of the K18-hACE2 model following SARS-CoV-2 infection, highlighting its significance to further study the mechanisms of SARS-CoV-2 neuropathogenesis. The Dual Mechanisms of Cognitive Control dataset, a theoretically-guided within-subject task fMRI battery Cognitive control is a critical higher mental function, which is subject to considerable individual variation, and is impaired in a range of mental health disorders. We describe here the initial release of Dual Mechanisms of Cognitive Control (DMCC) project data, the DMCC55B dataset, with 55 healthy unrelated young adult participants. Each participant performed four well-established cognitive control tasks (AX-CPT, Cued Task-Switching, Sternberg Working Memory, and Stroop) while undergoing functional MRI scanning. The dataset includes a range of state and trait self-report questionnaires, as well as behavioural tasks assessing individual differences in cognitive ability. The DMCC project is on-going and features additional components (e.g., related participants, manipulations of cognitive control mode, resting state fMRI, longitudinal testing) that will be publicly released following study completion. This DMCC55B subset is released early with the aim of encouraging wider use and greater benefit to the scientific community. The DMCC55B dataset is suitable for benchmarking and methods exploration, as well as analyses of task performance and individual differences.