A common SEC23B missense mutation in congenital dyserythropoietic anemia type II leads to growth restriction and chronic pancreatitis in mice Human loss-of-function mutations in SEC23B result in congenital dyserythropoietic anemia type II (CDAII). Complete deficiency of SEC23B in mice leads to perinatal death caused by massive degeneration of professional secretory tissues. Functions of SEC23B in postnatal mice are unclear. In this study, we generated knockout (KO) mice with deletion of exons 5 and 6 of Sec23b and knockin (KI) mice with the E109K mutation, the most common CDAII missense mutation. Despite decreased SEC23B protein level, Sec23bki/ki mice showed no obvious abnormalities. Hemizygous (Sec23bki/ko) mice exhibit a partial lethal phenotype, with only half of expected mice survive past weaning. Sec23bki/ko pancreas had exocrine insufficiency and chronic pancreatitis histology, as well as increased ER stress and apoptosis. Moreover, Sec23bki/ko mice exhibited severe growth restriction accompanied by growth hormone (GH) insensitivity, reminiscent of the Laron syndrome. Growth restriction is not associated with hepatocyte-specific Sec23b deletion, suggesting non-liver origin of the phenotype. Inflammation associated with chronic pancreatic deficiency may explain GH insensitivity in Sec23bki/ko mice. Our results indicate that phenotype severities are linked to the residual functions of SEC23B, demonstrating the importance of SEC23B in pancreatic acinar function in adult mice. The Sec23bki/ko mice provide a novel model of chronic pancreatitis and growth retardation.
Endometriotic Organoids: A Novel In Vitro Model of Endometriotic Lesion Development. The development and progression of endometriotic lesions are poorly understood, but immune cell dysfunction and inflammation are closely associated with the pathophysiology of endometriosis. A lack of suitable 3D in vitro models permitting the study of interactions between cell types and the microenvironment is a contributing factor. To address this limitation, we developed endometriotic organoids (EO) to explore the role of epithelial-stromal interactions and model peritoneal cell invasion associated with lesion development. Using a non-adherent microwell culture system, spherical organoids were generated with endometriotic epithelial cells (12Z) combined with immortalized endometriotic stromal cells (iEc-ESC) or immortalized uterine stromal cells (iHUF). Organoids self-organized with stromal cells occupying the center and epithelial cells on the periphery of the organoid. Endometriotic organoids (EO), containing iEc-ESC, resulted in the development of stratified 12Z epithelial cells compared to those with iHUF where the 12Z cells developed as a single layered epithelium. Transcriptomic analysis found 4,522 differentially expressed genes (DEG) between EO and 12Z/iHUF organoids, and the top DEG included increased expression of interleukins and prostaglandin synthase enzymes. An overlap of the EO DEG with baboon endometriotic lesions was highly significant. Finally, to mimic invasion of endometrial tissue into the peritoneum, a model was developed using EO and extracellular matrix containing human peritoneal mesothelial cells (LP9). Invasion of EO into the extracellular matrix-LP9 layer was increased in presence of estrogen or THP1-derived proinflammatory macrophages. Taken together, our results strongly support the concept that EO are an appropriate model for dissecting mechanisms that contribute to endometriotic lesion development.

