Experimental evolution reveals post-transcriptional regulation as a novel driver of Leishmania fitness gain The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network driving parasite fitness, with genome instability causing highly reproducible, gene dosage-dependent changes in protein abundance linked to post-transcriptional regulation. These in turn were associated with a gene dosage-independent reduction in abundance of flagellar transcripts and a coordinated increase in abundance of coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of fitness-adapted ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings.
A diagnostic primer pair to distinguish between wMel and wAlbB Wolbachia infection - A simplified real-time PCR approach for use in field collected Aedes aegypti in dengue reduction trials Detection of the Wolbachia endosymbiont in Aedes aegypti mosquitoes through real-time PCR assays is widely used during and after Wolbachia releases in dengue reduction trials involving the wMel and wAlbB strains. However, primers applied in current successful Wolbachia releases cannot be used in a single assay to distinguish between these strains. Here, we developed a new diagnostic primer pair wMA which can detect wMel or wAlbB infection in the same assay. We also tested current Wolbachia primers and show that there is variation in their performance when they are used to assess the relative density of Wolbachia. The new wMA primers provide an accurate and efficient estimate of the presence and density of both Wolbachia infections.
An alternatively spliced TREM2 isoform lacking the ligand binding domain is expressed in human brain BackgroundGenetic variants in TREM2 are strongly associated with Alzheimers Disease (AD) risk but alternative splicing in TREM2 transcripts has not been comprehensively described.

ObjectiveRecognizing that alternative splice variants can result in reduced gene expression and/or altered function, we sought to fully characterize splice variation in TREM2.

MethodsHuman blood and anterior cingulate autopsy tissue from 61 donors were used for end-point and quantitative PCR and Western blotting to identify and quantify novel TREM2 isoforms.

ResultsIn addition to previously described transcripts lacking exon 3 or exon 4, or retaining part of intron 3, we identified novel isoforms lacking exon 2, along with isoforms lacking multiple exons. Isoforms lacking exon 2 were predominant at approximately 10% of TREM2 mRNA in the brain. Expression of TREM2 and frequency of exon 2 skipping did not differ between AD samples and non-AD controls (p = 0.1268 and p = 0.4909, respectively). Further, these novel splice isoforms were also observed across multiple tissues with similar frequency (range 5.3 - 13.0%). We found that the exon 2 skipped isoform D2-TREM2 is translated to protein and localizes similarly to full-length TREM2 protein, that both proteins are primarily retained in the Golgi complex, and that D2-TREM2 is expressed in AD and non-AD brain.

ConclusionSince the TREM2 ligand binding domain is encoded by exon 2, and skipping this exon retains reading frame while conserving localization, we hypothesize that D2-TREM2 acts as an inhibitor of TREM2 and targeting TREM2 splicing may be a novel therapeutic pathway for AD. STING activation promotes robust immune response and NK cell-mediated tumor regression in glioblastoma models Immunotherapy has had a tremendous impact on cancer treatment in the past decade, with hitherto unseen responses at advanced and metastatic stages of the disease. However, the aggressive brain tumor glioblastoma (GBM) is highly immunosuppressive and remains largely refractory to current immunotherapeutic approaches. The cGAS-STING cytoplasmic double stranded DNA (dsDNA) sensing pathway has emerged as a next-generation immunotherapy target with potent local immune stimulatory properties.

Here, we investigated the status of the STING pathway in GBM and the modulation of the brain tumor microenvironment (TME) with the STING agonist ADU-S100. Our data reveal the presence of STING in human GBM specimens, where it stains strongly in the tumor vasculature. We show that human GBM explants can respond to STING agonist treatment by secretion of inflammatory cytokines. In murine GBM models, we show a profound shift in the tumor immune landscape after STING agonist treatment, with massive infiltration of the tumor-bearing hemisphere with innate immune cells including inflammatory macrophages, neutrophils and NK populations. Treatment of established murine intracranial GL261 and CT-2A tumors by biodegradable ADU-S100-loaded intracranial implants demonstrated a significant increase in survival in both models and long-term survival with immune memory in GL261. Responses to treatment were abolished by NK cell depletion. This study reveals therapeutic potential and deep remodeling of the TME by STING activation in GBM and warrants the further examination of STING agonists alone or in combination with other immunotherapies such as cancer vaccines, CAR T cells, NK therapies or immune checkpoint blockade.

