Genomic databanks and targeted assays help characterise domestic mosquito incursions Biosecurity strategies aiming to restrict the spread of invasive pests can benefit from knowing where new incursions have come from. This knowledge can be acquired using genomic databanks, where genetic variation in incursion samples is compared against variation in reference samples. Here, we use genomic databanks to investigate domestic incursions of two mosquito species in Australia. We used a deep learning method to trace a 2021 invasion of Aedes aegypti in Tennant Creek, Northern Territory, to Townsville, Queensland, and to trace Ae. albopictus incursions between several invaded and uninvaded islands in the Torres Strait. We observed high precision of tracing despite 30-70 generations separating incursions and source samples. Targeted assays also improved precision, in this case by comparing Wolbachia infection data to mitochondrial DNA variation in Ae. aegypti. For Ae. albopictus, incursions detected in 2021 were traced to different islands than incursions detected in 2019. Patterns of relatedness and inbreeding indicated that Tennant Creek was likely invaded by Ae. aegypti from one family, whereas Torres Strait incursions involved distinct kinship groups. Our results highlight the value of repeated tracing with genomic databanks that remain informative over years, and demonstrate how additional targeted assays (e.g. Wolbachia) can improve inferences.
Host Variation in Type I Interferon Signaling Genes (MX1), CCR5{triangleup}32, and MHC Class I Alleles in Treated HIV+ Non-Controllers Predict Viral Reservoir Size Prior genomewide association studies have identified variation in MHC Class I alleles and CCR5{Delta}32 as genetic predictors of viral control, especially in "elite" controllers, individuals who remain virally suppressed in the absence of therapy. We analyzed custom host whole exome sequencing and direct HLA typing from 202 ART-suppressed HIV+ non-controllers in relation to four measures of the peripheral CD4+ T cell reservoir: HIV intact DNA, total (t)DNA, unspliced (us)RNA, and RNA/DNA. Linear mixed models were adjusted for potential covariates including age, sex, nadir CD4+ T cell count, pre-ART HIV RNA, timing of ART initiation from date of detected HIV infection, and duration of ART suppression. Previously reported HLA "protective" B*57:01 (in untreated elite controllers) was significantly associated with lower usRNA (p=4.2x10-4, Benjamini-Hochberg false discovery rate [FDR] q=3.3x10-3). The favorable CCR5{Delta}32 mutation was also associated with lower tDNA and usRNA, primarily in Europeans (p=4.3x10-3 and 8.7x10-3). Several SNPs in interferon stimulated gene, MX1, were associated with HIV total DNA (p=1.3x10-7; q=0.02), and the direction of these associations paralleled predicted eQTL (expression quantitative trait) MX1 gene expression. Among Europeans, gene set enrichment analysis identified genes involved in "negative regulation of host cell viral entry" to be significantly associated with usRNA (q=0.03). These findings suggest a potential role for interferon signaling, MHC class I alleles, and CCR5{Delta}32 in maintaining the reservoir size in HIV non-controllers. These data need to be further validated in functional studies and replication cohorts.
Sexual signals persist over deep time: ancient co-option of bioluminescence for courtship displays in cypridinid ostracods Although the the diversity, beauty, and intricacy of sexually selected courtship displays command the attention of evolutionists, the longevity of these traits in deep time is poorly understood. Population-based theory suggests sexual selection could either lower or raise extinction risk, resulting in high or low persistence of lineages with sexually selected traits. Furthermore, empirical studies that directly estimate longevity of sexually selected traits are uncommon. Sexually selected signals - including bioluminescent courtship - originated multiple times during evolution, allowing empirical study of their longevity after careful phylogenetic and divergence time analyses. Here, we estimate the first transcriptome-based molecular phylogeny and divergence times of Cypridinidae. We report extreme longevity of bioluminescent courtship, a trait important in mate choice and probably under sexual selection. Our relaxed-clock estimates of divergence times coupled with stochastic character mapping show luminous courtship evolved only once in Cypridinidae in a Sub-Tribe we name Luxorina at least 151 Million Years Ago (Ma) from cypridinid ancestors that used bioluminescence only in anti-predator displays, defining a Tribe we name Luminini. This time-calibrated molecular phylogeny of cypridinids will serve as a foundation for integrative and comparative studies on the biochemistry, molecular evolution, courtship, diversification, and ecology of cypridinid bioluminescence. The persistence of luminous courtship for hundreds of millions of years indicates that rates of speciation within the group exceeded extinction risk, allowing for the persistence of a diverse clade of signalling species and that sexual selection did not cause rapid loss of associated traits. Male age and Wolbachia dynamics: Investigating how fast and why bacterial densities and cytoplasmic incompatibility strengths vary Endosymbionts can influence host reproduction and fitness to favor their maternal transmission. For example, endosymbiotic Wolbachia bacteria often cause cytoplasmic incompatibility (CI) that kills uninfected embryos fertilized by Wolbachia-modified sperm. Infected females can rescue CI, providing them a relative fitness advantage. Wolbachia-induced CI strength varies widely and tends to decrease as host males age. Since strong CI drives Wolbachia to high equilibrium frequencies, understanding how fast and why CI strength declines with male age is crucial to explaining age-dependent CIs influence on Wolbachia prevalence. Here, we investigate if Wolbachia densities and/or CI gene (cif) expression covary with CI-strength variation and explore covariates of age-dependent Wolbachia-density variation in two classic CI systems. wRi CI strength decreases slowly with Drosophila simulans male age (6%/ day), but wMel CI strength decreases very rapidly (19%/ day), yielding statistically insignificant CI after only three days of D. melanogaster emergence. Wolbachia densities and cif expression in testes decrease as wRi-infected males age, but both surprisingly increase as wMel-infected males age, and CI strength declines. We then tested if phage lysis, Octomom copy number (which impacts wMel density), or host immune expression covary with age-dependent wMel densities--only host immune expression correlated with density. Together, our results identify how fast CI strength declines with male age in two model systems and reveal unique relationships between male age, Wolbachia densities, cif expression, and host immunity. We discuss new hypotheses about the basis of age-dependent CI strength and its contributions to Wolbachia prevalence.