One Sentence SummaryEndometriotic organoids are an appropriate model to study epithelial-stromal interactions and model cell invasion associated with lesion development.
Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. The mutation effect reaction norm (Mu-RN) highlights environmentally dependent mutation effects and genetic interactions Since the modern synthesis, the fitness effects of mutations and epistasis have been central yet provocative concepts in evolutionary and population genetics. Studies of how the interactions between parcels of genetic information can change as a function of environmental context have added a layer of complexity to these discussions. Here I introduce the "mutation effect reaction norm" (Mu-RN), a new instrument through which one can analyze the phenotypic consequences of mutations and interactions across environmental contexts. It embodies the fusion of measurements of genetic interactions with the reaction norm, a classic depiction of the performance of genotypes across environments. I demonstrate the utility of the Mu-RN through the signature of a "compensatory ratchet" mutation that undermines reverse evolution of antimicrobial resistance. More broadly, I argue that the mutation effect reaction norm may help us resolve the dynamism and unpredictability of evolution, with implications for theoretical biology, genetic modification technology, and public health. A bacterial extracellular vesicle-based intranasal vaccine against SARS-CoV-2 protects against disease and elicits neutralizing antibodies to wild-type virus and Delta variant Several vaccines have been introduced to combat the coronavirus infectious disease-2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Current SARS-CoV-2 vaccines include mRNA-containing lipid nanoparticles or adenoviral vectors that encode the SARS-CoV-2 Spike (S) protein of SARS-CoV-2, inactivated virus, or protein subunits. Despite growing success in worldwide vaccination efforts, additional capabilities may be needed in the future to address issues such as stability and storage requirements, need for vaccine boosters, desirability of different routes of administration, and emergence of SARS-CoV-2 variants such as the Delta variant. Here, we present a novel, well-characterized SARS-CoV-2 vaccine candidate based on extracellular vesicles (EVs) of Salmonella typhimurium that are decorated with the mammalian cell culture-derived Spike receptor-binding domain (RBD). RBD-conjugated outer membrane vesicles (RBD-OMVs) were used to immunize the golden Syrian hamster (Mesocricetus auratus) model of COVID-19. Intranasal immunization resulted in high titers of blood anti-RBD IgG as well as detectable mucosal responses. Neutralizing antibody activity against wild-type and Delta variants was evident in all vaccinated subjects. Upon challenge with live virus, hamsters immunized with RBD-OMV, but not animals immunized with unconjugated OMVs or a vehicle control, avoided body mass loss, had lower virus titers in bronchoalveolar lavage fluid, and experienced less severe lung pathology. Our results emphasize the value and versatility of OMV-based vaccine approaches. Characterization of methylation profiles in spontaneous preterm birth placental villous tissue Preterm birth is a global public health crisis which results in significant neonatal and maternal mortality. Yet little is known regarding the molecular mechanisms of idiopathic spontaneous preterm birth, and we have few diagnostic markers for adequate assessment of placental development and function. Previous studies of placental pathology and our transcriptomics studies suggest a role for placental maturity in idiopathic spontaneous preterm birth. It is known that placental DNA methylation changes over gestation. We hypothesized that if placental hypermaturity is present in our samples, we would observe a unique idiopathic spontaneous preterm birth DNA methylation profile potentially driving the gene expression differences we previously identified in our placental samples. Our results indicate the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern suggesting hypermaturity. Only seven significant differentially methylated regions fitting the idiopathic spontaneous preterm birth specific (relative to the controls) profile were identified, indicating unusually high similarity in DNA methylation between idiopathic spontaneous preterm birth and term birth samples. We identified an additional 1,718 significantly methylated regions in our gestational age matched controls were the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern, again indicating a striking level of similarity between the idiopathic spontaneous preterm birth and term birth samples. Pathway analysis of these regions revealed differences in genes within the WNT and Cadherin signaling pathways, both of which are essential in placental development and maturation. Taken together, these data demonstrate that the idiopathic spontaneous preterm birth samples are molecularly more mature than expected given their respective gestational age which likely impacts birth timing. FGF8-mediated signaling regulates tooth developmental pace during odontogenesis The developing human and mouse teeth constitute an ideal model system to study the regulatory mechanism underlying organ growth control due to the fact that their teeth share highly conserved and well-characterized developmental processes and their developmental tempo vary notably. In the current study, we manipulated heterogenous recombination between human and mouse dental tissues and demonstrate that the dental mesenchyme dominates the tooth developmental tempo and FGF8 could be a critical player during this developmental process. Forced activation of FGF8 signaling in the dental mesenchyme of mice promoted cell proliferation, prevented cell apoptosis via p38 and perhaps PI3K-Akt intracellular signaling, and impelled the transition of the cell cycle from G1- to S-phase in the tooth germ, resulting in the slowdown of the tooth developmental pace. Our results provide compelling evidence that extrinsic signals can profoundly affect tooth developmental tempo and the dental mesenchymal FGF8 could be a pivotal factor in controlling developmental pace in a non-cell-autonomous manner during mammalian odontogenesis. An image-computable model on how endogenous and exogenous attention differentially alter visual perception Attention alters perception across the visual field. Typically, endogenous (voluntary) and exogenous (involuntary) attention similarly improve performance in many visual tasks, but they have differential effects in some tasks. Extant models of visual attention assume that the effects of these two types of attention are identical and consequently do not explain differences between them. Here, we develop a model of spatial resolution and attention that distinguishes between endogenous and exogenous attention. We focus on texture-based segmentation as a model system because it has revealed a clear dissociation between both attention types. For a texture for which performance peaks at parafoveal locations, endogenous attention improves performance across eccentricity, whereas exogenous attention improves performance where the resolution is low (peripheral locations) but impairs it where the resolution is high (foveal locations) for the scale of the texture. Our model emulates sensory encoding to segment figures from their background and predict behavioral performance. To explain attentional effects, endogenous and exogenous attention require separate operating regimes across visual detail (spatial frequency). Our model reproduces behavioral performance across several experiments and simultaneously resolves three unexplained phenomena: (1) the parafoveal advantage in segmentation, (2) the uniform improvements across eccentricity by endogenous attention and (3) the peripheral improvements and foveal impairments by exogenous attention. Overall, we unveil a computational dissociation between each attention type and provide a generalizable framework for predicting their effects on perception across the visual field. Atomic Force Microscopy Reveals Distinct Protofilament-scale Structural Dynamics in Depolymerizing Microtubule Arrays The dynamic reorganization of microtubule-based cellular structures, such as the spindle and the axoneme, fundamentally depends on the dynamics of individual polymers within multimicrotubule arrays. A major class of enzymes, implicated in both the complete demolition and fine size control of microtubule-based arrays, are depolymerizing kinesins. How different depolymerases differently remodel microtubule arrays is poorly understood. A major technical challenge in addressing this question is that existing optical or electron-microscopy methods lack the spatial-temporal resolution to observe the dynamics of individual microtubules within larger arrays. Here we use Atomic Force Microscopy (AFM) to image depolymerizing arrays at single microtubule and protofilament resolution. We discover previously unseen modes of microtubule array destabilization by conserved depolymerases. We find that the kinesin-13 MCAK mediates asynchronous protofilament depolymerization and lattice-defect propagation, whereas the kinesin-8 Kip3p promotes synchronous protofilament depolymerization. Unexpectedly, MCAK can depolymerize the highly stable axonemal doublets, but Kip3p cannot. We propose that distinct protofilament-level activities underlie the functional dichotomy of depolymerases, resulting in either large-scale destabilization or length regulation of microtubule arrays. Our work establishes AFM as a powerful strategy to visualize microtubule dynamics within arrays and reveals how nanometer-scale substrate specificity leads to differential remodeling of micron-sized cytoskeletal structures.

