Integrated quality control of allele-specific copy numbers, mutations and tumour purity from cancer whole genome sequencing assays Cancer genomes contain thousands of somatic point mutations, chromosome copy alterations and more complex structural variants, which contribute to tumour growth and therapy response. Whole genome sequencing is a well established approach for somatic variant identification, but its broad application comes with complications, particularly in how proposed calls are quality assessed. To address this issue, we present CNAqc, a quantitative framework to quality control somatic mutations and allele-specific copy numbers, both in clonal and subclonal settings while accounting for variations in tumour purity, as commonly seen in bulk sampling. We test the model via extensive simulations, validate it using low-pass single-cell data, and apply it to 2778 single-sample PCAWG whole-genomes, 10 in-house multi-region whole-genomes and 48 TCGA whole-exomes. CNAqc is compatible with common bioinformatic pipelines and designed to support automated parameterization processes that are crucial in the era of large-scale whole genome sequencing.
Network Analysis Reveals Different Cellulose Degradation Strategies across Trichoderma harzianum Strains Associated with XYR1 and CRE1 Trichoderma harzianum, whose gene expression is tightly controlled by the transcription factors (TFs) XYR1 and CRE1, is a potential candidate for hydrolytic enzyme production. Here, we performed a network analysis of T. harzianum IOC-3844 and T. harzianum CBMAI-0179 to explore how the regulation of these TFs varies between these strains. In addition, we explored the evolutionary relationships of XYR1 and CRE1 protein sequences among Trichoderma spp. The results of the T. harzianum strains were compared with those of Trichoderma atroviride CBMAI-0020, a mycoparasitic species. Although transcripts encoding carbohydrate-active enzymes (CAZymes), TFs, transporters, and proteins with unknown functions were coexpressed with cre1 or xyr1, other proteins indirectly related to cellulose degradation were identified. The enriched GO terms describing the transcripts of these groups differed across all strains, and several metabolic pathways with high similarity between both regulators but strain-specific differences were identified. In addition, the CRE1 and XYR1 subnetworks presented different topology profiles in each strain, likely indicating differences in the influences of these regulators according to the fungi. The hubs of the cre1 and xyr1 groups included transcripts not yet characterized or described as being related to cellulose degradation. The first-neighbor analyses confirmed the results of the profile of the coexpressed transcripts in cre1 and xyr1. The analyses of the shortest paths revealed that CAZymes upregulated under cellulose degradation conditions are most closely related to both regulators, and new targets between such signaling pathways were discovered. Although the evaluated T. harzianum strains are phylogenetically close and their amino acid sequences related to XYR1 and CRE1 are very similar, the set of transcripts related to xyr1 and cre1 differed, suggesting that each T. harzianum strain used a specific regulation strategy for cellulose degradation. More interestingly, our findings may suggest that XYR1 and CRE1 indirectly regulate genes encoding proteins related to cellulose degradation in the evaluated T. harzianum strains. An improved understanding of the basic biology of fungi during the cellulose degradation process can contribute to the use of their enzymes in several biotechnological applications and pave the way for further studies on the differences across strains of the same species.
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. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

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

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

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



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


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. Neural mechanism underlying task-specific enhancement of motor learning by concurrent transcranial direct current stimulation The optimal protocol for neuromodulation by transcranial direct current stimulation (tDCS) remains unclear. Using rotarod paradigm, we found that mouse motor learning was enhanced by anodal tDCS (3.2 mA/cm2) during but not before or after task performance. Dual-task experiments showed that motor learning enhancement was specific to the task accompanied by concurrent anodal tDCS. Studies using stroke model mice induced by middle cerebral artery occlusion (MCAO) showed that concurrent anodal tDCS restored motor learning capability in a task-specific manner. Transcranial in vivo calcium imaging further showed that anodal and cathodal tDCS elevated and suppressed neuronal activity in the primary motor cortex (M1), respectively. Anodal tDCS specifically promoted the activity of task-related M1 neurons during task performance, suggesting that elevated Hebbian synaptic potentiation in task-activated circuits accounts for motor learning enhancement. Thus, application of tDCS concurrent with the targeted behavioral dysfunction could represent a more effective approach for treating brain disorders. Marine gregarine genomes reveal the breadth of apicomplexan diversity and provide new insights on gliding motility Our current view of the evolutionary history, coding and adaptive capacities of Apicomplexa, protozoan parasites of a wide range of metazoan, is currently strongly biased toward species infecting humans, as data on early diverging apicomplexan lineages infecting invertebrates is extremely limited. Here, we characterized the genome of the marine eugregarine Porospora gigantea, intestinal parasite of Lobsters, remarkable for the macroscopic size of its vegetative feeding forms (trophozoites) and its gliding speed, the fastest so far recorded for Apicomplexa. Two highly syntenic genomes named A and B were assembled. Similar in size (~9 Mb) and coding capacity (~5300 genes), A and B genomes are 10.8% divergent at the nucleotide level, corresponding to 16-38 My in divergent time. Orthogroup analysis across 25 (proto)Apicomplexa species, including Gregarina niphandrodes, showed that A and B are highly divergent from all other known apicomplexan species, revealing an unexpected breadth of diversity. Phylogenetically these two species branch sister to Cephaloidophoroidea, and thus expand the known crustacean gregarine superfamily. The genomes were mined for genes encoding proteins necessary for gliding, a key feature of apicomplexans parasites, currently studied through the molecular model called glideosome. Sequence analysis shows that actin-related proteins and regulatory factors are strongly conserved within apicomplexans. In contrast, the predicted protein sequences of core glideosome proteins and adhesion proteins are highly variable among apicomplexan lineages, especially in gregarines. These results confirm the importance of studying gregarines to widen our biological and evolutionary view of apicomplexan species diversity, and to deepen our understanding of the molecular bases of key functions enabling parasitism, such as the glideosome. A flexible electronic strain sensor for the real-time monitoring of tumor regression Assessing the efficacy of cancer therapeutics in mouse models is a critical step in treatment development. However, low resolution measurement tools and small sample sizes make determining drug efficacy in vivo a difficult and time-intensive task. Here we present a commercially scalable wearable electronic sensor that automates the in vivo testing of cancer therapeutics by continuously monitoring micrometer-scale tumor volume progression or regression in real-time. Histology, caliper measurements, and bioluminescence imaging over a one-week treatment period validated the sensors recordings, but the sensors higher time and length scale resolutions provided the only measurement technique capable of continuously monitoring the immediate hour-scale pharmacodynamic response of a given drug. In mice with two subcutaneously implanted tumor models our sensors recorded significant volume reductions in tumors just 5 hours after small molecule or immunotherapy treatment initiation. We anticipate that real-time tumor regression datasets could help expedite and automate the process of screening cancer therapies in vivo. What can cold-induced transcriptomes of Arctic Brassicaceae tell us about the evolution of cold tolerance? O_LILittle is known about the evolution of cold tolerance in polar plant species and how they differ from their temperate relatives. To gain insight into their biology and the evolution of cold tolerance, we compared the molecular basis of cold response in three Arctic Brassicaceae species.
C_LIO_LIWe conducted a comparative time series experiment to examine transcriptional responses to low temperature. RNA was sampled at 22 {degrees}C, and after 3h, 6h, and 24h at 2 {degrees}C. We then identified sets of genes that were differentially expressed in response to cold and compared them between species, as well as to published data from the temperate Arabidopsis thaliana.
C_LIO_LIMost differentially expressed genes were species-specific, but a significant portion of the cold response was also shared among species. Among thousands of differentially expressed genes, [~]200 were shared among the three Arctic species and A. thaliana, while [~]100 were exclusively shared among the three Arctic species.
C_LIO_LIOur results show that cold response differs markedly between Arctic Brassicaceae species, but likely builds on a conserved basis found across the family. They also confirm that highly polygenic traits such as cold tolerance may show little repeatability in their patterns of adaptation.
C_LI Naming of human diseases on the wrong side of history BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. The evolution of stomatal traits along the trajectory towards C4 photosynthesis C4 photosynthesis optimizes plants carbon and water relations, allowing high photosynthetic rate with low stomatal conductance. Stomata have long been believed as a part of C4 syndrome. However, it remains unclear how stomata traits evolved along the path from C3 to C4. Stomatal patterning was examined in Flaveria genus, a model for studying C4 evolution. Comparative, transgenic and semi-in-vitro experiments were used to study molecular basis that underlies stomatal traits along C4 evolution. Novel results: the evolution from C3 to C4 species through intermediate species is accompanied by a stepwise rather than an abrupt change in the stomatal traits. The initial change occurs near Type II and dramatic change occurs at the C4-like species. On the road to C4, stomata become less in number but bigger in size and changes in stomatal density dominates changes in maximum stomatal conductance (gsmax). The reduction of FSTOMAGEN expression underlies altered gsmax between Flaveria species with different photosynthetic pathways and likely occurs in other C4 lineages. Our study provides insight into the pattern, mechanism and role of stomatal evolution along the road towards C4. This work highlights the stomatal traits in the current C4 evolutionary model and the co-evolution of photosynthetic pathway and stomata. Production of Offspring from Azoospermic Mice with Meiotic Failure: Precise Biparental Meiosis within Halved Oocytes While the large volume of mammalian oocytes is necessary for embryo development, it can lead to error-prone chromosomal segregation during meiosis. Conversely, we hypothesized that smaller oocytes would have a great unidentified potential to stabilize unstable meiosis and improve the development of the resultant embryos. Here, we show that reducing ooplasmic volume can rescue highly error-prone fertilization using primary spermatocytes by preventing segregation errors of chromosomes during biparental meiosis. High-resolution live-imaging analysis revealed that erroneous chromosome segregation occurred in most (90%) spermatocyte-injected oocytes of normal size, but could be ameliorated to 40% in halved oocytes. The birth rate improved remarkably from 1% to 19% (P < 0.0001). Importantly, this technique enabled the production of offspring from azoospermic mice with spermatocyte arrest caused by STX2 deficiency, an azoospermia factor also found in humans. Thus, contrary to popular opinion, oocytes inherently have a strong potential for precise meiotic divisions, which can be evoked by reduction of the ooplasmic volume. Their potential might help rescue cases of untreatable human azoospermia with spermatocyte arrest. Novel Apelin-expressing gCap Endothelial Stem-like Cells Orchestrate Lung Microvascular Repair QuestionWe sought to define the mechanism underlying lung microvascular regeneration in a severe acute lung injury (ALI) model induced by selective lung endothelial cell ablation.

