NECAB2 participates in an endosomal pathway of mitochondrial stress response at striatal synapses Synaptic signaling depends on ATP generated by mitochondria. Due to extensive connectivity, the striatum is especially vulnerable to mitochondrial dysfunction and thus requires efficient mitochondrial quality control and repair. We found that global knockout of the neuronal calcium-binding protein 2 (NECAB2) in the mouse causes loss of striatal synapses and behavioral phenotypes related to striatal dysfunction such as reduced motivation and altered sensory gating. Striatal mitochondria from Necab2 knockout mice are more abundant and smaller. They are characterized by increased respiration and superoxide production resulting in oxidative stress. This accumulation of dysfunctional mitochondria is caused by a defective assembly of mitochondria with early endosomes in a pathway that involves the small GTPase Rab5 and its guanine nucleotide exchange factor Alsin/ALS2. NECAB2 therefore participates in an endosomal pathway of mitochondrial stress response and repair important for striatal function.
GWAS-associated variants, non-genetic factors, and transient transcriptome in Multiple Sclerosis etiopathogenesis: a colocalization analysis A clinically actionable understanding of multiple sclerosis (MS) etiology goes through GWAS interpretation, prompting research on new gene regulatory models. We previously suggested a stochastic etiologic model where small-scale random perturbations could reach a threshold for MS development. The recently described mapping of the transient transcriptome (TT), including intergenic and intronic RNAs, seems appropriate to verify this model through a rigorous colocalization analysis. We show that genomic regions coding for the TT were significantly enriched for MS-associated GWAS variants and DNA binding sites for molecular transducers mediating putative, non-genetic, etiopathogenetic factors for MS (e.g., vitamin D deficiency, Epstein Barr virus latent infection, B cell dysfunction). These results suggest a model whereby TT-coding regions are hotspots of convergence between genetic ad non-genetic factors of risk/protection for MS, and plausibly for other complex disorders. Our colocalization analysis also provides a freely available data resource (www.mscoloc.com) for future research on MS transcriptional regulation.
Genome-scale metabolic modeling reveals SARS-CoV-2-induced metabolic changes and antiviral targets Tremendous progress has been made to control the COVID-19 pandemic caused by the SARS-CoV-2 virus. However, effective therapeutic options are still rare. Drug repurposing and combination represent practical strategies to address this urgent unmet medical need. Viruses, including coronaviruses, are known to hijack host metabolism to facilitate viral proliferation, making targeting host metabolism a promising antiviral approach. Here, we describe an integrated analysis of 12 published in vitro and human patient gene expression datasets on SARS-CoV-2 infection using genome-scale metabolic modeling (GEM), revealing complicated host metabolism reprogramming during SARS-CoV-2 infection. We next applied the GEM-based metabolic transformation algorithm to predict anti-SARS-CoV-2 targets that counteract the virus-induced metabolic changes. We successfully validated these targets using published drug and genetic screen data and by performing an siRNA assay in Caco-2 cells. Further generating and analyzing RNA-sequencing data of remdesivir-treated Vero E6 cell samples, we predicted metabolic targets acting in combination with remdesivir, an approved anti-SARS-CoV-2 drug. Our study provides clinical data-supported candidate anti-SARS-CoV-2 targets for future evaluation, demonstrating host metabolism-targeting as a promising antiviral strategy. Stress fiber strain recognition by the LIM protein testin is cryptic and mediated by RhoA The actin cytoskeleton is a key regulator of mechanical processes in cells. The family of LIM domain proteins have recently emerged as important mechanoresponsive cytoskeletal elements capable of sensing strain in the actin cytoskeleton. The mechanisms regulating this mechanosensitive behavior, however, remain poorly understood. Here we show that the LIM domain protein testin is peculiar in that despite the full-length protein primarily appearing diffuse in the cytoplasm, the C-terminal LIM domains alone recognize focal adhesions and strained actin while the N-terminal domains alone recognize stress fibers. Phosphorylation mutations in the dimerization regions of testin, however, reveal its mechanosensitivity and cause it to relocate to focal adhesions and sites of strain in the actin cytoskeleton. Finally, we demonstrate activated RhoA causes testin to adorn stress fibers and become mechanosensitive. Together, our data show that testins mechanoresponse is regulated in cells and provide new insights into LIM domain protein recognition of the actin cytoskeleton mechanical state. 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. Gradients in cellular metabolism regulate development of tonotopy along the cochlea. In vertebrates with elongated auditory organs, mechanosensory hair cells (HCs) are organised such that complex sounds are broken down into their component frequencies along the basal-to-apical long (tonotopic) axis. To generate the frequency-specific characteristics required at the appropriate positions, nascent HCs must determine their tonotopic properties during development. This relies on complex signalling between the developing HC and its local niche along the axis within the auditory organ. Here we apply NAD(P)H fluorescence lifetime imaging (FLIM) and live imaging of mitochondria to reveal metabolic gradients along the tonotopic axis of the developing cochlea. We further show that re-shaping these gradients during development, by inhibiting cytosolic glucose metabolism, alters normal Bmp7 and Chordin like-1 signalling leading to flattening of tonotopic patterning along the axis. Our work supports a causal link between morphogen signalling and metabolic reprogramming in specifying tonotopic morphologies in auditory HCs. Fecal immunoglobulin A (IgA) and its subclasses in systemic lupus erythematosus patients are nuclear antigen reactive and this feature correlates with gut permeability marker levels Systemic lupus erythematosus (SLE) is characterized by the production of anti-nuclear autoantibodies. Here, for the first time, we show that the abundances of gut permeability marker Zonulin and IgA1- and IgA2-subclasses are significantly higher in the fecal samples of SLE patients compared to HCs. Importantly, IgA-total, IgA1- and IgA2-subclasses from SLE patients showed higher nAg reactivity titers. Notably, we found that not only the nuclear antigen (nAg) reactive fecal IgA1:IgA2 ratio is higher in SLE patients, but also the abundance and nAg reactivity of fecal IgA and subclasses, IgA1 particularly, correlate with the fecal levels of Zonulin, which is produced primarily in the small intestine. These observations suggest that higher amounts of nAg-reactive IgA are produced, particularly, in the proximal gut as indicated by higher IgA1:IgA2 abundance and nAg reactivity ratios and correlation of these features with the Zonulin levels. Analysis of Drosophila cardiac hypertrophy by micro-computerized tomography for genetic dissection of heart growth mechanisms Heart failure is often preceded by pathological cardiac hypertrophy, a thickening of the heart musculature driven by complex gene regulatory and signaling processes. The Drosophila heart has great potential as a genetic model for deciphering the underlying mechanisms of cardiac hypertrophy. However, current methods for evaluating hypertrophy of the Drosophila heart are laborious and difficult to carry out reproducibly. Here we demonstrate that micro-computerized tomography (microCT) is an accessible, highly reproducible method for non-destructive, quantitative analysis of Drosophila heart morphology and size. To validate our microCT approach for analyzing Drosophila cardiac hypertrophy, we show that expression of constitutively active Ras (Ras85DV12), previously shown to cause hypertrophy of the fly heart, results in significant thickening of both adult and larval heart walls when measured from microCT images. We then show using microCT analysis that genetic upregulation of store-operated Ca2+ entry (SOCE) driven by expression of constitutively active Stim (StimCA) or Orai (OraiCA) proteins also results in significant hypertrophy of the Drosophila heart, through a process that specifically depends on Orai Ca2+ influx channels. Intravital imaging of heart contractility revealed significantly reduced end diastolic dimensions in StimCA and OraiCA expressing hearts, consistent with the hypertrophic phenotype. These results demonstrate that increased SOCE activity is an important driver of hypertrophic cardiomyocyte growth, and demonstrate how microCT analysis combined with tractable genetic tools in Drosophila can be used to delineate molecular signaling processes that underlie cardiac hypertrophy and heart failure.