SIGNIFICANCE STATEMENTOne cannot help but marvel at the precise organization of microtubule polymers in cellular structures such as the axoneme and the spindle. However, our understanding of the biochemical mechanisms that sculpt these arrays comes largely from in vitro experiments with a small number (one or two) of microtubules. This is somewhat akin to studying the architecture of multilane highways by studying one-lane streets. Here we directly visualize depolymerizing microtubule arrays at individual microtubule and protofilament resolution using Atomic Force Microscopy. Our results reveal differences in microtubule depolymerase activity and provide insights into how these differences in enzymatic activity on the nanometer scale can result in the differential remodeling of multi-microtubule arrays on the micron-length scale. Genetic architecture of natural variation of cardiac performance in flies Deciphering the genetic architecture of human cardiac disorders is of fundamental importance but their underlying complexity is a major hurdle. We investigated the natural variation of cardiac performance in the sequenced inbred lines of the Drosophila Genetic Reference Panel (DGRP)1. Genome Wide Associations Studies (GWAS) identified genetic networks associated with natural variation of cardiac traits which were used to gain insights as to the molecular and cellular processes affected. Non-coding variants that we identified were used to map potential regulatory non-coding regions, which in turn were employed to predict Transcription Factors (TFs) binding sites. Cognate TFs, many of which themselves bear polymorphisms associated with variations of cardiac performance, were also validated by heart specific knockdown. Additionally, we showed that the natural variations associated with variability in cardiac performance affect a set of genes overlapping those associated with average traits but through different variants in the same genes. Furthermore, we showed that phenotypic variability was also associated with natural variation of gene regulatory networks. More importantly, we documented correlations between genes associated with cardiac phenotypes in both flies and humans, which supports a conserved genetic architecture regulating adult cardiac function from arthropods to mammals. Specifically, roles for PAX9 and EGR2 in the regulation of the cardiac rhythm were established in both models, illustrating that the characteristics of natural variations in cardiac function identified in Drosophila can accelerate discovery in humans. Co-infection of Candidatus Piscichlamydia trichopodus (order Chlamydiales) and Henneguya sp (Myxosporea, Myxobolidae) in snakeskin gourami Trichopodus pectoralis (Regan 1910) The present study describes a simultaneous infection of a novel Chlamydia-like organism (CLO) with a Myxozoa parasite, Henneguya sp. in snakeskin gourami Trichopodus pectoralis in Thailand. A new CLO is proposed "Candidatus Piscichlamydia trichopodus" (CPT) based on 16S rRNA phylogenetic analysis. Systemic intracellular CPT infection was confirmed by histological examination, in situ hybridization, PCR assay, and sequencing of 16S rRNA. This novel pathogen belongs to the order Chlamydiales but differs in certain aspects from other species. The histopathological changes associated with CPT infection were different from the typical pathological lesions of epitheliocystis caused by previously known CLO. Unlike other CLO, CPT localized in the connective tissue rather than in the epithelial cells and formed smaller clumps of intracellular bacteria that stained dark blue with hematoxylin. On the other hand, typical myxospores of the genus Henneguya with tails were observed in the gill sections. Infection with Henneguya sp. resulted in extensive destruction of the gill filaments, most likely leading to respiratory distress. Due to the frequency of co-infections and the unavailability of culture methods for CLO and Henneguya sp, it was difficult to determine which pathogens were directly responsible for the associated mortality. However, co-infections may increase the negative impact on the host and the severity of the disease. Given the commercial importance of the snakeskin gourami and its significant aquaculture potential, the findings of this study are important for further studies on disease prevention. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. A novel nongenetic murine model of hyperglycemia and hyperlipidemia-associated aggravated atherosclerosis ObjectiveAtherosclerosis, the main pathology underlying cardiovascular diseases is accelerated in diabetic patients. Genetic mouse models require breeding efforts which are time-consuming and costly. Our aim was to establish a new nongenetic model of inducible metabolic risk factors that mimics hyperlipidemia, hyperglycemia, or both and allows the detection of phenotypic differences dependent on the metabolic stressor(s).

