Revisiting regulatory decoherence and phenotypic integration: accounting for temporal bias in co-expression analyses Environment can alter the degree of phenotypic variation and covariation, potentially influencing evolutionary trajectories. However, environment-driven changes in phenotypic variation remain understudied. In an effort to exploit the abundance of RNASequencing data now available, an increasing number of ecological studies rely on population-level correlation to characterize the plastic response of the entire transcriptome and to identify environmentally responsive molecular pathways. These studies are fundamentally interested in identifying groups of genes that respond in concert to environmental shifts. We show that population-level differential co-expression exhibits biases when capturing changes of regulatory activity and strength in rice plants responding to elevated temperature. One possible cause of this bias is regulatory saturation, the observation that detectable co-variance between a regulator and its target may be low as their transcript abundances are induced. This phenomenon appears to be particularly acute for rapid-onset environmental stressors. However, our results suggest that temporal correlations may be a reliable means to detect transient regulatory activity following rapid onset environmental perturbations such as temperature stress. Such temporal bias is likely to confound the studies of phenotypic integration, where high-order organismal traits are hypothesized to be more integrated with strong correlation under stressful conditions, while recent transcriptome studies exhibited weaker coexpression between genes under stressful conditions. Collectively, our results point to the need to account for the nuances of molecular interactions and the possibly confounding effects that these can introduce into conventional approaches to study transcriptome datasets.
Quantitative analysis of gene expression on the inactive human X chromosome The human X chromosome exists in active (Xa) and inactive (Xi) states. To quantify Xi gene expression and explore its implications, we conducted RNA sequencing in two cell types from individuals with one Xa and zero to three Xis. Linear models fit the data well, indicating equal contributions from each Xi and quantifying significant Xi-driven expression changes for 124 of 406 expressed X-chromosome genes. As revealed by allele-specific expression analyses, these Xi-driven changes reflect both i) Xi transcription and ii) Xi modulation - both positive and negative - of Xa transcript levels. Incorporating population and evolutionary measures of constraint on expression levels, we identified 13 X-chromosome genes most likely to drive the phenotypic consequences of natural variation in Xi copy number. These analytics and the underlying datasets offer a foundation for understanding Turner, Klinefelter, and other sex chromosome aneuploidy syndromes, and sex-chromosome effects on health and disease in euploid males and females.
Critical re-evaluation of the slope factor of the operational model of agonism Background and PurposeAlthough being a relative term, agonist efficacy is a cornerstone in the proper assessment of agonist selectivity and signalling bias. The operational model of agonism (OMA) has become successful in the determination of agonist efficacies and ranking them. In 1985, Black et al. introduced the slope factor to the OMA to make it more flexible and allow for fitting steep as well as flat concentration-response curves.

Experimental ApproachFunctional analysis of OMA demonstrates that the slope factor implemented by Black et al. affects relationships among parameters of the OMA, leading to wrong estimates of operational efficacy and affinity of an agonist. In contrast, functional analysis of OMA with Hill coefficient as a slope factor has not this flaw.

Key ResultsFitting of the OMA with Black et al. slope factor to concentration-response curves of several models (homotropic allosteric modulation, substrate inhibition and non-competitive auto-inhibition) gives wrong estimates of operational efficacy and affinity. In contrast, fitting of the OMA modified by the Hill coefficient to concentration-response curves of these several models gives correct estimates of operational efficacy and affinity.