NEW AND NOTEWORTHYGenetic analysis of cardiac hypertrophy in Drosophila holds immense potential for the discovery of new therapeutic targets to prevent and treat heart failure. However, this potential has been hindered by a lack of rapid and effective methods for analysis of heart size in flies. Here we demonstrate that analysis of the Drosophila heart with micro-computerized tomography yields accurate and highly reproducible heart size measurements that can be used to efficiently analyze heart growth and cardiac hypertrophy in Drosophila. Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. Loci for prediction of penicillin and tetracycline susceptibility in Neisseria gonorrhoeae: a genome wide association study BackgroundWhile Neisseria gonorrhoeae poses an urgent public health threat because of increasing antimicrobial resistance, much of the circulating population remains susceptible to historical treatment regimens. Point-of-care diagnostics that report susceptibility could allow for reintroduction of these regimens, but development of such diagnostics has been limited to ciprofloxacin, for which susceptibility can be predicted from a single locus.

MethodsWe assembled a dataset of 12,045 N. gonorrhoeae genomes with phenotypic resistance data for tetracycline (n = 3,611) and penicillin (n = 6,935). Using conditional genome wide association studies (GWAS), we sought to define genetic variants associated with susceptibility to penicillin and tetracycline. We evaluated the sensitivity and specificity of these variants for predicting susceptibility and non-resistance in our collection of gonococcal genomes.