Methods and ResultsWild-type mice were injected with gain-of-function PCSK9D377Y (proprotein convertase subtilisin/kexin type 9) mutant adeno-associated viral particles (AAV) and streptozotocin and fed either a high-fat diet (HFD) for 12 or 20 weeks or a high-cholesterol/high-fat diet (Paigen diet, PD) for 8 weeks. To evaluate atherosclerosis, two different vascular sites (aortic sinus and the truncus of the brachiocephalic artery) were examined in the mice. Combined hyperlipidemic and hyperglycemic (HGHCi) mice fed a HFD or PD displayed characteristic features of aggravated atherosclerosis when compared to hyperlipidemia (HCi HFD or PD) mice alone. Atherosclerotic plaques of HGHCi HFD animals were larger, showed a less stable phenotype (measured by the increased necrotic core area, reduced fibrous cap thickness, and less -SMA-positive area) and had more inflammation (increased plasma IL-1{beta} level, aortic pro-inflammatory gene expression, and MOMA-2-positive cells in the BCA) after 20 weeks of HFD. Differences between the HGHCi and HCi HFD models were confirmed using RNA-seq analysis of aortic tissue, revealing that significantly more genes were dysregulated in mice with combined hyperlipidemia and hyperglycemia than in the hyperlipidemia-only group. The HGHCi-associated genes were related to pathways regulating inflammation (increased Cd68, iNos, and Tnfa expression) and extracellular matrix degradation (Adamts4 and Mmp14). When comparing HFD with PD, the PD aggravated atherosclerosis to a greater extent in mice and showed plaque formation after 8 weeks. Hyperlipidemic and hyperglycemic mice fed a PD (HGHCi PD) showed less collagen (Sirius red) and increased inflammation (CD68-positive cells) within aortic plaques than hyperlipidemic mice (HCi PD). HGHCi-PD mice represent a directly inducible hyperglycemic atherosclerosis model compared with HFD-fed mice, in which atherosclerosis is severe by 8 weeks.