Significance statementModulation of the immune microenvironment is critical for immunosuppressive and therapy refractory tumors like glioblastoma. Activation of the STING pathway deeply remodels the brain tumor environment and attracts innate immune cells and natural killer cell populations, producing a robust antitumor effect with long-term immune memory. We further show that human glioblastoma tissue can respond to the therapy and lay the foundations for combined intracranial immunotherapies by using crosslinked biodegradable brain implants. Social behavioral profiling by unsupervised deep learning reveals a stimulative effect of dopamine D3 agonists on zebrafish sociality It has been a major challenge to systematically evaluate and compare how pharmacological perturbations influence social behavioral outcomes. Although some pharmacological agents are known to alter social behavior, precise description and quantification of such effects have proven difficult. The complexity of brain functions regulating sociality makes it challenging to predict drug effects on social behavior without testing in live animals, and most existing behavioral assays are low-throughput and provide only unidimensional readouts of social function. To achieve richer characterization of drug effects on sociality, we developed a scalable social behavioral assay for zebrafish named ZeChat based on unsupervised deep learning. High-dimensional and dynamic social behavioral phenotypes are automatically classified using this method. By screening a neuroactive compound library, we found that different classes of chemicals evoke distinct patterns of social behavioral fingerprints. By examining these patterns, we discovered that dopamine D3 agonists possess a social stimulative effect on zebrafish. The D3 agonists pramipexole, piribedil, and 7-hydroxy-DPAT-HBr rescued social deficits in a valproic acid-induced zebrafish autism model. The ZeChat platform provides a promising approach for dissecting the pharmacology of social behavior and discovering novel social-modulatory compounds. The compound of baicalein, wogonin, and oroxylin-A inhibiting EMT in A549 cell line via PI3K/AKT-TWIST1-glycolysis pathway Non-small cell lung cancer (NSCLC) is a worldwide disease with high morbidity and mortality, which is most derived from its metastasis. Some researches show that the epithelial-mesenchymal transition (EMT) process promotes lung cancer cells migration and invasion, leading to NSCLC metastasis. Total Flavonoid Aglycones Extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds (baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)), the reconstructed TFAE (reTFAE) present inhibitory effect on the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC. Alpha oscillations shape sensory representation and perceptual sensitivity Alpha activity (8-14 Hz) is the dominant rhythm in the awake brain, and thought to play an important role in setting the brains internal state. Previous work has associated states of decreased alpha power with enhanced neural excitability. However, evidence is mixed on whether and how such excitability enhancement modulates sensory signals of interest versus noise differently, and what, if any, the consequences are for subsequent perception. Here, human subjects (male and female) performed a visual detection task in which we manipulated their decision criteria in a block-wise manner. While our manipulation led to substantial criterion shifts, these shifts were not reflected in pre-stimulus alpha-band changes. Rather, lower pre-stimulus alpha power in occipital-parietal areas improved perceptual sensitivity and enhanced information content decodable from neural activity patterns. Additionally, oscillatory alpha phase immediately before stimulus presentation modulated accuracy. Together, our results suggest that alpha-band dynamics modulate sensory signals of interest more strongly than noise.