FindingsIn our conditional penicillin GWAS, the presence of a genetic variant defined by a non-mosaic penA allele without an insertion at codon 345 was significantly associated with penicillin susceptibility and had the highest negative effect size of significant variants (p = 5.0 x 10-14, {beta} = -2.5). In combination with the absence of blaTEM, this variant predicted penicillin susceptibility with high specificity (99.8%) and modest sensitivity (36.7%). For tetracycline, the wild type allele at rpsJ codon 57, encoding valine, was significantly associated with tetracycline susceptibility (p = 5.6 x 10-16, {beta} = -1.6) after conditioning on the presence of tetM. The combination of rpsJ codon 57 allele and tetM absence predicted tetracycline susceptibility with high specificity (97.2%) and sensitivity (88.7%).

InterpretationAs few as two genetic loci can predict susceptibility to penicillin and tetracycline in N. gonorrhoeae with high specificity. Molecular point-of-care diagnostics targeting these loci have the potential to increase available treatments for gonorrhea.

FundingNational Institute of Allergy and Infectious Diseases, the National Science Foundation, and the Smith Family Foundation

Research in ContextO_ST_ABSEvidence before this studyC_ST_ABSWe searched PubMed with the terms "Neisseria gonorrhoeae" and "diagnostic" or "assay" plus "penicillin" or "tetracycline" for reports in any language published up to July 1, 2021. We additionally searched for "Neisseria gonorrhoeae" and "genome wide association study". We found that previously proposed molecular diagnostics for penicillin and tetracycline susceptibility either exclusively focused on plasmid-mediated resistance (i.e., targeting blaTEM or tetM) or did not include variants in genes encoding antibiotic targets (e.g., did not include penA or rpsJ). Targets for molecular surveillance have focused on resistance-associated alleles rather than susceptibility-associated alleles. We did not find any previous penicillin or tetracycline conditional genome wide association studies (GWAS) in N. gonorrhoeae.

Added value of this studyTo identify targets for molecular diagnostics that predict penicillin and tetracycline susceptibility, we conducted GWAS conditioning on the presence of plasmid-mediated resistance determinants to detect chromosomal loci with the highest association with susceptibility. We discovered a sequence (penA_01) that differentiates susceptible isolates from those with a resistance-associated insertion at codon 345 and from those with mosaic penA alleles, which is associated with penicillin susceptibility. We also found that rpsJ codon 57 was the chromosomal locus contributing the most to tetracycline susceptibility. The combination of these chromosomal loci and the absence of plasmid encoded determinants predicts penicillin and tetracycline susceptibility with high specificity in both a large global collection of N. gonorrhoeae and a validation dataset consisting of recently published genomes from CDCs Gonococcal Isolate Surveillance Program (GISP) surveillance collected in 2018.