Conclusion and ImplicationsThe OMA modified by the Hill coefficient can be used for correct ranking of agonism and subsequent analysis, like quantification of signalling bias. An Essential Role of UBXN3B in B Lymphopoiesis Hematopoiesis is finely regulated to enable timely production of the right numbers and types of mature immune cells to maintain tissue homeostasis. Dysregulated hematopoiesis may compromise antiviral immunity and/or exacerbate immunopathogenesis. Herein, we report an essential role of UBXN3B in maintenance of hematopoietic homeostasis and restriction of immunopathogenesis during respiratory viral infection. Ubxn3b deficient (Ubxn3b-/-) mice are highly vulnerable to SARS-CoV-2 and influenza A infection, characterized by more severe lung immunopathology, lower virus-specific IgG, significantly fewer B cells, but more myeloid cells than Ubxn3b+/+ littermates. This aberrant immune compartmentalization is recapitulated in uninfected Ubxn3b-/- mice. Mechanistically, UBXN3B controls precursor B-I (pre-BI) transition to pre-BII and subsequent proliferation in a cell-intrinsic manner, by maintaining BLNK protein stability and pre-BCR signaling. These results reveal an essential role of UBXN3B for the early stage of B cell development. The Digital 3D-Atlas MAKER (DAMAKER): a dynamic and expandable digital 3D-tool for monitoring the temporal changes in tissue growth during hindbrain morphogenesis Reconstruction of prototypic three-dimensional (3D) atlases at the scale of whole tissues or organs requires specific methods to be developed. We have established a digital 3D-atlas maker (DAMAKER) and built a digital 3D-temporal atlas to monitor the changes in the growth of the tissue and neuronal differentiation domain in the zebrafish hindbrain. DAMAKER integrates spatial and temporal data from cell populations, neuronal differentiation and brain morphogenesis, through in vivo imaging techniques paired with image analyses and segmentation tools. First, we generated a 3D-reference from several imaged hindbrains and segmented them using a trainable tool; these were aligned using rigid registration, revealing distribution of neuronal differentiation growth patterns along the axes. Second, we quantified the dynamic growth of the neuronal differentiation domain vs. the progenitor domain, and by in vivo neuronal birthdating experiments we generated a digital 3D-temporal map of the neuronal growth in the whole hindbrain, revealing the spatiotemporal dynamics of neuronal differentiation upon morphogenesis. Last, we applied it to glutamatergic and GABAergic neurons, as proof-of-concept that the digital 3D-temporal map could be used as a proxy to infer neuronal birthdate. As this protocol uses open-access tools and algorithms, it can be shared for standardized, accessible, tissue-wide cell population atlas construction. The activity and functions of soil microbial communities in the Finnish sub-Arctic vary across vegetation types Increased microbial activity in high-latitude soils due to climate change might lead to higher greenhouse gas (GHG) emissions. However, mechanisms of microbial GHG production and consumption in tundra soils are not thoroughly understood. We analyzed 116 soil metatranscriptomes from 73 sites in the Finnish sub-Arctic to investigate how the diversity and functional potential of bacterial and archaeal communities vary across vegetation types and soil layers. Soils differed in physicochemical conditions, with meadow soils being characterized by higher pH and low soil organic matter (SOM) and carbon/nitrogen ratio whereas dwarf shrub-dominated ecosystems with high SOM and low pH. Actinobacteria, Acidobacteria, Alphaproteobacteria, and Planctomycetes predominated all communities but there were significant differences on genus level between vegetation types, as plant polymer degrading groups were more active in shrub-dominated soils compared to meadows. Given that climate change scenarios predict expansion in dwarf shrubs at high latitudes, our results indicate that the rate of carbon turnover in tundra soils may increase in the future. Additionally, transcripts of methanotrophs were detected in the mineral layer of all soils, potentially moderating methane fluxes from deeper layers. In all, this study provides new insights into possible shifts in tundra microbial diversity and activity with climate change. 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 interplay of chromatin phase separation and lamina interactions in nuclear organization The genetic material of eukaryotes is segregated into transcriptionally active euchromatin and silent heterochromatin compartments. The spatial arrangement of chromatin compartments evolves over the course of cellular life in a process that remains poorly understood. The latest nuclear imaging experiments reveal a number of dynamical signatures of chromatin which are reminiscent of active multi-phase liquids. This includes the observations of viscoelastic response, coherent motions, Ostwald ripening and coalescence of chromatin compartments. There is also growing evidence that liquid-liquid phase separation of protein and nucleic acid components is the underlying mechanism for the dynamical behavior of chromatin. In order to dissect the organizational and dynamical implications of chromatins liquid behavior, we have devised a phenomenological field-theoretic model of nucleus as a multi-phase condensate of liquid chromatin types. Employing the liquid chromatin model of Drosophila nucleus, we have carried out an extensive set of simulations with an objective to shed light on the dynamics and chromatin patterning observed in the latest nuclear imaging experiments. Our simulations reveal the emergence of experimentally detected mesoscale chromatin channels and spheroidal droplets which arise from the dynamic interplay of chromatin type to type interactions and intermingling of chromosomal territories. We also quantitatively reproduce coherent motions of chromatin domains observed in displacement correlation spectroscopy measurements which are explained within the framework of our model by phase separation of chromatin types operating within constrained intra-chromosomal and inter-chromosomal boundaries. Finally, we illuminate the role of heterochromatin-lamina interactions in the nuclear organization by showing that these interactions enhance the mobility of euchromatin and indirectly introduce correlated motions of heterochromatin droplets.