ImportanceWolbachia are the most common animal-associated endosymbionts due in large part to their manipulation of host reproduction. Many Wolbachia cause cytoplasmic incompatibility (CI) that kills uninfected host eggs. Infected eggs are protected from CI, favoring Wolbachia spread in natural systems and in transinfected mosquito populations where vector-control groups use strong CI to maintain pathogen-blocking Wolbachia at high frequencies for biocontrol of arboviruses. CI strength varies considerably in nature and declines as males age for unknown reasons. Here, we determine that CI strength weakens at different rates with age in two model symbioses. Wolbachia density and CI gene expression covary with wRi-induced CI strength in Drosophila simulans, but neither explain rapidly declining wMel-induced CI in aging D. melanogaster males. Patterns of host immune gene expression suggest a candidate mechanism behind age-dependent wMel densities. These findings inform how age-dependent CI may contribute to Wolbachia prevalence in natural systems and potentially in transinfected systems. Linking null models to hypotheses to improve nestedness analysis O_LINestedness is a common pattern in metacommunities and interaction networks, whose causes are still subject to investigation. Nestedness inference is challenging because, beyond calculating an index, we need to compare observed and expected values using null models. There are many different null models and the choice between them strongly affects test outcomes. Furthermore, there is no theoretical basis to guide this choice.
C_LIO_LIHere, we propose a different look at the meaning of nestedness that improves our understanding of its causes and directly links null models to hypotheses. The nestedness of a matrix is the combination of marginal sum inequality and high overlap. The higher the overlap, the more predictable the cell values by marginal sums. Therefore, the only way a process can increase nestedness is by promoting marginal sum inequalities, without concomitantly introducing preferences. We propose a combination of null models to test for different topological hypotheses. In these null models, interactions are randomized following a probability matrix. The probability matrix of the equiprobable model excludes marginal sum inequalities and consequently all nestedness-generating mechanisms, providing a distribution of expected values for nestedness significance tests. The probability matrix of the proportional model, by conserving marginal sums and maximizing overlap, conserves nestedness-generating mechanisms and excludes nestedness-disrupting mechanisms, thereby yielding the expected nestedness values for fully nested matrices.
C_LIO_LIWe advance a new protocol to test for nestedness in binary and weighted matrices, in which those two null models represent different topological hypotheses. Our protocol helps to tell apart cases in which nestedness is different from expected by chance and cases in which nestedness represents the best topological archetype. We also evaluate the efficiency of applying this protocol using some common nestedness indices. Analyses based on weighted data proved to be highly efficient and informative, while binary analyses induced severe distortions. Finally, we empirically illustrate our approach using a plant-pollinator network.
C_LIO_LIThrough a shift of perspective, our approach reconciliates contradictions in null model analysis and provides a guideline for nestedness interpretation.
C_LI Towards a unified model of naive T cell dynamics across the lifespan Naive CD4 and CD8 T cells are cornerstones of adaptive immunity, but the dynamics of their establishment early in life and how their kinetics change as they mature following release from the thymus are poorly understood. Further, due to the diverse signals implicated in naive T cell survival, it has been a long-held and conceptually attractive view that they are sustained by active homeostatic control as thymic activity wanes. Here we employ multiple experimental systems to identify a unified model of naive CD4 and CD8 T cell population dynamics across the mouse lifespan. We infer that both subsets divide rarely and progressively increase their survival capacity with cell age. Strikingly, this simple model captures naive CD4 T cell dynamics throughout life. In contrast, we find that newly generated naive CD8 T cells are lost more rapidly during the first 3-4 weeks of life, likely due to increased recruitment into memory. We find no evidence for elevated division rates in neonates, or for feedback regulation of naive T cell numbers at any age. We show how confronting mathematical models with diverse datasets can reveal a quantitative and remarkably simple picture of naive T cell dynamics from birth into old age. Brain expressed FKBP5 delineates a therapeutic subtype of severe mental illness Deducing genes capable of classifying biologically-distinct psychiatric subtypes, and their targets for treatment, is a priority approach in psychiatry. FKBP5 is one such gene with strong evidence of utility to delineate a trans-diagnostic psychiatric subtype. Yet how brain-expressed FKBP5 is affected in psychiatric disorders in humans is not fully understood and critical for propelling FKBP5-targeting treatment development. We performed a large-scale postmortem study (n=895) of FKBP5 using dorsolateral prefrontal cortex samples derived from individuals with severe psychiatric disorders with a comprehensive battery of bulk/single-cell omics and histological analyses. We observed consistently heightened FKBP5 mRNA and protein in psychopathology, moderated by genotype and age, and accompanied by DNA methylation changes in key enhancers. These effects were most prominent in superficial-layer pyramidal cells. Heightened FKBP5 was also differentially associated with downstream pathways according to age, being specifically associated with synaptic transmission in early adulthood, and neuroinflammation and neurodegeneration in later life. One-carbon pathway metabolites are altered in the plasma of subjects with Down syndrome: relation to chromosomal dosage Down syndrome (DS) is the most common chromosomal disorder and it is caused by trisomy of chromosome 21 (Hsa21). Subjects with DS show a large heterogeneity of phenotypes and the most constant clinical features present are typical facies and intellectual disability (ID). Several studies demonstrated that trisomy 21 causes an alteration in the metabolic profile, involving among all one-carbon cycle. We performed enzyme-linked immunosorbent assays (ELISAs) to identify the concentration of 5 different intermediates of the one-carbon cycle, namely tetrahydrofolate (THF), 5-methyltetrahydrofolate (5-methyl-THF), 5-formyltetrahydrofolate (5-formyl-THF), S-adenosyl-homocysteine (SAH) and S-adenosyl-methionine (SAM) in plasma samples obtained from a total of 164 subjects with DS compared to 54 euploid subjects. Results highlight specific alterations of THF with a median concentration ratio DS/control of 2:3, a decrease of a necessary molecule perfectly consistent with a chromosomal dosage effect. Moreover, SAM and SAH show a ratio of 2:1 and 3.6:1, respectively. The relevance of these results for the biology of intelligence and its impairment in trisomy 21 is discussed, leading to the final proposal of 5-methyl-THF as the best candidate for a clinical trial aimed at restoring the dysregulation of one-carbon cycle in trisomy 21, possibly improving cognitive skills of subjects with DS. Royal knifefish generate powerful suction feeding through large neurocranial elevation and high epaxial muscle power Suction feeding in ray-finned fishes involves powerful buccal cavity expansion to accelerate water and food into the mouth. Previous XROMM studies in largemouth bass (Micropterus salmoides), bluegill sunfish (Lepomis macrochirus), and channel catfish (Ictalurus punctatus) have shown that more than 90% of suction power in high performance strikes comes from the axial musculature. Thus, the shape of the axial muscles and skeleton may impact suction feeding mechanics. Royal knifefish (Chitala blanci) have an unusual postcranial morphology, with a ventrally flexed vertebral column and relatively large mass of epaxial muscle. Based on their body shape, we hypothesized that royal knifefish would generate high power strikes by utilizing large neurocranial elevation, vertebral column extension, and epaxial shortening. As predicted, C. blanci generated high suction expansion power compared to the other three species studied to date (up to 160 W), which was achieved by increasing both the rate of volume change and the intraoral subambient pressure. The large epaxial muscle (25% of body mass) shortened at high velocities to produce large neurocranial elevation and vertebral extension (up to 41 deg, combined), as well as high muscle mass-specific power (up to 800 W kg-1). For the highest power strikes, axial muscles generated 95% of the power, and 64% of the axial muscle mass consisted of the epaxial muscles. The epaxial-dominated suction expansion of royal knifefish supports our hypothesis that postcranial morphology may be a strong predictor of suction feeding biomechanics.