Implications of all the available evidenceThe chromosomal loci penA_01 and rpsJ codon 57 in combination with plasmid loci blaTEM and tetM are candidates for the development of point-of-care molecular diagnostics for penicillin and tetracycline susceptibility. The loci may be combined with the currently available ciprofloxacin susceptibility diagnostics to predict susceptibility to multiple antibiotics. Additionally, our study suggests that conditional GWAS focused on variants associated with susceptibility may be a promising approach to identify minimal sets of loci for molecular diagnostics and surveillance. Classification of daily crop phenology in PhenoCams using deep learning and hidden markov models Near surface cameras, such as those in the PhenoCam network, are a common source of ground truth data in modelling and remote sensing studies. Despite having locations across numerous agricultural sites, few studies have used near surface cameras to track the unique phenology of croplands. Due to management activities, crops do not have a natural vegetation cycle which many phenological extraction methods are based on. For example, a field may experience abrupt changes due to harvesting and tillage throughout the year. A single camera can also record several different plants due to crop rotations, fallow fields, and cover crops. Current methods to estimate phenology metrics from image time series compress all image information into a relative greenness metric, which discards a large amount of contextual information. This can include the type of crop present, whether snow or water is present on the field, the crop phenology, or whether a field lacking green plants consists of bare soil, fully senesced plants, or plant residue. Here we developed a modelling workflow to create a daily time series of crop type and phenology, while also accounting for other factors such as obstructed images and snow covered fields. We used a mainstream deep learning image classification model, VGG16. Deep learning classification models do not have a temporal component, so to account for temporal correlation among images our workflow incorporates a hidden markov model in the post-processing. The initial image classification model had out of sample F1 scores of 0.83-0.85, which improved to 0.86-0.91 after all post-processing steps. The resulting time series show the progression of crops from emergence to harvest, and can serve as a daily, local scale dataset of field states and phenological stages for agricultural research. Common Dynamic Determinants Govern Quorum Quenching Activity in N-terminal Serine Hydrolases Growing concerns about microbial antibiotic resistance have motivated extensive research into ways of overcoming antibiotic resistance. Quorum quenching (QQ) processes disrupt bacterial communication via quorum sensing, which enables bacteria to sense the surrounding bacterial cell density and markedly affects their virulence. Due to its indirect mode of action, QQ is believed to exert limited pressure on essential bacterial functions and may thus avoid inducing resistance. Although many enzymes display QQ activity against various bacterial signaling molecules, their mechanisms of action are poorly understood, limiting their potential optimization as QQ agents. Here we evaluate the capacity of three N-terminal serine hydrolases to degrade N-acyl homoserine lactones that serve as signaling compounds for Gram-negative bacteria. Using molecular dynamics simulations of the free enzymes and their complexes with two signaling molecules of different lengths, followed by quantum mechanics/molecular mechanics molecular dynamics simulations of their initial catalytic steps, we clarify the molecular processes underpinning their QQ activity. We conclude that all three enzymes degrade bacterial signaling molecules via similar reaction mechanisms. Moreover, we experimentally confirmed the activity of two penicillin G acylases from Escherichia coli (ecPGA) and Achromobacter spp. (aPGA), adding these biotechnologically well-optimized enzymes to the QQ toolbox. We also observed enzyme- and substrate-dependent differences in the catalytic actions of these enzymes, arising primarily from the distinct structures of their acyl-binding cavities and the dynamics of their molecular gates. As a consequence, the first reaction step catalyzed by ecPGA with a longer substrate had an elevated energy barrier because its shallow acyl binding site could not accommodate a productive substrate-binding configuration. Conversely, aPGA in complex with both substrates exhibited unfavorable energetics in both reaction steps due to the dynamics of the residues gating the acyl binding cavity entrance. Finally, the energy barriers of the second reaction step catalyzed by Pseudomonas aeruginosa acyl-homoserine lactone acylase with both substrates were higher than in the other two enzymes due to the unique positioning of Arg297{beta} in this enzyme. The discovery of these dynamic determinants will guide future efforts to design robust QQ agents capable of selectively controlling virulence in resistant bacterial species. Global membrane tension is independent of polyacrylamide substrate stiffness Most animal cells are surrounded by a cell membrane and an underlying actomyosin cortex. Both structures are linked, and they are under tension. In-plane membrane tension and cortical tension both influence many cellular processes, including cell migration, division, and endocytosis. However, while actomyosin tension is regulated by substrate stiffness, how membrane tension responds to mechanical substrate properties is currently poorly understood. Here, we probed the effective membrane tension of neurons and fibroblasts cultured on glass and polyacrylamide substrates of varying stiffness using optical tweezers. In contrast to actomyosin-based traction forces, both peak forces and steady state tether forces of cells cultured on hydrogels were independent of substrate stiffness and did not change after blocking myosin II activity using blebbistatin, indicating that tether and traction forces are not directly linked. Peak forces in fibroblasts on hydrogels were about twice as high as those in neurons, indicating stronger membrane-cortex adhesion in fibroblasts. Steady state tether forces were generally higher in cells cultured on hydrogels than on glass, which we explain by a mechanical model. Our results provide new insights into the complex regulation of effective membrane tension and pave the way for a deeper understanding of the biological processes it instructs. Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. Modulation of RNA polymerase processivity affects double-strand break repair in the presence of a DNA end-binding protein Homologous recombination is the predominant double-strand DNA break repair pathway in Escherichia coli. A recent report of non-homologous end joining involving the Ku-like Gam protein of bacteriophage Mu (MuGam) has sparked interest in understanding the functions of DNA end-binding proteins in bacteria. MuGam binds to DNA ends, but how it interferes with DNA repair or transcription in live bacteria remains elusive. In E. coli, RNA polymerase secondary channel interactors, such as the Gre factors, play a role in coordinating transcription with DNA replication and break repair. Here we show that MuGam inhibits break repair by slowing down DNA resection and impeding recombination in living cells. Loss of GreA restores break repair in the presence of MuGam by allowing for increased DNA resection due to a potential release of MuGam from the DNA. Using MuGam as a DNA break sensor, we found that breaks are generated when translation is inhibited, more so in the presence of GreA, supporting the model where transcription-translation uncoupling increases transcription/replication collisions. Significantly, this work reveals that modulation of RNA polymerase can impact DNA break repair in presence of a Ku-like protein. Supergenes origin and maintenance in Atlantic cod Supergenes are sets of genes that are inherited as a single marker and encode complex phenotypes through their joint action. They are identified in an increasing number of organisms, yet their origins and evolution remain enigmatic. In Atlantic cod, four megabase-scale supergenes have been identified and linked to migratory lifestyle and environmental adaptations. Here, we investigate the origin and maintenance of these four supergenes through analysis of whole-genome-sequencing data, including a new long-read-based genome assembly for a non-migratory Atlantic cod individual. We corroborate that chromosomal inversions underlie all four supergenes, and show that they originated at different times between 0.40 and 1.66 million years ago. While we found no evidence for a role of introgression in the origin of the four supergenes, we reveal gene flux between supergene haplotypes with derived and ancestral arrangements, occurring both through gene conversion and double crossover. Our results suggest that supergenes can be maintained over long timescales in the same way as hybridizing species, through the purging of introduced genetic variation. Cis-Regulatory Hubs: a new 3D Model of Complex Disease Genetics with an Application to Schizophrenia The 3-dimensional (3D) conformation of the chromatin creates complex networks of noncoding regulatory regions (distal elements) and promoters impacting gene regulation. Despite the importance of the role of noncoding regions in complex diseases, little is known about their interplay within regulatory hubs and implication in multigenic diseases like schizophrenia. Here we show that cis-regulatory hubs (CRHs) in neurons highlight functional interactions between distal elements and promoters, providing a model to explain epigenetic mechanisms involved in complex diseases. CRHs represent a new 3D model, where distal elements interact to create a complex network of active genes. In a disease context, CRHs highlighted strong enrichments in schizophrenia-associated genes, schizophrenia-associated SNPs and schizophrenia heritability compared to equivalent structures. Finally, CRHs exhibit larger proportions of genes differentially expressed in schizophrenia compared to promoter-distal element pairs or TADs. CRHs thus capture causal regulatory processes improving the understanding of complex disease etiology such as schizophrenia. These multiple lines of genetic and statistical evidence support CRHs as 3D models to study dysregulation of gene expression in complex diseases more generally.