MethodsChanges in lung cell populations and gene expression profiles were determined in transgenic mice expressing human diphtheria toxin (DT) receptor targeted to ECs using single-cell RNA sequencing at baseline (day 0) and days 3, 5 and 7 after lung EC ablation.

ResultsEight distinct endothelial clusters were resolved, including alveolar aerocytes (aCap) ECs expressing apelin at baseline, and general capillary (gCap) ECs expressing the apelin receptor. Intratracheal instillation of DT resulted in ablation of >70% of lung ECs, producing severe ALI with near complete resolution by 7 days. At 3 days post injury, a novel gCap population emerged characterized by de novo expression of apelin, together with the stem cell marker, protein C receptor. These stem-like cells transitioned to proliferative ECs, expressing apelin receptor together with the pro-proliferative transcription factor, FoxM1. This progenitor-like cell population was responsible for the rapid replenishment of all depleted EC populations by 7 days post injury, including aerocytes which play a critical role in re-establishment of the air-blood barrier. Treatment with an apelin receptor antagonist prevented recovery and resulted in excessive mortality, consistent with a central role for apelin signaling in EC regeneration and microvascular repair.

ConclusionThe lung has a remarkable capacity for microvasculature EC regeneration which is orchestrated by signaling between newly emergent apelin-expressing gCap endothelial stem-like cells and highly proliferative, apelin receptor positive endothelial progenitors.