SUMMARY STATEMENTRoyal knifefish rely on their distinct postcranial morphology--with a curved vertebral column and large dorsal body muscles--to produce large neurocranial elevation and powerful suction feeding. Different contribution of the monkey prefrontal and premotor dorsal cortex in decision making during a transitive inference task Several studies have reported similar neural modulations between brain areas of the frontal cortex, such as the dorsolateral prefrontal (DLPFC) and the premotor dorsal (PMd) cortex, in tasks requiring encoding of the abstract rules for selecting the proper action. Here we compared the neuronal modulation of the DLPFC and PMd of monkeys trained to choose the higher rank from a pair of abstract images (target item), selected from an arbitrarily rank-ordered set (A>B>C>D>E>F) in the context of a transitive inference task. Once acquired by trial-and-error, the ordinal relationship between pairs of adjacent images (i.e., A>B; B>C; C>D; D>E; E>F), monkeys were tested in indicating the ordinal relation between items of the list not paired during learning. During these decisions, we observed that the choice accuracy increased and the reaction time decreased as the rank difference between the compared items enhanced. This result is in line with the hypothesis that after learning, the monkeys built an abstract mental representation of the ranked items, where rank comparisons correspond to the items position comparison on this representation. In both brain areas, we observed higher neuronal activity when the target item appeared in a specific location on the screen with respect to the opposite position and that this difference was particularly enhanced at lower degrees of difficulty. By comparing the time evolution of the activity of the two areas, we observed that the neural encoding of target item spatial position occurred earlier in the DLPFC than in the PMd. Load adaptation of endocytic actin networks Clathrin-mediated endocytosis (CME) robustness under elevated membrane tension is maintained by actin assembly-mediated force generation. However, whether more actin assembles at endocytic sites in response to increased load, as has been observed in lamellipodia, has not previously been investigated. Here actin network ultrastructure at CME sites was examined under low and high membrane tension. Actin and N-WASP spatial organization indicate that actin polymerization initiates at the base of clathrin-coated pits and that the network then grows away from the plasma membrane. Actin network height at individual CME sites was not coupled to coat shape, raising the possibility that local differences in mechanical load feedback on assembly. By manipulating membrane tension and Arp2/3 complex activity we tested the hypothesis that actin assembly at CME sites increases in response to elevated load. Indeed, in response to elevated membrane tension, actin grew higher, resulting in greater coverage of the clathrin coat, and CME slowed. When membrane tension was elevated and the Arp2/3 complex was inhibited, shallow clathrin-coated pits accumulated, indicating that this adaptive mechanism is especially crucial for coat curvature generation. We propose that actin assembly increases in response to increased load to ensure CME robustness over a range of plasma membrane tensions. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. Kidney failure alters parathyroid Pin1 phosphorylation and parathyroid hormone mRNA binding proteins leading to secondary hyperparathyroidism Parathyroid hormone (PTH) regulates calcium metabolism and bone strength. Chronic kidney disease leads to secondary hyperparathyroidism (SHP) which increases morbidity and mortality. In experimental SHP, the increased PTH gene expression is due to enhanced PTH mRNA stability, mediated by changes in its interaction with stabilizing AUF1 and destabilizing KSRP. Pin1 isomerizes target proteins and leads to KSRP dephosphorylation. In SHP, Pin1 isomerase activity is decreased and phosphorylated KSRP fails to bind PTH mRNA, resulting in high PTH mRNA stability and levels. However, the up- and down-stream mechanisms by which kidney failure stimulates the parathyroid remain elusive. We now reveal a pathway where kidney failure induces parathyroid Pin1 phosphorylation, linking Pin1, KSRP and PTH mRNA stability as key players for the high PTH expression in SHP. We identified by mass-spectrometry, changes in rat parathyroid proteome and phosphoproteome profiles induced by impaired renal function, including KSRP phosphorylation at Pin1 target sites. Furthermore, both acute and chronic kidney failure led to parathyroid-specific Pin1 Ser16 and Ser71 phosphorylation, which disrupts Pin1 activity. Accordingly, pharmacologic Pin1 inhibition, that mimics the decreased Pin1 activity in SHP, increased PTH expression ex-vivo in parathyroid organ cultures and in transfected cells, through the PTH mRNA protein-interacting element and KSRP phosphorylation at potential Pin1-binding motifs. Therefore, kidney failure leads to loss of parathyroid Pin1 activity by inducing Pin1 phosphorylation. This predisposes parathyroids to increase PTH production through impaired PTH mRNA decay that is dependent on KSRP phosphorylation. Pin1 and KSRP phosphorylation and the Pin1-KSRP-PTH mRNA axis thus drive SHP. Activation of the type 3 secretion system of enteropathogenic E. coli leads to remodeling of its membrane composition The cell envelope of gram-negative bacteria is a complex structure, essential for bacterial survival and for resistance to many antibiotics. Channels that cross the bacterial envelope and the host cell membrane form secretion systems that are activated upon attachment to host, enabling bacteria to inject effector molecules into the host cell, required for bacterial-host interaction. The type III secretion system (T3SS) is critical for the virulence of several pathogenic bacteria, including enteropathogenic E. coli (EPEC). The EPEC T3SS activation is associated with repression of carbon storage regulator (CsrA), resulting in gene expression remodeling, which is known to affect EPEC central carbon metabolism and contributes to the adaptation to a cell-adherent lifestyle in a poorly understood manner. We reasoned that the changes in bacterial envelope upon attachment to host and the activation of a secretion system may involve a modification of the lipid composition of bacterial envelope. Accordingly, we performed a lipidomics analysis on mutant strains that simulate T3SS activation. We saw a shift in glycerophospholipid metabolism towards the formation of lysophospholipids, attributed to corresponding upregulation of the phospholipase pldA and the acyltransferase ygiH upon T3SS activation in EPEC. We also detected a shift from menaquinones and ubiquinones to undecaprenyl lipids, concomitant to abnormal synthesis of O-antigen. The remodeling of lipid metabolism is mediated by CsrA and associated with increased bacteria cell size and Zeta potential, and a corresponding alteration in EPEC permeability to vancomycin, increasing the sensitivity of T3SS-activated strains and of adherent wild type EPEC to the antibiotic.