Summary BlurbGenes and their regulatory elements are organized within 3D networks in neurons which model functional structures and explain schizophrenia genetic etiology. Chromosome breaks in breast cancers occur near herpes tumor virus sequences and explain why the cancer comes back One explanation for why breast cancer has a relentless tendency to come back after treatment is that breast cancer phenomes are fluid, so it is difficult for therapy to target new phenotypes that arise. In breast cancer, viral interactions cause chromosome breakage and fusion, producing anaphase bridges containing two centromeres. Cell division then leads to more chromosome breaks, more mutations, and new cancer phenotypes. Analyses of public data from about 2100 breast cancers implicate Epstein-Barr virus (EBV or Human Herpes Virus 4) in the cycles that break chromosomes and produce anaphase bridges. EBV generates breakage cycles with highly fragmented chromosomes (chromothripsis) as an inherent feature. These episodes continually produce new cancer driver mutations. DNA at chromothripsis boundaries in breast cancers consistently matches EBV tumor variants. Breakpoints are often near known EBV binding sites, and breast cancer shares an EBV methylation signature with other known EBV cancers. Breast cancer breakpoints on all chromosomes cluster around the same positions as in cancers definitively associated with EBV infection (nasopharyngeal cancer and Burkitts lymphoma). The status of EBV in the microbiome and the presence of anaphase bridges are phenotypic markers of breast cancer that are not explicitly targeted by chemotherapy. The breast cancer phenome continually supplies new cancer driver genes that foster breast cancer recurrence.