Significance statementThe internal state of our brain fluctuates, giving rise to variability in perception and action. Neural oscillations, most prominently in the alpha-band, have been suggested to play a role in setting this internal state. Here, we show that ongoing alpha-band activity in occipital-parietal regions predicts the quality of visual information decodable in neural activity patterns, and subsequently human observers sensitivity in a visual detection task. Our results provide comprehensive evidence that visual representation is modulated by ongoing alpha-band activity, and advance our understanding on how, when faced with unchanging external stimuli, internal neural fluctuations influence perception and behavior. Naming Human Diseases: Ethical Principles of Curating Exclusive Substitute for Inopportune Nosology 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. Sablefish (Anoplopoma fimbria Pallas, 1814) plasma biochemistry and hematology reference intervals including blood cell morphology Plasma biochemistry and hematology reference intervals are integral health assessment tools in all medical fields, including aquatic animal health. As sablefish (Anoplopoma fimbria) are becoming aquaculturally and economically more important, this manuscript provides essential reference intervals (RI) for their plasma biochemistry and hematology along with reference photomicrographs of blood cells in healthy, fasted sablefish. Blood cell morphology can differ between fish species. In addition, blood cell counts and blood chemistry can vary between fish species, demographics, water conditions, seasons, diets, and culture systems, which precludes the use of RIs from other fish species. For this study, blood was collected for plasma biochemistry and hematology analysis between June 20 and July 18, 2019, from healthy, yearling sablefish, hatched and reared in captivity on a commercial diet. Overnight fast of 16-18 hours did not sufficiently reduce lipids in the blood, which led to visible lipemia and frequent rupture of blood cells during analysis. Therefore, sablefish should be fasted for 24 to 36 hours before blood is collected to reduce hematology artifacts or possible reagent interference in plasma biochemistry analysis. Lymphocytes were the most dominant leukocytes (98%), while eosinophils were rare, and basophils were not detected in sablefish. Neutrophils were very large cells with Dohle bodies. In mammals and avian species, Dohle bodies are usually signs of toxic change from inflammation, but no such association was found in these fish. In conclusion, lipemia can interfere with sablefish blood analysis, and available removal methods should be evaluated as fasting for up to 36 h might not always be feasible. Also, more studies are required to establish RI for different developmental stages and rearing conditions. 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. Collective polymerase dynamics emerge from DNA supercoiling during transcription All biological processes ultimately come from physical interactions. The mechanical properties of DNA play a critical role in transcription. RNA polymerase can over or under twist DNA (referred to as DNA supercoiling) when it moves along a gene resulting in mechanical stresses in DNA that impact its own motion and that of other polymerases. For example, when enough supercoiling accumulates, an isolated polymerase halts and transcription stops. DNA supercoiling can also mediate non-local interactions between polymerases that shape gene expression fluctuations. Here, we construct a comprehensive model of transcription that captures how RNA polymerase motion changes the degree of DNA supercoiling which in turn feeds back into the rate at which polymerases are recruited and move along the DNA. Surprisingly, our model predicts that a group of three or more polymerases move together at a constant velocity and sustain their motion (forming what we call a polymeton) whereas one or two polymerases would have halted. We further show that accounting for the impact of DNA supercoiling on both RNA polymerase recruitment and velocity recapitulates empirical observations of gene expression fluctuations. Finally, we propose a mechanical toggle switch whereby interactions between genes are mediated by DNA twisting as opposed to proteins. Understanding the mechanical regulation of gene expression provides new insights into how endogenous genes can interact and informs the design of new forms of engineered interactions.

PACS numbers: Mechanical bistability of the mesoderm epithelium facilitates mesoderm invagination during Drosophila gastrulation Apical constriction driven by non-muscle myosin II ("myosin") provides a well-conserved mechanism to mediate epithelial folding. It remains unclear how contractile forces near the apical surface of a cell sheet drive out-of-the-plane bending of the sheet and whether myosin contractility is required throughout folding. By optogenetic-mediated acute inhibition of myosin, we find that during Drosophila mesoderm invagination, myosin contractility is critical to prevent tissue relaxation during the early, "priming" stage of folding but is dispensable for the actual folding step after the tissue passes through a stereotyped transitional configuration. The binary response suggests that the mesoderm is mechanically bistable during gastrulation. Combined modeling analysis and experimental measurements suggest that the observed mechanical bistability may arise from apicobasal shrinkage of the surrounding ectoderm, which promotes mesoderm invagination by facilitating a buckling transition. Our results suggest that Drosophila mesoderm invagination requires a joint action of local myosin contractility and mechanical bistability of the epithelium to trigger epithelial buckling. Assigning metabolic rate measurements to torpor and euthermy in heterothermic endotherms: 'torpor', a new package for R. Torpor is a state of controlled reduction of metabolic rate (M) in endotherms. Assigning measurements of M to torpor or euthermy can be challenging, especially when the difference between euthermic M and torpid M is small, in species defending a high minimal body temperature in torpor, in thermolabile species, and slightly below the thermoneutral zone (TNZ). Here, we propose a novel method for distinguishing torpor from euthermy. We use the variation in M measured during euthermic rest and torpor at varying ambient temperatures (Ta) to objectively estimate the lower critical temperature (Tlc) of the TNZ and to assign measurements to torpor, euthermic rest or rest within TNZ. In addition, this method allows the prediction of M during euthermic rest and torpor at varying Ta, including resting M within the TNZ. The present method has shown highly satisfactory results using 28 published sets of metabolic data obtained by respirometry on 26 species of mammals. Ultimately, this novel method aims to facilitate analysis of respirometry data in heterothermic endotherms. Finally, the development of the associated R-package (torpor) will enable widespread use of the method amongst biologists.