SIGNIFICANCE STATEMENTThe latest super-resolution imaging experiments have revealed a surprisingly dynamic and stochastic nature of chromatin in eukaryotic nuclei which is reminiscent of multi-phase fluid behavior. As a result, the concepts from the theory of complex fluids such as phase separation, viscoelasticity, and droplet nucleation have found widespread utility in understanding salient features of nuclear organization. In order to understand and disentangle the complex interplay of forces that contribute to the emergent patterns of organization and dynamics, we have devised a phenomenological field-theoretic model of nucleus as a multi-phase condensate of liquid chromatin types. Armed with a mesoscopic model of nuclear chromatin, we have shed light on the distinct dynamical and structural contributions of chromatin type interactions intermingling of chromosomal territories and lamina binding. We also shed light on the dynamical heterogeneity and coherent motions of chromatin domains which are fully captured by an interplay of micro-phase separation of chromatin types and lamina binding. Enhancing CDK4/6 inhibitor therapy for medulloblastoma using nanoparticle delivery and scRNA-seq-guided combination with sapanisertib CDK4/6 inhibitors hold promise for brain tumor treatment, but efficacy has been limited by recurrence in both preclinical models and clinical trials. To address recurrence, we tested a nanoparticle formulation of the CDK4/6 inhibitor palbociclib (POx-palbo) in mice genetically-engineered to develop SHH-driven medulloblastoma. We then analyzed medulloblastomas in mice receiving palbociclib treatment, and compared the efficacy of combining palbociclib with specific inhibitors suggested by our analysis. POx-Palbo showed reduced toxicity compared to conventional palbociclib, was tolerable in parenteral administration, improved CNS pharmacokinetics, and extended survival of mice with medulloblastoma. Recurrence, however, remained problematic as fractions of tumor cells proliferated during therapy. ScRNA-seq identified a gene expression pattern unique to proliferating medulloblastoma cells in POx-Palbo-treated mice, marked by up-regulation of the glutamate transporter Slc1a2 and down-regulation of diverse ribosomal genes. Reduced mTORC1 signaling, suggested by ribosomal suppression in POx-Palbo-treated tumors was confirmed by decreased 4EBP1 phosphorylation (p4EBP1). Further reducing mTORC1 activity by combining POx-Palbo with the mTORC1 inhibitor sapanisertib produced mutually enhancing effects, with increased suppression of both pRB and p4EBP1, and prolonged mouse survival compared to either agent alone. In contrast, targeting cell cycle progression by combining POx-Palbo with the SHH-pathway inhibitor vismodegib, or with the replication-targeting agents gemcitabine or etoposide, failed to enhance efficacy. Our data show the potential of nanoparticle formulation and scRNA-seq analysis of resistance to improve brain tumor treatment, and identify POx-palbo plus sapanisertib as effective combinatorial therapy for SHH medulloblastoma. This combination may be appropriate for testing in patients with recurrence, who need new options. Computation of the mitochondrial age distribution along the axon length We describe a compartmental model of mitochondrial transport in axons, which we apply to compute mitochondrial age at different distances from the soma. The model predicts that the average mitochondrial age at the tip of an axon with a length of 1 cm is approximately 22 hours. The mitochondria are youngest closest to the soma and their age scales approximately linearly with distance from the soma. To the best of the authors knowledge, this is the first attempt to predict the spatial distribution of mitochondrial age within an axon. A sensitivity study of the mean age of mitochondria to various model parameters is also presented. Differential regulation of degradation and immune pathways underlies adaptation of the ectosymbiotic nematode Laxus oneistus to oxic-anoxic interfaces Eukaryotes may experience oxygen deprivation under both physiological and pathological conditions. Because oxygen shortage leads to a reduction in cellular energy production, all eukaryotes studied so far conserve energy by suppressing their metabolism. However, the molecular physiology of animals that naturally and repeatedly experience anoxia is underexplored. One such animal is the marine nematode Laxus oneistus. It thrives, invariably coated by its sulfur-oxidizing symbiont Candidatus Thiosymbion oneisti, in anoxic sulfidic or hypoxic sand. Here, transcriptomics and proteomics showed that, whether in anoxia or not, L. oneistus mostly expressed genes involved in ubiquitination, energy generation, oxidative stress response, immune response, development, and translation. Importantly, ubiquitination genes were also upregulated when the nematode was subjected to anoxic sulfidic conditions, together with genes involved in autophagy, detoxification and ribosome biogenesis. We hypothesize that these degradation pathways were induced to recycle damaged cellular components (mitochondria) and misfolded proteins into nutrients. Remarkably, when L. oneistus was subjected to anoxic sulfidic conditions, lectin and mucin genes were also upregulated, potentially to promote the attachment of its thiotrophic symbiont. Furthermore, the nematode appeared to survive oxygen deprivation by using an alternative electron carrier (rhodoquinone) and acceptor (fumarate), to rewire the electron transfer chain. On the other hand, under hypoxia, genes involved in costly processes (e.g., amino acid biosynthesis, development, feeding, mating) were upregulated, together with the worms Toll- like innate immunity pathway and several immune effectors (e.g., Bacterial Permeability Increasing proteins, fungicides).