SummaryBreaks in breast cancer chromosomes and in fragments of anaphase bridges occur mainly near regions that resemble Epstein-Barr virus (EBV) tumor variants

Extensive breaks are an inherent outcome of anaphase bridges and the reassembly of highly fragmented chromosomes continually generates new cancer driver mutations and new cancer phenotypes.

Treatment does not explicitly target the underlying breakage and reassembly cycles so breast cancer often comes back.

The breast cancer phenome is not fixed and can continuously generate new phenotypes

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=154 HEIGHT=200 SRC="FIGDIR/small/467751v4_ufig1.gif" ALT="Figure 1">
View larger version (38K):
org.highwire.dtl.DTLVardef@6537cforg.highwire.dtl.DTLVardef@132fa50org.highwire.dtl.DTLVardef@1cfcb04org.highwire.dtl.DTLVardef@11979aa_HPS_FORMAT_FIGEXP  M_FIG C_FIG No impact of biocontrol agent's predation cues on development time or size of surviving Aedes albopictus under optimal nutritional availability Cyclopoid copepods have been applied successfully to limit populations of highly invasive Aedes albopictus mosquitoes that can transmit diseases of public health importance. However, there is concern that changes in certain mosquito traits, induced by exposure to copepod predation, might increase the risk of disease transmission. In this study, third instar Ae. albopictus larvae, hereafter referred to as "focal individuals", were exposed to Megacyclops viridis predator cues associated with both consumption of newly-hatched mosquito larvae and attacks on focal individuals. The number of newly-hatched larvae surrounding each focal larva was held constant to control for density effects on size, and the focal individuals day of pupation and wing length were recorded for each replicate.