ConclusionWe established a nongenetically inducible mouse model allowing comparative analyses of atherosclerosis in HCi and HGHCi conditions and its modification by diet, allowing analyses of multiple metabolic hits in mice. CryoEM of endogenous mammalian V-ATPase interacting with the TLDc protein mEAK-7 V-ATPases are rotary proton pumps that serve as signaling hubs with numerous proposed binding partners in cells. We used cryoEM combined with exhaustive focused classification to detect endogenous proteins that associate with V-ATPase from porcine kidney. A super-stoichiometric copy of subunit C was found in ~3% of complexes, while an additional ~1.6% of complexes bound mEAK-7, a protein with proposed roles in dauer formation in nematodes and mTOR signaling in mammals. High-resolution cryoEM of porcine kidney V-ATPase with recombinant mEAK-7 shows that mEAK-7s TLDc domain interacts with V-ATPases stator while its C-terminal  helix binds V-ATPases rotor. This crosslink would be expected to inhibit rotary catalysis. However, unlike inhibition of yeast V-ATPase by the TLDc protein Oxr1p, exogenous mEAK-7 does not inhibit V-ATPase and mEAK-7 overexpression in cells does not alter lysosomal or phagosomal pH. Instead, cryoEM suggests that interaction of mEAK-7 with V-ATPase is disrupted by ATP-induced rotation of the rotor. Comparison of Oxr1p and mEAK-7 binding explains this difference. Together, these results show that differences in V-ATPase binding by TLDc domain-containing proteins can lead to effects ranging from strong inhibition to formation of labile interactions that are sensitive to the enzymes activity. A northern range extension of a Canadian species of Special Concern, Dielis pilipes (Hymenoptera: Scoliidae), in the Okanagan Valley of British Columbia The only known Canadian records of the yellow scarab hunter wasp, Dielis pilipes (Saussure) (Hymenoptera: Scoliidae), are from the southern Okanagan and Similkameen valleys of British Columbia. We report a 25-km northern range extension of the species into the Ponderosa Pine Biogeoclimatic Ecosystem Classification zone, collected in an unmanaged agricultural field in Summerland, British Columbia. This finding is of conservation importance and has implications for natural biological control of ten-lined June beetles, Polyphylla decemlineata (Say) and P. crinita LeConte (Coleoptera: Scarabaeidae), incidental agricultural pests in the Okanagan. Naming of human diseases on the wrong side of history BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

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

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

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

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data. Migrators Within Migrators: Exploring Transposable Element Dynamics in the Monarch Butterfly, Danaus plexippus BackgroundLepidoptera (butterflies and moths) are an important model system in ecology and evolution. A high-quality chromosomal genome assembly is available for the monarch butterfly (Danaus plexippus), but it lacks an in-depth transposable element (TE) annotation, presenting an opportunity to explore monarch TE dynamics and the impact of TEs on shaping the monarch genome.