Summary statementThe presented method and its associated R-package (torpor) enable the assignment of metabolic rate measurements to torpor or euthermy, ultimately improving the standardization of respirometry analyses in heterotherms. Comparative proteome signatures of trace samples by multiplexed Data-Independent Acquisition Single cell transcriptomics has revolutionized our understanding of basic biology and disease. Since transcript levels often do not correlate with protein expression, it is crucial to complement transcriptomics approaches with proteome analyses at single cell resolution. Despite continuous technological improvements in sensitivity, mass spectrometry-based single cell proteomics ultimately faces the challenge of reproducibly comparing the protein expression profiles of thousands of individual cells. Here, we combine two hitherto opposing analytical strategies, DIA and Tandem-Mass-Tag (TMT)-multiplexing, to generate highly reproducible, quantitative proteome signatures from ultra-low input samples. While conventional, data-dependent shotgun proteomics (DDA) of ultra-low input samples critically suffers from the accumulation of missing values with increasing sample-cohort size, data-independent acquisition (DIA) strategies do usually not take full advantage of isotope-encoded sample multiplexing. We developed a novel, identification-independent proteomics data-analysis pipeline that allows to quantitatively compare DIA-TMT proteome signatures across hundreds of samples independent of their biological origin, and to identify cell types and single protein knockouts. We validate our approach using integrative data analysis of different human cell lines and standard database searches for knockouts of defined proteins. These data establish a novel and reproducible approach to markedly expand the numbers of proteins one detects from ultra-low input samples, such as single cells. Sharpened and mechanically durable carbon fiber electrode arrays for neural interfacing Bioelectric medicine treatments target disorders of the nervous system unresponsive to pharmacological methods. While current stimulation paradigms effectively treat many disorders, the underlying mechanisms are relatively unknown, and current neuroscience recording electrodes are often limited in their specificity to gross averages across many neurons or axons. Here, we develop a novel, durable carbon fiber electrode array adaptable to many neural structures for precise neural recording. Carbon fibers were sharpened using a reproducible blowtorch method that uses the reflection of fibers against the surface of a water bath. The arrays were developed by partially embedding carbon fibers in medical-grade silicone to improve durability. We recorded acute spontaneous electrophysiology from the rat cervical vagus nerve (CVN), feline dorsal root ganglia (DRG), and rat brain. Blowtorching resulted in fibers of 72.3 {+/-} 33.5 degree tip angle with 146.8 {+/-} 17.7 m exposed carbon. Silicone-embedded carbon fiber arrays were durable with 87.5% of fibers remaining after 50,000 passes. Observable neural clusters were recorded using sharpened carbon fiber electrodes from rat CVN (41.8 Vpp), feline DRG (101.1 Vpp), and rat brain (80.7 Vpp). Recordings from the feline DRG included physiologically-relevant signals from increased bladder pressure and cutaneous brushing. These results suggest that this carbon fiber array is a uniquely durable and adaptable neural recording device. In the future, this device may be useful as a bioelectric medicine tool for diagnosis and closed-loop neural control of therapeutic treatments and monitoring systems. Regulatory perturbations of ribosome allocation reshape the growth proteome with a trade-off in adaptation capacity Bacteria regulate their cellular resource allocation to enable fast growth-adaptation to a variety of environmental niches. We studied the ribosomal allocation, growth and expression profiles of two sets of fast-growing mutants of Escherichia coli K-12 MG1655 in glucose minimal medium. Mutants with only 3 of the seven copies of ribosomal RNA operons grew faster than the wild-type strain in minimal media and show similar phenotype to previously studied rpoB mutants. Higher growth rates due to increased ribosome content affected resource allocation. Expression profiles of fast-growing mutants shared downregulation of hedging functions and upregulated growth functions. Mutants showed longer diauxic shifts and reduced activity of gluconeogenic promoters during glucose-acetate shifts, suggesting reduced availability of the RNA Polymerase for expressing hedging proteome. These results show that the regulation of ribosomal allocation underlies the growth/hedging phenotypes obtained from laboratory evolution experiments. We show how two different regulatory perturbations (rRNA promoters or rpoB mutations) reshape the proteome for growth with a concomitant fitness cost

HighlightsMutants with only 3 ribosomal operons grow faster than wild-type in minimal medium

 {Delta}4 rrn and rpoB mutants share phenotypic traits

Faster growth of mutants is achieved by increased ribosome content

Fast-growing mutants display reduced hedging expression and adaptation trade-offs