Take-Home messageUsing sublethal lung endothelial cell (EC) ablation, we show for the first that EC regeneration and resolution of acute lung injury is orchestrated by novel apelin-expressing, gCap endothelial stem-like cells by a mechanism requiring apelin signaling.

Graphical Abstract O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY

A schematic representation of EC populations contributing to microvascular repair. At baseline (Day 0), there are two main alveolar groups of capillary ECs: larger apelin positive aCap ECs, termed aerocytes, that play a key structural role in forming the air-blood barrier; and smaller apelin receptor (Aplnr) expressing gCap ECs, which are found in the thicker regions at the corners of the alveoli. After DT-induced EC ablation, there is a marked depletion of both EC populations and the appearance of novel transitional and transient populations. At Day 3, there is the appearance of stem-like gCap ECs that paradoxically express apelin, but not its receptor, and are characterized by various stem and progenitor cell markers but show no evidence of proliferation. By Day 5, these transition to ECs expressing Aplnr which have a strong proliferative phenotype, as evidenced by FoxM1 and Ki67 expression, and then rapidly replenish depleted EC pools, including aCap ECs, by Day 7. This transition is orchestrated by the interaction of apelin with its receptor as a critical mechanism in lung microvascular regeneration after EC injury. AT1 = alveolar type -1 epithelial cell; AT2 = alveolar type-2 epithelial cell; APLNR = apelin receptor; ANGPT2 = angiopoietin 2; EPCR = Endothelial protein C receptor. Task-independent metrics of computational hardness predict human cognitive performance The survival of human organisms depends on our ability to solve complex tasks in the face of limited cognitive resources. However, little is known about the factors that drive the complexity of those tasks. Here, building on insights from computational complexity theory, we quantify the computational hardness of cognitive tasks using a set of task-independent metrics related to the computational resource requirements of individual instances of a task. We then examine the relation between those metrics and human behavior and find that they predict both time spent on a task as well as accuracy in three canonical cognitive tasks. Our findings demonstrate that performance in cognitive tasks can be predicted based on generic metrics of their inherent computational hardness.