ResultsWe find 6.21% of the monarch genome is comprised of TEs, a reduction of 6.85% compared to the original TE annotation performed on the draft genome assembly. Monarch TE content is low compared to two closely related species with available genomes, Danaus chrysippus (33.97% TE) and Danaus melanippus (11.87% TE). The biggest TE contributions to genome size in the monarch are LINEs and Penelope-like elements, and three newly identified families, r2-hero_dPle (LINE), penelope-1_dPle (Penelope-like), and hase2-1_dPle (SINE), collectively contribute 34.92% of total TE content. We find evidence of recent TE activity, with two novel Tc1 families rapidly expanding over recent timescales (tc1-1_dPle, tc1- 2_dPle). LINE fragments show signatures of genomic deletions indicating a high rate of TE turnover. We investigate associations between TEs and wing colouration and immune genes and identify a three-fold increase in TE content around immune genes compared to other host genes.

ConclusionsWe provide a detailed TE annotation and analysis for the monarch genome, revealing a considerably smaller TE contribution to genome content compared to two closely related Danaus species with available genome assemblies. We identify highly successful novel DNA TE families rapidly expanding over recent timescales, and ongoing signatures of both TE expansion and removal highlight the dynamic nature of repeat content in the monarch genome. Our findings also suggest that insect immune genes are promising candidates for future interrogation of TE-mediated host adaptation. The initial response of females towards congeneric males matches the propensity to hybridise in Ophthalmotilapia. AO_SCPLOWBSTRACTC_SCPLOWCichlid radiations often harbour closely related species with overlapping niches and distribution ranges. Such species sometimes hybridise in nature, which raises the question how can they coexist. This also holds for the Tanganyika mouthbrooders Ophthalmotilapia ventralis and O. nasuta. Earlier studies found indications of asymmetrical hybridisation with females of O. ventralis accepting males of O. nasuta, but not the other way around. We hypothesised that this was due to differences in the capacity for species recognition. Given the higher propensity of O. ventralis females towards hybridisation, we expect a reduced ability for species recognition in O. ventralis females, compared to O. nasuta females. We staged two experiments, one focusing on 22 female O. nasuta and one on 21 female O. ventralis. These fish were placed in one half of a tank and briefly exposed to a conspecific or a heterospecific male, a conspecific female, or nothing (control). Female response was evaluated by scoring six tracking parameters and by noting the occurrence of ten discrete behaviours before and during the encounter. Females always responded to the presence of another fish by approaching it. Remarkably, for both O. nasuta and O. ventralis, we did not find a different response between encounters with conspecific males and females. However, in agreement with our hypothesis, O. nasuta females behaved differently towards conspecific or heterospecific males, whereas O. ventralis females did not. When presented with a heterospecific male, O. nasuta females performed a lower number of  ram behaviours. Additionally, they never displayed the  flee behaviour, a component of the species mating repertoire that was seen in all but one of the presentations with a conspecific male. Our findings show that differences in species recognition at first encounter predict to a large degree the outcome of the mating process, even in the absence of mating behaviour. 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. In vitro pectate lyase activity and carbon uptake assays and whole genome sequencing of Bacillus amyloliquefaciens subsp. plantarum strains for a pectin defective pathway The pectin lyase activity of 59 Bacillus amyloliquefaciens subsp. plantarum (Bap) strains was tested in vitro on Pectate Agar (PA) and Tris-Spizizen Salts (TSS) medium. Bap strains were cultured on TSA medium and washed three times with sterile water before the inoculation on PA media. Higher and lower pectate lyase activity were observed in six (AP193, AP203, AP299, AP80, AP102, and AP52) and four (AP 194, AP214, AP215, and AP305) Bap strains compared to other Bap strains. A total of 12 Bap strains (AP67, AP71, AP77, AP78, AP85, AP102, AP108, AP135, AP143, AP189, AP192, and AP193) grew vigorously on TSS medium. A total of six Bap strains (AP194, AP204, AP214, AP216, AP219, and HD73) had lower growth compared to other Bap strains. Pectin (1%) were used for in vitro PA and TSS medium. Pectate lyase and utilization activity were not found in Bacillus thuringiensis subsp. kurstaki strain HD73 compared to Bap strains. A draft genome sequence for strains AP194 and AP214 that were negative for pectin utilization were generated using an Illumina MiSeq. RAST analysis revealed that the pectin-associated gene altronate hydrolase (uxaA) absent in AP 214 strain. Multiple amino acid alignments of exuT and uxuB gene sequence showed dissimilarities among AP194, AP214, and reference Bap strains. From fitting the average to fitting the individual: A cautionary tale for mathematicalmodelers An outstanding challenge in the clinical care of cancer is moving from a one-size-fits-all approach that relies on population-level statistics towards personalized therapeutic design. Mathematical modeling is a powerful tool in treatment personalization, as it allows for the incorporation of patient-specific data so that treatment can be tailor-designed to the individual. Herein, we work with a mathematical model of murine cancer immunotherapy that has been previously-validated against the average of an experimental dataset. We ask the question: what happens if we try to use this same model to perform personalized fits, and therefore make individualized treatment recommendations? Typically, this would be done by choosing a single fitting methodology, and a single cost function, identifying the individualized best-fit parameters, and extrapolating from there to make personalized treatment recommendations. Our analyses show the potentially problematic nature of this approach, as predicted personalized treatment response proved to be sensitive to the fitting methodology utilized. We also demonstrate how a small amount of the right additional experimental measurements could go a long way to improve consistency in personalized fits. Finally, we show how quantifying the robustness of the average response could also help improve confidence in personalized treatment recommendations.