In conclusion, we hypothesize that, in anoxic sulfidic sand, L. oneistus upregulates degradation processes, rewires oxidative phosphorylation and by reinforces its coat of bacterial sulfur-oxidizers. In upper sand layers, instead, it appears to produce broad-range antimicrobials and to exploit oxygen for biosynthesis and development. Community Evaluation of Glycoproteomics Informatics Solutions Reveals High-Performance Search Strategies of Serum N- and O-Glycopeptide Data Glycoproteome profiling (glycoproteomics) is a powerful yet analytically challenging research tool. The complex tandem mass spectra generated from glycopeptide mixtures require sophisticated analysis pipelines for structural determination. Diverse software aiding the process have appeared, but their relative performance remains untested. Conducted through the HUPO Human Proteome Project - Human Glycoproteomics Initiative, this community study, comprising both developers and users of glycoproteomics software, evaluates the performance of informatics solutions for system-wide glycopeptide analysis. Mass spectrometry-based glycoproteomics datasets from human serum were shared with all teams. The relative team performance for N- and O-glycopeptide data analysis was comprehensively established and validated through orthogonal performance tests. Excitingly, several high-performance glycoproteomics informatics solutions were identified. While the study illustrated that significant informatics challenges remain, as indicated by a high discordance between annotated glycopeptides, lists of high-confidence (consensus) glycopeptides were compiled from the standardised team reports. Deep analysis of the performance data revealed key performance-associated search variables and led to recommendations for improved "high coverage" and "high accuracy" glycoproteomics search strategies. This study concludes that diverse software for comprehensive glycopeptide data analysis exist, points to several high-performance search strategies, and specifies key variables that may guide future software developments and assist informatics decision-making in glycoproteomics. 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. Full-likelihood genomic analysis clarifies a complex history of species divergence and introgression: the example of the erato-sara group of Heliconius butterflies AO_SCPLOWBSTRACTC_SCPLOWIntrogressive hybridization plays a key role in adaptive evolution and species diversification in many groups of species. However, frequent hybridization and gene flow between species makes estimation of the species phylogeny and key population parameters challenging. Here, we show that by accounting for phasing and using full-likelihood analysis methods, introgression histories and population parameters can be estimated reliably from whole-genome sequence data. We employ full-likelihood methods under the multispecies coalescent (MSC) model with and without gene flow to analyze the genomic data from six members of the erato-sara clade of Heliconius butterflies and infer the species phylogeny and cross-species introgression events. The methods naturally accommodate random fluctuations in genealogical history across the genome due to deep coalescence. To avoid heterozygote phasing errors in haploid sequences commonly produced by genome assembly methods, we process and compile unphased diploid sequence alignments and use analytical methods to average over uncertainties in heterozygote phase resolution. There is robust evidence for introgression across the genome, both among distantly related species deep in the phylogeny and between sister species in shallow parts of the tree. We obtain chromosome-specific estimates of key population parameters such as introgression directions, times and probabilities, as well as species divergence times and population sizes for modern and ancestral species. We confirm ancestral gene flow between the sara clade and an ancestral population of H. telesiphe, a likely hybrid speciation origin for H. hecalesia, and gene flow between sister species H. erato and H. himera. Inferred introgression among ancestral species also explains the history of two chromosomal inversions deep in the phylogeny of the group. This study illustrates how a full-likelihood approach based on the multispecies coalescent makes it possible to extract rich historical information of species divergence and gene flow from genomic data. Feed efficiency of lactating Holstein cows was not as repeatable across diets as within diet over subsequent lactation stages BackgroundImproving feed efficiency has become a common target for dairy farmers to meet the requirement of producing more milk with fewer resources. To improve feed efficiency, a prerequisite is to ensure that the cows identified as most or least efficient will remain as such, independently of diet composition. Therefore, the current research analysed the ability of lactating dairy cows to maintain their feed efficiency while changing the energy density of the diet by changing its concentration in starch and fibre. A total of 60 lactating Holstein cows, including 33 primiparous cows, were first fed a high starch diet-low fibre (diet S+F-), then switched over to a low starch diet-high fibre (diet S-F+). To know if diet affect feed efficiency, we compared the ability of feed efficiency to be maintained within a diet over subsequent lactation stages, known as repeatability of feed efficiency, with its ability to be maintained across diets, known as reproducibility of feed efficiency. To do so we used two indicators: the estimation of the error of repeatability/reproducibility, which is commonly used in metrology, and the coefficient of correlation of concordance (CCC), which is used in biology. The effect of diet change could also lead to a change in cows sorting behaviour which could potentially affect feed efficiency if for example the most efficient cows select more concentrate than the least efficient. We therefore analysed the relationship between the differences in individual feed refusals composition and the differences in feed efficiency. To do so, the composition of each feed refusal was described with its near infra-red (NIR) spectroscopy and was performed on each individual feed ingredient, diet and refusals and used as composition variable. The variability of the NIR spectra of the refusals was described with its principal components thanks to a principal component analysis (PCA). The Pearson correlation was estimated to check the relationship between feed efficiency and refusals composition, i.e. sorting behaviour.