TeaserThe ability of humans to solve cognitive tasks is affected by generic mathematical properties of problems related to their computational complexity. Reconstructing SNP Allele and Genotype Frequencies from GWAS Summary Statistics The emergence of genomewide association studies (GWAS) has led to the creation of large repositories of human genetic variation, creating enormous opportunities for genetic research and worldwide collaboration. Methods that are based on GWAS summary statistics seek to leverage such records, overcoming barriers that often exist in individual-level data access while also offering significant computational savings. Here, we propose a novel framework that can reconstruct allelic and genotypic counts/frequencies for each SNP from case-control GWAS summary statistics. Our framework is simple and efficient without the need of any complicated underlying assumptions. Illustrating the great potential of this framework we also propose three summary-statistics-based applications implemented in a new software package (ReACt): GWAS meta-analysis (with and without sample overlap), case-case GWAS, and, for the first time, group polygenic risk score (PRS) estimation. We evaluate our methods against the current state-of-the-art on both synthetic data and real genotype data and show high performance in power and error control. Our novel group PRS method based on summary statistics could not be achieved prior to our proposed framework. We demonstrate here the potential applications and advantages of this approach. Our work further highlights the great potential of summary-statistics-based methodologies towards elucidating the genetic background of complex disease and opens up new avenues for research. Chimeric mRNA based COVID-19 vaccine induces protective immunity against Omicron and Delta The emerging SARS-CoV-2 variants of concern (VOCs) exhibit enhanced transmission and immune escape, reducing the efficacy and effectiveness of the two FDA-approved mRNA vaccines. Here, we explored various strategies to develop novel mRNAs vaccines to achieve safer and wider coverage of VOCs. Firstly, we constructed a cohort of mRNAs that feature a furin cleavage mutation in the spike (S) protein of predominant VOCs, including Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1) and Delta (B.1.617.2). Not present in the mRNA vaccines currently in use, the mutation abolished the cleavage between the S1 and S2 subunits, potentially enhancing the safety profile of the immunogen. Secondly, we systematically evaluated the induction of neutralizing antibodies (nAb) in vaccinated mice, and discovered that individual VOC mRNAs elicited strong neutralizing activity in a VOC-specific manner. Thirdly, the IgG produced in mice immunized with Beta-Furin and Washington (WA)-Furin mRNAs showed potent cross-reactivity with other VOCs, which was further corroborated by challenging vaccinated mice with the live virus of VOCs. However, neither WA-Furin nor Beta-Furin mRNA elicited strong neutralizing activity against the Omicron variant. Hence, we further developed an Omicron-specific mRNA vaccine that restored protection against the original and the sublineages of Omicron variant. Finally, to broaden the protection spectrum of the new Omicron mRNA vaccine, we tested the concept of bivalent immunogen. Instead of just fusing two RBDs head-to-tail, we for the first time constructed an mRNA-based chimeric immunogen by introducing the RBD of Delta variant into the entire S antigen of Omicron. The resultant chimeric mRNA was capable of inducing potent and broadly acting nAb against Omicron (both BA.1 and BA.2) and Delta, which paves the way to develop new vaccine candidate to target emerging variants in the future. 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. Betacoronavirus binding dynamics relevant to the functional evolution of the highly transmissible SARS-CoV-2 variant N501Y Comparative functional analysis of the dynamic interactions between various Betacoronavirus mutant strains and broadly utilized target proteins such as ACE2 and CD26, is crucial for a more complete understanding of zoonotic spillovers of viruses that cause diseases such as COVID-19. Here, we employ machine learning to replicated sets of nanosecond scale GPU accelerated molecular dynamics simulations to statistically compare and classify atom motions of these target proteins in both the presence and absence of different endemic and emergent strains of the viral receptor binding domain (RBD) of the S spike glycoprotein. Machine learning was used to identify functional binding dynamics that are evolutionarily conserved from bat CoV-HKU4 to human endemic/emergent strains. Conserved dynamics regions of ACE2 involve both the N-terminal helices, as well as a region of more transient dynamics encompassing K353, Q325 and a novel motif AAQPFLL 386-92 that appears to coordinate their dynamic interactions with the viral RBD at N501. We also demonstrate that the functional evolution of Betacoronavirus zoonotic spillovers involving ACE2 interaction dynamics are likely pre-adapted from two precise and stable binding sites involving the viral bat progenitor strains interaction with CD26 at SAMLI 291-5 and SS 333-334. Our analyses further indicate that the human endemic strains hCoV-HKU1 and hCoV-OC43 have evolved more stable N-terminal helix interactions through enhancement of an interfacing loop region on the viral RBD, whereas the highly transmissible SARS-CoV-2 variants (B.1.1.7, B.1.351 and P.1) have evolved more stable viral binding via more focused interactions between the viral N501 and ACE2 K353 alone. Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. 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. Intermittent compression increases flow through a gel: a mechanical model perhaps relevant to flow of cerebrospinal fluid Exchange of molecules between cerebrospinal fluid (CSF) and brain cells contributes to brain function and protection from dementia. Despite widespread acceptance of the  glymphatic theory, the route by which CSF is brought close enough to the neural tissue for solutes to be exchanged by extracellular diffusion is not entirely clear. Exogenous molecules injected into CSF are observed to reach the basement lamina that surrounds the dense capillary network. Transport of solutes by diffusion along the basement lamina, a gel of macromolcules about 100 nm thick, would be too slow; bulk flow in a static geometry would require unphysiologically high pressures. However, it is known that the pulsation of blood aids transport of CSF, and we hypothesized that this is because the pulsation intermittently squeezes the pericapillary lamina. In a primitive mimicry, we have tested whether intermittent squeezing increases flow through an agar gel. In all but one of 216 tests, pulsation caused a reversible increase, sometimes by a factor of 100 or more. The enhancement was greatest for frequencies 5-11 Hz and, over the tested range of pressure heads (20 - 50 cmH2O), was greatest for the lowest pressure. The results suggest a reason why exercise slows the aging of the brain. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. Usutu virus escapes langerin-induced restriction to productively infect human Langerhans cells, unlike West Nile virus. Usutu virus (USUV) and West Nile virus (WNV) are phylogenetically close emerging arboviruses transmitted by mosquitoes, and constitute a global public health threat. Since USUV and WNV enter the body through the skin, the first immune cells they encounter are skin-resident dendritic cells, the most peripheral outpost of immune defense. This unique network is composed of Langerhans cells (LCs) and dermal DCs, which reside in the epidermis and the dermis, respectively.