Author summaryAs we enter the era of healthcare where personalized medicine becomes a more common approach to treating cancer patients, harnessing the power of mathematical models will only become more essential. Using a preclinical dataset on cancer immunotherapy, we explore the challenges and limitations that arise when trying to move from fitting and making predictions for the population-level average, to fitting and making predictions for an individual. We find that the standard of approach of picking a single fitting methodology and a single cost function may end up having limited predictive value when applied to individual data. We also show how having a small amount of the right additional experimental data, and establishing the robustness of average treatment response, can help improve confidence in personalized model predictions. Differences in gait stability and acceleration characteristics between healthy young and older females Our aim was to evaluate differences in gait acceleration intensity, variability and stability of feet and trunk between older females and young females using inertial sensors. Twenty older females (OF; mean age 68.4, SD 4.1 years) and eighteen young females (YF; mean age 22.3, SD 1.7 years) were asked to walk straight for 100 meters at their preferred speed, while wearing inertial sensors on heels and lower back. We calculated spatiotemporal measures, foot and trunk acceleration characteristics and their variability, as well as trunk stability using the local divergence exponent (LDE). Two-way analysis of variance (including the factors foot and age), Students t-test, and Mann-Whitney U test were used to compare statistical differences of measures between groups. Cohens d effects were calculated for each variable. Foot maximum vertical acceleration and amplitude, trunk-foot vertical acceleration attenuation, as well as their variability were significantly smaller in OF than in YF. In contrast, trunk mediolateral acceleration amplitude, maximum vertical acceleration, and amplitude, as well as their variability were significantly larger in OF than in YF. Moreover, OF showed lower stability (i.e., higher LDE values) in mediolateral acceleration, mediolateral and vertical angular velocity of the trunk. Even though we measured healthy older females, these participants showed lower vertical foot accelerations with higher vertical trunk acceleration, lower trunk-foot vertical acceleration attenuation, less gait stability, and more variability of the trunk, and hence, were more likely to fall. These findings suggest that instrumented gait measurements may help for early detection of changes or impairments in gait performance, even before this can be observed by clinical eye or gait speed. Drug development of an affinity enhanced, broadly neutralizing heavy chain-only antibody that restricts SARS-CoV-2 in rodents We have identified camelid single-domain antibodies (VHHs) that cross-neutralize SARS-CoV-1 and -2, such as VHH72, which binds to a unique highly conserved epitope in the viral receptor-binding domain (RBD) that is difficult to access for human antibodies. Here, we establish a protein engineering path for how a stable, long-acting drug candidate can be generated out of such a VHH building block. When fused to human IgG1-Fc, the prototype VHH72 molecule prophylactically protects hamsters from SARS-CoV-2. In addition, we demonstrate that both systemic and intranasal application protects hACE-2-transgenic mice from SARS-CoV-2 induced lethal disease progression. To boost potency of the lead, we used structure-guided molecular modeling combined with rapid yeast-based Fc-fusion prototyping, resulting in the affinity-matured VHH72_S56A-Fc, with subnanomolar SARS-CoV-1 and -2 neutralizing potency. Upon humanization, VHH72_S56A was fused to a human IgG1 Fc with optimized manufacturing homogeneity and silenced effector functions for enhanced safety, and its stability as well as lack of off-target binding was extensively characterized. Therapeutic systemic administration of a low dose of VHH72_S56A-Fc antibodies strongly restricted replication of both original and D614G mutant variants of SARS-CoV-2 virus in hamsters, and minimized the development of lung damage. This work led to the selection of XVR011 for clinical development, a highly stable anti-COVID-19 biologic with excellent manufacturability. Additionally, we show that XVR011 is unaffected in its neutralizing capacity of currently rapidly spreading SARS-CoV-2 variants, and demonstrate its unique, wide scope of binding across the Sarbecovirus clades. Swine H1N1 influenza virus variants with enhanced polymerase activity and HA stability promote airborne transmission in ferrets Understanding how animal influenza A viruses (IAVs) acquire airborne transmissibility in humans and ferrets is needed to prepare for and respond to pandemics. Here, we investigated in ferrets the replication and transmission of swine H1N1 isolates P4 and G15, whose majority population had decreased polymerase activity and poor HA stability, respectively. For both isolates, a minor variant was selected and transmitted in ferrets. Polymerase-enhancing variant PA-S321 airborne-transmitted and propagated in one ferret. HA-stabilizing variant HA1-S210 was selected in all G15-inoculated ferrets and was transmitted by contact and airborne routes. With an efficient polymerase and a stable HA, the purified minor variant G15-HA1-S210 had earlier and higher peak titers in inoculated ferrets and was recovered at a higher frequency after airborne transmission than P4 and G15. Overall, HA stabilization played a more prominent role than polymerase enhancement in the replication and transmission of these viruses in ferrets. The results suggest pandemic risk-assessment studies may benefit from deep sequencing to identify minor variants with human-adapted traits.