Exposing late instar Ae. albopictus to predation decreased their chances of surviving to adulthood, and three focal larvae that died in the predator treatment showed signs of melanization, indicative of wounding. Among surviving focal Ae. albopictus, no significant difference in either pupation day or wing length was observed due to copepod predation. The absence of a significant sublethal impact from M. viridis copepod predation on surviving later-stage larvae in this analysis supports the use of M. viridis as a biocontrol agent against Ae. albopictus.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=106 SRC="FIGDIR/small/460697v2_ufig1.gif" ALT="Figure 1">
View larger version (9K):
org.highwire.dtl.DTLVardef@1abaabborg.highwire.dtl.DTLVardef@1561651org.highwire.dtl.DTLVardef@d325a7org.highwire.dtl.DTLVardef@4d1830_HPS_FORMAT_FIGEXP  M_FIG C_FIG 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. Royal knifefish generate powerful suction feeding through large neurocranial elevation and high epaxial muscle power Suction feeding in ray-finned fishes involves powerful buccal cavity expansion to accelerate water and food into the mouth. Previous XROMM studies in largemouth bass (Micropterus salmoides), bluegill sunfish (Lepomis macrochirus), and channel catfish (Ictalurus punctatus) have shown that more than 90% of suction power in high performance strikes comes from the axial musculature. Thus, the shape of the axial muscles and skeleton may impact suction feeding mechanics. Royal knifefish (Chitala blanci) have an unusual postcranial morphology, with a ventrally flexed vertebral column and relatively large mass of epaxial muscle. Based on their body shape, we hypothesized that royal knifefish would generate high power strikes by utilizing large neurocranial elevation, vertebral column extension, and epaxial shortening. As predicted, C. blanci generated high suction expansion power compared to the other three species studied to date (up to 160 W), which was achieved by increasing both the rate of volume change and the intraoral subambient pressure. The large epaxial muscle (25% of body mass) shortened at high velocities to produce large neurocranial elevation and vertebral extension (up to 41 deg, combined), as well as high muscle mass-specific power (up to 800 W kg-1). For the highest power strikes, axial muscles generated 95% of the power, and 64% of the axial muscle mass consisted of the epaxial muscles. The epaxial-dominated suction expansion of royal knifefish supports our hypothesis that postcranial morphology may be a strong predictor of suction feeding biomechanics.

SUMMARY STATEMENTRoyal knifefish rely on their distinct postcranial morphology--with a curved vertebral column and large dorsal body muscles--to produce large neurocranial elevation and powerful suction feeding. Accuracy and Power Analysis of Social Networks Built From Count Data O_LIPower analysis is used to estimate the probability of correctly rejecting a null hypothesis for a given statistical model and dataset. Conventional power analyses assume complete information, but the stochastic nature of behavioural sampling can mean that true and estimated networks are poorly correlated. Power analyses do not currently take the effect of sampling into account. This could lead to inaccurate estimates of statistical power, potentially yielding misleading results.
C_LIO_LIHere we develop a method for computing network correlation: the correlation between an estimated social network and its true network, using a Gamma-Poisson model of social event rates for networks constructed from count data. We use simulations to assess how the level of network correlation affects the power of nodal regression analyses. We also develop a generic method of power analysis applicable to any statistical test, based on the concept of diminishing returns.
C_LIO_LIWe demonstrate that our network correlation estimator is both accurate and moderately robust to its assumptions being broken. We show that social differentiation, mean social event rate, and the harmonic mean of sampling times positively impacts the strength of network correlation. We also show that the required level of network correlation to achieve a given power level depends on many factors, but that 0.80 network correlation usually corresponds to around 0.80 power for nodal regression in ideal circumstances.
C_LIO_LIWe provide guidelines for using our network correlation estimator to verify the accuracy of networks built from count data, and to conduct power analysis. This can be used prior to data collection, in post hoc analyses, or even for subsetting networks in dynamic network analysis. The network correlation estimator and custom power analysis methods have been made available as an R package.
C_LI Synthetic periphyton as a model system to understand species dynamics in complex microbial freshwater communities Phototrophic biofilms, also known as periphyton, are microbial freshwater communities that drive crucial ecological processes in streams and lakes. Gaining a deep mechanistic understanding of the biological processes occurring in natural periphyton remains challenging due to the high complexity and variability of such communities. To address this challenge, we rationally developed a workflow to construct a synthetic community by co-culturing 26 phototrophic species (i.e., diatoms, green algae and cyanobacteria) that were inoculated in a successional sequence to create a periphytic biofilm on glass slides. We show that this community is diverse, stable and highly reproducible in terms of microbial composition, function and 3D spatial structure of the biofilm. We also demonstrate the ability to monitor microbial dynamics at the single species level during periphyton development and how their abundances are impacted by stressors such as increased temperature and a herbicide, singly and in combination. Overall, such a synthetic periphyton, grown under controlled conditions, can be used as a model system for theory testing through targeted manipulation. Microorganisms Characterization in Anaerobic Solventogenic Processing of Volatile Fatty Acids 1This work evaluates four different physicochemical pretreatments (acidic, thermal, acidic-thermal and thermal-acidic) on an anaerobic inoculum used for alcohol production from acetate and butyrate. All experiments were conducted in single batches using acetate and butyrate as substrates at 30 {degrees}C and with a pressurized headspace of pure H2 at 2.15 atm (218.2 MPa). Thermal and acidic-thermal pretreatments lead to higher production of both ethanol and butanol. Mathematical modelling shows that the highest attainable concentrations of ethanol and butanol produced were 122 mg L-1 and 97 mg L-1 for the thermal pretreatment (after 17.5 days) and 87 mg L-1 and 143 mg L-1 for the acidic-thermal pretreatment (after 18.9 days). Acetate was produced in all assays. Thermodynamic data indicated that a high H2 partial pressure favoured solventogenic metabolic pathways. Finally, sequencing data showed that both thermal and acidic-thermal pretreatments selected mainly the bacterial genera Pseudomonas, Brevundimonas and Clostridium. The acidic-thermal pretreatment selected a bacterial community more adapted to the conversion of acetate and butyrate into ethanol and butanol, respectively. Thermal-acidic pretreatment was unstable, showing significant variability between replicates. Acidic pretreatment showed the lowest alcohol production. Towards a cumulative science of vocal markers of autism: a cross-linguistic meta-analysis-based investigation of acoustic markers in American and Danish autistic children Acoustic atypicalities in speech production are argued to be potential markers of clinical features in Autism Spectrum Disorder (ASD). A recent meta-analysis highlighted shortcomings in the field, in particular small sample sizes and study heterogeneity (Fusaroli et al., 2017). We showcase a cumulative (i.e., explicitly building on previous studies both conceptually and statistically) yet self-correcting (i.e., critically assessing the impact of cumulative statistical techniques) approach to prosody in ASD to overcome these issues.