ImportanceThe characterization of EPEC membrane lipid metabolism upon attachment to host is an important step towards a better understanding the shift of EPEC, a notable human pathogen, from a planktonic to adherent life style. It may also apply to other pathogenic bacteria that use this secretion system. We predict that upon attachment to host cells the lipid remodeling upon T3SS activation contributes to bacterial fitness and promotes host colonization, and show that it is associated with increased cell permeability and higher sensitivity to vancomycin.

To the best of our knowledge, this is the first demonstration of a bacterial lipid remodeling due to activation of a secretion system. Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: the Multi-Targeting Drug DREAM Challenge A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ( polypharmacology). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Author SummaryMany modern drugs are developed with the goal of modulating a single cellular pathway or target. However, many drugs are, in fact,  dirty; they target multiple cellular pathways or targets. This phenomenon is known as multi-targeting or polypharmacology. While some strive to develop  cleaner therapeutics that eliminate secondary targets, recent work has shown that multi-targeting therapeutics have key advantages for a variety of diseases. However, while multi-targeting drugs that affect a precisely-defined profile of targets may be more effective, it is difficult to computationally predict which molecules have desirable target profiles. Here, we report the results of a competitive crowdsourcing project (the Multi-Targeting Drug DREAM Challenge), where we challenged participants to predict chemicals that have desired target profiles for cancer and neurodegenerative disease. Enhancers regulate 3' end processing activityto control expression of alternative 3'UTR isoforms Multi-UTR genes are widely transcribed but show cell type-specific 3'UTR isoform expression. As transcriptional enhancers regulate mRNA expression, we investigated if they also regulate mRNA isoform expression. Deletion of an endogenous enhancer of a multi-UTR gene did not impair transcript production but prevented a switch in 3'UTR isoform expression. Also, the same enhancers are able to increase transcript production in the context of single-UTR gene promoters, but they increase 3' end processing activity when paired with multi-UTR gene promoters. We show that transcription factors regulate processing activity of weak polyadenylation sites to control cell type-specific alternative 3'UTR isoform expression of widely expressed genes. Cis-Regulatory Element Hijacking by Structural Variants Overshadows Higher-Order Topological Changes in Primary Prostate Cancer Prostate cancer is a heterogeneous disease whose progression is linked to genome instability. However the impact of this instability on the three-dimensional chromatin organization and how this drives progression is unclear. Using primary benign and tumour tissue, we find a high concordance in the higher-order three-dimensional genome organization across normal and prostate cancer cells. This concordance argues for constraints to the topology of prostate tumour genomes. Nonetheless, we identify changes to focal chromatin interactions and show how structural variants can induce these changes to guide cis-regulatory element hijacking. Such events result in opposing differential expression on genes found at antipodes of rearrangements. Collectively, our results argue that cis-regulatory element hijacking from structural variant-induced altered focal chromatin interactions overshadows higher-order topological changes in the development of primary prostate cancer. MetaTrass: High-quality metagenome assembly on the human gut microbiome by co-barcoding sequencing reads With the development of sequencing technologies and computational analysis in metagenomics, the genetic diversity of non-conserved regions has been receiving intensive attention to unravel the human gut microbial community. However, it remains a challenge to obtain enough microbial draft genomes at a high resolution from a single sample. In this work, we presented MetaTrass with a strategy of binning first and assembling later to assemble high-quality draft genomes based on metagenomics co-barcoding reads and the public reference genomes. We applied the tool to the single tube long fragment reads datasets for four human faecal samples, and generated more high-quality draft genomes with longer contiguity and higher resolution than the common combination strategies of genome assembling and binning. A total of 178 high-quality genomes was successfully assembled by MetaTrass, but the maximum of 58 was generated by the optimal common combination strategy in our tests. These high-quality genomes paved the way for genetic diversity and lineage analysis among different samples. With the high capability of assembling high-quality genomes of metagenomics datasets, MetaTrass will facilitate the study of spatial characters and dynamics of complex microbial communities at high resolution. The open-source code of MetaTrass is available at https://github.com/BGI-Qingdao/MetaTrass. Multiplex knockout of trichome-regulating MYB duplicates in hybrid poplar using a single gRNA As the focus for CRISPR edited plants moves from proof-of-concept to real world applications, precise gene manipulation will increasingly require concurrent multiplex editing for polygenic traits. A common approach for editing across multiple sites is to design one gRNA per target; however, this complicates construct assembly and increases the possibility of off-target mutations. In this study, we utilized one gRNA to target MYB186, a known positive trichome regulator, as well as its paralogs MYB138 and MYB38 at a consensus site for mutagenesis in Populus tremula x P. alba INRA 717-1B4. Unexpected duplications of MYB186 and MYB138 resulted in a total of eight alleles for the three targeted genes in the hybrid poplar. Deep sequencing and PCR analyses confirmed editing across all eight targets in nearly all of the resultant glabrous mutants, ranging from small indels to large genomic dropouts, with no off-target activity detected at four potential sites. This highlights the effectiveness of a single gRNA targeting conserved exonic regions for multiplex editing. Additionally, cuticular wax and whole leaf analyses showed a complete absence of triterpenes in the trichomeless mutants, hinting at a previously undescribed role for the non-glandular trichomes of poplar.