Using human skin explants, we show that while both viruses can replicate in keratinocytes, they can also infect resident DCs with distinct tropism, since WNV preferentially infects DCs in the dermis, whereas USUV has a greater propensity to infect LCs. Using both purified human epidermal LCs (eLCs) and monocyte derived LCs (MoLCs), we confirm that LCs sustain a faster and more efficient replication of USUV compared with WNV and that this correlates with a more intense innate immune response to USUV compared with WNV. Next, we show that ectopic expression of the LC-specific C-type lectin receptor (CLR), langerin, in HEK293T cells allows WNV and USUV to bind and enter, but supports the subsequent replication of USUV only. Conversely, blocking or silencing langerin in MoLCs or eLCs made them resistant to USUV infection, thus demonstrating that USUV uses langerin to enter and replicate in LCs. Altogether, our results demonstrate that LCs constitute privileged target cells for USUV in human skin, because langerin favors its entry and replication. Intriguingly, this suggests that USUV efficiently escapes the antiviral functions of langerin, which normally safeguards LCs from most viral infections. Limited selection towards lighter seeds in Danish grasslands over an eight-year period Time-series analysis of cover data from 236 Danish grassland sites demonstrated community selection on seed mass. Across four grassland habitat types and during an eight-year period there was significant community selection towards lighter seeds. In the analysis, the continuous seed size variable was used to group plant species into functional types. This method was chosen in order to account for the sampling process of the abundance data, and is in contrast to most other analyses of trait selection, where the community weighted mean of the traits is used as the dependent variable. Generally, little is known on the underlying causes of possible plant trait community selection in semi-natural grasslands that are not subjected to successional processes. DNA-controlled Spatio-temporal Patterning of a Cytoskeletal Active Gel Living cells move and change their shape because signaling chemical reactions modify the state of their cytoskeleton; an active gel that converts chemical energy into mechanical forces. To create life-like materials, it is thus key to engineer chemical pathways that drive active gels. Here, we describe the preparation of DNA-responsive surfaces that control the activity of a cytoskeletal active gel comprised of microtubules: a DNA signal triggers the release of molecular motors from the surface into the gel bulk, generating forces that structure the gel. Depending on the DNA sequence and concentration, the gel forms a periodic band pattern or contracts globally. Finally, we show that the structuration of the active gel can be spatially controlled in the presence of a gradient of DNA concentration. We anticipate that such DNA-controlled active matter will contribute to the development of life-like materials with self-shaping properties. Structure-Functional-Selectivity Relationship Studies on A-86929 Analogs and Small biaryl Fragments toward Discovery of biased D1 agonists Dopamine regulates normal functions such as movement, reinforcement learning, and cognition, and its dysfunction has been implicated in multiple psychiatric and neurological disorders. Dopamine acts through the D1- (D1R and D5R) and D2-class (D2R, D3R and D4R) of seven transmembrane receptors, and activates both G-protein- and {beta}-arrestin-dependent signaling pathways, to mediate its physiological effects. Current dopamine receptor-based therapies are used to ameliorate motor deficits in Parkinsons disease, or as antipsychotic medications for schizophrenia. These drugs show efficacy for ameliorating only some symptoms caused by dopamine dysfunction and are plagued by debilitating side-effects. Studies in primates and rodents have shown that shifting the balance of dopamine receptor signaling towards the arrestin pathway can be beneficial for inducing normal movement, while reducing motor side-effects such as dyskinesias, and can be efficacious at enhancing cognitive function compared to balanced agonists. Several structure-activity-relationship (SAR) studies have embarked on discovering {beta}-arrestin-biased dopamine agonists, focused on D2 partial agonists, non-catechol D1 agonists, and mixed D1/D2R dopamine receptor agonists. Here, we describe an SAR study to identify novel D1R {beta}-arrestin biased ligands using A-86929, a high-affinity D1R catechol agonist, as a core scaffold. Previously described and novel analogs of A-86929 were synthesized and screened in vitro for structure-functional-selectivity relationships (SFSR) studies to identify chemical motifs responsible for {beta}-arrestin biased activity at both D1 and D2Rs. Most of the A-86929 analogs screened were G protein biased but none of them were exclusively arrestin-biased. Additionally, various catechol aryl fragments were designed and synthesized. Other compounds surveyed included hydroxyl and chloro analogs of dopamine to test for hydrogen bonding and ionic interactions. Some of these small molecular probes displayed weak bias towards the {beta}-arrestin pathway. Continued in-depth SFSR studies informed by structure determination, molecular modeling, and mutagenesis studies will facilitate discovery of potent and efficacious arrestin-biased dopamine receptor ligands. Poly ADP-ribosylation of SET8 leads to aberrant H4K20 methylation domains in mammalian cells In mammalian cells, SET8 mediated Histone H4 Lys 20 monomethylation (H4K20me1) has been implicated in regulating mitotic condensation, DNA replication, DNA damage response, and gene expression. Here we show SET8, the only known enzyme for H4K20me1 is post-translationally poly ADP-ribosylated by PARP1 on lysine residues. PARP1 interacts with SET8 in a cell cycle-dependent manner. Poly ADP-ribosylation on SET8 renders it catalytically compromised, and degradation via ubiquitylation pathway. Knockdown of PARP1 led to an increase of SET8 protein levels, leading to aberrant H4K20me1 and H4K20me3 domains genome-wide. H4K20me1 is associated with higher gene transcription levels while the increase of H4K20me3 levels was predominant in DNA repeat elements correlating with modulation of chromatin states. Hence, SET8 mediated chromatin remodeling in mammalian cells are modulated by poly ADP-ribosylation by PARP1.