We relied on the recommendations contained in the meta-analysis to build and analyze a cross-linguistic corpus of multiple speech productions in 77 autistic and 72 neurotypical children and adolescents (>1000 recordings in Danish and US English). We used meta-analytically informed and skeptical priors, with informed priors leading to more generalizable inference. We replicated findings of a minimal cross-linguistically reliable distinctive acoustic profile for ASD (higher pitch and longer pauses) with moderate effect sizes. We identified novel reliable differences between the two groups for normalized amplitude quotient, maxima dispersion quotient, and creakiness. However, the differences were small, and there is likely no one acoustic profile characterizing all autistic individuals. We identified reliable relations of acoustic features with individual differences (age, gender), and clinical features (speech rate and ADOS sub-scores).

Besides cumulatively building our understanding of acoustic atypicalities in ASD, the study shows how to use systematic reviews and meta-analyses to guide the design and analysis of follow-up studies. We indicate future directions: larger and more diverse cross-linguistic datasets, focus on heterogeneity, self-critical cumulative approaches and open science.

Lay SummaryAutistic individuals are reported to speak in distinctive ways. Distinctive vocal production can affect social interactions and social development and could represent a noninvasive way to support the assessment of ASD. We systematically checked whether acoustic atypicalities highlighted in previous articles could be actually found across multiple recordings and two languages. We find a minimal acoustic profile of ASD: higher pitch, longer pauses, increased hoarseness and creakiness of the voice. However, there is much individual variability (by age, sex, language, and clinical characteristics). This suggests that the search for one common "autistic voice" might be naive and more fine-grained approaches are needed. The first sheep graph-based pan-genome reveals the spectrum of structural variations and their effects on tail phenotypes Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here, we assembled 26 haplotype-resolved genome assemblies from 13 genetically diverse sheep using PacBio HiFi sequencing. We constructed a graph-based ovine pan-genome and discovered 142,422 biallelic insertions and deletions, 7,028 divergent alleles and 13,419 multiallelic variations. We then used a graph-based approach to genotype the biallelic SVs in 684 individuals from 45 domestic breeds and two wild species. Integration with RNA-seq data allows to identify candidate expression-associated SVs. We demonstrate a direct link of SVs and phenotypes by localizing the putative causative insertion in HOXB13 gene responsible for the long-tail trait and identifying multiple large SVs associated with the fat-tail. Beyond generating a benchmark resource for ovine structural variants, our study highlights that animal genetic research will greatly benefit from using a pan-genome graph rather than a single reference genome.