ONE SENTENCE SUMMARYTargeting conserved sequences with a single gRNA allowed efficient mutagenesis of a multigene family and the recovery of trichomeless and triterpene-free poplar mutants. Gradients in cellular metabolism regulate development of tonotopy along the cochlea. In vertebrates with elongated auditory organs, mechanosensory hair cells (HCs) are organised such that complex sounds are broken down into their component frequencies along the basal-to-apical long (tonotopic) axis. To generate the frequency-specific characteristics required at the appropriate positions, nascent HCs must determine their tonotopic properties during development. This relies on complex signalling between the developing HC and its local niche along the axis within the auditory organ. Here we apply NAD(P)H fluorescence lifetime imaging (FLIM) and live imaging of mitochondria to reveal metabolic gradients along the tonotopic axis of the developing cochlea. We further show that re-shaping these gradients during development, by inhibiting cytosolic glucose metabolism, alters normal Bmp7 and Chordin like-1 signalling leading to flattening of tonotopic patterning along the axis. Our work supports a causal link between morphogen signalling and metabolic reprogramming in specifying tonotopic morphologies in auditory HCs. Effect of sublethal prenatal endotoxaemia on murine placental transport systems and lipid homeostasis Infection alters the expression of transporters that mediate the placental exchange of xenobiotics, lipids and cytokines. We hypothesized that lipopolysaccharide (LPS) modifies the expression of placental transport systems and lipid homeostasis. LPS (150 g/kg; i.p.) treatments were administered for 4 h or 24 h, animals were euthanized at gestational days (GD) 15.5 or 18.5, and maternal blood, foetuses and placentae were collected. Increased rates of foetal demise were observed at GD15.5 following LPS treatment, whereas at GD18.5, high rates of early labour occurred and were associated with distinct proinflammatory responses. LPS did not alter ABC transporter mRNA expression but decreased Fabppm at GD15.5 (LPS-4 h) and increased Fat/Cd36 lipid transporter mRNA at GD18.5 (LPS-4 h). At the protein level, breast cancer-related protein (BCRP) and Abcg1 levels were decreased in the placental labyrinth zone (Lz) at GD15.5, whereas P-glycoprotein (P-gp) and Bcrp Lz-immunostaining was decreased at GD18.5. In the placental junctional zone (Jz), P-gp, Bcrp and Abcg1 levels were higher at GD18.5. Specific maternal plasma and placental changes in triacylglycerol, free fatty acid, cholesterol, cholesterol ester and monoacylglycerol levels were detected in a gestational age-dependent manner. In conclusion, LPS-induced foetal death and early labour were associated with altered placental ABC and lipid transporter expression and deranged maternal plasma and placental lipid homeostasis. These changes likely modify foetal xenobiotic exposure and placental lipid exchange in cases of bacterial infection. Production and characterization of virus-free, CRISPR-CAR T cells capable of inducing solid tumor regression BackgroundChimeric antigen receptor (CAR) T cells traditionally harbor viral vector-based sequences that encode the CAR transgene in the genome. These T cell products have yet to show consistent anti-tumor activity in patients with solid tumors. Further, viral vector manufacturing is resource intensive, suffers from batch-to-batch variability, and includes several animal components, adding regulatory and supply chain pressures.

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

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