Despite similar ribosomal content in rich medium the mutants present growth defects Network centrality is the most important variable determining gene sequence evolution rates in Schizosaccharomyces pombe Which variables determine the rate of gene and protein sequence evolution is a central question in evolutionary genomics. In the model organism fission yeast (Schizosaccharomyces pombe), the determinants of the rate of sequence evolution have yet to be determined. Previous studies in other organisms have typically found gene expression levels to be most significant, with numerous other variables identified as having a smaller impact. Here, partial least squares regression (PLS) and partial correlation analysis are used to model sequence evolution rates in the fission yeast genome by a range of variables. Variable importance in projection (VIP) scores as well as partial correlation coefficients are calculated for each variable, and used as estimates of the influence of each independent variable on sequence evolution rate. Unlike many previous studies in other organisms, centrality in the PPI network is shown to be the most important variable, and gene expression found to be less influential. Considerable heterogeneity is found in the influence of different gene ontology terms as well as amino acid composition. However, the majority of variance in constraint in fission yeast remains unexplained by this study, indicating that variables not yet considered and stochastics probably have considerable impact on the rate of molecular evolution. The diverse effects of phenotypic dominance on hybrid fitness When divergent populations interbreed, their alleles are brought together in hybrids. These hybrids may express novel phenotypes, not previously exposed to selection. In the initial F1 cross, most divergent alleles are present as heterozygotes. Therefore, F1 fitness can be influenced by dominance effects that first appear together in the hybrids, and so could not have been selected to function well together. We present a systematic study of these F1 dominance effects by introducing variable phenotypic dominance into Fishers geometric model. We show that dominance often reduces hybrid fitness, which can lead to patterns of optimal outbreeding and a steady decline in F1 fitness at high levels of divergence. We also show that "lucky" beneficial effects sometimes arise by chance, which might be especially important when hybrids can access novel environments. We then explore the interaction of phenotypic dominance with uniparental inheritance, showing that dominance can lead to violations of Haldanes Rule (reduced fitness of the heterogametic sex) while strengthening Darwins Corollary (fitness differences between cross directions). Taken together, our results show that dominance could play an important role in the outcomes of hybridisation after secondary contact, and thus to the maintenance or collapse of isolating barriers. Nevertheless, the telltale signs of dominance are relatively few and subtle. Results also suggest that dominance effects are smaller than the cost of segregation variance, implying that simple additive models may still give good predictions for later-generation recombinant hybrids, even when dominance qualitatively alters outcomes for the F1. Including random effects in statistical models in ecology: fewer than five levels? As generalized linear mixed-effects models (GLMMs) have become a widespread tool in ecology, the need to guide the use of such tools is increasingly important. One common guideline is that one needs at least five levels of a random effect. Having such few levels makes the estimation of the variance of random effects terms (such as ecological sites, individuals, or populations) difficult, but it need not muddy ones ability to estimate fixed effects terms - which are often of primary interest in ecology. Here, I simulate ecological datasets and fit simple models and show that having too few random effects terms does not influence the parameter estimates or uncertainty around those estimates for fixed effects terms. Thus, it should be acceptable to use fewer levels of random effects if one is not interested in making inference about the random effects terms (i.e. they are  nuisance parameters used to group non-independent data). I also use simulations to assess the potential for pseudoreplication in (generalized) linear models (LMs), when random effects are explicitly ignored and find that LMs do not show increased type-I errors compared to their mixed-effects model counterparts. Instead, LM uncertainty (and p values) appears to be more conservative in an analysis with a real ecological dataset presented here. These results challenge the view that it is never appropriate to model random effects terms with fewer than five levels - specifically when inference is not being made for the random effects, but suggest that in simple cases LMs might be robust to ignored random effects terms. Given the widespread accessibility of GLMMs in ecology and evolution, future simulation studies and further assessments of these statistical methods are necessary to understand the consequences of both violating and blindly following simple guidelines. Unravelling Rubber Tree Growth by Integrating GWAS and Biological Network-Based Approaches Hevea brasiliensis (rubber tree) is a large tree species of the Euphorbiaceae family with inestimable economic importance. Rubber tree breeding programs currently aim to improve growth and production, and the use of early genotype selection technologies can accelerate such processes, mainly with the incorporation of genomic tools, such as marker-assisted selection (MAS). However, few quantitative trait loci (QTLs) have been used successfully in MAS for complex characteristics. Recent research shows the efficiency of genome-wide association studies (GWAS) for locating QTL regions in different populations. In this way, the integration of GWAS, RNA-sequencing (RNA-Seq) methodologies, coexpression networks and enzyme networks can provide a better understanding of the molecular relationships involved in the definition of the phenotypes of interest, supplying research support for the development of appropriate genomic based strategies for breeding. In this context, this work presents the potential of using combined multiomics to decipher the mechanisms of genotype and phenotype associations involved in the growth of rubber trees. Using GWAS from a genotyping-by-sequencing (GBS) Hevea population, we were able to identify molecular markers in QTL regions with a main effect on rubber tree plant growth under constant water stress. The underlying genes were evaluated and incorporated into a gene coexpression network modelled with an assembled RNA-Seq-based transcriptome of the species, where novel gene relationships were estimated and evaluated through in silico methodologies, including an estimated enzymatic network. From all these analyses, we were able to estimate not only the main genes involved in defining the phenotype but also the interactions between a core of genes related to rubber tree growth at the transcriptional and translational levels. This work was the first to integrate multiomics analysis into the in-depth investigation of rubber tree plant growth, producing useful data for future genetic studies in the species and enhancing the efficiency of the species improvement programs. 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. Distinct repertoire development history during affinity maturation of antigen-responding VHH antibody Time-resolved antibody repertoire analysis is essential to elucidate the antibody maturation process during immunization. However, most analyses are merely snapshot comparisons of unimmunized and fully immunized animals. In this study, an alpaca was immunized with multiple antigen injections and its repertoire development was traced by high-throughput next-generation sequencing every week for 3 months. The use of alpaca-derived, single-domain antibody circumvents the problems associated with conventional antibodies including heavy-chain and light-chain matching and labor-intensive production. The sequences were processed into clusters and the antibodies were generated to determine whether the clusters included antigen-binding antibodies. For most of the antigen-responding clusters, repertoire development exhibited distinct time-dependent patterns compared with clusters unresponsive to antigens. Moreover, the sequences at the tips and root of the phylogenetic tree had strong and weak antigen-binding affinity, respectively. The foregoing characteristics could be exploited in the prediction of antigen-responding clusters and the physical properties of antibodies.