IMPORTANCEDiverse IAVs circulate in animals, yet few acquire the viral traits needed to start a human pandemic. A stabilized HA and mammalian-adapted polymerase have been shown to promote the adaptation of IAVs to humans and ferrets (the gold-standard model for IAV replication, pathogenicity, and transmissibility). Here, we used swine IAV isolates of the gamma lineage a model to investigate the importance of HA stability and polymerase activity in promoting replication and transmission in ferrets. These are emerging viruses that bind to both -2,6- and -2,3-linked receptors. Using isolates containing mixed populations, a stabilized HA was selected within days in inoculated ferrets. An enhanced polymerase was also selected and propagated after airborne transmission to a ferret. Thus, HA stabilization was a stricter requirement, yet both traits promoted transmissibility. Knowing the viral traits needed for pandemic potential, and the relative importance of each, will help identify emerging viruses of greatest concern. Assessing sequence-based protein-protein interaction predictors for use in therapeutic peptide engineering Engineering peptides to achieve a desired therapeutic effect through the inhibition of a specific target activity or protein interaction is a non-trivial task. Few of the existing in silico peptide design algorithms generate target-specific peptides. Instead, many methods produce peptides that achieve a desired effect through an unknown mechanism. In contrast with resource-intensive high-throughput experiments, in silico screening is a cost-effective alternative that can prune the space of candidates when engineering target-specific peptides. Using a set of FDA-approved peptides we curated specifically for this task, we assess the applicability of several sequence-based protein-protein interaction predictors as a screening tool within the context of peptide therapeutic engineering. We show that similarity-based protein-protein interaction predictors are more suitable for this purpose than the state-of-the-art deep learning methods publicly available at the time of writing. We also show that this approach is mostly useful when designing new peptides against targets for which naturally-occurring interactors are already known, and that deploying it for de novo peptide engineering tasks may require gathering additional target-specific training data. Taken together, this work offers evidence that supports the use of similarity-based protein-protein interaction predictors for peptide therapeutic engineering, especially peptide analogs.