ConclusionsThis proof-of-principle study leveraging virus-free genome editing technology could enable flexible manufacturing of clinically relevant, high-quality CAR T cells to treat cancers, including solid tumors. 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. 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. 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. G Protein-Coupled Receptor-Ligand Dissociation Rates and Mechanisms from {tau}RAMD Simulations There is a growing appreciation of the importance of drug-target binding kinetics for lead optimization. For G protein-coupled receptors (GPCRs), which mediate signaling over a wide range of timescales, the drug dissociation rate is often a better predictor of in vivo efficacy than binding affinity, although it is more challenging to compute. Here, we assess the ability of the {tau}-Random Acceleration Molecular Dynamics ({tau}RAMD) approach to reproduce relative residence times and reveal dissociation mechanisms and the effects of allosteric modulation for two important membrane-embedded drug targets: the {beta}2-adrenergic receptor and the muscarinic acetylcholine receptor M2. The dissociation mechanisms observed in the relatively short RAMD simulations (in which molecular dynamics (MD) simulations are performed using an additional force with an adaptively assigned random orientation applied to the ligand) are in general agreement with much more computationally intensive conventional MD and metadynamics simulations. Remarkably, although decreasing the magnitude of the random force generally reduces the number of egress routes observed, the ranking of ligands by dissociation rate is hardly affected and agrees well with experiment. The simulations also reproduce changes in residence time due to allosteric modulation and reveal associated changes in ligand dissociation pathways.

TABLE OF CONTENTS GRAPHIC

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=100 SRC="FIGDIR/small/449151v2_ufig1.gif" ALT="Figure 1">
View larger version (43K):
org.highwire.dtl.DTLVardef@743dcdorg.highwire.dtl.DTLVardef@1857381org.highwire.dtl.DTLVardef@62359aorg.highwire.dtl.DTLVardef@1ad38a3_HPS_FORMAT_FIGEXP  M_FIG C_FIG