ResultsThe error of reproducibility of feed efficiency across diets was 2.95 MJ/d. This error was significantly larger than the errors of repeatability estimated within diet, which were 2.01 MJ/d within diet S-F+ and 2.40 MJ/d within diet S+F-. The CCC was 0.64 between feed efficiency estimated within diet S+F- and feed efficiency estimated within diet S-F+. This CCC was smaller than the one observed for feed efficiency estimated within diet between two subsequent lactation stages (CCC = 0.72 within diet S+F- and 0.85 within diet S-F+). Feed efficiency was poorly correlated to the first two principal components, which explained 90% of the total variability of the NIR spectra of the individual refusals. This suggests that feed sorting behaviour did not explain differences in feed efficiency.

ConclusionsFeed efficiency was significantly less reproducible across diets than repeatable within the same diet over subsequent lactation stages, but cows ranking for feed efficiency was not significantly affected by diet change. This loss in repeatability across diets could be due to a more pronounced feed sorting subsequent to the change in diet composition. However, the differences in sorting behaviour between cows were not associated to feed efficiency differences in this trial neither with the S+F- diet nor with the S-F+ diet. Those results have to be confirmed on diets having different forage to concentrate ratios to ensure that the least and most efficient cows will not change. Principal component analysis reveals multiple consistent responses to naturalistic stimuli in children and adults Functional MRI (fMRI) study of naturalistic conditions, e.g. movie watching, usually focuses on shared responses across subjects. However, individual differences have been attracting increasing attention in search of group differences or associations with behavioral outcomes. Individual differences are typically studied by directly modeling the pair-wise intersubject correlation matrix or projecting the relations onto a single dimension. We contend that it is critical to examine whether there are one or more consistent responses underlying the whole sample, because multiple components, if exist, may undermine the intersubject relations using the previous methods. We propose to use principal component analysis (PCA) to examine the heterogeneity of brain responses across subjects and project the individual variability into higher dimensions. By analyzing an fMRI dataset of children and adults watching a cartoon movie, we showed evidence of two consistent responses in the supramarginal gyrus and other regions. While the first components in many regions represented a response pattern mostly in older children and adults, the second components mainly represented the younger children. The second components in the supramarginal network resembled a delayed version of the first PCs for 4 seconds (2 TR), indicating slower responses in the younger children than the older children and adults. The analyses highlight the importance of identifying multiple consistent responses in responses to naturalistic stimuli. This PCA-based approach could be complementary to the commonly used intersubject correlation to analyze movie watching data. Robust retrieval of data stored in DNA by de Bruijn graph-based de novo strand assembly DNA data storage is a rapidly developing technology with great potential due to its high density, long-term durability, and low maintenance cost. The major technical challenges include various errors, such as the strand breaks, rearrangements, and indels that frequently arise during DNA synthesis, amplification, sequencing, and preservation. Through a de novo assembly strategy, we developed an algorithm based on the de Bruijn graph and greedy path search (DBGPS) to address these issues. DBGPS shows distinct advantages in handling DNA breaks, rearrangements, and indels. The robustness of DBGPS is demonstrated by accelerated aging, multiple independent data retrievals, deep error-prone PCR, and large data scale simulations. Remarkably, 6.8 MB of data can be retrieved accurately from a seriously corrupted sample that has been treated at 70 {degrees}C for 70 days. With DBGPS, we were able to achieve a logical density of 1.30 bits/cycle and a physical density of 295 PB/g.