TeaserClusters of antigen-responding, single-domain antibodies exhibited distinct time-dependent repertoire development history A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. 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. From GWAS Variant to Function: a Study of ~148,000 Variants for Blood Cell Traits. Genome-wide association studies (GWAS) have identified hundreds of thousands of genetic variants associated with complex diseases and traits. However, most variants are noncoding and not clearly linked to genes, making it challenging to interpret these GWAS signals. We present a systematic variant-to-function study, prioritizing the most likely functional elements of the genome for experimental follow-up, for >148,000 variants identified for hematological traits. Specifically, we developed VAMPIRE: Variant Annotation Method Pointing to Interesting Regulatory Effects, an interactive web application implemented in R Shiny (http://shiny.bios.unc.edu/vampire/). This tool efficiently integrates and displays information from multiple complementary sources, including epigenomic signatures from blood cell relevant tissues or cells, functional and conservation summary scores, variant impact on protein and gene expression, chromatin conformation information, as well as publicly available GWAS and phenome-wide association study (PheWAS) results. Leveraging data generated from independently performed functional validation experiments, we demonstrate that our prioritized variants, genes, or variant-gene links are significantly more likely to be experimentally validated. This study not only has important implications for systematic and efficient revelation of functional mechanisms underlying GWAS variants for hematological traits, but also provides a prototype that can be adapted to many other complex traits, paving the path for efficient variant to function (V2F) analyses. Rescue of stalled clathrin-mediated endocytosis by asymmetric Arp2/3-mediated actin assembly Actin assembly facilitates vesicle formation in several trafficking pathways. Clathrin-mediated endocytosis (CME) shows elevated actin assembly dependence under high membrane tension. Why actin assembly at CME sites occurs heterogeneously even within the same cell, and how assembly forces are harnessed, are not fully understood. Here, endocytic dynamics, actin presence, and geometry of CME proteins from three different functional modules, were analyzed using three-dimensional (3D) super-resolution microscopy, live-cell imaging, and machine-learning-based computation. When hundreds of CME events were compared, sites with actin assembly showed a distinct signature, a delay between completion of coat expansion and vesicle scission, indicating that actin assembly occurs preferentially at stalled CME sites. N-WASP is recruited to one side of CME sites where it is positioned to stimulate asymmetric actin assembly. We propose that asymmetric actin assembly rescues stalled CME sites by pulling vesicles into the cell much like a bottle opener pulls off a bottle cap.