One-Sentence SummaryA de Bruijn graph-based de novo assembly algorithm for DNA data storage enables fast and robust data readouts even with DNA samples that have been severely corrupted. Manduca sexta experience high parasitoid pressures in the field but minor fitness costs of consuming plant secondary compounds O_LIPlant-herbivore co-evolutionary interactions have led to a range of plant defenses that minimize insect damage and a suite of counter-adaptations that allow herbivores to feed on defended plants. Consuming plant secondary compounds results in herbivore growth and developmental costs but can have beneficial effects such as deterrence or harm of parasitoid enemies. Therefore, the role of secondary compounds on herbivore fitness must be considered in the context of the abundance and level of harm from natural enemies and the costs herbivores incur feeding on plant secondary compounds.
C_LIO_LIIn this study, I combined field measurements of Cotesia congregata wasp parasitism pressure with detailed measurements of the costs of plant secondary compounds across developmental stages in the herbivore host, Manduca sexta.
C_LIO_LII show that C. congregata parasitoids exert large negative selective pressures, killing 31-57% of M. sexta larvae in the field. Manduca sexta developed fastest during instars most at risk for parasitoid oviposition but growth was slowed by consumption of plant secondary compounds. The negative effects of consuming plant secondary compounds as larvae influenced adult size traits but there were no immune, survival, or fecundity costs.
C_LIO_LIThese results suggest that developmental costs experienced by M. sexta herbivores consuming defensive compounds are minor in comparison to the strong negative survival pressures from abundant parasitoid enemies.
C_LI Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. Using Hierarchical Similarity To Examine The Genetics of Behcet's disease Behcets disease (BD) is a multisystem inflammatory disease that affects patients along the historic silk road. Thus far, the pathogenesis of the disease has proved elusive due to the complex genetic interactions and unknown environmental or viral triggering factors of the disease. In this paper, we seek to clarify the genetic factors of the disease while also uncovering other diseases of interest that present with a similar genotype as BD. To do this, we employ a convergent functional genomics approach by leveraging the hierarchical similarity tool available in Geneweaver. Through our analysis, we were able to ascertain 7 BD consensus genes and 16 autoimmune diseases with genetic overlap with BD. The results of our study will inform further research into the pathogenesis of Behcets Disease. Cryptic fertility costs of viral infection in a model social insect Declining insect populations emphasize the importance of understanding the drivers underlying reductions in insect fitness. Here, we investigated viruses as a threat to social insect reproduction, using honey bees as a model species. We report that in a sample of N = 93 honey bee (Apis mellifera) queens from nine beekeeping operations across a wide geographic range, high levels of natural viral infection are associated with decreased ovary mass. We confirmed this finding in an independent sample of N = 54 queens. Failed (poor quality) queens displayed higher levels of viral infection, reduced sperm viability, smaller ovaries, and altered ovary protein composition compared to healthy queens. We experimentally infected queens with Israeli acute paralysis virus (IAPV) and found that the ovary masses of IAPV-injected queens were significantly smaller than control queens, demonstrating a causal relationship between viral infection and ovary size. Queens injected with IAPV also had significantly lower expression of vitellogenin, the main source of nutrition deposited into developing oocytes, and higher levels of heat-shock proteins (HSPs), which are part of the honey bees antiviral response. This work together shows that viral infections occurring naturally in the field are compromising queen reproductive success. TSC22D4 interacts with Akt1 in response to metabolic and stress signals Transforming Growth Factor {beta} 1 Stimulated Clone 22 D4 (TSC22D4) is an intrinsically disordered protein that regulates cellular and physiological processes such as cell proliferation, cellular senescence as well as hepatic glucose and lipid metabolism. The molecular mechanism of TSC22D4 action in these cellular and metabolic functions, however, remains largely elusive. Here, we identified TSC22D4 as a novel protein kinase B/Akt1 interacting protein, a critical mediator of insulin/PI3K signaling pathway implicated in diverse set of diseases including type 2 diabetes, obesity and cancer. TSC22D4 interacts with Akt1 not constitutively but rather in a regulatory manner. While glucose and insulin stimulation of cells or refeeding of mice impair the hepatic TSC22D4-Akt1 interaction, inhibition of mitochondria and oxidative stress, promote it; indicating that extra- and intra-cellular cues play a key role in controlling TSC22D4-Akt1 interaction. Our results also demonstrate that together with its dimerization domain, i.e. the TSC box, TSC22D4 requires its intrinsically disordered region (D2 domain) to interact with Akt1. To understand regulation of TSC22D4 function further, we employed tandem mass spectrometry and identified 15 novel phosphorylation sites on TSC22D4. Similar to TSC22D4-Akt1 interaction, TSC22D4 phosphorylation also responds to environmental signals such as starvation, mitochondrial inhibition and oxidative stress. Interestingly, 6 out of the 15 novel phosphorylation sites lie within the TSC22D4 D2 domain, which is required for TSC22D4-Akt1 interaction. Characterization of the regulation and function of these novel phosphorylation sites, in the future, will shed light on our understanding of the role of TSC22D4-Akt1 interaction in both cell biological and physiological functions. Overall, our findings postulate a model whereby TSC22D4 acts as an environmental sensor and interacts with Akt1 to regulate cell proliferation, cellular senescence as well as maintain metabolic homeostasis. The extracellular matrix controls stem cell specification and crypt morphology in the developing and adult gut The rapid renewal of the epithelial gut lining is fuelled by stem cells that reside at the base of intestinal crypts. In recent years, the signal transduction pathways and morphogens that regulate intestinal stem cell self-renewal and differentiation have been extensively characterised. In contrast, although extracellular matrix (ECM) components form an integral part of the intestinal stem cell niche, their direct influence on the cellular composition is less well understood. Here, we set out to systematically compare the effect of two major ECM classes, fibrillar collagen type I and the basement membrane, on the intestinal epithelium. We found that both collagen I and laminin-containing cultures allow growth of small intestinal epithelial cells with all differentiated and undifferentiated cell types present in both cultures, albeit at different ratios. Specific to the collagen culture was a subset of cells with a fetal-like gene expression program. In contrast, laminin, but not collagen IV, increased Lgr5+ stem cells and Paneth cells, and induced crypt-like morphology changes. The transition from a collagen culture to a laminin culture, resembles the gut development in vivo. Here, the ECM is dramatically remodelled by mesenchymal cells, which is accompanied by a specific and local expression of the laminin receptor ITGA6 in the crypt-forming epithelium. This laminin:ITGA6 signalling is essential for the stem cell induction and crypt formation in vitro. Importantly, deletion of laminin in the adult mouse results in a fetal-like epithelium with a marked reduction of adult intestinal stem cells. Overall, our data support the hypothesis that the formation of intestinal crypts is induced by an increased laminin concentration in the ECM. Long term analysis of social structure: evidence of age-based consistent associations in male Alpine ibex Despite its recognized importance for understanding the evolution of animal sociality as well as for conservation, long term analysis of social networks of animal populations is still relatively uncommon. We investigated social network dynamics in males of a gregarious mountain ungulate (Alpine ibex, Capra ibex) over ten years focusing on groups, sub-groups and individuals, exploring the dynamics of sociality over different scales. Despite the social structure changing between seasons, the Alpine ibex population was highly cohesive: fission-fusion dynamics lead almost every male in the population to associate with each other male at least once. Nevertheless, we found that male Alpine ibex showed preferential associations that were maintained across seasons and years. Age seemed to be the most important factor driving preferential associations while other characteristics, such as social status, appeared less crucial. We also found that centrality measures were influenced by age and were also related to individual physical condition. The multi-scale and long-term frame of our study helped us show that ecological constrains, such as resource availability, may play a role in shaping associations in a gregarious species, but they cannot solely explain sociality and preferential association that are likely also to be driven by life-history linked physiological and social needs. Our results highlight the importance of long-term studies based on individually recognizable subjects to help us build on our understanding of the evolution of animal sociality. Genetic variation at mouse and human ribosomal DNA influences associated epigenetic states BackgroundRibosomal DNA (rDNA) displays substantial inter-individual genetic variation in human and mouse. A systematic analysis of how this variation impacts epigenetic states and expression of the rDNA has thus far not been performed.

ResultsUsing a combination of long- and short-read sequencing, we establish that 45S rDNA units in the C57BL/6J mouse strain exist as distinct genetic haplotypes that influence the epigenetic state and transcriptional output of any given unit. DNA methylation dynamics at these haplotypes are dichotomous and life-stage specific: at one haplotype, the DNA methylation state is sensitive to the in utero environment, but refractory to post-weaning influences, whereas other haplotypes entropically gain DNA methylation during ageing only. On the other hand, individual rDNA units in human show limited evidence of genetic haplotypes, and hence little discernible correlation between genetic and epigenetic states. However, in both species, adjacent units show similar epigenetic profiles, and the overall epigenetic state at rDNA is strongly positively correlated with total rDNA copy number. Analysis of different mouse inbred strains reveals that in some strains, such as 129S1/SvImJ, rDNA copy number is only approximately 150 copies per diploid genome and DNA methylation levels are <5%.

ConclusionsOur work demonstrates that rDNA-associated genetic variation has a considerable influence on rDNA epigenetic state and consequently rRNA expression outcomes. In the future, it will be important to consider the impact of inter-individual rDNA (epi)genetic variation on mammalian phenotypes and diseases. The desert green algae Chlorella ohadii thrives at excessively high light intensities by exceptionally enhancing the mechanisms that protect photosynthesis from photoinhibition Although light is the driving force of photosynthesis, excessive light can be harmful. One of the main processes that limits photosynthesis is photoinhibition, the process of light-induced photodamage. When the absorbed light exceeds the amount that is dissipated by photosynthetic electron flow and other processes, damaging radicals are formed that mostly inactivate photosystem II (PSII). Damaged PSII must be replaced by a newly repaired complex in order to preserve full photosynthetic activity. Chlorella ohadii is a green micro-alga, isolated from biological desert soil crusts, that thrives under extreme high light and is highly resistant to photoinhibition. Therefore, C. ohadii is an ideal model for studying the molecular mechanisms underlying protection against photoinhibition. Comparison of the thylakoids of C. ohadii cells that were grown under low light versus extreme high light intensities, found that the alga employs all three known photoinhibition protection mechanisms: i) massive reduction of the PSII antenna size; ii) accumulation of protective carotenoids; and iii) very rapid repair of photo-damaged reaction center proteins. This work elucidated the molecular mechanisms of photoinhibition resistance in one of the most light-tolerant photosynthetic organisms and shows how photoinhibition protection mechanisms evolved to marginal conditions, enabling photosynthesis-dependent life in severe habitats.

One Sentence HighlightAnalysis of the photosynthetic properties of a desert algae that thrives at extreme high light intensities revealed protection from photoinhibition driven by the remarkable enhancement of three protection mechanisms.